CN103156647A - 超声波诊断装置、图像处理装置及图像处理方法 - Google Patents
超声波诊断装置、图像处理装置及图像处理方法 Download PDFInfo
- Publication number
- CN103156647A CN103156647A CN2012105185486A CN201210518548A CN103156647A CN 103156647 A CN103156647 A CN 103156647A CN 2012105185486 A CN2012105185486 A CN 2012105185486A CN 201210518548 A CN201210518548 A CN 201210518548A CN 103156647 A CN103156647 A CN 103156647A
- Authority
- CN
- China
- Prior art keywords
- time series
- seasonal effect
- image
- blood flow
- dimensional
- 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
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
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/466—Displaying means of special interest adapted to display 3D data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
Abstract
实施方式涉及一种超声波诊断装置、图像处理装置及图像处理方法。准确且时间分辨率良好地观测二维或三维关注区域的血流动态。实施方式的超声波诊断装置包括处理部、图像生成部、检测部、控制部。处理部从由超声波收发而收集到的二维或三维的反射波数据,取得扫描范围的按照时间序列的二维或三维的血流信息,超声波收发是在由多条扫描线形成的扫描范围中进行的。图像生成部根据扫描范围的按照时间序列的二维或三维的血流信息,生成按照时间序列的血流图像。检测部在按照时间序列的血流图像之间,检测预先设定的二维或三维的关注区域内的斑点的按照时间序列的移动信息。控制部使规定的显示部显示基于按照时间序列的移动信息的数据即移动信息数据。
Description
本申请主张享有2011年12月8日申请的日本专利申请号2011-268884的优先权,并在本申请中引用上述日本专利申请的全部内容。
技术领域
实施方式涉及超声波诊断装置、图像处理装置及图像处理方法。
背景技术
以往,为了观察、诊断生物体的血流,而广泛应用超声波诊断装置。超声波诊断装置利用基于多普勒(Doppler)效应的多普勒法,使用超声波的反射波,生成并显示血流信息。作为由超声波诊断装置生成并显示的血流信息,有彩色(color)多普勒图像、多普勒波形(多普勒频谱(spectrum))等。
彩色多普勒图像是通过彩色血流图(CFM:Color Flow Mapping)法而拍摄的超声波图像。在CFM法中,针对包含观察部位、诊断部位的区域(二维区域或三维区域),在多条扫描线上收发超声波。而且,在CFM法中,通过MTI(Moving Target Indicator:动目标指示)滤波器(filter),从反射波数据(data)中,去除因组织运动引发的频率成分,提取血流成分的数据,并利用自相关法对血流成分的数据进行频率解析,由此运算出血流的速度、血流的分散、血流的能量(power)。彩色多普勒图像是将所述运算结果的分布以二维方式彩色显示所得的超声波图像。而且,在CFM法中,已知有用于提高低流速的血流的检测能力的被称为交替扫描(scan)的方法、能够以高帧频(frame rate)拍摄彩色多普勒图像的方法。
另一方面,多普勒波形是通过连续波(CW:Continuous Wave)多普勒法、脉冲(pulse)波(PW:Pulsed Wave)多普勒法而收集的数据。在CWD法、PWD法中,在包含观察部位、诊断部位的一根扫描线中,超声波的收发是依照时间序列进行的。在CWD法中,利用具有与MTI滤波器相同滤波器特性的HPF(High Pass Filter:高通滤波器),从按照时间序列的反射波数据中,提取扫描线上的血流成分的数据。而且,在PWD法中,是从在一根扫描线上所设定的取样容积(sample volume)(也被称为距离选通(range gate))内按照时间序列的反射波数据中,提取扫描线上的取样容积(例如1点)的血流成分的数据。
在CWD法、PWD法中,是通过快速傅里叶变换(FFT:Fast FourierTransform)法,对该血流成分的数据进行频率解析,由此运算出血流信息(血流的速度、血流的分散、血流的能量)。在CWD法中,输出的是扫描线上的平均血流信息,在PW多普勒法中,输出的是扫描线上的取样容积的血流信息。多普勒波形是依照时间序列将该运算结果的信息描绘(Plot)而成的波形的图像。多普勒波形的显示也被称为FFT显示。
在FFT显示中,由于扫描线为一根,因此能够准确且时间分辨率良好地观测一维关注区域内的血流动态。但是,在FFT显示中,无法观察到二维或三维关注区域的血流动态。
发明内容
本发明所要解决的课题在于,提供一种能够准确且时间分辨率良好地观测二维或三维关注区域内的血流动态的超声波诊断装置、图像处理装置及图像处理方法。
实施方式的超声波诊断装置具备处理部、图像生成部、检测部、控制部。处理部从由超声波收发而收集到的二维或三维的反射波数据,取得扫描范围的按照时间序列的二维或三维的血流信息,所述超声波收发是在由多条扫描线形成的扫描范围中进行的。图像生成部根据所述扫描范围的按照时间序列的二维或三维的血流信息,生成按照时间序列的血流图像。检测部在所述按照时间序列的血流图像之间,检测预先设定的二维或三维的关注区域内的斑点(speckle)的按照时间序列的移动信息。控制部使规定的显示部显示基于所述按照时间序列的移动信息的数据即移动信息数据。
根据实施方式的超声波诊断装置,可以准确且时间分辨率良好地观测二维或三维关注区域内的血流动态。
附图说明
图1是用于说明第1实施方式涉及的超声波诊断装置的构成例的图。
图2A、图2B及图2C是用于说明第1实施方式涉及的超声波诊断装置执行的各种彩色血流图法的扫描形态的图。
图3是用于说明第1实施方式涉及的CFM处理部的处理的一个例子的图。
图4是用于说明第1实施方式涉及的检测部的处理的一个例子的图。
图5是用于说明第1实施方式涉及的检测部的处理的一个例子的图。
图6A及图6B是用于说明第1实施方式涉及的检测部的处理的一个例子的图。
图7A、图7B是用于说明通过第1实施方式涉及的控制部而显示的移动信息数据的具体例子的图。
图8A、图8B用于说明通过第1实施方式涉及的控制部而显示的移动信息数据的具体例子的图。
图9是用于说明通过第1实施方式涉及的控制部而显示的移动信息数据的具体例子的图。
图10用于说明通过第1实施方式涉及的控制部而显示的移动信息数据的具体例子的图。
图11是用于说明在第1实施方式中设定多个关注区域时显示的移动信息数据的一个例子的图。
图12是用于说明第1实施方式涉及的超声波诊断装置的处理的流程图(flowchart)。
图13是用于说明生成多普勒波形时的PDW处理部及图像生成部的处理的图。
图14是用于说明第2实施方式涉及的图像处理部执行的处理次序的一个例子的图。
图15是用于说明第2实施方式涉及的图像处理部执行的处理次序的一个例子的图。
图16是用于说明第2实施方式涉及的图像处理部执行的处理次序的一个例子的图。
图17是用于说明第2实施方式涉及的图像处理部执行的处理次序的一个例子的图。
图18是用于说明第2实施方式涉及的超声波诊断装置的处理的流程图。
图19是用于说明第3实施方式的图。
图20是用于说明第3实施方式的图。
图21是用于说明第3实施方式涉及的超声波诊断装置的处理的流程图。
具体实施方式
以下,参照附图,详细说明超声波诊断装置的实施方式。
(第1实施方式)
首先,说明第1实施方式涉及的超声波诊断装置的构成例。图1是用于说明第1实施方式涉及的超声波诊断装置的构成例的图。如图1所示,第1实施方式涉及的超声波诊断装置具有超声波探头(probe)1、监视器(monitor)2、输入装置3、装置主体10。
超声波探头1具有多个压电振子,这些多个压电振子基于从下述的装置主体10具有的收发部11供给的驱动信号,而产生超声波。而且,超声波探头1接收来自被检测体P的反射波,将其变换成电信号。而且,超声波探头1具有设于压电振子的匹配层及声透镜(lens)、以及防止超声波从压电振子向后方传播的背衬(backing)材料等。超声波探头1装卸自由地与装置主体10连接。
若从超声波探头1向被检测体P发送超声波,则发送的超声波因被检测体P的体内组织的声阻抗(impedance)的非连续面而依次反射,并作为反射波信号被超声波探头1具有的多个压电振子接收。接收的反射波信号的振幅依赖于反射超声波的非连续面的声阻抗的差。此外,发送的超声波脉冲在由移动的血流、心脏壁等的表面反射时的反射波信号,由于多普勒效应,依赖于移动体针对超声波发送方向的速度成分,而受到频率偏移(多普勒偏移)。
此外,第1实施方式还可以应用于以下情况:利用作为将多个压电振子配置成一排的一维超声波探头的超声波探头1,二维地对被检测体P进行扫描(scan)的情况;利用作为使一维超声波探头的多个压电振子机械地摆动的超声波探头1或将多个压电振子格子状二维配置的二维超声波探头的超声波探头1,三维地对被检测体P进行扫描的情况。在此,一维超声波探头也可以利用一根扫描线,一维地对被检测体P进行扫描。而且,二维超声波探头还可以通过将超声波会聚后发送,而二维地扫描被检测体P。
而且,第1实施方式涉及的超声波诊断装置如后所述可以进行彩色多普勒图像的拍摄。而且,第1实施方式涉及的超声波诊断装置如后所述可以对B模式(mode)图像、彩色多普勒图像、B模式图像的一部分上重叠彩色多普勒图像而成的图像上设定的取样容积内的多普勒波形进行收集。因此,超声波探头1有时候会根据收集的图像的类别,例如将用于进行彩色血流图(CFM:Color Flow Mapping)法的超声波探头1,替换成用于进行连续波(CW:Continuous Wave)多普勒法或脉冲波(PW:Pulsed Wave)多普勒法的超声波探头1。此外,所述取样容积有时候也被称为距离选通。
输入装置3具有鼠标(mouse)、键盘(keyboard)、按钮(button)、面板开关(panel switch)、触控指令屏幕(touch command screen)、脚踏开关(foot switch)、轨迹球(trackball)等,从超声波诊断装置的操作者处接收各种设定请求,并将接收到的各种设定请求发送给装置主体10。
例如,输入装置3从操作者处接收后面所述的图像处理部15用于进行图像处理的关注区域(ROI:Region Of Interest)的设定。此外,在第1实施方式中,关于输入装置3接收的关注区域将在下文进行详细叙述。
监视器2显示GUI(Graphical User Interface:图形用户界面),或者显示在装置主体10上生成的超声波图像等,该GUI是用于由超声波诊断装置的操作者利用输入装置3输入各种设定请求。
装置主体10是基于超声波探头1接收到的反射波而生成超声波图像的装置。如图1所示,装置主体10具有收发部11、帧缓冲器(frame buffer)12、B模式处理部13、多普勒处理部14、图像处理部15、图像存储器(memory)16、控制部17、内部存储部18。
收发部11具有触发(trigger)产生电路、发送延迟电路及脉冲产生(pulsar)电路等,并向超声波探头1供给驱动信号。脉冲产生电路重复产生速率脉冲(rate pulse),该速率脉冲用于以规定的重复频率(PRF:Pulse Repetition Frequency)形成发送超声波。此外,PRF也被称为速率频率。而且,发送延迟电路对于脉冲产生电路产生的各速率脉冲,赋予将由超声波探头产生的超声波会聚成束状并决定发送指向性所需要的每个压电振子的发送延迟时间。而且,触发产生电路以基于速率脉冲的定时(timing),向超声波探头1施加驱动信号(驱动脉冲)。即,发送延迟电路通过改变对各速率脉冲赋予的发送延迟时间,任意调整来自压电振子面的发送方向。
此外,收发部11具有如下功能,即,基于后面所述的控制部17的指示,为了执行规定的扫描序列(scan sequence),而可瞬时变更发送频率、发送驱动电压等。尤其是,发送驱动电压的变更可以通过能够瞬时切换其值的线性放大器(linear amplifier)型的振荡电路、或可电切换多个电源单元(unit)的机构来实现。
而且,收发部11具有放大器(amplifier)电路、A/D(analog/digital)变换器、接收延迟电路、加法器、正交检波电路等,对超声波探头1接收到的反射波信号执行各种处理,生成反射波数据。放大器电路针对各通道(channel)放大反射波信号,进行增益(gain)校正处理。A/D变换器对增益校正后的反射波信号进行A/D变换。接收延迟电路对数字数据(digital data)赋予用于决定接收指向性所需的接收延迟时间。加法器对由接收延迟电路处理后的反射波信号进行加法处理。通过加法器的加法处理,而强调反射波信号的来自与接收指向性相应的方向的反射成分。然后,正交检波电路将加法器的输出信号变换成基频(baseband)频带的同相信号(I信号、I:In-pahse)和正交信号(Q信号、Q:Quadrature-phase)。然后,正交检波电路将I信号及Q信号(以下记做I/Q信号)作为反射波数据而保存到后级的帧缓冲器12中。这样,收发部11对超声波收发时的发送指向性和接收指向性进行控制。
帧缓冲器12是临时存储收发部11生成的反射波数据(I/Q信号)的缓冲器。具体来说,帧缓冲器12保持与存储容量相应的数据量的I/Q信号。例如,帧缓冲器12是FIFO(First-In/First-Out:先入先出)存储器,存储规定帧大小的I/Q信号,在收发部11新生成1个帧的I/Q信号时,废弃生成时间最久的1个帧的I/Q信号,存储新生成的1个帧的I/Q信号。
此外,所谓1个帧的I/Q信号,例如是指用于生成一张超声波图像的反射波数据,收发部11在多条扫描线(扫描线(scan line))形成的扫描范围内使超声波探头1进行超声波收发,从而生成1个帧的I/Q信号。
B模式处理部13从帧缓冲器12中,读出收发部11生成的反射波数据(I/Q信号),对读出的反射波数据进行对数放大、包络线检波处理、对数压缩等,生成以亮度的明亮度表现信号强度的数据(B模式数据)。
多普勒处理部14从帧缓冲器12中,读出收发部11生成的反射波数据(I/Q信号),通过对读出的反射波数据进行频率解析,而提取多普勒偏移(多普勒偏移(shift)频率),并通过使用多普勒偏移,提取基于多普勒效应的血流、组织、造影剂回声成分,从而生成关于多个点或1点提取平均速度、分散、能量等移动体信息而成的数据(多普勒数据)。
具体来说,如图1所示,多普勒处理部14具有CFM处理部14a和PWD处理部14b。CFM处理部14a是通过CFM法,生成用于生成彩色多普勒图像的多普勒数据的处理部。CFM处理部14a通过自相关法取得扫描范围内的血流的移动体信息(血流信息)。
而且,PWD处理部14b是通过PWD法生成用于生成多普勒波形的多普勒数据的处理部。PWD处理部14b通过进行频率解析,取得取样容积内的血流的移动体信息(血流信息)。例如,PWD处理部14b通过利用快速傅里叶变换(FFT:Fast Fourier Transform)法进行频率解析,而取得处于取样容积内的血流的移动体信息(血流信息)。此外,PWD处理部14b可以作为利用CWD法生成用于生成多普勒波形的多普勒数据的处理部发挥功能。而且,只要能够取得可生成多普勒波形的多普勒数据,则PWD处理部14b进行的频率解析也可以是快速傅里叶变换法以外的方法。
此外,关于CFM处理部14a进行的处理、PWD处理部14b进行的处理,将于下文进行详细叙述。而且,多普勒处理部14也可以是具有用于进行组织多普勒法的处理部的情形。
图像处理部15是如下处理部,即,使用B模式处理部13及多普勒处理部14生成的数据,生成显示用的图像数据,并对生成的图像数据进行图像处理。图1所示的图像处理部15具有图像生成部15a及检测部15b。
图像生成部15a利用B模式处理部13及CFM处理部14a生成的数据,生成超声波图像。即,图像生成部15a利用B模式处理部13生成的B模式数据,生成以亮度表现反射波的强度的B模式图像。而且,图像生成部15a利用CFM处理部14a生成的多普勒数据,生成作为表示移动体信息(血流信息)的速度图像、分散图像、能量图像、或者,这些图像的组合图像的彩色多普勒图像。例如,图像生成部15a生成与能量值相应地以红色系使色调变化的能量图像。此外,图像生成部15a除了可以生成彩色显示用的彩色多普勒图像以外,例如还可以生成与能量值相应地以灰度(grayscale)方式使亮度变化的灰度的能量图像。
以下,将图像生成部15a利用CFM处理部14a生成的数据而生成的彩色多普勒图像等图像记做“血流图像”。
进而,图像生成部15a利用PWD处理部14b生成的多普勒数据,生成将血流的速度信息依照时间序列进行描绘而成的多普勒波形。具体来说,图像生成部15a生成以取样容积内的血流的速度为纵轴、以时间为横轴的时间变化曲线。然后,图像生成部15a根据取样容积内的血流的分散值,设定纵轴方向的宽度,并根据取样容积内的血流的能量值来设定亮度值,从而生成多普勒波形。
这里,图像生成部15a一般来说是将超声波扫描的扫描线信号流,变换(扫描转换(scan convert))成以电视(television)等为代表的视频格式(video format)的扫描线信号流,生成作为显示用图像的超声波图像(B模式图像、血流图像)。具体来说,图像生成部15a通过进行与基于超声波探头1的超声波的扫描形态相应的坐标变换,生成作为显示用图像的超声波图像。而且,图像生成部15a除了进行扫描转换以外,作为各种图像处理,例如还进行使用扫描转换后的多个图像帧,重新生成亮度的平均值图像的图像处理(平滑化处理)、在图像内使用微分滤波器的图像处理(边缘(edge)增强处理)等。
此外,图像生成部15a在以二维方式进行了超声波收发时,通过进行坐标变换,而生成作为显示用图像的二维B模式图像或二维血流图像。而且,图像生成部15a在以三维方式进行了超声波收发时,生成体数据(volume data)(三维B模式图像或三维血流图像),并通过各种绘制(rendering)处理,利用体数据生成用于显示在监视器2上的二维图像。
而且,图像生成部15a在各种图像上合成各种参数(parameter)的文字信息、刻度、体位标志(body mark)等,生成合成图像。而且,图像生成部15a生成B模式图像和彩色多普勒图像的重叠图像等重叠了各种图像而成的重叠图像,或者生成用于并列显示各种图像的图像。
图1所示的检测部15b是对图像生成部15a生成的图像数据进行图像处理的处理部。具体来说,检测部15b进行对图像内的斑点进行追踪的斑点追踪(Speckle Tracking)。此外,在第1实施方式中,关于检测部15b进行的处理内容将于下文详细叙述。
图像存储器16是存储图像生成部15a生成的各种数据的存储器。而且,图像存储器16也可以存储B模式处理部13、多普勒处理部14生成的数据(原始数据)。而且,图像存储器16还可以视需要存储帧缓冲器12所保持的数据。
内部存储部18存储用于进行超声波收发、图像处理及显示处理的控制程序(program)、诊断信息(例如患者ID、医生意见等)、诊断协议(protocol)、各种体位标志等各种数据。而且,内部存储部18还可以视需要被用于图像存储器16所存储的数据的保管等中。而且,内部存储部18所存储的数据可以通过未图示的接口(interface),转送给外部的周边装置。
控制部17控制超声波诊断装置的全体处理。具体来说,控制部17基于操作者通过输入装置3输入的各种设定请求、从内部存储部18读取的各种控制程序及各种数据,控制收发部11、B模式处理部13、多普勒处理部14、图像处理部15的处理。而且,控制部17以在监视器2上显示图像存储器16存储的数据、用于操作者指定各种处理的GUI等的方式进行控制。此外,利用PWD法、DWD法而显示多普勒波形的操作,也被称为“FFT显示”。而且,以下,将利用CFM法显示血流图像(彩色多普勒图像)的操作,记做彩色多普勒显示。
以上,说明了第1实施方式涉及的超声波诊断装置的全体构成。基于该构成,第1实施方式涉及的超声波诊断装置通过CFM法而生成血流图像。在此,第1实施方式涉及的超声波诊断装置可以利用图2A、图2B及图2C例示的各种扫描形态,执行各种CFM法。图2A、图2B及图2C是用于说明第1实施方式涉及的超声波诊断装置执行的各种彩色血流图法的扫描形态的图。
在以往的彩色多普勒显示中,对于一根扫描线(扫描线)要进行两次以上的超声波收发。以下,将以往的彩色多普勒显示中进行的CFM法记做“以往的CFM法”。在以往的CFM法中,在同一方向上重复多次超声波收发,将这些多次的数据(反射波数据)设为一个分组(packet),在分组内进行封闭处理,并关于1分组输出1个数据(多普勒数据)。若将1帧内的扫描线数设为“M”、将与在每一根扫描线进行的收发次数对应的分组大小(packet size)设为“N”、将重复频率设为“PRF”,则即便在不进行B模式用扫描的情况下,在以往的CFM法中,帧频“Fr”也能由以下的式(1)表示。
若设“PRF=5000(Hz)、M=50、N=10”,根据式(1),帧频变成10fps(frame per second)。即,在以往的CFM法中,随着分组大小“N”增大,帧频将变慢。
图2A是表示设“M=13、N=8”时,以往的CFM法进行的扫描形态的图。此外,以下将图2A、图2B及图2C所示的13根扫描线从右向左依次设为L1~L13。在图2A所示的一例中,以L1进行“第1次、第2次、…、第8次”的超声波收发,以L2进行“第9次、第10次、…、第16次”的超声波收发,以L3进行“第17次、第18次、…、第24次”的超声波收发,以L4进行“第25次、第26次、…、第32次”的超声波收发。并且,在图2的(A)所示的一例中,以L5进行“第33次、第34次、…、第40次”的超声波收发,以L6进行“第41次、第42次、…、第48次”的超声波收发,以L7进行“第49次、第50次、...、第56次”的超声波收发,以L8进行“第57次、第58次、...、第64次”的超声波收发。同样地,进行L9~L13的超声波收发。
另一方面,作为与以往的CFM法相比使低流速检测能力提高的CFM法,已实用化的有被称为交替扫描的方法。在此,若将超声波的音速设为“C”、将接收中心频率设为“f0”,则作为不会引起混叠的最大可检测流速(混叠速度)的“Vmax”通过以下的式(2)表示。
在使用相同特性的MTI(Moving Target Indicator)滤波器的情况下,混叠速度小则低流速检测能力高。在交替扫描法中,将多条扫描线设为一个块(block),在块内交替进行超声波收发。若将块内的扫描线数(光栅(raster)数)设为“L”,则用于输出块内的所有分组的数据的时间成为以往的CFM法的“L”倍。因此,在光栅数为“L”的交替扫描法中,“Vmax”是通过以下的式(3)表示。
图2B是表示设“M=13、N=8、L=4”时利用交替扫描法进行的扫描形态的图。在“N=8、L=4”时的交替扫描法中,针对4根扫描线的每根扫描线的1次超声波收发,通过在4根扫描线中交替重复8次,由此每一个块合计进行32次超声波收发。例如,设“L1~L4”为“块1”、“L5~L8”为“块2”进行说明。该情况下,在“块1”中,第1次至第4次的超声波收发是按照“L1、L2、L3、L4”的顺序进行,第5次至第8次的超声波收发是按照“L1、L2、L3、L4”的顺序进行,最后第29次至第32次的超声波收发是按照“L1、L2、L3、L4”的顺序进行。另外,在“块2”中,第33次至第36次的超声波收发是按照“L5、L6、L7、L8”的顺序进行,第37次至第40次的超声波收发是按照“L5、L6、L7、L8”的顺序进行,最后第61次至第64次的超声波收发是按照“L5、L6、L7、L8”的顺序进行。
其中,若“PRF、M、N”的值相同,则在以往的CFM法和交替扫描法中帧频是相同的。
以下,再次说明所述以往的CFM法、交替扫描法中实际进行的处理,说明帧频为式(1)的理由。在CFM法中,通过MTI滤波器,除去组织运动所致的频率成分,从而从反射波数据中提取血流成分的数据。并且,在CFM法中,利用自相关法根据提取到的数据运算“速度·分散·能量”,并显示运算结果的二维分布。在以往的CFM法、交替扫描法中,对同一光栅进行多次收发。
同一光栅上的同一地点的接收数据(反射波数据)的集合被称为分组,CFM处理部14a对该分组进行MTI滤波器处理、自相关运算、“速度·分散·能量推断”这样的处理。在此,分组的大小通常被设为“6~20”。CFM处理部14a最终在分组内生成1组“速度·分散·能量”的数据。CFM处理部14a在整个二维空间进行该处理,从而生成二维的多普勒数据。而且,例如在利用机械式扫描探头进行三维扫描的情况下,CFM处理部14a在多张二维空间中进行该二维空间的处理,从而生成三维的多普勒数据。借此,显示血流图像。
这样,在以往的CFM法、交替扫描法中,以分组为单位而生成1组数据。进而,在进行CFM法时,超声波探头1还进行用于在B模式图像上重叠血流图像的B模式用扫描。因这两个理由,血流的帧频变成5~30fps这样的较慢帧频。在以往的CFM法、交替扫描法中,由于受到这种时间分辨率的制约,所以难以生成显示描绘出瞬时血流动态的血流图像。
因此,开发了一种用于提高帧频的CFM法。以下,将该CFM法记做“高帧频法”。在高帧频法中,并不是以分组为单位进行MTI滤波器处理、自相关运算处理、“速度·分散·能量推断”处理,而是进行与B模式用扫描相同的超声波扫描。在高帧频法中,例如,如图2C所示,在形成1帧的扫描范围的13根扫描线(L1~L13)的各个中,逐次进行超声波收发。并且,在高帧频法中,在帧方向上对各帧的相同位置的数据流进行处理。由此,在高帧频法中,MTI滤波器处理可以从对分组这种有限长度的数据进行处理变成对无限长度的数据进行处理,从而可提高MTI滤波器的性能,同时能够以与扫描帧频相同的帧频显示血流信息。即,在高帧频法中,能够将MTI滤波器处理变成无限冲击响应型滤波器(IIR滤波器、IIR:Infinite Impulse Response)处理。由此,在高帧频法中,相对于式(1)而言可以实现“N”倍的帧频。
这样,第1实施方式涉及的超声波诊断装置可以通过各种扫描形态来执行CFM法。例如,可以通过高帧频法进行彩色多普勒显示,从而操作者能准确且时间分辨率良好地观测瞬时的血流动态。
相对于此,在FFT显示中,由于扫描线为一根,因此时间分辨率高于彩色多普勒显示。而且,FFT的频率解析是FFT处理,精度高于自相关法的频率解析。由此,为了准确且时间分辨率良好地观测瞬时血流动态,FFT显示优于彩色多普勒显示。其中,若将超声波束与血流方向所成的角度设为“θ”,观测到的多普勒偏移频率变成真正的多普勒偏移频率的“cosθ”倍的值,因此在CWD法、PWD法中,通过赋予“θ”而对测量速度施加校正(角度校正)。
但是,FFT显示如扫描线上、例如1点的取样容积这样,只能显示一维的血流动态。即,在FFT显示中可观察到的血流动态为一维的关注区域内的血流动态。
另一方面,在彩色多普勒显示中可以在二维空间内显示血流动态。而且,在进行三维扫描时,在彩色多普勒显示中可以在三维空间内显示血流动态。然而,彩色多普勒显示如上所述般,和FFT显示相比,帧频不够快。而且,以二维或三维进行角度校正时,操作者必须对各场所设定“θ”,因此在彩色多普勒显示中通常不进行角度校正。
因这些理由,彩色多普勒显示被用于观察二维或三维的关注区域的血流动态,并不能代替FFT显示。但是,如上述那样,FFT显示可观察到的血流动态为一维的关注区域内的血流动态,要观察二维或三维的关注区域的血流动态,需要例如在二维空间或三维空间设定多个取样容积,针对多个取样容积进行PDW法的多普勒波形的收集。
因此,在第1实施方式中,为了准确且时间分辨率良好地观测二维或三维关注区域内的血流动态,而进行以下处理。即,CFM处理部14a利用在通过多条扫描线形成的扫描范围内进行的超声波收发而收集到的二维或三维的反射波数据,取得扫描范围的按照时间序列的二维或三维的血流信息。具体来说,CFM处理部14a利用自相关法而取得扫描范围的按照时间序列的二维或三维血流信息。图像生成部15a利用扫描范围的按照时间序列的二维或三维血流信息,生成按照时间序列的血流图像。检测部15b在按照时间序列的血流图像之间,检测预先设定的二维或三维关注区域内的斑点的按照时间序列的移动信息。控制部17在监视器2上显示基于按照时间序列的移动信息的数据即移动信息数据。
以下,详细说明所述处理的一个例子。此外,以下说明的情况是,收发部11通过高帧频法的扫描形态,执行超声波收发,CFM处理部14a基于高帧频法而取得血流信息。即,以下,CFM处理部14a利用通过在各扫描线进行一次扫描范围内的超声波收发的扫描形态而收集到的二维或三维反射波数据,取得扫描范围的按照时间序列的二维或三维血流信息。其中,第1实施方式中,根据扫描范围的大小,例如如果能将帧频设为60fps以上的高速,则即使在执行所述以往的CFM法或交替扫描法的情况下也是可以的。
首先,CFM处理部14a从帧缓冲器12中,读出通过在多条扫描线形成的扫描范围内进行的超声波收发而收集到的二维或三维反射波数据。然后,CFM处理部14a针对二维或三维反射波数据进行自相关法,从而取得扫描范围的按照时间序列的二维或三维血流信息。例如,根据控制部17指定的扫描序列,收发部11使超声波探头1执行高帧频法的扫描形态,即,在各扫描线进行一次二维的扫描范围内的超声波收发,并接收1帧的反射波信号。由此,收发部11生成二维的反射波数据,并将其保存到帧缓冲器12内。CFM处理部14a针对二维的反射波数据执行自相关法,从而按照时间序列取得扫描范围内的二维血流信息。在此,血流信息为血流的速度、分散及能量。图3是用于说明第1实施方式涉及的CFM处理部及图像生成部的处理的一个例子的图。
如图3所示,CFM处理部14a的处理包含“MTI滤波器处理”、“自相关运算处理”及“速度·分散·能量推断处理”这3个块。在图3所示的“MTI滤波器处理”块中,例如,作为MTI滤波器,使用4阶的IIR滤波器。
在此,将当前帧为“第n帧”中的“某个位置”的反射波数据设为“x(n)”、将相同位置的前1帧为“第(n-1)帧”的反射波数据设为“x(n-1)”、将相同位置的前2帧为“第(n-2)帧”的反射波数据设为“x(n-2)”。而且,将相同位置的前3帧为“第(n-3)帧”的反射波数据设为“x(n-3)”、将相同位置的前4帧为“第(n-4)帧”的反射波数据设为“x(n-4)”。而且,将MTI滤波器的滤波器系数设为“ak”及“bk”。该情况下,从“MTI滤波器处理”块输出的“x(n)”的血流成分即“y(n)”是通过以下的式(4)表示。
“MTI滤波器处理”块的输出结果被输入至图3所示的“自相关运算处理”块。在“自相关运算处理”块中,进行滞后(lag)0及滞后1的自相关值的运算。若将滞后0设为“C0(n)”、滞后1设为“C1(n)”,则“C0(n)”通过以下的式(5)算出,“C1(n)”通过以下的式(6)算出。
此外,在式(5)及式(6)中,上标星号(asterisk)“*”表示复共轭(complex conjugate),“N”表示自相关运算中设定的帧方向的加法幅度。“C0(n)”和能量等价。此外,图3中虽未表示,但通过对式(5)及式(6)获得的滞后0及滞后1,在超声波的波束方向进行移动平均等平滑化,而可牺牲距离分辨率以提高信噪比(S/N)。
“自相关运算处理”块的输出结果被输入至图3所示的“速度·分散·能量推断处理”块。在“速度·分散·能量推断处理”块中,利用“C0(n)”及“C1(n)”,通过以下所示的式(7),算出“第n帧”的“某个位置”的血流的速度“V”、血流的速度的分散“T”和血流的能量“P”。
此外,在式(7)的第1式中,“atan2”是输出“-π~+π”的角度的“反正切(arc tangent)函数”,“image”是从复数中仅输出虚部的函数,“real”是从复数中仅输出实部的函数。即,在式(7)的第1式中,算出“C1(n)”的虚部和实部的相位变化作为血流的速度“V”。这样,能量通过自相关运算作为滞后0而被求出,速度及分散则通过对自相关运算的结果进行频率解析而被求出。
通过所述处理,CFM处理部14a针对每一帧取得扫描范围内的各扫描线上设定的多个点各自的血流信息,并将其输出至图像生成部15a。图像生成部15a利用扫描范围的按照时间序列的二维或三维血流信息,生成按照时间序列的血流图像。在第1实施方式中,图像生成部15a利用二维的血流信息而生成按照时间序列的血流图像。具体来说,如图3所示,图像生成部15a通过进行“坐标变换处理”,而利用按照时间序列的二维的血流信息,生成按照时间序列的二维血流图像。
例如,图像生成部15a生成以灰度描绘能量值的灰度的能量图像。或者,图像生成部15a生成如下彩色的能量图像,即,仅将速度符号作为方向来使用,例如在符号为正(plus)的某个位置上根据能量值使红色系的明亮度变化,在符号为负(minus)的位置上根据能量值使蓝色系的明亮度变化。或者,图像生成部15a生成以二维图(map)变换速度和能量而得的彩色的能量图像。此外,图像生成部15a不仅生成能量显示用的血流图像,有时也可以生成速度显示用血流图像、速度-分散显示用血流图像。其中,作为以下说明的检测部15b的处理中使用的血流图像,理想的是能量显示用血流图像,因此第1实施方式涉及的图像生成部15a即使在仅要求进行速度显示、速度-分散显示的情况下,也会生成能量显示用血流图像。
图1所示的检测部15b在按照时间序列的血流图像之间,检测预先设定的二维或三维关注区域内的斑点的按照时间序列的移动信息。具体来说,检测部15b检测关注区域内的多个点各自的按照时间序列的移动信息。首先,说明关注区域的设定例子。
例如,控制部17根据操作者的ROI设定请求,在监视器2上显示图像生成部15a生成的一个血流图像。操作者参照监视器2,使用输入装置3,设定希望观察血流动态的关注区域。在生成显示二维血流图像的情况下,操作者设定二维的关注区域。而且,在生成三维血流图像并显示三维血流图像的一截面的情况下,操作者设定二维关注区域。或者,在生成三维血流图像并显示例如三维血流图像的正交3截面的情况下,操作者通过在各截面设定二维关注区域,而设定三维关注区域。在第1实施方式中,说明了操作者对二维血流图像设定二维关注区域的情况。图4、图5、图6A及图6B是用于说明第1实施方式涉及的检测部的处理的一个例子的图。
例如,如图4所示,操作者对灰度的能量图像设定矩形的ROI。此外,第1实施方式中,ROI的形状可以任意设定。而且,第1实施方式中,也可以在包含彩色多普勒显示用扫描范围的扫描范围内的B模式图像上,设定用于检测部15b的处理的ROI。
并且,如图5所示,检测部15b主要进行“斑点追踪处理”和“速度向量(vector)计算处理”这两个处理。在“斑点追踪处理”块中,检测部15b利用相互相关法,在连续的帧间、或相隔数帧的帧间进行斑点的追踪(tracking),从而计算ROI内的多个点的运动向量。并且,检测部15b在“速度向量计算处理”块中,使用帧间时间差,将运动向量变换成速度向量。
例如,如图6A所示,在ROI内设定35个追踪点。检测部15b通过追踪“第(n-1)帧”的追踪点中的斑点在“第n帧”移动到哪个位置,而计算运动向量。并且,检测部15b在图5所示的“速度向量计算处理”中,使用作为相互相关处理对象的两个帧的帧间时间差,利用运动向量计算出速度向量。例如,如图6B所示,检测部15b计算出“第n帧”的ROI内的35点的各自的速度向量。通过在血流图像间重复进行该处理,检测部15b按照时间序列而计算出ROI内的多个点各自的速度向量。
在此,在远远小于发送超声波波长的反射体密集的情况下,反射波信号会相互干涉。干涉的强弱会变成反射波信号的振幅的强弱,基于该振幅信息的超声波图像内会产生点状的伪影(artifact)的斑点。能量图像是将表示作为反射体的血液的量的程度的能量值图像化而成的血流图像。因此,从能量图像的斑点的移动向量变换的速度向量变成与血流的速度向量对应的信息。即,在第1实施方式中,具有能量信息的血流图像的ROI内的多个点各自的基于斑点的速度向量,变成表示血流动态的“移动信息”。
控制部17在监视器2上显示作为基于按照时间序列的移动信息的数据的移动信息数据。以下,针对通过第1实施方式涉及的控制部17的控制而显示的移动信息数据,使用图7~图10等进行具体说明。图7~图10是用于说明通过第1实施方式涉及的控制部而显示的移动信息数据的具体例子的图。
具体来说,控制部17使图像生成部15a,生成按照时间序列将多个点各自的按照时间序列的移动信息的代表值描绘而成的时间变化曲线,作为移动信息数据。并且,控制部17在监视器2上显示时间变化曲线。移动信息的代表值为平均速度。在此,作为代表值的平均速度例如是多个点各自的速度向量的标量(scalar)的平均值。图7A中表示图像生成部15a生成以帧序号(frame number)为横轴、以平均速度(单位:cm/s)为纵轴的平均速度时间变化曲线时的一个例子。而且,图7B表示图像生成部15a生成以使用帧间时间差而从帧序号换算来的时间(单位:sec(秒))为横轴、以平均速度(单位:cm/s)为纵轴的平均速度时间变化曲线时的一个例子。用于横轴的信息是操作者可以任意变更的参数。
此外,作为代表值的平均速度既可以由检测部15b计算出,也可以由控制部17计算出。而且,第1实施方式中,还可以将作为代表值的平均速度设为多个点各自的速度向量的设定方向成分的标量的平均值。此外,设定方向成分例如是方位方向、深度方向、血管的行走方向等。
而且,作为代表值的平均速度,有时候也可以是平均速度向量的标量、平均速度向量的设定方向成分的标量。此外,设定方向成分例如为方位方向、深度方向、血管的行走方向等。
而且,代表值也可以是最大流速。作为代表值的最大流速是速度向量的标量的最大值、或速度向量的设定方向成分的标量的最大值。此外,设定方向成分例如为方位方向、深度方向、血管的行走方向等。
而且,代表值也可以是流速的中央值。作为代表值的中央值是速度向量的标量的中央值、或速度向量的设定方向成分的标量的中央值。此外,设定方向成分例如为方位方向、深度方向、血管的行走方向等。
而且,代表值也可以是分散值。作为代表值的分散值是速度向量的标量的分散值、或速度向量的设定方向成分的标量的分散值。此外,设定方向成分例如为方位方向、深度方向、血管的行走方向等。
而且,代表值还可以是多个点各自的速度向量的标量的直方图(histogram)中成为最大频度、“第n高的频度”的区间的值。例如,在“50cm/s~60cm/s”的区间为最大频度的情况下,代表值可以是“55cm/s”。或者,为了避免使用因噪音而变成最大频度的区间来计算出代表值,也可以利用“第3高的频度”的区间来计算出代表值。在利用直方图计算出代表值的情况下,有时候也可以使用多个点各自的速度向量的设定方向成分的标量的直方图。此外,设定方向成分例如为方位方向、深度方向、血管的走行方向等。
此外,控制部17在从所述各种代表值中选择多个代表值的情况下,使图像生成部15a生成被选择出的每多个代表值的时间变化曲线,并显示在监视器2上。
在此,由于图7A、图7B所示的平均速度的时间变化曲线类似于多普勒波形,因此能够用于代替FFT显示。于是,控制部17有时也可以使图像生成部15a生成将时间变化曲线图像化而成的时间变化图像,作为移动信息数据。即,如图8A所示,图像生成部15a生成绘制了利用速度向量计算出的ROI内的速度的平均值的时间变化曲线(参照图中的实线),根据ROI内的速度的分散值而设定纵轴方向的宽度(参照图中的虚线)。并且,如图8A所示,图像生成部15a根据ROI内的能量的平均值,设定以速度的平均值为中心且纵轴方向的长度为速度的分散值的线段的亮度,由此生成时间变化图像。
由此,图像生成部15a生成图8B例示的时间变化图像,控制部17在监视器2上显示时间变化图像。图8B所示的平均速度的时间变化图像是比平均速度的时间变化曲线更类似于多普勒波形的图像。
或者,控制部17使图像生成部15a生成在血流图像的关注区域内重叠了表示多个点各自的移动信息的规定图形而成的重叠图像,作为移动信息数据。例如,如图9所示,图像生成部15a生成在能量图像的ROI内,重叠了表示多个点各自的速度向量的箭头而成的重叠图像。
或者,控制部17还可以使图像生成部15a生成如下图像作为移动信息数据,该图像是将代表值及表示作为该代表值的计算依据的多个点各自的移动信息的值的分布的直方图中的至少一个重叠在重叠图像的规定位置上而形成的图像。例如,如图9所示,图像生成部15a在重叠图像内的ROI上合成作为代表值的平均速度“67cm/s”。另外,例如图像生成部15a如图9所示,在代表值左侧合成表示作为平均速度的计算依据的35个追踪点各自的速度向量的标量的分布的直方图。此外,代表值、直方图被合成的位置除了在重叠图像上以外,还可以在重叠图像的上方。
在此,血流图像是按照时间序列而生成了多个,因此能够动画显示。而且,在能量图像的ROI内所重叠的箭头是表示该能量图像中计算出的多个点各自的速度向量的箭头。而且,对重叠图像进一步重叠的代表值、直方图是与该重叠图像对应的能量图像的ROI内的代表值和直方图。即,重叠图像、或重叠了代表值、直方图的重叠图像均是按照时间序列而生成了多个,因此能够动画显示。
于是,控制部17在与时间变化曲线或时间变化图像并列地动画显示重叠图像的情况下,会在时间变化曲线或时间变化图像中,显示与显示的重叠图像的时相相对应的位置。例如,控制部17在与图7A所示的平均速度的时间变化曲线并列地动画显示重叠了代表值及直方图的重叠图像的情况下,进行图10所例示的显示控制。
即,如图10所示,控制部17将表示计算出监视器2上显示的重叠图像的平均速度的时刻的位置的线标记(line marker)a,合成到平均速度的时间变化曲线上予以显示。在此,图10所示的线标记a的位置根据所显示的血流图像的帧而从左向右移动。此外,在第1实施方式中,控制部17在与时间变化曲线或时间变化图像并列地动画显示速度显示、速度-分散显示的彩色多普勒图像的情况下,也可以将与所显示的血流图像的时相对应的位置,显示在时间变化曲线或时间变化图像中。
在所述一个例子中说明了基于在二维的ROI中检测到的移动信息的各种移动信息数据被显示的情况,但第1实施方式中,也可以显示基于在三维的ROI中检测到的移动信息的各种移动信息数据。
而且,由于血流图像是二维图像,因此操作者可以将二维的ROI设定在任意位置上。而且,在三维拍摄血流图像的情况下,血流图像是以拍摄空间的任意截面二维显示,因此操作者可以将二维或三维的ROI设定在任意位置上。所以,在第1实施方式中,也可以设定多个ROI。在设定多个ROI的情况下,通过控制部17的控制,检测部15b及图像生成部15a并行地进行针对多个ROI各自的处理。图11是用于说明第1实施方式中设定多个关注区域时所显示的移动信息数据的一个例子的图。
例如,在希望更细微地观察图5所例示的ROI内的血流动态的情况下,如图11所示,操作者将图5所例示的ROI分割为两部分而设定ROI-A及ROI-B。这种情况下,检测部15b计算出ROI-A内的多个点各自的速度向量、和ROI-B内的多个点各自的速度向量。并且,例如,如图11所示,图像生成部15a生成ROI-A及ROI-B各自的平均速度的时间变化曲线。此外,图11所例示的ROI-A及ROI-B的时间变化曲线的横轴为“时间(单位:sec(秒))”。
并且,如图11所示,图像生成部15a将表示ROI-A内的多个点各自的速度向量的箭头、和表示ROI-B内的多个点各自的速度向量的箭头,重叠于能量图像。另外,如图11所示,图像生成部15a将ROI-A内的平均速度和速度的直方图、以及ROI-B内的平均速度和速度的直方图重叠。并且,如图11所示,图像生成部15a生成合成图像,该合成图像是使两个ROI的信息重叠于同一能量图像的重叠图像、以及两个时间变化曲线合成而成。该合成图像通过控制部17而被监视器2显示。此外,图11中虽未图示,但在动画显示重叠图像的情况下,也可以将表示与所显示的重叠图像的时相对应的位置的线标记,分别合成到两个时间变化曲线上。
此外,移动信息数据的类别设定、设定多个类别的移动信息数据时的显示形态既可以由操作者事先设定,也可以初始设定。
接下来,使用图12,说明第1实施方式涉及的超声波诊断装置的处理。图12是用于说明第1实施方式涉及的超声波诊断装置的处理的流程图。此外,在图12中,说明的是如下情况,即,在收集血流图像后进行ROI设定,收集血流图像后,在血流图像的电影(cine)再生过程中并行进行所述检测部15b及控制部17的处理。
如图12所示,第1实施方式涉及的超声波诊断装置的控制部17判定是否从操作者接受了关注区域设定(步骤(step)S101)。在此,在未接受到关注区域的设定的情况下(步骤S101,否),控制部17待机至接受到关注区域的设定为止。
另一方面,在接受到关注区域的设定的情况下(步骤S101,是),控制部17将关注区域的位置信息通知给检测部15b,检测部15b在血流图像的各帧中计算出关注区域内的多个点的运动向量(步骤S102),并将计算出的运动向量换算成速度向量(步骤S103)。
并且,通过控制部17的控制,图像生成部15a生成附线标记的时间变化曲线、附线标记的时间变化图像、重叠图像等移动信息数据(步骤S104)。并且,通过控制部17的控制,监视器2显示移动信息数据(步骤S105),并结束处理。
此外,在第1实施方式中,说明的是如下情况,即,收集血流图像后进行ROI设定,并且收集血流图像后,在血流图像的电影再生过程中进行所述检测部15b及控制部17的处理。但是,第1实施方式也可以是如下情况,即,收集血流图像前进行ROI设定,且与收集血流图像并行地实时进行所述检测部15b及控制部17的处理。
如上所述,在第1实施方式中,检测血流图像内设定的ROI内的速度向量,并将检测到的速度向量的信息,以时间变化曲线或时间变化图像、重叠图像等形态来显示。即,在第1实施方式中,通过在血流图像内任意位置设定任意形状的ROI,能够将与在多普勒显示下所观察的一维血流动态相关的信息同样的信息,以二维或三维来进行显示。因此,在第1实施方式中,可以准确且时间分辨率良好地观测二维或三维关注区域内的血流动态。而且,在第1实施方式中,可以设定多个关注区域。因此,在第1实施方式中,可以同时观测二维或三维的多个关注区域内的血流动态。此外,移动信息的检测优选是以能量图像为对象来进行的,因此血流信息有时也可以仅为通过自相关运算作为滞后0而被求出的能量。即,第1实施方式也可以是如下情况,即,CFM处理部14a利用自相关法仅取得能量,图像生成部15a仅生成按照时间序列的能量图像。
(第2实施方式)
在第2实施方式中,说明了与第1实施方式中说明的移动信息数据一起,另行地显示通过多普勒显示用扫描而生成的多普勒波形。
即,在第2实施方式中,PDW处理部14b针对通过在第1实施方式中说明的关注区域所含的扫描线利用另行地进行的由超声波收发所收集到的反射波数据,利用例如基于快速傅里叶变换法的频率解析,从而取得扫描线上的取样容积内的按照时间序列的血流信息。并且,在第2实施方式中,图像生成部15a根据通过频率解析取得的按照时间序列的血流信息,而生成多普勒波形。
首先,PDW处理部14b从帧缓冲器12中,读取通过脉冲波收发而生成的一维反射波数据之中、与扫描线上的取样容积的位置相对应的一维反射波数据。或者,也可以是如下情况,即,收发部11使用较短的时间门,仅将与扫描线上的取样容积的位置相对应的一维反射波数据保存到帧缓冲器12内,PDW处理部14b读出与扫描线上的取样容积的位置相对应的一维反射波数据。图13是用于说明生成多普勒图像时的PDW处理部及图像生成部的处理的图。
首先,PDW处理部14b在图13所示的“HPF处理”块内,对一维反射波数据进行具有与MTI滤波器相同特性的HPF(High Pass Filter)处理,从而除去低频成分(杂波:clutter),提取血流成分。并且,PDW处理部14b在图13所示的“FFT处理”块内,通过快速傅里叶变换法,利用血流成分进行血流信息(速度·分散·能量)的运算。此外,在图13所示的“FFT处理”块内,还进行所述角度校正。并且,图像生成部15a在图13所示的“多普勒波形生成处理”块内,生成多普勒波形。
图14~图17是用于说明第2实施方式涉及的图像处理部执行的处理次序的一个例子的图。例如,与第1实施方式同样地,在第2实施方式中,如图14所示,对能量图像设定ROI。并且,在第2实施方式中,图像生成部15a使用检测部15b的处理结果,如图14所示,生成例如平均速度的时间变化曲线等移动信息数据。
另外,在第2实施方式中,例如,在图15的右图所示的B模式图像上,设定ROI内所含的取样容积(参照图中的虚线内)。并且,在第2的实施方式中,在与血流图像的拍摄不同的时期,在取样容积上的扫描线中进行多普勒显示用扫描。由此,图像生成部15a生成图15的左图例示的多普勒波形。
并且,第2实施方式涉及的控制部17将多普勒波形和移动信息数据,以能够明确表示两数据的同一时相的形态,显示在监视器2上。例如,在并列显示多普勒波形和时间变化曲线的情况下,如图16所示,控制部17判定多普勒波形中波形为峰值的时期、和时间变化曲线中代表值为峰值的时期对应于大致同一时相。控制部17为了让操作者容易视觉辨认两者达到峰值(peak)的时相,例如,如图16所示,显示将两者的峰值连结的箭头。或者,控制部17以达到峰值的时相在横轴方向上大体一致的方式,变更例如多普勒波形的横轴的显示比例尺(scale)、或变更时间变化曲线的横轴的显示比例尺。
而且,第2实施方式涉及的控制部17在与多普勒波形并列地动画显示重叠图像的情况下,在多普勒波形中显示与所显示的重叠图像的时相相对应的位置。而且,第2实施方式涉及的控制部17在与多普勒波形及时间变化曲线并列地动画显示重叠图像的情况下,在多普勒波形及时间变化曲线中显示与所显示的重叠图像的时相相对应的位置。
例如,如图17所示,控制部17在与多普勒波形及平均速度的时间变化曲线并列地动画显示重叠有作为代表值的平均速度及直方图的重叠图像的情况下,与第1实施方式同样地,在平均速度的时间变化曲线上合成表示计算出监视器2所显示的重叠图像的平均速度的时刻的位置的线标记a。进而,如图17所示,控制部17在多普勒波形上合成表示计算出监视器2所显示的重叠图像的平均速度的时刻的位置的线标记b。线标记a及线标记b与重叠图像的再生的进展同时地进行移动。
此外,合成有线标记的移动信息数据有时也可以是时间变化图像。而且,匹配时相的处理如上所述既可以使用峰值来进行,也可以例如若与数据收集并行地收集ECG(electrocardiograph:心电图),则使用ECG的信息来进行。
此外,第2实施方式中,若设定大体包含发送连续波的所有扫描线在内的关注区域作为移动信息的检测对象,则也可以应用CWD法来进行所述处理。这种情况下,针对通过连续波收发而生成的一维的反射波数据,作为CDW处理部的PDW处理部14b取得按照扫描线的时间序列的一维的血流信息,图像生成部15a生成多普勒波形。并且,通过CWD法所生成的多普勒波形和移动信息数据一起显示。
接下来,使用图18说明第2实施方式涉及的超声波诊断装置的处理。图18是用于说明第2实施方式涉及的超声波诊断装置的处理的流程图。此外,在图18中,说明的是如下情况,即,生成移动信息数据后,开始进行多普勒显示用扫描,与多普勒显示用扫描并行,实时地执行多普勒波形和移动信息数据的并列显示。
如图18所示,第2实施方式涉及的超声波诊断装置的控制部17判定是否从操作者接受到在移动信息数据显示用的关注区域所包含的扫描线上对多普勒显示用的取样容积的设定(步骤S201)。在此,在未接受到取样容积的设定的情况下(步骤S201,否),控制部17待机至接受取样容积的设定为止。
另一方面,在接受到取样容积的设定的情况下(步骤S201,是),超声波探头1通过收发部11的控制,在包含取样容积的扫描线上进行超声波收发,收发部11利用超声波探头1接收到的反射波信号,生成扫描线上的反射波数据,并将其保存到帧缓冲器12内。并且,PWD处理部14b从帧缓冲器12中读出取样容积的反射波数据(步骤S202),执行HPF处理及FFT处理,从而取得取样容积的血流信息(步骤S203)。并且,图像生成部15a生成多普勒波形(步骤S204)。并且,通过控制部17的控制,监视器2使多普勒波形和移动信息数据的时相相匹配进行显示(步骤S205),然后结束处理。
此外,第2实施方式中,也可以是如下情况,即,收集多普勒波形后,进行多普勒波形和移动信息数据的并列显示。而且,在第2实施方式中,多普勒波形的收集也可以在收集血流图像之前进行。而且,在第2实施方式中,也可以设定多个移动信息数据显示用的ROI。这种情况下,多普勒显示用的扫描要针对多个ROI分别进行。而且,在第2实施方式中,移动信息数据显示用的ROI也可以是三维的。
如上所述,在第2实施方式中,与移动信息数据显示一起进行多普勒显示。而且,在PWD法中,可以将混叠速度设定得比CFM法高。因此,关于高流速的血流而言,相比CFM法,PWD法的FFT显示的可靠性更高。因此,如第2实施方式那样,将利用PWD法生成的多普勒波形,和高帧频的血流图像、利用该血流图像生成的移动信息数据并列显示,由此操作者可以确认通过血流图像、移动信息数据所观察的与血流动态相关的信息的可靠性。
(第3实施方式)
在第3实施方式中,使用图19及图20,来说明利用快速傅里叶变换法生成显示移动信息数据的情况。图19及图20是用于说明第3实施方式的图。
在第2实施方式中,多普勒波形和高帧频血流图像是分别进行收集的,用于生成高帧频血流图像的I/Q信号例如是被保持在帧缓冲器12内。因此,PWD处理部14b可以取得高帧频血流图像的截面上的任意位置的I/Q信号。
因此,在第3实施方式中,PWD处理部14b针对从帧缓冲器12取得的关注区域内的反射波数据,如使用图13说明的那样,例如通过快速傅里叶变换法进行频率解析,从而取得关注区域的按照时间序列的血流信息。即,如图19所示,PWD处理部14b取得高帧频血流图像上在任意位置设定的ROI内的按照时间序列的I/Q信号,并在帧方向上对取得的I/Q信号进行HPF处理然后进行FTT处理。由此,PWD处理部14b可以取得用于生成与以往多普勒显示中所使用的多普勒波形相同的多普勒波形的血流信息。换句话说,在第3实施方式中,将在第1实施方式中说明的ROI、或者在该ROI内设定的ROI作为取样容积。
控制部17在监视器2中显示基于关注区域的按照时间序列的血流信息的按照时间序列的数据。例如,在第3实施方式中,控制部17可以利用通过FFT处理取得的ROI内的血流信息,生成显示多普勒波形,或者生成显示时间变化曲线,或者生成显示重叠图像。
在此,PWD处理部14b可以在ROI内取得为了生成高帧频血流图像而设定的多个点各自的“速度、分散及能量”作为血流信息。或者,PWD处理部14b可以利用ROI内的多个点各自的“速度、分散及能量”,取得ROI内的“平均速度、平均分散及平均能量”作为平均血流信息。或者,PWD处理部14b通过将ROI分割成多个ROI,而可取得多个ROI的各自的平均血流信息。
因此,在第3实施方式中,作为多普勒波形,控制部17可以生成显示ROI的多普勒波形、ROI内的多个点各自的多普勒波形、或者ROI内设定的多个ROI各自的多普勒波形。而且,在第3实施方式中,作为时间变化曲线,控制部17可以生成显示ROI的时间变化曲线、ROI内的多个点各自的时间变化曲线、或者ROI内设定的多个ROI各自的时间变化曲线。而且,在第3实施方式中,作为重叠图像,控制部17可以生成显示以一个箭头表示ROI的平均血流信息的重叠图像、以多个箭头表示ROI内的多个点各自的血流信息的重叠图像、以一个箭头表示ROI内设定的多个ROI各自的血流信息的重叠图像、或者以多个箭头表示ROI内设定的多个ROI各自的血流信息的重叠图像。
此外,在第3实施方式中,也可以使PWD处理部14b和CFM处理部14a及检测部15b对同一反射波数据进行处理。即,如图20所示,控制部17也可以使得在CFM处理部14a及PWD处理部14b分别取得ROI内的血流信息。假设,若将CFM处理部14a的输出数据设为第1血流信息、将PWD处理部14b的输出数据设为第2血流信息,则控制部17显示时间变化曲线、时间变化图像、重叠图像等,作为基于第1血流信息的移动信息数据。而且,控制部17如上所述那样显示多普勒波形、时间变化曲线、重叠图像,作为基于第2血流信息的移动信息数据。
此外,第2血流信息的ROI可以是和第1血流信息的ROI相同的区域,也可以被设定在第1血流信息的ROI内。而且,第2血流信息的ROI也可以在第1血流信息的ROI内设定多个。而且,在第3实施方式中,也可以进一步进行FFT显示用扫描。而且,在第3实施方式中,若设定大体包含发送连续波的所有扫描线在内的关注区域作为移动信息的检测对象,则也可以应用CWD法进行所述的处理。
此外,在第3实施方式中,动画显示重叠图像的情况下,也可以在时间变化曲线等中显示表示同一时相的线标记。
接下来,使用图21,说明第3实施方式涉及的超声波诊断装置的处理。图21是用于说明第3实施方式涉及的超声波诊断装置的处理的流程图。此外,图21中说明的是如下情况,即,收集血流图像后进行ROI设定,并在收集血流图像后,进行所述PWD处理部14b及控制部17的处理。
如图21所示,第3实施方式涉及的超声波诊断装置的控制部17判定是否从操作者接受到对作为取样容积的关注区域的设定(步骤S301)。在此,在未接受到取样容积的设定的情况下(步骤S301,否),控制部17待机至接受到取样容积的设定为止。
另一方面,在接受到取样容积的设定的情况下(步骤S301,是),PWD处理部14b利用血流图像用所保持的反射波数据而取得取样容积的反射波数据(步骤S302),进行HPF处理及FFT处理,取得取样容积的血流信息(步骤S303)。然后,图像生成部15a生成多普勒波形(步骤S304)。然后,通过控制部17的控制,监视器2显示多普勒波形,或者将多普勒波形和移动信息数据的时相匹配地显示(步骤S305),然后结束处理。此外,第3实施方式也可以是如下情况,即,收集反射波数据之前进行ROI设定,在血流图像的拍摄过程中,实时显示基于利用FFT处理而从血流图像用所保持的反射波数据中取得的血流信息的多普勒波形等数据。
如上所述,在第3实施方式中,使用与高帧频血流像相同的信号进行FFT显示,因此与第2实施方式相比,能够以更准确的位置和时相观察血流动态。然而,与如第2实施方式那样利用PWD法等另行地收集反射波数据的情况相比,在第3实施方式中PRF不得不变低,所以混叠速度变低。即,从混叠速度的观点出发,若在心脏或颈动脉这样流速快的部位进行FFT显示,则第2实施方式中说明过的方法可靠性更高。因此,使用与高帧频血流像相同的信号而进行FFT显示的第3实施方式的方法,理想的是被用于观察速度比较慢的血流的情况中。在观察速度比较慢的血流的情况下,通过使用第3实施方式的方法,从而无须另外进行PDW用扫描,所以能够简化检查方法,从而缩短检查时间。
此外,所述第1实施方式~第3实施方式中说明的图像处理,也可以由与医用图像诊断装置独立地设置的图像处理装置来执行。具体来说,也可以是如下情况,即,在管理各种医用图像的数据的系统(system)即PACS(Picture Archiving and Communication Systems:影像归档和传输系统)的数据库(database)、或管理附有医用图像的电子病历(chart)的电子病历系统的数据库等中,存储帧缓冲器12所保持的数据,具有多普勒处理部14、图像处理部15、控制部17等功能的图像处理装置从这些数据库中接收反射波数据,进行所述图像处理。而且,在进行第1实施方式说明的图像处理的情况下,该图像处理装置例如也可以从数据库接收以60fps以上的帧频所生成的血流图像,并进行所述图像处理。
而且,图示的各装置的各构成要素为功能性概念,并非必须物理性地构成图示那样。即,各装置的分散·合并的具体方式并不限定于图示,其全部或者一部分可以根据各种负载、使用状况等,而以任意单位功能性或物理性地分散·合并而构成。进而,在各装置所进行的各处理功能的全部或任意一部分,可以由CPU及该CPU解析执行的程序来实现,或者能作为基于布线逻辑(wired logic)的硬件(hardware)来实现。
而且,第1实施方式~第3实施方式中说明的图像处理方法可以通过由个人计算机(personal computer)或工作站(workstation)等计算机,执行预先准备的图像处理程序而实现。该程序可以通过因特网(internet)等网络(network)进行发布。而且,该图像处理程序也可以被记录在硬盘(hard disk)、软盘(flexible disk)(FD)、CD-ROM、MO、DVD等计算机可读取的记录介质内,由计算机从记录介质读出后执行。
以上,如说明所示,根据第1实施方式~第3实施方式,可以准确且时间分辨率良好地观测二维或三维关注区域内的血流动态。
虽然说明了本发明的几个实施方式,但这些实施方式是作为例子而提示的,并不意图限定本发明的范围。这些实施方式能够以其他的各种方式来实施,在不脱离发明的要旨的范围内,能够进行各种的省略、置换、变更。这些实施方式或其变形与包含于发明的范围或要旨中一样,包含于权利要求书记载的发明及其均等的范围中。
Claims (14)
1.一种超声波诊断装置,其特征在于,包括:
处理部,从由超声波收发而收集到的二维或三维的反射波数据,取得扫描范围的按照时间序列的二维或三维的血流信息,所述超声波收发是在由多条扫描线形成的所述扫描范围中进行的;
图像生成部,根据所述扫描范围的按照时间序列的二维或三维的血流信息,生成按照时间序列的血流图像;
检测部,在所述按照时间序列的血流图像之间,检测预先设定的二维或三维的关注区域内的斑点的按照时间序列的移动信息;以及
控制部,使规定的显示部显示基于所述按照时间序列的移动信息的数据即移动信息数据。
2.根据权利要求1所述的超声波诊断装置,其特征在于:
所述处理部通过自相关法,取得所述扫描范围的按照时间序列的二维或三维的血流信息。
3.根据权利要求1所述的超声波诊断装置,其特征在于:
所述处理部从通过以各扫描线进行一次所述扫描范围内的超声波收发的扫描形态收集到的二维或三维的反射波数据,取得所述扫描范围的按照时间序列的二维或三维的血流信息。
4.根据权利要求1所述的超声波诊断装置,其特征在于:
所述检测部检测所述关注区域内的多个点各自的按照时间序列的移动信息,
所述控制部使所述图像生成部生成将所述多个点各自的按照时间序列的移动信息的代表值按照时间序列绘制而成的时间变化曲线、或者使所述时间变化曲线图像化而成的时间变化图像,作为所述移动信息数据。
5.根据权利要求4所述的超声波诊断装置,其特征在于:
所述控制部使所述图像生成部生成在所述血流图像的关注区域内重叠了表示所述多个点各自的移动信息的规定图形而成的重叠图像,作为所述移动信息数据。
6.根据权利要求5所述的超声波诊断装置,其特征在于:
所述控制部使所述图像生成部生成将所述代表值及直方图中的至少一个重叠到所述重叠图像的规定位置而成的图像,作为所述移动信息数据,所述直方图表示成为该代表值的计算依据的所述多个点各自的移动信息的值的分布。
7.根据权利要求6所述的超声波诊断装置,其特征在于:
在与所述时间变化曲线或所述时间变化图像并列地动画显示所述重叠图像的情况下,所述控制部使与所显示的重叠图像的时相对应的位置显示在所述时间变化曲线或所述时间变化图像中。
8.根据权利要求1所述的超声波诊断装置,其特征在于:
所述处理部针对通过以所述关注区域所包含的扫描线另行进行的超声波收发而收集到的反射波数据,进行频率解析,从而取得所述扫描线的按照时间序列的血流信息、或该扫描线上的取样容积中的按照时间序列的血流信息,
所述图像生成部利用通过所述频率解析取得的按照时间序列的血流信息,生成按照时间序列的多普勒波形,
所述控制部使所述多普勒波形和所述移动信息数据,以能够明确表示两者的数据的同一时相的方式,显示于所述规定的显示部。
9.根据权利要求8所述的超声波诊断装置,其特征在于:
所述检测部对所述关注区域内的多个点各自的按照时间序列的移动信息进行检测,
所述控制部使所述图像生成部生成在所述血流图像的关注区域内重叠了表示所述多个点各自的移动信息的规定图形而成的重叠图像,作为所述移动信息数据,进而在与所述多普勒波形并列地动画显示所述重叠图像的情况下,使与所显示的重叠图像的时相对应的位置显示在所述多普勒波形中。
10.根据权利要求1所述的超声波诊断装置,其特征在于,还包括:
存储部,该存储部保持所述二维或三维的反射波数据,
所述处理部针对从所述存储部取得的所述关注区域内的反射波数据进行频率解析,从而取得所述关注区域的按照时间序列的血流信息,
所述控制部进一步使基于通过所述频率解析取得的所述关注区域的按照时间序列的血流信息的数据,作为所述移动信息数据而显示于所述规定的显示部。
11.根据权利要求1所述的超声波诊断装置,其特征在于:
设定多个所述关注区域。
12.一种超声波诊断装置,其特征在于,包括:
存储部,保持通过超声波收发而收集到的二维或三维的反射波数据,该超声波收发是在由多条扫描线形成的扫描范围中进行的;
处理部,从所述存储部取得预先设定的二维或三维的关注区域内的反射波数据,对该取得的反射波数据进行频率解析,从而取得所述关注区域的按照时间序列的血流信息;以及
控制部,使基于所述关注区域的按照时间序列的血流信息的按照时间序列的数据显示于规定的显示部。
13.一种图像处理装置,其特征在于,包括:
检测部,在按照时间序列的血流图像之间,检测预先设定的二维或三维的关注区域内的斑点的按照时间序列的移动信息,该按照时间序列的血流图像是根据从二维或者三维的反射波数据取得的扫描范围的按照时间序列的二维或者三维的血流信息生成的,所述二维或者三维的反射波数据是通过在由多条扫描线形成的上述扫描范围中进行的超声波收发而收集到的;以及
控制部,使基于按照所述时间序列的移动信息的数据即移动信息数据显示于规定的显示部。
14.一种图像处理方法,其特征在于,包括:
检测部在按照时间序列的血流图像之间,检测预先设定的二维或三维的关注区域内的斑点的按照时间序列的移动信息,该按照时间序列的血流图像是根据从二维或者三维的反射波数据取得的扫描范围的按照时间序列的二维或者三维的血流信息生成的,所述二维或者三维的反射波数据是通过在由多条扫描线形成的所述扫描范围中进行的超声波收发而收集到的;以及
控制部使基于按照所述时间序列的移动信息的数据即移动信息数据显示于规定的显示部。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011268884A JP5972561B2 (ja) | 2011-12-08 | 2011-12-08 | 超音波診断装置、画像処理装置及び画像処理プログラム |
JP2011-268884 | 2011-12-08 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103156647A true CN103156647A (zh) | 2013-06-19 |
CN103156647B CN103156647B (zh) | 2016-03-16 |
Family
ID=48572636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210518548.6A Active CN103156647B (zh) | 2011-12-08 | 2012-12-06 | 超声波诊断装置、图像处理装置及图像处理方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9345446B2 (zh) |
JP (1) | JP5972561B2 (zh) |
CN (1) | CN103156647B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104188687A (zh) * | 2013-10-25 | 2014-12-10 | 中国科学院深圳先进技术研究院 | 基于超声回波射频信号的多普勒血流速度估测方法和系统 |
CN105310726A (zh) * | 2014-06-18 | 2016-02-10 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
WO2016172890A1 (zh) * | 2015-04-29 | 2016-11-03 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN106923864A (zh) * | 2015-11-03 | 2017-07-07 | 东芝医疗系统株式会社 | 超声波诊断装置、图像处理装置及图像处理程序 |
CN108523931A (zh) * | 2017-02-17 | 2018-09-14 | 通用电气公司 | 用于诊断医学成像的空间彩色流的方法和系统 |
CN110037740A (zh) * | 2018-01-15 | 2019-07-23 | 通用电气公司 | 用于超声流体成像的系统和方法 |
CN110403681A (zh) * | 2018-04-27 | 2019-11-05 | 佳能医疗系统株式会社 | 超声波诊断装置以及穿刺辅助方法 |
CN110477955A (zh) * | 2019-08-22 | 2019-11-22 | 电子科技大学 | 一种基于iq数据的血管自动识别方法 |
CN110507356A (zh) * | 2018-05-22 | 2019-11-29 | 佳能医疗系统株式会社 | 超声波诊断装置 |
CN112336375A (zh) * | 2019-08-08 | 2021-02-09 | 株式会社日立制作所 | 超声波诊断装置以及超声波图像处理方法 |
CN112512398A (zh) * | 2018-07-27 | 2021-03-16 | 富士胶片株式会社 | 医疗图像处理装置 |
CN113951930A (zh) * | 2021-09-16 | 2022-01-21 | 李世岩 | 一种三维颈部超声自动扫查与评估系统及方法 |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102081697B (zh) | 2009-11-27 | 2013-12-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种在超声成像空间中定义感兴趣容积的方法及其装置 |
JP6494903B2 (ja) * | 2013-02-14 | 2019-04-03 | ソニー株式会社 | 分析システム、分析プログラム及び分析方法 |
KR102243022B1 (ko) * | 2014-03-05 | 2021-04-21 | 삼성메디슨 주식회사 | 선택 정보에 기초하여 관심 영역에 포함된 혈류에 대한 정보를 출력하는 방법, 장치 및 시스템. |
EP2989986B1 (en) * | 2014-09-01 | 2019-12-18 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus and method of operating the same |
KR20160087221A (ko) * | 2015-01-13 | 2016-07-21 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 그 동작방법 |
CN107613878B (zh) * | 2015-03-30 | 2021-04-06 | 毛伊图像公司 | 用于检测物体运动的超声成像系统和方法 |
CN106333707B (zh) * | 2015-07-09 | 2020-12-01 | 深圳迈瑞生物医疗电子股份有限公司 | 超声多普勒图谱校正方法、装置及超声诊断系统 |
JP6674289B2 (ja) * | 2016-03-07 | 2020-04-01 | 国立大学法人旭川医科大学 | 蓄積一括事後処理による3次元体積ドプライメージング装置、3次元体積ドプライメージング方法、およびプログラム |
EP3616595B1 (en) * | 2016-05-23 | 2023-03-29 | Leica Instruments (Singapore) Pte. Ltd. | Medical observation device such as a microscope or an endoscope, and method for displaying medical images |
WO2017223563A1 (en) * | 2016-06-24 | 2017-12-28 | Duke University | Systems and methods for estimating cardiac strain and displacement using ultrasound |
EP3629936A1 (en) | 2017-05-25 | 2020-04-08 | Koninklijke Philips N.V. | Systems and methods for automatic detection and visualization of turbulent blood flow using vector flow data |
JP6829656B2 (ja) * | 2017-06-12 | 2021-02-10 | 株式会社日立製作所 | 超音波診断装置 |
KR20190085741A (ko) * | 2018-01-11 | 2019-07-19 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그의 제어 방법 |
JP7282564B2 (ja) * | 2018-03-30 | 2023-05-29 | キヤノンメディカルシステムズ株式会社 | 医用診断装置、医用画像処理装置、及び画像処理プログラム |
JP7215053B2 (ja) | 2018-10-02 | 2023-01-31 | コニカミノルタ株式会社 | 超音波画像評価装置、超音波画像評価方法および超音波画像評価プログラム |
US20200273181A1 (en) * | 2019-02-26 | 2020-08-27 | Canon Medical Systems Corporation | Analysis apparatus and ultrasound diagnosis apparatus |
KR102112699B1 (ko) * | 2019-06-12 | 2020-05-19 | 주식회사 뷰노 | 시계열적 신호 시각화 방법 및 이를 이용한 장치 |
US11760412B2 (en) * | 2021-10-07 | 2023-09-19 | Denso Corporation | Systems and methods for controlling a vehicle localization routine |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06114059A (ja) * | 1992-10-02 | 1994-04-26 | Toshiba Corp | 超音波カラードプラ断層装置 |
CN101352354A (zh) * | 2007-07-27 | 2009-01-28 | 株式会社东芝 | 超声波诊断装置 |
CN101563626A (zh) * | 2006-11-03 | 2009-10-21 | 皇家飞利浦电子股份有限公司 | 用于最佳散斑跟踪的双路径处理 |
CN101732036A (zh) * | 2008-11-05 | 2010-06-16 | 东芝医疗系统株式会社 | 医用图像解析装置以及图像解析用控制程序 |
US20110033098A1 (en) * | 2009-08-07 | 2011-02-10 | Medinol Ltd. | Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images |
CN102038522A (zh) * | 2009-10-21 | 2011-05-04 | 株式会社东芝 | 超声波诊断装置以及超声波数据取得方法 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5357867A (en) * | 1976-11-04 | 1978-05-25 | Hewlett Packard Yokogawa | Signal analyzer |
JPH0613031B2 (ja) | 1987-08-12 | 1994-02-23 | 株式会社東芝 | 超音波血流イメ−ジング装置 |
US5622174A (en) * | 1992-10-02 | 1997-04-22 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and image displaying system |
JP3724846B2 (ja) | 1995-06-15 | 2005-12-07 | 株式会社東芝 | 超音波診断装置 |
JPH11164831A (ja) * | 1997-12-03 | 1999-06-22 | Aloka Co Ltd | 超音波診断装置 |
US6048312A (en) * | 1998-04-23 | 2000-04-11 | Ishrak; Syed Omar | Method and apparatus for three-dimensional ultrasound imaging of biopsy needle |
JP4481386B2 (ja) * | 1999-06-07 | 2010-06-16 | 株式会社東芝 | 超音波診断装置 |
US6618493B1 (en) * | 1999-11-26 | 2003-09-09 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for visualization of motion in ultrasound flow imaging using packet data acquisition |
US6780152B2 (en) * | 2002-06-26 | 2004-08-24 | Acuson Corporation | Method and apparatus for ultrasound imaging of the heart |
JP4245428B2 (ja) * | 2003-07-29 | 2009-03-25 | 東芝医用システムエンジニアリング株式会社 | 超音波ドプラ診断装置及び画像データ生成方法 |
JP4801912B2 (ja) * | 2004-03-24 | 2011-10-26 | 株式会社東芝 | 超音波診断装置 |
JP4660126B2 (ja) * | 2004-06-18 | 2011-03-30 | 株式会社東芝 | 超音波血流イメージング装置 |
JP5570877B2 (ja) | 2010-06-04 | 2014-08-13 | 株式会社東芝 | 超音波診断装置 |
-
2011
- 2011-12-08 JP JP2011268884A patent/JP5972561B2/ja active Active
-
2012
- 2012-11-29 US US13/688,556 patent/US9345446B2/en active Active
- 2012-12-06 CN CN201210518548.6A patent/CN103156647B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06114059A (ja) * | 1992-10-02 | 1994-04-26 | Toshiba Corp | 超音波カラードプラ断層装置 |
CN101563626A (zh) * | 2006-11-03 | 2009-10-21 | 皇家飞利浦电子股份有限公司 | 用于最佳散斑跟踪的双路径处理 |
CN101352354A (zh) * | 2007-07-27 | 2009-01-28 | 株式会社东芝 | 超声波诊断装置 |
CN101732036A (zh) * | 2008-11-05 | 2010-06-16 | 东芝医疗系统株式会社 | 医用图像解析装置以及图像解析用控制程序 |
US20110033098A1 (en) * | 2009-08-07 | 2011-02-10 | Medinol Ltd. | Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images |
CN102038522A (zh) * | 2009-10-21 | 2011-05-04 | 株式会社东芝 | 超声波诊断装置以及超声波数据取得方法 |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104188687B (zh) * | 2013-10-25 | 2016-08-24 | 中国科学院深圳先进技术研究院 | 基于超声回波射频信号的多普勒血流速度估测方法和系统 |
CN104188687A (zh) * | 2013-10-25 | 2014-12-10 | 中国科学院深圳先进技术研究院 | 基于超声回波射频信号的多普勒血流速度估测方法和系统 |
CN105310726B (zh) * | 2014-06-18 | 2018-07-03 | 东芝医疗系统株式会社 | 超声波诊断装置、图像处理装置以及图像处理方法 |
CN105310726A (zh) * | 2014-06-18 | 2016-02-10 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
WO2016172890A1 (zh) * | 2015-04-29 | 2016-11-03 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN106102587A (zh) * | 2015-04-29 | 2016-11-09 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN106102587B (zh) * | 2015-04-29 | 2019-06-14 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN110013272A (zh) * | 2015-04-29 | 2019-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN110013273A (zh) * | 2015-04-29 | 2019-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN110013274A (zh) * | 2015-04-29 | 2019-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
US11944497B2 (en) | 2015-04-29 | 2024-04-02 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Ultrasonic blood flow imaging display method and ultrasonic imaging system |
CN110013273B (zh) * | 2015-04-29 | 2022-07-29 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像系统 |
CN106923864A (zh) * | 2015-11-03 | 2017-07-07 | 东芝医疗系统株式会社 | 超声波诊断装置、图像处理装置及图像处理程序 |
CN106923864B (zh) * | 2015-11-03 | 2021-01-08 | 东芝医疗系统株式会社 | 超声波诊断装置、图像处理装置及图像处理程序 |
CN108523931A (zh) * | 2017-02-17 | 2018-09-14 | 通用电气公司 | 用于诊断医学成像的空间彩色流的方法和系统 |
CN110037740A (zh) * | 2018-01-15 | 2019-07-23 | 通用电气公司 | 用于超声流体成像的系统和方法 |
CN110403681A (zh) * | 2018-04-27 | 2019-11-05 | 佳能医疗系统株式会社 | 超声波诊断装置以及穿刺辅助方法 |
CN110403681B (zh) * | 2018-04-27 | 2023-12-19 | 佳能医疗系统株式会社 | 超声波诊断装置以及图像显示方法 |
CN110507356A (zh) * | 2018-05-22 | 2019-11-29 | 佳能医疗系统株式会社 | 超声波诊断装置 |
CN112512398A (zh) * | 2018-07-27 | 2021-03-16 | 富士胶片株式会社 | 医疗图像处理装置 |
US11881190B2 (en) | 2018-07-27 | 2024-01-23 | Fujifilm Corporation | Medical image processing apparatus |
CN112512398B (zh) * | 2018-07-27 | 2024-03-29 | 富士胶片株式会社 | 医疗图像处理装置 |
CN112336375A (zh) * | 2019-08-08 | 2021-02-09 | 株式会社日立制作所 | 超声波诊断装置以及超声波图像处理方法 |
CN112336375B (zh) * | 2019-08-08 | 2024-04-12 | 富士胶片医疗健康株式会社 | 超声波诊断装置以及超声波图像处理方法 |
CN110477955B (zh) * | 2019-08-22 | 2021-05-11 | 电子科技大学 | 一种基于iq数据的血管自动识别方法 |
CN110477955A (zh) * | 2019-08-22 | 2019-11-22 | 电子科技大学 | 一种基于iq数据的血管自动识别方法 |
CN113951930A (zh) * | 2021-09-16 | 2022-01-21 | 李世岩 | 一种三维颈部超声自动扫查与评估系统及方法 |
CN113951930B (zh) * | 2021-09-16 | 2023-11-28 | 杭州影想未来科技有限公司 | 一种三维颈部超声自动扫查与评估系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5972561B2 (ja) | 2016-08-17 |
JP2013118984A (ja) | 2013-06-17 |
US9345446B2 (en) | 2016-05-24 |
US20130150717A1 (en) | 2013-06-13 |
CN103156647B (zh) | 2016-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103156647B (zh) | 超声波诊断装置、图像处理装置及图像处理方法 | |
CN103826541B (zh) | 超声波诊断装置以及控制方法 | |
JP4745133B2 (ja) | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム | |
JP3746115B2 (ja) | 超音波診断装置 | |
CN1883397B (zh) | 超声波诊断装置及超声波计测方法 | |
CN101361665B (zh) | 超声波诊断装置、超声波图像处理装置以及方法 | |
CN101292879B (zh) | 超声波诊断装置及其控制方法 | |
CN101797165B (zh) | 超声波诊断装置及超声波诊断装置的控制方法 | |
JPH08229038A (ja) | 3次元超音波画像作成方法および装置 | |
CN101601593B (zh) | 超声波诊断装置 | |
JP2007518512A (ja) | 心筋灌流を表示するための画像分割 | |
CN104640506A (zh) | 超声波诊断装置、医用图像处理装置和医用图像处理方法 | |
EP2610641B1 (en) | Ultrasound and system for forming a Doppler ultrasound image | |
CN100421627C (zh) | 超声波成像装置和超声波成像方法 | |
JP6113892B2 (ja) | 超音波診断装置、画像処理装置及び画像処理プログラム | |
CN102309339A (zh) | 超声波诊断装置、图像处理装置和图像处理方法 | |
US5123417A (en) | Apparatus and method for displaying ultrasonic data | |
JP4945277B2 (ja) | 超音波診断装置 | |
JP4381028B2 (ja) | 超音波診断装置 | |
US6425868B1 (en) | Ultrasonic imaging system | |
JP6356276B2 (ja) | 超音波診断装置及び制御方法 | |
JP2008253663A (ja) | 超音波診断装置およびその制御処理プログラム | |
JP4537754B2 (ja) | 超音波診断装置及びパルスドプラ計測装置 | |
JP2013099386A (ja) | 超音波診断装置および医用画像処理装置 | |
JP6444519B2 (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 | ||
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: 20160706 Address after: Japan Tochigi Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |