CN115869005B - 超声图像的干扰抑制方法、超声设备和存储介质 - Google Patents
超声图像的干扰抑制方法、超声设备和存储介质 Download PDFInfo
- Publication number
- CN115869005B CN115869005B CN202111131090.4A CN202111131090A CN115869005B CN 115869005 B CN115869005 B CN 115869005B CN 202111131090 A CN202111131090 A CN 202111131090A CN 115869005 B CN115869005 B CN 115869005B
- Authority
- CN
- China
- Prior art keywords
- signal
- ultrasonic
- empty
- ultrasound
- echo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000001629 suppression Effects 0.000 title claims abstract description 19
- 230000000737 periodic effect Effects 0.000 claims abstract description 56
- 238000002604 ultrasonography Methods 0.000 claims description 80
- 238000001914 filtration Methods 0.000 claims description 33
- 238000004458 analytical method Methods 0.000 claims description 16
- 239000000523 sample Substances 0.000 claims description 14
- 238000010219 correlation analysis Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 18
- 238000004590 computer program Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- 241001203771 Eudonia echo Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本申请涉及超声图像处理技术领域,提供一种超声图像的干扰抑制方法、超声设备和存储介质,用于解决相关技术中一些干扰信号存在的频段通常会覆盖超声设备的工作频段而传统的滤波器无法有效对滤除该类干扰信号的问题。在本申请中,通过在超声信号扫描序列中添加空扫描的方式能够检测并获取周期性干扰信号,然后在超声信号扫描序列中滤除该干扰信号,能够达到对周期性干扰信号的检测和抑制的目的。
Description
技术领域
本申请涉及超声图像处理技术领域,尤其涉及一种超声图像的干扰抑制方法、超声设备和存储介质。
背景技术
通过控制探头发射超声信号(亦称之为超声波束)和接受超声波来得到超声图像。接收到的超声回波信号在前端经过滤波放大后通过波束合成即可得到原始的超声图像。
在医院的工作环境中通常会有像X光,CT,高频电刀等各种医疗设备同时运行,这些都会对超声设备产生电磁干扰。干扰信号经前端放大后通过波束合成会最终显示在超声图像上,进而对病灶部位产生噪声影响医生的诊断。
而,一些干扰信号存在的频段通常会覆盖超声设备的工作频段,因此传统的滤波器无法有效对滤除该类干扰信号。
发明内容
本申请实施例提供一种超声图像的干扰抑制方法、超声设备和存储介质,用于解决相关技术中一些干扰信号存在的频段通常会覆盖超声设备的工作频段而传统的滤波器无法有效对滤除该类干扰信号的问题。
第一方面,本申请提供一种超声图像的干扰抑制方法,所述方法包括:
获取多条第一空扫描线的回波数据;
对所述多条第一空扫描线的回波数据进行信号分析;
若分析到周期性干扰信号,则从超声信号的超声回波信号中滤除所述周期性干扰信号。
可选的,所述对所述多条第一空扫描线的回波数据进行信号分析,包括:
对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差;
若相邻第一空扫描线的延迟差在指定范围内波动,则确定分析到所述周期性干扰信号。
可选的,所述从超声信号的超声回波信号中滤除所述周期性干扰信号之前,所述方法还包括:
确定所述多条第一空扫描线的回波信号的信号强度高于预设强度阈值。
可选的,所述从超声信号的超声回波信号中滤除所述周期性干扰信号,包括:
针对每条超声信号的超声回波信号,基于所述第一空扫描线的回波信号的延迟差对所述第一空扫描线的回波信号进行调整,得到与所述超声回波信号对应的干扰信号;
将所述干扰信号从所述超声回波信号中滤除。
可选的,一帧超声图像对应多条超声回波信号;所述从超声信号的超声回波信号中滤除所述周期性干扰信号:
在一帧图像的超声扫描序列中添加第二空扫描线,得到所述第二空扫描线的回波数据作为干扰信号;其中所述超声扫描序列中至少一条超声扫描线对应一条所述第二扫描线;
基于所述周期性干扰信号的延迟差对所述干扰信号进行调整,得到待滤除信号;
将所述一帧图像的与所述干扰信号对应的所述至少一条超声扫描线的回波信号进行滤波处理,得到滤除所述周期性干扰信号的超声回波信号。
可选的,针对B型模式,每条超声信号对应一条空扫描线;
针对C型模式,一组超声信号对应一条空扫描线。
可选的,所述一组超声信号包括A线、B线和C线,所述一组超声信号对应的一条扫描线为D线。
可选的,所述对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差,包括:
采用频域自相关分析方法对相邻第一空扫描线的回波信号进行相关性分析,得到自相关结果信号;所述自相关结果信号用于表述延迟差和自相关值之间的对应关系;
从所述自相关结果信号中搜索自相关值的最大值;
获取所述最大值对应的延迟差作为所述相邻第一空扫描线之间的延迟差。
第二方面,本申请还提供一种超声设备,包括:处理器、存储器、显示单元和探头;
探头,用于发射超声信号;
显示单元,用于显示超声图像;
处理器,分别与所述探头以及所述显示单元相连接,被配置为执行:
获取多条第一空扫描线的回波数据;
对所述多条第一空扫描线的回波数据进行信号分析;
若分析到周期性干扰信号,则从超声信号的超声回波信号中滤除所述周期性干扰信号。
可选的,执行所述对所述多条第一空扫描线的回波数据进行信号分析,所述处理器具体被配置为执行:
对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差;
若相邻第一空扫描线的延迟差在指定范围内波动,则确定分析到所述周期性干扰信号。
可选的,执行所述从超声信号的超声回波信号中滤除所述周期性干扰信号之前,所述处理器还被配置为执行:
确定所述多条第一空扫描线的回波信号的信号强度高于预设强度阈值。
可选的,执行所述从超声信号的超声回波信号中滤除所述周期性干扰信号,所述处理器具体被配置为执行:
针对每条超声信号的超声回波信号,基于所述第一空扫描线的回波信号的延迟差对所述第一空扫描线的回波信号进行调整,得到与所述超声回波信号对应的干扰信号;
将所述干扰信号从所述超声回波信号中滤除。
可选的,一帧超声图像对应多条超声回波信号;执行所述从超声信号的超声回波信号中滤除所述周期性干扰信号,所述处理器具体被配置为执行:
在一帧图像的超声扫描序列中添加第二空扫描线,得到所述第二空扫描线的回波数据作为干扰信号;其中所述超声扫描序列中至少一条超声扫描线对应一条所述第二扫描线;
基于所述周期性干扰信号的延迟差对所述干扰信号进行调整,得到待滤除信号;
将所述一帧图像的与所述干扰信号对应的所述至少一条超声扫描线的回波信号进行滤波处理,得到滤除所述周期性干扰信号的超声回波信号。
可选的,针对B型模式,每条超声信号对应一条空扫描线;
针对C型模式,一组超声信号对应一条空扫描线。
可选的,所述一组超声信号包括A线、B线和C线,所述一组超声信号对应的一条扫描线为D线。
可选的,执行所述对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差,所述处理器具体被配置为执行:
采用频域自相关分析方法对相邻第一空扫描线的回波信号进行相关性分析,得到自相关结果信号;所述自相关结果信号用于表述延迟差和自相关值之间的对应关系;
从所述自相关结果信号中搜索自相关值的最大值;
获取所述最大值对应的延迟差作为所述相邻第一空扫描线之间的延迟差。
第三方面,本申请一实施例还提供了一种计算机可读存储介质,当所述计算机可读存储介质中的指令由电子设备的处理器执行时,使得电子设备能够执行如本申请第一方面中提供的任一方法。
第四方面,本申请一实施例提供了一种计算机程序产品,包括计算机程序,所述计算机程序被处理器执行时实现如本申请第一方面中提供的任一方法。
本申请的实施例提供的技术方案至少带来以下有益效果:通过在超声信号扫描序列中添加空扫描的方式能够检测并获取周期性干扰信号,然后在超声信号扫描序列中滤除该干扰信号,能够达到对周期性干扰信号的检测和抑制的目的。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本公开。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例中所需要使用的附图作简单地介绍,显而易见地,下面所介绍的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的超声设备的框架示意图;
图2为本申请一实施例提供的超声设备实现超声图像的干扰抑制的原理示意图;
图3为本申请一实施例提供的超声图像的干扰抑制方法;
图4为本申请一实施例提供的空扫描线的示意图;
图5为本申请一实施例提供的进行延迟调整的示意图;
图6为本申请一实施例提供的添加空扫描线方式的另一示意图;
图7为本申请一实施例提供的添加空扫描线方式的又一示意图;
图8为本申请一实施例提供的添加空扫描线方式的另一示意图;
图9为本申请一实施例提供的添加空扫描线方式的又一示意图。
具体实施方式
为了使本领域普通人员更好地理解本申请的技术方案,下面将结合附图,对本申请实施例中的技术方案进行清楚、完整地描述。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应所述理解这样使用的数据在适当情况下可以互换,以便这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
以下,对本申请实施例中的部分用语进行解释说明,以便于本领域技术人员理解。
(1)本申请实施例中术语“多个”是指两个或两个以上,其它量词与之类似。
(2)“和/或”,描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
(3)超声设备,本申请实施例中超声设备指支持B型模式和C型模式的超声成像设备。用户在显示界面中可以通过***菜单或者控制按键选择解剖B型模式或C型模式。
鉴于相关技术中难以过滤掉超声统的工作频段内的干扰信号的问题,本申请实施例提供了一种超声图像的干扰抑制方法、超声设备和存储介质。
对超声***产生的干扰的信号分为周期性干扰与非周期性干扰。周期性干扰通常由高频电刀等医疗设备引起,这类干扰信号的功率通常较大,对超声图像的影响较为明显。并且干扰一般不会是单一频率,通常由基频信号和其各次谐波构成,这些频段往往会覆盖超声工作频段因此传统滤波器难以滤除。
例如,超声回波信号在形成最终超声图像之前会经历一系列的滤波操作,通常包含前端电路中的低通以及高通滤波器。在经过解调后还有动态滤波的操作。在图像处理领域还会经历一系列噪声抑制及图像增强操作。这些流程都会对干扰信号产生一定的效果,但如果干扰信号所包含成分与超声信号处于同一频带那么这些操作都不会对干扰有抑制作用,然而这又是非常常见的情况。
为了能够滤除在超声***频段内的周期性干扰信号,本申请实施例中提出在正常扫描序列(即发出超声信号的扫描序列)中添加空扫描线来获取干扰信号的信息,最后可通过相位调整的手段进行滤除。其中,空扫描线即可控制探头不发射超声波束但是仍接收回波信号。也可理解为对探头的激励很小甚至没有,导致超声探头达到不发射超声波束,但是能接收回波信号的效果。
下面结合附图对本申请实施例提供的超声设备和超声图像的干扰抑制方法进行说明。
参见图1所示,为本申请实施例提供的超声设备的结构框图。
应该理解的是,图1所示超声设备100仅是一个范例,并且超声设备100可以具有比图1中所示的更多的或者更少的部件,可以组合两个或多个的部件,或者可以具有不同的部件配置。图中所示出的各种部件可以在包括一个或多个信号处理和/或专用集成电路在内的硬件、软件、或硬件和软件的组合中实现。
图1中示例性示出了根据示例性实施例中超声设备100的硬件配置框图。
如图1所示,超声设备100例如可以包括:处理器110、存储器120、显示单元130和探头140;其中,
探头140,用于发射超声波束;
显示单元130,用于显示解剖M型图像;
存储器120被配置为存储用于超声成像所需的数据,可包括软件程序,应用界面数据等;
处理器110,分别与所述探头140、所述显示单元130和存储器120相连接,被配置为执行:
获取多条第一空扫描线的回波数据;
对所述多条第一空扫描线的回波数据进行信号分析;
若分析到周期性干扰信号,则从超声信号的超声回波信号中滤除所述周期性干扰信号。
可选的,所述对所述多条第一空扫描线的回波数据进行信号分析,包括:
对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差;
若相邻第一空扫描线的延迟差在指定范围内波动,则确定分析到所述周期性干扰信号。
可选的,所述从超声信号的超声回波信号中滤除所述周期性干扰信号之前,所述方法还包括:
确定所述多条第一空扫描线的回波信号的信号强度高于预设强度阈值。
可选的,所述从超声信号的超声回波信号中滤除所述周期性干扰信号,包括:
针对每条超声信号的超声回波信号,基于所述第一空扫描线的回波信号的延迟差对所述第一空扫描线的回波信号进行调整,得到与所述超声回波信号对应的干扰信号;
将所述干扰信号从所述超声回波信号中滤除。
可选的,一帧超声图像对应多条超声回波信号;所述从超声信号的超声回波信号中滤除所述周期性干扰信号:
在一帧图像的超声扫描序列中添加第二空扫描线,得到所述第二空扫描线的回波数据作为干扰信号;其中所述超声扫描序列中至少一条超声扫描线对应一条所述第二扫描线;
基于所述周期性干扰信号的延迟差对所述干扰信号进行调整,得到待滤除信号;
将所述一帧图像的与所述干扰信号对应的所述至少一条超声扫描线的回波信号进行滤波处理,得到滤除所述周期性干扰信号的超声回波信号。
可选的,针对B型模式,每条超声信号对应一条空扫描线;
针对C型模式,一组超声信号对应一条空扫描线。
可选的,所述一组超声信号包括A线、B线和C线,所述一组超声信号对应的一条扫描线为D线。
可选的,所述对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差,包括:
采用频域自相关分析方法对相邻第一空扫描线的回波信号进行相关性分析,得到自相关结果信号;所述自相关结果信号用于表述延迟差和自相关值之间的对应关系;
从所述自相关结果信号中搜索自相关值的最大值;
获取所述最大值对应的延迟差作为所述相邻第一空扫描线之间的延迟差。
图2为根据本申请一个实施例的应用原理的示意图。其中,该部分可由图1所示超声设备的部分模块或功能组件实现,下面将仅针对主要的部件进行说明,而其它部件,如存储器、控制器、控制电路等,此处将不进行赘述。
如图2所示,应用环境中可以包括用户界面210、用于显示所述用户界面的显示单元220以及处理器230。
显示单元220可以包括显示面板221、背光组件222。其中,显示面板321被配置为对超声图像进行显示,背光组件222位于显示面板221背面,背光组件222可以包括多个背光分区(图中未示出),各背光分区可以发光,以点亮显示面板221。
处理器230可以被配置为控制背光组件222中各背光分区的背光源亮度,以及控制探头发射超声波束和接收回波信号。
其中,处理器230可以包括空扫描线处理单元和超声信号处理单元。其中,空扫描线处理单元用于获取空扫描线的回波数据并对其进行分析,得到周期性干扰信号。信号处理单元用于对超声图像的超声回波信号进行滤波处理以滤除超声回波信号中的周期性干扰信号。
如图3所示,为本申请实施例中超声图像的干扰抑制方法的流程示意图,包括以下步骤:
在步骤301中,获取多条第一空扫描线的回波数据。
在一些实施例中,可以每帧超声图像之前设置两条甚至更多条空扫描线,并得到每条空扫描线的回波数据以便于进行分析。
在另一些实施例中,还可以在每次超声成像之前,设置两条甚至更多条空扫描线,并得到每条空扫描线的回波数据以便于进行分析。也即,可以多帧超声图像共同对应同一第一空扫描线。
在步骤302中,对所述多条第一空扫描线的回波数据进行信号分析,以便于确定第一空扫描线的回波数据是否为周期性干扰信号。
在一些实施例中,所述对所述多条第一空扫描线的回波数据进行信号分析,可实施为基于周期信号的延迟差基本不变的特点进行分析,例如可实施为:
对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差;若相邻第一空扫描线的延迟差在指定范围内波动,则确定分析到所述周期性干扰信号。而后,针对每条超声信号的扫描线分别调整空扫描线的超声回波信号的延迟后进行相减即可将干扰信号滤除。
例如,如图4所示,可以设置3条扫描线,分别为1号、2号和3号扫描线。然后计算1号扫描线和2号扫描线之间的延迟差1,并计算2号扫描线和3号扫描线之间的延迟差2。比较延迟差1和延迟差2,若两个延迟差之间相同或差距小于预设差距,则确定延迟差基本不变。由此,可知得到的干扰信号为周期性干扰信号。由此,可以后续从超声图像的超声回波信号中滤除得到的周期性干扰信号。
在一些实施例中,可采用频域自相关分析方法对相邻第一空扫描线的回波信号进行相关性分析,得到自相关结果信号;所述自相关结果信号用于表述延迟差和自相关值之间的对应关系;然后,从所述自相关结果信号中搜索自相关值的最大值;最后,获取所述最大值对应的延迟差作为所述相邻第一空扫描线之间的延迟差。
如,通过空扫描得到的数据后,需要分析计算各第一空扫描线之间的延迟差来得到。由于干扰信号经常包含基波及其各次谐波分量,因此通过计算频谱的方式较为复杂。本神器在计算延迟时采用频域自相关的方式可有效降低计算复杂度,以相邻的1号扫描线和2号扫描线为例延迟差的计算方式如公式(1)所示:
Rcorr(τ)=IFFT(FFT(X1)*CONJ(FFT(x2))) (1)
在公式(1)中,X1为1号空扫描线的回波数据,X2为2号空扫描线的回波数据,Rcorr(τ)表示τ索引的延迟差对应的自相关值。
其他空扫描线以此类推。延迟差通过自相关结果最大值的索引得到,在若干条第一空扫描线之间,若延迟差基本相同且达到一定数目时则说明当前环境存在周期性干扰。针对每一条常规扫描信号的超声回波信号,在得到该扫描回波信号时理论上有其对应的周期性干扰信号。
实施时,可以利用延迟差信息调整常规扫描线的干扰信号,并与其相减即可消除常规扫描信号中所存在的干扰信号。
在一些实施例中,由于信号强度较高的周期性干扰信号对超声图像的干扰较大,需要抑制,故此在实施时,可针对信号强度较高的信号进行抑制。可实施为对在抑制干扰信号前,可确定所述多条第一空扫描线的回波信号的信号强度是否高于预设强度阈值。若高于预设强度阈值,则进一步分析该回波信号是否为周期性干扰信号。由此可以有针对性的进行干扰抑制。
在步骤303中,若分析到周期性干扰信号,则从超声信号的超声回波信号中滤除所述周期性干扰信号。
一种可能的实施方式为基于第一空扫描线的回波信号预测每条常规扫描线对应的干扰信号然后进行抑制。如可实施为:
针对每条超声信号的超声回波信号,基于所述第一空扫描线的回波信号的延迟差对所述第一空扫描线的回波信号进行调整,得到与所述超声回波信号对应的干扰信号;然后将所述干扰信号从所述超声回波信号中滤除。
如图5所示,假设一帧超声图像需要n条超声信号扫描线,在第一帧超声图像之前添加m条空扫描线,得到这m条空扫描线各自的回波信号。然后针对其中相邻的任意两条空扫描线,确定其回波信号之间的延迟差。若多个延迟差基本不变说明为周期性干扰信号,若信号强度高于预设强度阈值,则该周期性干扰信号对超声图像影响较大需要抑制。而后,采用m条空扫描线的回波信号(即干扰信号)及其延迟差来估计1-n条超声信号扫描线中每条扫描线对应的干扰信号,然后分别进行干扰抑制。如基于延迟差得到第1条超声信号扫描线的干扰信号1(如图6所示),则从第1条超声信号扫描线的超声回波信号中滤除干扰信号1,实施时直接相减即可。类似的,如图6所示,得到第2条超声信号扫描线的干扰信号2,则从第2条超声信号扫描线的超声回波信号中滤除干扰信号2,实施时直接相减即可。以此类推,能够得到一帧图像中各条超声信号扫描线对应的干扰信号并进行抑制。然后采用干扰抑制后的超声信号扫描线的回波信号进行成像。即采用图5中超声信号扫描线1-n的回波信号成像。
实施时,可以如图5所示每一帧超声图像之前都添加m条空扫描线,也可以多帧超声图像之前添加m条空扫描线,均适用于本申请实施例。当然,也可以在每帧超声图像之后,或多帧超声图像之后添加m条空扫描线。
在另一种可能的实施方式中,干扰信号也会受电磁环境的影响发生一些变动。本申请实施例中为了提高抑制干扰信号的精度可以在一帧图像的扫描序列中添加多条空扫描线来实时的获取成像过程中的干扰信号。
可实施为,一帧超声图像对应多条超声回波信号;其中所述超声扫描序列中至少一条超声扫描线对应一条所述第二扫描线。如图7所示,其中虚线为空扫描线,实线为超声信号扫描线。除第一个扫描位置(超声图像之前)两条扫描线皆为空扫描外,其他位置超声信号扫描线与空扫描线交替。其中,一条超声信号扫描线和其对应的一条空扫描线构成一组扫描线。在这样的扫描方式下可以使用前两条空扫描线的回波信号来计算干扰信号的延迟差,利用得到的延迟差对每一组扫描线中的空扫描线的回波信号进行延迟调整,调整之后与该组内的超声信号扫描线的超声回波信号进行干扰抑制。如相减即可实现对干扰信号的有效抑制,无需要计算出干扰信号的准确周期。
实施时,也可以多条超声信号扫描线构成一组超声信号扫描线,一组超声信号扫描线和一条空扫描线构成一组扫描线。如图8所示,前两条扫描线均为空扫描线用于确定干扰信号的延迟差,而后没两条超声信号扫描线对应一条空扫描线构成一组扫描线。一组扫描线内空扫描线用于对超声信号扫描线进行干扰抑制。
针对C型模式,一组超声信号可以为一次循环扫描所需的扫描线。例如如图9所示,针对一帧超声图像,第一个位置均为空扫描线用于计算延迟差,其余位置中一组超声信号包括A线、B线和C线,然后在D线进行空扫描。
需要说明的是,由于不同扫描线的信号来源于不同的阵元且干扰信号存在一定波动,因此直接相减的方式可能无法完整的滤除干扰。本专利滤除干扰的方式还可以使用并不限于自适应滤波等方式。
综上所述,本申请实施例中,能够在扫描序列中添加空扫描对周期性干扰信号进行检测和抑制,实施方式简单易行,能够有效进行周期性干扰信号的抑制。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。
Claims (9)
1.一种超声图像的干扰抑制方法,其特征在于,所述方法包括:
获取多条第一空扫描线的回波数据;
对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差,若相邻第一空扫描线的延迟差在指定范围内波动,则确定分析到周期性干扰信号;
从超声信号的超声回波信号中滤除所述周期性干扰信号。
2.根据权利要求1所述的方法,其特征在于,所述从超声信号的超声回波信号中滤除所述周期性干扰信号之前,所述方法还包括:
确定所述多条第一空扫描线的回波信号的信号强度高于预设强度阈值。
3.根据权利要求1所述的方法,其特征在于,所述从超声信号的超声回波信号中滤除所述周期性干扰信号,包括:
针对每条超声信号的超声回波信号,基于所述第一空扫描线的回波信号的延迟差对所述第一空扫描线的回波信号进行调整,得到与所述超声回波信号对应的干扰信号;
将所述干扰信号从所述超声回波信号中滤除。
4.根据权利要求1所述的方法,其特征在于,一帧超声图像对应多条超声回波信号;所述从超声信号的超声回波信号中滤除所述周期性干扰信号:
在一帧图像的超声扫描序列中添加第二空扫描线,得到所述第二空扫描线的回波数据作为干扰信号;其中所述超声扫描序列中至少一条超声扫描线对应一条所述第二空扫描线;
基于所述周期性干扰信号的延迟差对所述干扰信号进行调整,得到待滤除信号;
将所述一帧图像的与所述干扰信号对应的所述至少一条超声扫描线的回波信号进行滤波处理,得到滤除所述周期性干扰信号的超声回波信号。
5.根据权利要求4所述的方法,其特征在于,
针对B型模式,每条超声信号对应一条空扫描线;
针对C型模式,一组超声信号对应一条空扫描线。
6.根据权利要求5所述的方法,其特征在于,所述一组超声信号包括A线、B线和C线,所述一组超声信号对应的一条扫描线为D线。
7.根据权利要求1所述的方法,其特征在于,所述对所述多条第一空扫描线的回波数据的延迟差进行分析,得到相邻第一空扫描线之间的延迟差,包括:
采用频域自相关分析方法对相邻第一空扫描线的回波信号进行相关性分析,得到自相关结果信号;所述自相关结果信号用于表述延迟差和自相关值之间的对应关系;
从所述自相关结果信号中搜索自相关值的最大值;
获取所述最大值对应的延迟差作为所述相邻第一空扫描线之间的延迟差。
8.一种超声设备,其特征在于,包括:处理器、存储器、显示单元和探头;
探头,用于发射超声信号;
显示单元,用于显示超声图像;
处理器,分别与所述探头以及所述显示单元相连接,被配置为执行如权利要求1-7中任一所述的方法。
9.一种计算机可读存储介质,其特征在于,当所述计算机可读存储介质中的指令由电子设备的处理器执行时,使得所述电子设备能够执行如权利要求1-7中任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111131090.4A CN115869005B (zh) | 2021-09-26 | 2021-09-26 | 超声图像的干扰抑制方法、超声设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111131090.4A CN115869005B (zh) | 2021-09-26 | 2021-09-26 | 超声图像的干扰抑制方法、超声设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115869005A CN115869005A (zh) | 2023-03-31 |
CN115869005B true CN115869005B (zh) | 2024-06-04 |
Family
ID=85762656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111131090.4A Active CN115869005B (zh) | 2021-09-26 | 2021-09-26 | 超声图像的干扰抑制方法、超声设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115869005B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0584510A2 (de) * | 1992-08-25 | 1994-03-02 | Leuze electronic GmbH + Co. | Aus einem Sender und einem Empfänger bestehende Einrichtung zum Erfassen von Gegenständen |
CN111374701A (zh) * | 2018-12-29 | 2020-07-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声诊断设备及其降低超声回波数据干扰的方法 |
CN111714155A (zh) * | 2020-06-30 | 2020-09-29 | 深圳开立生物医疗科技股份有限公司 | 一种超声抗干扰***、方法及超声设备 |
CN111915515A (zh) * | 2020-07-21 | 2020-11-10 | 无锡声亚医疗科技有限公司 | 超声图像中噪声去除的方法、超声设备及存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6761691B2 (en) * | 2000-07-21 | 2004-07-13 | Fuji Photo Film Co., Ltd. | Image forming method used in ultrasonic diagnosis, ultrasonic diagnostic apparatus, signal processing apparatus, and recording medium for recording signal processing program |
US9172476B2 (en) * | 2012-03-09 | 2015-10-27 | The United States Of America As Represented By The Secretary Of The Army | Method and system for removal of noise in signal |
JP6758864B2 (ja) * | 2016-03-03 | 2020-09-23 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置及び医用画像処理方法 |
US20180028153A1 (en) * | 2016-07-26 | 2018-02-01 | Toshiba Medical Systems Corporation | Ultrasound diagnostic apparatus and ultrasound imaging method |
-
2021
- 2021-09-26 CN CN202111131090.4A patent/CN115869005B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0584510A2 (de) * | 1992-08-25 | 1994-03-02 | Leuze electronic GmbH + Co. | Aus einem Sender und einem Empfänger bestehende Einrichtung zum Erfassen von Gegenständen |
CN111374701A (zh) * | 2018-12-29 | 2020-07-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声诊断设备及其降低超声回波数据干扰的方法 |
CN111714155A (zh) * | 2020-06-30 | 2020-09-29 | 深圳开立生物医疗科技股份有限公司 | 一种超声抗干扰***、方法及超声设备 |
CN111915515A (zh) * | 2020-07-21 | 2020-11-10 | 无锡声亚医疗科技有限公司 | 超声图像中噪声去除的方法、超声设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN115869005A (zh) | 2023-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11635514B2 (en) | Imaging methods and apparatuses for performing shear wave elastography imaging | |
US6368277B1 (en) | Dynamic measurement of parameters within a sequence of images | |
CA2340246C (en) | Adaptive cancellation of ring-down artifact in ivus imaging | |
JP5797364B1 (ja) | 超音波観測装置、超音波観測装置の作動方法、及び超音波観測装置の作動プログラム | |
US8864671B2 (en) | Methods and systems for color flow imaging | |
KR102286299B1 (ko) | 초음파 이미지 디스플레이 장치 및 초음파 이미지를 디스플레이하기 위한 방법 | |
US11071525B2 (en) | Ultrasonic diagnostic apparatus and method | |
US10456115B2 (en) | Ultrasound system and clutter filtering method thereof | |
US20120136248A1 (en) | Ultrasound diagnosis apparatus, image generating method, and image processing apparatus | |
CN113040878B (zh) | 超声穿刺针的位置信息处理方法、超声设备及存储介质 | |
JP2017038837A (ja) | 超音波画像診断装置、超音波画像処理方法及び超音波画像処理プログラム | |
CN113379664B (zh) | 超声图像中超声穿刺针的增强方法、超声设备及存储介质 | |
JP6997508B2 (ja) | 超音波診断装置及びその制御プログラム | |
WO2012090658A1 (ja) | 超音波診断装置および画像処理方法 | |
US20110190632A1 (en) | Ultrasonic diagnostic appratus and ultrasonic image processing method | |
CN115869005B (zh) | 超声图像的干扰抑制方法、超声设备和存储介质 | |
CN111434311B (zh) | 超声波摄像装置以及图像处理方法 | |
JP2011000470A (ja) | 超音波診断装置 | |
CN110604595B (zh) | 一种脂肪肝定量分析方法及脂肪肝定量分析*** | |
CN115666400A (zh) | 辅助用户执行医学超声检查 | |
CN113545806A (zh) | ***弹性成像方法和超声弹性成像*** | |
JP2019010189A (ja) | 超音波装置及びその制御プログラム | |
JP2019072438A (ja) | 超音波診断装置及びその制御プログラム | |
WO2021232192A1 (zh) | 超声造影成像方法、装置和存储介质 | |
CN112702955A (zh) | 一种检查模式的确定方法及超声设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |