CN101646390A - 用于超声波成像技术领域的方法和仪器 - Google Patents

用于超声波成像技术领域的方法和仪器 Download PDF

Info

Publication number
CN101646390A
CN101646390A CN200880010582A CN200880010582A CN101646390A CN 101646390 A CN101646390 A CN 101646390A CN 200880010582 A CN200880010582 A CN 200880010582A CN 200880010582 A CN200880010582 A CN 200880010582A CN 101646390 A CN101646390 A CN 101646390A
Authority
CN
China
Prior art keywords
frequency
prf
max
peak
minimum
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
Application number
CN200880010582A
Other languages
English (en)
Other versions
CN101646390B (zh
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.)
Hitachi Ltd
Original Assignee
Arockar K K
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 Arockar K K filed Critical Arockar K K
Publication of CN101646390A publication Critical patent/CN101646390A/zh
Application granted granted Critical
Publication of CN101646390B publication Critical patent/CN101646390B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • 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/50Systems of measurement, based on relative movement of the target
    • G01S15/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S15/582Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets
    • G01S15/584Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets with measures taken for suppressing velocity ambiguities, i.e. anti-aliasing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

获得多普勒频谱中的最大频率并用作混淆检测器。当混淆发生时,大于频率界限的频率从一个频率区域变化到另一个频率区域。当检测到混淆时,移动零频率基线以防止以后的混淆。

Description

用于超声波成像技术领域的方法和仪器
技术领域
本发明主要涉及超声波成像领域。特别的,本发明的具体实施例涉及用于频谱图像的方法和系统。
背景技术
超声波被用于对各种器官、心脏、肝脏、胎儿和血管进行成像。对于心血管疾病的诊断,频谱多普勒通常被用于测量血流速度。由于与没有空间辨别能力的连续波(CW)技术相比,脉冲频谱多普勒技术具有对血管中的速度进行采样的空间采样能力,而且可以沿着超声波束对所有的信号进行采样,因此脉冲频谱多普勒技术通常得到使用。
多普勒技术中,超声波以脉冲重复频率(PRF)进行传输,并且随着在接收到的超声波信号中的频移(多普勒移频)而检测血流速度。接收到的超声波与和传输超声波频率相同频率的同相(0度)和正交(90度)的基准信号混合。经过低通滤波高频部件(即二次谐波)之后,只得到了基带信号。壁滤波(wall-filtering)(如高通滤波)应用于基带信号以从组织去除强杂波噪声和缓慢移动如血管壁的组织,从而产生复杂I-Q多普勒信号。
通常,I-Q多普勒信号被输入到诸如快速傅立叶变换(FFT)的频谱分析器,以得到代表血流速度的多普勒频谱。多普勒移频与血流速度存在如下关系:
Δf = 2 f t v cos θ c - - - ( 1 )
此处Δf是多普勒移频,ft是传输频率,v是血流速度,θ是超声波束方向和速度矢量间的夹角,c是声速。
通常使用128点、256点或512点快速傅立叶变换(FFT)。因为多普勒信号是通过脉冲超声波(以及采样)技术得到的,所以采样理论规定了最大频率界限。最大频率通常是脉冲重复频率(PRF)或fPRF的一半。因为在复杂I-Q多普勒信号上执行FFT,所以负方向的血流速度出现于负频域。因此,多普勒频谱FFT输出具有相应于负速度的负频率。这样,多普勒频谱通常具有从
Figure G2008800105824D00022
的频率范围。但是,负频率范围可分配为表示大于
Figure G2008800105824D00023
且一直到fPRF的正频率。在相反的情形中,正频率范围可分配为表示小于且一直到-fPRF的负频率。多普勒频谱模型中,上述情形是通过基线移动执行的。基线移动在正频率或负频率方向上移动零频率基线的位置。这样,在极端的情形中由于基线移动的作用,多普勒频谱可具有从-fPRF到0或从0到fPRF的范围。全部频率范围恒为fPRF
通常在心血管应用方面,血流速度可能超过这些最大速度,从而导致混淆(aliasing)。混淆发生时,频谱可环绕在正的最大频率处,超过最大界限的频率出现于负频率;或者可环绕在负的最大频率处,超过负的最大界限的频率出现于正频率。混淆使得确定血流速度变得困难。
相反地,fPRF可能太大以至于不能准确测量血流速度。最大血流速度(最大频率)可能只有最大频率界限的大约十分之一,造成所显示的频谱太小以至于不能准确测量。
大多数超声波应用中,用户手动调节相应于血流速度的PRF,和/或手动调节相应于频谱标度中的零速度而作为零频率位置的基线。但是,用户在调节这些设定中所花费的时间更应该用于诊断上面。
这就存在解决这些问题的必要。
发明内容
本发明的发明人发现人们期望存在一种获得多普勒频谱中的最大频率并用作混淆检测器的系统和方法。当混淆发生时,最大频率从正频率卷绕到负频率,或从负频率卷绕到正频率。当检测到混淆时,基线就被移动以适应在正确的频率极性中被卷绕频率的大小。
本发明的一个方案提供了一种用于多普勒频谱中探测和校正混淆的方法。根据本发明的该方案的方法包括随时间接收多普勒频谱信号;从多普勒频谱计算最大频率fmax和最小频率fmin;随时间追踪最大频率fmax和最小频率fmin;如果在正频率区域中的频率变化(卷绕)到负频率区域,那么检测混淆是否发生自最大频率fmax;或如果在负频率区域中的负频率变化(卷绕)到正频率区域,那么检测混淆是否发生自最小频率fmin;并且如果检测到混淆,那么根据最大频率偏差fa在正或负方向上移动分隔多普勒频谱的负和正频率区域的零频率基线。
本发明的另一方案提供了一种用于确定超声波系统的脉冲重复频率的方法。根据本发明的该方案的方法包括随时间接收多普勒频谱信号;从多普勒频谱计算最大频率fmax,从多普勒频谱计算最小频率fmin;随时间追踪最大频率fmax和最小频率fmin;捕捉被追踪的最大频率fmax的最高值high fmax和最小频率fmin的最低值low fmin;将最高值high fmax和最低值low fmin作比较以确定最大频率fmax和最小频率fmin是否是双极的、还是负或正的单极的;如果最大频率fmax和最小频率fmin是双极的,则:基于最高最大频率high fmax和最低最小频率low fmin之间的差值确定频带宽度,将频带宽度与当前PRF设定值作比较,如果频带宽度大于当前PRF设定值,则增大PRF设定值,如果频带宽度小于当前PRF设定值的预定分数,则减小PRF设定值,如果频带宽度小于当前PRF设定值但同时大于当前PRF的预定分数,则使用当前PRF设定值;如果最大频率fmax和最小频率fmin是正的单极的,则:将当前正的最大频率界限b1fPRF与最高最大频率high fmax作比较,其中如果最高最大频率highfmax大于当前正的最大频率界限b1fPRF,则增大当前PRF设定值到相应于最高最大频率high fmax的设定,如果最高最大频率high fmax小于当前正的最大频率界限b1fPRF,则将低级阈值b2b1fPRF与最高最大频率high fmax作比较,其中如果最高最大频率high fmax小于低级阈值b2b1fPRF,那么减小PRF直到与最高最大频率high fmax相等;如果最大频率fmax和最小频率fmin是负的单极的,则:将当前负的最大频率界限-(1-b1)fPRF的绝对值与最低最小频率low fmin的绝对值作比较,其中如果最低最小频率low fmin的绝对值大于当前负的最大频率界限-(1-b1)fPRF的绝对值,那么增大当前PRF设定值到相应于最低最小频率lowfmin的绝对值的设定,如果最低最小频率low fmin的绝对值小于当前负的最大频率界限-(1-b1)fPRF的绝对值,那么将低级阈值-b2(1-b1)fPRF的绝对值与最低最小频率low fmin的绝对值作比较,其中如果最低最小频率low fmin的绝对值小于低级阈值-b2(1-b1)fPRF的绝对值,那么减小PRF到与最低最小频率low fmin的绝对值相等。
本发明的另一方案提供了一种用于确定超声波系统的脉冲重复频率的方法。根据本发明的该方案的方法包括设定初始脉冲重复频率;随时间接收多普勒频谱信号;从多普勒频谱计算最大频率fmax,从多普勒频谱计算最小频率fmin;随时间追踪最大频率fmax和最小频率fmin;捕捉被追踪的最大频率fmax的最高值high fmax和最小频率fmin的最低值low fmin;将最低最小频率low fmin的绝对值与最高最大值high fmax的绝对值作比较,以确定是正频区域还是负频率区域占优势;如果最高最大值high fmax较大,那么正频率区域占优势并且计算正的低级阈值b2b1fPRF,并且将正的最大频率界限b1fPRF和正的低级阈值b2b1fPRF与最高最大频率high fmax作比较,其中如果最高最大频率high fmax小于正的低级阈值b2b1fPRF,那么减小PRF直到正的最大频率界限b1fPRF与最高最大频率high fmax相等,或混淆开始发生在负的最大频率界限-(1-b1)fPRF处,无论哪个先发生,而且其中如果最高最大频率high fmax大于正的最大频率界限b1fPRF,那么增大PRF到与最高最大频率high fmax相等,如果最低最小频率low fmin的绝对值较大,那么负频率区域占优势而且计算低级阈值-b2(1-b1)fPRF,将负的最大频率界限-(1-b1)fPRF的绝对值和低级阈值-b2(1-b1)fPRF的绝对值与最低最小频率low fmin的绝对值作比较,其中如果最低最小频率low fmin的绝对值小于低级阈值-b2(1-b1)fPRF的绝对值,那么减小PRF直到负的最大频率界限-(1-b1)fPRF的绝对值为最低最小频率low fmin的绝对值或混淆开始发生于正频率界限处,无论哪个先发生,而且其中如果最低最小频率lowfmin的绝对值大于负的最大频率界限-(1-b1)fPRF的绝对值,那么增大PRF到与最低最小频率low fmin相等。
本发明另一方案提供了一种用于多普勒频谱中检测和校正混淆的系统。根据本发明的该方面的系统包括用于随时间接收多普勒频谱信号的装置;用于从多普勒频谱计算最大频率fmax和最小频率fmin的装置;用于随时间追踪最大频率fmax和最小频率fmin的装置;如果在正频率区域中的频率变化(卷绕)到负频率区域,用于检测混淆是否发生自最大频率fmax的装置;如果在负频率区域中的负频率fmin变化(卷绕)到正频率区域,用于检测混淆是否发生自最小频率fmin的装置;以及如果检测到混淆,根据最大频率偏差fa在正或负方向移动分隔多普勒频谱的负和正频率区域的零频率基线的装置。
本发明另一方案提供了一种用于确定超声波系统的脉冲重复频率的系统。根据本发明的该方面的系统包括用于设定初始脉冲重复频率的装置;用于随时间接收多普勒频谱信号的装置;用于从多普勒频谱计算最大频率fmax的装置,用于从多普勒频谱计算最小频率fmin的装置;用于随时间追踪最大频率fmax和最小频率fmin的装置;用于捕捉被追踪的最大频率fmax的最高值highfmax和最小频率fmin的最低值low fmin的装置;用于将最低最小频率low fmin的绝对值与最高最大值high fmax的绝对值作比较以确定正或负频率区域哪个占优势的装置,如果最高最大值high fmax较大,那么正频率区域占优势而且计算正的低级阈值b2b1fPRF;和用于将正的最大频率界限b1fPRF和正的低级阈值b2b1fPRF与最高最大频率high fmax作比较的装置,其中如果最高最大频率highfmax小于正的低级阈值b2b1fPRF,则不管哪个先发生减小PRF直到正的最大频率界限b1fPRF与最高最大频率high fmax相等,或者混淆开始发生于负的最大频率界限-(1-b1)fPRF处,而且其中如果最高最大频率high fmax大于正的最大频率界限b1fPRF,那么增大PRF到等于最高最大频率high fmax,如果最低最小频率low fmin的绝对值较大,那么负频率区域占优势而且计算低级阈值-b2(1-b1)fPRF;用于将负的最大频率界限-(1-b1)fPRF的绝对值和低级阈值-b2(1-b1)fPRF的绝对值与最低最小频率low fmin的绝对值作比较的装置,其中如果最低最小频率low fmin的绝对值小于低级阈值-b2(1-b1)fPRF的绝对值,那么减小PRF直到负的最大频率界限-(1-b1)fPRF的绝对值为最低最小频率low fmin的绝对值,或者混淆就开始发生于正频率界限处,无论哪一个先发生,而且其中如果最低最小频率low fmin的绝对值大于负的最大频率界限-(1-b1)fPRF的绝对值,那么增大PRF直到与最低最小频率low fmin相等。
在附图和下面的描述中将对本发明的一个或多个具体实施例的细节予以阐述。本发明的其它特点、目的和优点将从说明书、附图和权利要求的阐述中变得清晰。
附图说明
图1是表示展示出混淆的最大多普勒频率的示意图。
图2A是表示百分位数频谱的最大多普勒频率的示意图。
图2B是表示百分位数频谱的最小多普勒频率的示意图。
图3是表示校正的基线移动后的最大多普勒频率的示意图。
图4是表示频谱的最小、平均和最大频率的示意图。
图5A是表示多普勒频谱的双极最大和最小频率的示意图。
图5B是表示多普勒频谱的单极正的最大和最小频率的示意图。
图5C是表示多普勒频谱的单极负的最大和最小频率的示意图。
图6是描述自动基线移动方法的示范性流程图。
图7是描述自动PRF设定和基线移动方法的示范性流程图。
图8是描述带有固定基线的自动PRF设定方法的示范性流程图。
图9是带有自动基线移动和PRF设定的示范性超声波系统。
图10是表示随时间变化的最大和最小频率的示范性多普勒频谱。
具体实施方式
下面将参照附图对本发明的具体实施例进行描述,其中相同的数字始终代表相同的元件。在详细解释本发明的具体实施例之前,需要理解的是:本发明并不仅限应用于下面描述中所阐述或附图中所阐明的例子的细节。本发明还能够以多种应用和多种方式实施或执行其它实施例。而且,需要理解的是:此处所使用的措词和术语是为了描述的目的,而不应被认为是进行限制。这里使用的“包含”、“包括”或“具有”及其变形意味着包括后面所列的项目及其等价物以及附加项目。术语“安装的”、“连接的”和“结合的”被广泛地使用而且包括直接和间接的安装、连接和结合。甚至,“连接的”和“结合的”并不限于物理或机械的连接或结合。
应该注意的是,本发明并不限于所描述的或在图中所暗示的任何特殊的软件语言。本领域中的任一普通技术人员将会理解可使用多种可选择的软件语言来实施本发明。同样可以理解的是,其中的一些部件和项目被当作硬件阐明和描述,因为这在本领域内是常用做法。但是,基于对上述详细描述的阅读,本领域中的任一普通技术人员将理解在至少一个实施例中,本方法和系统中的部件可以以软件或硬件实施。
图9表示带有自动基线移动和PRF设定的超声波系统901。图6、图7和图8表示了描述系统901所使用的各种方法的流程图。超声波信号从通过发送/接收开关907的发送器905驱动的超声波探头903发送。接收器909接收通过开关907从探头903接收到的超声波信号,并且处理信号911。
被处理的信号913结合到多普勒频谱处理器915、色流处理器(color flowprocessor)921和B型图像处理器923。多普勒频谱处理器915包括多普勒信号处理器917和频谱分析器919,并且处理多普勒流速信号,并计算且输出多普勒频谱925。色流处理器921处理接收到的信号913,并计算且输出速度、功率和变化信号927。B型图像处理器923处理接收到的信号913,并计算且输出B型图像929或者通过振幅检测的信号振幅。
多普勒频谱信号925、色流处理器信号(速度、功率和变化)927和B型处理器信号929结合到可将信号转换成扫描转换信号的扫描转换器931。扫描转换器931的输出结合到用于显示超声波图像的显示监视器933。
被处理的信号913结合到用于计算时域中的多普勒流信号的多普勒信号处理器917。多普勒流速信号结合到频谱分析器919,该频谱分析器919将时域多普勒信号转换成其频谱频率分量925。频率分量或频谱925间接地结合到脉冲重复频率(PRF)发生器935。依靠来自通过开关939结合到PRF生成器935的手动用户输入937或来自自动基线移动和PRF设定处理器941的输入,PRF发生器935产生脉冲重复频率(PRF)。自动基线移动和PRF设定处理器941包括PRF设定装置943、基线位置装置945和可实施为DSP、FPGA、ASIC或作为分立部件的处理器947。处理器947得到基线移动和/或结合到PRF发生器935的PRF设定。基线移动或者受控于通过开关959的用户输入961,或者自动地受控于通过开关959的基线位置装置945。开关959允许用户在用户输入模式或自动模式之间进行选择。
处理器947包括引擎,所述引擎可以计算最大频率和最小频率949,检测混淆和偏差951,以及从多普勒频谱925中追踪最大频率953、最小频率955和平均频率957。处理器947通过分析多普勒频谱925来最优化成像,并且如果必要就产生PRF设定943和基线零频率移动945。
参照图6,在使用中,超声波系统901可使用用于如心脏、颈动脉或肝脏成像的特殊应用的缺省PRF,以观察血流多普勒频谱(步骤602)。最大PRF是超声波系统的最高频率范围或最高速度范围。
如图10所示,多普勒频谱图像输出925典型地是随时间变化的频谱,或是以功率作为亮度相对于时间(横轴)变化的频率(竖轴)。多普勒频谱的亮度表示在频率处的频谱功率。最大多普勒频率949计算自多普勒频谱925,而且最大多普勒频率949作为如图10所示的最大频率的曲线随时间被追踪。
最大频率引擎949将最大频率计算为百分位数频率。通过对所有频率中的功率进行积分而首先得到多普勒频谱的总面积,如下面表达式的分母所示,
∫ 0 f max pdf ∫ pdf = 0.999 - - - ( 2 )
此处p是频谱功率(或频谱振幅谱a、或提高到ab的功率,此处b是实数或任何获得自振幅的信号)。如99%或99.9%的百分位数应用于产生百分位数面积的全部面积(即公式(2)的分母)。第二个积分(公式(2)的分子)始于零频率,且当此积分达到百分位数面积时终止。最大频率是积分停止处的频率。在频谱发生混淆的情况下,即使积分(公式(2)的分子)达到最大频率范围,公式(2)也可能得不到满足。在这种情况下,积分继续到负的最大频率范围并且在负频率范围中继续进行到零频率,直到公式(2)得到满足。
图2A表示在给定时间作为频率相对功率的多普勒频谱。图2A所示的多普勒频谱表示了,对于介于正、负频率范围界限之间、即-(1-b1)fPRF到b1fPRF之间的频谱样本,99%的频率代表了最大频率值fmax(步骤604),此处b1为介于0和1之间的分数而且确定了零频率基线的位置,这样就确定了从0到b1fPRF的正频率范围和从-(1-b1)fPRF到0的负频率范围。如果b1=1/2,那么正频率范围和负频率范围相等。对于每个频谱样本的最大频率值fmax就如同曲线的最大频率值那样随时间而被追踪。
可使用降噪技术来降低来自多普勒频谱925中的噪声。多普勒频谱功率可由降噪增益控制来抑制。功率谱可由振幅谱a或提高到ab的功率来代替,此处b是实数或任何来自振幅的信号。
图1表示了被混淆的最大频率fmax曲线101。最大频率曲线101可相对于零频率基线103在正或负频率方向上移动。
但是,如果最大频率fmax超过PRF频率范围界限,即正的最大频率界限b1fPRF或负的最大频率界限-(1-b1)fPRF,那么如b1fPRF处所示,大于频率界限的频率就改变(卷绕)到相反的最大频率区域。所述突然的极性改变由混淆检测器和偏差引擎951检测为混淆(步骤606,610)。在频率由正转变到负105的基线附近,可以自然地发生极性变化而没有混淆。
当检测到混淆时,由偏差引擎951计算最大频率偏差fa,该最大频率偏差fa相应于从最大正频率范围界限b1fPRF或最大负频率范围界限-(1-b1)fPRF卷绕的频率的振幅。图1中,计算来自负的最大频率范围-(1-b1)fPRF的最大偏差fa。当PRF太小而且混淆发生时,一个以上的频率极值可能发生混淆(频率卷绕),如fa1、fa2、fa3等。混淆检测器和偏差引擎951检测每个混淆(频率卷绕),并且将所有被混淆的频率作比较以发现在观察期间内的最大频率偏差fa
取决于正或负频率是否被混淆,最大频率偏差fa被用于在正或负频率方向上偏移基线103。预定频率安全余量fs可被增加到最大频率偏差fa,以确保在实施基线103的移动之后,将不会有大于最大正频率界限b1fPRF或最大负频率界限-(1-b1)fPRF的频率。基线移动由下式决定:
基线移动=±(fa+fs)               (3)
式(3)中的符号显示了基线移动的方向。负号显示基线在负频率方向移动,而正号显示基线在正频率方向的移动。
图3表示基线移动到图1中被混淆的最大频率fmax曲线101的结果。基线移动301在正或负频率方向调节基线以得到没有被混淆的最大频率fmax曲线303。因为图1中的最大频率偏差fa在负频率范围被检测到,所以式(3)中的方向是负的,并且基线103由被计算出的基线移动303所代替,该基线移动303包括最大频率偏差fa和预定的频率安全余量fs(式(3))(步骤608)。图6中的方法调节了基线而且维持了恒定的PRF设定。
当基线被移动时,正和负频率范围随着基线移动而改变。在基线移动之后,正的最大频率界限变成b1fprf+fa+fs1,而负的最大频率界限变成-(1-b1)fprf+fa+fs1。例如,如果用于图1的被计算出的基线偏移为
Figure G2008800105824D00091
(3),那么图3中的基线301在负频率方向被移动
Figure G2008800105824D00092
如果当前PRF分数b1为1/2,那么这意味着负和正频率范围是从
Figure G2008800105824D00093
到0和从0到
Figure G2008800105824D00094
新的负频率范围变成从
Figure G2008800105824D00095
到0而且新的正频率范围变成从0到
Figure G2008800105824D00096
基线移动调节PRF分数b1
b1 new fPRF=b1 current fprf-基线移动                (4)
以上描述了使用最大频率fmax的基线移动来校正在正频率界限处所经历的混淆。这种方法通过使用最小频率fmin而应用于负频率界限处。最小频率fmin被计算为百分位数值。如下的表达式的分母所示,首先通过对所有频率的功率积分得到多普勒频谱的总面积,
∫ f min 0 pdf ∫ pdf = 0.999 - - - ( 5 )
此处p是频谱功率(或频谱振幅谱a、或提高到ab的功率,此处b为实数或从振幅导出的任何信号)。如99%或99.9%的百分位数被应用于产生百分位数面积的总面积。第二个积分(式(5)的分子)始于零频率且当积分达到如图2B所示的百分位数面积时终止。最大频率是积分停止处的频率。在混淆涉及负的最大频率的情况下,将使用最大频率的基线移动简单地通过最小频率转换成基线移动。当最小频率从负的最大频率界限改变(卷绕)到正的最大频率界限时,检测到位于负频率范围处的混淆。被混淆的部分将由前述混淆在正的最大频率范围处相反方向上的基线移动所校正。
此外,最大频率和最小频率可由如下的交替的方法而得到。
首先,平均频率fmean由下式得到:
f mean = ∫ fpdf ∫ pdf - - - ( 6 )
然后,最大频率fmax和最小频率fmin由下式计算:
∫ f mean f max pdf ∫ pdf = 0.499 , 及                                                        (7)
∫ f min f mean pdf ∫ pdf = 0.499 - - - ( 8 )
此处f是频率,p是多普勒频谱功率(或频谱振幅谱a,或提高到ab的功率,此处b是实数或从振幅导出的任何信号)。
图7表示描述了同样包括调节PRF设定的基线移动的变量的流程图。最大PRF首先可用于观察不会发生混淆的血流多普勒频谱(步骤702)。可选地,首先可使用预置的PRF。
与上面情况类似,当计算多普勒最大频率fmax、最小多普勒频率fmin和平均多普勒频率fmean时,通过最大频率引擎953、最小频率引擎955和平均频率引擎957计算。图4表示多普勒功率谱,识别该频谱的已被计算出的最大频率fmax、最小频率fmin和平均频率fmean的值。对于每个频谱样本,其最大频率fmax、最小频率fmin和平均频率fmean的值像曲线一样被随时间追踪。
可如下首先从频谱925计算平均频率fmean作为一阶矩:
f mean = ∫ fpdf ∫ pdf - - - ( 9 )
此处f是频率,且p是多普勒频谱功率(或频谱振幅谱a,或提高到ab的功率,此处b是实数或从振幅导出的任何信号)。
在平均频率fmean从频谱被计算出来后,最大多普勒频率fmax和最小多普勒频率fmin也被计算出。
最大频率fmax和最小频率fmin被计算为从被计算的平均频率fmean值的频谱的百分位数。例如,从平均频率fmean,49.9%的最大频率fmax可在从平均频率fmean开始的正频率方向上被计算出来。相似地,最小频率fmin在负方向上被计算出来。
最大频率fmax和最小频率fmin一起设定了全部频谱功率的99.8%的合并边界,如下:
∫ f mean f max pdf ∫ pdf = 0.499 , 及                                                                (10)
∫ f min f mean pdf ∫ pdf = 0.499 - - - ( 11 )
因为平均频率fmean值是频谱的加权平均频率,只要百分值小于50%,那么最大频率fmax和最小频率fmin的值就通过使用公式(10)和(11)由最大引擎953和最小引擎955计算出来(步骤704)。可选地,最大频率和最小频率的值可分别使用公式(2)和(5)来计算。
图5A、图5B和图5C表示了随时间变化的被计算出的最大频率fmax501和最小频率fmin503的值。所述曲线设定了高的多普勒频谱边界high fmax505和低的多普勒频谱边界low fmin507。最大多普勒频率fmax曲线的最高值highfmax505和最小多普勒频率fmin曲线的最低值low fmin被捕捉并被记录。
如果最大频率fmax曲线或最小频率fmin曲线在观察期间经历混淆(如图1所示),那么混淆检测器和偏差引擎951通过将每个被卷绕的频率所经历的偏差加到其各自的被削减的峰值(clipped peak)而继续追踪最大频率fmax曲线和最小频率fmin曲线。如果在正的最大频率范围和负的最大频率这两个范围都检测到削减(clipping),那么当前PRF的设定就太小了。
如果所有频率分量都位于可能包括校正过的混淆频率(如果频谱曾被混淆)的正频率区域中,那么频谱就是单极正的。如果所有频率分量都位于可能包括校正过的混淆频率(如果频谱曾被混淆)的负频率区域中,那么频谱就是单极负的。在校正混淆(如果频谱曾被混淆)之后,如果频率分量位于正和负频率两个区域中,那么频谱就是双极的。
图5A表示了双极的频谱。最高最大频率high fmax505和最低最小频率low fmin507之间的频带宽度509被计算并被用于确定用于基于观察期间的最佳图像显示的新PRF。所述频带宽度为,
频带宽度=(high fmax)-(low fmin)           (12)
其可被认为用于记录被观察的血流的最小PRF。可添加频率安全余量fs1和fs2以调节确保在频谱和最大频率范围之间的余量充足的频带宽度509。
调节过的频带宽度=((hign fmax)-(low fmin))+fs1+fs2         (13)
调节过的频带宽度与当前PRF设定作比较(步骤706)。如果调节过的频带宽度大于当前PRF设定943,
调节过的频带宽度>当前PRF                        (14)
那么当前PRF设定943由处理器947增加到相应于调节过的频带宽度的设定,并且被输出到PRF发生器935(步骤718)。如果调节过的频带宽度小于当前PRF设定,那么可能不发生混淆但是当前PRF设定可能会太大。
进一步,将当前PRF设定的分数与调节过的频带宽度作比较,以将PRF减小到产生最佳图像显示的值。如果PRF设定对于血流速度太大,那么多普勒频谱925显示将会太小以至于不能准确表示血流速度。
当前PRF的分数被用作低级阈值。介于0和1之间的预定的分数,如
Figure G2008800105824D00131
可用作分数。
(分数)(当前PRF)<调节过的频带宽度<当前PRF    (15)
如果调节过的频带宽度小于分数PRF,那么多普勒频谱图像需要增大尺寸(步骤708)。因此,PRF 943被减小到调节过的频带宽度并且被输出到PRF发生器935(步骤716)。减小或增大PRF设定,直到调节过的频带宽度小于当前PRF设定而大于分数PRF。
图5B表示的是单极正频谱。此情形中,最高最大频率high fmax505加上频率安全余量fs1被用于确定新的PRF。将最高最大频率high fmax505加上频率安全余量fs1与当前正的最大频率界限b1fPRF作比较。如果最高最大频率high fmax505加上频率安全余量fs1大于当前正的最大频率界限b1fPRF 943,
(high fmax+fs1)>b1fPRF                (16)
那么当前PRF设定943由处理器947增大到对应于最高最大频率highfmax505加上频率安全余量fs1的设定,并且被输出到PRF发生器935。如果最高最大频率high fmax505加上频率安全余量fs1小于当前正的最大频率界限b1fPRF,那么可能不发生混淆但是当前PRF可能被设定得太大。
进一步,将最高最大频率high fmax505加上频率安全余量fs1与当前正的最大频率界限b1fPRF的分数作比较以减小PRF到产生最佳图像显示的值。如果对于血流速度的测量而言PRF设定太小,那么就将发生混淆。但是,如果对于血流速度而言PRF设定太大,那么多普勒频谱925显示将太小以至于不能准确表示血流速度。
正的低级阈值b2b1fPRF(此处b2是介于0和1之间的分数)被计算并且与最高最大频率high fmax505加上频率安全余量fs1作比较。
b2b1fPRF<(high fmax+fs1)                (17)
如果最高最大频率high fmax505加上频率安全余量fs1小于当前正的最大频率界限b1fPRF,那么多普勒频谱图像就需要增大尺寸。因此,PRF 943被减小到最高最大频率加上频率安全余量high fmax+fs1而且被输出到PRF发生器935。减小或者增大PRF设定,直到最高最大频率high fmax505加上频率安全余量fs1小于当前正的最大频率界限b1fPRF,同时大于正的低级阈值b2b1fPRF
图5C表示的是单极负频谱。此情形中,最低最小频率low fmin507加上频率安全余量fs2被用于确定新的PRF。将最低最小频率low fmin507加上频率安全余量fs2与当前负的最小频率界限-(1-b1)fPRF作比较。如果最低最小频率low fmin507的绝对值加上频率安全余量fs2大于当前负的最大频率界限-(1-b1)fPRF 943的绝对值,
(|low fmin|+fs2)>(1-b1)fPRF    (18)
当前PRF设定943由处理器947增加到对应于最低最小频率low fmin507的绝对值加上频率安全余量fs2的设定,并被输出到PRF发生器935。如果最低最小频率low fmin507的绝对值加上频率安全余量fs2小于当前负的最大频率界限-(1-b1)fPRF的绝对值,那么可能不发生混淆但是当前PRF设定可能会太大。
进一步,将最低最小频率low fmin507的绝对值加上频率安全余量fs2与当前负的最大频率界限-(1-b1)fPRF的绝对值的分数作比较,以将PRF减小到产生最佳图像显示的值。如果对于测量血流速度来说PRF设定太小,那么将会发生混淆。但是,如果对于血流速度来说PRF设定太大,那么多普勒频谱925显示将会太小以至于不能准确表示血流速度。
负的低级阈值-b2(1-b1)fPRF(此处b2是介于0和1之间的分数)被计算并且与最低最小频率low fmin507加上频率安全余量fs2作比较,最低最小频率low fmin507加上频率安全余量fs2又与当前负的最大频率界限-(1-b1)fPRF作比较。
b2(1-b1)fPRF<|lowfmin|+fs                (19)
如果最低最小频率low fmin507的绝对值加上频率安全余量fs2小于当前负的最大频率界限-(1-b1)fPRF的绝对值的分数,那么多普勒频谱图像就需要增大尺寸。因此,PRF 943被减小到最低最小频率low fmin507的绝对值加上频率安全余量fs2,并且被输出到PRF发生器935上。减小或者增大PRF设定,直到最低最小频率low fmin507的绝对值加上频率安全余量fs2小于当前负的最大频率界限-(1-b1)fPRF的绝对值,同时大于负的低级阈值-b2(1-b1)fPRF的绝对值。
如果不管频谱是双极的、还是正或负的单极,在调节PRF后检测到混淆,那么如上所述通过基线移动校正(步骤710、720、712、714)。因为频谱不必位于频率范围的中心,所以即使在PRF被确定的观察期间混淆没有发生,在调节PRF之后混淆还是可能发生。在减小PRF之后,高的最大频率或者低的最小频率可能超过相应的界限。
图8表示的流程图描述了调节PRF设定但没有进行基线移动的变量。基线可被固定在正的最大频率范围和负的最大频率范围之间的任何预定的位置。最初,PRF被设定为缺省的PRF值,或者最大PRF(步骤802)。在所述PRF处被输超声波,而且执行多普勒频谱925处理以得到多普勒频谱。
如上述公式(10)和(11),计算最大多普勒频率fmax和最小多普勒频率fmin。在观察期间(例如,至少一个心搏周期,心跳,或小于一个心搏周期)监控最大多普勒频率fmax和最小多普勒频率fmin,并且记录最大多普勒频率fmax曲线的最高值high fmax和最小多普勒频率fmin曲线的最低值low fmin
可以将频率安全余量fs1、fs2加到最高最大频率high fmax的绝对值和最低最小频率low fmin的绝对值,
|high fmax|+fs1,及               (20)
|low fmin|+fs2                    (21)
公式(20)和(21)被用于发现最佳PRF设定。
将最高最大频率high fmax加上频率安全余量fs1与最大正的频率界限b1fPRF作比较。如果最高最大频率high fmax加上频率安全余量fs1大于正的最大频率界限b1fPRF,那么将PRF增大到最高最大频率high fmax加上频率安全余量fs1的水平。相反地,将最低最小频率low fmin的绝对值加上频率安全余量fs2与负的最大频率界限-(1-b1)fPRF作比较。如果最低最小频率low fmin的绝对值加上频率安全余量fs2大于负的最大频率界限-(1-b1)fPRF的绝对值,那么将PRF增大到最低最小频率low fmin的绝对值加上频率安全余量fs2(步骤806,818)。
如果最高最大频率high fmax加上频率安全余量fs1小于正的最大频率界限b1fPRF,而且最低最小频率low fmin加上频率安全余量fs2小于负的最大频率界限-(1-b1)fPRF,那么将最高最大频率high fmax的绝对值与最低最小频率low fmin的绝对值作比较以确定哪侧的频率分量占优势(步骤808)。
所述比较确定正频率区域还是负频率区域占优势。如果
(|high fmax|+fs1)>(|low fmin|+fs2)        (22)
成立,那么正频率区域占优势而且计算正的低级阈值b2b1fPRF,此处b2是介于0和1之间的分数。
将最高最大频率high fmax加上频率安全余量fs1与正的低级阈值b2b1fPRF作比较(步骤820)
b2b1fPRF<(high fmax+fs1)       (23)
如果条件(23)满足,那么PRF设定就完成(步骤814)。如果最高最大频率high fmax加上频率安全余量fs1(12)小于低级阈值b2b1fPRF,那么将PRF减小以满足条件(23),与此同时维持在负频率范围没有混淆发生(步骤816)。如果混淆开始发生,甚至在满足上述条件(23)之前停止减小PRF。
如果条件(22)不满足,那么负频率区域占优势并且计算负的低级阈值-b2(1-b1)fPRF(步骤808)。
最低最小频率low fmin的绝对值加上频率安全余量fs2与负的低级阈值-b2(1-b1)fPRF作比较(步骤822)。
(|low fmin|+fs2)>b2(1-b1)fPRF    (24)
如果条件(24)满足,那么PRF设定就完成(步骤814)。如果最低最小频率low fmin的绝对值加上频率安全余量fs2小于低级阈值-b2(1-b1)fPRF的绝对值,那么将PRF减小以满足条件(24),与此同时在正频率一侧维持没有混淆发生。如果混淆开始发生,那么甚至在满足上述条件(24)之前停止减小PRF。
一个测试确定最高最大频率high fmax加上频率安全余量fs1是否大于用于混淆的最大正频率界限b1fPRF,或者,最低最小频率low fmin的绝对值加上频率安全余量fs2是否大于用于混淆的最小负频率界限-b2(1-b1)fPRF的绝对值。
如果最高最大频率high fmax加上频率安全余量fs1小于最大正频率界限b1fPRF,而且,如果最低最小频率low fmin的绝对值加上频率安全余量fs2小于最小负频率界限-(1-b1)fPRF的绝对值,那么就进行另一个测试。
另一个测试确定,如果正频率占优势(或者说条件(22)成立),那么最高最大频率high fmax加上频率安全余量fs1是否大于正的低级阈值b2b1fPRF;或者,如果负频率占优势(或者说条件(22)不成立),那么最低最小频率low fmin的绝对值加上频率安全余量fs2是否大于负的低级阈值-b2(1-b1)fPRF的绝对值。所述测试确保多普勒频谱对于显示足够大。如果PRF太高,那么多普勒频谱显示受损,而且这对于准确的医疗诊断来说是不可接受的。在所述变量中,基线103被固定且不能作基线移动。
因为基线没有被移动,所以在没有占优势的频率区域中减小的PRF可能在频谱中导致混淆。例如,如果正频率占优势,那么上述有条件的测试就基于正频率最大值调节当前PRF而且相应地调节PRF。减小PRF的过程中,联合了频谱的负的部分可能开始被混淆。当频谱的负的部分开始混淆,停止减小PRF。
已经描述了本发明的一个或多个实施例。不过,可以理解的是,在不脱离本发明的精神和范围的情况下,可以实施各种改变。因此,其它实施例包括于随后的权利要求的范围内。

Claims (24)

1、一种用于检测和校正多普勒频谱中的混淆的方法,其包括:
随时间接收多普勒频谱信号;
由所述多普勒频谱计算最大频率fmax和最小频率fmin
随时间追踪所述最大频率fmax和最小频率fmin
如果正频率区域中的频率改变(卷绕)到负频率区域,则检测混淆是否发生自所述最大频率fmax,或如果负频率区域中的频率改变(卷绕)到正频率区域,则检测混淆是否发生自所述最小频率fmin;及
如果检测到混淆,则根据最大频率偏差fa在正或负方向上对分隔多普勒频谱的负频率区域和正频率区域的零频率基线进行移动。
2、如权利要求1所述的方法,其中,所述多普勒频谱信号来自由振幅谱a、功率谱a2或提高到ab的功率组成的组,这里b是实数。
3、如权利要求1所述的方法,其中,所述最大频率fmax是所述多普勒频谱的百分值。
4、如权利要求1所述的方法,其中,所述最小频率fmin是所述多普勒频谱的百分值。
5、如权利要求1所述的方法,进一步包括由自一个(正或负的)区域卷绕到另一个(负或正的)区域的频率的大小确定最大频率偏差fa
6、如权利要求1所述的方法,进一步包括:
将频率安全余量fs加到所述最大频率偏差fa;及
将所述频谱基线移动所述最大频率偏差加上安全余量±(fa+fs)。
7、一种用于确定超声波系统的脉冲重复频率的方法,其包括:
随时间接收多普勒频谱信号;
由所述多普勒频谱计算最大频率fmax
由所述多普勒频谱计算最小频率fmin
随时间追踪所述最大频率fmax和最小频率fmin
捕捉被追踪的所述最大频率fmax的最高值high fmax和所述最小频率fmin的最低值low fmax
将所述最高值high fmax和所述最低值low fmin作比较,以确定所述最大频率fmax和最小频率fmin是双极、负的单极还是正的单极的;
如果是双极的情况:
基于所述最高最大频率high fmax和所述最低最小频率low fmin之间的差而确定频带宽度;
将所述频带宽度与当前脉冲重复频率设定值作比较;
如果所述频带宽度大于所述当前脉冲重复频率设定值,则增大所述脉冲重复频率设定值;
如果所述频带宽度小于预定分数的当前脉冲重复频率设定值,则减小所述脉冲重复频率设定值;及
如果所述频带宽度小于所述当前脉冲重复频率设定值且大于所述预定分数的当前脉冲重复频率设定值,则使用所述当前脉冲重复频率设定值;
如果是正的单极的情况;
将所述最高最大频率high fmax与当前正的最大频率界限b1fPRF作比较,其中如果所述最高最大频率high fmax大于所述当前正的最大频率界限b1fPRF,则将当前脉冲重复频率设定值增大到对应于所述最高最大频率highfmax的设定;
如果所述最高最大频率high fmax小于当前正的最大频率界限b1fPRF,则将所述最高最大频率high fmax与低级阈值b2b1fPRF作比较,其中如果所述最高最大频率high fmax小于所述低级阈值b2b1fPRF,则减小所述脉冲重复频率直到等于所述最高最大频率high fmax;及
如果是负的单极的情况:
将所述最低最小频率low fmin的绝对值与当前负的最大频率界限-(1-b1)fPRF的绝对值作比较,其中如果所述最低最小频率low fmin的绝对值大于所述当前负的最大频率界限-(1-b1)fPRF的绝对值,则将所述当前脉冲重复频率设定值增大到对应于所述最低最小频率low fmin的绝对值的设定;
如果所述最低最小频率low fmin的绝对值小于所述当前负的最大频率界限-(1-b1)fPRF的绝对值,则将所述最低最小频率low fmin的绝对值与低级阈值-b2(1-b1)fPRF的绝对值作比较,其中如果所述最低最小频率low fmin的绝对值小于所述低级阈值-b2(1-b1)fPRF的绝对值,则减小所述脉冲重复频率到等于所述最低最小频率low fmin的绝对值。
8、如权利要求7所述的方法,其中,所述多普勒频谱信号来自由振幅谱a、功率谱a2或提高到ab的功率组成的组,这里b是实数。
9、如权利要求7所述的方法,其中,所述最大频率fmax是所述多普勒频谱的百分位数频率。
10、如权利要求7所述的方法,其中,所述最小频率fmin是所述多普勒频谱的百分位数频率。
11、如权利要求7所述的方法,进一步包括:
如果在正频率区域中的频率卷绕到负频率区域,则检测混淆是否发生自所述最大频率fmax
如果在所述负频率区域中的负频率卷绕到所述正频率区域,则检测混淆是否发生自所述最小频率fmin;以及
如果检测到混淆,则根据最大频率偏差fa在正或负方向上对分隔多普勒频谱的所述负频率区域和正频率区域的零频率基线进行移动。
12、如权利要求7所述的方法,进一步包括:
如果在正频率区域中的频率卷绕到负频率区域,则检测混淆是否发生自所述最大频率fmax
如果在所述负频率区域中的频率卷绕到所述正频率区域,则检测混淆是否发生自所述最小频率fmin;以及
如果检测到混淆,则将被混淆的频率偏差加到所述被混淆的最大频率以及加到被混淆的最小频率。
13、如权利要求11所述的方法,其中,所述频率偏差还包括:
由自一个(正或负的)区域卷绕到另一个(负或正的)区域的频率的大小确定所述最大偏差fa,以计算所述最大频率fmax和所述最小频率fmin
14、如权利要求13所述的方法,进一步包括:
将频率安全余量加到所述最大频率偏差fa;及
将所述频谱基线移动所述最大频率偏差fa加安全余量。
15、如权利要求7所述的方法,进一步包括:
将频率安全余量加到所述频带宽度、最高最大频率high fmax或最低最小频率low fmin的绝对值;及
与所述当前脉冲重复频率作比较。
16、一种确定用于超声波系统的脉冲重复频率的方法,其包括:
设定初始脉冲重复频率;
随时间接收多普勒频谱信号;
由所述多普勒频谱计算最大频率fmax
由所述多普勒频谱计算最小频率fmin
随时间追踪所述最大频率fmax和最小频率fmin
捕捉被追踪的所述最大频率fmax的最高值high fmax和所述最小频率fmin的最低值low fmin
将所述最高最大值high fmax的绝对值与所述最低最小频率low fmin的绝对值作比较,以确定所述正频率区域和负频率区域哪个占优势;
如果所述最高最大值high fmax较大,则所述正频率区域占优势且计算正的低级阈值b2b1fPRF;及
将所述最高最大频率high fmax与所述正的最大频率界限b1fPRF和所述正的低级阈值b2b1fPRF作比较,其中如果所述最高最大频率high fmax小于所述正的低级阈值b2b1fPRF,则执行减小所述PRF直到所述正的最大频率界限b1fPRF等于所述最高最大频率high fmax,或者混淆开始发生于所述负的最大频率界限-(1-b1)fPRF处,这里执行先到来的那个,而且其中如果所述最高最大频率high fmax大于所述正的最大频率界限b1fPRF,则增大所述脉冲重复频率直到等于所述最高最大频率high fmax
如果所述最低最小频率low fmin的绝对值较大,则所述负频率区域占优势并且计算负的低级阈值-b2(1-b1)fPRF
将所述最低最小频率low fmin的绝对值与所述负的最大频率界限-(1-b1)fPRF的绝对值和所述低级阈值-b2(1-b1)fPRF作比较,其中如果所述最低最小频率low fmin的绝对值小于低级阈值-b2(1-b1)fPRF的绝对值,则执行减小所述脉冲重复频率直到所述负的最大频率界限-(1-b1)fPRF的绝对值为所述最低最小频率low fmin的绝对值,或者混淆开始发生于所述正的频率界限处,这里执行先到来的那个,而且其中如果所述最低最小频率low fmin的绝对值大于所述负的最大频率界限-(1-b1)fPRF的绝对值,则增大所述脉冲重复频率直到等于所述最低最小频率low fmin
17、如权利要求16所述的方法,进一步包括将频率安全余量加到所述最高最大频率high fmax和所述最低最小频率low fmin的绝对值。
18、如权利要求16所述的方法,其中,所述正的最大频率界限b1fPRF和所述负的最大频率界限-(1-b1)fPRF由所述脉冲重复频率和固定的零频率基线位置所确定。
19、如权利要求16所述的方法,其中,所述最大频率fmax是所述多普勒频谱的百分位数频率。
20、如权利要求16所述的方法,其中,所述最小频率fmin是所述多普勒频谱的百分位数频率。
21、如权利要求16所述的方法,其中,观察期能够小于心搏周期,或小于至少一个心搏周期的长周期。
22、如权利要求16所述的方法,其中,由带有或不带有降噪的多普勒频谱计算所述最大频率fmax和所述最小频率fmin
23、一种用于检测和校正多普勒频谱中的混淆的系统,其包括:
用于随时间接收多普勒频谱信号的装置;
用于由所述多普勒频谱计算最大频率fmax和最小频率fmin的装置;
用于随时间追踪所述最大频率fmax和最小频率fmin的装置;
如果在正频率区域中的频率改变(卷绕)到负频率区域,检测混淆是否发生自所述最大频率fmax的装置;
如果在所述负频率区域中的频率改变(卷绕)到所述正频率区域,检测混淆是否发生自所述最小频率fmin的装置;以及
如果检测到混淆,根据最大频率偏差fa对在正或负方向上分隔所述多普勒频谱的负频率区域和正频率区域的零频率基线进行移动的装置。
24、一种用于确定用于超声波系统的脉冲重复频率的系统,其包括:
用于设定初始脉冲重复频率的装置;
用于随时间接收多普勒频谱信号的装置;
用于由所述多普勒频谱计算最大频率fmax的装置;
用于由所述多普勒频谱计算最小频率fmin的装置;
用于随时间追踪所述最大频率fmax和最小频率fmin的装置;
用于捕捉被追踪的所述最大频率fmax的最高值high fmax和所述最小频率fmin的最低值low fmin的装置;
用于将所述最高最大频率high fmax的绝对值与所述最低最小频率low fmin的绝对值作比较,以确定所述正频率区域和负频率区域哪个占优势的装置;
如果所述最高最大频率high fmax较大,则所述正频率区域占优势而且计算正的低级阈值b2b1fPRF;及
用于将所述最高最大频率high fmax与所述正的最大频率界限b1fPRF和所述正的低级阈值b2b1fPRF作比较的装置,其中如果最高最大频率high fmax小于所述正的低级阈值b2b1fPRF,则执行减小脉冲重复频率直到所述正的最大频率界限b1fPRF等于所述最高最大频率high fmax,或者混淆开始发生于所述负的最大频率界限-(1-b1)fPRF处,这里执行先到来的那个,而且其中如果所述最高最大频率high fmax大于所述正的最大频率界限b1fPRF,则增大所述脉冲重复频率到等于所述最高最大频率high fmax
如果所述最低最小频率low fmin的绝对值较大,则所述负的频率区域占优势而且计算低级阈值-b2(1-b1)fPRF
用于将所述最低最小频率low fmin的绝对值与所述负的最大频率界限-(1-b1)fPRF的绝对值和所述低级阈值-b2(1-b1)fPRF的绝对值作比较的装置,其中如果所述最低最小频率low fmin的绝对值小于所述低级阈值-b2(1-b1)fPRF的绝对值,则执行减小所述脉冲重复频率直到所述负的最大频率界限-(1-b1)fPRF的绝对值为所述最低最小频率low fmin的绝对值,或者混淆开始发生于所述正频率界限处,这里执行先到来的那个,而且其中如果所述最低最小频率low fmin的绝对值大于所述负的最大频率界限-(1-b1)fPRF的绝对值,则增大所述脉冲重复频率到等于所述最低最小频率low fmin
CN2008800105824A 2007-03-29 2008-03-28 用于超声波成像技术领域的方法和仪器 Expired - Fee Related CN101646390B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US92108907P 2007-03-29 2007-03-29
US60/921,089 2007-03-29
US11/926,251 US20080242994A1 (en) 2007-03-29 2007-10-29 Methods and apparatus for ultrasound imaging
US11/926,251 2007-10-29
PCT/JP2008/056728 WO2008123596A1 (en) 2007-03-29 2008-03-28 Methods and apparatus for ultrasound imaging

Publications (2)

Publication Number Publication Date
CN101646390A true CN101646390A (zh) 2010-02-10
CN101646390B CN101646390B (zh) 2012-03-21

Family

ID=39795586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008800105824A Expired - Fee Related CN101646390B (zh) 2007-03-29 2008-03-28 用于超声波成像技术领域的方法和仪器

Country Status (5)

Country Link
US (1) US20080242994A1 (zh)
EP (1) EP2129295A4 (zh)
JP (2) JP5385790B2 (zh)
CN (1) CN101646390B (zh)
WO (1) WO2008123596A1 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102551791A (zh) * 2010-12-17 2012-07-11 深圳迈瑞生物医疗电子股份有限公司 一种超声成像方法和装置
CN102764140A (zh) * 2012-08-16 2012-11-07 无锡祥生医学影像有限责任公司 一种用于触摸屏超声诊断仪的多普勒频谱优化方法及其装置
CN105476663A (zh) * 2016-01-13 2016-04-13 北京悦琦创通科技有限公司 调整多普勒频谱图的方法和设备
CN105476661A (zh) * 2016-01-13 2016-04-13 北京悦琦创通科技有限公司 调整多普勒频谱图的方法和设备
CN108472013A (zh) * 2016-01-22 2018-08-31 奥林巴斯株式会社 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序
CN110477955A (zh) * 2019-08-22 2019-11-22 电子科技大学 一种基于iq数据的血管自动识别方法
CN111265251A (zh) * 2020-01-22 2020-06-12 青岛海信医疗设备股份有限公司 一种多普勒谱校正方法和设备
CN111329522A (zh) * 2018-12-19 2020-06-26 株式会社日立制作所 超声波摄像装置及其控制方法

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8142361B2 (en) * 2007-08-28 2012-03-27 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method and apparatus for automatic optimization of doppler imaging parameters
KR100969537B1 (ko) * 2008-01-10 2010-07-12 주식회사 메디슨 도플러 모드 영상을 형성하는 초음파 시스템 및 방법
CN101884551B (zh) * 2009-05-15 2014-10-15 深圳迈瑞生物医疗电子股份有限公司 提高超声多普勒成像自动调整性能的方法及其超声系统
US8882675B2 (en) 2009-07-16 2014-11-11 Hitachi Aloka Medical, Ltd. Methods and apparatus for ultrasound imaging
CN102469982B (zh) * 2009-08-11 2014-04-23 日立阿洛卡医疗株式会社 用于超声波成像的方法及装置
CN101991437B (zh) * 2009-08-24 2013-09-04 深圳迈瑞生物医疗电子股份有限公司 一种信号处理方法和装置及多普勒超声系统
US9155519B1 (en) 2010-04-05 2015-10-13 Hitachi Aloka Medical, Ltd. Methods and apparatus for ultrasound imaging
JP5950271B2 (ja) * 2011-07-27 2016-07-13 株式会社日立製作所 超音波診断装置
JP6391912B2 (ja) * 2013-02-26 2018-09-19 キヤノンメディカルシステムズ株式会社 超音波診断装置
JP6039088B2 (ja) 2013-09-09 2016-12-07 パイオニア株式会社 流速検出装置及び流速検出方法
JP6515205B2 (ja) * 2016-01-04 2019-05-15 株式会社日立製作所 超音波撮像装置
CN106236148B (zh) * 2016-07-29 2019-06-28 飞依诺科技(苏州)有限公司 一种脉冲重复频率确定方法及装置
US11213278B2 (en) 2017-05-19 2022-01-04 Mayo Foundation For Medical Education And Research System and method for visualization of tissue microvascular using ultrasound
WO2019064471A1 (ja) * 2017-09-29 2019-04-04 本多電子株式会社 超音波診断装置及びそれにおける信号処理方法、信号処理プログラムを格納した記録媒体
JP7334494B2 (ja) * 2019-06-25 2023-08-29 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0651035B2 (ja) * 1985-08-20 1994-07-06 松下電器産業株式会社 超音波パルスドツプラ血流計
US5152291A (en) * 1990-09-07 1992-10-06 Hewlett-Packard Company Acoustic fiber measurement of intravascular blood
JP2693575B2 (ja) * 1989-05-24 1997-12-24 株式会社東芝 超音波ドプラ血流測定装置
JP3698173B2 (ja) * 1995-05-24 2005-09-21 東芝医用システムエンジニアリング株式会社 超音波診断装置
JP3892538B2 (ja) * 1997-07-18 2007-03-14 株式会社東芝 超音波ドプラ診断装置
JP2000005177A (ja) * 1998-06-17 2000-01-11 Matsushita Electric Ind Co Ltd 超音波ドプラ血流計
JP4740436B2 (ja) * 1998-12-31 2011-08-03 ゼネラル・エレクトリック・カンパニイ ドップラー超音波スペクトル写真のための速度スケールおよびパルス繰返し数の自動調整
KR100381874B1 (ko) * 2000-07-08 2003-04-26 주식회사 메디슨 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
US6663566B2 (en) * 2002-02-19 2003-12-16 Ge Medical Systems Global Technology Company, Llc Method and apparatus for automatic control of spectral doppler imaging
JP2004194911A (ja) * 2002-12-18 2004-07-15 Ge Medical Systems Global Technology Co Llc 超音波診断装置
JP2006520619A (ja) * 2003-02-13 2006-09-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 超音波カラーフロードプラ情報から合成されるフロースペクトログラム
US6733454B1 (en) * 2003-02-26 2004-05-11 Siemens Medical Solutions Usa, Inc. Automatic optimization methods and systems for doppler ultrasound imaging
US6976960B2 (en) * 2003-06-10 2005-12-20 Siemens Medical Solutions Usa, Inc. Automatic velocity anti-aliased ultrasound methods and systems
JP4458825B2 (ja) * 2003-11-27 2010-04-28 株式会社日立メディコ 超音波診断装置
US7288068B2 (en) * 2003-12-15 2007-10-30 Siemens Medical Solutions Usa, Inc. Automatic optimization for ultrasound medical imaging
KR20060124824A (ko) * 2005-05-26 2006-12-06 주식회사 메디슨 초음파 스펙트럼 영상을 처리하는 방법 및 초음파 진단시스템

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102551791A (zh) * 2010-12-17 2012-07-11 深圳迈瑞生物医疗电子股份有限公司 一种超声成像方法和装置
CN102764140A (zh) * 2012-08-16 2012-11-07 无锡祥生医学影像有限责任公司 一种用于触摸屏超声诊断仪的多普勒频谱优化方法及其装置
CN102764140B (zh) * 2012-08-16 2014-07-02 无锡祥生医学影像有限责任公司 一种用于触摸屏超声诊断仪的多普勒频谱优化方法及其装置
CN105476663A (zh) * 2016-01-13 2016-04-13 北京悦琦创通科技有限公司 调整多普勒频谱图的方法和设备
CN105476661A (zh) * 2016-01-13 2016-04-13 北京悦琦创通科技有限公司 调整多普勒频谱图的方法和设备
CN108472013A (zh) * 2016-01-22 2018-08-31 奥林巴斯株式会社 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序
CN111329522A (zh) * 2018-12-19 2020-06-26 株式会社日立制作所 超声波摄像装置及其控制方法
CN110477955A (zh) * 2019-08-22 2019-11-22 电子科技大学 一种基于iq数据的血管自动识别方法
CN110477955B (zh) * 2019-08-22 2021-05-11 电子科技大学 一种基于iq数据的血管自动识别方法
CN111265251A (zh) * 2020-01-22 2020-06-12 青岛海信医疗设备股份有限公司 一种多普勒谱校正方法和设备
CN111265251B (zh) * 2020-01-22 2023-01-24 青岛海信医疗设备股份有限公司 一种多普勒谱校正方法和设备

Also Published As

Publication number Publication date
JP2013116372A (ja) 2013-06-13
US20080242994A1 (en) 2008-10-02
EP2129295A4 (en) 2012-07-11
JP5385790B2 (ja) 2014-01-08
EP2129295A1 (en) 2009-12-09
WO2008123596A1 (en) 2008-10-16
JP2010522581A (ja) 2010-07-08
JP5462966B2 (ja) 2014-04-02
CN101646390B (zh) 2012-03-21

Similar Documents

Publication Publication Date Title
CN101646390B (zh) 用于超声波成像技术领域的方法和仪器
US7850609B2 (en) Ultrasound diagnostic apparatus
US9895138B2 (en) Ultrasonic diagnostic apparatus
JP3746115B2 (ja) 超音波診断装置
US7513872B2 (en) Ultrasonic doppler measuring apparatus and control method therefor
US4961427A (en) Ultrasonic doppler diagnostic apparatus
WO2014013839A1 (ja) 超音波診断装置及び画像処理装置
CN106333707B (zh) 超声多普勒图谱校正方法、装置及超声诊断系统
CN105534549B (zh) 一种超声骨密度仪探头位置监测系统及其监测方法
US20110092819A1 (en) Ultrasonic diagnosis apparatus and ultrasoinc data acquisition method
JP2010259799A (ja) クラッタ信号をフィルタリングするための固有ベクトル設定方法およびこれを用いた超音波システム
KR20130081626A (ko) 벡터 도플러를 이용한 초음파 장치 및 초음파 영상 생성 방법
JPH06327672A (ja) 超音波ドプラ診断装置
JP3578680B2 (ja) 超音波診断装置
US8189427B2 (en) Clutter signal filtering for doppler signal
JP2008237647A (ja) 超音波診断装置、ドプラ計測装置、及びドプラ計測プログラム
JP4357260B2 (ja) 加速度脈波測定装置
JPH05337111A (ja) 超音波診断装置
JPS6237980B2 (zh)
JP5909076B2 (ja) 超音波診断装置及び制御プログラム
KR101221288B1 (ko) 적응적 클러터 필터링을 수행하는 초음파 시스템 및 방법
JP5235110B2 (ja) 超音波診断装置
JP2007054226A (ja) 超音波診断装置
Stegman et al. A feasibility study for measuring accurate tendon displacements using an audio-based Fourier analysis of pulsed-wave Doppler ultrasound signals
JPH0556973A (ja) 超音波血流イメージング装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB02 Change of applicant information

Address after: Tokyo, Japan

Applicant after: Hitachi Aloka Medical Ltd.

Address before: Tokyo, Japan

Applicant before: Arockar K. K.

COR Change of bibliographic data

Free format text: CORRECT: APPLICANT; FROM: AROCKAR K.K. TO: HITACHI ALOKA MEDICAL LTD.

C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20170216

Address after: Tokyo, Japan, Japan

Patentee after: Hitachi Ltd.

Address before: Tokyo, Japan

Patentee before: Hitachi Aloka Medical Ltd.

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120321

Termination date: 20190328

CF01 Termination of patent right due to non-payment of annual fee