CN101347343A - 超声波图像取得装置 - Google Patents
超声波图像取得装置 Download PDFInfo
- Publication number
- CN101347343A CN101347343A CNA2008101286852A CN200810128685A CN101347343A CN 101347343 A CN101347343 A CN 101347343A CN A2008101286852 A CNA2008101286852 A CN A2008101286852A CN 200810128685 A CN200810128685 A CN 200810128685A CN 101347343 A CN101347343 A CN 101347343A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- flow velocity
- velocity waveforms
- actual measurement
- waveforms
- 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
- 238000001228 spectrum Methods 0.000 claims abstract description 43
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 238000002604 ultrasonography Methods 0.000 claims abstract description 16
- 239000012530 fluid Substances 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims description 72
- 238000003384 imaging method Methods 0.000 claims description 59
- 230000017531 blood circulation Effects 0.000 claims description 47
- 238000003860 storage Methods 0.000 claims description 17
- 238000003745 diagnosis Methods 0.000 claims description 14
- 230000005856 abnormality Effects 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 8
- 230000003595 spectral effect Effects 0.000 claims description 7
- 230000000737 periodic effect Effects 0.000 claims description 4
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 239000008280 blood Substances 0.000 abstract description 10
- 210000004369 blood Anatomy 0.000 abstract description 10
- 230000000694 effects Effects 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 25
- 230000006870 function Effects 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 230000009471 action Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 210000004204 blood vessel Anatomy 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000010247 heart contraction Effects 0.000 description 2
- 210000003709 heart valve Anatomy 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 210000001715 carotid artery Anatomy 0.000 description 1
- 230000004087 circulation Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000803 paradoxical effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Hematology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (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
本发明提供一种超声波图像取得装置,其中,发送接收部向诊断部位发送接收具有与速度范围对应的重复频率的超声波。多普勒信号处理部根据通过超声波的发送接收而得到的信号,生成表示运动流体的速度的多普勒频谱图像。存储部预先存储有基于使基准的血流速度波形与ECG波形相关联的模型的模型化后的值。实测波形计算部根据特定患者的频谱图像来求出实测的血流速度波形。ECG波形取得部取得与实测的血流速度波形对应的定时的ECG波形。波形推定部根据实测的血流速度波形、ECG波形、以及模型化后的值来推定患者的血流速度波形。插值部对实测的血流速度波形和推定的血流速度波形进行比较,从实测的血流速度波形中除去其差超过阈值的部分,对该部分进行插值。
Description
技术领域
本发明涉及一种生成多普勒频谱图像的超声波图像取得装置。更具体地说涉及一种自动地调整多普勒频谱图像的速度范围的超声波图像取得装置。
背景技术
以往,如下的超声波图像取得装置是已知的:并用超声波脉冲反射法和超声波多普勒法,通过使用了一个超声波探测器的超声波操作来取得诊断部位的断层像和其血流信息,并且实时地至少显示该血流信息。在该超声波图像取得装置中,由于对体内的有血液等的流的诊断部位发送接收的超声波的多普勒效应,接收频率相对于发送频率稍微偏移,根据该偏移频率(多普勒偏移频率)与血流速度成比例的超声波多普勒法的原理,来进行多普勒偏移频率的频率解析,根据其结果取得血流信息。
在上述的超声波图像取得装置中,将多普勒频率的频谱图像作为对象来进行在诊断中使用的项目(参数)的计测处理,该多普勒频率的频谱图像是以频率f为纵轴、以时间t为横轴、将角频率分量的功率(强度)作为亮度(灰度),将对所取得的多普勒信号的快速傅立叶变换(FFT)的频率分析的结果进行频谱显示而得到的。
按顺序对该计测处理的流程进行说明。(1)在多普勒频率的频谱图像上,求出与该频率f轴方向的频率分布内的最大频率对应的最大流速Vp(Vpeak)以及与平均频率对应的平均速度Vm(Vmean)的位置。(2)在该时间轴t方向上描绘出该最大流速Vp、平均流速Vm的1个时间变化。(3)在示出该Vp、Vm的时间性位置变化曲线的描绘波形上,在每个心周期(1次心跳)进行心脏收缩期的波形峰值PS(Peak of Systolic,心脏收缩峰值)以及心脏舒张期的波形峰值ED(End of Diastolic,舒张期末)的同步检测。(4)根据该PS·ED的信息,执行如下处理:计测血管内的血流量、脉动流的HR(HeartRate:心率)、PI(Pulsatility Index,脉动指数)、以及RI(ResistanceIndex,阻力指数)等与诊断相关的各种参数(指标),并显示其计测值(参数计测处理)。
上述的Vp、Vm的描绘波形检测处理、PS·ED的峰值检测处理、以及PI、RI等参数计测处理是以凝固图像为对象基本进行人工操作。而且,近年来,以实时图像为对象并以自动操作进行的超声波图像取得装置也得到了普及。
但是,在脉冲多普勒(PW)法中,发送具有预定的重复频率的脉冲,利用预定的采样频率来进行所接收到的信号的频率解析。在用于该频率解析的采样频率fs低于多普勒偏移频率的情况下,引起混迭(Aliasing)现象(折叠)。因此,为了防止该现象,需要提高脉冲的重复频率(RPF:Pulse Repetition Frequency),缩短每次的观测时间间隔。在该情况下,如果指定了所希望测定的位置,则必然地决定了最大的PRF,如果决定了PRF,则也决定了可以计测的最高血流速度。将该可以计测的最高血流速度称为速度范围。
例如,在希望对30cm左右的血流速度进行测量时,如果设为10cm左右的速度范围,则发生混迭,而无法对血流进行测量。因此,在该情况下,需要将速度范围设定为50cm左右。
在进行了多普勒频谱显示的情况下,如果速度范围过小,则如上述产生折叠部分。在这样的情况下,通过操作者手动地将多普勒速度范围设定得较大,而使折叠部分落于奈奎斯特频率(Nyquist rate)(PRF的一半)内,从而得到在显示上连接良好的多普勒频谱图像。
相反,如果速度范围过大,则频谱的波形变小而难以观测。在这样的情况下,操作者通过将速度范围设定得较小,可以得到有效地使用了显示画面的上下的易于观测的多普勒频谱图像。
另外,在超声波多普勒法中,对血流的方向流向超声波探测器的血流附加正的符号。另外,对远离超声波探测器的血流附加负的符号。如果向某血管施加超声波探测器,则在该血管为动脉的情况下血流的速度随着脉动而变动,但跨越正负的变化少,而通常偏向正负的某一侧。
例如,在进行多普勒频谱显示时产生了折叠部分的情况下,操作者有时对基线移位开关进行操作来使多普勒频谱图像的基线(BL=0)移位。将其称为速度偏差的调整。通过使该基线移位-0.25(基线移位量=-0.25),折叠部分超过奈奎斯特频率地移动,从而得到在显示上连接良好的多普勒频谱。
此处,参照图1来对超声波图像取得装置取得的多普勒频谱图像进行说明。图1是示出由超声波图像取得装置取得的断层像和多普勒频谱图像的图。此处,对将颈动脉作为诊断部位而取得并显示断层像和多普勒频谱图像的情况进行说明。
例如,在画面110中,如果利用用于指定取得多普勒频谱图像的位置的范围门(range gate)102来指定了重叠有B模式断层像100和彩色多普勒图像101的图像上所显示出的血管,则可以得到表示该位置处的血流速度分布的时间变化的多普勒频谱图像并显示在画面上。在画面111中,显示出决定速度范围(可以测定的最高血流速度)的PRF为7.1(kHz)、速度偏差(BLS:Baseline Shift,基线移位)为0(Hz)时的多普勒频谱图像103(画面110的用虚线表示的部分)。
进而,操作者为了根据多普勒频谱图像103的自动描绘的峰值形状来进行血流状态的计测,以在速度范围(纵轴)的中心以一定的比率显示的方式,在画面120中调整脉冲重复频率(PRF)和速度偏差(BLS)。例如,通过将PRF变更成5(kHz),并将BLS移位到负侧,来如画面121所示放大多普勒频谱图像103来进行显示(画面120的用虚线表示的部分)。
在使用超声波图像取得装置来对血流速度等进行测定时,由于被检者的疾病和身体状况、探测器的施加方式(角度)、利用PW多普勒时血管内的范围门位置和宽度、诊断部位的不同,要测定的血流速度较大地变化。因此,以往操作者每次都进行调整装置的速度范围或者使基线移位那样的最优化,根据放大后的波形,来进行HR、PI、RI的计测。但是,象这样每当要诊断的血流的状态变化时,都调整与速度范围对应的PRF和速度偏差(BLS)的作法很麻烦。
因此,提出了如下的超声波图像取得装置:使多普勒频谱图像的速度范围以及速度偏差的调整操作自动化,从而使血流计测的操作性提高,操作者无需注意装置设定就可以进行诊断(例如,日本特开2005-185731号公报)。在以往技术的方法中,通过根据在预定时间(例如1次心跳)内取得的多普勒波形计算出角速度(频率)的频度来生成表示速度的频度分布的直方图,根据该直方图,决定速度范围来进行反馈,以在显示区域的纵向的α%(例如70%)以内显示出多普勒波形。
这里存在的课题是,当自动地调整多普勒频谱图像的速度范围时,多普勒频谱的计测的稳定性和可靠性变得重要。但是,在循环器(心脏)中存在用于使血流等循环的阀,产生该阀动作时的信号(以下称为“阀信号”)。因此,在以往的超声波图像取得装置中,特别在循环器(心脏)的多普勒血流诊断中,除了血流信号以外还混入高功率的阀信号,从而难以进行仅以血流成分为对象的自动的速度范围调整。
发明内容
本发明的目在于提供一种超声波图像取得装置,推定与期望的心时相相当、且没有受到阀信号的波形的影响的波形,根据该推定的波形从实测的血流速度波形中自动地除去心时相中的阀信号混入区域,通过对该除去后的部分进行插值来生成除去了阀信号的影响的实测的血流波形。
另外,本发明的目的在于提供一种超声波图像取得装置,根据所推定的波形从实测的血流速度波形中自动地除去心时相中的阀信号混入区域,根据除去了该阀信号后的实测的血流速度波形来进行自动的速度范围调整。
本发明的第1方式的超声波图像取得装置具有发送接收部、多普勒信号处理部、存储部、实测波形计算部、ECG波形取得部、波形推定部、以及插值部。而且,发送接收部向包含被检体内的运动流体的诊断部位发送接收具有与表示可测定的速度的速度范围对应的重复频率的超声波。多普勒信号处理部根据通过上述超声波的发送接收而得到的信号,生成表示上述运动流体的速度的多普勒频谱图像。存储部预先存储有基于使基准的血流速度波形与ECG波形相关联的模型的模型化后的值。实测波形计算部根据特定的患者的上述频谱图像来求出实测的血流速度波形。ECG波形取得部取得与上述实测的血流速度波形对应的定时的ECG波形。波形推定部根据上述实测的血流速度波形、上述ECG波形、以及上述模型化后的值来推定上述患者的血流速度波形。插值部对上述实测的血流速度波形和上述推定的血流速度波形进行比较,从上述实测的血流速度波形中除去其差超过阈值的部分,对上述实测的血流速度波形的该除去的部分进行插值。
根据该第1方式,可以自动地生成除去了阀信号的影响的血流速度波形。由此,可以生成不受阀信号影响的易于描绘的多普勒频谱图像,可以使血流量等的计测精度提高。
另外,本发明的第2方式的超声波图像取得装置具有发送接收部、多普勒信号处理部、存储部、实测波形计算部、ECG波形取得部、波形推定部、以及速度范围设定部。而且,发送接收部向包含被检体内的运动流体的诊断部位发送接收具有与表示可测定的速度的速度范围对应的重复频率的超声波。多普勒信号处理部根据通过上述超声波的发送接收而得到的信号,生成表示上述运动流体的速度的多普勒频谱图像。存储部预先存储有基于使基准的血流速度波形与ECG波形相关联的模型的模型化后的值。实测波形计算部根据特定的患者的上述频谱图像来求出实测的血流速度波形。ECG波形取得部取得与上述实测的血流速度波形对应的定时的ECG波形。波形推定部根据上述实测的血流速度波形、上述ECG波形、以及上述模型化后的值来推定上述患者的血流速度波形。速度范围设定部对上述实测的血流速度波形和上述推定的血流速度波形进行比较,根据除去了其差超过阈值的部分后的上述实测的血流速度波形来求出上述速度范围。
根据该第2方式,可以自动地求出与除去了阀信号的影响的血流速度波形对应的最佳的速度范围。由此,可以生成对于操作者而言易于观察的多普勒频谱图像。
本发明的第3方式的超声波图像取得装置具有发送接收部、多普勒信号处理部、存储部、以及异常判定部。而且,发送接收部向包含被检体内的运动流体的诊断部位发送接收具有基于表示可测定的速度的速度范围的重复频率的超声波。多普勒信号处理部根据通过上述超声波的发送接收而得到的信号,生成表示上述运动流体的速度的多普勒频谱图像。存储部预先存储有根据个别的ECG波形以及与该ECG波形对应的基准的血流速度波形的统计数据而模型化成一般的基准的血流速度波形的模型后的值。异常判定部根据上述个别的理想的血流速度波形的统计数据,来生成作为包含可视为正常状态的血流速度波形的区域的状态空间,当来自该状态空间的上述实测的血流信息的马哈拉诺比斯距离超过预定的阈值时,判定为血流速度异常。
根据该第3方式,可以自动地判定血流速度的异常。由此,可以减少血流速度的异常的遗漏,并且,可以早期地发现异常。
附图说明
图1是示出由超声波图像取得装置取得的断层像和多普勒频谱图像的图。
图2是第1实施方式的超声波图像取得装置的框图。
图3是用于说明第1实施方式的超声波图像取得装置的一系列动作的流程图。
图4A是用于使用本发明的超声波图像取得装置来求出血流速度波形的曲线图。
图4B是在以往的超声波图像取得装置中以手动操作设定了血流速度波形的速度范围时的曲线图。
图4C是示出实施例中的血流速度波形的理想的速度范围的曲线图。
图4D是使用第1实施方式的超声波图像取得装置来设定了血流速度波形的速度范围时的曲线图。
图5是用于说明系统识别中的参数模型的图。
图6是第2实施方式的超声波图像取得装置的框图。
具体实施方式
(第1实施方式)
以下对本发明的第1实施方式的超声波图像取得装置进行说明。图2是示出本发明的第1实施方式的超声波图像取得装置的概略结构的框图。
第1实施方式的超声波图像取得装置1可以按照显示超声波断层像(B模式断层像)的B模式、显示血流信息的多普勒模式(脉冲多普勒(PW)或连续波多普勒(CW))、二维地显示血流信息的CFM(Color Flow Mapping,彩色血流成像)模式等已知的模式来动作。
在超声波探测器2中,使用多个超声波振子沿着预定方向(扫描方向)排列成1列的一维超声波探测器、矩阵状地配置超声波振子的二维超声波探测器。
发送接收部3构成为具备:发送部(未图示),向超声波探测器2供给电信号而发生超声波;以及接收部(未图示),接收来自超声波探测器2的信号。
发送接收部3的发送部具备未图示的时钟发生电路、发送延迟电路、以及脉冲发生器电路。时钟发生电路是发生决定超声波信号的发送定时、发送频率的时钟信号的电路。发送延迟电路是在超声波的发送时施加延迟来实施发送聚焦(focus)的电路。脉冲发生器电路内置有与各超声波振子对应的各个路径数量的脉冲发生器,在被施加了延迟的发送定时发生驱动脉冲,供给到超声波探测器2的各超声波振子。
另外,发送接收部3的接收部具备未图示的前置放大器电路、A/D转换电路、以及接收延迟·加法电路。前置放大器电路针对每个接收通道对从超声波探测器2的各超声波振子输出的回波信号进行放大。A/D转换电路对放大后的回波信号进行A/D转换。接收延迟·加法电路向A/D转换后的回波信号提供决定接收指向性所必要的延迟时间并相加。通过该相加,来自与接收指向性对应的方向的反射成分被强调。另外,将由该发送接收部3进行加法处理后的信号称为“RF信号”。从发送接收部3输出的RF信号输入到B模式处理部4或多普勒处理部5。
进而,发送接收部3从速度范围设定部94接收到重复频率(PRF)后,按照重复频率使超声波探测器2发送接收超声波。
B模式处理部4进行回波的振幅信息的影像化,根据回波信号生成B模式超声波光栅数据。具体而言,B模式处理电路对RF信号进行带通滤波处理,之后,对输出信号的包络线进行检波,并对检波后的数据实施基于对数变换的压缩处理。由B模式处理部4生成的B模式超声波光栅数据输出到DSC 6。
多普勒处理部5构成为具备正交相位检波部51、范围门(RangeGate:RG)处理部52、壁滤波器53以及FFT运算部54。
正交相位检波部51将来自基准发送器的基准信号及其具有90度相位差的基准信号、与从发送接收部3输出的RF信号混合。范围门(RG)处理部52除去混合信号中的高频分量而取得仅由多普勒偏移频率分量构成的多普勒信号,并从其中提取出被检体内的期望深度的多普勒信号。壁滤波器53从来自范围门处理部52的、由范围门指定的被检体内的预定位置的多普勒信号中,除去动作比较慢的血管壁、心脏壁等不需要的低频多普勒信号,提取出应检测的血流的多普勒信号。FFT运算部54对由壁滤波器53提取出的多普勒信号进行频率解析而取得作为其解析结果的多普勒频谱信号,并输出到DSC(Digital Scan Converter,数字扫描转换器)6。由此,在显示部8上,例如与B模式断层像一起显示出多普勒频谱图像。
另外,FFT运算部54从速度范围设定部94接收到速度偏差(BLS)后,使FFT处理的读出地址变化该移动量来调整速度的偏差。
为了取得以正交坐标系表示的图像,DSC 6将超声波光栅数据转换成以正交坐标系表示的图像数据(扫描转换处理)。从DSC 6向显示部8输出图像数据,在显示部8上显示基于该图像数据的图像。例如,DSC 6根据B模式超声波光栅数据生成作为二维信息的断层像数据,并向显示部8输出该断层像数据。显示部8显示基于该断层像数据的断层像。
自动范围/自动BLS处理部9接收从FFT运算部54输出的多普勒频谱信号,计算出最适合于该多普勒频谱信号的速度范围。自动范围/自动BLS处理部9由实测波形计算部91、ECG波形取得部92、波形推定部93、插值部94以及速度范围设定部95构成。以下,对自动范围/自动BLS处理部9进行详细说明。
实测波形计算部91通过对从FFT运算部54输出的多普勒频谱信号的最大速度Vp进行检测并在时间方向上连线(例如,可以是检波,也可以是伪滤波)(以下称为“描绘”),而检测出Vp的实测的血流速度波形。由此,最大速度Vp的实测的血流速度波形成为描绘了多普勒频谱图像的最大速度Vp后得到的波形。然后,实测波形计算部91根据所检测出的实测的血流速度波形计算出表示实测的血流速度波形的函数Vp(t)(以下简称为“实测的血流速度波形Vp(t)”)。进而,实测波形计算部91向波形推定部93以及插值部94发送实测的血流速度波形Vp(t)。
ECG波形取得部91从心电计10接收信号,使其与实测波形计算部91求出的实测的血流速度波形同步,进而在时间方向上描绘最大振幅来生成ECG波形。进而,ECG波形取得部91根据所生成的ECG波形来生成表示该ECG波形的函数u(t)(以下简称为“ECG波形u(t)”)。然后,ECG波形取得部91向波形推定部93发送所生成的ECG波形u(t)。
存储部7根据事先针对每个年龄以及疾病收集的个别患者的ECG波形以及与该ECG波形对应的没有阀信号影响的个别患者的理想的血流速度波形的统计数据,使用ARX模型(Auto-Regressiveexogeneous model,外生自回归模型)来进行系统识别,模型化成一般的理想血流速度波形的模型,并预先存储该模型的系数序列。该没有阀信号影响的理想的血流速度波形相当于本发明中的“基准的血流速度波形”。此处,ARX模型是指使用于系统识别的线性时变的参数模型,是将当前的输出y(t)与有限个过去的输出数据y(t-k)和输入数据u(t-k)关联起来的模型(例如,参照“基于MATLAB的系统识别”足立修一著东京电机大学出版社)。
此处,说明上述的系数序列的求法。首先对年龄以及所具有的疾病进行分组,针对每个组收集几百个病历的ECG(Electrocardiogram,心电图)波形以及没有阀信号影响的理想的血流速度波形。此处,作为理想的血流速度波形的求法,有如下的方法:通过在循环器中的不存在相互阀的部分进行血流速度的测定而求出没有阀信号影响的血流速度波形,或者由经验丰富的医师根据经验从实测的血流速度波形手动地除去阀信号的影响而求出没有阀信号影响的血流速度波形。
接下来,根据所收集到的ECG波形以及理想的血流速度波形来进行基于ARX模型的系统识别。设表示ECG波形的函数为u(t)、表示一般的理想的血流速度波形的函数为yi(t),则
A(q)*yi(t)=B(q)*u(t-nk)+e(t)
e(t):残差,即期待值与实测值之差
nk:与和ECG波形对应的一般的理想的血流速度波形的时间延迟
A(q)=1+a1q-1+...+anaq-na
B(q)=b1+b2q-1+...+bnbq-nb+1
(A(q)、B(q)为不可约的移位算子q的多项式)
na、nb、nk:整数自变量。
因此,在存储部7中,针对每个被分类的组,预先存储一般的理想的血流速度波形的模型中的系数,即ai=(a1、a2、...、ana)、bj=(b1、b2、...、bnb)。
波形推定部93根据从实测波形计算部91接收到的实测的血流速度波形Vp(t)、从ECG波形取得部92接收到的ECG波形u(t)、以及存储在存储部7中的系数ai=(a1、a2、...、ana)、bj=(b1、b2、...、bnb),推定当前正在计测血流速度的患者的理想的血流速度波形。在此,对于在该推定中使用的系数,判断患者属于上述的哪个组来使用与该组对应的系数。以下,将该推定出的波形称为推定波形ye(t)。
该推定波形ye(t)可以表示成
ye(t)=-∑{ai*Vp(t-i)}+∑{bj*u(t-j)}。
波形推定部95向插值部94发送推定波形ye(t)。
插值部94存储有实测的血流速度波形Vp(t)与推定波形ye(t)之差的阈值。此处,该阈值是用于在实测的血流速度波形Vp(t)与推定波形ye(t)之差超过了该值的情况下,判断为在该时间实测的血流速度波形Vp(t)受到了阀信号的影响的值。这样,该阈值的目的在于除去非常高速的阀信号的影响,相反,如果较强地要求与推定波形的近似,则会除去大部分的值。因此,优选按照必须留下多大范围来作为希望除去影响的阀信号的影响范围、以及除去了阀信号的影响后的实测的血流速度波形Vp(t)的范围的要求,来设定该阈值。然后,插值部94求出所接收到的当前正在计测血流速度的患者的实测的血流速度波形Vp(t)与推定波形ye(t)之差,来判定超过了阈值的部分。
在本实施方式中,阈值被提供给正、负这两方。但是,该阈值的设定方法也可以是其他方法。例如,存在阀的影响的部分与推定波形相比值更大,所以在从实测的血流速度波形中减去推定波形的情况下,可以仅对正侧设定阈值的值。另外,在从推定波形中减去实测的血流速度波形的情况下,可以仅对负侧设定阈值的值。进而,也可以取实测的血流速度波形与推定波形之差的绝对值,针对该绝对值设定阈值。
插值部94从当前正在计测血流速度的患者的实测的血流速度波形Vp(t)中除去超过了阈值的部分。
进而,插值部94对实测的血流速度波形Vp(t)的所除去的部分适用推定波形ye(t)。然后,插值部94对该推定波形ye(t)和实测的血流速度波形Vp(t)的差距进行插值。该插值方法没有特别限定,例如可以是直线插值、样条插值、或根据实测的血流速度波形的峰值进行线性预测来进行插值的方法等任意方法。另外,在本实施方式中,在对所除去的部分适用了推定波形ye(t)后进行插值,但也可以不适用推定波形ye(t)而直接对实测的血流速度波形Vp(t)的所除去的部分进行插值。该插值方法也没有特别限定,例如可以是直线插值、样条插值、或根据实测的血流速度波形的峰值进行线性预测来进行插值的方法等。此处,直线插值由于处理简单,所以对插值部94的负荷较少即可。而在样条插值、根据实测的血流速度波形的峰值进行线性预测来进行插值的方法的情况下,处理烦杂,但可以使用平滑的曲线来进行插值,从而可以生成更接近于实际波形的波形。进而,样条插值虽然对要插值的区间短的情况有效,但根据实测的血流速度波形的峰值进行线性预测来进行插值的方法可以在要插值的区间长的情况下进行插值。以下,将该插值部94所生成的血流速度波形称为“范围设定用波形”。
插值部94向速度范围设定部95发送范围设定用波形。
速度范围设定部95进行统计计算处理来计算出速度范围的上限以及下限。此处,在统计计算处理中存在正则分布模型和平滑化后阈值处理模型等,其中正则分布模型从范围设定用波形生成直方图,求出速度分布,根据进行了基于该速度分布的加权后的范围设定用波形的分布计算出平均和方差,将平均±系数*σ设为速度范围的上限以及下限的推定值,而平滑化后阈值处理模型根据加权后的范围设定用波形的分布,将相当于峰值的系数%的值设为速度范围的上限以及下限。该速度范围的计算方法详细记载于日本特开2005-185731号公报。
然后,速度范围设定部95求出与所设定的速度范围对应的重复频率(PRF)。
另外,速度范围设定部95使用范围设定用波形,根据现行的基准位置(基线)上侧(正侧)的最大值求出最大速度范围,根据下侧(负侧)的最大值求出最小速度范围,进行最大速度范围和最小速度范围的比较,求出基线(基准位置)的移动量。例如,使最大速度范围和最小速度范围的平均值置于画面中央,通过求出其平均值与基线=0的距离,来求出移位量。进行作为基准位置(基线)的移动量的速度偏差(BLS)的计算。该速度偏差(BLS)的计算方法详细记载于日本特开2005-185731号公报。
如上所述决定了重复频率(PRF)以及速度偏差(BLS)后,速度范围设定部95向发送接收部3输出重复频率(PRF)。同时,速度范围设定部95向多普勒信号处理部5的FFT运算部54输出速度偏差(BLS)。
另外,在本实施方式中,自动范围/自动BLS处理部9可以由硬件构成,也可以构成为软件。例如,可以由CPU构成自动范围/自动BLS处理部9,从存储区域(未图示)读入程序并执行该程序,由此来执行实测波形计算部91、ECG波形取得部92、波形推定部93、插值部94以及速度范围设定部95的功能。
此处,虽然在本实施方式中一同进行了速度偏差的调整以及速度范围的调整,但本发明的超声波图像取得装置也可以仅通过速度范围的调整来动作。在该情况下,偏差不移动而总是以0Hz为基准,但进行速度范围的调整而得到减轻了折叠的多普勒频谱图像。
(动作)
接下来,参照图3来对本发明的第1实施方式的超声波图像取得装置的一系列动作进行说明。图3是示出本发明的第1实施方式的超声波图像取得装置的一系列动作的流程图。
(步骤S001)
首先,对被检体发送超声波,根据来自被检体的反射波来生成B模式断层像和多普勒频谱图像。由FFT运算部54生成的多普勒波形数据从FFT运算部54经由DSC 6输出到显示部8,与B模式断层像一起显示在显示部8上。进而,多普勒波形数据从FFT运算部54输出到自动描绘部7。
(步骤S002)
在步骤S001取得了多普勒波形后,实测波形计算部91在时间方向上描绘多普勒频谱图像的边缘部(最大速度Vp),对最大速度Vp的血流速度波形进行检测。进而,实测波形计算部91求出作为表示最大速度Vp的血流速度波形的函数的实测的血流速度波形Vp(t)。然后,实测波形计算部91向波形推定部93以及插值部94输出该实测的血流速度波形Vp(t)。
(步骤S003)
在步骤S003中,ECG波形取得部92描绘出从心电计10接收到的信号的最大振幅,求出作为表示该描绘出的波形的函数的ECG波形u(t)。然后,ECG波形取得部92向波形推定部93输出该ECG波形u(t)。
(步骤S004)
在步骤S004中,波形推定部93根据实测的血流速度波形Vp(t)、ECG波形u(t)、以及预先存储在存储部7中的一般的理想的血流速度波形的模型的系数来推定当前正在进行血流速度的测定的患者的理想的血流速度波形(推定波形ye(t))。然后,波形推定部93向插值部94输出推定波形ye(t)。
(步骤S005)
插值部94求出实测的血流速度波形Vp(t)以及推定波形ye(t)之差,从实测的血流速度波形Vp(t)中除去该差超过了所存储阈值的部分,进而通过进行该除去的部分的插值,从而求出范围设定用波形。
(步骤S006)
速度范围设定部95根据范围设定用波形进行统计值计算处理来求出重复频率(PRF)以及速度偏差(BLS)。
(步骤S007)
如上所述决定了重复频率(PRF)和速度偏差(BLS)后,速度范围设定部95向发送接收部3输出该重复频率(PRF)。同时,速度范围设定部95向多普勒信号处理部5内的FFT运算部54输出速度偏差(BLS)。发送接收部3按照由速度范围设定部95求出的重复频率(PRF)来使超声波探测器2发送接收超声波。另外,FFT运算部54按照由速度范围设定部95求出的速度偏差(BLS),使FFT处理的读出地址变化该移位量,来调整速度的偏差。由此,速度范围和速度偏差(BLS)被更新。
如上所述,本实施方式的超声波图像取得装置1根据通过扫描取得的多普勒频谱图像来自动地求出除去了阀信号的影响的血流速度波形的速度范围以及速度偏差。然后,可以使用该速度范围和速度偏差,追踪血流的状态偏移来改变多普勒速度范围等。
(实施例)
以下,参照图4A~D,利用通过使用本实施方式的超声波图像取得装置来除去了阀信号的影响的血流速度波形的实施例、和不使用本实施方式的超声波图像取得装置而使用以往方法的血流速度波形的比较例,来对血流速度波形的显示进行说明。图4A是用于使用本发明的超声波图像取得装置来求出血流速度波形的曲线图。图4B是在以往的超声波图像取得装置中以手动操作设定了血流速度波形的速度范围时的曲线图。图4C是示出本实施例中的血流速度波形的理想的速度范围的曲线图。图4D是使用本实施方式的超声波图像取得装置设定了血流速度波形的速度范围时的曲线图。在图4的各曲线中,纵轴表示速度范围(kHz),横轴表示时间(sec)。
图4A所示的曲线501是示出除去阀的影响之前的实测的血流速度波形Vp(t)的曲线,曲线502是示出推定波形ye(t)的曲线,曲线503是示出与曲线501的血流速度波形对应的ECG波形ECG(t)的曲线,曲线504是使用本实施方式的超声波图像取得装置除去了阀的影响后的血流速度波形的曲线。
在本实施例中,首先,通过对所实测的多普勒信号进行多普勒自动描绘而生成曲线501示出的实测的血流速度波形Vp(t)。然后,根据曲线503示出的ECG波形ECG(t)和统计求出的理想的血流速度波形,生成曲线502示出的推定波形ye(t)。然后,通过从实测的血流速度波形Vp(t)除去与该推定波形ye(t)的差超过了阈值的部分并对该部分进行插值,来生成曲线504示出的除去了阀的影响的血流速度波形。
此处,在不从实测波形Vp(t)除去阀的影响而显示血流速度波形的情况下,速度范围需要是5.6(kHz)。而在利用以往技术并通过手动操作从实测波形Vp(t)除去了阀的影响的情况下,可以如图4B所示抑制成3.9(kHz)的速度范围。但是,如果根据考虑要极力抑制阀的影响的情况而统计学地计算出的理想的波形来求出速度范围,则在本实施例中的血流速度波形中,理想的是将速度范围抑制成3.2(kHz)。因此,使用本实施方式的超声波图像取得装置,在使用了曲线504示出的除去了阀的影响的血流速度波形的情况下,可以如图4D所示将速度范围抑制成3.2(kHz)。
如上所述,难以通过手动操作接近于理想的速度范围,但可以通过利用本实施方式的超声波图像取得装置除去阀的影响,来使该速度范围大致近似于理想的速度范围。
此处,在本实施方式中在进行系统识别时使用ARX模型,但只要是系统识别中的参数模型(例如参照“用于控制的上级系统识别”足立修一著东京电机大学),则也可以使用其他数学模型。作为该参数模型,除了ARX模型以外,例如还有FIR(Finite Impulse Response,有限脉冲响应)模型、ARMAX(Auto Regressive Moving AverageeXogenious,自回归移动平均外生)模型、OE(Output Error,输出误差)模型、BJ(Box and Jenkins,博克斯和詹金斯)模型等。系统识别中的参数模型如图5所示。图5是用于说明系统识别中的参数模型的图。此处在图5中,A(z)*y(k)={B(z)/F(z)}*u(k)+{C(z)/D(z)}*w(k)
A(z)=1+a1*z-1+...+an*z-n
B(z)=b1*z-1+b2*z-2+...+bm*z-m
C(z)=c1*z-1+c2*z-2+...+cp*z-p
D(z)=1+d1*z-1+d2*z-2+...+dq*z-q
F(z)=1+f1*z-1+f2*z-2+...+fr*z-r
于是,该系统识别中的参数模型表示为:
e(k)+d1*e(k-1)+...+dq*e(k-q)=c1*w(k-1)+c2*w(k-2)+cp*w(k-p)
x(k)+f1*x(k-1)+...+fq*x(k-q)=b1*u(k-1)+b2*u(k-2)+bm*w(k-m)
y(k)+a1*y(k-1)+...+aq*y(k-q)=e(k)+x(k)
另外,在本实施方式中在进行了插值后进行速度范围的调整,但也可以仅生成进行插值并除去了阀的影响的血流速度波形。利用该结构可以取得除去了阀的影响的血流速度波形,从而可以易于进行血流的流量的计测等而使计测精度提高。
进而,也可以使用不进行插值而仅除去了阀的影响的血流速度波形来进行速度范围的调整。在该情况下,速度范围设定部95存储阈值。而且,速度范围设定部95进行推定波形和实测的血流速度波形的比较来求出超过阈值的部分。利用该结构,速度范围调整的精度虽然降低几%,但由于计算简单,所以可以使速度范围的计算速度提高。
(第2实施方式)
接下来对本发明的第2实施方式的超声波图像取得装置进行说明。第2实施方式的超声波图像取得装置1A构成为在第1实施方式的超声波图像取得装置1中还具备异常判定部11。图6是示出第2实施方式的超声波图像取得装置的功能的框图。以下对超声波图像取得装置1A中的异常的判定进行说明。
异常判定部11根据作为统计数据而存储在存储部7中的每个年龄以及疾病的、没有阀信号的影响的个别的理想的血流速度波形来生成状态空间。此处,状态空间是指,可视为正常状态的血流速度波形包含在其中的空间。
异常判定部11从实测波形计算部91取得实测的血流速度波形Vp(t)。
异常判定部11根据具有该血流速度波形Vp(t)的患者的年龄以及疾病来选择所对应的状态空间,求出该状态空间与血流速度波形Vp(t)的马哈拉诺比斯(Mahalanobis)距离。此处,马哈拉诺比斯距离是使用了变量间的相关的指标,是表示离开作为对象的基准空间的距离的距离尺度。
异常判定部11在所计算出的马哈拉诺比斯距离超过了预先存储的针对马哈拉诺比斯距离的阈值的情况下判断为血流速度异常。然后,异常判定部11例如使用显示部8等向操作者通知血流速度异常的信息。
由此,在发生了超过阈值那样的血流速度异常的情况下。自动地向操作者通知,所以可以减少血流速度的异常的遗漏,从而可以向患者提供更迅速且正确的诊断。
Claims (14)
1.一种超声波图像取得装置,其特征在于,具有:
发送接收部,向包含被检体内的运动流体的诊断部位发送接收具有与表示可测定的速度的速度范围对应的重复频率的超声波;
多普勒信号处理部,根据通过上述超声波的发送接收而得到的信号,生成表示上述运动流体的速度的多普勒频谱图像;
存储部,预先存储有基于使基准的血流速度波形与ECG波形相关联的模型的模型化后的值;
实测波形计算部,根据特定的患者的上述频谱图像来求出实测的血流速度波形;
ECG波形取得部,取得与上述实测的血流速度波形对应的定时的ECG波形;
波形推定部,根据上述实测的血流速度波形、上述ECG波形、以及上述模型化后的值来推定上述患者的血流速度波形;以及
插值部,对上述实测的血流速度波形和上述推定的血流速度波形进行比较,从上述实测的血流速度波形中除去其差超过阈值的部分,对上述实测的血流速度波形的该除去的部分进行插值。
2.根据权利要求1所述的超声波图像取得装置,其特征在于,上述基准的血流速度波形是没有阀信号的影响的血流速度波形。
3.根据权利要求1所述的超声波图像取得装置,其特征在于,上述基准的血流速度波形的模型是使用了系统识别中的参数模型的模型,上述模型化后的值是进行基于上述系统识别中的参数模型的系统识别而求出的系数。
4.根据权利要求1所述的超声波图像取得装置,其特征在于,上述插值部使用直线插值、样条插值、或利用根据上述实测的血流速度波形中的峰值的线性预测的插值中的任意一个来进行上述插值。
5.根据权利要求1所述的超声波图像取得装置,其特征在于,上述插值部利用上述推定的血流速度波形来置换上述推定的血流速度波形与上述实测的血流速度波形之差超过阈值的上述实测的血流速度波形的部分。
6.一种超声波图像取得装置,其特征在于,具有:
发送接收部,向包含被检体内的运动流体的诊断部位发送接收具有与表示可测定的速度的速度范围对应的重复频率的超声波;
多普勒信号处理部,根据通过上述超声波的发送接收而得到的信号,生成表示上述运动流体的速度的多普勒频谱图像;
存储部,预先存储有基于使基准的血流速度波形与ECG波形相关联的模型的模型化后的值;
实测波形计算部,根据特定的患者的上述频谱图像来求出实测的血流速度波形;
ECG波形取得部,取得与上述实测的血流速度波形对应的定时的ECG波形;
波形推定部,根据上述实测的血流速度波形、上述ECG波形、以及上述模型化后的值来推定上述患者的血流速度波形;以及
速度范围设定部,对上述实测的血流速度波形和上述推定的血流速度波形进行比较,根据除去了其差超过阈值的部分后的上述实测的血流速度波形来求出上述速度范围。
7.根据权利要求6所述的超声波图像取得装置,其特征在于,上述基准的血流速度波形是没有阀信号的影响的血流速度波形。
8.根据权利要求6所述的超声波图像取得装置,其特征在于,上述基准的血流速度波形的模型是使用了系统识别中的参数模型的模型,上述模型化后的值是进行基于上述系统识别中的参数模型的系统识别而求出的系数。
9.根据权利要求6所述的超声波图像取得装置,其特征在于,上述速度范围设定部针对每个预定观测时间,根据对除去了上述超过阈值的部分后的上述实测的血流速度的最大流速以及平均流速中的至少一方的分布的统计值,求出速度范围。
10.根据权利要求6所述的超声波图像取得装置,其特征在于,上述速度范围设定部还具备速度偏差调整部,该速度偏差调整部根据除去了上述超过阈值的部分后的上述实测的血流速度波形的最大流速的正侧最大值以及负侧最大值,求出多普勒频谱图像显示中的速度方向的基准位置的移动量,来调整速度偏差。
11.根据权利要求6所述的超声波图像取得装置,其特征在于,还具备插值部,该插值部从上述实测的血流速度波形中除去上述推定的血流速度波形与上述实测的血流速度波形之差超过阈值的部分,并对上述实测的血流速度波形的该除去的部分进行插值。
12.根据权利要求11所述的超声波图像取得装置,其特征在于,上述插值部使用直线插值、样条插值、或利用根据上述实测的血流速度波形中的峰值的线性预测的插值中的任意一个来进行上述插值。
13.根据权利要求11所述的超声波图像取得装置,其特征在于,上述插值部利用上述推定的血流速度波形来置换上述推定的血流速度波形与上述实测的血流速度波形之差超过阈值的上述实测的血流速度波形的部分。
14.一种超声波图像取得装置,其特征在于,具有:
发送接收部,向包含被检体内的运动流体的诊断部位发送接收具有基于表示可测定的速度的速度范围的重复频率的超声波;
多普勒信号处理部,根据通过上述超声波的发送接收而得到的信号,生成表示上述运动流体的速度的多普勒频谱图像;
存储部,预先存储有根据个别的ECG波形以及与该ECG波形对应的基准的血流速度波形的统计数据而模型化成一般的基准的血流速度波形的模型后的值;以及
异常判定部,根据上述个别的理想的血流速度波形的统计数据,来生成作为包含可视为正常状态的血流速度波形的区域的状态空间,当来自该状态空间的上述实测的血流信息的马哈拉诺比斯距离超过预定的阈值时,判定为血流速度异常。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007185882A JP4960162B2 (ja) | 2007-07-17 | 2007-07-17 | 超音波診断装置 |
JP2007185882 | 2007-07-17 | ||
JP2007-185882 | 2007-07-17 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101347343A true CN101347343A (zh) | 2009-01-21 |
CN101347343B CN101347343B (zh) | 2010-12-22 |
Family
ID=40265414
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101286852A Expired - Fee Related CN101347343B (zh) | 2007-07-17 | 2008-06-23 | 超声波图像取得装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8523776B2 (zh) |
JP (1) | JP4960162B2 (zh) |
CN (1) | CN101347343B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101933816A (zh) * | 2009-06-25 | 2011-01-05 | 株式会社东芝 | 三维超声波诊断装置 |
CN102100567A (zh) * | 2009-12-21 | 2011-06-22 | 株式会社东芝 | 彩色多普勒超声波诊断装置 |
CN102551791A (zh) * | 2010-12-17 | 2012-07-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像方法和装置 |
CN103126724A (zh) * | 2011-11-30 | 2013-06-05 | 株式会社东芝 | 超声波诊断装置以及图像处理方法 |
CN104220005A (zh) * | 2012-04-18 | 2014-12-17 | 日立阿洛卡医疗株式会社 | 超声波摄像装置以及超声波摄像方法 |
CN105212963A (zh) * | 2015-08-28 | 2016-01-06 | 飞依诺科技(苏州)有限公司 | 实现超声彩色多普勒血流成像的壁滤波方法及系统 |
CN105380680A (zh) * | 2014-09-01 | 2016-03-09 | 三星麦迪森株式会社 | 超声诊断设备和操作该超声诊断设备的方法 |
CN105476661A (zh) * | 2016-01-13 | 2016-04-13 | 北京悦琦创通科技有限公司 | 调整多普勒频谱图的方法和设备 |
CN107123112A (zh) * | 2017-01-23 | 2017-09-01 | 上海联影医疗科技有限公司 | 血流状态分析系统及方法 |
CN107920802A (zh) * | 2015-08-27 | 2018-04-17 | 皇家飞利浦有限公司 | 利用自适应采样窗口尺寸的谱多普勒处理 |
CN108169331A (zh) * | 2017-12-04 | 2018-06-15 | 北京星航机电装备有限公司 | 薄板栅格翼结构焊缝相控阵超声检测装置及检测方法 |
CN116849702A (zh) * | 2023-06-01 | 2023-10-10 | 南方科技大学医院(深圳市南山区西丽人民医院) | 一种基于三维超声心动图的肾脏健康状况的评价方法及系统 |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5295684B2 (ja) * | 2008-08-21 | 2013-09-18 | 株式会社東芝 | 超音波診断装置及び診断パラメータ自動計測方法 |
US8792876B1 (en) * | 2008-12-12 | 2014-07-29 | Cisco Technology, Inc. | System and method for provisioning flows in a WiMAX network environment |
JP5558727B2 (ja) * | 2009-02-27 | 2014-07-23 | 株式会社東芝 | 超音波診断装置および超音波診断装置のデータ処理プログラム |
JP5481108B2 (ja) | 2009-06-26 | 2014-04-23 | 株式会社東芝 | 超音波診断装置及び自動診断支援装置 |
KR100977367B1 (ko) * | 2009-10-28 | 2010-08-20 | (주)메디슨 | 3차원 pw 스펙트럼 초음파 진단 장치 및 3차원 pw 스펙트럼 데이터 생성 방법 |
CN102053051A (zh) * | 2009-10-30 | 2011-05-11 | 西门子公司 | 一种体液分析系统和用于体液分析的图像处理设备、方法 |
WO2011135482A1 (en) * | 2010-04-28 | 2011-11-03 | Koninklijke Philips Electronics N.V. | Property determining apparatus for determining a property of an object |
GB2488092B (en) * | 2010-11-03 | 2014-10-29 | Kittiwake Developments Ltd | A sensor based means of monitoring the mechanical condition of rotating machinery that operates intermittently |
WO2013057999A1 (ja) * | 2011-10-20 | 2013-04-25 | 日立アロカメディカル株式会社 | 超音波撮像装置及び方法 |
US9890870B2 (en) * | 2011-12-29 | 2018-02-13 | Schneider Electric Buildings Llc | Valve flow control optimization via customization of an intelligent actuator |
KR101462174B1 (ko) | 2013-01-29 | 2014-11-20 | 삼성전자주식회사 | 영상 표시 장치 및 영상 표시 방법 |
JP6391912B2 (ja) * | 2013-02-26 | 2018-09-19 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置 |
EP4220221A1 (en) | 2013-03-15 | 2023-08-02 | BFLY Operations, Inc. | Monolithic ultrasonic imaging devices, systems and methods |
CA2919183A1 (en) | 2013-07-23 | 2015-01-29 | Butterfly Network, Inc. | Interconnectable ultrasound transducer probes and related methods and apparatus |
JP6552599B2 (ja) | 2014-04-18 | 2019-07-31 | バタフライ ネットワーク,インコーポレイテッド | 単一基板超音波撮像装置の構造、関連装置、及び方法 |
WO2016176855A1 (zh) * | 2015-05-07 | 2016-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | 血流成像方法及系统 |
CN106419890B (zh) * | 2016-11-14 | 2024-04-30 | 佛山科学技术学院 | 一种基于时空调制的血流速度测量装置及方法 |
CN108814649B (zh) * | 2018-07-05 | 2024-04-05 | 河南省计量测试科学研究院 | 多普勒超声诊断仪血流波形模体及校准方法 |
US11365992B2 (en) * | 2019-12-10 | 2022-06-21 | Schneider Electric Systems Usa, Inc. | Systems and methods for validating a vortex fluid flow signal |
CN111860636A (zh) * | 2020-07-16 | 2020-10-30 | 无锡祥生医疗科技股份有限公司 | 测量信息提示方法和超声训练方法 |
CN113066083B (zh) * | 2021-04-25 | 2022-08-02 | 青岛海信医疗设备股份有限公司 | 流体的多普勒参数确定方法和电子设备 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0420085B1 (de) * | 1989-09-27 | 1998-07-01 | Pyrotec Limited | Messung der Herzfunktion |
JP3109090B2 (ja) * | 1990-08-31 | 2000-11-13 | 日本電気株式会社 | 半導体記憶装置の製造方法 |
US5634465A (en) * | 1995-06-09 | 1997-06-03 | Advanced Technology Laboratories, Inc. | Continuous display of cardiac blood flow information |
US5735281A (en) * | 1996-08-09 | 1998-04-07 | Hewlett-Packard Company | Method of enhancing and prolonging the effect of ultrasound contrast agents |
JP3892538B2 (ja) * | 1997-07-18 | 2007-03-14 | 株式会社東芝 | 超音波ドプラ診断装置 |
US6024701A (en) * | 1998-08-27 | 2000-02-15 | T.A.O. Medical Technologies Ltd. | Method of and system for estimating placenta and fetus well being using system identification techniques |
JP4240633B2 (ja) * | 1999-02-10 | 2009-03-18 | 株式会社東芝 | 超音波ドプラ診断装置 |
US6340346B1 (en) * | 1999-11-26 | 2002-01-22 | T.A.O. Medical Technologies Ltd. | Method and system for system identification of physiological systems |
US6692438B2 (en) * | 2001-12-18 | 2004-02-17 | Koninklijke Philips Electronics Nv | Ultrasonic imaging system and method for displaying tissue perfusion and other parameters varying with time |
US6733454B1 (en) * | 2003-02-26 | 2004-05-11 | Siemens Medical Solutions Usa, Inc. | Automatic optimization methods and systems for doppler ultrasound imaging |
JP4245428B2 (ja) * | 2003-07-29 | 2009-03-25 | 東芝医用システムエンジニアリング株式会社 | 超音波ドプラ診断装置及び画像データ生成方法 |
JP4393137B2 (ja) * | 2003-08-25 | 2010-01-06 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
WO2005039418A1 (en) * | 2003-10-23 | 2005-05-06 | Koninklijke Philips Electronics, N.V. | Ultrasound imaging method and apparatus |
JP4476619B2 (ja) | 2003-12-26 | 2010-06-09 | 株式会社東芝 | ドプラ超音波診断装置 |
US8303507B2 (en) * | 2004-09-07 | 2012-11-06 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnostic apparatus and measuring method of diagnostic parameter |
US7513872B2 (en) * | 2004-10-18 | 2009-04-07 | Kabushiki Kaisha Toshiba | Ultrasonic doppler measuring apparatus and control method therefor |
JP5085170B2 (ja) * | 2007-03-28 | 2012-11-28 | 株式会社東芝 | 超音波診断装置 |
-
2007
- 2007-07-17 JP JP2007185882A patent/JP4960162B2/ja not_active Expired - Fee Related
-
2008
- 2008-06-17 US US12/140,607 patent/US8523776B2/en active Active
- 2008-06-23 CN CN2008101286852A patent/CN101347343B/zh not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101933816B (zh) * | 2009-06-25 | 2013-12-04 | 株式会社东芝 | 三维超声波诊断装置 |
CN101933816A (zh) * | 2009-06-25 | 2011-01-05 | 株式会社东芝 | 三维超声波诊断装置 |
CN102100567A (zh) * | 2009-12-21 | 2011-06-22 | 株式会社东芝 | 彩色多普勒超声波诊断装置 |
US9907531B2 (en) | 2010-12-17 | 2018-03-06 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd | Methods and systems for ultrasonic imaging |
CN102551791B (zh) * | 2010-12-17 | 2016-04-27 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像方法和装置 |
CN102551791A (zh) * | 2010-12-17 | 2012-07-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像方法和装置 |
CN103126724A (zh) * | 2011-11-30 | 2013-06-05 | 株式会社东芝 | 超声波诊断装置以及图像处理方法 |
CN104220005A (zh) * | 2012-04-18 | 2014-12-17 | 日立阿洛卡医疗株式会社 | 超声波摄像装置以及超声波摄像方法 |
CN105380680A (zh) * | 2014-09-01 | 2016-03-09 | 三星麦迪森株式会社 | 超声诊断设备和操作该超声诊断设备的方法 |
CN107920802A (zh) * | 2015-08-27 | 2018-04-17 | 皇家飞利浦有限公司 | 利用自适应采样窗口尺寸的谱多普勒处理 |
CN105212963B (zh) * | 2015-08-28 | 2018-02-06 | 飞依诺科技(苏州)有限公司 | 实现超声彩色多普勒血流成像的壁滤波方法及系统 |
CN105212963A (zh) * | 2015-08-28 | 2016-01-06 | 飞依诺科技(苏州)有限公司 | 实现超声彩色多普勒血流成像的壁滤波方法及系统 |
CN105476661A (zh) * | 2016-01-13 | 2016-04-13 | 北京悦琦创通科技有限公司 | 调整多普勒频谱图的方法和设备 |
CN107123112A (zh) * | 2017-01-23 | 2017-09-01 | 上海联影医疗科技有限公司 | 血流状态分析系统及方法 |
CN107123112B (zh) * | 2017-01-23 | 2020-08-04 | 上海联影医疗科技有限公司 | 血流状态分析系统及方法 |
CN108169331A (zh) * | 2017-12-04 | 2018-06-15 | 北京星航机电装备有限公司 | 薄板栅格翼结构焊缝相控阵超声检测装置及检测方法 |
CN116849702A (zh) * | 2023-06-01 | 2023-10-10 | 南方科技大学医院(深圳市南山区西丽人民医院) | 一种基于三维超声心动图的肾脏健康状况的评价方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
JP4960162B2 (ja) | 2012-06-27 |
US8523776B2 (en) | 2013-09-03 |
JP2009022364A (ja) | 2009-02-05 |
CN101347343B (zh) | 2010-12-22 |
US20090024037A1 (en) | 2009-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101347343B (zh) | 超声波图像取得装置 | |
US11635514B2 (en) | Imaging methods and apparatuses for performing shear wave elastography imaging | |
JP4864547B2 (ja) | 超音波診断装置およびその制御処理プログラム | |
DE102012101395B4 (de) | Mehrstrahliger Spektral-Doppler in der medizinischen diagnostischen Ultraschallbildgebung | |
JP4676334B2 (ja) | 生体信号モニタ装置 | |
EP2022404B1 (en) | Ultrasonograph, medical image processing device, and medical image processing program | |
US8303507B2 (en) | Ultrasonic doppler diagnostic apparatus and measuring method of diagnostic parameter | |
JP3746115B2 (ja) | 超音波診断装置 | |
US8265358B2 (en) | Ultrasonic image processing apparatus and method for processing ultrasonic image | |
JP5100343B2 (ja) | 超音波診断装置、及び超音波診断装置の制御プログラム | |
US7175598B2 (en) | Ultrasound diagnosis apparatus that adjusts a time phase between a plurality of image series | |
US11471130B2 (en) | Method and ultrasound system for shear wave elasticity imaging | |
CN101516270B (zh) | 超声波诊断装置 | |
US8888703B2 (en) | Ultrasound imaging apparatus and method for processing image with motion tracking and error correction | |
EP1637081B1 (en) | Ultrasonograph | |
JP5231768B2 (ja) | 超音波診断装置および超音波診断装置のデータ処理プログラム | |
JP4763588B2 (ja) | 超音波診断装置 | |
US10624608B2 (en) | Ultrasonic diagnostic apparatus | |
DE60036056T2 (de) | Ultraschalldiagnosegerät | |
CN102038522A (zh) | 超声波诊断装置以及超声波数据取得方法 | |
US11039777B2 (en) | Ultrasonic diagnostic apparatus and control method | |
CN104665877A (zh) | 颈动脉血管局部脉搏波传播速度测量方法 | |
US20130158403A1 (en) | Method for Obtaining a Three-Dimensional Velocity Measurement of a Tissue | |
JP2010046284A (ja) | 超音波診断装置及び診断パラメータ自動計測方法 | |
CN106659470B (zh) | 超声波诊断装置 |
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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160712 Address after: Japan Tochigi Patentee after: TOSHIBA MEDICAL SYSTEMS Corp. Address before: Tokyo, Japan Patentee before: Toshiba Corp. Patentee before: TOSHIBA MEDICAL SYSTEMS Corp. |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101222 |