CN101170948A - 超声波诊断装置及其图像处理方法 - Google Patents

超声波诊断装置及其图像处理方法 Download PDF

Info

Publication number
CN101170948A
CN101170948A CNA200680014972XA CN200680014972A CN101170948A CN 101170948 A CN101170948 A CN 101170948A CN A200680014972X A CNA200680014972X A CN A200680014972XA CN 200680014972 A CN200680014972 A CN 200680014972A CN 101170948 A CN101170948 A CN 101170948A
Authority
CN
China
Prior art keywords
profile
frame
speckle
left ventricle
diagnostic ultrasound
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CNA200680014972XA
Other languages
English (en)
Other versions
CN101170948B (zh
Inventor
长野智章
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Healthcare Corp
Original Assignee
Hitachi Medical Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Publication of CN101170948A publication Critical patent/CN101170948A/zh
Application granted granted Critical
Publication of CN101170948B publication Critical patent/CN101170948B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • 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/48—Diagnostic techniques
    • A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
    • 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/52087—Details related to the ultrasound signal acquisition, e.g. scan sequences using synchronization techniques
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • 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/52074—Composite displays, e.g. split-screen displays; Combination of multiple images or of images and alphanumeric tabular information
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Cardiology (AREA)
  • Quality & Reliability (AREA)
  • Physiology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)

Abstract

一种超声波诊断装置，包括向被检体内收发超声波并对该被检体的动态图像进行拍摄的单元，还包括：斑点测量单元，其对所述动态图像的各帧，测量各帧上出现的斑点的大小及/或形状；和平滑化单元，其根据该测量的斑点的大小及/或形状，对各帧的图像数据实施平滑化处理。

Description

超声波诊断装置及其图像处理方法
技术领域
本发明涉及超声波诊断装置及其图像处理方法，尤其涉及考虑斑点的形状以及大小而能获得有用的超声波图像的超声波诊断装置及其图像处理方法。
背景技术
在由超声波诊断装置获得的超声波图像中，混有称为斑点噪声的噪音。认为由于与超声波的波长相比非常小的生物体组织内的反射体群反射的散射波以各种相位产生并干涉，从而出现该斑点噪声。
一般，历来斑点噪声对于图像诊断而言是不需要的噪声，因此认为是应降低的部分。例如，根据专利文献1记载的现有技术，具备判定斑点噪声并将其除去的电路。
专利文献1：特开平9-94248号公报
但是，本发明者们认为超声波图像中的斑点未必是在诊断被检体时不需要的信息。即，若实施与图像上出现的斑点的形状以及大小相应的滤波处理，则认为可不设置斑点噪声的除去电路而获得对诊断更为有用的图像。
发明内容
本发明的目的在于，提供一种考虑斑点的形状以及大小来进行滤波处理，从而可获得更有用的超声波图像的超声波诊断装置及其图像处理方法。
更具体而言，提供一种有助于心脏病变的适宜的诊断的超声波诊断装置及其图像处理方法，尤其是，提供一种可获得心脏的多个部位、即左心室、心肌、左心房、右心室、右心房等四腔的良好质量的图像，并对这些功能进行评价的超声波诊断装置及其图像处理方法。
为了实现上述目的，本发明的超声波诊断装置，包括向被检体内收发超声波并对该被检体的动态图像进行拍摄的单元，还包括：斑点测量单元，其对所述动态图像的各帧，测量各帧上出现的斑点的大小及/或形状；和平滑化单元，其根据该测量的斑点的大小及/或形状，对各帧的图像数据实施平滑化处理。
另外，本发明的超声波图像处理方法，包括(1)向被检体内收发超声波并对该被检体的动态图像进行拍摄的工序，该方法还包括：(2)对所述动态图像的各帧，测量各帧上出现的斑点的大小及/或形状的工序；和(3)根据该测量的斑点的大小及/或形状，对各帧的图像数据实施平滑化处理的工序。
附图说明
图1是本发明的实施例1的超声波诊断装置的系统构成图；
图2是本发明的实施例1的超声波诊断装置1的处理次序；
图3是表示在超声波图像(B模式)上设定有朝超声波束方向的窗口的情况的图；
图4是表示理想的斑点的轮廓的图；
图5是以椭圆形显示斑点的图；
图6是表示斑点彼此的间隔窄、对比度饱和的情况的图；
图7是表示二维高斯滤波器的特性的一例的图；
图8是表示完成了手动跟踪(manual trace)的状态的图；
图9是表示通过步骤26进行了轮廓点83～85的凹凸和间隔的偏差的修正后的结果的图；
图10是说明辛普森(Simpson)法的图；
图11是以曲线图表示各参数随着帧的更新而如何变化(时间变化)的图；
图12是表示实施例2中的显示例的图；
图13是表示实施例3中的显示例的图；
图14是表示显示装置上的其他显示例的图。
具体实施方式
下面，参照附图进行说明。
(实施例1)
图1是本发明的实施例1的超声波诊断装置的系统构成图。
在图1中，本发明的实施例1中的超声波诊断装置1是利用超声波进行心脏功能的测量用的装置，采用至少其一部分中包括公知的超声波诊断装置的装置构成。
超声波诊断装置1构成为包括：探头2、发送部3、接收部4、收发分离部、调相加法运算部6、信号处理部7、A/D转换器8、帧存储器9a、视频存储器(cine memory)9b、控制器10、输入设备11、接口12、结果存储部13、显示电路部14、显示装置15、心电图16。在图1中，仅表示了超声波诊断装置的主要功能。下面，对图1所示的各构成进行说明。
探头2构成为：向诊断部位(在此为心脏)送出超声波，并能接受其反射波。在探头2的内部设置有未图示的多个振动器，这些振动器成为超声波的产生源(送波源)，另外还能接受反射波。发送部3构成为：能对探头2进行驱动来生成送出超声波用的送波脉冲信号。另一方面，接收部4构成为：可接收由探头2接受并被转换为电信号的回声信号。
收发分离部5构成为：能够在发送时向探头2送出来自发送部3的送波脉冲信号，在接收时向接收部4送出来自探头2的回声信号。调相加法运算部6构成为：能够对来自接收部4的多个回声信号进行调相加法运算，生成接收射束信号。
信号处理部7、A/D转换器8、帧存储器9a和视频存储器9b构成为：作为根据所述接收射束信号来获得诊断部位的浓淡断层像(黑白断层像)用的信号处理部而发挥功能。即，信号处理部7输入来自调相加法运算部6的接收射束信号，进行增益修正、对数压缩、检波、增强轮廓、滤波处理等信号处理。A/D转换器8构成为：可将从信号处理部7输出的信号变换为数字信号。帧存储器9a构成为：可将从A/D转换器输出的接收射束信号以图像帧单位进行存储。另外，视频存储器9b构成为：可存储连续拍摄的多个图像帧。此外，存储到帧存储器9a和视频存储器9b的图像与由心电计16测量的ECG波形的相位信息对应。
帧存储器9a中存储的上述断层帧数据被基于控制器10的控制信号以电视同步方式读出。另外，控制器10构成为可进行如下各种处理：基于控制程序进行各构成要素的控制；将从视频存储器9b读出的断层帧数据加工成超声波断层像数据；进行后述的轮廓点或各轮廓线的数据生成或向显示装置15的输出控制；进行同样后述的心脏功能测量所涉及的容积计算和距离计算等规定的运算；进行同样后述的修正、平滑化处理或组织追踪。
控制器10构成为具有所谓的微型计算机的功能。控制器10具有如后述的运算单元、运算结果输出单元、平滑化处理单元、组织追踪单元等功能。
输入设备11通过接口12与控制器10连接。作为输入设备11，可举出鼠标或跟踪球作为一例。输入设备11为了由技术人员(操作者)一边参照显示装置15上显示的超声波图像，一边在超声波图像上手动跟踪心脏的左心室、心肌、以及左心房的各轮廓而设置。输入设备11和控制器10具有如后述的跟踪单元以及修正单元的功能。
结果存储部13具有作为存储后述的轮廓点的坐标数据、由控制器运算出的运算结果等的存储器的功能。结果存储部13中存储的坐标数据和运算结果等基于控制器10的控制信号而被读出。
显示电路部14构成为：基于来自控制器10输出的控制信号而工作。显示电路部14构成为：将来自控制器10的超声波断层像数据或后述的轮廓点、各轮廓线的数据变换为模拟信号，来生成显示用的影像信号。在显示电路部14中未特别进行图示，但设置有D/A变换器或影像信号变换电路等。显示装置15构成为：输入从显示电路部14输出的影像信号，显示超声波图像。作为显示装置15，例如使用电视监视器。
接着，参照图2的流程图，对本发明的实施例1的超声波诊断装置1的处理次序进行说明。图2是表示控制器10中的处理的流程图。这里说明的处理中还包含由用户利用输入设备11和显示装置15所进行的各种输入等的次序等。另外，下述所示的图2的流程图的各步骤说明还随时参照图3～图11进行。
(步骤21)
基于操作者利用输入设备11输入的信号，从影片存储器9b读出第一帧的图像，在显示装置15上显示超声波的动态图像的第一帧。
(步骤22)
对步骤21中显示的第一帧图像通过后述的方法实施滤波处理，改善画质。
下面，详细描述在该步骤中执行的滤波处理。
这里，在步骤32中，详细描述用于使进行组织追踪等时的运算(微分运算等)更适合的滤波处理。该滤波处理由步骤22a和步骤22b构成。
一般，为了适合进行微分运算，对图像数据实施图像浓度的平滑化滤波。作为平滑化滤波器之一有二维高斯滤波器。在超声波图像中，由于在超声波束的收发方向和扫描方向(与该收发方向交叉的方向)上分辨率不同，因此，需要实施与该分辨率的差异相对应的二维高斯滤波处理。更具体而言，在超声波图像上，存在称为斑点的浓淡的凹凸(例如，参照特开平7-51270号公报)，但在超声波图像中，该斑点并非完全的圆形，而是椭圆形状(长径方向和短径方向分别为超声波束的扫描方向或发送方向)。因此，本发明者鉴于斑点为椭圆形状，发明了二维高斯滤波器以各向异性方式实施的方法。
下面，对步骤22a和22b进行说明。
(步骤22a)
图3是在超声波图像(B模式)上设定有超声波束方向上的窗口41的情况。首先，根据该窗口41内的图像数据，求出窗口内部的斑点的平均大小及/或形状。作为具体的次序，直接取出窗口内的像素值，按照使射束的发送方向在画面上为垂直方向的方式进行仿射变换(affinetransformation)，然后，在表示为42的窗口内，对水平方向和垂直方向的浓度共生矩阵(例如，参照特开平5-123318)进行运算，求出对比度的特征量。当斑点的大小及/或形状为理想情况时，如图4所示，对比度最高的像素位置(51)与对比度最低的像素位置(52)之间的距离(53)成为斑点的大小(54)的二分之一。因此，通过计算水平方向和垂直方向的对比度最高的像素位置(51)与对比度最低的像素位置(52)之间的距离(53)，求出以椭圆近似时(图5)的斑点的大小的尺寸(短径A以及长径B)。
或者，另外，图像的某线段上(图3的窗口42内的线段上)的轮廓也并非如图4所示那样理想。例如，图6是对浓度共生矩阵进行运算的结果，61是像素间距离，62是对比度，63是横向对比度，64是纵向对比度，由A、B表示的距离65是斑点的短径A以及长径B。在图6中，斑点彼此的间隔窄，浓淡共生矩阵的对比度特征量饱和。在如图6的情况下，检测像素值达到最大的距离(图6中的A以及B)，作为斑点的尺寸(短径A以及长径B)。
(步骤22b)
在该步骤中，利用在步骤22a中求出的斑点的尺寸(短径A以及长径B)，进行与其特性相应的二维高斯滤波器的应用。
图7是二维高斯滤波器(71)的特性的一例。二维高斯滤波器如图7所示，将无论取X轴方向和Y轴方向上的哪个截面都为正态分布的函数作为滤波器，但若在XY平面上截断图7所示的函数，则其截面成为椭圆，与在步骤22a中求出的斑点同样。在该步骤中，将步骤22a中的短径A和长径B的大小调整为图7所示的二维高斯滤波器的X轴方向和Y轴方向的标准偏差，从而使平滑化处理最佳化。
关于相对在步骤22a中求出的短径A和长径B如何调整二维高速滤波器的X轴方向和Y轴方向的标准偏差σA和σB，为了不产生不需要的噪声，设定为比短径A和长径B的长度稍长是重要的。另外，若使标准偏差σA和σB的设定过大，则存在过于平滑化导致失去图像的特征的问题，因此需要调整为适当的值。
通过上述的滤波处理(步骤22a以及22b)，可描绘到心脏的左心房(或右心房)。
(步骤23)
操作者一边观察显示装置15上显示的第一帧超声波图像，一边利用鼠标或跟踪球所构成的输入设备11，开始心脏四腔的手动跟踪。在本实施例中，通过步骤22a和22b，可清晰地描绘到心脏的左心房，因此，到左心房为止的手动跟踪容易。这里，手动跟踪是指由操作者在超声波图像上以点(轮廓点)来描绘左心室、心肌、左心房的轮廓(更具体而言，左心室内膜、左心室外膜、左心房的轮廓)。另外，在本实施例中，将一对阀环(左心室与左心房的结合点)作为对左心室和左心房进行手动跟踪时的公共部分。
这里的手动跟踪的具体次序的一例为：在一方的阀环的位置(例如，图8中的81)配置轮廓点，从该轮廓点沿左心室内膜依次配置多个轮廓点。在沿左心室内膜配置了多个轮廓点之后，在另一方的阀环的位置(例如，图8中的82)配置轮廓点。同样，在一方的阀环的位置配置轮廓点，从该轮廓点沿左心室外膜依次配置多个轮廓点。在沿左心室外膜配置了多个轮廓点之后，在另一方的阀环的位置配置轮廓点。进而，在一方的阀环的位置配置轮廓点，从该轮廓点沿左心房依次配置多个轮廓点。在沿左心房配置了多个轮廓点之后，在另一方的阀环的位置配置轮廓点。
这里所示的手动跟踪的次序作为一例，从哪个轮廓开始手动跟踪均可。另外，进行画面的右向旋转和左向旋转的哪个朝向的手动跟踪(轮廓点的配置)均可。另外，也可仅对左心房进行手动跟踪，而对左心室和心肌利用现有技术(例如，特开平8-206117号公报所公开的技术等)进行自动跟踪。
(步骤24)
通过步骤23由输入设备11输入而配置的轮廓点以叠加在超声波图像上的形式显示于显示装置15，并伴随于此存储到结果存储部13中。
图8表示完成了手动跟踪的状态的图，配置有多个轮廓点。图8中的81、82表示一对阀环的位置。另外，83是由多个轮廓点形成的左心室内膜的轮廓线，84是左心室外膜的轮廓线，85是左心房的轮廓线。从图8可知，在通过手动跟踪配置了轮廓点的情况下，轮廓线83～85凹凸较多，轮廓点的间隔存在偏差。
(步骤25)
基于控制器10的控制，对手动跟踪的三条轮廓线83～85的凹凸或间隔的偏差进行自动修正。更具体而言，例如通过样条曲线等进行拟合，按照使轮廓点达到预先确定的个数和间隔的方式进行重新配置即可。图9是通过步骤25进行了轮廓线83～85的凹凸或间隔的偏差的修正后的结果，轮廓线变得更圆滑。
(步骤26)
通过步骤25进行了自动修正的结果，当存在操作者明确认为是拟合错误的轮廓点时，利用输入设备11进行手动修正。更具体而言，通过对各个轮廓点进行点击、拖动来手动修正。手动修正后的各轮廓点的坐标数据被重新存储到结果存储部13中。
此外，在上述的步骤23～25的手动跟踪以及手动跟踪后的轮廓点的修正中，阀环的位置可能会因轮廓线83～85而稍微偏移。在这种情况下，可以通过某一轮廓线的阀环的位置来公用，也可求出多个轮廓线的阀环的位置坐标的平均坐标将其作为阀环的位置而公用。由此，左心室以及左心房由一条线连接，从而可对左心室左心房之间流动的血流量等无遗漏地进行测量。另外，由左心室内膜和左心室外膜包围的区域是心肌，因此，也可对心肌区域的面积等无遗留地进行测量。
(步骤27)
判断是否恰当地实施了最初帧的手动跟踪，当判断为已恰当地实施的情况下转移到步骤28，当判断为未恰当实施的情况下转移到步骤21。
(步骤28)
基于到步骤26为止求出的轮廓线，对心脏的各部位的容积以及大小(距离等)进行测量。在本实施例中，例如，为了求出容积而利用辛普森法，通过图10进行其具体说明。首先，求出阀环间的中点101，并检索左心室内膜、左心室外膜、左心房各个轮廓线上距该点最远的点，通过连结该远点和中点101求出轴102、103、104，从而应用辛普森法。利用辛普森法进行脏器的求积的方法例如在特开平7-289545号公报中公开。通过如在特开平7-289545号公报中公开的方法，计算左心室内膜、左心室外膜、左心房、左心室与左心房的和、心肌(左心室外膜容积与左心室内膜容积之差)的各容积。另外，计算轴102的长度(连接在轴102上相对于中点101最向图上侧的与左心室外膜相交的点、和在轴102上相对于中点101最向图下侧延长时与左心房相交的点的长度)、左心室以及左心房的壁间距离(将构成左心室内膜以及左心房的轮廓线的一对阀环间连接起来的线段方向的宽度105、106、心肌的内外膜间的距离107)。另外，还计算轮廓线方向的轮廓点间距离108。
(步骤29)
判断是否存在下一帧。当存在下一帧时前进到步骤30，当不存在下一帧时前进到步骤33。
(步骤30)
控制器10从视频存储器9b读入下一帧图像数据，存储到结果存储部13中。
(步骤31)
对在步骤30中显示的第二帧图像与步骤22同样实施滤波处理，进行画质改善。步骤31由步骤31a和步骤31b构成，分别进行与步骤22a、步骤22b同样的处理。
(步骤32)
在该步骤中，自动追踪从第一帧向第二帧(或者，在步骤30中根据相继实施的帧的读入，从第n帧向第n+1帧)转移时产生的各脏器的轮廓线的变化。这里，将各脏器的轮廓线的变化(变动)的追踪称为组织追踪的处理。本实施例中的组织追踪的处理的具体方法采用鲁棒性高的算法，以便在画质差的情况下也能应用。例如，可采用光学流(optical flow)法，块匹配法、梯度法、粒子追踪法。在梯度法中，具体而言，利用图像浓度的梯度，以解析方式求出速度矢量。由于向图像的访问只要微分计算，因此可高速求出速度矢量。尤其是，在膜的部分可获得大至某一程度的微分值，因此能稳定进行组织追踪。当组织追踪结束时，转移到步骤28，对第二帧(第n+1帧)计算左心室内膜、左心室外膜、左心房、左心室与左心房的和、心肌(左心室外膜容积与左心室内膜容积之差)的各容积、102的长度、左心室以及左心房的壁间距离、心肌的内外膜间的距离、轮廓线方向的轮廓点间距离(以下将这些由步骤28求出的值称为参数)。
(步骤33)
当针对所有帧的处理结束时，在显示装置15上以曲线图等形式显示各参数随着帧的更新如何变化(时间变化)。此外，关于这里的曲线图显示，是在横轴上表示时间或帧编号，在纵轴上表示各参数的计算值，例如如图11所示。
在图11的显示中，显示左心室内膜、左心室外膜、左心房、左心室与左心房的和、心肌(左心室外膜容积-左心室内膜容积)的容积的时间变化，可一边相互参照包括左心房在内的心脏的各部的容积变化和ECG(心电图)，一边进行诊断(在图11中，最上面的线表示左心室外膜、接下来依次是左心室与左心房的和、心肌、左心室内膜、左心房，最下面的线表示ECG(心电图))。另外，左心室内膜、左心室外膜、左心房的各轴(图10中的102、103、104)的长度也可同样地进行曲线化。进而，左心室、左心房的壁间距离105、106、心肌的内外膜间距离107、轮廓线方向的轮廓点间距离108也可进行曲线化。容积、轴长、壁间距离是评价心功能的重要指标，与左心室心肌、左心房膜的运动性能密切相关。如上所述，通过进行曲线化，可观察时段与心功能异常的关系，还可观察左心室、左心房间的心功能的差异。
根据上述实施例，对时间上连续获得的多帧所构成的动态图像中的最初的一帧，通过手动跟踪等确定各脏器的轮廓线，进而在对与最初一帧相连的多帧通过组织追踪来追踪各脏器的轮廓线时，考虑在超声波图像上特征表现的矢量的大小以及形状，进行各图像数据的滤波处理，由此，还可追踪以往不清晰的左心房等的轮廓。另外，基于各帧的各脏器的轮廓线，求取供被检体的诊断用的各种参数，还可显示其时间上的变化，从而可提供适宜的医用图像诊断装置以及方法。
(实施例2)
本实施例是本发明的显示到显示装置15上的显示例的另一例。
如图12所示，还可对心室、左心房、心肌的剖面(121)以时间序列进行三维排列显示。由此，从视觉上可观察心脏的形状的变化。
(实施例3)
本实施例是本发明的显示到显示装置15上的显示例的又一例。
如图14所示，可将心室、左心房、心肌划分为几个部分进行三维显示(131)，并且，通过对帧进行连续显示，可使三维显示的心室、左心房、心肌在时间上变化的情况可视化。
本发明并不限定于上述实施例，在不脱离本发明的宗旨的范围内可进行各种变形来实时。例如，上述步骤22、步骤31中对各帧实施的滤波处理也可安排在步骤21之前对所有帧实施，而不是在步骤21之后、步骤30之后进行。
另外，在步骤22b、31b中实施的二维高斯滤波也可采用不是利用高斯函数的滤波处理，而采用其他函数。
另外，本发明不限于观察心脏的动态，还可应用于观察其他脏器。例如，可用于观察脖颈的颈动脉的脉搏。不仅是具有动态的脏器，对于通常的基于超声波的摄像，通过考虑斑点的大小/及形状来进行平滑化处理，也存在画质改善效果，本发明还可应用于通常的超声波诊断装置以及方法中。
另外，斑点的大小及/或形状在一帧图像数据内也会因位置而不同，因此，可根据基于位置引起的大小及/或形状的变化来改变基于高斯滤波的平滑化处理。
另外，图11～13所示的显示例也可不单独在显示装置15上显示，而是与B模式像排列显示。例如，将B模式像(141)和图11(142)组合的显示例成为图14所示。在图14中，由143表示的线表示图14上侧的B模式像是在下侧的显示例的横轴的时间轴的哪个时刻的像。

Claims (20)

1.一种超声波诊断装置，包括向被检体内收发超声波并对该被检体的动态图像进行拍摄的单元，还包括：
斑点测量单元，其对所述动态图像的各帧，测量各帧上出现的斑点的大小及/或形状；和
平滑化单元，其根据该测量的斑点的大小及/或形状，对各帧的图像数据实施平滑化处理。
2.根据权利要求1所述的超声波诊断装置，其特征在于，
所述斑点的大小及/或形状由具有长轴和短轴的椭圆来近似，所述斑点测量单元通过求取由所述椭圆来近似的斑点的长轴和短轴，由此对所述斑点的大小及/或形状进行测量。
3.根据权利要求1所述的超声波诊断装置，其特征在于，
所述斑点测量单元基于对所述动态图像的各帧的图像数据实施了浓度共生矩阵运算后的结果，测量所述斑点的大小及/或形状。
4.根据权利要求2所述的超声波诊断装置，其特征在于，
所述平滑化处理通过高斯滤波器处理进行，所述高斯滤波器的垂直的两个轴上的标准偏差基于由所述斑点测量单元求出的斑点的长轴和短轴来确定。
5.根据权利要求1所述的超声波诊断装置，其特征在于，
还包括：
参数测量单元，其基于所述动态图像的各帧中的各部位的轮廓，对表示所述运动部位的形状的参数进行测量；和
显示单元，其显示所述参数在时间上的变化。
6.根据权利要求1所述的超声波诊断装置，其特征在于，
所述运动部位是被检体的心脏，所述轮廓是心脏的左心室内膜、左心室外膜、左心房、右心室内膜、右心室外膜、右心房的轮廓。
7.根据权利要求6所述的超声波诊断装置，其特征在于，
还包括如下单元，该单元使作为所述心脏的左心室内膜和左心室外膜和左心房、或者右心室内膜和右心室外膜和右心房的接合点的阀环，在各轮廓中公用。
8.根据权利要求5所述的超声波诊断装置，其特征在于，
所述参数是所述心脏的四腔的容积或轴长、构成心脏的膜间的距离、或心肌的厚度。
9.根据权利要求5所述的超声波诊断装置，其特征在于，
所述测量单元通过辛普森法求取参数。
10.根据权利要求1所述的超声波诊断装置，其特征在于，
还包括如下单元，该单元通过对所述运动部位的剖面进行三维排列来显示所述轮廓在时间上的变化。
11.根据权利要求1所述的超声波诊断装置，其特征在于，
还包括如下单元，该单元通过使所述运动部位的轮廓的三维图像在时间上变化，来显示所述轮廓在时间上的变化。
12.一种超声波图像处理方法，包括(1)向被检体内收发超声波并对该被检体的动态图像进行拍摄的工序，
该方法还包括：
(2)对所述动态图像的各帧，测量各帧上出现的斑点的大小及/或形状的工序；和
(3)根据该测量的斑点的大小及/或形状，对各帧的图像数据实施平滑化处理的工序。
13.根据权利要求12所述的超声波图像处理方法，其特征在于，
还包括：
(4)对所述动态图像的任一帧提取所述运动部位的轮廓的工序；和
(5)对所述动态图像的其他帧检测所述轮廓的运动的工序。
14.根据权利要求12所述的超声波图像处理方法，其特征在于，
所述工序(2)通过对各帧的图像数据实施浓度共生矩阵，来检测所述斑点的大小及/或形状。
15.根据权利要求12所述的超声波图像处理方法，其特征在于，
所述轮廓是心脏的左心室内膜、左心室外膜、左心房、右心室内膜、右心室外膜、右心房的轮廓。
16.根据权利要求15所述的超声波图像处理方法，其特征在于，
所述工序(2)包括：
(6)利用输入单元，由所述动态图像的一帧输入心脏的左心室内膜、左心室外膜、左心房的轮廓点的工序；
(7)将心脏的左心室内膜、左心室外膜、左心房的轮廓点相交的部分作为阀环并使其公用的工序；
(8)圆滑地连接所述轮廓点来导出轮廓线的工序；和
(9)对所述轮廓线的导出进行修正的工序。
17.根据权利要求13所述的超声波图像处理方法，其特征在于，
包括：
(10)检测所述工序(9)中作为轮廓线导出的轮廓在动态图像的各帧中如何运动的工序；
(11)根据所述轮廓的运动，算出与所述运动部位的运动有关的参数的工序；和
(12)显示该参数在时间上的变化的工序。
18.根据权利要求13所述的超声波图像处理方法，其特征在于，
包括：
(13)进行用于适当显示所述轮廓的运动的信号处理的工序；和
(14)基于由所述信号处理获得的结果，对所述轮廓的运动进行显示的工序。
19.根据权利要求18所述的超声波图像处理方法，其特征在于，
在所述工序(14)中，通过对所述运动部位的剖面进行三维排列来显示所述轮廓在时间上的变化。
20.根据权利要求18所述的超声波图像处理方法，其特征在于，
在所述工序(14)中，通过使所述运动部位的轮廓的三维图像在时间上变化，来显示所述轮廓在时间上的变化。
CN200680014972XA 2005-05-19 2006-05-18 超声波诊断装置及其图像处理方法 Active CN101170948B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2005146656 2005-05-19
JP146656/2005 2005-05-19
PCT/JP2006/309902 WO2006123729A1 (ja) 2005-05-19 2006-05-18 超音波診断装置及びその画像処理方法

Publications (2)

Publication Number Publication Date
CN101170948A true CN101170948A (zh) 2008-04-30
CN101170948B CN101170948B (zh) 2010-09-08

Family

ID=37431301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200680014972XA Active CN101170948B (zh) 2005-05-19 2006-05-18 超声波诊断装置及其图像处理方法

Country Status (5)

Country Link
US (1) US20090216124A1 (zh)
EP (1) EP1882450B1 (zh)
JP (1) JP5138369B2 (zh)
CN (1) CN101170948B (zh)
WO (1) WO2006123729A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101912273A (zh) * 2010-07-23 2010-12-15 深圳市理邦精密仪器股份有限公司 超声动态图像处理方法及系统
CN102274051A (zh) * 2011-05-27 2011-12-14 深圳市理邦精密仪器股份有限公司 一种超声图像膀胱容积自动测量方法及系统
CN103732134A (zh) * 2010-12-29 2014-04-16 迪亚卡帝奥有限公司 自动左心室功能评价
CN112043259A (zh) * 2015-07-10 2020-12-08 深圳迈瑞生物医疗电子股份有限公司 一种监护系统、方法及装置
CN112750099A (zh) * 2019-10-16 2021-05-04 深圳市理邦精密仪器股份有限公司 卵泡测量方法、超声设备以及计算机可读存储介质

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101051555B1 (ko) * 2007-11-20 2011-07-22 삼성메디슨 주식회사 개선된 3차원 초음파 영상을 형성하는 초음파 영상 장치 및방법
US20100113929A1 (en) * 2008-10-31 2010-05-06 National Tsing Hua University High-frequency ultrasonic imaging system and method
JP5209026B2 (ja) * 2010-10-27 2013-06-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
JP6274495B2 (ja) * 2013-10-25 2018-02-07 キヤノンメディカルシステムズ株式会社 画像処理装置および超音波診断装置
JP6418881B2 (ja) * 2014-10-08 2018-11-07 ザイオソフト株式会社 パラメータ処理装置、パラメータ処理方法、及びパラメータ処理プログラム
WO2019156139A1 (ja) * 2018-02-08 2019-08-15 興和株式会社 画像処理装置、画像処理方法及び画像処理プログラム
JP7258538B2 (ja) * 2018-12-14 2023-04-17 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用情報処理装置、医用情報処理プログラム
JP7328156B2 (ja) * 2020-01-22 2023-08-16 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用画像処理装置、および医用画像処理プログラム

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0734798B2 (ja) * 1991-11-08 1995-04-19 アロカ株式会社 超音波画像表示装置
US5479926A (en) * 1995-03-10 1996-01-02 Acuson Corporation Imaging system display processor
CA2326596A1 (en) * 1998-03-30 1999-10-07 Stanley Zietz Echocardiography workstation
US6674879B1 (en) * 1998-03-30 2004-01-06 Echovision, Inc. Echocardiography workstation
US6155978A (en) * 1998-12-09 2000-12-05 General Electric Company Three-dimensional imaging by projecting morphologically filtered pixel data
JP4749592B2 (ja) * 2000-05-01 2011-08-17 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー パルス・ドプラ超音波イメージングにおいてサンプル・ゲートを自動設定する方法及び装置
US6312385B1 (en) * 2000-05-01 2001-11-06 Ge Medical Systems Global Technology Company, Llc Method and apparatus for automatic detection and sizing of cystic objects
JP4907798B2 (ja) 2001-08-24 2012-04-04 株式会社東芝 超音波診断装置
US6676605B2 (en) * 2002-06-07 2004-01-13 Diagnostic Ultrasound Bladder wall thickness measurement system and methods
JP2004141514A (ja) * 2002-10-28 2004-05-20 Toshiba Corp 画像処理装置及び超音波診断装置
US8021301B2 (en) * 2003-12-26 2011-09-20 Fujifilm Corporation Ultrasonic image processing apparatus, ultrasonic image processing method and ultrasonic image processing program

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101912273A (zh) * 2010-07-23 2010-12-15 深圳市理邦精密仪器股份有限公司 超声动态图像处理方法及系统
CN103732134A (zh) * 2010-12-29 2014-04-16 迪亚卡帝奥有限公司 自动左心室功能评价
CN103732134B (zh) * 2010-12-29 2016-08-17 迪亚卡帝奥有限公司 用于自动左心室功能评价的系统、装置、设备和方法
CN102274051A (zh) * 2011-05-27 2011-12-14 深圳市理邦精密仪器股份有限公司 一种超声图像膀胱容积自动测量方法及系统
CN102274051B (zh) * 2011-05-27 2012-12-26 深圳市理邦精密仪器股份有限公司 一种超声图像膀胱容积自动测量方法及系统
CN112043259A (zh) * 2015-07-10 2020-12-08 深圳迈瑞生物医疗电子股份有限公司 一种监护系统、方法及装置
CN112750099A (zh) * 2019-10-16 2021-05-04 深圳市理邦精密仪器股份有限公司 卵泡测量方法、超声设备以及计算机可读存储介质

Also Published As

Publication number Publication date
JP5138369B2 (ja) 2013-02-06
EP1882450B1 (en) 2013-01-02
EP1882450A4 (en) 2010-05-26
CN101170948B (zh) 2010-09-08
WO2006123729A1 (ja) 2006-11-23
JPWO2006123729A1 (ja) 2008-12-25
EP1882450A1 (en) 2008-01-30
US20090216124A1 (en) 2009-08-27

Similar Documents

Publication Publication Date Title
CN101170948B (zh) 超声波诊断装置及其图像处理方法
CN103222879B (zh) 用于识别超声成像的最佳图像帧的系统和方法
CN103169500B (zh) 超声波诊断装置、医用图像诊断装置以及医用图像处理方法
JP6160487B2 (ja) 超音波診断装置およびその制御方法
CN101066211B (zh) 用于在超声波系统中显示信息的方法
CN104114102A (zh) 超声波诊断装置、图像处理装置以及图像处理方法
CN104883982A (zh) 针对照护点的解剖学智能回波心动描记
CN102551813B (zh) 超声诊断设备和用于追踪组织移动的方法
JP2002224116A (ja) 超音波診断装置及び画像処理装置
US11931201B2 (en) Device and method for obtaining anatomical measurements from an ultrasound image
KR20110128197A (ko) 심장 ｍ-모드 뷰들의 자동 분석
CN109770943A (zh) 一种利用计算机视觉定位的超声自动优化方法
JP2012090820A (ja) 超音波診断装置
CN107427279A (zh) 使用具有用户控制的心脏模型腔室的心脏功能的超声诊断
EP3977938B1 (en) Method and ultrasound system for shear wave elasticity imaging
US20180192987A1 (en) Ultrasound systems and methods for automatic determination of heart chamber characteristics
CN106030657B (zh) 医学4d成像中的运动自适应可视化
US20210345986A1 (en) Automatic evaluation of ultrasound protocol trees
CN115330681A (zh) 一种超声图像测量方法、系统、智能终端及存储介质
US20180049718A1 (en) Ultrasonic diagnosis of cardiac performance by single degree of freedom chamber segmentation
EP4265191A1 (en) Ultrasound imaging
CN114271850B (zh) 超声检测数据的处理方法及超声检测数据的处理装置
CN116269496A (zh) 基于隐式神经表示的心脏三维超声成像及心功能评估系统
CN101283916B (zh) 超声波诊断装置及超声波诊断方法
US20220071599A1 (en) Ultrasound control unit

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
TR01 Transfer of patent right

Effective date of registration: 20170323

Address after: Tokyo, Japan, Japan

Patentee after: Hitachi Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Hitachi Medical Corporation

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211117

Address after: Chiba County, Japan

Patentee after: Fujifilm medical health Co.,Ltd.

Address before: Tokyo, Japan

Patentee before: Hitachi, Ltd.

TR01 Transfer of patent right