CN101563626A - 用于最佳散斑跟踪的双路径处理 - Google Patents
用于最佳散斑跟踪的双路径处理 Download PDFInfo
- Publication number
- CN101563626A CN101563626A CNA2007800408907A CN200780040890A CN101563626A CN 101563626 A CN101563626 A CN 101563626A CN A2007800408907 A CNA2007800408907 A CN A2007800408907A CN 200780040890 A CN200780040890 A CN 200780040890A CN 101563626 A CN101563626 A CN 101563626A
- Authority
- CN
- China
- Prior art keywords
- speckle
- line data
- scan
- adopt
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8995—Combining images from different aspect angles, e.g. spatial compounding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明总体上涉及结合了用于散斑跟踪的增强和减弱技术的改进的系统和方法,其用于随着时间的推移获得诸如组织的目标的运动的一系列图像。所述方法包括以下步骤:将声波发射到人体内并输出这些声波的回声;接收所述回声并对其进行射束成形处理以生成扫描线数据;采用减少散斑的方法对扫描线数据进行处理以显示解剖学信息;采用不减少散斑的方法或过程处理扫描线数据;以及在一个扫描序列中,同时获取两种扫描线数据,即,经处理减少了散斑的数据和经处理未减少散斑的数据。
Description
本发明总体上涉及超声成像,更具体地,涉及采用增强的超声散斑图和减弱的超声散斑图两者的超声成像。
在过去的十年间,包括频率复合(compounding)和空间复合(SonoCT)技术的先进复合技术使得超声图像质量得到了显著提高。这些技术通过减弱超声散斑而起作用,所述超声散斑是与由Raleigh散射回声获得的相长/相消干涉图相关的人为噪声图。
散斑是由与包含在超声射束的分辨单元中的很多小的解剖学目标相关的随机相长和相消干涉引起的。根据定义,这些目标或者Raleigh散射体比询问声波的波长短得多。所传输的声束倾向于是宽带的,其包含这样的概念,即,该射束含有具有各种波长的声波。如为本领域技术人员所已知的,不同的波长具有不同的相长和相消干涉图,从而具有不同的散斑图。与棱镜将白光分成其构成波长(颜色)的方式非常类似,正交带通滤波器将返回的声音回声分离成两个组,一组具有较短的波长,另一组具有较长的波长。因此,这两个组将具有不同的干涉图,并从而具有不同的散斑图。
近来,出现了这样的需求,即,(在1D、2D和3D空间内)跟踪血液和组织的速度及其位移两者。由于通过超声成像获得的散斑图倾向于对组织和组织位移进行短距离跟踪,因而可以通过对在空间上获得的散斑图和在时间上获得的类似的散斑图进行交叉相关来计算对组织速度和位移的准确测量。在业内将这些技术称为2D散斑跟踪和3D散斑跟踪。尽管通过减弱返回回声的超声散斑而获得最佳的黑白(BW)图像质量,但是在增强超声散斑时获得最佳的散斑跟踪(位移和速度)。
本发明总体上涉及结合了用于散斑跟踪的增强和减弱技术的改进的系统和方法,其用于随着时间的推移获得诸如组织的目标的运动的一系列图像。
所述方法包括以下步骤:将声波发射到人体内并输出这些声波的回声;接收所述回声并对其进行射束成形处理以生成扫描线数据;采用减少散斑的方法对扫描线数据进行处理以显示解剖学信息;采用不减少散斑的方法或过程处理扫描线数据;以及在一个扫描序列中,同时获取两种扫描线数据,即经处理减少了散斑的数据和经处理未减少散斑的数据。
因此,根据下述说明和权利要求,本发明的上述和其他目的、方面、特征和优点将变得更加明显。
在附图中,类似的附图标记总体上指所有不同的图示中的相同的部分。同样,附图不必按比例绘制,而是附图的重点总体上在于对本发明的原理进行说明。
图1是由低频正交带通滤波器生成的患者(组织)的一部分的示例性散斑图像;
图2是患者(组织)的一部分的示例性合成模型或真实图像;
图3是由低频正交带通滤波器生成的患者(组织)的一部分的另一示例性散斑图像;
图4示出了在后来的时间处的图3的组织;
图5是患者(组织)的一部分的另一示例性合成模型或真实图像;
图6是在后来的时间处图5所示的组织的真实图像;
图7是被配置为减少散斑图的现有技术超声成像系统的示例性示意图;
图8是被配置为用于最佳散斑跟踪的现有技术超声成像系统的示例性示意图;
图9是根据本发明的一个实施例的被配置为用于最佳散斑跟踪的超声成像系统的示例性示意图;
图10是根据本发明的另一实施例的被配置为用于最佳散斑跟踪的超声成像系统的示例性示意图。
本发明涉及采用增强的超声散斑图和减弱的超声散斑图两者的超声成像。复制对应于单条扫描线的返回数字化回声,并将其发送至两个独立的处理路径。针对黑白(BW)图像质量对一条路径进行优化,即,减少散斑。针对散斑跟踪优化另一条路径,即,增强散斑。
在采用超声波照射目标(例如,人体组织)时,所述目标能够与超声信号发生相长或相消干涉。使目标组织的图像看起来呈粒状,或者看起来具有纹理。将这种粒状外观称为散斑。散斑不影响图像中的底层数据。散斑只是数据中随组织运动而变化的任意隆起(bump)或噪声。因而,对散斑进行跟踪,即,随着时间的推移捕获散斑数据能够随着时间的推移对组织移动和/或位移或者血流进行跟踪。
例如,可以如下所述地采用散斑跟踪经过一个心动周期的心脏搏动或运动。在血液正常地流经组织时,组织是软的。在血液没有正常地流经组织时,组织将变硬。心脏是海绵样的,具有可收缩特性。随着心脏搏动,其收缩并还原。但是,死亡的或者受损的组织将不再收缩或运动。因此,随着时间的推移由心脏的超声成像跟踪散斑图能够实现对心脏搏动或运动或其缺失的跟踪。
图1示出了由低频正交带通滤波器生成的患者的一部分,即,组织的散斑图像。图2示出了患者的一部分,即,组织的合成模型或真实图像。在图2中,已经去除了所有的人为散斑(artificial speckle),从而更加容易地看到位于左侧的点目标、位于左上方的小的黑色脉管以及背景灰度级的细微变化(右下)。从而,可以将图2看作是2D和/或解剖学显示中的最佳的。但是,如果组织运动,并且想要探测该位置上的这一位移相对变化,那么图2就没用了,这是由于其缺少任何显著的“纹理”(尤其在右下方)。因而,采用图2进行运动检测将是非常困难的。
在图3中,识别由中间的灰框示出的组织的任意区域以用于跟踪。所述灰框示出了被称为感兴趣区域(ROI)的区域。图4示出了相同的组织,连同所述ROI,但是其位于后来的时间处。如所述ROI所示的,组织已经从图3所示的其初始位置发生了运动。但是,更为重要地,在图3和图4当中,也被称为散斑或颗粒的纹理是相同的。正是这一纹理允许各种“散斑跟踪”方法确定任何给定的组织运动了多远。
图5示出了与图3中相同的组织,并且位于相同的时间处。但是,在这种情况下,已经消除了所有的散斑。并且,已经识别了用于跟踪的具体的组织区域(ROI,灰框)。通过采用与图5中所采用的相同的散斑减少技术,图6示出了位于稍后的时间处的相同的组织。然而,在图6中,由于已经去除了所有的散斑,因而对于任意散斑跟踪方法不能确定图5的ROI中的预期组织运动了多远。从而,消除所有的散斑妨碍了对组织运动的跟踪。
参照图7,其示出了被配置为减少散斑图的现有技术超声成像系统100的示意图。成像系统100包括超声换能器(XD)105、扫描仪110、第一正交带通滤波器(QBP1)115、第二正交带通滤波器(QBP2)120、LogDetect125、LogDetect 130、平均器件135、多速低通滤波器(LPF)140、SonoCT145和显示器150。在优选实施例中,预期扫描仪110已经对返回回声进行数字化,使得采用数字硬件或者采用作为CPU的部分的软件来处理接下来的处理步骤。所述平均器件135可以简单到将LogDetect 125和130的两个输出相加并将所得结果除2。
在操作中,超声换能器(XD)105是将电信号转换成声波并执行反向转换的超声压电换能器。XD 105对对象(患者)进行扫描,生成超声波并将所述超声波输出至扫描仪110,所述扫描仪110是用于引导超声射束的方向并使其聚焦的相位到波(phase to wave)射束形成器。将扫描仪110的输出输入至QBP1 115和QBP2120。QBP1 115和QBP2120是每者包括希尔伯特变换器的带通滤波器(1-3MHz)。QBP1 115的中心为2MHz,QBP2 120的中心为3MHz。QBP1 115和QBP2 120中的每者均输出被称为IQ信号的复解析信号,其具有实同相信号(I)和复正交信号(Q)。通过对平方和取平方根,能够计算出所述回声的包络为:LogDetect 125和LogDetect 130分别接收来自QBP1 115和QBP2 120的复信号,并检测所接收到的复信号的包络,之后对所探测到的结果取对数。注意,对所检测到的来自不同频率的带通滤波器的信号进行组合的方法被称为“频率复合”,这是超声领域的一种成熟完善的技术。
平均器件135接收来自LogDetect 125和LogDetect 130的对数化包络。来自LogDetect 125和LogDetect 130的对数化包络源自于两个不同频率(例如,分别为2MHz和3MHz)。散斑根据频率发生变化,同时底层信号则保持相同。在将所述对数化包络一起取平均时,散斑得到平均。之后,将平均信号输入至多速低通滤波器140并输出至SonoCT 145。由于散斑比底层平均信号变化得更快,因而对这一数据进行低通滤波将进一步减少散斑变化。所述多速低通滤波器140还减少了高空间频率信息,由此实现了对所述信号的抽取(decimate)。其将每扫描线的采样数量从几千个减少到了仅几百个。具有更少的采样减小了下游处理操作的计算负担。
SonoCt 145是一种复合成像设备,其从不同的视角获得图像,之后将其组合成单幅图像。散斑图随着视角的不同而发生变化。将平均器件135的输出输入到多速低通滤波器140和SonoCt 145内的目的在于从超声图像中进一步去除散斑。之后,将SonoCt 145的输出输入到诸如监视器的显示器150中。
参照图8,其示出了被配置为用于最佳散斑跟踪的现有技术超声成像系统200的示意图。成像系统200包括超声换能器(XD)105、扫描仪110、正交带通滤波器(QBP1)115、LogDetect 125、多速低通滤波器202、散斑跟踪器205和显示器150。
在操作中,XD 105生成超声波,并将其输出至扫描仪110。将扫描仪110的输出输入至QBP1 115。QBP1 115输出如上所述的IQ信号。LogDetect125接收来自QBP1 115的复信号,并检测所接收到的复信号的包络。之后,将所述包络输入至多速低通滤波器202,并输出至散斑跟踪器205。与图7中采用的用于最佳散斑减少的多速低通滤波器140不同,这一多速低通滤波器202提供了更少的平滑,并潜在地提供了更少的抽取。为了实现最佳散斑跟踪,希望增强散斑,因而用于屏蔽散斑的现有技术现在变得不利了。散斑跟踪器205是一种交叉相关设备,其跟踪处于不同时间点上的散斑,即,随着目标(例如,组织)的运动而记录图像数据,以获得散斑的变化。通过使处于不同的时间点上的散斑交叉相关,散斑跟踪器能够计算组织位移、组织运动和组织收缩。之后,将散斑跟踪器205的输出输入至显示器150。
存在各种已公开的“散斑跟踪”方法(例如,Chen等人的美国专利No.5876342)。典型的用于对行动和速度二者进行跟踪的方法是“归一化交叉相关”法。其被定义如下:
其中:
NCC是归一化交叉相关函数
dx、dy是用于确定散斑运动了多远的搜索空间
x∈ROI对感兴趣区域(ROI)中的x&y求和
u1是位于时间1处的图像
u2是位于时间2处的图像
该方程的应用如下:
1、首先,识别第一图像中的所选择的用于进行跟踪的感兴趣区域(ROI)。注意,可以选择多个ROI,并且可以选择每一像素(或者3D体积中的每一体素)进行跟踪。其定义了ROI以及x和y在第一图像u1中的范围。
2、接下来,改变dx和dy,从而对稍后的时间处观察到的图像u2中的相同大小的ROI进行位移。
3、对于每一dx和dy,估计归一化交叉相关(NCC)函数。
4、重复步骤2和3,直到观察到NCC的峰最大值为止。NCC值1.0表示最大相关。这一峰值处的dx和dy的值表示所述ROI中的期望组织移动了多远。
如对本领域的技术人员显而易见地,在源ROI(u1中)或者发生位移的ROI(u2中)内缺乏任何纹理或者散斑变化将使得NCC搜索算法失败。因为对于dx和dy的所有位移值而言都观察到了1.0的相关值,因而不能识别出峰值。
本发明提供了改进的系统和方法,其用于将从图像增强超声信号路径获得的数据和从散斑增强超声信号路径获得的数据相结合,从而随着时间的推移获得组织运动的一系列图像。
参照图9,其示出了被配置为用于最佳散斑跟踪的超声成像系统300的优选实施例的示意图。所述成像系统300包括超声换能器(XD)105、扫描仪/射束形成器110、第一正交带通滤波器(QBP1)115、第二正交带通滤波器(QBP2)120、LogDetect 125、LogDetect 130、平均器件135、第一多速低通滤波器305、第二多速低通滤波器310、散斑跟踪器205、SonoCT145和显示器150。
在操作中,扫描仪110向超声换能器XD 105发送电信号,超声换能器XD 105将这一电信号转换成声波。将这些声波传播到身体内,并使其在各解剖学结构上发生反射。通过同一超声换能器XD 105将返回的声波回声转换回电信号,之后,将其发送回扫描仪110。之后,扫描仪110对这些信号进行处理,以隔离来自特定扫描方向和深度的回声,由此确定处于这些位置的解剖学结构。
将扫描仪110的输出输入至QBP1 115和QBP2 120。在一个实施例中,QBP1 115的中心为2MHz,QBP2 120的中心为3MHz。QBP1 115和QBP2120每者输出IQ信号,所述IQ信号是从中去除了信号噪声的复信号。LogDetect 125和LogDetect 130分别接收来自QBP1 115和QBP2 120的复信号,并检测所接收到的复信号的包络。平均器件135经由信号路径320接收来自LogDetect 125的信号包络,并接收来自LogDetect 130的信号包络,并由所述图像求得平均噪声(散斑),如上文所述。
之后,将所述平均信号输入到多速低通滤波器310内。将多速低通滤波器310的输出输入至SonoCT 145,SonoCT 145获得来自不同视角的图像,之后将其组合成单幅图像。之后,将SonoCt 145的输出输入至显示器150。
还将来自LogDetect 125的信号包络经由信号路径315输入至多速低通滤波器305。将多速低通滤波器305的输出输入至散斑跟踪器205,所述散斑跟踪器跟踪处于不同时间点上的散斑。如上所述地,通过使处于不同的时间点上的散斑交叉相关,散斑跟踪器能够计算组织位移、组织运动和组织收缩。之后,将散斑跟踪器205的输出输入至显示器150。
显示器150同时获得来自散斑跟踪器205的散斑数据和来自SonoCT145的图像数据。可以使这一散斑数据或“功能信息”作为曲线图或辅助图像与解剖学图像数据并排显示。在优选实施例中,例如,可以采用不同于解剖学图像的颜色使这一功能信息覆盖或叠加在解剖学图像数据之上。在超声领域,通常将这样的图像称为“参数图像”。
因而,可以将散斑数据叠加到图像数据上,以创建参数图像,其允许对所成像的组织的运动的观察。可以基于散斑数据的值用各种颜色显示所述散斑数据。例如,在一个实施例中,将指示运动的组织的“不断变化”的散斑数据显示为绿色,将指示非运动的组织的“非变化”的散斑数据显示为灰色。有利地,在将“着色”的散斑数据叠加到同时获得的图像数据上时,能够观察到运动的组织和非运动的组织。除了组织运动之外,还可以采用所获得的散斑数据和图像数据观察血流。随着血液流动,组织随时间扩张收缩,从而引起不断变化的散斑数据。如果不存在血流,那么所获得的散斑数据将不会发生变化。
散斑跟踪器205的直接输出提供了针对所询问的解剖学结构的运动和位移信息。可以采用这一信息确定很多功能属性。在一个示例中,可以关于时间对位移场(field)求微分,以确定不同结构的速度。在另一示例中,可以采用位移的空间差异计算局部应变。可以采用这样的应变量度将心肌的健康、收缩的部分和缺血、坏死、非收缩的部分区分开。在另一示例中,可以采用运动场进行定时分析,以确定心脏的不同部分何时收缩。在正常的健康心脏中,左心室的所有部分都倾向于同时收缩。但是,在非同步收缩的患有疾病的心脏中,心肌的不同部分在不同的时间收缩,从而导致有效泵送(pumping)减少。
可以采用专用硬件或者在计算机中运行的软件计算所有的上述导出量度。而且,有可能实时(当正在获取声波时)或者非实时(获取后)地导出这样的量度。
可以将本发明的上述系统和方法用于检测乳房组织中的肿瘤。只有当肿瘤被较不致密的组织包绕时,诸如,在四十到五十岁的妇女身上时,诸如乳房X线照相术的现有方法才有效。而本发明则不管包绕组织密度如何都能有效地检测到肿瘤,即,不存在血流或组织移动的区域,因而能够检测到二十到四十岁的妇女体内的肿瘤。
本发明的上述方法还可以起到发现心脏的梗死(infracted)区域的作用。这样的区域已经受损并且血流减少,因此运动减弱,可以对其运行进行跟踪和观察。
此外,与包括电离辐射或者放射性染料的引入的当前诊断方法相比,本发明更加快速、安全,并且介入性更低。
图9所示的实施例的一个关键限制在于,散斑减少图像质量路径和最佳散斑跟踪路径两者共享QBP-滤波器-LogDetect处理组(bank)之一(例如,QBP滤波器115和LogDetect 125)。然而这一共享由于只需要两个QBP-滤波器-LogDetect处理体,因而能够得到更低成本的实现,但是其潜在地损害了散斑减少图像质量路径和最佳散斑跟踪路径两者的性能。例如,可能希望将所述路径之一配置成用于基频操作(QBP滤波器具有接近传输频率的中心频率),而将另一路径配置为用于组织谐波成像(QBP滤波器具有二倍于传输频率的中心频率)。
参照图10,其示出了在解决图9的性能限制的备选实施例中被配置成用于最佳散斑跟踪的超声成像系统400的示意图。成像系统400包括超声换能器(XD)105、扫描仪/射束形成器110、第一正交带通滤波器(QBP1)115、第二正交带通滤波器(QBP2)120、第三正交带通滤波器(QBP3)405、LogDetect 125、LogDetect 130、LogicDetect 410、平均器件135、第一多速低通滤波器305、第二低速低通滤波器310、散斑跟踪器205、SonoCT 145和显示器150。
在操作中,XD 105将超声波转换成电信号,并将其输出至扫描仪110。将扫描仪110的输出输入至QBP1 115、QBP2 120和QBP3 405。在一个实施例中,QBP1 115的中心为2MHz,QBP2 120的中心为3MHz。其可能涉及这样一种情形,即,传输频率以2.5MHz为中心,而QBP1 115和QBP2 120试图在接近所述传输频率的基频上执行频率复合。这些频率可以是针对最佳图像质量和最佳散斑减少而选定的。在这一相同情形下,可以推断应当采用组织谐波成像执行最佳散斑跟踪(参照美国专利No.5879303)。在这种情况下,使QBP3 405的中心位于5MHz处是适当的,该频率是所传输的声波的频率的二倍。QBP1 115、QBP2 120和QBP3 405中的每者输出IQ信号,所述IQ信号是从中去除了信号噪声的复信号。LogDetect 125和LogDetect 130分别接收来自QBP1 115和QBP2 120的复信号,并检测所接收到的复信号的包络。LogDetect 410接收来自QBP3405的复信号,并检测所接收到的复信号的包络。
平均器件135接收来自LogDetect 125和LogDetect 130的信号包络,并对来自图像的散斑求平均。之后,将所述平均信号输入至多速低通滤波器310。将多速低通滤波器310的输出输入至SonoCT 145。之后,将SonoCt145的输出输入至诸如监视器的显示器150。
同时,将来自LogDetect 410的信号包络输入至多速低通滤波器305。将多速低通滤波器305的输出输入至散斑跟踪器205。之后,将散斑跟踪器205的输出输入至显示器150。
所有的实施例和方框图描述了处理相同扫描线的不同方法,从而针对最佳图像质量和散斑减少对一条路径进行优化,并针对散斑跟踪和散斑增强对第二条路径进行优化。将扫描线定义为询问人体内的特定视线的单个声射束,其具有轴向深度尺寸(例如,以mm为单位)。根据如何对这种扫描线进行排序,能够获得不同的成像模式和显示。在一个实施例中,所述扫描线可以询问同一条视线(其被称为M模式)。在第二实施例中,所述扫描线可以贯穿身体内的层析成像片层排序,其被称为2D或B模式操作。在又一实施例中,所述扫描线可以根据方位(横向)和高度尺寸两者发生变化,从而进行体积扫描(其被称为3D或4D成像)。
同样,对本领域技术人员显而易见的,本发明适于任何类型的超声换能器,包括但不限于单元件机械换能器、相控阵列、线性装置(linears)、弯曲线性(curved-linear)阵列(CLA)、2D矩阵阵列和相控阵列摇动器(wobbler)。
在本发明的另一实施例中,假设使所述并行处理路径时分复用,从而在逐条线的基础上改变单个处理路径,使得在一接收扫描事件的过程中,针对散斑跟踪对路径进行优化,而在针对另一接收扫描事件的过程中,所述线可以是相同的视线,可以将该路径优化为获得具有减弱的散斑的最佳图像质量。
在本发明的另一实施例中,用于散斑跟踪的处理包括将RF滤波器用于带通滤波器,而不是具有LogDetect。本发明的另一实施例包括通过将检测后的低通滤波器限定在低于散斑跟踪路径中所采用截止频率的截止频率上来减少散斑。
在不脱离本发明的精神和范围的情况下,本领域技术人员可以想到文中描述的内容的各种变化、修改和其他实现。相应地,本发明不仅受前面的示意性说明的限定。
Claims (22)
1、一种超声系统,包括:
用于将声波发射到人体内并输出所述声波的回声的器件;
用于接收所述回声并对其进行射束成形处理以生成至少一个扫描线数据的器件;
用于处理所述扫描线数据之一以显示解剖学信息的第一器件,所述用于处理的第一器件包括减少散斑;
用于处理所述扫描线数据之一的第二器件,所述用于处理的第二器件不包括减少散斑;以及
用于在一个扫描序列中同时获取采用所述第一器件处理的所述扫描线数据之一和采用所述第二器件处理的所述扫描线数据之一的器件。
2、根据权利要求1所述的系统,其中,所述第一处理器件和所述第二处理器件处理相同的扫描线数据。
3、根据权利要求1所述的系统,其中,所述第一处理器件和所述第二处理器件处理来自一个扫描序列的不同扫描线数据。
4、根据权利要求1所述的系统,其中,所述第一处理器件包括RF带通滤波器。
5、根据权利要求1所述的系统,其中,所述第二处理器件包括RF带通滤波器。
6、根据权利要求1所述的系统,其中,所述第一处理器件包括用于检测所述回声的包络的器件。
7、根据权利要求6所述的系统,其中,所述第一处理器件包括用于对所检测到的包络求对数的器件。
8、根据权利要求1所述的系统,其中,采用频率复合实现减少散斑。
9、根据权利要求8所述的系统,还包括两个或更多滤波器组,每一滤波器组包括探测器和具有针对每一滤波器组的唯一响应的带通滤波器。
10、根据权利要求1所述的系统,其中,通过空间复合实现减少散斑。
11、根据权利要求1所述的系统,其中,所述扫描序列询问单条视线、平面和体积的至少其中之一。
12、根据权利要求1所述的系统,其中,所述的用于发射声波的器件选自由相控阵列、线性装置、弯曲线性阵列、机械摇动器和3D摇动器构成的组。
13、根据权利要求1所述的系统,其中,所述第一处理器件和所述第二处理器件包括RF带通滤波器。
14、根据权利要求1所述的系统,其中,可以采用专用硬件或者采用在CPU内运行的软件实现所述第一处理器件。
15、根据权利要求1所述的系统,其中,可以采用专用硬件或者采用在CPU内运行的软件实现所述第二处理器件。
16、根据权利要求1所述的系统,其中,可以采用参数成像显示技术将采用所述第二处理器件处理的所述扫描线数据之一叠加到所述解剖学信息上。
17、根据权利要求1所述的系统,其中,使采用所述第二处理器件处理的所述扫描线数据之一与从先前扫描序列获取的数据交叉相关。
18、根据权利要求17所述的系统,其中,可以采用所述经过交叉相关的数据求得应变、应变率、弹性图像、壁增厚和收缩定时。
19、根据权利要求1所述的系统,其中,通过将检测后低通滤波器限定在低于散斑跟踪路径中采用的截止频率的截止频率上实现减少散斑。
20、根据权利要求1所述的系统,其中,可以重复一个扫描序列,以确定组织随着时间的推移而发生的空间位移。
21、一种用于执行散斑跟踪的方法,所述方法包括如下步骤:
将声波发射到人体内并输出所述声波的回声;
接收所述回声并对其进行射束成形处理以生成至少一个扫描线数据;
对所述扫描线数据之一进行处理以显示解剖学信息,所述处理包括减少散斑;
对所述扫描线数据之一进行额外处理,所述额外处理不包括减少散斑;以及
在一个扫描序列中同时获取经过处理以显示解剖学信息的所述扫描线数据之一和采用额外处理处理过的所述扫描线数据之一。
22、一种具有用于在计算机上运行以执行散斑跟踪的计算机可读程序代码的计算机可读介质,所述计算机可读程序代码包括:
将声波发射到人体内并输出所述声波的回声;
接收所述回声并对其进行射束成形处理以生成至少一个扫描线数据;
对所述扫描线数据之一进行处理以显示解剖学信息,所述处理包括减少散斑;
对所述扫描线数据之一进行额外处理,所述额外处理不包括减少散斑;以及
在一个扫描序列中同时获取经过处理以显示解剖学信息的所述扫描线数据之一和采用额外处理处理过的所述扫描线数据之一。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US86425906P | 2006-11-03 | 2006-11-03 | |
US60/864,259 | 2006-11-03 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101563626A true CN101563626A (zh) | 2009-10-21 |
Family
ID=39301157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007800408907A Pending CN101563626A (zh) | 2006-11-03 | 2007-11-02 | 用于最佳散斑跟踪的双路径处理 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20100004540A1 (zh) |
EP (1) | EP2082261A2 (zh) |
JP (1) | JP5627890B2 (zh) |
KR (1) | KR101407425B1 (zh) |
CN (1) | CN101563626A (zh) |
WO (1) | WO2008053457A2 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103156647A (zh) * | 2011-12-08 | 2013-06-19 | 株式会社东芝 | 超声波诊断装置、图像处理装置及图像处理方法 |
CN103747742A (zh) * | 2011-04-14 | 2014-04-23 | 明尼苏达大学评议会 | 使用超声成像的脉管表征 |
CN104125801A (zh) * | 2011-10-03 | 2014-10-29 | 海浪科技有限公司 | 利用双带脉冲复合体的非线性成像 |
CN107303186A (zh) * | 2016-04-19 | 2017-10-31 | 美国西门子医疗解决公司 | 弹性成像中的频率复合 |
CN109803588A (zh) * | 2016-09-02 | 2019-05-24 | 挪威科技大学 | 流体路径的增强分辨率超声成像 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010068450A1 (en) * | 2008-11-25 | 2010-06-17 | Mayo Foundation For Medical Education And Research | System and method for analyzing carpal tunnel using ultrasound imaging |
KR101120820B1 (ko) | 2009-11-19 | 2012-03-22 | 삼성메디슨 주식회사 | 초음파 공간 합성 영상을 제공하는 초음파 시스템 및 방법 |
EP2385391A3 (en) | 2010-05-04 | 2012-08-01 | Sony Corporation | Active imaging device and method for speckle noise reduction |
US10231712B2 (en) | 2010-06-09 | 2019-03-19 | Regents Of The University Of Minnesota | Dual mode ultrasound transducer (DMUT) system and method for controlling delivery of ultrasound therapy |
JP2012176232A (ja) * | 2011-02-04 | 2012-09-13 | Toshiba Corp | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
TWI446897B (zh) * | 2011-08-19 | 2014-08-01 | Ind Tech Res Inst | 超音波影像對齊裝置及其方法 |
TWI487402B (zh) * | 2012-08-10 | 2015-06-01 | Mstar Semiconductor Inc | 可用於一無線通訊系統的搜尋方法 |
WO2014109392A1 (ja) * | 2013-01-11 | 2014-07-17 | 日立アロカメディカル株式会社 | 超音波撮像装置 |
WO2014136016A1 (en) | 2013-03-05 | 2014-09-12 | Koninklijke Philips N.V. | Consistent sequential ultrasound acquisitions for intra-cranial monitoring |
CN109044407A (zh) | 2013-07-23 | 2018-12-21 | 明尼苏达大学评议会 | 使用多频率波形的超声图像形成和/或重建 |
WO2016120763A1 (en) * | 2015-01-29 | 2016-08-04 | Koninklijke Philips N.V. | Evaluation of cardiac infarction by real time ultrasonic strain imaging |
KR102387708B1 (ko) | 2015-01-30 | 2022-04-19 | 삼성메디슨 주식회사 | 향상된 hprf 도플러 영상을 위한 가이드를 제공하는 방법 및 초음파 시스템 |
CN104586433B (zh) * | 2015-02-02 | 2016-08-24 | 声泰特(成都)科技有限公司 | 基于变频的基波/谐波融合与空间复合相结合的成像方法 |
WO2018051265A1 (en) | 2016-09-15 | 2018-03-22 | Koninklijke Philips N.V. | Ultrasonic elastographic pre-load measurement and display |
US11458337B2 (en) | 2017-11-28 | 2022-10-04 | Regents Of The University Of Minnesota | Adaptive refocusing of ultrasound transducer arrays using image data |
US11596812B2 (en) | 2018-04-06 | 2023-03-07 | Regents Of The University Of Minnesota | Wearable transcranial dual-mode ultrasound transducers for neuromodulation |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1715360A2 (en) * | 2005-03-31 | 2006-10-25 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus and ultrasound image processing program |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2005208A (en) * | 1934-08-20 | 1935-06-18 | Schuchat Jonas | Wrapping |
EP0321717B1 (en) * | 1987-11-30 | 1995-07-05 | Aloka Co. Ltd. | Ultrasonic speckle velocity measurement method and apparatus |
US5474070A (en) * | 1989-11-17 | 1995-12-12 | The Board Of Regents Of The University Of Texas System | Method and apparatus for elastographic measurement and imaging |
US5879303A (en) * | 1996-09-27 | 1999-03-09 | Atl Ultrasound | Ultrasonic diagnostic imaging of response frequency differing from transmit frequency |
US5735797A (en) * | 1996-12-30 | 1998-04-07 | General Electric Company | Method and apparatus for combining topographic flow power imagery with a B-mode anatomical imagery |
US5876342A (en) * | 1997-06-30 | 1999-03-02 | Siemens Medical Systems, Inc. | System and method for 3-D ultrasound imaging and motion estimation |
US6213946B1 (en) * | 1998-12-24 | 2001-04-10 | Agilent Technologies, Inc. | Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging |
US6142942A (en) * | 1999-03-22 | 2000-11-07 | Agilent Technologies, Inc. | Ultrasound imaging system and method employing an adaptive filter |
US6464637B1 (en) * | 2000-06-23 | 2002-10-15 | Koninklijke Philips Electronics N.V. | Automatic flow angle correction by ultrasonic vector |
JP2004073620A (ja) * | 2002-08-21 | 2004-03-11 | Toshiba Medical System Co Ltd | 超音波診断装置 |
CA2445137A1 (en) * | 2002-10-15 | 2004-04-15 | Matsushita Electric Industrial Co., Ltd. | Image processing apparatus, method and program |
US20050053305A1 (en) * | 2003-09-10 | 2005-03-10 | Yadong Li | Systems and methods for implementing a speckle reduction filter |
JP4590256B2 (ja) * | 2004-05-20 | 2010-12-01 | 富士フイルム株式会社 | 超音波撮像装置、超音波画像処理方法、及び、超音波画像処理プログラム |
US20050288589A1 (en) * | 2004-06-25 | 2005-12-29 | Siemens Medical Solutions Usa, Inc. | Surface model parametric ultrasound imaging |
US7678050B2 (en) * | 2004-08-24 | 2010-03-16 | General Electric Company | Method and apparatus for detecting cardiac events |
-
2007
- 2007-11-02 KR KR1020097008907A patent/KR101407425B1/ko active IP Right Grant
- 2007-11-02 US US12/447,969 patent/US20100004540A1/en not_active Abandoned
- 2007-11-02 WO PCT/IB2007/054466 patent/WO2008053457A2/en active Application Filing
- 2007-11-02 EP EP07826970A patent/EP2082261A2/en not_active Withdrawn
- 2007-11-02 JP JP2009535182A patent/JP5627890B2/ja not_active Expired - Fee Related
- 2007-11-02 CN CNA2007800408907A patent/CN101563626A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1715360A2 (en) * | 2005-03-31 | 2006-10-25 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus and ultrasound image processing program |
Non-Patent Citations (1)
Title |
---|
ABD EL-MONEM EL-SHARKAWY等: "improved ultrasound speckle motion tracking using nonlinear diffusion filtering", 《MEDICAL IMAGING 2001:ULTRASONIC IMAGING AND SIGNAL PROCESSING(PROCEEDING VOLUME)》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103747742A (zh) * | 2011-04-14 | 2014-04-23 | 明尼苏达大学评议会 | 使用超声成像的脉管表征 |
CN103747742B (zh) * | 2011-04-14 | 2016-04-06 | 明尼苏达大学评议会 | 使用超声成像的脉管表征 |
CN104125801A (zh) * | 2011-10-03 | 2014-10-29 | 海浪科技有限公司 | 利用双带脉冲复合体的非线性成像 |
CN103156647A (zh) * | 2011-12-08 | 2013-06-19 | 株式会社东芝 | 超声波诊断装置、图像处理装置及图像处理方法 |
CN103156647B (zh) * | 2011-12-08 | 2016-03-16 | 株式会社东芝 | 超声波诊断装置、图像处理装置及图像处理方法 |
US9345446B2 (en) | 2011-12-08 | 2016-05-24 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method |
CN107303186A (zh) * | 2016-04-19 | 2017-10-31 | 美国西门子医疗解决公司 | 弹性成像中的频率复合 |
US10675007B2 (en) | 2016-04-19 | 2020-06-09 | Siemens Medical Solutions Usa, Inc. | Frequency compounding in elasticity imaging |
CN109803588A (zh) * | 2016-09-02 | 2019-05-24 | 挪威科技大学 | 流体路径的增强分辨率超声成像 |
CN109803588B (zh) * | 2016-09-02 | 2021-09-28 | 挪威科技大学 | 流体路径的增强分辨率超声成像 |
Also Published As
Publication number | Publication date |
---|---|
WO2008053457A3 (en) | 2008-07-03 |
JP5627890B2 (ja) | 2014-11-19 |
US20100004540A1 (en) | 2010-01-07 |
JP2010508881A (ja) | 2010-03-25 |
EP2082261A2 (en) | 2009-07-29 |
WO2008053457A2 (en) | 2008-05-08 |
KR20090084840A (ko) | 2009-08-05 |
KR101407425B1 (ko) | 2014-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101563626A (zh) | 用于最佳散斑跟踪的双路径处理 | |
JP6934497B2 (ja) | 多開口超音波撮像を用いた点音源送信及び音速補正 | |
US10338203B2 (en) | Classification preprocessing in medical ultrasound shear wave imaging | |
JP6063553B2 (ja) | 超音波イメージング方法及び超音波イメージング装置 | |
US10959704B2 (en) | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method | |
EP1579244B1 (en) | Segmentation tool for identifying flow regions in an imaging system | |
CN101951839B (zh) | 利用超声造影剂的呼吸门控治疗评估 | |
US11471130B2 (en) | Method and ultrasound system for shear wave elasticity imaging | |
US20100286516A1 (en) | High pulse repetition frequency for detection of tissue mechanical property with ultrasound | |
EP2022404A1 (en) | Ultrasonograph, medical image processing device, and medical image processing program | |
JP6063552B2 (ja) | 超音波イメージング方法及び超音波イメージング装置 | |
US9173629B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus | |
CN102458260A (zh) | 测量心脏收缩功能的方法和装置 | |
EP3769690B1 (en) | Method for tissue characterization by ultrasound wave attenuation measurements and ultrasound system for tissue characterization | |
US9161737B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus | |
CN107303186B (zh) | 弹性成像中的频率复合 | |
US11006928B2 (en) | Sound speed imaging using shear waves | |
US9448100B2 (en) | Signal processing apparatus | |
CN105025807B (zh) | 用来对静止流体池进行识别和定界的扫描 | |
EP2138102B1 (en) | Formation of an enhanced elastic image in an ultrasound system | |
US20190239861A1 (en) | Ultrasonic diagnostic apparatus | |
CN108135573A (zh) | 超声波诊断系统及超声波诊断方法 | |
CN114246610A (zh) | 用于生成增强图像的超声成像系统和方法 | |
Onyia | Displacement Data Processing for ARFI Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20091021 |