CN1886094A - 多普勒速度检测装置以及使用该装置的超声波诊断装置 - Google Patents
多普勒速度检测装置以及使用该装置的超声波诊断装置 Download PDFInfo
- Publication number
- CN1886094A CN1886094A CNA2004800349981A CN200480034998A CN1886094A CN 1886094 A CN1886094 A CN 1886094A CN A2004800349981 A CNA2004800349981 A CN A2004800349981A CN 200480034998 A CN200480034998 A CN 200480034998A CN 1886094 A CN1886094 A CN 1886094A
- Authority
- CN
- China
- Prior art keywords
- expansion coefficient
- signal
- time
- plural
- velocity
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 54
- 238000004458 analytical method Methods 0.000 claims abstract description 64
- 230000005540 biological transmission Effects 0.000 claims abstract description 15
- 230000017531 blood circulation Effects 0.000 claims description 41
- 238000002604 ultrasonography Methods 0.000 claims description 36
- 230000007246 mechanism Effects 0.000 claims description 32
- 230000011514 reflex Effects 0.000 claims description 25
- 230000033001 locomotion Effects 0.000 claims description 23
- 239000000203 mixture Substances 0.000 claims description 21
- 239000000523 sample Substances 0.000 claims description 19
- 210000001835 viscera Anatomy 0.000 claims description 13
- 230000003068 static effect Effects 0.000 claims description 8
- 238000000034 method Methods 0.000 abstract description 11
- 230000006870 function Effects 0.000 description 54
- 230000014509 gene expression Effects 0.000 description 31
- 238000012545 processing Methods 0.000 description 15
- 238000012360 testing method Methods 0.000 description 15
- 238000010606 normalization Methods 0.000 description 11
- 230000004044 response Effects 0.000 description 11
- 238000001914 filtration Methods 0.000 description 8
- 238000003745 diagnosis Methods 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 238000012937 correction Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 229940126655 NDI-034858 Drugs 0.000 description 2
- 241000290929 Nimbus Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 210000003494 hepatocyte Anatomy 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 101100063942 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) dot-1 gene Proteins 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
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
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8981—Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
-
- 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/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种多普勒速度检测装置,其特征在于,具备:对受检体多次收发脉冲波的机构、以接收回波信号为基础分析与所述受检体内的运动反射物有关的速度的速度分析机构,所述速度分析机构以下述方式构成,即,以虚数单位为系数,当将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,从各复数展开系数的大小和各复数展开系数间的大小的比得到与所述受检体内的运动反射物有关的带符号的速度信号。由此,本发明实现在很好地对杂波信号进行严加区别的同时可以进行已抑制了误差的速度检测/分析的多普勒速度检测技术,另外,还提供使用了该技术的超声波诊断装置。
Description
技术领域
本发明涉及通过对生物体收发超声波检测出生物体内部的血流等运动的反射物来显示速度的分布或速度的空间分布的医疗诊断用的超声波多普勒血流计或血流描绘装置等超声波诊断技术。
背景技术
在脉冲/多普勒装置中,向对象物多次收发脉冲状的超声波或电磁波等波,对于已接收的图1那样的多个回波信号,将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的着眼时刻的信号而成时间序列信号,通过对该时间序列信号进行分析,得到与运动的物体有关的速度等信息。作为最普通的信号处理方法,是对这些接收信号进行正交检波,作为复数时间序列信号而进行分析。这样的脉冲/多普勒装置被广泛用作对生物体内部的血流等进行检测/描绘的超声波诊断装置、对雨云等进行检测/描绘的气象雷达、对飞行物体进行检测的航空雷达等。
在反射物只存在一个的情况下,可以从上述的时间序列信号的相位旋转速度λΔΦ/Δt即带符号的角频率,按照下式容易地求出反射物接近或远离收发两用机的运送速度v。
v=λΔΦ/Δt/2π (1)
在这里,λ是波长,Δt是脉冲发送的时间间隔,ΔΦ是相位旋转角。在接近的情况下,v的符号为正,在远离的情况下,v的符号为负。
然而,在如上所述的实际的脉冲/多普勒装置的接收回波中,基于静止反射物的回波信号或所谓杂波信号的强度比基于目标运动反射物的回波信号的强度通常大几个数量级以上。为此,只通过单纯将(1)式的处理应用于上述的时间序列信号,无法检测出运动反射物。因此,实际的脉冲/多普勒装置抑制基于静止反射物的回波信号,在进行相对强调基于运动反射物的回波信号的被称之为MTI(Moving Target Indicator)处理的信号处理之后,进行速度的检测或分析。
另外,在IEEE Transactions on Ultrasonics,Ferroelectrics,and FrequencyControl,Vol.42,927-937页(1995年)中,提出有多项回归滤波器(Polynomial Regression Filter)。
这是通过在时间序列信号上顺次最小二乘拟合0次式(常数)、1次式、……、M次式,并消除拟合成分,除去原来的时间序列信号所具有的漂移(drift)成分的处理。
作为MTI处理,最了解的处理是使用在时域以卷积积和表现的通常的低频截止滤波器的处理。该处理有下述的缺点。
1)当使用输入N1个时间序列数据点的低频截止滤波器时,输入到后段的速度检测/分析处理部的数据点损耗(N1-1)个。
2)难以获得锐截止特性的滤波器
实际上,基于静止反射物的回波信号在时间轴上也并不是完全静止的,当中途的介质摇动或是医疗用超声波诊断装置时,因为血流以外的静止脏器的缓慢运动,进行漂移。因此,基于静止反射物的回波信号的频率成分不只是相位旋转速度为零的直流成分,还具有相位旋转速度不为零的低频成分,在保持基于运动反射物的回波信号的同时,需要将其除去的截止特性,上述2)成为问题。当从通过N次的收发信号得到的N点的时间序列信号进行速度检测/分析处理时,上述1)是指,如果使用该种类的滤波器,实际能用于速度检测/分析运算的数据点减小至(N-N1+1)。这在像医疗诊断用超声波血流描绘装置那样实时性非常重要的应用中,并不优选。
作为解决该问题的处理,提出有上述的Polynomial Regression Filter(多项回归滤波器)。该处理是通过N点的输入时间序列信号乘以N×N的矩阵来表示。因此,作为输出信号,得到相同N点的时间序列信号,
1)没有数据点的损耗。
另外,其低频截止特性基于拟合的最大次数M的大小,
2)与具有同等的截止频率的上述以往型滤波器相比,具有相当大的锐截止特性。
不过,在通过多项式拟合/滤波器的处理的后段,当通过用上述(1)式表示的种类的以往方式,进行运动物体的速度的检测/分析时,问题是在截止频率的附近,速度运算误差比以往型滤波器大。为此,当对作为检测/分析结果的速度的精度有某种程度的要求时,通过多项式拟合/滤波器的锐截止特性,结果无法使用位于已延伸至低通的部分的频带的运动物体的信号。即,此时无法充分发挥上述2)的优势。
发明内容
因此,本发明鉴于现有的脉冲/多普勒装置具有的如上所述问题,实现在很好地对杂波信号进行严加区别的同时可以进行已抑制了误差的速度检测/分析的多普勒速度检测技术,另外,还提供使用了该技术的超声波诊断装置。
本发明的目的在于从本质上解决上述课题,首先,关于截止杂波信号的MTI滤波器的功能和其必要性,试着进行深入考察。作为具体例子,考虑从N点的时间序列复数信号进行速度分析的情况。这样的作为脉冲/多普勒速度分析而使用的典型的时间序列复数信号,接收的回波信号乘以90°相位的不同的2个输送频率信号而得到一组信号,以虚数单位为系数线性结合上述得到的一组信号并作为复数信号,按照发送时刻的顺序排列该信号的以各脉冲发送时刻为基准的时相相等的部分而得到。
在这样的离散时间序列信号所具有的频率或相位的分析中的常用方法是离散性傅立叶展开。具体而言,用表示频率的指数k=-N/,……,0,……,N/2的复数正弦波系展开为排列成N=1,……,N的时间序列信号。
Cs(n,k)=cos[kπ(2n-N-1)/N]+jsin[kπ(2n-N-1)/N]
………(2)
j是虚数单位。负的k对应于远离方向上的速度,正的k对应于接近方向上的速度。K=±N/2是乃奎斯特极限,因此,没有正负区别,所以由(2)表示的(N+1)个复数正弦波函数中独立的函数为N个,它们形成了正交函数。
假设反射物只存在1个且其在乃奎斯特极限范围的一定速度内运动的情况,作为运动速度的函数,计算通过上述展开得到的离散傅立叶展开系数,关于N=8的情况,对其绝对值绘图如图2所示。横轴是以脉冲发送的时间间隔的倒数PRF(pulse Repetition Frequency)的2π倍为单位表示与运动速度对应的相位旋转速度。在图2中,重叠表示展开系数的8个函数进行图示,其中,通过区别线在图3和图4中表示其中的各4个。
通过发送时间间隔的N倍的时间,以相位进行整数次旋转那样的运动速度,不是零的展开系数只有1个,其他的展开系数全都为零,非常好地进行速度分析。不过,以除此之外的相位旋转不是整数次的一般运动速度,各自的展开系数即使在距离横轴上峰值相当远的运动速度中,也具有峰值的1/10左右的大小的绝对值。这意味着作为速度分析计具有-20dB左右的远程串扰。当振幅比作为检测/分析对向的运动反射物回波更大的杂波信号进行漂移时,其出现严重问题。
该问题是,通过应用基于海明(Hanning)函数那样缓慢上升下降的窗口函数的权重,可以很好地进行抑制,但不能从本质上解决。另外,该方法与大幅度损耗时间序列信号的数量一样,是出现不需要的效果、尤其是不适于时间序列信号的数量原本较少的情况下的方法。
作为具体例,有作为以乃奎斯特极限的2/3的恒定速度运动的检测/分析对向的运动反射物,考虑具有其1000倍的回波信号的振幅的杂波信号以乃奎斯特极限速度的1/100进行偏移的情况。在图5中,表示作为此时的上述速度分析计的输出的速度频谱。作为检测/分析对象的运动反射物的波谱,被来自在零速度附近具有峰值的杂波信号的串扰成分所覆盖,在其输出波谱中不能检测出来。
这是需要截止杂波信号的MTI滤波器的最本质原因。作为MTI而被众所周知的、在时域内以卷积积和表示的通常的低频截止滤波器,在图6中用实线表示其低通振幅特性的典型例子。当输入这样的以n1=1,……,N+M-1排列的时间序列信号S1(n1)时,作为输出,得到以n0=1,……,N排列的时间序列信号S0(n0)。
S0(n0)=∑F(m)S1(n0+m-1)
………………………(3)
在这里,S表示m=1,……M的和。在图6中,用实线表示作为低频截止滤波器的最简单的M=3、F(1)=-1、F(2)=2、F(3)=-1的情况。在与图5相同的输入信号通过该低频截止滤波器之后,与图5一样进行速度分析,此时,如图7中的实线所示。通过低频截止滤波器的偏移将杂波信号的振幅抑制在1/2000倍左右,其结果,可以看到在乃奎斯特极限的0.3~0.4倍的速度处具有峰值的运动反射物的波谱。
作为相同目的的MTI滤波器,在图6中用虚线表示设计上述多项式拟合/滤波器时的通过振幅特性的典型例子。在该例中,对0~3次式进行拟合并减去。比上述的卷积滤波器更能看到锐低频截止特性。另外,在与图5相同的输入信号通过该低频截止滤波器之后,用虚线将与图5一样进行速度分析的结果显示于图7中。由于低频截止特性的锐度,可以完全抑制杂波信号。多项式拟合/滤波器特有的锐低频截止特性可以理解成如下所示。
关于狭义勒让德多项式,例如,在“岩波数学公式III”(特殊函数,82~85页)中有记述。当对其普遍化时,偶数次(2n)勒让德多项式可以定义成与次数比2n低的所有勒让德多项式正交的最高次数为2n的偶函数,奇数次(2n+1)勒让德多项式可以定义成与次数比2n+1低的所有勒让德多项式正交的最高次数为2n+1的奇函数。作为例子,关于由一系列的8点构成的时间序列信号,计算0~7次的离散勒让德函数,将偶数次显示于图8,将奇数次显示于图9。用实线表示0次、1次,用虚线表示2次、3次,用点划线表示4次、5次,用点线表示6次、7次。对振幅进行归一化并使RMS(均方根)值为1。将符号规定成最低次项的符号为正,即对于偶数次多项式,为常数项,对于奇数次多项式,一次项的符号为正。
图10是其频谱的绘图,用各个勒让德展开系数的最大值对频率成分的振幅的绝对值进行归一化并显示。在图10中,重叠8个函数并图示,但其中,通过区别线在图11和图12中表示其中的各4个。图中,用实线表示0次、7次的勒让德展开系数A(0)、A(7),用虚线表示1次、6次的勒让德展开系数A(1)、A(6),用点划线表示2次、5次的勒让德展开系数A(2)、A(5),用点线表示3次、4次的勒让德展开系数A(3)、A(4)。
与图2~4一样,图10~12表示针对输入复数正弦波的频率响应。将峰值视为主要响应,将其以外的部分视为次要响应。勒让德函数与图2~4的傅立叶展开系相比,高频侧的次要响应大于主要响应,但在低频侧没有次要响应,次数越大,越在频率原点附近急剧收敛成零。在这里,虽然省略分析证明,但作为频率的函数的m次勒让德函数,在频率原点处从0次到m-1次导数全部成为0。
上述的多项式拟合/滤波器是以其离散的勒让德函数为基础来展开输入信号,除外0~M次,将M+1次以上的成分的和作为输出信号的滤波器。因此,输出信号以在频率原点附近与频率的M+1次方成比例收敛成零。像在图6和图7中用虚线表示的例子那样,多项式拟合/滤波器的低频截止特性出色的信号源自于这样的勒让德函数的本质特性。如果进行如此考虑,其结论是,利用多项式拟合/滤波器或勒让德函数展开的方式,在本质上适合于抑制称之为杂波的在频率原点附近具有非常大的振幅的干扰成分的影响而检测目标运动物体,或者进行必须分析其速度的多普勒速度检测。
众所周知,在傅立叶展开系中,通过以复数正弦波为基础,可以像图2~图5或图7所示的那样进行带符号的速度分析,这样的分析方法得到广泛使用。然而,勒让德展开系尽管如上所述具有极其适合于多普勒速度检测的性质,但如果直接使用,从图10~12可知,绝对值虽相等,但对符号不同的速度显示出相等的响应,不能进行带符号的速度分析。
因此,在本发明中,以虚数单位j为系数,线性结合次数只差一个的偶数次和奇数次的勒让德函数,形成复数勒让德函数,由此可以进行带符号的速度分析。注意到偶数次的勒让德函数显示出类似余弦波的变化,奇数次的勒让德函数显示出类似正弦波的变化,通过以虚数单位j为系数线性结合余弦函数和正弦函数,可以得到相位角的j倍的指数函数、即复数正弦函数,可以与从该相位角的增减得到带符号的速度一样,从复数勒让德函数得到带符号的速度。具体而言,首先,以0次~N-1次的离散勒让德函数为基础,展开由一系列的N点构成的时间序列信号,得到A(0)~A(N-1)的展开系数,以此为基础,如下述式所示计算出一系列的复数勒让德展开系数,
C(±(2n-1))=A(2n-2)±jA(2n-1) (1≤2n-1≤N-1)
C(±(2n))=A(2n)±jA(2n-1) (2≤2n≤N-1;n:自然数,±号选取顺序一致)
…………(4)
通过该系数的相对大小进行带符号的速度分析。
在图13中,关于与图10~12相同n=8的情况,表示作为与复数正弦波输入相对的输出的复数勒让德展开系数的绝对值,作为与带符号的速度相对应的相位旋转速度的函数。复数勒让德系数的绝对值以各自的最大值进行归一化并绘图。在图13中,是交叠14个函数而进行图示,但通过区别线显示于图14~图17中。用实线表示C(1)、C(3)、C(5)、C(7),用虚线表示C(2)、C(4)、C(6),用点划线表示C(-2)、C(-4)、C(-6),用点线表示C(-1)、C(-3)、C(-5)、C(-7)。可知2N-2个复数勒让德展开系数相对于速度区分成正负而响应。已确认具有这种响应的复数勒让德系可以在N≤35的范围内进行构建。
在图18中,将由构成2N-2个复数勒让德展开系数的偶数次和奇数次勒让德系数的绝对值比决定的相位角,作为与带符号的速度相对应的相位旋转速度的函数进行显示。图中,从内侧开始,为C(±1)、C(±2)、C(±3)、C(±4)、C(±5)、C(±6)、C(±7)的偶奇相位角。各个复数勒让德展开系数的偶奇相位角作为速度的函数发生单调的变化,所以如果利用这一点,即使使用相对少数的复数勒让德展开系数,也可以进行作为连续函数的速度检测/分析。
与作为检测/分析对象的运动反射物回波相比,存在具有非常大的振幅的杂波信号,当其发生漂移时,可以按照其漂移程度的增大,使用以按顺序除掉C(±1)、C(±2)、……、和次数小的系数的复数勒让德展开系数系,进行运动反射物的速度检测/分析。
下面,对本发明的代表性构成例进行描述。
(1)本发明是一种多普勒速度检测装置,其特征在于,是对速度测量对象物多次收发脉冲波、以接收回波信号为基础进行与所述速度测量对象物有关的速度分析的多普勒速度检测装置,以下述方式构成:将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为勒让德多项式的各成分进行展开,根据各展开系数的大小得到与所述速度测量对象物有关的速度信号。
(2)是在上述(1)记载的多普勒速度检测装置中,其特征在于,以下述方式构成:以虚数单位为系数,当将所述接收回波时间序列信号作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比得到带符号的速度信号。
(3).本发明是一种多普勒速度检测装置,具备对受检体多次收发脉冲波的机构、以接收回波信号为基础进行与所述受检体内的运动反射物有关的速度分析的速度分析机构,
所述速度分析机构以下述方式构成:以虚数单位为系数,当将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比得到与所述受检体内的运动反射物有关的带符号的速度信号。
(4).是在上述(3)记载的多普勒速度检测装置中,其特征在于,设置显示机构,其对有关所述运动反射物的速度信号连同来自所述受检体内的静止反射物的回波信号一起进行二维或三维显示。
(5).本发明是一种超声波诊断装置,其特征在于,具备超声波探头、从所述超声波探头对受检体多次收发超声波脉冲的机构、以来自所述受检体的接收回波信号为基础进行所述受检体内的运动反射物的速度分析的速度分析机构,
所述速度分析机构以下述方式构成:将以发送时刻的次序排列始自所述超声波脉冲各发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为勒让德多项式的各成分进行展开,根据各展开系数的大小得到与所述受检体内的运动反射物有关的速度信号。
(6)是在上述(5)记载的多普勒速度检测装置中,其特征在于,
所述速度分析机构以下述方式构成:以虚数单位为系数,当将所述接收回波时间序列信号作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比得到带符号的速度信号。
(7)本发明是一种超声波诊断装置,其特征在于,具备超声波探头、从所述超声波探头对受检体多次收发超声波脉冲的机构、以来自所述受检体的接收回波信号为基础进行血流的速度分析的速度分析机构,所述的血流在所述受检体内的运动的脏器中流动,所述速度分析机构以下述方式构成:以虚数单位为系数,当将以发送时刻的次序排列始自所述超声波脉冲各发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比,得到与所述受检体内的血流有关的带符号的速度信号。
(8)是在上述(7)记载的多普勒速度检测装置中,其特征在于,设置显示机构,其对所述脏器的运动速度或其空间变化的分布图像和同时取得的血流分布图像进行重叠或并列显示。
(9)本发明是一种超声波诊断装置,其特征在于,从多个超声波探头对受检体多次收发超声波脉冲的超声波收发机构;对所述受检体内的超声波脉冲的发送焦点位置进行控制的发送射波束形成装置;对所述受检体内的接收焦点位置进行控制的接收射波束形成装置;对所述超声波收发机构、所述发送射波束形成装置和所述接收射波束形成装置进行控制的控制系统;和以来自所述受检体的接收回波信号为基础进行所述受检体内的运动反射物的速度分析的速度分析机构,
而且,所述速度分析机构以下述方式构成,即,以虚数单位为系数,当将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比,得到与所述受检体内的运动反射物有关的带符号的速度信号。
附图说明
图1是并列表示由N次脉冲发送所产生的的接收信号的图。
图2是将离散傅立叶展开系数的绝对值作为与速度对应的带符号频率的函数进行显示的图。
图3是将图2中的其他4个离散傅立叶展开系数的绝对值作为与速度对应的带符号频率的函数进行显示的图。
图4是将图2中的其他4个离散傅立叶展开系数的绝对值作为与速度对应的带符号频率的函数进行显示的图。
图5是表示通过使用了离散傅立叶展开的速度分析的速度频谱系数的例子。
图6是表示以往的低频截止型MTI滤波器(实线)和多项式拟合/滤波器(虚线)所具有的通过振幅特性的典型例子的图。
图7是表示分别通过以往的低频截止型MTI滤波器(实线)和多项式拟合/滤波器(虚线)之后,进行使用了离散傅立叶展开的速度分析时得到的速度频谱系数的例子的图。
图8是表示与由一系列的8点构成的信号有关的离散勒让德函数的偶数次的图。
图9是表示与由一系列的8点构成的信号有关的离散勒让德函数的奇数次的图。
图10表示将勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图11是表示将图10中的4个勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图12是表示将图10中的其他4个勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图13表示将复数勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图14是表示将图13中的±6次、±7次的复数勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图15是表示将图13中的±4次、±5次的复数勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图16是表示将图13中的±2次、±3次的复数勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图17是表示将图13中的±1次的复数勒让德展开系数的归一化绝对值作为与速度对应的带符号频率的函数的图。
图18是表示将复数勒让德展开系数的偶奇相位角作为与速度对应的带符号频率的函数的图。
图19是表示成为本发明的实施例的带血流检测/描绘功能的超声波诊断装置的构成的框图。
图20是表示将使用多个复数勒让德系数绝对值之间的比而形成的速度校正曲线作为与速度对应的带符号频率的函数的图。
图21是表示已对原点附近进行修正的复数勒让德展开系数的相位和速度之间的校正曲线的图。
图22是表示本发明的多普勒速度检测算法的例子的图。
图23是表示基于本发明的方式的速度检测结果的例子(1)的图。
图24是表示基于以往方式的速度检测结果的例子(1)的图。
图25是表示基于多项式拟合/滤波方式的速度检测结果的例子(1)的图。
图26是表示基于本发明的方式的速度检测结果的例子(2)的图。
图27是表示基于以往方式的速度检测结果的例子(2)的图。
图28是表示基于多项式拟合/滤波方式的速度检测结果的例子(2)的图。
图29是表示基于本发明的方式的速度检测结果的例子(3)的图。
图30是表示基于以往方式的速度检测结果的例子(3)的图。
图31是表示基于多项式拟合/滤波方式的速度检测结果的例子(3)的图。
具体实施方式
下面,关于本发明的实施例,参照附图进行详述。
图19是用于实施本发明而构成的带有血流描绘功能的超声波诊断装置框图的典型例子。构成超声波探头1的各元件,通过切换开关组2与发送射波束形成装置3和接收射波束形成装置10连接。在发送射波束形成装置3中,按照基于收发顺序控制部6的控制,使用通过发送波形选择部4从发送波形存储器5选择并读出的波形,通过各元件发送,此时生成具有指向性的超声波脉冲那样的信号。该信号通过超声波探头1的各元件变换成超声波脉冲并被发送给生物体,在生物体中反射或散射并返回至超声波探头1的超声波回波信号,由各元件接收并变换成电信号。在接收射波束形成装置10中,按照收发顺序控制部6的控制,为了生成具有指向性的接收灵敏度,付与各接收信号延迟时间并相互相加。通过延迟相加得到的时间序列信号,仍然按照收发顺序控制部6的控制,临时写入到利用接收存储选择部11选择的接收信号存储器12中的一个存储体(bank),在用于进行多普勒信号分析的N个时间序列信号齐备之后,读出,在13-16中进行用于检测/分析速度的信号处理。
读出的N个时间序列接收信号,在按照发送时刻的次序排列将为了获得各信号而进行的脉冲发送的时刻作为基准的时相相等的部分。首先,在相位旋转检测器13中,从含有杂波信号状态下的时间序列信号求出相位旋转。检测出相位旋转量Pn的最典型方法,基于以n=1,……,N排列的复数时间序列信号S(n)进行计算,
P(n)=S(n+1)S(n)*/||S(n+1)||/||S(n)|| (n=1,···,N-1)
………………………… (5)从P(n)的均值Pa的相位求出平均相位旋转速度。在这里,S*是S的复共轭,||S||是S的绝对值。通常,杂波/回波信号与血流回波信号相比,其信号振幅非常大,所以Pa可以认为大致是杂波/回波信号的相位旋转平均值。在下行混频器14(down mixer)中,使用在这里求得的Pa或Pa的周边区域的空间均值,进行使回波信号的相位旋转均值成为零、换言之使杂波信号的中心频率大致成为零那样的混频处理。即,通过S(n)和Pa得到下述式。
Sd(n)=S(n)Pa*n (n=1,···,N)
………………………(6)
在这里,Pa*是Pa的复共轭。
通过进行该下行混频处理,可以更有效地进行抑制杂波信号的后段处理。以得到的复数时间序列信号Sd(n)为基础,在相位旋转检测器15中,与13一样按照下述式计算下行混频处理后的相位旋转量Pd(n),求出该相位的绝对值的最大值。
Pd(n)=Sd (n+1)Sd(n)*/||Sd(n+1)||/||Sd(n)|| (n=1,··,N-1)
………………………(7)
按照该相位旋转最大值的大小,控制血流信号检测/分析部16的截止特性。即,当相位旋转最大值较大时,与较小时相比,通过进行更高设定勒让德系数的截止次数M的控制,有效抑制杂波成分。
在这里,关于相当于本发明的核心的血流信号检测/分析部16的运行,稍加详细描述。输入的时间序列信号Sd(1),……,Sd(N),首先从图9和图9的例子那样的0开始展开成N-1次的勒让德函数,求出其展开系数A(0),……,A(N-1)。该计算可以通过如下所示的阵列运算容易地进行。用行向量L(n)表示n次勒让德函数,将其从n=0排列至N-1,生成N×N阵列LL,此时,用于从以Sd(1),……,Sd(N)为元素(element)的列向量Sd得到以A(0),……,A(N-1)为元素的列向量A的阵列,可以以下述式求出。
FF=(LL·tLL)-1LL
……………………(8)
在这里,tLL表示LL的转置矩阵,LL-1表示LL的逆矩阵。如果预先准备该矩阵,可以通过下式表示的矩阵运算迅速地求出勒让德展开系数。
A=FF·Sd
…………(9)
按照式(4),从其展开系数A(0),……,A(N-1)得到复数勒让德展开系数C(±1),……,C(±(N-1)),其中,按照相位旋转最大值决定截止次数M,放弃自C(±1)至C(±(M+1))的低次数复数勒让德系数。接着,求得各复数勒让德系数的绝对值,对其进行校正,然后选择最大系数的复数勒让德展开系数。如在图13-17中例示的那样,首先,关于作为输入一定振幅的复数正弦波时的输出的复数勒让德系数绝对值,求出对与带符号的速度对应的频率进行扫描时的最大值,使其该值校正系数的绝对值。从所选择的复数勒让德系数的偶奇相位角求得与血流速度对应频率。关于偶奇相位角,如在图18中例示的那样,首先,关于作为输入一定振幅的复数正弦波时的输出的复数勒让德系数的偶奇相位角,对与带符号的速度对应的相位旋转速度进行扫描,求出速度校正曲线。
作为用于速度校正的曲线,除了使用图18那样的各复数勒让德展开系数的偶奇相位角之外,还可以使用多个复数勒让德展开系数绝对值之间的比。将其例子显示于图20。这是用于在复数勒让德展开系数的绝对值当中,当||C(n)||为最大时,使用||C(n+1)/C(n-1)||或者||C(n-1)/C(n+1)||的反正切函数作为校正曲线的例子。当||C(±7)||为最大时,与图18一样,分别使用C(±7)的偶奇相位角。另外,当||C(±1)||为最大时,使用||C(±2)/C(±1)||(±号选取顺序一致)的反正切函数。
得到的血流信号连同通过回波振幅检测器17和回波振幅压缩器18得到的来自静止脏器的回波信号被输入到扫描变换器19。在扫描变换器19中,适当重叠所输入的多个信号并用显示部进行2~3维显示,进行信号的生成/控制。如上所示,放弃C(±1)~C(±(M+1))的低次数复数勒让德展开系数,通过C(±(M+1))~C(±(N-1))的复数勒让德展开系数进行速度分析,此时,对于绝对值比最大值为勒让德系数A(M)的速度大的速度,难以进行正确的速度分析,所以需要采取一些对策。
接着,将本发明的多普勒速度检测/分析装置的运行例子与以往方式进行比较并显示。杂波回波信号的振幅为血流回波信号的振幅的300倍,关于杂波速度,假定初期值为零且以一定的加速度上升的情况。关于时间序列信号数N=8的情况,进行数值计算模拟并比较。作为以往的方式,使用图6的情况下的低频截止滤波器作为MTI滤波器,在其后段,采用与使用式(5)得到相位旋转量均值Pa的方法相同的方法,选择在血流信号检测/分析部16进行计算速度的处理的例子。进而,还与将MTI滤波器置换成如图5例示那样的多项式拟合/滤波器的情况进行比较。此时,截止次数M与本发明的方法一样,昭相位旋转检测器15的输出信号进行控制。
对于绝对值比最大值为A(M)的速度小的速度范围内的上述问题,对图20例示的速度校正曲线的原点附近进行修正,以使该范围的速度输出全部为零,使用该校正曲线进行速度检测/分析。在图21中,作为例子,显示M=3的情况下的已被修正的校正曲线。因为放弃C(±3),在||C(±4)||为最大时使用的校正曲线不是||C(±5)/C(±3)||,被改写成||C(±5)/C(±4)||的反正切函数。在这里,当从C(±(M+2))到C(±(N-1))的绝对值的最大值不到基于假定的血流回波信号的绝对值的1/100时,控制成速度输出为零。在以往方式和多项式拟合/滤波器的情况下,作为各自的MTI滤波器的输出的信号振幅,在不到假定的血流回波信号的绝对值的1/100时,控制成速度输出为零。
将以上说明的速度检测相加汇总于图22。在该图中,记录了使用绝对值最大的复数勒让德系数的相位检测速度的相加,但为了进一步进行稳定的速度检测,还有不仅使用绝对值最大系数的相位、还使用在从多个系数的相位求得的速度上加权系数绝对值的权重而得到的均值的方法。例如,当nm(>0)次的复数勒让德系数的绝对值为最大时,使用在从nm-1次至N-1次的系数的相位求得的速度上加权系数绝对值的权重而得到的均值,当-nm(<0)次的复数勒让德系数的绝对值为最大时,使用在从-nm+1次至-N+1次的系数的相位求得的速度上加权系数绝对值的权重而得到的均值。
在图23、24、25中,关于杂波的最终到达速度为乃奎斯特极限速度的0.8%的情况,分别表示基于本发明的方式、以往方式、多项式拟合/滤波方式的速度检测结果。横轴是作为输入的血流速度,纵轴是作为输出而检测的速度,用实线对有效速度检测范围进行绘图,用点线对范围外进行绘图。另外,用点划线表示理想的情况。如此,当杂波速度较低时,任何方式都能给以大致准确的速度分析结果,但基于本发明的方式的误差最小。
在图26、27、28中,同样地,关于杂波的最终到达速度为乃奎斯特极限速度的3%的情况,分别表示基于本发明的方式、以往方式、多项式拟合/滤波方式的速度检测结果。在该杂波速度域中,基于以往方式的速度检测大致都有问题。另外,多项式拟合/滤波方式尽管其有效速度检测范围内的误差不大,但在血流速度比其小的情况下的误差较大,只要不采取对策就无法使用。对此,本发明的方式不仅在有效速度检测范围内能给出准确的速度检测结果,而且在血流速度比其小的情况下,直接输出检测速度为零。
如此,在有效速度检测范围内输出准确的速度检测结果,而且,在血流速度比其小的情况下,直接显示输出检测速度为零或者不能进行速度检测,这可以说是已实施本发明的多普勒速度检测装置的特征。
在图29、30、31中,进而,关于杂波的最终到达速度为乃奎斯特极限速度的20%的情况,分别表示基于本发明的方式、以往方式、多项式拟合/滤波方式的速度检测结果。迄今为止,当杂波速度增大时,基于以往方式的速度检测都有问题。另外,多项式拟合/滤波方式,其有效速度检测范围内的误差增大,同时在血流速度比其小的情况下的误差也显著增大。对此,本发明的方式尽管其有效速度检测范围变得稍窄,当血流速度为乃奎斯特极限速度的60%以上的范围时,不仅能给出大致准确的速度检测结果,而且在血流速度比其小的情况下的运行也仍然比较自然。
当对描绘如下所述的低速血流的情况进行具体描述时,通过本发明,在脏器向超声波探头的方向的运动速度以1mm/s改变的情况下,可以实时地描绘向超声波探头的方向的速度成分3mm/s以上的血流,其中,所述的低速血流在通过呼吸等有软组织运动的肝脏和肾脏中流过。这在以往的方式中是根本无法实现的血流速度检测性能。
如此通过本发明,可以将血流和脏器的运动严加区别并描绘。为了进一步发挥该优势,在图20的诊断装置中,将通过相位旋转检测器13检测出来的杂波信号即脏器回波的相位旋转速度的信号输入到扫描变换器19中,可以将表示脏器的运动速度和其空间微分的分布的图像与血流图像重叠或者并列显示。以肝癌的情况为例说明该构成的有用性。在肝癌的边缘部,新生血管发达,血流动态与周边的正常肝不同。另外,因为与周围的正常肝的硬度不同,所以组织的运动也与周边不同。因此,除了血流图像之外,表示由脏器的运动速度的场所导致的变化,可以提供对肝癌的诊断极其有用的图像。
如上所述,通过本发明,就生物体内部的血流等运动的反射物的速度而言,除掉与其回波信号相比振幅非常大的杂波回波信号的影响,可以对其进行准确检测。具体而言,在脏器向超声波探头的方向的运动速度以1mm/s改变的情况下,可以实时地描绘向超声波探头的方向的速度成分3mm/s以上的血流。
如此,根据本发明,通过可以区别正负速度的复数勒让德展开系,提供抑制由振幅比血流大2~3数量级的体动引起的杂波而可以检测血流速度的新方式,另外,可以提供能进行准确医疗诊断的带血流检测/描绘功能的超声波诊断装置。即,已实施本发明的装置在医用诊断上的有用性极大,因此,本发明在支持医用诊断的工业中的意义也非常大。
进而,本发明的方法可以飞跃性地提高通过收发电磁波检测雨云等运动的反射物并进行描绘的气象雷达、对飞行物体进行检测的航空雷达、或者、对接近的物体进行检测的碰撞防止雷达等脉冲/多普勒/雷达装置等其他脉冲/多普勒装置的运动反射体检测能力,本发明的工业和社会意义在这一层面也非常大。
工业上的可利用性
通过本发明,就生物体内部的血流等运动的反射物的速度而言,可以抑制与其回波信号相比振幅非常大的杂波回波信号的影响,而准确检测出反射物的速度。
Claims (10)
1.一种多普勒速度检测装置,是对速度测量对象物多次收发脉冲波、以接收回波信号为基础进行与所述速度测量对象物有关的速度分析的多普勒速度检测装置,
以下述方式构成:将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为勒让德多项式的各成分进行展开,根据各展开系数的大小得到与所述速度测量对象物有关的速度信号。
2.如权利要求1所述的多普勒速度检测装置,其特征在于,
以下述方式构成:以虚数单位为系数,当将所述接收回波时间序列信号作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比得到带符号的速度信号。
3.一种多普勒速度检测装置,具备对受检体多次收发脉冲波的机构、以接收回波信号为基础进行与所述受检体内的运动反射物有关的速度分析的速度分析机构,
所述速度分析机构以下述方式构成:以虚数单位为系数,当将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比得到与所述受检体内的运动反射物有关的带符号的速度信号。
4.如权利要求3所述的多普勒速度检测装置,其特征在于,
设置显示机构,其对有关所述运动反射物的速度信号连同来自所述受检体内的静止反射物的回波信号一起进行二维或三维显示。
5.一种超声波诊断装置,具备超声波探头、从所述超声波探头对受检体多次收发超声波脉冲的机构、以来自所述受检体的接收回波信号为基础进行所述受检体内的运动反射物的速度分析的速度分析机构,
所述速度分析机构以下述方式构成:将以发送时刻的次序排列始自所述超声波脉冲各发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为勒让德多项式的各成分进行展开,根据各展开系数的大小得到与所述受检体内的运动反射物有关的速度信号。
6.如权利要求5所述的超声波诊断装置,其特征在于,
所述速度分析机构以下述方式构成:以虚数单位为系数,当将所述接收回波时间序列信号作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比得到带符号的速度信号。
7.一种超声波诊断装置,具备超声波探头、从所述超声波探头对受检体多次收发超声波脉冲的机构、以来自所述受检体的接收回波信号为基础进行血流的速度分析的速度分析机构,所述的血流在所述受检体内的运动的脏器中流动,
所述速度分析机构以下述方式构成:以虚数单位为系数,当将以发送时刻的次序排列始自所述超声波脉冲各发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比,得到与所述受检体内的血流有关的带符号的速度信号。
8.如权利要求7所述的超声波诊断装置,其特征在于,
设置显示机构,其对所述脏器的运动速度或其空间变化的分布图像和同时取得的血流分布图像进行重叠或并列显示。
9.如权利要求7或者8所述的超声波诊断装置,其特征在于,
以下述方式构成:在所述脏器向所述超声波探头的方向的运动速度以1mm/秒变化当中,检测并显示向所述超声波探头的方向的速度成分3mm/秒以上的血流。
10.一种超声波诊断装置,具备:从多个超声波探头对受检体多次收发超声波脉冲的超声波收发机构;对所述受检体内的超声波脉冲的发送焦点位置进行控制的发送射波束形成装置;对所述受检体内的接收焦点位置进行控制的接收射波束形成装置;对所述超声波收发机构、所述发送射波束形成装置和所述接收射波束形成装置进行控制的控制系统;和以来自所述受检体的接收回波信号为基础进行所述受检体内的运动反射物的速度分析的速度分析机构,
而且,所述速度分析机构以下述方式构成,即,以虚数单位为系数,当将以发送时刻的次序排列始自各脉冲发送时刻的经过时间相等的接收回波信号而成的接收回波时间序列信号,作为从0次开始的勒让德多项式的各成分进行展开时,对此时的偶数次项的展开系数和次数与其差1个的奇数次项的展开系数进行线性结合,得到复数展开系数,根据各复数展开系数的大小和各复数展开系数间的大小的比,得到与所述受检体内的运动反射物有关的带符号的速度信号。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003397074 | 2003-11-27 | ||
JP397074/2003 | 2003-11-27 | ||
PCT/JP2004/010157 WO2005051203A1 (ja) | 2003-11-27 | 2004-07-09 | ドプラ速度検出装置及びそれを用いた超音波診断装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1886094A true CN1886094A (zh) | 2006-12-27 |
CN1886094B CN1886094B (zh) | 2010-12-01 |
Family
ID=34631530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2004800349981A Expired - Fee Related CN1886094B (zh) | 2003-11-27 | 2004-07-09 | 多普勒速度检测装置以及使用该装置的超声波诊断装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7942821B2 (zh) |
EP (1) | EP1688093B1 (zh) |
JP (1) | JP4369427B2 (zh) |
CN (1) | CN1886094B (zh) |
WO (1) | WO2005051203A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176859A (zh) * | 2008-10-13 | 2011-09-07 | 埃西斯创新有限公司 | 对象的物理特性的研究 |
CN106232016A (zh) * | 2014-04-21 | 2016-12-14 | 株式会社日立制作所 | 超声波诊断装置 |
JP2017153592A (ja) * | 2016-02-29 | 2017-09-07 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、信号処理装置および信号処理プログラム |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4634738B2 (ja) * | 2004-04-30 | 2011-02-16 | 株式会社東芝 | 超音波診断装置 |
US7590184B2 (en) * | 2005-10-11 | 2009-09-15 | Freescale Semiconductor, Inc. | Blind preamble detection for an orthogonal frequency division multiplexed sample stream |
US7740611B2 (en) | 2005-10-27 | 2010-06-22 | Patented Medical Solutions, Llc | Method and apparatus to indicate prior use of a medical item |
US7623599B2 (en) * | 2005-11-21 | 2009-11-24 | Freescale Semiconductor, Inc. | Blind bandwidth detection for a sample stream |
US7675844B2 (en) * | 2006-02-24 | 2010-03-09 | Freescale Semiconductor, Inc. | Synchronization for OFDM signals |
US20080025197A1 (en) * | 2006-07-28 | 2008-01-31 | Mccoy James W | Estimating frequency error of a sample stream |
JP5593116B2 (ja) * | 2009-04-30 | 2014-09-17 | 三星メディソン株式会社 | 適応クラッターフィルタリングを行う超音波システムおよび方法 |
JP5498346B2 (ja) * | 2010-10-07 | 2014-05-21 | 株式会社東芝 | 超音波診断装置 |
EP2654569B1 (en) * | 2010-12-23 | 2020-04-15 | Koninklijke Philips N.V. | Wall filter for ultrasonic mitral regurgitation analysis |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2563656B2 (ja) | 1990-07-30 | 1996-12-11 | 松下電器産業株式会社 | 超音波ドプラ映像装置 |
US5228009A (en) * | 1992-04-10 | 1993-07-13 | Diasonics, Inc. | Parametric clutter elimination |
US5622174A (en) * | 1992-10-02 | 1997-04-22 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and image displaying system |
US5315562A (en) * | 1992-10-23 | 1994-05-24 | Rowe, Deines Instruments Inc. | Correlation sonar system |
US5570691A (en) * | 1994-08-05 | 1996-11-05 | Acuson Corporation | Method and apparatus for real-time, concurrent adaptive focusing in an ultrasound beamformer imaging system |
-
2004
- 2004-07-09 JP JP2005515729A patent/JP4369427B2/ja not_active Expired - Fee Related
- 2004-07-09 US US10/580,225 patent/US7942821B2/en not_active Expired - Fee Related
- 2004-07-09 CN CN2004800349981A patent/CN1886094B/zh not_active Expired - Fee Related
- 2004-07-09 EP EP04747623A patent/EP1688093B1/en not_active Expired - Fee Related
- 2004-07-09 WO PCT/JP2004/010157 patent/WO2005051203A1/ja not_active Application Discontinuation
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176859A (zh) * | 2008-10-13 | 2011-09-07 | 埃西斯创新有限公司 | 对象的物理特性的研究 |
CN102176859B (zh) * | 2008-10-13 | 2014-01-15 | 埃西斯创新有限公司 | 对象的物理特性的研究 |
CN106232016A (zh) * | 2014-04-21 | 2016-12-14 | 株式会社日立制作所 | 超声波诊断装置 |
CN106232016B (zh) * | 2014-04-21 | 2019-09-24 | 株式会社日立制作所 | 超声波诊断装置 |
JP2017153592A (ja) * | 2016-02-29 | 2017-09-07 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、信号処理装置および信号処理プログラム |
Also Published As
Publication number | Publication date |
---|---|
US20070167791A1 (en) | 2007-07-19 |
EP1688093B1 (en) | 2013-01-02 |
EP1688093A4 (en) | 2011-04-13 |
JPWO2005051203A1 (ja) | 2007-06-14 |
WO2005051203A1 (ja) | 2005-06-09 |
US7942821B2 (en) | 2011-05-17 |
JP4369427B2 (ja) | 2009-11-18 |
CN1886094B (zh) | 2010-12-01 |
EP1688093A1 (en) | 2006-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Perrot et al. | So you think you can DAS? A viewpoint on delay-and-sum beamforming | |
CN103181780B (zh) | 信号分析方法、超声波图像分析方法以及超声波成像系统 | |
US6309356B1 (en) | Method and apparatus for forming medical ultrasound images | |
KR101586998B1 (ko) | 의료용 초음파 이미징에서의 전단파 정보의 해석 | |
CN105816205B (zh) | 声辐射力脉冲成像中的稀疏跟踪 | |
CN103462643B (zh) | 一种剪切波速度测量方法、装置和系统 | |
Shin et al. | Estimation of average speed of sound using deconvolution of medical ultrasound data | |
US20130258805A1 (en) | Methods and systems for producing compounded ultrasound images | |
WO2010104863A9 (en) | Method for ultrasound vibrometry using orthogonal basis functions | |
CN1886094A (zh) | 多普勒速度检测装置以及使用该装置的超声波诊断装置 | |
US9883851B2 (en) | System and method for shear wave generation with steered ultrasound push beams | |
JP2014515980A (ja) | 超音波アレイ用の分離可能なビーム形成 | |
CN107714091A (zh) | 经颅低频超声线性调频脉冲逆转微泡成像方法 | |
EP2610641A2 (en) | Ultrasound and system for forming a Doppler ultrasound image | |
US6517489B1 (en) | Method and apparatus for forming medical ultrasound images | |
CN105919624B (zh) | 一种高脉冲重复扫查频率信号的编码和装置 | |
Cacko et al. | Shear wave elastography implementation on a portable research ultrasound system: Initial results | |
KR101193710B1 (ko) | 적응적 클러터 필터링을 수행하는 초음파 시스템 및 방법 | |
Peng et al. | Enhanced ultrasound strain imaging using chirp-coded pulse excitation | |
Gray et al. | Compensation of array lens effects for improved co-registration of passive acoustic mapping and B-mode images for cavitation monitoring | |
Bottenus | Implementation of constrained swept synthetic aperture using a mechanical fixture | |
CN111789632A (zh) | 超声波诊断装置、信号处理装置以及存储介质 | |
JP5355327B2 (ja) | 超音波診断装置 | |
CN108784744A (zh) | 用于剪切波成像的可变焦点 | |
Liang et al. | Continuous wave ultrasonic Doppler tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101201 Termination date: 20160709 |
|
CF01 | Termination of patent right due to non-payment of annual fee |