CN101595396A - 高分辨率超声位移测量的装置与方法 - Google Patents

高分辨率超声位移测量的装置与方法 Download PDF

Info

Publication number
CN101595396A
CN101595396A CNA2007800308186A CN200780030818A CN101595396A CN 101595396 A CN101595396 A CN 101595396A CN A2007800308186 A CNA2007800308186 A CN A2007800308186A CN 200780030818 A CN200780030818 A CN 200780030818A CN 101595396 A CN101595396 A CN 101595396A
Authority
CN
China
Prior art keywords
displacement
echo
signal
gcs
clutter
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
Application number
CNA2007800308186A
Other languages
English (en)
Inventor
彭旗宇
张立群
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Rehab Tek
Original Assignee
Rehab Tek
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Rehab Tek filed Critical Rehab Tek
Publication of CN101595396A publication Critical patent/CN101595396A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/06Systems determining the position data of a target
    • G01S15/08Systems for measuring distance only
    • G01S15/10Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
    • G01S15/102Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics
    • G01S15/104Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications

Abstract

高分辨率位移测量的超声波方法包括发射一对Golay互补序列或其他的互补序列,接收物体(或物体内某一区域)的回波,再对这对回波序列进行脉冲压缩。压缩脉冲输出的主瓣信号周围的剩余杂波信号被用来确定物体(或物体内某一区域)在这两个序列发射之间的位移。此外,物体(或物体内某一区域)的运动速度,物体厚度,物体应变,物体的弹性刚度,及粘性阻尼也可随后确定。

Description

高分辨率超声位移测量的装置与方法
相关专利申请的交叉引用
本申请要求2006年6月20日提交的美国专利申请No.11/471,129的优先权。
技术领域
本发明涉及位移,速度,厚度和应变的非侵入测量,及其在医疗和工业应用。
传统的用于测量位移的脉冲回波(echo)技术中,发射超声波(雷达或声纳),并从目标接收回波。通过比较两个回波之间的时移、频率变化或相位变化来测量位移。例如,传统的测量角膜厚度(约500纳米)超声方法根据来自角膜两侧的回波的两个峰值之间的时间间隔来获得角膜厚度。传统的超声血流速度测量方法将测量时域或频域中的多普勒频率变化。传统的超声波弹性的测量和成像通过时域或频域中的相关性来检测回波的时移,由此获得物体内的应变或位移。
基于传统的脉冲-回波技术的位移测量的分辨率受发射波的中心频率,回波的采样率和测量系统中的电子噪声所限制。对于现有技术而言,有必要克服这些限制,从而提高精确度,分辨率,响应率,灵敏度,抗噪音和速度。
发明内容
在以下所提供的详细阐述将展示本发明的适用的其他领域。以下提供的具体描述和例子,仅以阐述本发明为目的。本发明的实际应用并不仅限于此。
概括的说,本发明提供了一种以高分辨率确定物体或物体中的期望区域的位移的系统,并且所述测量对于噪声具有高鲁棒性。此系统包括发射一对Golay互补序列(GCS),接收来自同一物体或物体的相同区域的两个回波,压缩脉冲和去除主瓣。通过压缩脉冲输出的主瓣周围的剩余杂波信号来得到两次发射之间的物体的位移。
发明的具体实现可以包括一个或几个以下特征。信号发生器电路产生Golay互补序列以控制波发射器装置。波发射器可以是超声波换能器、雷达天线或声纳发射器。来自物体目标区域的回波随后被接收器接收。在某些情况下,波发射器装置和接收装置可以是相同的装置。信号调整电路执行前置放大,回波衰减补偿,回波的带通滤波。模数转换器(ADC)将模拟回波采样为数字回波。数字信号处理装置或计算机执行脉冲压缩,主瓣去除,杂波收集和位移计算(图1)。
在另一个方面,本发明提供了一种确定物体或者物体中的期望区域的速度的系统。该系统包括发射一对GCS序列(A和B)。接收来自相同的物体或物体的相同区域的两个相应的回波。执行脉冲压缩和主瓣去除。然后,根据压缩脉冲输出的主瓣周围的剩余杂波来获得两个发射之间的物体或期望区域的位移。物体或期望区域的速度可从这两个发射之间的位移和时间间隔得到。
在另一个方面,本发明提供了一种确定物体或物体中的期望区域的厚度的系统。该系统包括发射一对GCS序列。然后接收来自同一物体或相同区域两侧的两个相应的回波,并随后进行脉冲压缩和主瓣去除。随后,通过压缩脉冲输出的主瓣周围的剩余杂波信号得到物体或期望区域的厚度。
在另一个方面,本发明提供了一种以高分辨率确定物体或物体中期望区域的应变的系统,并且所述测量具有高鲁棒性。该系统包括发射一对Golay互补序列。然后接收来自同一物体或物体的相同区域的两个相应的回波,并随后进行脉冲压缩和主瓣去除。通过压缩脉冲输出的主瓣周围的剩余杂波信号获得两个Golay发射之间的物体的区域或部分的应变。
在另一个方面中,本发明提供了一种以高分辨率和高鲁棒性的测量来对物体或物体的期望区域进行2维或3维位移映射的系统。该系统包括发射一对Golay互补序列。接收来自相同的物体或物体的相同区域的两个回波。然后,执行脉冲压缩和主瓣去除。然后,根据压缩脉冲输出的主瓣周围的剩余杂波信号获得两发射之间的物体区域的位移。随后,根据在2维或3维空间中不同位置处测量的位移获得位移的2维或3维映射。
在另一个方面中,本发明提供了对物体或物体的期望区域进行2维或3维速度映射的系统。该系统包括发射一对Golay互补序列。接收来自相同的物体或相同区域的相应的两次回波。然后,执行脉冲压缩和主瓣去除。然后,从压缩脉冲输出的主瓣周围的剩余杂波信号获得两发射之间物体或期望区域的位移。物体或期望区域的速度可随后根据这两个发射之间的位移和时间间隔导出。根据在2维或3维空间中的不同位置处测量的速度来获得速度的2维或3维映射。
在另一个方面中,本发明提供了以高分辨率和鲁棒性的测量对物体或物体的期望区域进行2维或3维应变映射的系统。该系统包括发射一对Golay互补序列。来自相同的物体或物体的相同区域的两个相应的回波接收。然后,执行脉冲压缩和主瓣去除。根据压缩脉冲输出的主瓣周围的剩余杂波信号获得两发射之间的散射体位移。物体或期望区域的应变随后由物体或期望区域的长度和位移导出。
一种数字信号处理装置或计算机,实施脉冲压缩,主瓣去除,杂波收集和位移计算。通过2维或3维空间中不同位置处测量的应变得出应变的2维或3维映射。
发明的各种不同方面可能提供一个或更多个以下特点。能够通过物体或物体中区域的位移,速度,厚度,应变导出物体或物体的某一区域的机械属性,如弹性刚度和粘性阻尼。
发明的上述以及其他特点,连同发明本身一起,可以通过下列图示,详细的描述和权利要求得到充分了解。
附图说明
图1是以Golay互补序列的杂波信号为基础的位移测量的一个简单图示。
图2是使用GCS脉冲压缩的示例。图3是物体运动对压缩脉冲的影响:(a)GCS A和GCS B之间物体不发生移动;(b)GCS A和GCS B之间物体移动7.5微米;(c)GCS A和GCS B之间物体移动15微米;(d)GCS A和GCS B之间物体移动30微米。
图4是在没有进行对数压缩的超声图像中组织运动的效应。图像尺寸大约是5毫米×5毫米。(a)GCS A和GCS B之间组织不发生移动;(b)GCS A和GCS B之间组织移动7.5微米;(c)GCS A和GCS B之间组织移动15微米;(d)GCS A和GCS B之间组织移动30微米。
图5是两个GCS激励之间的组织运动和被压缩脉冲的SCR之间的关系。用于模拟的换能器的中心频率f0和-3分贝带宽η分别是10MHz和70%。GCS的长度N和宽度Tc分别是256和100纳秒。数字回波的采样率fs是100MHz。
图6(a)示出了来自GCS A(实线表示)和GCS B(虚线表示)的回波的自相关的部分波形。虚线波形移动了GCS码宽度5/10。图6(b)示出了由不同时移所产生的杂波(GCS码宽度的1/10,2/10,3/10,4/10和5/10)。CGS的长度N是256。X轴是自相关函数的采样点n。
图7示出了不同回波噪声水平下的Sac-位移曲线。仿真参数设定与图3相同。在每条曲线中,当位移d等于±7.5微米,±15微米,…,75微米时的Sac值用星号标记。回波幅度被归一化为1。(a):没有噪声;(b)~(f):随机噪声均匀分布在间隔(-0.2,0.2),(-0.4,0.4),(-0.6,0.6),(-0.8,0.8)和(-1,1)中。
图8示出了计算单个物体位移的算法流程。
图9示出了实际系统的分辨率。(a)实际的位移和基于杂波信号的测量之间的关系。(b)测量误差的分布的直方图。
图10示出了物体中的期望区域中的位移测量算法。
图11示出了预白化滤波器的系数(GCS的长度是256)。
图12示出了仿真中使用的散射体的分布。
图13示出了(a)GCS A的回波。(b)GCS B的回波。(c)被压缩的脉冲。(d)主瓣去除后获得的杂波信号。(e)解卷积的结果。(f)幅度归一和解调的结果。
图14示出了杂波(幅度归一后的结果)的幅度和直线电机控制的实际位移之间的关系。
图15示出了仿真散射体的设定。
图16示出了(a)GCS A的回波。(b)GCS B的回波。(c)被压缩的脉冲。(d)主瓣去除后获得的杂波信号。(e)解卷积的结果。(f)幅度归一的结果。
图17示出了(a)实际轴位移设定的绝对值,图16(f)所示波形的包络.(c)位移的仿真2维分布的图像。
具体实施方式
以下对优选实施例的描述本质上只是示意性的,并不是对本发明及其应用的限制。
本发明的实施例为医疗和工业应用中位移、速度、厚度以及应变提供了非侵入检测的技术。比如,发明中的系统可以在内部外部扰动下在活的有机体内检测内部的位移或应变,所述位移或应变与组织的弹性模量或刚度有关。为了测量位移,用一对互补序列激发超声换能器。在所示的实施例中,所使用的是Golay互补序列(GCS),且受控的时间间隔与待测组织的位移以及超声在组织中的速度相关。来自组织的两个回波随后由相同或另一个超声换能器所接收。在一些情形中,回波需要被预先放大以补偿在组织中的能量损失。然后,模拟的回波由模-数转换器采样并传送到电子设备(微处理器、数字信号处理器或者FPGA)。计算机进行脉冲压缩,主瓣去除,杂波积累和位移计算的处理。
见图1,用于测量位移的系统10,包括待测对象11,波发射器和回波接收装置(换能器)12,编码激发的信号产生器13,前置放大器、滤波器和时间增益补偿(TGC)装置14,模/数转换器ADC 15,PC计算机或微处理器16。功能块16可包括编码激励,TGC和ADC时序控制17,回波A和B的脉冲压缩18,主瓣去除19,杂波收集20和位移计算21。
波发射器和回波接收装置(换能器)12是向目标产生超声波(或雷达,声纳)并接收来自目标的回波的换能器。用于编码激励的信号发生器13被配置为产生GCS对从而控制换能器。前置放大器、滤波器和TGC装置14被配置为进行前置放大、滤波和对衰减产生的回波能量损耗进行补偿。ADC 15被配置为将模拟回波转换成数字回波。PC计算机或微处理器16被配置为控制编码激励的时序、TGC和ADC17,对回波A和B 18进行脉冲压缩,主瓣去除19,杂波累积20,以及检测位移21。
Golay互补序列是指一对序列A(a0,a1,…,aN-1)和B(b0,b1,…,bN-1)。A和B的长度均为N。A和B的元素是-1或1。在所描述的实施例中,GCS具有如下特征:
cj+dj=2N j=0
                    (1)
cj+dj=0  j≠0
其中cj和dj是GCS A和GCS B的自相关:
c j = Σ i = j i = N - 1 a i a i - j j = 0 , . . . , N - 1 Σ i = 0 i = N - 1 + j a i a i - j j = - N + 1 , . . . , - 1 - - - ( 2 )
d j = Σ i = j i = N - 1 b i b i - j j = 0 , . . . , N - 1 Σ i = 0 i = N - 1 + j b i b i - j j = - N + 1 , . . . , - 1 - - - ( 3 )
脉冲压缩是为了获得来自目标的陡峭(sharp)的回波所进行的数字信号处理环节。GCS A和GCS B被用来顺序控制换能器时,将来自同一物体的两个回波信号的自相关相加,就会获得一个非常陡峭的、通常称为压缩脉冲或主瓣的脉冲(见图2)。用于获得窄脉冲的自相关操作和相加操作被称为脉冲压缩。
经过脉冲压缩后可以获得陡峭的脉冲(也称为主瓣)。物体在GCS A和GCS B之间的时间间隔期间的位移在主瓣周围导致噪声(也称为杂波)并导致主瓣幅度减小。
图3显示了目标位移对脉冲压缩产生的影响。图3(a)显示了在GCS A和GCS B之间物体未发生移动时的压缩脉冲(两个自相关的相加)。在图3(b)、(c)和(d)中,物体分别在GCS A和GCS B之间移动了7.5、15和30微米。仿真中采用的换能器的中心频率f0和-3分贝的带宽η分别为10MHz和70%。GCS的长度N和宽度Tc分别为256和100纳秒。回波采样率fs为100MHz。
在2维超声图像中,杂波会在图像中导致图像失真。图4显示了未经过对数压缩的超声图像中的组织运动的影响。在图4(a)中,组织在GCS A和GCS B间没有移动。图4(b),(c)和(d)显示了当组织在GCS A和GCS B之间分别移动了7.5、15和30微米时回波的求和后的互相关的幅度。图像大小约为5毫米×5毫米。
图5显示了在GCS A和GCS B激励之间物体运动和压缩脉冲的信杂比(SCR)之间的关系。在GCS A和B激励之间的物体的移动与压缩脉冲的SCR是相关联的。仿真中采用的换能器中心频率f0为10MHz,-3分贝的带宽η为70%。GCS长度N和宽度Tc分别为256和100纳秒。回波采样率fs为100MHz。
当两次回波(分别来自GCS A和GCS B)之间存在一定的时移时,这两个回波的自相关和就会产生杂波。杂波幅度与时移成正比。时移越大,杂波越强。图6(a)显示了来自GCS A和GCS B的回波自相关的部分波形,分别由实线和虚线表示。虚线波形移动了GCS码的5/10。图6(b)显示了不同时移量(GCS码宽度的1/10、2/10、3/10、4/10和5/10)所产生的杂波。GCS长度N为256。X轴为自相关函数的采样点n。
在传统的脉冲-回波方法中,一般通过两个回波之间的时移得到位移。计算时移的精度受噪声水平、换能器带宽和ADC采样率限制。
在本发明中,杂波(在传统的方法中作为噪声对待)幅度,取代了传统方法中的时移,被用来获得位移。换能器带宽和ADC的采样率对杂波幅度没有直接影响。由于脉冲压缩去掉了大部分白噪声,噪声对杂波也几乎没有影响。所以杂波的幅度是比时移更加鲁棒性的检测物体位移的方法。
物体位移与利用Golay编码激励(见图5)的压缩脉冲的SCR(信杂比)相关。然而,在实际系统中,SCR由于采样精度和噪声的限制而并不是可靠的物体位移测量方法。因此,使用新的参数Sac来取代SCR来测量物体位移。Sac是所有杂波绝对值的和。一般主瓣宽度与单脉冲的宽度相等,因此,
S ac = Σ t = 0 t = L / 2 - l | P c ( t ) | + Σ t = L / 2 + l t = L | P c ( t ) |
其中,Pc(t)是压缩脉冲。L是整个Golay序列/码的长度。l是单脉冲的宽度。
图7显示了不同噪声水平下物体位移与Sac之间的关系,并且如果没有噪声,则位移测量分辨率几乎是无限的。白噪声只能对位移和Sac之间关系产生偏置。图7显示了对于回波中不同噪声水平的Sac-位移曲线。仿真中参数设定与图3相同。在位移等于d=±7.5微米、±15微米、…、±75微米时的Sac值由每条曲线中的星号标记出来。回波幅度被归一化为1。图7(a)显示无噪声情况。图7(b)~(f)中随机噪声均匀地分别在间隔(-0.2,0.2)、(-0.4,0.4)、(-0.6,0.6)、(-0.8,0.8)和(-1,1)之间。
图8显示了单个目标的位移测量算法,包括以下9个步骤:(1)回波A和B的脉冲压缩;(2)以范围从-10到10倍变化的采样间隔(步长为1)对回波A和B进行脉冲压缩;(3)计算杂波幅度绝对值的和;(4)找到杂波幅度绝对值和的最小值(P0);(5)对于回波A和B脉冲压缩进行10倍插值;(6)以偏置从-1(P1)到1(P2)倍的采样间隔(步长为0.1)对回波A和B进行脉冲压缩;(7)找到杂波幅度绝对值和的最小值(P’0);(8)对回波A和B的脉冲压缩进行100倍插值;(9)计算位移。
在单一目标的测量中,根据Sac-位移曲线中直接得到位移。位移测量的有效范围由Sac-位移曲线中的单调区间决定(图7中0~75微米)。可以通过将某一时移加到GCS B的回波从而将超出范围的位移调节到Sac-位移曲线的单调区间中。
基于上述算法(图8),用图1所示系统获取来自物体的实验数据并进行分析。实验中实际位移从1纳米到6微米。所有测量的误差的均值和STD为-5.76纳米±36.27纳米。在图9(a)中基于杂波信号的实际位移与超声测量的位移进行了比较。图9(b)为测量误差的直方图。
图10显示了对物体中期望区域进行位移测量的算法。该算法被分为三个部分:主瓣去除,匹配滤波和幅度归一化。
主瓣去除:主瓣去除的目的是把杂波信号从压缩脉冲中提取出来。如果GCS A和GCS B之间的时间间隔很短,则物体的期望区域中的位移可以被忽略。因此间隔时间很短的压缩脉冲只有主瓣,可以从具有长时间间隔的压缩脉冲中减去该间隔时间很短的压缩脉冲以去除主瓣。
解卷积:解卷积的目的是累积杂波能量并去除杂波混叠。采用预白化滤波器来提高解卷积的SNR。预白化的目的在于使信号包含所有频率的强度相等的分量。预白化滤波器使非白色信号转化成接近白色信号。这是通过预测器实现的。白化信号的方式是尝试根据前面样本的信息对样本n进行预测。如果我们把这一预测从实际样本n中减掉,则样本n只剩下了与其它样本不相关的样本部分。
预测器是使用AR模型设计的:
Figure A20078003081800151
ak是AR模型的系数。ak可以依照Yule-Walker方程计算:
R 0 R - 1 · · · R 1 - p R 1 R 0 R 2 - p · · · · · · · · · R p - 1 R p - 2 · · · R 0 a 1 a 2 · · · a p = R 1 R 2 · · · R p
预白化滤波器的传递函数是:
H ( z ) = 1 - Σ k = 1 p a k z - k
图11显示了GCS(N=256)的预白化滤波器的系数。经过预白化后,匹配滤波器随即被应用于白色信号以产生解卷信号。
幅度归一化:幅度归一化的目的是把来自不同物体/散射体的回波的幅度进行归一化。“散射体”是散射传输波的分界面,由此至少返回一些回波。“目标”是指在物体中所选的研究区域,可以包含一个或多个散射体。杂波的幅度与位移和回波幅度成比例。通过幅度归一化,杂波幅度仅与位移成比例。
幅度归一化的方法是:(1)采用Hilbert变换对匹配滤波器的输出(信号M)和具有主瓣的原始压缩脉冲(信号C)进行解调;(2)将解调的信号M除以信号C。
进行两个计算机仿真以对图10中的算法进行验证。
在第一仿真中,采用10MHz,90%带宽的换能器。发射频率为10MHz。采样率为100MHz。码长度是128。发射128位码后开始对回波进行采样。在第一码和最后码的发射之间的回波是不可用的。由此产生了一段短的、无法检测的区域,被称为死区。本仿真中的死区为9.6毫米。沿着超声波束放置十个散射体(图12)。如图12所示,散射体之间的距离分别是:1.5毫米、2.25毫米、3.0毫米、3.75毫米、4.5毫米、5.25毫米、6.0毫米、6.75毫米和7.5毫米。从左至右,散射子的位移分别为0.75微米、1.5微米、2.25微米、3微米、3.75微米、4.5微米、5.25微米、6微米、6.75微米和7.5微米。来自各散射子的回波幅度分别为1.1、1、0.9、0.8、0.7、0.6、0.5、0.4、0.3和0.2。
图13显示了各个步骤仿真的结果。图13分别示出了:(a)GCS A的回波,(b)GCS B的回波,(c)压缩脉冲,(d)主瓣去除后获得的杂波信号,(e)解卷积的结果,(f)幅度归一化和解调的结果。
图14比较了杂波的幅度(幅度归一化后)和实际位移。图14表明幅度归一化的结果与实际位移成正比。
在第二仿真中,散射子被均匀地沿着扫描线分布,且散射子之间的间距为225μm。来自散射子的回波幅度与散射子的深度成指数关系。散射子的位移具有正弦分布,如图15所示。GCS编码激励和超声换能器的参数设置与第一仿真相同。
图16显示了各步骤仿真的结果。图16分别显示了:(a)GCS A的回波,(b)GCS B的回波,(c)压缩脉冲,(d)主瓣去除后获得的杂波信号,(e)解卷积的结果,(f)幅度归一化的结果。
图17中,对实际轴向位移设置的绝对值和图16(f)中所示波形的包络进行了比较。结果表明,图10中描述的算法可以测量物体期望原因(desired reason)的很小的位移(±0.75微米)。
如以上描述以及相应的插图所示,对示意性的实施例可以在未偏离本发明保护范围的前提下做出各种修改,因此以上描述和插图中所包含的内容应视为示意性而非限制。故当前发明的广度和范围不受以上任何示意性实施例的限制,而只由以下各声明及其等价物所定义。

Claims (45)

1.一种用超声测量物体期望区域的位移的方法,包括:
发射两个互补的序列;
接收所述互补序列的回波;
对所述两个回波进行脉冲压缩;以及
根据压缩脉冲输出的主瓣周围的剩余杂波信号来获得两个发射之间的散射体的位移。
2.如权利要求1的方法,所发射的信号是一对Golay互补序列。
3.如权利要求1的方法,进一步包含自相关所述互补序列的步骤。
4.如权利要求1中的方法,其中根据所获得的位移获得以下至少其中之一:厚度,速度,应变,弹性刚度,和粘性阻尼。
5.如权利要求1的方法,其中通过散射体的运动产生所获得的位移。
6.如权利要求1的方法,其中回波是来自散射体的反射或散射波。
7.如权利要求1的方法,其中脉冲压缩包含回波GCS A(B)与GCS码A(B)的相关以及这两个相关的和。
8.如权利要求1方法,其中主瓣是压缩脉冲的主峰。
9.如权利要求1方法,其中剩余杂波信号是由两个回波的时移中的一个或由在两个GCS发射期间的散射体位移所产生的。
10如权利要求1的方法,进一步包含去除主瓣的步骤。
11.如权利要求9的方法,进一步包含杂波累积的步骤。
12.如权利要求11中的方法,其中杂波累积是累积杂波信号中的能量。
13.如权利要求11中的方法,其中使用预白化滤波器进行杂波累积。
14.如权利要求11中的方法,进一步包含匹配滤波以累积杂波信号的能量的步骤。
15.如权利要求1的方法,其中位移检测是根据杂波信号的幅度和时移之间的关系获得的。
16.如权利要求16的方法,其中通过重复插值来提高杂波信号的幅度和时移之间的关系的分辨率。
17.如权利要求1的方法,其中通过幅度归一化获得位移检测。
18.如权利要求17的方法,其中通过去除来自不同散射体的回波的幅度的影响来进行幅度归一化。
19.一种位移测量装置,用于测量选定物体中期望目标的位移,所述装置包括:
波发射器,被配置用于向所述目标发射能量;
接收器,被配置用于接收从所述目标反射的信号;
处理器,所述处理器被配置用于向所述波发射器提供两个互补的序列,以及所述处理器被配置用于根据来自所述两个互补的序列的反射或散射信号的主瓣周围的剩余杂波信号来测量目标的位移。
20.如权利要求19的装置,其中所述两个互补的序列是Golay互补序列。
21.如权利要求19的装置,进一步包含以下至少其中之一:信号产生器,前置放大器,滤波器,时间增益补偿电路和模拟数字转换电路,其与所述发射器或所述接收器中的至少一个操作地通信。
22.如权利要求21的装置,其中前置放大器、滤波器和时间增益补偿电路被配置用来放大回波,过滤噪声和补偿回波的能量损失。
23.如权利要求19的装置,其中所述处理器被配置用来控制换能器/接收器、前置放大器、滤波器、时间增益补偿电路和模拟数字转换电路的时序,并计算物体或物体的期望区域的位移。
24.如权利要求20的装置,其中所发射的信号是一对Golay互补序列A和B或者其他互补序列,以及互补序列A和B的自相关的和是德尔塔函数。
25.如权利要求19的装置,其中所述处理器进一步被配置为获得以下至少其中之一:厚度、速度、应变、弹性刚度和粘性阻尼。
26.如权利要求19的装置,其中所述处理器进一步被配置为获得由目标运动而产生的维度参数。
27.如权利要求19的装置,其中位移计算包括分析压缩脉冲输出的主瓣周围的剩余杂波信号。
28.如权利要求19的装置,进一步包含脉冲压缩。
29.如权利要求28的装置,其中所述脉冲压缩包含至少两个回波A(B)与GCS码A(B)之间的相关以及所述相关之和。
30.如权利要求19的装置,其中主瓣是压缩脉冲的主峰。
31.如权利要求19的装置,其中杂波是由两回波的时移产生的信号。
32.如权利要求19的装置,其中杂波是在两Golay互补序列发射期间的目标位移所产生的信号。
33.如权利要求19的装置,其中杂波分析包括主瓣去除,杂波累积,以及位移检测。
34.如权利要求19的装置,其中通过去除主瓣获得杂波信号。
35.如权利要求19的装置,其中进一步包括杂波累积。
36.如权利要求35的装置,其中杂波积聚包括预白化滤波器。
37.如权利要求35的装置,进一步包括匹配滤波以累积杂波信号的能量。
38.如权利要求33的装置,其中位移检测来自于累积的杂波信号。
39.如权利要求19的装置,其中位移检测来自于时移与杂波的关系。
40.如权利要求39的装置,其中通过重复插值来提高杂波信号的幅度与时移之间的关系的分辨率。
41.如权利要求19的装置,其中通过幅度归一化进行位移检测。
42.如权利要求41的装置,其中幅度归一化去除了来自不同物体或散射体的回波的幅度的影响。
43.如权利要求1的方法,进一步包括在脉冲压缩前时移GCS B的回波。
44.如权利要求1的方法,其中所述发射是雷达发射和声纳发射之一。
45.如权利要求19的装置,其中所述发射是雷达发射和声纳发射之一。
CNA2007800308186A 2006-06-20 2007-06-20 高分辨率超声位移测量的装置与方法 Pending CN101595396A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US11/471,129 2006-06-20
US11/471,129 US7535797B2 (en) 2006-06-20 2006-06-20 High-resolution ultrasound displacement measurement apparatus and method

Publications (1)

Publication Number Publication Date
CN101595396A true CN101595396A (zh) 2009-12-02

Family

ID=38834350

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2007800308186A Pending CN101595396A (zh) 2006-06-20 2007-06-20 高分辨率超声位移测量的装置与方法

Country Status (3)

Country Link
US (1) US7535797B2 (zh)
CN (1) CN101595396A (zh)
WO (1) WO2007149918A2 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103123265A (zh) * 2011-11-21 2013-05-29 柯宜京 一种提高超声测量或检测仪器轴向分辨力的方法和装置
CN103908304A (zh) * 2014-03-14 2014-07-09 中瑞科技(常州)有限公司 一种超声弹性成像系统
CN109343044A (zh) * 2018-10-11 2019-02-15 中国人民解放军国防科技大学 一种提高格雷互补波形多普勒分辨率的方法

Families Citing this family (54)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7535797B2 (en) * 2006-06-20 2009-05-19 Rehabtek High-resolution ultrasound displacement measurement apparatus and method
US9867530B2 (en) 2006-08-14 2018-01-16 Volcano Corporation Telescopic side port catheter device with imaging system and method for accessing side branch occlusions
US9596993B2 (en) 2007-07-12 2017-03-21 Volcano Corporation Automatic calibration systems and methods of use
WO2009009799A1 (en) 2007-07-12 2009-01-15 Volcano Corporation Catheter for in vivo imaging
WO2009009802A1 (en) 2007-07-12 2009-01-15 Volcano Corporation Oct-ivus catheter for concurrent luminal imaging
US11141063B2 (en) 2010-12-23 2021-10-12 Philips Image Guided Therapy Corporation Integrated system architectures and methods of use
US11040140B2 (en) 2010-12-31 2021-06-22 Philips Image Guided Therapy Corporation Deep vein thrombosis therapeutic methods
WO2013033489A1 (en) 2011-08-31 2013-03-07 Volcano Corporation Optical rotary joint and methods of use
US9420997B2 (en) * 2012-06-14 2016-08-23 Siemens Medical Solutions Usa, Inc. Motion artifact suppression in ultrasound diagnostic imaging
US11272845B2 (en) 2012-10-05 2022-03-15 Philips Image Guided Therapy Corporation System and method for instant and automatic border detection
US9292918B2 (en) 2012-10-05 2016-03-22 Volcano Corporation Methods and systems for transforming luminal images
US9324141B2 (en) 2012-10-05 2016-04-26 Volcano Corporation Removal of A-scan streaking artifact
US9286673B2 (en) 2012-10-05 2016-03-15 Volcano Corporation Systems for correcting distortions in a medical image and methods of use thereof
US10568586B2 (en) 2012-10-05 2020-02-25 Volcano Corporation Systems for indicating parameters in an imaging data set and methods of use
US9367965B2 (en) 2012-10-05 2016-06-14 Volcano Corporation Systems and methods for generating images of tissue
US9307926B2 (en) 2012-10-05 2016-04-12 Volcano Corporation Automatic stent detection
JP2015532536A (ja) 2012-10-05 2015-11-09 デイビッド ウェルフォード, 光を増幅するためのシステムおよび方法
US9858668B2 (en) 2012-10-05 2018-01-02 Volcano Corporation Guidewire artifact removal in images
US10070827B2 (en) 2012-10-05 2018-09-11 Volcano Corporation Automatic image playback
US9840734B2 (en) 2012-10-22 2017-12-12 Raindance Technologies, Inc. Methods for analyzing DNA
EP2931132B1 (en) 2012-12-13 2023-07-05 Philips Image Guided Therapy Corporation System for targeted cannulation
US10942022B2 (en) 2012-12-20 2021-03-09 Philips Image Guided Therapy Corporation Manual calibration of imaging system
US10939826B2 (en) 2012-12-20 2021-03-09 Philips Image Guided Therapy Corporation Aspirating and removing biological material
WO2014099899A1 (en) 2012-12-20 2014-06-26 Jeremy Stigall Smooth transition catheters
US11406498B2 (en) 2012-12-20 2022-08-09 Philips Image Guided Therapy Corporation Implant delivery system and implants
US9709379B2 (en) 2012-12-20 2017-07-18 Volcano Corporation Optical coherence tomography system that is reconfigurable between different imaging modes
JP2016506276A (ja) 2012-12-20 2016-03-03 ジェレミー スティガール, 血管内画像の位置の特定
US9486143B2 (en) 2012-12-21 2016-11-08 Volcano Corporation Intravascular forward imaging device
US10993694B2 (en) 2012-12-21 2021-05-04 Philips Image Guided Therapy Corporation Rotational ultrasound imaging catheter with extended catheter body telescope
US9383263B2 (en) 2012-12-21 2016-07-05 Volcano Corporation Systems and methods for narrowing a wavelength emission of light
US10191220B2 (en) 2012-12-21 2019-01-29 Volcano Corporation Power-efficient optical circuit
CA2895993A1 (en) 2012-12-21 2014-06-26 Jason Spencer System and method for graphical processing of medical data
WO2014099672A1 (en) 2012-12-21 2014-06-26 Andrew Hancock System and method for multipath processing of image signals
US9612105B2 (en) 2012-12-21 2017-04-04 Volcano Corporation Polarization sensitive optical coherence tomography system
CN103006274A (zh) * 2012-12-21 2013-04-03 深圳大学 一种超声检测角膜粘弹性的方法和系统
US10058284B2 (en) 2012-12-21 2018-08-28 Volcano Corporation Simultaneous imaging, monitoring, and therapy
US10413317B2 (en) 2012-12-21 2019-09-17 Volcano Corporation System and method for catheter steering and operation
US10166003B2 (en) 2012-12-21 2019-01-01 Volcano Corporation Ultrasound imaging with variable line density
JP6243453B2 (ja) 2013-03-07 2017-12-06 ボルケーノ コーポレイション 血管内画像におけるマルチモーダルセグメンテーション
US10226597B2 (en) 2013-03-07 2019-03-12 Volcano Corporation Guidewire with centering mechanism
CN105228518B (zh) 2013-03-12 2018-10-09 火山公司 用于诊断冠状微脉管疾病的系统和方法
US20140276923A1 (en) 2013-03-12 2014-09-18 Volcano Corporation Vibrating catheter and methods of use
US11026591B2 (en) 2013-03-13 2021-06-08 Philips Image Guided Therapy Corporation Intravascular pressure sensor calibration
US9301687B2 (en) 2013-03-13 2016-04-05 Volcano Corporation System and method for OCT depth calibration
US10758207B2 (en) 2013-03-13 2020-09-01 Philips Image Guided Therapy Corporation Systems and methods for producing an image from a rotational intravascular ultrasound device
US10292677B2 (en) 2013-03-14 2019-05-21 Volcano Corporation Endoluminal filter having enhanced echogenic properties
US20160030151A1 (en) 2013-03-14 2016-02-04 Volcano Corporation Filters with echogenic characteristics
US10219887B2 (en) 2013-03-14 2019-03-05 Volcano Corporation Filters with echogenic characteristics
ES2544581B1 (es) * 2014-01-31 2016-06-08 Vicente Diaz Fuente Método y aparato para la detección de distancia, velocidad y dirección de objetos
EP3268768B1 (en) * 2015-03-11 2023-08-02 Raytheon Company Wideband channel equalization for signals propagated in lossy transmission media
WO2017040267A1 (en) * 2015-08-28 2017-03-09 Soneter, Inc. Flow meter configuration and calibration
US11116415B2 (en) * 2017-06-07 2021-09-14 Bragi GmbH Use of body-worn radar for biometric measurements, contextual awareness and identification
EP3435116A1 (en) 2017-07-24 2019-01-30 Koninklijke Philips N.V. An ultrasound probe and processing method
IT201900008994A1 (it) * 2019-06-14 2020-12-14 Aresys S R L Sistema radar per misurare lo spessore di un oggetto e metodo per misurare lo spessore di un oggetto.

Family Cites Families (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4862177A (en) * 1974-11-25 1989-08-29 Hughes Aircraft Company Processor for discriminating between ground and moving targets
US4257275A (en) * 1977-06-10 1981-03-24 Yokogawa Electric Works, Ltd. Velocity detecting apparatus
JPS58191504A (ja) 1982-05-04 1983-11-08 Mitsubishi Electric Corp 受信装置
JPS60169782A (ja) * 1984-02-14 1985-09-03 Nec Corp 移動目標表示装置
US5014712A (en) 1989-12-26 1991-05-14 General Electric Company Coded excitation for transmission dynamic focusing of vibratory energy beam
US5457818A (en) * 1990-06-08 1995-10-10 Butler; James A. Detection threshold adjustment method for linear matched filter receivers
US6809681B1 (en) * 1992-08-25 2004-10-26 Raytheon Company Random-modulation radar signal-induced interference cancellation method and apparatus
US5738635A (en) 1993-01-22 1998-04-14 Technomed Medical Systems Adjustable focusing therapeutic apparatus with no secondary focusing
US5573497A (en) 1994-11-30 1996-11-12 Technomed Medical Systems And Institut National High-intensity ultrasound therapy method and apparatus with controlled cavitation effect and reduced side lobes
US6009046A (en) 1995-03-02 1999-12-28 Acuson Corporation Ultrasonic harmonic imaging system and method
US6193663B1 (en) 1997-12-18 2001-02-27 Acuson Corporation Diagnostic ultrasound imaging method and system with improved frame rate
US5964706A (en) 1998-03-18 1999-10-12 General Electric Company Method and apparatus for pulsed doppler imaging using coded excitation on transmit and pulse compression on receive
US6095977A (en) 1998-03-26 2000-08-01 Hall; Anne Lindsay Method and apparatus for color flow imaging using Golay-coded excitation on transmit and pulse compression on receive
US5938611A (en) 1998-03-26 1999-08-17 General Electric Company Method and apparatus for color flow imaging using coded excitation with single codes
US6186949B1 (en) 1998-03-31 2001-02-13 General Electric Company Method and apparatus for three-dimensional flow imaging using coded excitation
US6210332B1 (en) 1998-03-31 2001-04-03 General Electric Company Method and apparatus for flow imaging using coded excitation
US6312384B1 (en) 1998-03-31 2001-11-06 General Electric Company Method and apparatus for flow imaging using golay codes
US5984869A (en) 1998-04-20 1999-11-16 General Electric Company Method and apparatus for ultrasonic beamforming using golay-coded excitation
US6050947A (en) * 1998-04-20 2000-04-18 General Electric Company Method and apparatus for harmonic tissue imaging and contrast imaging using coded transmission
US5961463A (en) 1998-08-24 1999-10-05 General Electric Company Nonlinear imaging using orthogonal transmit and receive codes
US6113545A (en) 1998-04-20 2000-09-05 General Electric Company Ultrasonic beamforming with improved signal-to-noise ratio using orthogonal complementary sets
US6239741B1 (en) 1998-07-20 2001-05-29 Multispectral Solutions, Inc. UWB dual tunnel diode detector for object detection, measurement, or avoidance
US6123670A (en) 1998-12-15 2000-09-26 General Electric Company Ultrasound imaging with optimal image quality in region of interest
US6010456A (en) 1998-12-30 2000-01-04 General Electric Company Method and apparatus for acoustic subtraction imaging using linear and nonlinear ultrasonic images
US6146328A (en) 1999-02-23 2000-11-14 General Electric Company Method and apparatus for ultrasonic beamforming using golay-coded excitation
US6155980A (en) 1999-03-16 2000-12-05 General Electric Company Ultrasonic imaging system with beamforming using unipolar or bipolar coded excitation
US6213947B1 (en) 1999-03-31 2001-04-10 Acuson Corporation Medical diagnostic ultrasonic imaging system using coded transmit pulses
US6241674B1 (en) 1999-03-31 2001-06-05 Acuson Corporation Medical ultrasound diagnostic imaging method and system with nonlinear phase modulation pulse compression
US6179780B1 (en) 1999-08-06 2001-01-30 Acuson Corporation Method and apparatus for medical diagnostic ultrasound real-time 3-D transmitting and imaging
US6056693A (en) 1999-08-16 2000-05-02 General Electric Company Ultrasound imaging with synthetic transmit focusing
US6375618B1 (en) 2000-01-31 2002-04-23 General Electric Company Enhanced tissue-generated harmonic imaging using coded excitation
US6491631B2 (en) 2001-01-11 2002-12-10 General Electric Company Harmonic golay-coded excitation with differential pulsing for diagnostic ultrasound imaging
US6760486B1 (en) 2000-03-28 2004-07-06 General Electric Company Flash artifact suppression in two-dimensional ultrasound imaging
US6790182B2 (en) 2000-09-05 2004-09-14 Koninklijke Philips Electronics N.V. Ultrasound system and ultrasound diagnostic apparatus for imaging scatterers in a medium
US6487433B2 (en) 2001-01-08 2002-11-26 General Electric Company Method and apparatus using golay-coded excitation for echocardiology
JP3800966B2 (ja) * 2001-02-13 2006-07-26 東海ゴム工業株式会社 空気圧制御式制振器
US6666825B2 (en) 2001-07-05 2003-12-23 General Electric Company Ultrasound transducer for improving resolution in imaging system
US6796944B2 (en) 2002-05-17 2004-09-28 Ge Medical Systems Global Technology, Llc Display for subtraction imaging techniques
US7094204B2 (en) 2002-08-23 2006-08-22 Siemens Medical Solutions Usa, Inc. Coded excitation imaging for use with bipolar, unipolar and other waveforms
US20050096544A1 (en) 2003-10-30 2005-05-05 Xiaohui Hao Method and apparatus for single transmission Golay coded excitation
US9310475B2 (en) 2003-11-21 2016-04-12 General Electric Company Method and apparatus for transmitting multiple beams
US20060074320A1 (en) 2004-08-27 2006-04-06 Yoo Yang M Home ultrasound system
US7535797B2 (en) * 2006-06-20 2009-05-19 Rehabtek High-resolution ultrasound displacement measurement apparatus and method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103123265A (zh) * 2011-11-21 2013-05-29 柯宜京 一种提高超声测量或检测仪器轴向分辨力的方法和装置
CN103123265B (zh) * 2011-11-21 2016-11-02 柯宜京 一种提高超声测量或检测仪器轴向分辨力的方法和装置
CN103908304A (zh) * 2014-03-14 2014-07-09 中瑞科技(常州)有限公司 一种超声弹性成像系统
CN109343044A (zh) * 2018-10-11 2019-02-15 中国人民解放军国防科技大学 一种提高格雷互补波形多普勒分辨率的方法

Also Published As

Publication number Publication date
US20070291591A1 (en) 2007-12-20
WO2007149918A2 (en) 2007-12-27
WO2007149918A3 (en) 2008-12-04
US7535797B2 (en) 2009-05-19

Similar Documents

Publication Publication Date Title
CN101595396A (zh) 高分辨率超声位移测量的装置与方法
O'Donnell et al. Internal displacement and strain imaging using ultrasonic speckle tracking
US8062223B2 (en) Using pulsed-wave ultrasonography for determining an aliasing-free radial velocity spectrum of matter moving in a region
CN103269639B (zh) 利用质心估计剪切波速度
CN1022135C (zh) 傅里叶处理脉冲多普勒信号的方法和装置
CN105997148B (zh) 脉冲多普勒超高谱分辨率成像处理方法及处理系统
JPH0246213B2 (zh)
KR100381874B1 (ko) 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
CN110954885B (zh) 基于sbl的捷变频雷达自适应目标重构方法
CN105844645A (zh) 一种瞬时弹性成像编码检测方法及系统
US7363177B2 (en) Apparatus and method for performing the time delay estimation of signals propagating through an environment
US20200163649A1 (en) Estimating Phase Velocity Dispersion in Ultrasound Elastography Using a Multiple Signal Classification
Wilhjelm et al. Target velocity estimation with FM and PW echo ranging Doppler systems I. Signal analysis
CN106022369A (zh) 一种基于超声背向散射信号成像的骨质评价方法
CN103690194A (zh) 基于宽带信号的血流流速估计方法和系统
US6509866B2 (en) Fast chirp transform
US20220341885A1 (en) Free-encoder positioning system using acoustic features and imu
CN105676225B (zh) 不透明浑浊液体中的超声测距系统和方法
EP2366997B1 (en) Method and device for determining the structural organization of an object with ultrasounds
Murino et al. A confidence-based approach to enhancing underwater acoustic image formation
JP3881078B2 (ja) 周波数推定方法、周波数推定装置、ドップラソナーおよび潮流計
Leetang et al. Evaluation of ultrasonic target detection by alternate transmission of different codes in M-sequence pulse compression
US9918698B2 (en) System and method for gradient-based k-space search for shear wave velocity dispersion estimation
Zedel Noise in coherent Doppler sonar velocity measurements: Where does it come from and what can you do about it?
Martin et al. Range discrimination in ultrasonic vibrometry: Theory and experiment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20091202