CN104640506B - 超声波诊断装置、医用图像处理装置和医用图像处理方法 - Google Patents
超声波诊断装置、医用图像处理装置和医用图像处理方法 Download PDFInfo
- Publication number
- CN104640506B CN104640506B CN201380001735.XA CN201380001735A CN104640506B CN 104640506 B CN104640506 B CN 104640506B CN 201380001735 A CN201380001735 A CN 201380001735A CN 104640506 B CN104640506 B CN 104640506B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- arrival time
- shearing wave
- displacement
- region
- 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.)
- Active
Links
Classifications
-
- 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/485—Diagnostic techniques involving measuring strain or elastic properties
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52019—Details of transmitters
- G01S7/5202—Details of transmitters for pulse systems
- G01S7/52022—Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
- G01S7/52095—Details related to the ultrasound signal acquisition, e.g. scan sequences using multiline receive beamforming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
根据本实施方式的超声波诊断装置(1),包括:超声波探头(3)、经由超声波探头(3)向第1区域发送在被检体内产生剪切波的第1超声波,向上述被检体内的第2区域发送第2超声波的发送部(11);基于第2超声波产生接收信号的接收部(13);用接收信号计算伴随着剪切波的向第2区域的传播的组织位移量的位移量计算部(19);基于与第2区域的各位置有关的位移量的时间变化确定剪切波到达各位置的到达时间的到达时间确定部(21);以及基于到达时间和根据到达时间预先设定的像素值产生根据到达时间分配了像素值的剪切波到达图像的图像产生部(23)。
Description
技术领域
本发明的实施方式涉及通过扫描被检体内而把被检体内部的脏器图像化并诊断被检体的疾病等的超声波诊断装置、医用图像处理装置和医用图像处理方法。
背景技术
在超声波诊断中,超声波诊断装置仅通过从体表触抵接声波探头这么简单的操作就可以实时地显示心脏的博动、胎儿的活动的样子。而且,超声波诊断装置因安全性高而反复地用于检查。而且,超声波诊断装置的系统的规模比X射线诊断装置、X射线计算机断层摄影(Computed Tomography:以下称为CT)装置、磁共振成像装置(Magnetic ResonanceImaging:以下称为MRI)等其它医用图像诊断装置小。因此,超声波诊断装置还具有通过向床侧移动而可以容易地在检查中使用等的简便性。另外,超声波诊断装置由于没有像使用了X射线的医用图像诊断装置等那样的对被检体的照射的影响,也被使用在产科、居家医疗等中。
近年来,具有评价被检体内的组织的硬度的功能的超声波诊断装置十分普及。作为硬度的评价方法大致分两种。第一种是,通过用超声波探头从体表按压和释放被检体内的组织,根据按压和释放时观测的断面内的各点的畸曲的大小把被检体内的组织的相对硬度可视化的方法。第二种是,通过从体表向组织赋予声辐射力、机械振动,随时间变化地观测断面内的各点中的组织的位移。具体地,是通过求出由声辐射力、机械振动产生的剪切波(shear wave)的传播速度,求出诊断对象组织的弹性模量的方法。
在第一种方法中,有局部的畸曲的大小依赖于操作者移动超声波探头的大小的特征。另外,在第一种方法中,操作者关注的区域(以下称为关注区域)的硬度评价通过与关注区域的周边的区域(以下称为周边区域)的硬度比较而进行。即,在第一种方法中,关注区域的硬度是与周边区域的硬度比较的相对硬度。
在第二种方法中,可以求出关注区域的绝对的弹性模量。具体地,第二种方法,一般地,时序地取得因产生的剪切波通过而造成的被检体内的各位置处的组织的位移。另外,有时除了剪切波造成的组织的位移以外在组织的位移中包含被检体的体动造成的组织整体的位移,所以一般执行消除组织整体的位移的处理。然后,推定各位置处的位移变成最大的时间作为剪切波到达的到达时间。基于推定的到达时间和从剪切波产生位置到各位置的距离,求出在被检体内传播的剪切波的传播速度。
但是,在组织整体的移动的方向和速度像短时间变化的那样的不规则时,不能正确地消除组织整体的位移。另外,由于关注区域的位置和组织的硬度,剪切波造成的组织的位移很微小,不能精度良好地推定到达时间。而且,在组织的边界面剪切波被反射、折射时,和在关注区域存在剪切波不传播的血管时,会有错误地求出到达时间的可能性。
在不能正确地消除组织整体的位移时,或剪切波被反射、折射时,或在关注区域存在有剪切波不传播的区域时,从剪切波的产生位置到关注区域的各位置的剪切波的到达时间成为从到达时间的整体的趋势偏离的点。因此,从整体的趋势偏离的点在剪切波的传播速度的计算过程中被排除,或被判断为算不出剪切波的传播速度,而以不被绘图的方式输出。因此,还提出了示出针对计算出的剪切波的传播速度的可靠性的程度的方法。
但是,有在剪切波的传播速度等的计测结果为未输出时,操作者不能具体地知道计测结果为未输出的理由的问题。而且,即使在剪切波的传播速度等的计测结果为输出了时,也存在操作者不能具体地知道计测结果是具有何种程度的可靠性的结果的问题。另外,计测结果低时,还存在如果改变与计测有关的区域,则操作者不能判断是否得到有可靠性的计测结果的问题。因此,操作者盲目地设定关心区域,根据求出的计测结果(数值)、触诊和其它临床所见综合地判断,判断求出的数值是否妥当。从以上可知,在现有的超声波诊断装置中,存在如果操作者不预先对被检体内的组织的硬度进行一定程度的分析,则不能判断上述方法获得的计测结果的妥当性的问题。
发明内容
(发明要解决的问题)
希望有即使操作者不能预测被检体内的组织的硬度时,也可以直观地、容易地把握计测结果的妥当性的显示方法。例如,操作者只要可以确认作为硬度计测的一次信息的剪切波在被检体内怎样传播(例如,像操作者期待的那样在方位(Lateral,横向)方向上均匀地传播,产生剪切波的反射和折射,或剪切波带有某种趋势地传播等)就可以。作为操作者确认剪切波的传播的样子的方法,有例如以动画显示剪切波的传播的样子的方法,但在临床中,每次计测时需要显示确认剪切波的传播的动画并确认剪切波的传播的时间,所以有不实用的问题。
本实施方式的目的在于,提供能够以检查者能够简便地确认的方式把被检体内产生的剪切波的传播可视化的超声波诊断装置。
(用来解决问题的方案)
根据本实施方式的超声波诊断装置,包括:具有多个振子的超声波探头;经由上述超声波探头,向第1区域发送在被检体内产生剪切波的第1超声波,向上述被检体内的第2区域发送用于观测预先确定的位置上的位移的第2超声波的发送部;基于上述第2超声波产生接收信号的接收部;用上述接收信号,计算伴随着上述剪切波向上述第2区域的传播的组织位移量的位移量计算部;基于与上述第2区域的各位置有关的上述位移量的时间变化,确定上述剪切波到达上述各位置的到达时间的到达时间确定部;以及基于上述到达时间和根据上述到达时间预先设定的像素值,产生根据上述到达时间分配上述像素值的剪切波到达图像的图像产生部。
(发明的效果)
如果采用本实施方式,则能够以检查者能够简便地确认的方式把被检体内产生的剪切波的传播可视化。
附图说明
图1是示出根据实施方式1的超声波诊断装置的构成的一例的构成图。
图2是示出根据实施方式1的、示出与位移量观测扫描功能有关的处理的步骤的流程图的一例的图。
图3是根据实施方式1的、把与产生剪切波的第1超声波有关的扫描线和与用来得到第1、第2接收信号的第2超声波有关的扫描线与关心区域一起示出的图。
图4是示出根据实施方式1的、示出与剪切波到达图像产生功能有关的处理的步骤的流程图的一例的图。
图5是示出根据实施方式1的、在有关心区域的位置组织位移量的时间变化的一例的图。
图6是示出根据实施方式1的、在有关心区域的位置除去了被检体的体动造成的组织位移量的位移量的时间变化的一例的图。
图7是示出根据实施方式1的、在有关心区域的位置组织位移量的时间变化的一例的图。
图8是示出根据实施方式1的、在有关心区域的位置除去了被检体的体动造成的组织位移量的位移量的时间变化的一例的图。
图9是根据实施方式1的、把产生的剪切波到达图像的一例与示出与剪切波到达时间对应的色相的图例的彩色映像表(color map)的一例一起示出的图。
图10是示出根据实施方式1的变形例的、示出与剪切波到达修正图像产生功能有关的处理的步骤的流程图的一例的图。
图11是根据实施方式1的变形例的、把剪切波到达修正图像的一例与剪切波到达图像的一例一起示出的图。
图12是示出根据实施方式2的超声波诊断装置的构成的一例的构成图。
图13是示出根据实施方式2的、示出剪切波传播速度图像产生处理的步骤的流程图的一例的图。
图14是示出根据实施方式2的第1变形例的、在位于4个不同的扫描线上的同一深度的4点处,远离超声波探头的方向的位移随时间的变化的一例的图。
图15是示出根据实施方式2的第1变形例的、在位于4个不同的扫描线上的同一深度的4点处,剪切波的到达时间与到剪切波的产生位置的距离的关系的一例的图。
(附图标记说明)
1:超声波诊断装置;3:超声波探头;5:装置本体;7:显示部;9:输入部;11:发送部;13:接收部;15:B模式处理部;17:多普勒处理部;19:位移量计算部;21:到达时间确定部;22:速度计算部;23:图像产生部;25:图像合成部;27:存储部;29:接口部;31:控制处理器(CPU);100:医用图像处理装置;101:医用图像处理装置
具体实施方式
以下,参照附图说明根据本实施方式的超声波诊断装置。另外,在以下的说明中,对具有基本相同的构成的构成要素赋予相同的附图标记,只在必要时进行重复说明。
(实施方式1)
图1是示出根据实施方式1的超声波诊断装置1的构成的构成图。像该图所示的那样,超声波诊断装置1具有:超声波探头3;装置本体5;显示部7;以及与装置本体5连接的用来把来自操作者的各种指示、命令、信息取入装置本体5的输入部9。而且,本超声波诊断装置1也可以经由后述的接口部29连接以心电仪、听诊器、脉波计、呼吸传感器为代表的未图示的生物体信号计测部和网络。
超声波探头3具有:多个压电振子、匹配层、和在多个压电振子的背面侧设置的背衬部件。多个压电振子是压电陶瓷等的声电可逆变换元件。多个压电振子并排安装在超声波探头3的前端。以下,说明一个压电振子构成一个信道的情况。压电振子响应从后述的发送部11供给的驱动信号而产生超声波。如果超声波经由超声波探头3发送到被检体P,则被发送的超声波(以下称为发送超声波)在被检体内的生物体组织中的声阻抗的不连续面处被反射。压电振子接收被反射的超声波,产生回波信号。回波信号的振幅取决于以与超声波的反射有关的不连续面为边界的声阻抗的差。另外,发送超声波在移动的血流和心脏壁等的表面处被反射了时的回波信号的频率因多普勒效应而与移动体(血流和心脏壁的表面)的超声波发送方向的速度分量相关地发生偏移。
以下,说明超声波探头3为用一维阵列进行二维扫描的探头的情况。另外,超声波探头3也可以是使一维阵列在与多个振子的排列方向正交的方向上摆动而执行三维扫描的机械四维探头。另外,超声波探头3不限于机械四维探头,也可以是二维阵列探头。
为了高效地对被检体P进行超声波的收发,在多个压电振子的超声波发射面侧设置匹配层。背衬部件防止超声波向压电振子的后方传播。
装置本体5具有:发送部11、接收部13、B模式处理部15、多普勒处理部17、位移量计算部19、到达时间确定部21、图像产生部23、图像合成部25、存储部27、接口部29和控制部(中央运算处理装置:Central Processing Unit,以下称为CPU)31。
发送部11具有:脉冲发生器11A、发送延迟电路11B和脉冲发生电路11C。脉冲发生器11A以预定的速度频率frHz(周期:1/fr秒)反复产生用来形成发送超声波的速度脉冲。产生的速度脉冲按信道数分配,送到发送延迟电路11B。
发送延迟电路11B针对多个信道的每一个使发送超声波收缩成束状,且对各速度脉冲赋予为了确定发送指向性所需的延迟时间(以下称为发送延迟时间)。发送超声波的发送方向或发送延迟时间(以下称为发送延迟模式)存储在后述的存储部27中。在后述的利用CPU3发送超声波时参照在存储部27中存储的发送延迟模式。
脉冲发生电路11C在基于该速度脉冲的定时,针对超声波探头3的每一个振子施加电压脉冲(驱动信号)。由此,超声波束被发送给被检体。
以下,详述产生剪切波的第1超声波的发送。发送部11经由超声波探头3向后述的第1区域发送在被检体内产生剪切波的推式脉冲(第1超声波)。发送部11在发送第1超声波的前后,向后述的第2区域发送第2超声波。为了得到在第1区域产生的剪切波造成的组织位移量,利用发送部11发送第2超声波。具体地,发送部11在经由后述的输入部9输入的关心区域(Region Of Interest:以下称为ROI)或操作者关注的区域(以下称为关注区域)等的预定区域的方位方向的边界附近(以下称为第1区域),以使焦点(以下称为剪切波产生位置)为预先设定的深度的方式,发送推式脉冲。换言之,剪切波产生位置设置在例如在预定区域外且在预定区域附近的扫描线上。预先设定的深度为例如预定区域的平均深度。
推式脉冲的频率与在B模式和多普勒模式中发送的超声波的频率基本相同。推式脉冲的波数比在B模式和多普勒模式中发送的超声波的波数多。具体地,发送部11利用后述的CPU31从存储部27读出例如用来在第1区域的深度方向的中点使超声波收缩成束状的发送延迟模式(以下称为剪切波产生发送延迟模式)。发送延迟电路11B把剪切波产生发送延迟模式造成的发送延迟时间赋予各速度脉冲。脉冲发生电路11C在基于该速度脉冲的定时针对每一个振子施加用来发生推式脉冲的驱动信号。
发送部11在发送第1超声波的前后,对预定区域(第2区域)发送第2超声波。
接收部13具有:前置放大器13A、未图示的模拟数字(Analog to digital,以下称为A/D)变换器、接收延迟电路13B和加法器13C。接收部13基于第2超声波和剪切波造成的反射波产生接收信号。
前置放大器13A针对每个信道把经由超声波探头3取入的来自被检体P的回波信号放大。A/D变换器把放大了的接收回波信号变换成数字信号。
接收延迟电路13B对变换成数字信号的接收回波信号赋予为了确定接收指向性所需的延迟时间(以下称为接收延迟时间)。回波信号的接收方向或接收延迟时间(以下称为接收延迟模式)被存储到后述的存储部27中。由后述的CPU31参照在存储部27中存储的接收延迟模式。加法器13C把被赋予了延迟时间的多个回波信号相加。通过该相加,接收部13产生强调了来自与接收指向性对应的方向的反射分量的接收信号(也称为RF(射频)信号)。利用该发送指向性和接收指向性确定超声波收发的综合的指向性。利用该综合的指向性确定超声波束(所谓的“超声波扫描线”)。具体地,接收部13接收与剪切波造成的物理的组织位移有关的第2超声波的反射波。
B模式处理部15具有未图示的包络线检波器、对数变换器等。包络线检波器对从接收部13输出的接收信号执行包络线检波。包络线检波器把进行了包络线检波的信号输出到后述的对数变换器。对数变换器对进行了包络线检波的信号进行对数变换,相对地强调弱的信号。B模式处理部15基于被对数变换器强调了的信号,产生各扫描线和各超声波收发中的每个深度的信号值(B模式数据)。
另外,在超声波探头3为机械四维探头时、为二维阵列探头时,B模式处理部15也可以产生三维B模式数据,该三维B模式数据由与被扫描区域中的方位(Azimuth)方向、高度(Elevation)方向、深度方向(以下称为距离(Range)方向)分别对应地排列的多个信号值构成。距离方向是扫描线上的深度方向。方位方向是例如沿一维超声波振子的排列方向的电子扫描方向。高度方向是一维超声波振子的机械的摆动方向。另外,三维B模式数据也可以是使多个像素值或多个亮度值等沿扫描线与方位方向、高度方向、距离方向分别对应地排列得到的数据。另外,三维B模式数据也可以是与在被扫描区域中预先设定的ROI有关的数据。另外,B模式处理部15也可以产生体数据来代替三维B模式数据。以下,把由B模式处理部15产生的数据总称为B模式数据。
多普勒处理部17具有:未图示的混相器、低通滤波器(Low Pass Filter:以下称为LPF)、速度/色散/功率运算器件等。混相器使具有与发送频率相同的频率f0的基准信号与从接收部17输出的接收信号相乘。通过该相乘,得到多普勒偏移频率fd的分量的信号和具有(2f0+fd)的频率分量的信号。LPF去除来自混相器的具有两种频率分量的信号中的、高频分量(2f0+fd)的信号。多普勒处理单元142通过去除高频分量(2f0+fd)的信号,产生具有多普勒偏移频率fd的分量的多普勒信号。
另外,多普勒处理部17也可以使用正交检波方式产生多普勒信号。此时,接收信号(RF信号)被正交检波而变换成IQ信号。多普勒处理单元142通过对IQ信号进行复数傅里叶变换而产生具有多普勒偏移频率fd的分量的多普勒信号。多普勒信号是例如血流、组织、造影剂造成的多普勒分量。
速度/色散/功率运算器件具有未图示的MTI(Moving Target Indicator,活动目标显示)滤波器、LPF滤波器、自相关运算器等。另外,也可以取代自相关运算器而具有互相关运算器。MTI滤波器对产生的多普勒信号除去脏器的呼吸性移动、搏动性移动等造成的多普勒分量(杂波分量)。为了从多普勒信号抽出与血流有关的多普勒分量(以下称为血流多普勒分量)而使用MTI滤波器。为了从多普勒信号抽出与组织的移动有关的多普勒分量(以下称为组织多普勒分量)而使用LPF。
自相关运算器对血流多普勒分量和组织多普勒分量算出自相关值。自相关运算器基于算出的自相关值算出血流和组织的平均速度值、色散值、多普勒信号的反射强度(功率)等。速度/色散/功率运算器件基于基于多个多普勒信号的血流和组织的平均速度值、色散值、多普勒信号的反射强度等,产生预定区域的各位置处的彩色多普勒数据。以下,把多普勒信号和彩色多普勒数据总称为多普勒数据。
以下,说明利用自相关运算器计算组织多普勒分量的相位差的处理(以下称为自相关处理)。首先,对一个扫描线方向,以预定的脉冲重复频率(Pulse RepetitionFrequency:以下称为PRF)多次(n次)收发超声波。然后,利用接收部13产生多个(n个)接收信号。从与多个接收信号分别对应的多个多普勒信号抽出多个组织多普勒分量。对组织多普勒分量附加表示一个扫描线方向的深度和收发次数的编号作为附带信息。
自相关运算器对一个扫描线方向的相同深度处的多个(n个)组织多普勒分量执行自相关处理。具体地,基于第i(1≤i≤n-1)个组织多普勒分量TD(i)和第(i+1)(2≤i+1≤n)个组织多普勒分量TD(i+1),计算组织多普勒分量的相位差。更详细地,自相关运算器通过对TD(i)与TD(i+1)的复数共轭积TD(i)×TD*(i+1)执行复数傅里叶变换,计算TD(i)与TD(i+1)的相位差。该计算相位差的步骤从i=1重复计算到i=n-1。即,用因一个扫描线方向上的n次超声波的收发产生的n个接收信号,针对每一个1/PRF和每一个深度计算(n-1)个相位差。另外,也可以在每次产生组织多普勒分量时进行相位差的计算。相位差与每一个1/PRF和每一个深度的组织位移量对应。通过累积从第1个到第i个相位差,可以求出以第1个组织多普勒分量作为基准的第i个组织位移量。
另外,也可以取代自相关处理而利用由未图示的互相关运算器进行的互相关处理进行相位差的计算。互相关处理是例如计算第1次超声波的收发造成的组织多普勒分量与第2次以后的超声波的收发造成的组织多普勒分量的相位差的处理。即,互相关处理是以第1次超声波的收发造成的组织多普勒分量为基准来计算相位差的处理。
互相关运算器对一个扫描线方向的相同深度处的多个(n个)组织多普勒分量执行互相关处理。具体地,基于第1个组织多普勒分量TD(1)和第j(2≤j≤n)个的组织多普勒分量TD(j)计算组织多普勒分量的相位差。该计算相位差的步骤从j=2重复计算到j=n。即,针对一个扫描线方向上的n次超声波的收发,针对每一个1/PRF和每一个深度计算(n-1)个相位差。互相关运算器的相位差示出以第1个组织多普勒分量为基准的第j个组织多普勒分量的相位。即,由互相关运算器输出的相位差示出以第1个组织多普勒分量为基准的组织位移量。
位移量计算部19基于推式脉冲发送之前由自相关器得到的组织位移量(以下称为第1位移数据)、和推式脉冲发送之后由自相关运算器得到的组织位移量(以下称为第2位移数据)中的第n次超声波收发即将结束时的组织位移量(第n次、第(n-1)次等),针对第2区域中的每一个深度把与被检体的体动有关的组织位移量(以下称为体动位移量)近似化。位移量计算部19通过从第1、第2位移数据减去近似得到的体动位移量,计算伴随着剪切波的传播的组织位移量(以下称为剪切波传播数据)。另外,剪切波传播数据也可以是沿扫描线方向远离超声波探头3的方向上的组织位移量。
以下,为了具体地说明,在推式脉冲发送之前对预定区域内的一个扫描线执行的超声波的收发次数(以下称为推式脉冲前收发次数)为15次。另外,在推式脉冲发送之后对预定区域内的一个扫描线执行的超声波的收发次数(以下称为推式脉冲后收发次数)n为65次。
具体地,位移量计算部19基于与推式脉冲前收发次数(15次)有关的第1位移数据、和分别与推式脉冲后收发次数(65次)中的第64次和第65次对应的多个位移量,在与第1、第2位移数据有关的超声波的整个收发期间,用例如2次多项式对体动位移量的时间变化进行近似。以下,把表示体动位移量的时间变化的数据称为体动位移数据。针对预定区域内的各位置(即,预定区域内的每一个扫描线和每一个深度)产生体动位移数据。位移量计算部19通过从第1、第2位移数据减去体动位移数据,计算示出剪切波造成的组织位移量的时间变化的剪切波传播数据。针对预定区域内的各位置产生剪切波传播数据。位移量计算部19把剪切波传播数据输出到后述的到达时间确定部21。位移量计算部19通过以上的处理计算与预定区域内的各位置对应的剪切波传播数据。
到达时间确定部21基于剪切波传播数据中的组织位移量的时间变化,确定到达预定区域内的各位置的剪切波的到达时间(以下称为剪切波到达时间)。具体地,到达时间确定部21针对预定区域内的各位置,以推式脉冲的发送开始时或推式脉冲的发送结束时为时刻0,确定与剪切波传播数据中的组织位移量的最大值对应的时间作为剪切波到达时间。到达时间确定部21把与预定区域内的各位置处的到达时间有关的数据(以下称为到达时间数据)输出到后述的图像产生部23。
另外,到达时间确定部21也可以确定与剪切波传播数据中的组织位移量的时间变化中的最大时间变化对应的时间,作为剪切波到达时间。另外,到达时间确定部21也可以把以推式脉冲的发送开始时或推式脉冲的发送结束时为基准而偏离了预定时间的时刻作为0,确定与剪切波传播数据中的组织位移量的最大值对应的时间作为剪切波到达时间。
另外,到达时间确定部21也可以基于剪切波传播数据中的组织位移量的时间变化中的极大值中的最大位移量确定到达时间。具体地,到达时间确定部21抽出剪切波传播数据中的组织位移量的时间变化中的、取极大值的数据。例如,在剪切波传播数据的波形具有多个峰时,也存在多个取极大值的点。到达时间确定部21抽出这些极大值中的、具有最大位移量的极大值。到达时间确定部21基于具有该极大值的位移量确定到达时间。
图像产生部23具有:未图示的数字扫描转换器(Digital Scan Converter:以下称为DSC)和图像存储器。图像产生部23对DSC执行坐标变换处理(再次取样)。坐标变换处理是例如把由B模式数据、多普勒数据和到达时间数据构成的超声波扫描的扫描线信号列变换成以电视等为代表的一般的视频格式的扫描线信号列的处理。图像产生部23通过坐标变换处理生成作为显示图像的超声波图像。具体地,图像产生部23基于B模式数据产生B模式图像。图像产生部23基于多普勒数据产生平均速度图像、色散图像、功率图像等的多普勒图像。
图像产生部23基于到达时间数据和在后述的存储部27中存储的到达时间色相对应表,产生在预定区域的各位置处根据到达时间分配了色相的剪切波到达图像。到达时间色相对应表是例如色相与到达时间的值的对应表。例如,到达时间为0时的色相是蓝色。色相例如伴随着到达时间的增大而以从蓝到蓝绿、绿、黄绿、黄、橙、红的顺序变化,最大到达时间被定义为红色。另外,在实施方式中,作为例子示出了分配色相的例子,但分配的参数不限于色相,只要是用来定义像素值的参数,什么都可以。例如,也可以在生成灰阶的剪切波到达图像时,到达时间色相对应表中没有色相而分配黑白的亮度值,产生在预定区域的各位置处根据到达时间分配了亮度值的灰阶的剪切波到达图像。
图像存储器存储与产生的超声波图像(B模式图像、平均速度图像、色散图像、功率图像、剪切波到达图像)对应的数据(以下称为图像数据)。利用经由后述的输入部9进行的操作者的指示读出在图像存储器中存储的图像数据。图像存储器是例如,保存与刚要定格(freeze)之前的多个帧对应的超声波图像的存储器。通过连续显示(电影显示)在该电影存储器(cine-memory)中存储的图像,也可以在后述的显示部7上显示超声波活动图像。
图像合成部25把各种参数的文字信息和刻度等合成到超声波图像。图像合成部25把合成后的超声波图像输出到后述的显示部7。图像合成部25产生将剪切波到达图像在B模式图像上定位并重叠的重叠图像。图像合成部23把产生的重叠图像输出到显示部7。
存储部27存储:聚焦深度不同的多个接收延迟模式、多个发送延迟模式、多个剪切波产生发送延迟模式、本超声波诊断装置1的控制程序、诊断规程、收发条件等的各种数据群、诊断信息(患者ID、医师观察等)、由接收部13产生的接收信号、由B模式处理部15产生的B模式数据、由多普勒处理部17产生的多普勒数据、预定区域内的各位置处的第1、第2位移数据、体动位移数据、剪切波传播数据、到达时间数据、到达时间色相对应表、B模式图像、平均速度图像、色散图像、功率图像、剪切波到达图像、与剪切波到达图像的产生有关的算法(以下称为剪切波到达图像产生算法)等。另外,上述的图像存储器也可以设置在存储部27中。
接口部29是与输入部9、网络、未图示的外部存储装置和生物体信号计测部有关的接口。由装置本体5得到的超声波图像等的数据和解析结果等可以经由接口部29和网络转送到其它装置。另外,接口部29也可以经由网络下载由未图示的其它医用图像诊断装置取得的与被检体有关的医用图像。
CPU31基于由操作者经由输入部9输入的对B模式、多普勒模式和剪切波到达图像显示模式的选择、帧频、被扫描深度、发送开始和结束,读出在存储部27中存储的发送延迟模式、接收延迟模式、剪切波产生发送延迟模式和装置控制程序,根据它们控制装置本体5。例如,CPU31从存储部27读出剪切波到达图像产生算法。CPU31根据读出的剪切波到达图像产生算法,控制发送部11、多普勒处理部17、位移量计算部19和到达时间确定部21。剪切波到达图像显示模式是产生剪切波到达图像,执行显示的模式。
具体地,如果经由输入部7被输入剪切波到达图像显示模式,则为了对预定区域内的扫描线方向发送用来得到第1位移数据的超声波,CPU31控制发送部11。为了产生第1位移数据,CPU31控制多普勒处理部17。
然后,为了产生推式脉冲,CPU31把从存储部27读出的剪切波产生发送延迟模式输出到发送部11。为了对预定区域内的扫描线方向发送用来得到第2位移数据的超声波,CPU31控制发送部11。为了产生第2位移数据,CPU31控制多普勒处理部17。为了基于第1、第2位移数据产生体动位移数据,CPU31控制位移量计算部19。为了基于第1、第2位移数据和体动位移数据产生剪切波传播数据,CPU31控制位移量计算部19。为了基于剪切波传播数据产生到达时间数据,CPU31控制到达时间确定部21。为了基于从存储部27读出的到达时间色相对应表和到达时间数据产生剪切波到达图像,CPU31控制图像产生部23。
显示部7基于来自图像合成部25的输出,显示B模式图像和多普勒图像等的超声波图像、剪切波到达图像、重叠图像等。另外,显示部7也可以对显示的图像执行亮度、对比度、动态容限、γ修正等的调整和彩色映像表的分配。
输入部9与接口部29连接,把来自操作者的各种指示、命令、信息、选择、设定取入装置本体5。输入部9具有未图示的轨迹球、开关按钮、鼠标、键盘等的输入器件。输入器件检测在显示画面上显示的光标的坐标,把检测到的坐标输出到后述的CPU31。另外,输入器件也可以是以覆盖显示画面的方式设置的触摸命令屏。此时,输入部9以电磁感应式、电磁失真式、感压式等的坐标读取原理检测触摸指示的坐标,把检测到的坐标输出到CPU31。另外,如果操作者操作输入部9的结束按钮或定格按钮,则超声波的收发结束,装置本体5成为暂时停止状态。
输入部9通过操作者的指示输入关心区域或操作者关注的预定区域。输入部9输入利用操作者的指示选择的模式。
(位移量观测扫描功能)
位移量观测扫描功能是与用来取得伴随着剪切波的传播的组织位移量的扫描步骤有关的功能。以下,说明与位移量观测扫描功能有关的处理(以下称为位移量观测扫描处理)。
图2示出位移量观测扫描处理的步骤的流程图。图3是在位移量观测扫描处理中,把产生剪切波的扫描线(以下称为剪切波产生扫描线)53和与用来观测剪切波造成的组织位移量的超声波的收发有关的扫描线(以下称为观测扫描线)54、55、56、57,与关心区域(或关注区域)52一起示意地示出的示意图。
在向被检体P发送超声波之前,利用操作者经由输入部9进行的指示设定关心区域52(步骤Sa1)。此时,示出观测扫描线的编号的变量k被初始化。另外,也可以取代关心区域52而设置操作者关注的区域。对关心区域52内的第k个扫描线(例如54)多次(例如15次)收发第2超声波,每次收发超声波时产生第1接收信号(步骤Sa2)。然后,基于关心区域52的位置,从存储部27读出剪切波产生发送延迟模式。用读出了的剪切波产生发送延迟模式对关心区域52附近的预定的扫描线54发送推式脉冲(第1超声波)(步骤Sa3)。发送推式脉冲之后,对第k个扫描线(例如54)多次(例如65次)收发超声波(第2超声波),每次收发超声波时产生第2接收信号(步骤Sa4)。在整个关心区域即扫描线55到57中,只要没有产生第1、第2接收信号(步骤Sa5),k就加一(步骤Sa6)。然后,重复步骤Sa2到步骤Sa4的处理。
位移量观测扫描处理是,对关心区域52内的多个观测扫描线的每一个,重复进行如下处理:在发送推式脉冲之前为了产生第1接收信号而收发超声波,然后,发送推式脉冲,在推式脉冲发送之后为了产生第2接收信号而进行超声波收发。即,位移量观测扫描处理是用一个推式脉冲执行针对一个扫描线的超声波收发的扫描处理。另外,也可以在发送一个推式脉冲之后,通过多个扫描线执行超声波收发。另外,也可以执行所谓的同时并列接收,即在发送一个推式脉冲之后发送1次超声波,执行针对多个扫描线方向的接收。
(剪切波到达图像产生功能)
剪切波到达图像产生功能是产生用色相示出预定区域的各位置处的剪切波到达时间的剪切波到达图像的功能。以下,说明根据剪切波到达图像产生功能的处理(以下称为剪切波到达图像产生处理)。
图4是示出剪切波到达图像产生处理的步骤的一例的流程图。
通过针对一个扫描线的多次超声波的收发,用推式脉冲发送之前产生的多个第1接收信号执行自相关处理。利用自相关处理产生与预定区域内的各位置有关的第1位移数据(步骤Sb1)。通过针对一个扫描线的多次超声波的收发,用推式脉冲发送后产生的多个第2接收信号执行自相关处理。通过自相关处理产生与预定区域内的各位置有关的第2位移数据(步骤Sb2)。
图5是示出在有关心区域的位置组织位移量的时间变化的一例的图。即,图5是汇总与有关心区域的位置有关的第1、第2位移数据得到的图。图5的横轴与除去了推式脉冲的超声波的发送次数对应。图5的纵轴表示组织位移量。针对从1到15的发送次数的组织位移量与第1位移数据对应。与从16到80的发送次数对应的组织位移量与第2位移数据对应。图5中的位移量是以针对第1次超声波的发送次数的组织位移量作为基准的位移量。即,图5中的位移量与由自相关处理产生的相位差的累积值、或由互相关处理产生的相位差对应。
图5示出组织因被检体的体动而在一个方向(靠近超声波探头3的方向)上大致等速地移动的情况。因被检体的体动,不能识别剪切波的传播造成的组织的位移。
基于第1位移数据和第2位移数据中的第n次超声波收发即将结束时的组织位移量(第n次、第(n-1)次等),产生示出被检体的体动位移量的体动位移数据(步骤Sb3)。通过从第1、第2位移数据减去体动位移数据,产生示出剪切波造成的组织位移量的时间变化的剪切波传播数据(步骤Sb4)。基于剪切波传播数据,确定从剪切波的产生位置到达与剪切波传播数据对应的关心区域内的各位置的剪切波到达时间(步骤Sb5)。
图6是示出从图5所示的组织位移量除去了被检体的体动造成的位移量的位移量的时间变化(剪切波传播数据)的一例的图。图6的横轴与除去了推式脉冲的超声波的发送次数对应。图6的纵轴表示组织位移量。图6的纵轴的刻度和符号与图5的纵轴不同。即,图6的纵轴表示远离超声波探头的方向,是与图5中的纵轴的符号相反的负。另外,剪切波造成的组织位移量的位数与体动造成的组织位移量的位数不同。
像图6所示的那样,在剪切波通过后(第22次的发送次数附近),组织位移量的绝对值最大。到达时间确定部21在图6中确定组织位移量的绝对值为最大的图中的48的时间(发送次数)作为剪切波到达时间。另外,在确定组织位移量的时间变化为负的最大的时间作为剪切波到达时间时,到达时间确定部21把发送次数16的附近确定为剪切波到达时间。
图7是示出在有与图5不同的关心区域的位置,组织位移量的时间变化的一例的图。图7的纵轴示出远离超声波探头的方向,是与图5中的纵轴的符号相反的负。图7示出组织因被检体的体动而在一个方向(远离超声波探头3的方向)上大致等速地移动的情况。因被检体的体动,不能识别剪切波的传播造成的组织的位移。
图8是示出从图7所示的组织位移量除去了被检体的体动造成的位移量的位移量的时间变化(剪切波传播数据)的一例的图。图8的纵轴和横轴与图6的纵轴和横轴分别相同。图8示出在得到剪切波传播数据的位置处,剪切波造成的组织位移量过小时或被检体的体动造成的组织位移量不规则时产生的剪切波传播数据的一例。到达时间确定部21确定发送次数24的附近作为剪切波到达时间。图8的曲线图与图6的曲线图不同,剪切波到达时间不明显。在得到图8那样的剪切波传播数据时,在剪切波到达图像中是与周围的色相相比例如大大不同的色相。
基于确定的剪切波到达时间和到达时间色相对应表,产生根据剪切波到达时间分配了色相的剪切波到达图像(步骤Sb6)。产生的剪切波到达图像与B模式图像重叠地显示在显示部7上。
图9把产生的剪切波到达图像的一例与示出与剪切波到达时间对应的色相的图例的彩色映像表的一例一起示出的图。图9的(a)、(b)是示出在B模式图像上重叠剪切波到达图像得到的重叠图像的一例的图。图9的(a)、(b)中的圆表示关注区域(目标)。图9的(c)是示出图9的(a)、(b)的彩色映像表的图例的图。
图9的(a)中的关注区域表示具有比关注区域的周围更柔软的组织的区域。图9的(b)中的关注区域表示具有比关注区域的周围更硬的组织的区域。图9的(a)、(b)都示出方位方向的色相的变化(由于图是淡色,所以是浓淡的变化)为面向图面时剪切波从左到右传播的样子。在图9的(a)中,在关注区域的下方区域显示斑点模样(马赛克模样)。即,在图9的(a)中,关注区域的下方区域示出例如像图8所示的那样,不能正确评价剪切波到达时间的区域。如果比较图9的(a)和图9的(b),则在图9的(a)的中央附近沿方位方向,色相(图中是浓淡)比图9的(b)更急剧地变化。因此,推定为在图9的(a)的中央附近,剪切波的传播速度慢,即,是比图9的(a)的关注区域周边的组织更柔软的组织。
在图9的(b)的中央附近,沿方位方向,色相(图中是浓淡)比图9的(a)更缓慢地变化。因此,推定为在图9的(b)的中央附近,剪切波的传播速度快,即,是比图9的(b)的关注区域周边的组织更硬的组织。
(变形例)
与实施方式1的不同之处在于,以利用操作者的指示在预定区域内指定的指定位置处的剪切波到达时间作为基准,确定对预定区域内的各位置处的剪切波到达时间进行修正后得到的修正到达时间,基于修正到达时间和到达时间色相对应表,产生对剪切波到达图像的色相进行修正后得到的剪切波到达修正图像。本变形例具有对用实施方式1产生的剪切波到达图像追加执行的处理。
输入部9把由操作者指定的指定位置输入装置本体5。
到达时间确定部21以指定位置处的剪切波到达时间作为基准,确定对预定区域内的各位置处的剪切波到达时间进行修正后得到的修正到达时间。具体地,到达时间确定部21通过从各位置处的剪切波到达时间减去指定位置处的剪切波到达时间,确定各位置处的修正到达时间。到达时间确定部21把确定了的修正到达时间输出到图像产生部23。
图像产生部23基于修正到达时间和到达时间色相对应表,产生对剪切波到达图像中的色相进行修正后得到的剪切波到达修正图像。具体地,图像产生部23用到达时间色相对应表确定与各位置处的修正到达时间的值对应的色相。图像产生部23通过把剪切波到达图像的各位置处的色相变更成与修正到达时间对应的色相,产生剪切波到达修正图像。
图像合成部25通过把在B模式图像上重叠的剪切波到达图像更新成剪切波到达修正图像,产生对重叠图像的色相进行修正后得到的修正重叠图像。
显示部7把显示的剪切波到达图像置换成剪切波到达修正图像进行显示。另外,显示部7在显示重叠图像时,也可以把显示的重叠图像置换成修正重叠图像进行显示。
(剪切波到达修正图像产生功能)
剪切波到达修正图像产生功能是以指定位置处的剪切波到达时间作为基准确定对预定区域内的各位置处的剪切波到达时间进行修正后得到的修正到达时间,基于修正到达时间和到达时间色相对应表产生对剪切波到达图像的色相进行修正后得到的剪切波到达修正图像的功能。以下,说明根据剪切波到达修正图像产生功能的处理(以下称为剪切波到达修正图像产生处理)。
图10是示出剪切波到达修正图像产生处理的步骤的一例的流程图。
在显示部7上显示的剪切波到达图像的关心区域内,利用操作者的指示输入指定位置(或指定区域)(步骤Sc1)。另外,指定位置可以是点,也可以是线段。作为指定位置的线段,是例如通过利用操作者的指示在关心区域中指定的位置与扫描线平行的线段。以指定位置处的剪切波到达时间作为基准,确定在各位置对剪切波到达时间进行修正后得到的修正到达时间(步骤Sc2)。从存储部27读出到达时间色相对应表。基于读出的到达时间色相对应表和修正到达时间,产生对剪切波到达图像中的色相进行修正后得到的剪切波到达修正图像(步骤Sc3)。
图11是把剪切波到达修正图像与剪切波到达图像一起示出的图。图11(a)是示出由实施方式1产生的剪切波到达图像的图。图11(b)示出由本变形例产生的剪切波到达修正图像的图。在图11中,面向图面时剪切波从左方向向右方向传播。图11中的圆形的区域91表示比周边组织剪切速度快的区域。图11中的圆形的区域92表示比周边组织剪切速度慢的区域。
在图11(a)中,在例如关注圆形的区域92地显示剪切波到达修正图像时,经由输入部9输入作为剪切波到达时间的基准的线段90(即,剪切波到达时间为0的位置)。到达时间确定部21从ROI内的同一深度的位置处的剪切波到达时间减去线段90(以下称为基准线)上的同一深度的位置处的剪切波到达时间。到达时间确定部21通过针对ROI内的各深度执行上述减法,确定修正到达时间。在到达时间为负时,到达时间确定部21使修正到达时间为0。图像产生部23基于修正到达时间和到达时间色相对应表产生对剪切波到达图像的色相进行修正后得到的剪切波到达修正图像。
图11(b)示出基于基准线上的剪切波到达时间对图11(a)的剪切波到达图像的色相进行修正后得到的剪切波到达修正图像的一例的图。像图11(b)所示的那样,在区域92中,方位方向的色相(图中是浓淡)急剧地变化。即,图11(b)表示区域92中的剪切波的传播速度比区域92的周边区域慢。
如果采用以上所述的构成,则可以得到以下的效果。
如果采用本实施方式的超声波诊断装置1,则可以基于除去了被检体的体动造成的组织位移量的组织的变量的时间变化,产生根据关心区域内的各位置处的剪切波到达时间分配了色相的剪切波到达图像。由此,可以把在被检体内产生的剪切波的传播的样子(例如,折射、反射等)作为一个图像可视化。由此,可以作为操作者用来把剪切波到达图像中的斑点模样的区域判断为因不能正确除去被检体的体动造成的位移量等而不能正确确定剪切波到达时间的区域的判断材料。另外,操作者可以在剪切波到达图像中从用来评价用来确定剪切波的传播速度的处理对象和组织的硬度的评价对象中除去斑点模样的区域。
另外,如果采用本实施方式的超声波诊断装置1,则可以抽出除去了被检体的体动造成的组织位移量的组织位移量的时间变化的极大值中的、具有最大位移量的值,基于该极大值具有的位移量确定到达时间,产生根据关心区域内的各位置处的剪切波到达时间分配了色相的剪切波到达图像。由此,例如,在位移量观测扫描处理中,即使超声波探头和被检体中的至少一个移动了时,也可以确定合适的到达时间。即,连续地产生整体移动时,有时表观上的组织位移量伴随着时间(发送次数)的经过而单调增加。此时如果简单地抽出最大位移量,则在单调增加的比例大时,会把因单调增加而增加的值作为最大位移量抽出。但是,使用因单调增加而增加的位移量,则使用产生了整体移动的时刻的位移量算出到达时间,作为结果会确定物理上不正确的到达时间。于是,通过抽出极大值中的具有最大位移量的值,基于该极大值具有的位移量确定到达时间,就可以不使用前面所述的由单调增加得到的最大值,而使用对于确定到达时间是合适的时刻的位移量。由此,在到达时间的确定中鲁棒性提高。
根据以上可知,如果采用本实施方式,则操作者即使不把将各时刻的位移量映像化了的图像作为动画观察,也可以用一幅剪切波到达图像把握剪切波的传播的样子。另外,操作者通过在剪切波到达图像中观察方位方向的色相的变化的程度(色相的变化的缓急),可以定性地把握预定区域中的各位置处的剪切波的传播速度,即各位置处的组织的硬度。
另外,如果采用本实施方式的变形例,则可以产生以操作者经由输入部9输入的指定位置处的剪切波到达时间作为基准,对剪切波到达时间的色相进行修正后得到的剪切波到达修正图像。由此,即使在操作者关注的关注区域从剪切波的产生位置离开时,也可以简便地显示在关注区域中传播的剪切波的样子。而且,如果采用本实施方式的变形例,则也可以把作为剪切波到达时间的基准的位置指定为线段。由此,在剪切波到达关注区域之前,即使因剪切波的传播速度的不均匀性或剪切波的折射等造成剪切波的波面紊乱,也可以把在关注区域中把波面的紊乱初始化。其结果,操作者可以直观地、容易地显示关注区域中的剪切波的传播的样子。
(实施方式2)
与实施方式1的不同之处在于,用预定区域(关心区域)内的各位置处的剪切波到达时间,计算预定区域内的各位置处的剪切波的传播速度,产生根据计算得到的传播速度分配了色相的剪切波传播速度图像。
图12是示出根据实施方式2的超声波诊断装置的构成的一例的构成图。
速度计算部22用预定区域内的各位置处的剪切波到达时间(到达时间数据)计算预定区域内的各位置处的剪切波的传播速度。具体地,速度计算部22在预定区域内的各位置处计算相邻的位置处的剪切波到达时间的差(以下称为到达时间差)。另外,速度计算部22同时计算相邻的位置间的距离。然后,速度计算部22计算预定区域内的各位置处的到达时间差的倒数与邻接的位置间的距离的积,作为剪切波传播速度。速度计算部22把预定区域的各位置处的剪切波传播速度的数据(以下称为传播速度数据)输出到图像产生部23。
图像产生部23基于传播速度数据和在存储部27中存储的传播速度色相对应表,产生在预定区域的各位置处根据剪切波的传播速度分配了色相的剪切波传播速度图像。传播速度色相对应表是例如传播速度的值与色相的对应表。例如,传播速度为0时的色相是蓝色。例如伴随着传播速度的增大,色相以从蓝到蓝绿、绿、黄绿、黄、橙、红的顺序变化,最大传播速度被定义为红色。
图像合成部25产生使剪切波传播速度图像在B模式图像上定位地重叠的传播速度重叠图像。图像合成部25把传播速度重叠图像或剪切波传播速度图像输出到显示部7。
存储部27存储传播速度数据、传播速度色相对应表、剪切波传播速度图像、传播速度重叠图像、与剪切波传播速度图像的产生有关的算法(以下称为剪切波传播速度图像产生算法)等。
CPU31基于由操作者经由输入部9输入的剪切波传播速度图像显示模式的选择,读出存储部27中存储的发送延迟模式、接收延迟模式、剪切波产生发送延迟模式和装置控制程序,根据它们控制装置本体5。例如CPU31从存储部27读出剪切波传播速度图像产生算法。CPU31根据读出的剪切波传播速度图像产生算法控制速度计算部22、图像产生部23等。剪切波传播速度图像显示模式是执行产生、显示剪切波传播速度图像的模式。
具体地,如果经由输入部9输入剪切波传播速度图像显示模式,则CPU31控制速度计算部22以基于到达时间数据产生传播速度数据。CPU31控制图像产生部23,以基于传播速度数据和从存储部27读出的传播速度色相对应表,产生剪切波传播速度图像。
显示部7基于来自图像合成部25的输出显示剪切波传播速度图像、传播速度重叠图像等。
(剪切波传播速度图像产生功能)
剪切波传播速度图像产生功能是产生以色相示出预定区域(第2区域)的各位置处的剪切波的传播速度的剪切波传播速度图像的功能。以下,说明产生剪切波传播速度图像的功能的处理(以下称为剪切波传播速度图像产生处理)。
图13是示出剪切波传播速度图像产生处理的步骤的一例的流程图。
基于预定区域内的各位置处的剪切波到达时间,计算各位置处的剪切波的传播速度(步骤Sd1)。基于各位置处的传播速度和传播速度色相对应表,产生在预定区域的各位置处根据剪切波的传播速度分配了色相的剪切波传播速度图像(步骤Sd2)。
(第1变形例)
与实施方式2的不同之处在于,计算利用操作者的指示在预定区域内指定的指定区域内的剪切波的平均传播速度(以下称为区域内平均传播速度),显示区域内平均传播速度。
输入部9向装置本体5输入由操作者在预定区域内指定的指定区域。
速度计算部22基于指定区域内的各位置处的剪切波到达时间和从剪切波产生位置到指定区域内的各位置的距离,计算区域内平均传播速度。例如,速度计算部22基于在指定区域内且在不同的扫描线上的同一深度的位置处的多个剪切到达时间和从剪切波产生位置到各位置的距离,通过线性回归分析计算区域内平均传播速度。
图像产生部23与剪切波到达图像一起产生显示区域内平均传播速度的图像。
图14是纵轴为位移量、横轴为时间地画出指定区域内的不同的4个扫描线上的同一深度的各位置处的剪切波传播数据得到的图。图14中的曲线64、65、66、67分别表示作为项的扫描线上的同一深度中的剪切波传播数据。图14中的曲线64、65、66、67例如与图3中的扫描线54、55、56、57对应。到达时间确定部21针对图14的曲线64、65、66、67分别确定远离超声波探头的方向(图14中的负的位移量)的最大值的时间,作为剪切波到达时间。
图15是示出剪切波到达时间与到剪切波产生位置的距离的关系的一例的图。图15中的点74表示针对对图14中的曲线64确定的剪切波到达时间的到剪切波产生位置的距离。图15中的点75表示针对对图14中的曲线65确定的剪切波到达时间的到剪切波产生位置的距离。图15中的点76表示针对对图14中的曲线66确定的剪切波到达时间的到剪切波产生位置的距离。图15中的点77表示针对对图14中的曲线67确定的剪切波到达时间的到剪切波产生位置的距离。图15中的直线78是由使用了点74、75、76、77的线性回归分析确定的近似直线。速度计算部21计算线性回归分析得到的直线的斜率,作为区域内平均传播速度。
(第2变形例)
与实施方式2中的第1变形例的不同之处在于,用计算得到的区域内平均传播速度计算指定区域内的平均弹性模量。
速度计算部22用计算得到的区域内平均传播速度计算平均弹性模量(杨氏模量)。具体地,速度计算部22以区域内平均传播速度(vs(m/s))的平方(vs 2)的3倍(3×vs 2)计算平均弹性模量(E(kPa))。即,速度计算部22计算用的公式为E=3×vs 2。
图像产生部23与剪切波到达图像一起产生用来表现平均弹性模量的图像。
显示部7与剪切波到达图像一起显示用来表现平均弹性模量的图像。
如果采用以上所述的构成,则可以得到以下的效果。
如果采用本实施方式的超声波诊断装置1,则用预定区域(关心区域)内的各位置处的剪切波到达时间计算预定区域内的各位置处的剪切波的传播速度,可以产生根据计算得到的传播速度分配了色相的剪切波传播速度图像。由此,可以把在被检体内产生的剪切波的传播速度的空间的分布可视化。另外,操作者,在剪切波的传播速度变成有杂音地不均匀的区域,可以判断为剪切波的传播没有顺利实现。而且,在剪切波到达图像上,指定剪切波传播的区域或正确实现了剪切波的传播的区域,作为用来解析剪切波的平均传播速度和平均弹性模量的区域,可以计算剪切波的平均传播速度和平均弹性模量。由此,可以提高对于剪切波的平均传播速度和平均弹性模量等的定量的解析结果的可靠性。
另外,作为实施方式1的变形例,在用医用图像处理装置100实现本超声波诊断装置1的技术思想时,具有例如图1的构成图中的虚线内的构成要素。另外,医用图像处理装置100也可以在存储部27中存储第1、第2位移数据来执行以后的处理。此时不需要位移量计算部19。另外,医用图像处理装置100也可以在存储部27中存储剪切波传播数据来执行以后的处理。此时不需要位移量计算部19和到达时间确定部21。
另外,作为实施方式2的变形例,在用医用图像处理装置101实现本超声波诊断装置1的技术思想时,具有例如图12的构成图中的虚线内的构成要素。另外,在医用图像处理装置101中,也可以读入从超声波诊断装置1输出的剪切波到达图像的数据等来执行上述处理。
(实施方式3)
与实施方式1、2的不同之处在于,多次执行与多个观测扫描线的每一个对应地发送的推式脉冲的发送次数。作为具体例,说明针对在第1区域内深度不同的多个焦点发送多个推式脉冲的情况。
发送部11在第1区域内发送与深度不同的多个焦点分别对应的多个推式脉冲。通过发送多个推式脉冲,在多个焦点的每一个处产生剪切波。在第2区域中,传播由在多个焦点处分别产生的剪切波重叠得到的重叠剪切波。发送部11例如在第1区域内对深度不同的4个焦点分别发送推式脉冲。具体地,在图2的步骤Sa3中,发送部11对4个焦点的每一个发送推式脉冲。
由多普勒处理部17的自相关器或互相关器得到的第2位移数据具有重叠剪切波造成的组织位移量。
位移量计算部19通过从第1、第2位移数据减去体动位移数据,计算示出重叠剪切波造成的组织位移量的时间变化的重叠剪切波传播数据。
到达时间确定部21基于重叠剪切波传播数据中的组织位移量的时间变化,确定剪切波到达时间。具体地,到达时间确定部21针对预定区域内的各位置,以任一个推式脉冲的发送开始时或任一个推式脉冲的发送结束时作为时刻0,确定与重叠剪切波传播数据中的组织位移量的最大值对应的时间,作为剪切波到达时间。到达时间确定部21把预定区域内的各位置处的到达时间数据输出到图像产生部23。另外,可以在任一个推式脉冲的发送开始时刻或发送结束时刻,任意选择剪切波到达时间的基准时刻。
如果采用以上所述的构成,则可以得到以下的效果。
如果采用本实施方式的超声波诊断装置1,则可以基于除去了被检体的体动造成的组织位移量的组织的变量的时间变化,产生根据关心区域内的各位置处的重叠剪切波的剪切波到达时间分配了色相的剪切波到达图像。由此,可以把在被检体内产生的剪切波的传播的样子(例如,折射、反射等)作为一个图像可视化。
经由网络与医用图像处理装置100、101连接的超声波诊断装置在超声波检查中执行数据的取得和保存即可。本医用图像处理装置100、101可以用于在超声波检查后详细观察剪切波到达图像。即,在进行超声波检查的检查者和解析由超声波诊断装置产生的医用图像的解析者不同时,显现出本实施方式的效果。即,解析者可以在超声波检查后,利用剪切波到达图像确认剪切波是否像解析者所希望的那样传播。解析者可以基于剪切波到达图像计算预定区域内的组织的弹性模量。因此,可以提高例如作为解析结果的弹性模量的可靠性。另外,还可以避免靠检查者的主观确定解析的范围的情况,可以减少误诊的可能性。
而且,实施方式中的各功能也可以通过把执行该处理的程序安装到工作站等的计算机中并把它们在存储器上展开来实现。此时,也可以把能够使计算机执行该方法的程序保存到磁盘(硬盘等)、光盘(CD-ROM、DVD等)、半导体存储器等的存储媒体中并分发。
虽然说明了本发明的几个实施方式,但这些实施方式都是作为例子提出的,并非用来限定本发明的范围。这些新的实施方式可以以其它的各种方式实施,在不脱离发明的主要构思的范围内,可以进行各种省略、改写、变更。这些实施方式及其变形都包含在发明的范围和主要构思内,且包含在权利要求书记载的发明及其等价的范围内。
Claims (11)
1.一种超声波诊断装置,包括:
超声波探头,具有多个振子;
发送部,经由上述超声波探头向第1区域发送在被检体内产生剪切波的第1超声波,向上述被检体内的关心区域发送第2超声波;
接收部,基于上述第2超声波产生接收信号;
位移量计算部,用上述接收信号,计算由上述剪切波造成的组织的位移量;
到达时间确定部,基于上述位移量的时间变化,关于上述关心区域内的多个位置的每一个位置确定上述剪切波到达的到达时间;以及
图像产生部,产生表示上述关心区域内的多个位置的每一个位置处的上述到达时间的剪切波到达图像。
2.如权利要求1所述的超声波诊断装置,其中,
上述到达时间确定部基于上述位移量中的最大位移量确定上述到达时间。
3.如权利要求1所述的超声波诊断装置,其中,
上述到达时间确定部基于上述位移量的时间变化中的极大值中的最大的位移量确定上述到达时间。
4.如权利要求1所述的超声波诊断装置,其中,
上述到达时间确定部基于上述位移量的时间变化中的最大的时间变化确定上述到达时间。
5.如权利要求1所述的超声波诊断装置,其中,
上述位移量计算部作为上述位移量计算沿远离上述超声波探头的方向的位移量。
6.如权利要求1所述的超声波诊断装置,其中,上述超声波诊断装置还包括:
输入部,向上述关心区域输入由操作者指定的指定区域,
上述到达时间确定部确定以上述指定区域中的到达时间为基准对上述关心区域内的多个位置的每一个位置处的到达时间进行修正后得到的修正到达时间,
上述图像产生部基于上述修正到达时间和根据上述关心区域内的多个位置的每一个位置处的到达时间分配的色相产生对上述剪切波到达图像中的上述色相进行修正后得到的剪切波到达修正图像。
7.如权利要求1所述的超声波诊断装置,其中,上述超声波诊断装置还包括:
速度计算部,用上述关心区域内的多个位置的每一个位置处的到达时间计算上述关心区域内的多个位置的每一个位置处的上述剪切波的传播速度,
上述图像产生部基于计算出的上述传播速度和根据上述传播速度预先设定的色相产生根据上述关心区域内的多个位置的每一个位置处的上述传播速度分配了上述预先设定的色相的剪切波传播速度图像。
8.如权利要求1所述的超声波诊断装置,其中,上述超声波诊断装置还包括:
输入部,向上述关心区域输入由操作者指定的指定区域;以及
速度计算部,用上述指定区域内的各位置处的上述到达时间计算在上述指定区域内传播的上述剪切波的传播速度。
9.如权利要求1所述的超声波诊断装置,其中,上述超声波诊断装置还包括:
输入部,向上述关心区域输入由操作者指定的指定区域;以及
速度计算部,用在上述指定区域内传播的上述剪切波的传播速度计算上述指定区域内的弹性模量。
10.一种医用图像处理装置,包括:
存储部,存储通过对被检体的关心区域发送接收超声波而得到的接收信号;
位移量计算部,用上述接收信号计算由在上述被检体内传播的剪切波造成的组织的位移量;
到达时间确定部,基于上述位移量的时间变化,关于上述关心区域内的多个位置的每一个位置确定上述剪切波到达的到达时间;以及
图像产生部,产生表示上述关心区域内的多个位置的每一个位置处的上述到达时间的剪切波到达图像。
11.一种医用图像处理方法,包括:
存储通过对被检体的关心区域发送接收超声波而得到的接收信号;
用上述接收信号计算由在上述被检体内传播的剪切波造成的组织的位移量;
基于上述位移量的时间变化,关于上述关心区域内的多个位置的每一个位置,确定上述剪切波到达的到达时间;以及
产生表示上述关心区域内的多个位置的每一个位置处的上述到达时间的剪切波到达图像。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2013/075207 WO2015040710A1 (ja) | 2013-09-18 | 2013-09-18 | 超音波診断装置、医用画像処理装置および医用画像処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104640506A CN104640506A (zh) | 2015-05-20 |
CN104640506B true CN104640506B (zh) | 2017-06-30 |
Family
ID=52668589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380001735.XA Active CN104640506B (zh) | 2013-09-18 | 2013-09-18 | 超声波诊断装置、医用图像处理装置和医用图像处理方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10959704B2 (zh) |
CN (1) | CN104640506B (zh) |
WO (1) | WO2015040710A1 (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10327738B2 (en) * | 2015-03-31 | 2019-06-25 | Samsung Medison Co., Ltd. | Ultrasound imaging apparatus and method of processing ultrasound image thereof |
JP6390516B2 (ja) * | 2015-05-27 | 2018-09-19 | コニカミノルタ株式会社 | 超音波診断装置及び超音波診断装置の制御方法 |
US10537305B2 (en) * | 2015-09-30 | 2020-01-21 | General Electric Company | Detecting amniotic fluid position based on shear wave propagation |
WO2017069068A1 (ja) | 2015-10-23 | 2017-04-27 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
US10561402B2 (en) * | 2016-07-29 | 2020-02-18 | Siemens Medical Solutions Usa, Inc. | Motion imaging with multiple parallel receive beams |
KR102646993B1 (ko) * | 2016-09-12 | 2024-03-14 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 초음파 영상 장치의 제어 방법 |
CN107440740B (zh) | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN107505232B (zh) * | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | 运动信息获取方法及装置 |
WO2020037674A1 (zh) * | 2018-08-24 | 2020-02-27 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像装置和方法、超声弹性检测装置和方法 |
CN109589139A (zh) * | 2018-12-06 | 2019-04-09 | 深圳开立生物医疗科技股份有限公司 | 一种超声测量生物量自动确认方法和超声诊断系统 |
JP7168474B2 (ja) * | 2019-01-31 | 2022-11-09 | 富士フイルムヘルスケア株式会社 | 超音波撮像装置、治療支援システム、及び、画像処理方法 |
KR20200110960A (ko) * | 2019-03-18 | 2020-09-28 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 그 제어 방법 |
JP7236312B2 (ja) * | 2019-04-04 | 2023-03-09 | 富士フイルムヘルスケア株式会社 | 超音波診断装置、信号処理装置、および、プログラム |
JP7354632B2 (ja) | 2019-07-12 | 2023-10-03 | コニカミノルタ株式会社 | 超音波診断装置、および、超音波診断装置の制御方法 |
US20210128100A1 (en) * | 2019-10-31 | 2021-05-06 | Canon Medical Systems Corporation | Analyzing apparatus and ultrasound diagnosis apparatus |
CN111772677B (zh) * | 2020-07-07 | 2023-06-13 | 意领科技有限公司 | 具有维度的生物组织弹性检测方法,检测系统和存储介质 |
CN113827278B (zh) * | 2021-10-27 | 2023-09-15 | 青岛海信医疗设备股份有限公司 | 剪切波传播速度的确定方法及装置 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1383374A (zh) * | 2000-04-26 | 2002-12-04 | 皇家菲利浦电子有限公司 | 用于切变波参数估算的超声波方法和系统 |
JP4564130B2 (ja) * | 2000-05-23 | 2010-10-20 | 株式会社東芝 | 超音波画像診断装置 |
JP4805669B2 (ja) * | 2005-12-27 | 2011-11-02 | 株式会社東芝 | 超音波画像処理装置及び超音波画像処理装置の制御プログラム |
US8187187B2 (en) | 2008-07-16 | 2012-05-29 | Siemens Medical Solutions Usa, Inc. | Shear wave imaging |
US9554770B2 (en) * | 2008-09-29 | 2017-01-31 | Siemens Medical Solutions Usa, Inc. | High pulse repetition frequency for detection of tissue mechanical property with ultrasound |
JP5455379B2 (ja) * | 2009-01-07 | 2014-03-26 | 株式会社東芝 | 医用画像処理装置、超音波診断装置、及び医用画像処理プログラム |
WO2011001776A1 (ja) * | 2009-07-02 | 2011-01-06 | 株式会社 日立メディコ | 超音波診断装置、せん断波の伝搬画像生成方法 |
JP5559788B2 (ja) | 2009-07-07 | 2014-07-23 | 株式会社日立メディコ | 超音波診断装置 |
BR112012012230B1 (pt) * | 2009-11-25 | 2021-04-20 | Koninklijke Philips N.V. | sistema de formação de imagem diagnóstica ultrassônica para análise de onda de cisalhamento e método para operar um sistema de formação de imagem diagnóstica ultrassônica para medir as ondas de cisalhamento |
JP5846755B2 (ja) * | 2010-05-14 | 2016-01-20 | 株式会社東芝 | 画像診断装置及び医用画像表示装置 |
JP5597455B2 (ja) * | 2010-06-25 | 2014-10-01 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
US8961418B2 (en) * | 2010-10-06 | 2015-02-24 | Siemens Medical Solutions Usa, Inc. | Solving for shear wave information in medical ultrasound imaging |
US10338203B2 (en) * | 2011-09-09 | 2019-07-02 | Siemens Medical Solutions Usa, Inc. | Classification preprocessing in medical ultrasound shear wave imaging |
US9468421B2 (en) * | 2012-02-16 | 2016-10-18 | Siemens Medical Solutions Usa, Inc. | Visualization of associated information in ultrasound shear wave imaging |
JP6462340B2 (ja) * | 2013-12-13 | 2019-01-30 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置、画像処理装置及び画像処理方法 |
CN104605890B (zh) * | 2014-12-18 | 2017-04-05 | 深圳开立生物医疗科技股份有限公司 | 一种剪切波峰值波形校正的方法、装置、系统及其应用 |
CN105232085B (zh) * | 2015-11-18 | 2017-10-31 | 中国人民解放军第三军医大学第三附属医院 | 基于动态孔径控制的超声剪切波弹性成像方法 |
-
2013
- 2013-09-18 CN CN201380001735.XA patent/CN104640506B/zh active Active
- 2013-09-18 WO PCT/JP2013/075207 patent/WO2015040710A1/ja active Application Filing
-
2014
- 2014-04-04 US US14/244,969 patent/US10959704B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US10959704B2 (en) | 2021-03-30 |
CN104640506A (zh) | 2015-05-20 |
WO2015040710A1 (ja) | 2015-03-26 |
US20150080730A1 (en) | 2015-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104640506B (zh) | 超声波诊断装置、医用图像处理装置和医用图像处理方法 | |
JP4745133B2 (ja) | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム | |
JP6058295B2 (ja) | 超音波診断装置、医用画像処理装置、医用画像処理方法、および医用画像処理プログラム | |
JP7299014B2 (ja) | 解析装置、及び解析プログラム | |
CN104905812B (zh) | 用于显示对象的多个不同图像的方法和设备 | |
EP3463098B1 (en) | Medical ultrasound image processing device | |
CN105407809A (zh) | 超声波诊断装置以及弹性评价方法 | |
JP5417048B2 (ja) | 超音波診断装置、及び超音波診断プログラム | |
CN103648400A (zh) | 超声波诊断装置以及方法 | |
JP3946815B2 (ja) | 超音波診断装置 | |
JP2008183063A (ja) | 医用画像診断装置、医用画像表示装置及びプログラム | |
US20100324420A1 (en) | Method and System for Imaging | |
JP2005087634A (ja) | 超音波診断装置及び計測データ表示方法 | |
JP2012030053A (ja) | 超音波診断装置、画像処理装置および画像処理方法 | |
KR20180070990A (ko) | 초음파 영상장치 및 그 제어방법 | |
EP3364881B1 (en) | Ultrasound imaging apparatus and controlling method for the same | |
KR20150118493A (ko) | 초음파 장치 및 그 제어 방법 | |
JP6358192B2 (ja) | 超音波診断装置、及び超音波診断装置の制御方法 | |
JP2019154571A (ja) | 解析装置、及び制御プログラム | |
JPWO2012090658A1 (ja) | 超音波診断装置および画像処理方法 | |
JP4350214B2 (ja) | 超音波診断装置 | |
JP2013141515A (ja) | 医用画像装置及び医用画像構成方法 | |
JP5443781B2 (ja) | 超音波診断装置 | |
JP7292896B2 (ja) | 解析装置、及び解析方法 | |
CN104869910B (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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20160712 Address after: Japan Tochigi Applicant after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Applicant before: Toshiba Corp Applicant before: Toshiba Medical System Co., Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |