CN104160423A - 用于抑制超声图像中的非稳定混响的滤波系统和方法 - Google Patents
用于抑制超声图像中的非稳定混响的滤波系统和方法 Download PDFInfo
- Publication number
- CN104160423A CN104160423A CN201380012973.0A CN201380012973A CN104160423A CN 104160423 A CN104160423 A CN 104160423A CN 201380012973 A CN201380012973 A CN 201380012973A CN 104160423 A CN104160423 A CN 104160423A
- Authority
- CN
- China
- Prior art keywords
- signal
- filtering
- group
- ultrasonoscopy
- data
- 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.)
- Pending
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 100
- 238000002604 ultrasonography Methods 0.000 title claims abstract description 42
- 230000001629 suppression Effects 0.000 title abstract 2
- 238000003384 imaging method Methods 0.000 claims description 82
- 230000004044 response Effects 0.000 claims description 18
- 238000013519 translation Methods 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 10
- 230000008520 organization Effects 0.000 claims description 10
- 230000001373 regressive effect Effects 0.000 claims description 10
- 230000007246 mechanism Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 43
- 238000010586 diagram Methods 0.000 description 26
- 230000008569 process Effects 0.000 description 12
- 230000033001 locomotion Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000007788 liquid Substances 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 239000000919 ceramic Substances 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 2
- 210000001124 body fluid Anatomy 0.000 description 2
- 239000010839 body fluid Substances 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 238000002386 leaching Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一般来说,本技术针对用于抑制超声图像中的混响伪像的滤波系统和方法。在一些实施例中,一种获取经滤波的超声图像的方法包括使用声头摄取目标组织的第一超声图像。移动声头的至少一部分,以使得混响伪像超声路径长度相对声头的第一位置变化。然后,摄取目标组织的第二超声图像。使用至少一种滤波方法合成第一和第二超声图像。所述滤波方法会减弱或去除合成的超声图像中的混响伪像。
Description
相关申请的交叉引用
本申请要求2012年2月2日提交的美国临时专利申请No.61/594,244的优先权益,通过引用将其完整地结合到本文中。
技术领域
一般来说,本技术针对用于抑制超声图像中的非稳定混响的滤波系统和方法。
背景技术
常规医疗超声成像系统从换能器元件阵列传送高频声音的短脉冲,然后在该声音传播过组织并反射回换能器阵列后接收它。在知道组织内声音速度的情况下(通常假设为恒定1540m/s),可基于超声的传送和接收之间的时间延迟创建组织的图像-假设声音直接从换能器阵列传播到反射组织并且直接回到换能器阵列。但是,如果声音多次反射,时间延迟将不再对应于超声反射器的物理位置,并且结构可能出现在超声图像中不同于其真实物理位置的位置。这些混响伪像可使超声图像中的组织特征模糊不清,并且可不利地影响超声图像的质量。伪像可以是“非稳定的”,因为混响伪像的表象位置会相对被成像组织/物质移动。
非稳定混响伪像可能出现于,例如,具有通过液体耦合到患者的独立成像和治疗换能器的超声治疗设备中,其中成像换能器的位置相对患者固定,而治疗换能器相对患者移动以改变治疗换能器的焦点以处理组织体。在这种系统中,由于在相同设备内稳定部件的超声波多次反射的结果,可能还存在稳定伪像。在这种系统中,来自治疗换能器或其它表面的二次反射可能在组织图像中产生混响伪像,由此使组织细节模糊不清。相应地,存在从超声图像中去除混响伪像的需要。
附图说明
图1是按照本技术的实施例配置的用于对非稳定混响伪像进行滤波的系统的示意图。
图2A是按照本技术的实施例采集第一位置的超声图像的成像换能器的示意图。
图2B是按照本技术的实施例的来自图2A的超声采集的混响伪像的超声图像的部分示意图。
图2C是按照本技术的实施例采集第二位置的超声图像的成像换能器的示意图。
图2D是按照本技术的实施例的来自图2C的超声采集的混响伪像的超声图像的部分示意图。
图2E是按照本技术的实施例滤出混响伪像之后从图2B和2D的超声图像合成的超声图像的部分示意图。
图3是示出按照本技术的实施例采集和滤波超声数据以去除非稳定混响伪像的方法的框图。
图4A是按照本技术的实施例进行滤波之前的超声图像的图示。
图4B是按照本技术的实施例采集图4A的超声图像期间的超声信号的速度分布的图示。图4C是按照本技术的实施例基于图4B的速度分布对图像进行滤波之后的图4A的超声图像的图示。
图5是示出按照本技术的实施例采集和滤波超声数据以去除非稳定混响伪像的另一种方法的框图。
图6是示出按照本技术的实施例采集和滤波超声数据以去除非稳定混响伪像的另一种方法的框图。
图7是示出按照本技术的实施例采集和滤波超声数据以去除非稳定混响伪像的另一种方法的框图。
具体实施方式
一般来说,本技术针对用于抑制超声图像中的混响伪像的滤波系统和方法。在一些实施例中,一种获取经滤波的超声图像的方法包括使用声头摄取目标组织的第一超声图像。移动声头的至少一部分,以使得混响伪像超声路径长度相对声头的第一位置变化。然后,摄取目标组织的第二超声图像。使用至少一种滤波方法合成第一和第二超声图像。所述滤波方法会减弱或去除合成的超声图像中的混响伪像。
下面将参照图1-7描述本技术的若干实施例的具体细节。以下公开没有阐述与超声成像相关的众所周知的结构和系统的其它细节,以避免不必要地使本技术的各种实施例的描述模糊不清。图中显示的许多细节、尺寸、角度和其它特征只是本技术的特定实施例的说明。相应地,只要没有背离本技术的精神或范围,其它实施例可具有其它细节、尺寸、角度和特征。因此,本领域的技术人员将相应地理解,本技术可能有包含其它元件的其它实施例,或者本技术可能有不包含以下参照图1-7示出或描述的若干特征的其它实施例。
图1是按照本技术的实施例配置的用于对非稳定混响伪像进行滤波的系统100的示意图。系统100可包括配置成摄取目标组织的超声图像的声头114。如下文更详细地描述,声头114可包含置于具有患者接口盖130的壳体140之内或邻近的成像换能器150和高强度聚焦超声(HIFU)换能器112。声头114还可包含配置成相对目标组织移动成像换能器150或HIFU换能器112中的至少一个的移动机构110。
HIFU来源或换能器112可配置成发射高能级超声能量。在若干实施例中,系统100包含共同或独立操作的多个换能器112。HIFU换能器112可从HIFU控制器118接收有关与HIFU处理相关的持续时间、频率、幅度和/或其它因素的指令。更具体地说,HIFU控制器118可与能够生成电波形以驱动HIFU换能器112的HIFU发射器122通信。
成像控制器102可控制系统100的操作。例如,成像控制器102可与HIFU控制器118通信(如发送和接收指令)。成像控制器102还可与能够控制或配置HIFU换能器112输出方面的波束形成器104通信和向其发送指令。例如,波束形成器104可设置HIFU换能器112中的波形延迟以引导超声波束。波束形成器104可与可包含脉冲和/或波形驱动电路的成像发射器106通信。例如,成像发射器106可将所需的波形驱动到成像换能器150。成像换能器150可从发射/接收开关108接收作为电信号的波形指令。在一些实施例中,发射/接收开关108可包含用于阻挡高幅度发射信号的级。成像换能器150可将来自发射/接收开关108的电信号转换成可由HIFU换能器112使用的声信号。
如上所述,HIFU换能器112可通过一个或多个移动机构110移动或调整。例如,在一些实施例中,移动机构110包含导向摇摆以沿角度θ调整HIFU换能器112的位置。在其它实施例中,移动机构110可包含高度调整以沿z轴修改HIFU换能器112和成像换能器150之一或两者的高度。在又一些实施例中,HIFU换能器112、成像换能器150或其它部件可被转动、平移、调整或以其它方式由其它机构或在其它方向或角度移动。移动机构110可从运动控制器116接收移动和/或定位指令,以将HIFU换能器112和/或成像换能器150调整到所需的位置。如下文中参照图2A-7更详细地讨论,在若干实施例中,运动控制器116可用于将HIFU换能器112和/或成像换能器150移动到序列的、偏移的位置,以获取迭代数据用作滤波过程的一部分以减弱非稳定混响伪像。在若干实施例中,运动控制器116可与HIFU控制器118通信,或可包含与HIFU控制器118相同的设备。
在若干实施例中,系统100包含可定位于HIFU换能器112与组织之间的换能器液体120。患者接口盖130可包含将换能器液体120留在壳体140中并邻近目标组织的薄膜接口。在一些实施例中,壳体140还封装或包含HIFU换能器112。虽然HIFU换能器112、成像换能器150、患者接口盖130、壳体140、移动机构110和/或其它部件被统称为“声头”114,但是将会理解,在其它实施例中,术语“声头”114可表示整个超声成像组件或其单个部件或子部件。例如,在一些实施例中,术语“声头”114可表示HIFU换能器112或成像换能器150中的一个。在又一些实施例中,术语“声头”114可表示成像换能器150、HIFU换能器112、壳体140或系统100的其它部分的反射表面部分或其它特征。在其它实施例中,“声头”包括与用户(如执业医生)或患者进行接口的系统100的任何方面。
随着向HIFU换能器112或壳体140内的其它表面提供声信号,成像换能器150可从HIFU换能器112接收声信号并且将该声信号转换成要传递到发射/接收开关108并由下游处理部件处理的电信号。例如,发射/接收开关108可将电信号传递到成像接收器160,它可包含用于从换能器元件接收电信号的放大器、以及滤波器和数据转换电路。成像接收器160可将接收的信号或数据传递到射频(RF)数据处理器162。更具体地说,如上所述,接收的信号(或时间或位置偏移信号系列)可传递到RF数据处理器162以进行滤波。RF数据处理器162可对来自RF数据的混响和/或其它伪像进行滤波。如以下将更详细地讨论,滤波可包括在图像渲染之前减弱非稳定混响信号的任何过程。RF数据处理器162可将数据传递到波束形成器164,它可执行RF信道数据处理以引导接收的数据。然后,可将该数据传递到另一RF数据处理器166以对数字信号进行滤波和/或处理,以便去除或减少混响和/或其它伪像。在其它实施例中,系统100可包含波束形成器164之前或之后的仅一个RF数据处理器,而在又一些实施例中,在系统100的这些或其它位置可能有超过两个RF数据处理器。然后,经滤波的数据可传递到图像形成器168以从波束形成的数据生成超声图像。将该图像发送到图像数据处理器170,它配置成对来自图像数据的混响和/或其它伪像进行滤波。在一些实施例中,成像接收器160、RF数据处理器162和166、波束形成器164、图像形成器168以及图像数据处理器170中的一个或多个与作为整个系统控制的组成部分的成像控制器102通信、从其接收指令和/或向其提供反馈。最后,经滤波的图像显示在输出设备或监视器上,例如图像显示面板172。将会理解,系统100可包含不同的特征和/或具有不同的布置。例如,在其它实施例中,可去除上述的系统100的一个或多个部件。在又一些实施例中,系统100的单独描述的部件的功能性可结合到单个设备或子部件中。
术语“滤波”在本文中被广泛使用,并可包括任何适用的滤波方法或技术。各种滤波方法可在信号处理通路中的各个点单独地或结合地实现。可能的滤波器包括但不限于有限脉冲响应(FIR)滤波器、无限脉冲响应(IIR)滤波器、回归滤波器及百分比滤波器。滤波可在波束形成之前(如在RF数据处理器162)、波束形成之后(如在RF数据处理器166)、或在图像形成之后在图像数据处理器170逐个信道完成。滤波器选择可由实现滤波方法所在的超声系统的硬件和软件能力来决定。下面提供关于若干适用滤波器的进一步的细节。
FIR/IIR滤波器-FIR和IIR滤波器可用于超声系统中,特别是用于多普勒血流成像以抑制来自组织的信号。在一些实施例中,低通或带阻滤波器可与一组固定的预先计算的滤波器系数一起使用。在另一实施例中,可基于混响伪像的速度(或表示运动的其它衡量标准)自适应地计算滤波器系数。在这个实施例中,会使用多普勒方法(或用于估计运动的其它已知超声方法)来估计速度,并且从这个速度估计,可自适应地调整滤波器系数以最优地抑制混响。在另一实施例中,多组滤波器系数可被预先计算,并且在脉冲传送/接收期间基于混响来源的速度的先验知识进行互换。
回归滤波器-回归滤波器可用于例如多普勒超声处理中,特别是用于彩色多普勒和能量多普勒。对于彩色多普勒和能量多普勒,超声数据可拟合到低阶多项式函数并且余数可用于血流计算。对于混响抑制,余数将由混响信号支配并可被丢弃,并且拟合将由组织信号支配。
平均滤波器-平均滤波器会对来自多个脉冲的信号求平均。这些滤波器是FIR滤波器的概念变体。平均可以是加权的或未加权的。在一些实施例中,在最强的信号是混响的操作假设下,每个脉冲可相对其幅度成反比加权。
百分比滤波器-百分比滤波器按幅度分类信号并保留处于、高于或低于某个阈值的信号。在一个实施例中,使用最小值滤波器,由此只在逐个像素的基础上保留具有最小幅度的脉冲。
将会理解,本文中描述的系统100和其它系统/方法可使用上述滤波器中的一种或多种或其它适用滤波器。
图2A是按照本技术的实施例采集第一位置的超声图像的成像换能器的示意图。图2B是按照本技术实施例的来自图2A的超声采集的混响伪像的超声图像的部分示意图。一起参照图2A和2B,成像换能器150定位于至少部分在壳体140中的深度z1。成像换能器150通过换能器液体120和患者接口盖130沿路径L1向患者组织200发送超声波。组织200包括组织特征202,它是组织200中产生超声图像回声的特征或者体。如图2B所示,超声图像204配置成向用户如执业医生显示。超声图像204包含伪像206,它是混响的声音波形的结果,该混响的声音波形具有在以d1表示的等效深度处显示的混响信号通路长度L3+L4+L5+L6。
图2C和2D示出相似的布置。例如,图2C是按照本技术的实施例采集第二位置的超声图像的成像换能器的示意图。但是,在这个布置中,成像换能器150定位于深度z2。图2D是按照本技术的实施例的来自图2C的超声采集的混响伪像的超声图像的部分示意图。如图2D所示,超声图像204包含伪像206,它是混响的声音波形的结果,该混响的声音波形具有在以d2表示的等效深度处显示的路径长度L3+L4+L5+L6。
在换能器液体120内出现声音波形的两个或更多反射时出现混响伪像206,并且全部或部分反射的声音波形由成像换能器150接收。常规成像换能器在陶瓷换能器元件与组织接口之间的物质内出现混响伪像条件,虽然在感兴趣成像深度之前能量通常减弱到足够的水平。成像换能器150在接收从组织反射的波形的同时接收混响的波形。这两个接收的波形由超声系统处理并显示在超声图像204上,以使得混响伪像206显示在距离换能器等效于一半混响路径长度((L2+L3+L4+L5)/2)的位置。虽然所示混响路径(L2+L3+L4+L5)属于简单的两个表面反射,声音波形由成像换能器150接收之前可从换能器液体120内的许多表面反射。在图像数据被接收并处理成超声图像时,成像换能器150的表面在距离成像换能器150等同于混响路径长度(d1或d2)的距离处显示为混响伪像206。
在成像换能器150或反射表面(如HIFU换能器112)移位以导致混响路径长度(L2+L3+L4+L5)改变时,以上参照图1描述的系统100可促使混响路径长度的移位改变。例如,在成像换能器150从位置z1移位到位置z2时,所产生的混响伪像206出现在各自的深度d1和d2。例如,图2E是滤出混响伪像之后从图2B和2D的超声图像合成的超声图像的部分示意图。如图2E所示,通过结合分别在z1和z2捕获的两个数据组,可产生其中不存在混响伪影206的第三数据组。虽然本文所述的技术可能表示处理两个数据组,但相应算法也可应用于超过两个数据组。在一些实施例中,换能器位置和各自数据组的数量可介于4到8的范围内。在其它实施例中,在使用复合成像技术时混响路径长度可改变,其中声音波形以不同的引导角度从陶瓷换能器阵列的不同陶瓷元件组传送。
虽然对于本文所述的一些技术来说,可能希望混响伪像206停留在深度d2,其中伪像206不与位于深度d1的伪像206重叠,但其它技术,例如促使伪像206的“速度”改变的那些技术不需要d2-d1的分离范围。在这些情况下,由于能量在时间和空间上的重叠,所接收能量的位置变化可能较小,虽然速度仍可被促使改变。然后,可应用当前已知的信号处理技术和算法,以滤出后续的或一系列的数据组之间的“移动”能量。
还应该注意到,虽然所述算法指示在多个位置或传送角度捕获数据和处理单个图像的多个数据组的一系列过程,但在其它实施例中,可采用使所显示的图像数量等于所采集的数据集数量的方式处理数据(如使用滑动窗口和并行处理)。
可采用各种方法来采集超声数据和对数据进行滤波以减弱非稳定混响伪像。在图3的框图所示的第一方法300中,由固定或可变时间间隔分隔的重复超声脉冲序列被获取并定义随后在图像形成之前滤波的脉冲的集合体,以减弱可归因于非稳定混响伪像的信号部分,例如移动的HIFU换能器。在框302,方法300包括在时间t(1)定位成像换能器,并将计数器k设置为等于2。在框304,在时间t(1)采集超声图像。在框306,方法300包括等待至时间t(k)。在框308采集后续超声图像k。在框310和312对迭代超声图像重复这个过程,直至过程到达最终时间t(K)。如上所述,超声脉冲之间的时间可以是固定的、可变的或随机的。
在框316,将计数器设置为指向图像数据中的第一个像素。在框318,在从图像1到图像K的像素n处生成像素数据矢量。在框320将滤波器(如以上参照图1提及的任何滤波器)应用于所述矢量。滤波器可用于从超出预定范围的速度移动的超声图像中去除伪像,例如比组织或体液移动得更快的混响伪像。在其它实施例中,可设置滤波器以使得移动组织和液体随混响伪像滤出;这种滤波器可能很有用,例如,如果图像的功能是表达稳定组织边界。在特定实施例中,滤波器是低通滤波器,如回归滤波器、FIR滤波器或IIR滤波器。备选地,如果混响伪像的表观速度是先验已知的,可使用带阻滤波器更有选择性地对其进行滤波。
在框322,滤波结果存储在图像1的像素n中。在框324和326,重复此滤波过程,直至到达图像文件的末尾。在一些实施例中,脉冲之间的时间间隔、集合体长度和/或滤波器阶数和系数可手动或自适应地改变,以提高滤波器性能。在对图像进行滤波后,在框328显示图像1。将会理解,在其它实施例中,方法300可包括另外的步骤和/或可去除上述一个或多个步骤。
图4A是按照本技术的实施例进行滤波之前的超声图像的图示。图4B是按照本技术的实施例采集图4A的超声图像期间的超声信号的速度分布的图示。图4C是按照以上参照图3所述的方法基于图4B的速度分布对图像进行滤波之后的图4A的超声图像的图示。在特定实施例中,可使用具有16个脉冲的集合体的回归滤波器。所产生的、经滤波的超声图像显示高速混响伪像的大大减少,提供更清晰的目标组织图像。
参照图3所述的方法300的有效性取决于组织速度和混响表观速度的差异。图5是示出包括平移成像换能器以创建或增强组织和混响的速度差异的方法500的框图。随着成像换能器被平移,组织的表观位置将被平移相同的量,并且更快移动的混响的表观位置将被平移更大的量。因此,使用此方法500,可更容易地识别稳定和非稳定伪像,以便减弱和去除。
例如,在框502,方法500包括将成像换能器定位于位置z(1),并将计数器k设置为等于2。在框504,在位置z(1)采集超声图像。在框506,方法500包括将成像换能器移到位置z(k)。在框508,在位置z(k)采集后续超声图像k。在框510和512对迭代超声图像重复这个过程,直至过程到达最终位置z(K)。在一些实施例中,成像换能器的位置由计算机控制,而在其它实施例中,成像换能器通过手动控制进行定位。
在框514,调准图像1和k,并且在框516将计数器n设置为指向图像数据中的第一个像素。在框518,在从图像1到图像K的像素n处生成像素数据矢量。在框520将滤波器(如以上参照图1提及的任何滤波器)应用于矢量。滤波器可用于从在两个或更多图像帧中落在两个或更多不同和独立位置的超声图像中去除伪像。备选地,滤波器可用于从超出预定范围的速度移动的超声图像中去除伪像,例如比组织或体液移动得更快的混响伪像。混响的速度可放大(并更容易滤出),因为成像换能器在多个位置采集数据。在框522,滤波结果存储在图像1的像素n中。在框524和526,重复此滤波过程,直至到达图像文件的末尾。在框528,显示经滤波的图像。
图6是示出采集超声数据并对其进行滤波以去除非稳定混响伪像的另一种方法600的框图。方法600考虑了混响对象的表观位置可取决于声音传播的方向。方法600使用在不同方向发射的脉冲序列。然后,对脉冲进行滤波以抑制并非所有成像方向共有的特征。为了减少由目标组织和/或HIFU系统的部件移动造成的伪像,可为特定聚焦区域完成序列脉冲。
例如,在框602,方法600包括编程波束形成器以角度θ(1)发射超声波形。在框604,发射波形并将计数器k设置为2;在框606采集超声图像1。迭代地,在框608-616,迭代地采集后续超声图像,直至最终发射角度θ(K)被采集。
在框618,调准图像1和k,并且在框620将计数器n设置为指向图像数据中的第一个像素。在框622,在从图像1到图像K的像素n处生成像素数据矢量。在框624将滤波器(如最小值滤波器、加权平均滤波器或其它滤波器)应用于所述矢量。除减弱混响对象之外,这些多个引导角度的空间复合可用于减少斑点。对各种引导脉冲的混响对象反射将具有比组织更高的差异。此差异可用于划分出混响对象。在框626,滤波结果存储在图像1的像素n中。在框628和630,重复此滤波过程,直至到达图像文件的末尾。在框632,显示经滤波的图像。
图7是示出采集超声数据并对其进行滤波以去除非稳定混响伪像的又一种方法700的框图。如图7所示,可对图像序列应用另外的滤波以进一步抑制混响。在单张图像中,可抑制混响的程度将随生成混响伪像的结构的速度而改变。如果图像之间的抑制程度差异较大,可对图像序列进行滤波或将其“缝合”到一起以进行另外的混响抑制。可能的滤波器包括但不限于应用于图像之间的均值、中值和最小值滤波器。
例如,在框702,方法700包括将成像换能器定位于位置z(1),并将计数器k设置为等于2。在框704,在位置z(1)采集超声图像。在框706,方法700包括将成像换能器移到位置z(k)。在框708,在位置z(k)采集后续超声图像k。在框710和712对迭代超声图像重复这个过程,直至换能器到达最终位置z(K)。在一些实施例中,成像换能器的位置由计算机控制,而在其它实施例中,成像换能器通过手动控制进行定位。
在框714,调准图像1和k,并且在框716将计数器n设置为指向图像数据中的第一个像素。在框718,在从图像1到图像K的像素n处生成像素数据矢量。在框720,对比图像像素亮度的差异,并应用滤波技术。在一些实施例中,如框722和724所示,如果差异大于预定阈值,则在图像1上保存最低(即最暗)的像素值。在框726,如果差异在阈值内,则在图像1上保存像素数据的平均值。在框728和730,对图像文件中的所有像素重复这个过程。在一个备选实施例中,只有[p(1,n),p(2,n),…,p(K,n)]中的最小值会存储在图像1的像素n中。在又一个实施例中,只有[p(1,n),p(2,n),…,p(K,n)]的平均值会存储在图像1的像素n中。在又一个实施例中,通过移动HIFU换能器100改变混响路径长度。在框732,显示经滤波的图像。
示例
以下示例说明本技术的若干实施例。
1.一种用于获取超声图像的方法,所述方法包括:
经由声头摄取目标组织的第一超声图像;
移动所述声头的至少一部分;
经由所述声头摄取所述目标组织的第二超声图像;
合成所述第一超声图像与所述第二超声图像;以及
消除合成的超声图像中的混响伪像。
2.如示例1所述的方法,其中,移动所述声头的至少一部分包括移动所述声头的反射表面。
3.如示例1所述的方法,其中,移动所述声头的至少一部分包括移动成像换能器。
4.如示例1所述的方法,其中,移动所述声头的至少一部分包括移动高强度聚焦的超声换能器。
5.如示例1所述的方法,其中,移动所述声头的至少一部分包括平移或调整所述声头的所述部分的角度其中的至少一个。
6.如示例1所述的方法,其中,消除所述混响伪像包括对所述第一超声图像和所述第二超声图像使用低通滤波器。
7.如示例1所述的方法,其中,消除所述混响伪像包括使用回归滤波器、有限脉冲响应滤波器、无限脉冲响应滤波器或带阻滤波器中的至少一个。
8.如示例1所述的方法,其中,消除所述混响伪像包括使用均值、中值或最小值滤波器中的至少一个。
9.一种抑制超声图像中的非稳定混响伪像的方法,所述方法包括:
通过成像换能器将多个超声脉冲传送到组织体,其中多个超声脉冲的传送以时间间隔分隔;
响应传送多个超声脉冲的每一个,从所述组织体接收反射信号,其中,所述反射信号包括并非源自所述组织体的混响信号;
对所述反射信号进行滤波以减弱可归因于所述混响信号的反射信号部分;以及
基于所述滤波信号生成所述组织体的超声图像。
10.如示例9所述的方法,其中,传送多个超声脉冲包括传送以固定或变化时间间隔分隔的脉冲。
11.如示例9所述的方法,其中,对所述反射信号进行滤波包括使用低通滤波器对所述信号进行滤波,并且其中所述低通滤波器是回归滤波器、有限脉冲响应滤波器、无限脉冲响应滤波器以及带阻滤波器中的至少一个。
12.如示例9所述的方法,还包括在多个超声脉冲的传送之间将所述成像换能器平移固定的量。
13.如示例9所述的方法,还包括:
通过所述成像换能器传送第二多个超声脉冲到所述组织体;
响应传送第二多个超声脉冲的每一个,从所述组织体接收第二反射信号,其中所述第二反射信号包括并非源自所述组织体的混响信号;
对所述第二反射信号进行滤波以减弱可归因于所述混响信号的第二反射信号的部分;
基于所述第二滤波信号生成所述组织体的第二超声图像;以及
对比所述第一和所述第二超声图像以进一步减弱可归因于所述混响信号的图像部分。
14.如示例13所述的方法,其中,对比所述第一和所述第二超声图像包括使用均值、中值和最小值滤波器中的至少一个对所述图像进行滤波。
15.一种减弱超声图像中的非稳定伪像的方法,所述方法包括:
通过在第一信号通路上应用第一信号来捕获包含伪像的第一组超声图像数据;
通过在第二信号通路上应用第二信号来捕获第二组超声图像数据;以及
结合所述第一组数据和所述第二组数据以消除伪像。
16.如示例15所述的方法,其中,结合所述第一组数据和所述第二组数据包括使用低通滤波器对所述第一组数据和所述第二组数据进行滤波。
17.如示例15所述的方法,其中,捕获所述第一组超声图像数据在第一时间间隔发生,并且捕获所述第二组超声图像数据在与所述第一时间间隔隔开的第二时间间隔发生。
18.如示例15所述的方法,其中,在第二信号通路上应用第二信号包括在具有与所述第一信号通路不同的通路长度或通路路线的第二信号通路上应用第二信号。
19.如示例15所述的方法,其中,结合所述第一组数据和所述第二组数据以消除伪像包括对所述第一组数据和所述第二组数据应用均值、中值或最小值滤波器中的至少一个。
20.如示例19所述的方法,还包括至少部分基于所述第一组数据与所述第二组数据之间的差异选择均值、中值或最小值滤波器中的至少一个。
21.一种其上存储了计算机可执行指令的物理计算机可读存储介质,如果所述指令在由计算系统执行时,使所述计算系统执行包括以下的操作:
指示成像换能器在第一信号通路上应用第一信号;
基于所述第一信号接收包含伪像的第一组超声图像数据;
指示所述成像换能器在第二信号通路上应用第二信号;
基于所述第二信号接收第二组超声图像数据;以及
结合所述第一组数据和所述第二组数据以消除伪像。
22.如示例21所述的物理计算机可读存储介质,其中,所述操作还包括使用均值、中值或最小值滤波器中的至少一个对所述第一组超声图像数据和所述第二组超声图像数据进行滤波。
23.如示例21所述的物理计算机可读存储介质,其中,所述操作还包括使用低通滤波器对所述第一组数据和所述第二组数据进行滤波。
24.如示例21所述的物理计算机可读存储介质,其中,指示所述成像换能器在第二信号通路上应用第二信号包括指示所述成像换能器在具有与所述第一信号通路不同的通路长度或通路路线的第二信号通路上应用第二信号。
25.如示例21所述的物理计算机可读存储介质,其中,所述操作还包括在应用所述第一信号和应用所述第二信号之间将所述成像换能器平移固定的量。
26.一种超声成像系统,包括:
成像换能器,配置成向组织体传递超声脉冲;
接收器,配置成接收来自所述组织体的反射信号,其中所述反射信号包括并非源自所述组织体的混响信号;以及
控制器,编程为对所述反射信号进行滤波以减弱可归因于所述混响信号的反射信号部分,并基于所述滤波信号生成所述组织体的超声图像。
27.如示例26所述的系统,其中,所述成像换能器配置成传递以固定或变化时间间隔分隔的超声脉冲。
28.如示例26所述的系统,其中,编程为对所述反射信号进行滤波的所述控制器包括编程为使用低通滤波器对所述反射信号进行滤波的控制器,并且其中所述低通滤波器包括回归滤波器、有限脉冲响应滤波器、无限脉冲响应滤波器以及带阻滤波器中的至少一个。
29.如示例26所述的系统,还包括配置成移动所述成像换能器并由此修改所述反射信号的长度的移动机构。
30.如示例26所述的系统,其中,编程为对所述反射信号进行滤波的所述控制器包括编程为使用均值、中值或最小值滤波器中的至少一个对所述反射信号进行滤波的控制器。
从以上所述将会理解,虽然为了说明的目的在本文中描述了本技术的具体实施例,但可进行各种修改,而不背离本技术的精神和范围。例如,本文所述的系统和方法可用于对非稳定混响伪像之外的伪像进行滤波。更具体地说,作为一种去除稳定伪像的手段,通过改变成像换能器相对患者组织接口的距离,非稳定混响伪像技术也可用于只有稳定伪像的设备中。此外,在特定实施例的上下文中描述的新技术的某些方面在其它实施例中可结合或消除。而且,虽然在这些实施例的上下文中描述了与本技术的某些实施例相关的优点,其它实施例也可展现此类优点,并且不是所有实施例都必须展现此类优点才落入本技术的范围。相应地,本公开和相关技术可包含本文中未明确示出或描述的其它实施例。因此,本公开仅受所附权利要求书限制。
Claims (30)
1.一种用于获取超声图像的方法,所述方法包括:
经由声头摄取目标组织的第一超声图像;
移动所述声头的至少一部分;
经由所述声头摄取所述目标组织的第二超声图像;
合成所述第一超声图像与所述第二超声图像;以及
消除合成的超声图像中的混响伪像。
2.如权利要求1所述的方法,其中,移动所述声头的至少一部分包括移动所述声头的反射表面。
3.如权利要求1所述的方法,其中,移动所述声头的至少一部分包括移动成像换能器。
4.如权利要求1所述的方法,其中,移动所述声头的至少一部分包括移动高强度聚焦的超声换能器。
5.如权利要求1所述的方法,其中,移动所述声头的至少一部分包括平移或调整所述声头的所述部分的角度其中的至少一个。
6.如权利要求1所述的方法,其中,消除所述混响伪像包括对所述第一超声图像和所述第二超声图像使用低通滤波器。
7.如权利要求1所述的方法,其中,消除所述混响伪像包括使用回归滤波器、有限脉冲响应滤波器、无限脉冲响应滤波器或带阻滤波器中的至少一个。
8.如权利要求1所述的方法,其中,消除所述混响伪像包括使用均值、中值或最小值滤波器中的至少一个。
9.一种抑制超声图像中的非稳定混响伪像的方法,所述方法包括:
通过成像换能器将多个超声脉冲传送到组织体,其中,多个超声脉冲的传送以时间间隔分隔;
响应传送多个超声脉冲的每一个,从所述组织体接收反射信号,其中,所述反射信号包括并非源自所述组织体的混响信号;
对所述反射信号进行滤波以减弱可归因于所述混响信号的反射信号部分;以及
基于所述滤波信号生成所述组织体的超声图像。
10.如权利要求9所述的方法,其中,传送多个超声脉冲包括传送以固定或变化时间间隔分隔的脉冲。
11.如权利要求9所述的方法,其中,对所述反射信号进行滤波包括使用低通滤波器对所述信号进行滤波,并且其中所述低通滤波器是回归滤波器、有限脉冲响应滤波器、无限脉冲响应滤波器以及带阻滤波器中的至少一个。
12.如权利要求9所述的方法,还包括在多个超声脉冲的传送之间将所述成像换能器平移固定的量。
13.如权利要求9所述的方法,还包括:
通过所述成像换能器传送第二多个超声脉冲到所述组织体;
响应传送第二多个超声脉冲的每一个,从所述组织体接收第二反射信号,其中,所述第二反射信号包括并非源自所述组织体的混响信号;
对所述第二反射信号进行滤波以减弱可归因于所述混响信号的第二反射信号部分;
基于所述第二滤波信号生成所述组织体的第二超声图像;以及
对比所述第一和所述第二超声图像以进一步减弱可归因于所述混响信号的图像部分。
14.如权利要求13所述的方法,其中,对比所述第一和所述第二超声图像包括使用均值、中值和最小值滤波器中的至少一个对所述图像进行滤波。
15.一种减弱超声图像中的非稳定伪像的方法,所述方法包括:
通过在第一信号通路上应用第一信号来捕获包含伪像的第一组超声图像数据;
通过在第二信号通路上应用第二信号来捕获第二组超声图像数据;以及
结合所述第一组数据和所述第二组数据以消除所述伪像。
16.如权利要求15所述的方法,其中,结合所述第一组数据和所述第二组数据包括使用低通滤波器对所述第一组数据和所述第二组数据进行滤波。
17.如权利要求15所述的方法,其中,捕获所述第一组超声图像数据在第一时间间隔发生,并且捕获所述第二组超声图像数据在与所述第一时间间隔隔开的第二时间间隔发生。
18.如权利要求15所述的方法,其中,在第二信号通路上应用第二信号包括在具有与所述第一信号通路不同的通路长度或通路路线的第二信号通路上应用第二信号。
19.如权利要求15所述的方法,其中,结合所述第一组数据和所述第二组数据以消除伪像包括对所述第一组数据和所述第二组数据应用均值、中值或最小值滤波器中的至少一个。
20.如权利要求19所述的方法,还包括至少部分基于所述第一组数据与所述第二组数据之间的差异选择均值、中值或最小值滤波器中的至少一个。
21.一种其上存储了计算机可执行指令的物理计算机可读存储介质,如果所述指令由计算系统执行时,使所述计算系统执行包括以下的操作:
指示成像换能器在第一信号通路上应用第一信号;
基于所述第一信号接收包含伪像的第一组超声图像数据;
指示所述成像换能器在第二信号通路上应用第二信号;
基于所述第二信号接收第二组超声图像数据;以及
结合所述第一组数据和所述第二组数据以消除伪像。
22.如权利要求21所述的物理计算机可读存储介质,其中,所述操作还包括使用均值、中值或最小值滤波器中的至少一个对所述第一组超声图像数据和所述第二组超声图像数据进行滤波。
23.如权利要求21所述的物理计算机可读存储介质,其中,所述操作还包括使用低通滤波器对所述第一组数据和所述第二组数据进行滤波。
24.如权利要求21所述的物理计算机可读存储介质,其中,指示所述成像换能器在第二信号通路上应用第二信号包括指示所述成像换能器在具有与所述第一信号通路不同的通路长度或通路路线的第二信号通路上应用第二信号。
25.如权利要求21所述的物理计算机可读存储介质,其中,所述操作还包括在应用所述第一信号和应用所述第二信号之间将所述成像换能器平移固定的量。
26.一种超声成像系统,包括:
成像换能器,配置成向组织体传递超声脉冲;
接收器,配置成接收来自所述组织体的反射信号,其中所述反射信号包括并非源自所述组织体的混响信号;以及
控制器,编程为对所述反射信号进行滤波以减弱可归因于所述混响信号的反射信号部分,并基于所述滤波信号生成所述组织体的超声图像。
27.如权利要求26所述的系统,其中,所述成像换能器配置成传递以固定或变化时间间隔分隔的超声脉冲。
28.如权利要求26所述的系统,其中,编程为对所述反射信号进行滤波的所述控制器包括编程为使用低通滤波器对所述反射信号进行滤波的控制器,并且其中所述低通滤波器包括回归滤波器、有限脉冲响应滤波器、无限脉冲响应滤波器以及带阻滤波器中的至少一个。
29.如权利要求26所述的系统,还包括配置成移动所述成像换能器并由此修改所述反射信号长度的移动机构。
30.如权利要求26所述的系统,其中,编程为对所述反射信号进行滤波的所述控制器包括编程为使用均值、中值或最小值滤波器中的至少一个对所述反射信号进行滤波的控制器。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261594244P | 2012-02-02 | 2012-02-02 | |
US61/594,244 | 2012-02-02 | ||
PCT/US2013/024663 WO2013116855A1 (en) | 2012-02-02 | 2013-02-04 | Filtering systems and methods for suppression of non-stationary reverberation in ultrasound images |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104160423A true CN104160423A (zh) | 2014-11-19 |
Family
ID=47750805
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380012973.0A Pending CN104160423A (zh) | 2012-02-02 | 2013-02-04 | 用于抑制超声图像中的非稳定混响的滤波系统和方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10713758B2 (zh) |
EP (1) | EP2810248A1 (zh) |
JP (1) | JP6161129B2 (zh) |
KR (1) | KR20140121859A (zh) |
CN (1) | CN104160423A (zh) |
WO (1) | WO2013116855A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110431443A (zh) * | 2017-03-13 | 2019-11-08 | 皇家飞利浦有限公司 | 用于对超声图像杂波进行滤波的方法和系统 |
CN110622033A (zh) * | 2017-05-11 | 2019-12-27 | 皇家飞利浦有限公司 | 用于控制复合超声图像的生成的方法和系统 |
CN110622034A (zh) * | 2017-05-11 | 2019-12-27 | 皇家飞利浦有限公司 | 超声诊断图像中的混响伪影消除 |
CN111226259A (zh) * | 2017-08-16 | 2020-06-02 | 皇家飞利浦有限公司 | 用于图像伪影消除的系统、方法和装置 |
CN111436971A (zh) * | 2020-04-20 | 2020-07-24 | 汕头市超声仪器研究所有限公司 | 一种超声混响伪像抑制方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112566559A (zh) * | 2018-07-11 | 2021-03-26 | 皇家飞利浦有限公司 | 具有像素外推图像增强的超声成像系统 |
JP7112588B2 (ja) * | 2019-03-19 | 2022-08-03 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法、及び超音波観測装置の作動プログラム |
US11986356B2 (en) * | 2019-11-21 | 2024-05-21 | Koninklijke Philips N.V. | Reduction of reverberation artifacts in ultrasound images and associated devices, systems, and methods |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060030776A1 (en) * | 2004-08-09 | 2006-02-09 | General Electric Company | Range dependent weighting for spatial compound imaging |
CN101305925A (zh) * | 2008-05-15 | 2008-11-19 | 上海理工大学 | 超声扫描探头 |
CN102136138A (zh) * | 2011-03-28 | 2011-07-27 | 深圳市蓝韵实业有限公司 | 一种医学图像降噪与增强的处理方法 |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3189293B2 (ja) | 1991-05-16 | 2001-07-16 | 株式会社日立製作所 | 超音波治療装置 |
JPH03114451A (ja) | 1989-09-29 | 1991-05-15 | Terumo Corp | 超音波診断装置 |
JPH03114452A (ja) | 1989-09-29 | 1991-05-15 | Terumo Corp | 超音波診断装置 |
US5465723A (en) | 1994-06-02 | 1995-11-14 | Vingmed Sound A/S | Method and apparatus for ultrasound imaging |
US5538004A (en) * | 1995-02-28 | 1996-07-23 | Hewlett-Packard Company | Method and apparatus for tissue-centered scan conversion in an ultrasound imaging system |
JPH09103434A (ja) | 1995-03-31 | 1997-04-22 | Toshiba Corp | 超音波治療装置 |
JP2969073B2 (ja) | 1995-12-28 | 1999-11-02 | アロカ株式会社 | 超音波診断装置 |
JP3732613B2 (ja) | 1997-04-15 | 2006-01-05 | フクダ電子株式会社 | 超音波診断装置 |
US6511426B1 (en) * | 1998-06-02 | 2003-01-28 | Acuson Corporation | Medical diagnostic ultrasound system and method for versatile processing |
US6210328B1 (en) * | 1998-10-01 | 2001-04-03 | Atl Ultrasound | Ultrasonic diagnostic imaging system with variable spatial compounding |
JP2003521341A (ja) | 2000-01-31 | 2003-07-15 | アー.ヤー. アンゲルセン、ビョルン | 医療用超音波イメージングにおける位相面収差およびパルス残響の補正 |
US6436041B1 (en) | 2000-12-22 | 2002-08-20 | Acuson Corporation | Medical ultrasonic imaging method with improved ultrasonic contrast agent specificity |
US6776759B2 (en) | 2002-02-27 | 2004-08-17 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for high strain rate rejection filtering |
US6905465B2 (en) | 2002-04-05 | 2005-06-14 | Angelsen Bjoern A. J. | Corrections for pulse reverberations and phasefront aberrations in ultrasound imaging |
US20050277835A1 (en) * | 2003-05-30 | 2005-12-15 | Angelsen Bjorn A | Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties |
US7311701B2 (en) * | 2003-06-10 | 2007-12-25 | Cierra, Inc. | Methods and apparatus for non-invasively treating atrial fibrillation using high intensity focused ultrasound |
US20050288587A1 (en) * | 2004-06-25 | 2005-12-29 | Yongrae Roh | Drive machanism for mechanically scanned ultrasound transducers |
CA2576068A1 (en) * | 2004-09-16 | 2006-03-23 | University Of Washington | Interference-free ultrasound imaging during hifu therapy, using software tools |
JP4615950B2 (ja) * | 2004-09-30 | 2011-01-19 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
JP2007236681A (ja) | 2006-03-09 | 2007-09-20 | Toshiba Corp | 医用画像診断装置及び医用画像診断装置の制御プログラム |
US8699765B2 (en) * | 2006-11-13 | 2014-04-15 | Siemens Medical Solutions Usa, Inc. | Reducing jittering in medical diagnostic ultrasound imaging |
JP5105966B2 (ja) | 2007-05-21 | 2012-12-26 | 一浩 飯沼 | 超音波検査装置 |
US7780601B2 (en) * | 2007-06-05 | 2010-08-24 | Siemens Medical Solutions Usa, Inc. | Adaptive clinical marker preservation in spatial compound ultrasound imaging |
CN102472814A (zh) | 2009-06-30 | 2012-05-23 | 皇家飞利浦电子股份有限公司 | 基于传播介质修改的回射信号消除 |
-
2013
- 2013-02-04 KR KR1020147023696A patent/KR20140121859A/ko not_active Application Discontinuation
- 2013-02-04 US US13/758,826 patent/US10713758B2/en active Active
- 2013-02-04 WO PCT/US2013/024663 patent/WO2013116855A1/en active Application Filing
- 2013-02-04 EP EP13706334.3A patent/EP2810248A1/en not_active Withdrawn
- 2013-02-04 JP JP2014555825A patent/JP6161129B2/ja not_active Expired - Fee Related
- 2013-02-04 CN CN201380012973.0A patent/CN104160423A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060030776A1 (en) * | 2004-08-09 | 2006-02-09 | General Electric Company | Range dependent weighting for spatial compound imaging |
CN101305925A (zh) * | 2008-05-15 | 2008-11-19 | 上海理工大学 | 超声扫描探头 |
CN102136138A (zh) * | 2011-03-28 | 2011-07-27 | 深圳市蓝韵实业有限公司 | 一种医学图像降噪与增强的处理方法 |
Non-Patent Citations (1)
Title |
---|
朱其刚等: "内窥镜超声图像伪像消除与灰度增强", 《微计算机信息》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110431443A (zh) * | 2017-03-13 | 2019-11-08 | 皇家飞利浦有限公司 | 用于对超声图像杂波进行滤波的方法和系统 |
CN110431443B (zh) * | 2017-03-13 | 2023-08-15 | 皇家飞利浦有限公司 | 用于对超声图像杂波进行滤波的方法和系统 |
CN110622033A (zh) * | 2017-05-11 | 2019-12-27 | 皇家飞利浦有限公司 | 用于控制复合超声图像的生成的方法和系统 |
CN110622034A (zh) * | 2017-05-11 | 2019-12-27 | 皇家飞利浦有限公司 | 超声诊断图像中的混响伪影消除 |
CN110622034B (zh) * | 2017-05-11 | 2023-10-20 | 皇家飞利浦有限公司 | 超声诊断图像中的混响伪影消除 |
CN110622033B (zh) * | 2017-05-11 | 2024-01-30 | 皇家飞利浦有限公司 | 用于控制复合超声图像的生成的方法和系统 |
CN111226259A (zh) * | 2017-08-16 | 2020-06-02 | 皇家飞利浦有限公司 | 用于图像伪影消除的系统、方法和装置 |
CN111226259B (zh) * | 2017-08-16 | 2024-03-12 | 皇家飞利浦有限公司 | 用于图像伪影消除的系统、方法和装置 |
CN111436971A (zh) * | 2020-04-20 | 2020-07-24 | 汕头市超声仪器研究所有限公司 | 一种超声混响伪像抑制方法 |
CN111436971B (zh) * | 2020-04-20 | 2022-06-10 | 汕头市超声仪器研究所股份有限公司 | 一种超声混响伪像抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2013116855A1 (en) | 2013-08-08 |
KR20140121859A (ko) | 2014-10-16 |
US20130204135A1 (en) | 2013-08-08 |
JP6161129B2 (ja) | 2017-07-12 |
US10713758B2 (en) | 2020-07-14 |
JP2015506258A (ja) | 2015-03-02 |
EP2810248A1 (en) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104160423A (zh) | 用于抑制超声图像中的非稳定混响的滤波系统和方法 | |
US20210085292A1 (en) | Ultrasound imaging with sparse array probes | |
KR102103137B1 (ko) | 가중 인자들을 적용함으로써 초음파 이미지 품질을 향상시키는 시스템들 및 방법들 | |
US9310475B2 (en) | Method and apparatus for transmitting multiple beams | |
JP6306012B2 (ja) | 超音波イメージングシステムのメモリアーキテクチャ | |
WO2019214127A1 (zh) | 一种经颅三维脑血管复合成像方法及系统 | |
CN110680380B (zh) | 超声波摄像装置及图像处理装置 | |
KR20160041516A (ko) | 빔포밍 장치 및 이를 포함하는 초음파 진단장치 | |
WO1995002993A1 (en) | Method and apparatus for performing imaging | |
WO1993019673A1 (en) | Method and apparatus for ultrasonic beam focusing | |
IL140280A (en) | Method and apparatus for harmonic imaging using multiple focal zones | |
EP2397867B1 (en) | Adaptive clutter filtering in an ultrasound system | |
JP2007512068A (ja) | ビーム走査反転による機械的プローブを用いた3次元超音波画像化 | |
US8235906B2 (en) | System and method for accelerated focused ultrasound imaging | |
KR101117545B1 (ko) | 평면파 중첩을 이용하여 탄성영상을 형성하는 초음파 진단 시스템 및 방법 | |
JP2007507271A (ja) | 超音波イメージングにおけるスモールアンサンブル長によるクラッタフィルタリング | |
KR20160085016A (ko) | 초음파 진단 장치 및 그 제어 방법 | |
CN106580365B (zh) | 超声设备及其控制方法 | |
KR101024851B1 (ko) | Bc-모드 영상을 형성하는 초음파 시스템 및 방법 | |
US10228462B2 (en) | Ultrasonic imaging apparatus and control method thereof | |
US20230296746A1 (en) | Systems and methods for scanning | |
WO2023177397A1 (en) | Systems and methods for scanning | |
KR101341092B1 (ko) | 가변 인덕터를 이용한 초음파 진단 방법 및 장치 | |
KR102324390B1 (ko) | 휴대용 하이브리드 초음파 진단장치 | |
AU2022407186A1 (en) | Technologies for ultrasound asynchronous resonance imaging (ari) for needle tip localization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141119 |
|
RJ01 | Rejection of invention patent application after publication |