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

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

Info

Publication number
CN106061397A
CN106061397A CN201480076573.0A CN201480076573A CN106061397A CN 106061397 A CN106061397 A CN 106061397A CN 201480076573 A CN201480076573 A CN 201480076573A CN 106061397 A CN106061397 A CN 106061397A
Authority
CN
China
Prior art keywords
waveform
candidate areas
correlation
stabilisation
heart rate
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.)
Pending
Application number
CN201480076573.0A
Other languages
English (en)
Inventor
中村雅志
村下贤
坂下肇
笠原英司
松下典义
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Ltd filed Critical Hitachi Ltd
Publication of CN106061397A publication Critical patent/CN106061397A/zh
Pending legal-status Critical Current

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/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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/02—Measuring pulse or heart rate
    • 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/0866—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
    • 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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
    • 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
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Gynecology & Obstetrics (AREA)
  • Pregnancy & Childbirth (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

基准帧选择部从表示胎儿的心脏的帧列中选择基准帧。候补区域组设定部对构成帧列的各个帧设定候补区域组。相关值运算部在基准帧与除此以外的各帧之间针对每个候补区域运算相关值。由此，生成与多个候补区域对应的多个相关值波形。稳定化波形部分确定部针对每个相关值波形确定稳定化波形部分。稳定化区域确定部在多个稳定化波形部分中确定稳定化度最高的稳定化波形部分(即，稳定化最高的候补区域)。心率运算部基于稳定化度最高的稳定化波形部分运算胎儿的心跳信息(心率等)。

Description

超声波诊断装置以及超声波图像处理方法
技术领域
本发明涉及超声波诊断装置，特别涉及获得进行周期性运动的脏器的周期信息的超声波诊断装置。
背景技术
针对胎儿的心脏，难以利用心电图记录仪等直接计测心率等。另一方面，通过利用超声波诊断装置能够获得心跳等的信息。
例如在专利文献1所公开的超声波诊断装置中，以表示胎儿的心脏的多个断层图像为对象，进行基准断层图像与除此以外的各断层图像之间的相关运算。根据表示该运算结果的相关值波形运算胎儿的心率。
另外，在专利文献2所公开的超声波诊断装置中，基于超声波图像分析胎儿的身体的动作以及心脏的动作，由此获得表示身体的变动的波形与表示心脏的运动的波形。基于减去表示身体的变动的波形的心脏的运动的波形，运算胎儿的心率。
现有技术文献
专利文献
专利文献1：日本特开2013-198635号公报
专利文献2：日本特开2013-198636号公报
发明内容
发明所要解决的课题
然而，在根据相关值波形运算心跳信息的情况下，如何设定成为相关值运算的对象的关心区域，对心跳信息的运算精度影响较大。例如，若在心脏中不稳定地进行周期运动的部分设定关心区域，则无法获得稳定的相关值波形，产生心跳信息的测定精度降低的问题。期望以适当的位置或者适当的尺寸来设定关心区域。
特别地，胎儿的心脏极小，容易移动。另外，超声波图像所表现的心脏的边界不清晰的情况较多。因此，用户手动地指定心脏中的稳定地进行周期运动的部位极其困难。
本发明的目的在于在超声波诊断装置中，针对进行周期性运动的脏器，提高周期信息的测定精度。或者，本发明的目的在于当在胎儿的心脏的剖面上设定心跳信息计测用的关心区域的情况下，减轻或者消除用户的负担。或者，本发明的目的在于能够使设定于胎儿的心脏的剖面上的心跳信息计测用的关心区域的位置以及尺寸的至少一方最佳化。
用于解决课题的手段
本发明的超声波诊断装置的特征在于，具有：帧列生成部，其基于通过对周期性地进行运动的脏器发送接收超声波而获得的信号生成帧列；候补区域组设定部，其针对上述帧列的每一个帧设定候补区域组；相关值运算部，其针对每个上述候补区域，在上述帧列中的基准帧与除此以外的各帧之间依次运算相关值，由此针对每个上述候补区域生成表示上述相关值的时间变化的相关值波形；稳定化波形部分确定部，其在每个上述候补区域的相关值波形中确定稳定化波形部分；最佳稳定化波形部分确定部，其从通过上述相关值运算部生成的多个相关值波形中被确定的多个稳定化波形部分中确定最佳稳定化波形部分；以及周期信息运算部，其基于从与上述最佳稳定化波形部分对应的候补区域获得的相关值波形，运算上述脏器的运动的周期信息。
根据上述的构成，生成与多个候补区域对应的多个相关值波形，针对各个相关值波形，确定其中的稳定化波形部分。即，各个相关值波形不整个成为评价对象，其中的稳定化波形部分成为评价对象。在该情况下，例如也可以将相关值波形上连续存在的偏差较少的部分确定为稳定化波形部分。或者，也可以将处于在相关值波形上离散地存在的偏差较少的关系的多个波形片段的集合确定为稳定化波形部分。若确定与多个候补区域对应的多个稳定化波形部分，则从其中确定最佳稳定化波形部分。该确定相当于多个候补区域中的稳定化区域(周期信息计测用的关心区域)的确定。因此，根据最佳稳定化波形部分，或者将其包含的相关值波形运算周期信息。若成为对象的脏器为心脏，则作为周期信息运算心率等心跳信息。
上述结构预先准备成为关心区域的候补的多个候补区域，评价针对这些候补区域而运算的多个相关值波形，从而选择最佳候补区域(或者所参照的波形部分)。因此，在评价相关值波形后，决定关心区域，从而提高关心区域的设定精度。另外，能够消除用户必须边预测或者考虑稳定性边设定关心区域的繁琐的问题。
根据经验，相关值波形的整体稳定的情况几乎不存在，多数情况下，在各个相关值波形中包含稳定的部分与不稳定的部分。特别是在胎儿的心脏的计测时，认为上述的趋势较强。根据本发明，在相关值波形的评价时，能够除去稳定化波形部分以外的不稳定化波形部分(例如值过于极端的部分)而评价相关值波形。因此，能够积极地利用有用或者优良的波形信息。与最佳稳定化波形部分对应的候补区域与脏器中的稳定地进行周期运动的部分对应。因此，根据本发明，能够从这样的稳定地进行周期运动的部分精度良好地测定周期信息。
优选，上述候补区域组由在帧区域的整体内或者一部分内相互具有非一致的关系而设定的多个候补区域构成。由此，能够确定适于运算周期信息的候补区域。帧列由在时间轴上排列的多个帧构成。各帧相当于组织中的计测对象剖面，具体而言，与波束扫描面或者断层图像对应。针对其全体设定候补区域组，或者针对其一部分设定候补区域组。优选预先准备具有相互不同的图案的多个候补区域组，手动或者与诊断部位等对应地自动设定任意的候补区域组。
优选，上述候补区域组包含在上述帧区域的整体内或者一部分内至少设定在不同的位置的多个候补区域。由此，能够使运算周期信息的区域的位置最佳化。
优选，上述候补区域组包含在上述帧区域的整体内或者一部分内至少具有不同的尺寸的多个候补区域。由此，能够使运算周期信息的区域的尺寸最佳化。
优选，上述稳定化波形部分确定部通过上述相关值波形的波形分析来确定上述稳定化波形部分。
优选，上述稳定化波形部分确定部包含：生成部，其在上述相关值波形中针对每个邻接峰值间隔运算伪周期信息，由此生成伪周期信息列；以及判定部，其在上述伪周期信息列中判定满足稳定化条件的多个伪周期信息，由此确定上述稳定化波形部分。由此，能够除去过于极端的伪周期信息来评价相关值波形，因此能够不受过于极端的伪周期信息的影响，而确定最佳稳定化波形部分。
优选，上述判定部包含：排序部，其根据排序条件对上述伪周期信息列进行排序；以及确定部，其从上述排序后的伪周期信息列中确定以排序方向排列的基准数的伪周期信息，作为上述多个伪周期信息。
优选，上述排序部按值从大到小的顺序或者从小到大的顺序对上述伪周期信息列进行排序，上述确定部将上述排序后的伪周期信息列的中间部分确定为上述基准数的伪周期信息。与中间部分相比，在除此以外的部分包含有过于极端的伪周期信息，在中间部分包含有与除此以外的部分相比稳定的伪周期信息。因此，通过将与中间部分所包含的伪周期信息对应的波形部分确定为稳定化波形部分，能够除去过于极端的伪周期信息来评价相关值波形。
优选，上述判定部包含：针对上述伪周期信息列设定多个偏差参照窗并运算多个偏差的功能；以及从上述多个偏差中确定最小的偏差，由此确定上述相关值波形中的上述稳定化波形部分的功能。在偏差成为最小的参照窗中，与除此以外的参照窗相比，包含有稳定的伪周期信息。根据该结构，能够除去过于极端的伪周期信息来评价相关值波形。
优选，上述最佳稳定化波形部分确定部将在上述多个稳定化波形部分中偏差成为最小的稳定化波形部分确定为上述最佳稳定化波形部分。在与偏差成为最小的最佳稳定化波形部分对应的候补区域中，与其他的候补区域相比周期运动稳定。因此，基于从该候补区域获得的相关值波形求得周期信息，从而周期信息的测定精度提高。
优选，上述周期信息运算部根据上述最佳稳定化波形部分运算上述周期信息。最佳稳定化波形部分与其他的波形部分相比稳定(过于极端的部分被除去)。因此，从最佳稳定化波形部分求得周期信息，从而周期信息的测定精度进一步提高。
另外，本发明的超声波图像处理方法的特征在于，包含：接受基于对周期性地进行运动的脏器发送接收超声波而获得的信号而生成的帧列，并针对上述帧列的每一个帧设定候补区域组的工序；针对每个上述候补区域，在上述帧列中的基准帧与除此以外的各帧之间依次运算相关值，由此针对每个上述候补区域生成表示上述相关值的时间变化的相关值波形的工序；针对每个上述候补区域，在相关值波形中确定稳定化波形部分的工序；从在生成的多个上述相关值波形中被确定的多个稳定化波形部分中确定最佳稳定化波形部分的工序；以及基于从与上述最佳稳定化波形部分对应的候补区域获得的相关值波形运算上述脏器的运动的周期信息的工序。
发明效果
根据本发明，能够在超声波诊断装置中提高进行周期性运动的脏器的周期信息的测定精度。
附图说明
图1是表示本发明的实施方式的超声波诊断装置的一个例子的框图。
图2是表示候补区域组的设定例的示意图。
图3是表示各候补区域中的相关值波形的一个例子的图。
图4是表示实施例1的处理的流程图。
图5A是用于对实施例1的处理进行说明的图。
图5B是用于对实施例1的处理进行说明的图。
图6是表示实施例2的处理的流程图。
图7A是用于对实施例2的处理进行说明的图。
图7B是用于对实施例2的处理进行说明的图。
图7C是用于对实施例2的处理进行说明的图。
图7D是用于对实施例2的处理进行说明的图。
图8A是表示变形例1的候补区域组的设定例的示意图。
图8B是表示变形例1的候补区域组的设定例的示意图。
图8C是表示变形例1的候补区域组的设定例的示意图。
图8D是表示变形例1的候补区域组的设定例的示意图。
图8E是表示变形例1的候补区域组的设定例的示意图。
图8F是表示变形例1的候补区域组的设定例的示意图。
图8G是表示变形例1的候补区域组的设定例的示意图。
图8H是表示变形例1的候补区域组的设定例的示意图。
图8I是表示变形例1的候补区域组的设定例的示意图。
图9A是表示变形例2的候补区域组的设定例的示意图。
图9B是表示变形例2的候补区域组的设定例的示意图。
图9C是表示变形例2的候补区域组的设定例的示意图。
图9D是表示变形例2的候补区域组的设定例的示意图。
图9E是表示变形例2的候补区域组的设定例的示意图。
图9F是表示变形例2的候补区域组的设定例的示意图。
图10是表示变形例3的候补区域组的设定例的示意图。
图11A是表示变形例3的候补区域组的设定例的示意图。
图11B是表示变形例3的候补区域组的设定例的示意图。
图11C是表示变形例3的候补区域组的设定例的示意图。
图11D是表示变形例3的候补区域组的设定例的示意图。
图11E是表示变形例3的候补区域组的设定例的示意图。
图11F是表示变形例3的候补区域组的设定例的示意图。
图11G是表示变形例3的候补区域组的设定例的示意图。
图11H是表示变形例3的候补区域组的设定例的示意图。
图11I是表示变形例3的候补区域组的设定例的示意图。
图11J是表示变形例3的候补区域组的设定例的示意图。
图11K是表示变形例3的候补区域组的设定例的示意图。
图11L是表示变形例3的候补区域组的设定例的示意图。
具体实施方式
图1表示本发明的实施方式的超声波诊断装置的一个例子。超声波诊断装置是设置于医院等医疗机构，通过对人体的超声波的发送接收而形成超声波图像的装置。本实施方式的超声波诊断装置如以下详述的那样，具备通过对胎儿发送接收超声波，计测胎儿的心跳信息的功能。周期性地进行运动的其他的组织也可以成为计测对象。
在图1中，探头10是对包含对象物的诊断区域发送接收超声波的收发波器。探头10具备发送接收超声波的多个振动元件。由多个振动元件形成超声波束。超声波束被反复电子扫描。由此，依次形成波束扫描面。作为电子扫描方式，公知有电子扇形扫描方式、电子线性扫描方式等。
收发部12在发送时，对探头10具备的多个振动元件输出被延迟处理后的多个发送信号。由此，从多个振动元件向生物体内输送发送波束。在接收时，若来自生物体内的反射波被多个振动元件接受，则从上述多个振动元件向收发部12输出多个接收信号。在收发部12对多个接收信号实施整相加法处理等，从而形成接收波束。即，收发部12输出整相加法处理后的接收信号(波束数据)。通过收发部12的作用，发送波束以及接收波束(两者合并称为超声波束)被电子地扫描。由此，构成上述的波束扫描面。波束扫描面相当于多个波束数据，它们构成接收帧(接收帧数据)。此外，各波束数据由在深度方向排列的多个回波数据构成。通过重复超声波束的电子扫描，从收发部12输出在时间轴上排列的多个接收帧。它们构成接收帧列。从收发部12输出的波束数据经由未图示的信号处理部后发送至图像形成部14。信号处理部具备检波电路、对数压缩电路等。此外，在超声波的收发时，也可以利用发送孔径综合等技术。
图像形成部14由具有坐标变换功能以及插值处理功能等的数字扫描转换器构成。图像形成部14基于接收帧列，形成由多个显示帧构成的显示帧列100。构成显示帧列100的各个显示帧是B模式断层图像的数据。显示帧列100被输出到监视器等显示部34来显示。由此，能够实时地将B模式断层图像显示为动态图像。在本实施方式中，显示帧列100被存储于帧列存储部18。
图像处理部16包含：帧列存储部18、基准帧选择部20、候补区域组设定部22、相关值运算部24、稳定化波形部分确定部26、稳定化区域确定部28以及心率运算部30。
基准帧选择部20从存储于帧列存储部18的显示帧列100中，选择成为相关值运算的基准的基准帧。基准帧选择部20例如与经由操作部32被输入的用户操作对应地选择基准帧。例如，通过显示部34显示存储于帧列存储部18的显示帧列100。用户边观察显示于显示部34的显示帧列100，边使用操作部32指定基准帧。此外，基准帧选择部20也可以从显示帧列100中将任意的显示帧选择为基准帧。也可以使基准帧的选择自动化。例如，也可以通过图像分析来确定基准帧。
候补区域组设定部22针对成为处理对象的显示帧列的每一个帧设定候补区域组110。例如，候补区域组110由具有非一致的关系而分散地设定的多个候补区域构成。具体而言，候补区域组设定部22针对显示帧列的每一个，在相互不同的位置设定多个候补区域。另外，候补区域组设定部22也可以设定尺寸相互不同的多个候补区域。在本实施方式中，候补区域组设定部22在显示帧列的每一个中，对胎儿的心脏设定候补区域组110。候补区域组设定部22例如与经由操作部32被输入的用户操作对应地设定候补区域组110。例如，通过显示部34显示基准帧。用户边观察显示于显示部34的基准帧，边使用操作部32指定候补区域组110的设定位置。针对该被指定的位置设定候补区域组110。候补区域组设定部22针对构成处理对象的显示帧列的各个显示帧，在与对基准帧设定的位置相同的位置设定候补区域组110。此外，候补区域组设定部22也可以对基准帧进行图像分析，在胎儿的心脏的区域设定候补区域组。候补区域组设定部22从帧列存储部18读出与各个候补区域对应的显示帧列120并将其输出至相关值运算部24。
图2表示候补区域组的设定例。在基准帧40中表现胎儿的身体42以及胎儿的心脏44。在图2所示的例子中设定了六个矩形状的候补区域(候补区域50A～50F)。候补区域50A～50E以局部包含心脏44的方式分别设定于不同的位置。候补区域50F被设定为包含心脏44的整体。候补区域50A～50E的尺寸相同。候补区域50A～50D相互不重叠地设定于将候补区域50F四等分而得的区域。候补区域50E被设定为与候补区域50A～50D局部重叠。在图2所示的例子中，候补区域50A～50F为矩形，但也可以为其他的多边形、圆形、椭圆形。另外，候补区域50A～50E可以为相同的大小，也可以为不同的大小。候补区域50A～50D也可以相互局部重叠地设定。另外，候补区域的个数也不限定于图2所示的例子，只要设定多个候补区域即可。候补区域的形状、尺寸、个数以及设定位置是任意的，例如也可以通过用户对操作部32的操作将它们指定。
返回图1进行说明。相关值运算部24针对每个候补区域，在基准帧与基准帧以外的各显示帧之间依次运算相关值。由此，针对每个候补区域生成表示相关值的时间变化的相关值波形130。列举具体例进行说明。假定显示帧F1、F2、F3、F4是作为处理对象的显示帧。在该情况下，相关值运算部24针对每个候补区域，运算基准帧(例如显示帧F1)与显示帧F2之间的相关值、基准帧与显示帧F3之间的相关值以及基准帧与显示帧F4之间的相关值。由此，针对每个候补区域获得表示相关值的时间变化的相关值波形。作为相关值，例如能够利用SSD(Sum of Square Difference：差的平方和)、SAD(Sum of Absolute Difference：差的绝对值的和)或者平均值的差等公知的方法。
图3表示与候补区域50A～50F对应的相关值波形的一个例子。图3中的横轴表示时间轴，纵轴表示相关值。相关值波形A是表示图2所示的候补区域50A中的相关值的时间变化的波形。相关值波形B是表示候补区域50B中的相关值的时间变化的波形。相关值波形C是表示候补区域50C中的相关值的时间变化的波形。相关值波形D是表示候补区域50D中的相关值的时间变化的波形。相关值波形E是表示候补区域50E中的相关值的时间变化的波形。相关值波形F是表示候补区域50F中的相关值的时间变化的波形。
若返回图1进行说明，则稳定化波形部分确定部26在每个候补区域的相关值波形130中确定稳定化波形部分140。即，稳定化波形部分确定部26针对每个候补区域，从相关值波形130除去值过于极端的部分，确定波形稳定的稳定化波形部分140。例如，稳定化波形部分确定部26针对每个候补区域，基于相关值波形130运算心脏的伪心率，基于伪心率从相关值波形130确定稳定化波形部分140。
稳定化区域确定部28将在多个相关值波形130中确定的多个稳定化波形部分140相互比较，从而从多个稳定化波形部分140中确定最佳稳定化波形部分。然后，稳定化区域确定部28将与最佳稳定化波形部分对应的候补区域确定为稳定化区域。例如，稳定化区域确定部28针对每个候补区域，运算与稳定化波形部分140对应的伪心率的偏差，确定伪心率的偏差为最小的稳定化波形部分(最佳稳定化波形部分)，将与该最佳稳定化波形部分对应的候补区域确定为稳定化区域。
心率运算部30基于从稳定化区域获得的相关值波形，运算胎儿的心跳信息。心跳信息例如为心率。
此外，图1所示的探头10以外的结构能够利用例如处理器、电路等硬件资源而实现，在该实现中也可以根据需要利用存储器等器件。另外，探头10以外的结构例如也可以通过计算机实现。换句话说，也可以通过计算机具备的CPU、存储器、硬盘等硬件资源与规定CPU等的动作的软件(程序)的协作，实现探头10以外的结构的全部或者一部分。该程序经由CD、DVD等记录介质，或者经由网络等通信路径，存储于未图示的存储装置。作为其他的例子，探头10以外的结构也可以通过DSP(Digital Signal Processor数字信号处理器)、FPGA(Field Programmable Gate Array现场可编程逻辑门阵列)等实现。
接下来，参照图4所示的流程图，对本实施方式的超声波诊断装置的处理的实施例1进行说明。首先，通过显示部34显示存储于帧列存储部18的显示帧列。用户使用操作部32从显示帧列中指定处理对象帧列(S01)。进而，用户使用操作部32从处理对象帧列中指定基准帧(S02)。接着，通过显示部34显示基准帧。用户边观察该基准帧，边使用操作部32指定候补区域组的设定位置。由此，通过候补区域组设定部22对处理对象帧列的每一个帧设定候补区域组(S03)。作为一个例子，如图2所示，候补区域组设定部22对处理对象帧列的每一个帧设定候补区域50A～50F。若设定候补区域组，则相关值运算部24针对每个候补区域在基准帧与各显示帧之间依次运算相关值，从而针对每个候补区域生成相关值波形(S04)。作为一个例子，如图3所示，相关值运算部24针对候补区域50A～50F生成相关值波形A～F。
而且，稳定化波形部分确定部26针对每个相关值波形依次运算伪心率(S05)。具体而言，稳定化波形部分确定部26针对每个候补区域，不断地探索相关值波形的峰值点(极大点或者极小点)，将相互邻接的峰值点(极大点或者极小点)的时间间隔分别运算为伪的一次心跳时间。而且，稳定化波形部分确定部26针对每个候补区域，基于多个伪的一次心跳时间，运算每单位时间的多个伪心率(bmp)。由此，运算在时间轴上排列的多个伪心率，它们构成伪心率列。若参照图3进行说明，则稳定化波形部分确定部26针对相关值波形A，将相互邻接的峰值点(例如极大值)的时间间隔T1～T9分别运算为伪的一次心跳时间。而且，稳定化波形部分确定部26根据时间间隔T1～T9运算每单位时间的伪心率R1～R9。在时间轴上排列的伪心率R1～R9构成伪心率列。稳定化波形部分确定部26针对相关值波形B～F也运算伪心率列。
接下来，稳定化波形部分确定部26针对每个候补区域，按值从大到小的顺序排序处理伪心率列(重排)(S06)。若以伪心率R1～R9为例进行说明，则如图5A所示，稳定化波形部分确定部26按值从大到小的顺序排序伪心率R1～R9。或者，如图5B所示，稳定化波形部分确定部26也可以按值从小到大的顺序排序伪心率R1～R9。稳定化波形部分确定部26针对相关值波形B～F也排序伪心率列。
然后，稳定化波形部分确定部26针对每个候补区域，运算排序后的伪心率列中的中央N个(配置于中央的N个)的平均值以及偏差(S07)。N为整数。作为一个例子，若N＝5，则在图5A所示的例子中，稳定化波形部分确定部26运算配置于中央的五个伪心率(伪心率R9、R6、R5、R7、R4)的平均值以及偏差。或者，如图5B所示，稳定化波形部分确定部26也可以运算按值从小到大的顺序排序的伪心率R1～R9中的中央N个的平均值以及偏差。稳定化波形部分确定部26针对相关值波形B～F也运算排序后的伪心率列中的中央N个的平均值以及偏差。此外，在图5A以及图5B所示的例子中，N＝5，但也可以使用除此以外的值。
此处，对偏差进行说明。若将N个伪心率的每一个的值设为xi(i＝1～N)，将它们的平均设为m，则通过以下的式(1)求得方差。
[数学式1]
σ 2 = 1 N Σ i = 1 N ( x i - m ) 2 ... ( 1 )
将该方差的正的平方根σ称为标准偏差。
另外，将标准偏差σ除以平均值m所得的值称为变动系数CV(Coefficient ofVariation)。变动系数CV由以下的式(2)表示。
CV＝σ/m····(2)
变动系数CV表示不取决于平均值的相对偏差。例如，即便标准偏差σ为相同的“20”，在平均值为“50”的情况与平均值为“200”的情况下，也认为后者(平均值＝200)的偏差较少(在平均值为“50”的情况下，CV＝0.4，在平均值为“200”的情况下，CV＝0.1)。稳定化波形部分确定部26针对相关值波形A～F，运算中央N个伪心率的偏差CV。
而且，稳定化区域确定部28将相关值波形A～F的偏差CV相互比较，从而确定偏差CV成为最小的相关值波形，确定与该相关值波形对应的候补区域(稳定化区域)(S08)。即，稳定化区域确定部28使用伪心率列中的中央N个的偏差CV评价相关值波形A～F的稳定化度，确定稳定化度最高(偏差CV成为最小)的相关值波形。作为一个例子，在相关值波形A的偏差CV在相关值波形A～F中为最小的情况下，稳定化区域确定部28将与相关值波形A对应的候补区域50A确定为稳定化区域。
此外，稳定化波形部分确定部26也可以针对每个相关值波形运算伪心率列中的中央N个的标准偏差σ，确定标准偏差σ为最小的相关值波形，将与该相关值波形对应的候补区域确定为稳定化区域。
若如以上那样确定稳定化区域，则心率运算部30运算稳定化区域的心率(S09)。例如，心率运算部30将稳定化区域的伪心率列(全部伪心率)的平均值运算为胎儿的心率(bpm)。心率运算部30也可以将排序后的伪心率列中的中央N个的平均值运算为胎儿的心率。胎儿的心率例如被输出到显示部34来显示。作为一个例子，在候补区域50A被确定为稳定化区域的情况下，心率运算部30将排序后的伪心率列中的中央N个的平均值(例如图5A所示的伪心率R9、R6、R5、R7、R4的平均值)运算为胎儿的心率。
在继续处理的情况下(S10，是)，更新作为处理对象的显示帧列(S11)，以更新后的显示帧列为对象，进行步骤S04～S09的处理。例如，若用户使用操作部32将在其他的时间段取得的显示帧列指定为处理对象，则以被指定的显示帧列为对象，进行步骤S04～S09的处理。在不继续处理的情况下(S10，否)，心率的计测结束。
如以上那样，在实施例1中，针对每个候补区域，按值从大到小的顺序或者从小到大的顺序排序伪心率列，运算排序后的伪心率列中的中央N个的偏差CV。而且，基于每个候补区域的偏差CV评价各候补区域的相关值波形。由此，能够从相关值波形除去稳定化波形部分以外的不稳定化波形部分(心率过于极端的部分)来评价相关值波形。其结果，能够确定得到不受不稳定化波形部分的影响，与其他的候补区域相比能够获得稳定的相关值波形的稳定化区域。
对该点详细地进行说明。异常值(过于极端的心率)不可避免地包含于相关值波形中。因此，若基于伪心率列所包含的全部伪心率的偏差CV评价相关值波形，则导致也包含异常值地进行评价。在该情况下，评价的精度降低。与此相对，在实施例1中，按值从大到小的顺序或者从小到大的顺序排序伪心率列，运算排序后的伪心率列中的中央N个的偏差CV。而且，基于该偏差CV评价相关值波形。在进行了排序处理的情况下，在中央N个以外的范围中包含有过于极端的伪心率。在中央N个的范围中与其他的范围相比，不包含过于极端的伪心率。因此，与中央N个伪心率对应的波形部分，与同中央N个以外的伪心率对应的波形部分相比稳定，从而相当于相关值波形中的稳定化波形部分。因此，通过使用排序后的伪心率列中的中央N个的偏差CV，能够在除去异常值的状态下评价相关值波形，而确定稳定化区域。在该实施例1中，排序处理了伪心率列，因此稳定化波形部分是在时间轴上未必连续的多个波形部分的集合。
稳定化区域的周期运动与其他候补区域的周期运动相比稳定。因此，基于从稳定化区域获得的相关值波形运算心率，从而心率的测定精度提高。另外，排序后的中央N个伪心率与稳定化波形部分对应。因此，根据该稳定化波形部分运算心率，从而心率的测定精度进一步提高。
接下来，参照图6所示的流程图，对本实施方式的超声波诊断装置的处理的实施例2进行说明。在实施例2中不进行伪心率列的排序处理，而运算在伪心率列中按时间顺序排列的连续N个的偏差CV。然后，基于该偏差CV确定稳定化区域。以下，对实施例2的处理进行详述。
首先，与实施例1相同地，通过用户从存储于帧列存储部18的显示帧列中指定处理对象帧列(S20)，从该处理对象帧列中指定基准帧(S21)。接着，通过候补区域组设定部22针对处理对象帧列的每一个帧设定候补区域组(S22)。然后，通过稳定化波形部分确定部26针对每个候补区域生成相关值波形(S23)，针对每个相关值波形运算伪心率列(S24)。伪心率列由在时间轴上排列的多个伪心率构成。作为一个例子，如图2所示，设定候补区域50A～50F，如图3所示，运算与候补区域50A～50F对应的相关值波形A～F以及伪心率列。
接下来，稳定化波形部分确定部26针对每个候补区域，对伪心率列设定闸门(时间窗)。在该闸门中包含有按时间顺序排列的连续N个的伪心率。而且，稳定化波形部分确定部26边使闸门在时间方向错开，边运算各闸门所包含的连续N个伪心率的平均值以及偏差CV(S25)。换句话说，稳定化波形部分确定部26针对每个候补区域运算伪心率列的移动平均以及偏差CV。然后，稳定化波形部分确定部26针对每个候补区域确定偏差CV为最小的最佳闸门(S26)。稳定化波形部分确定部26将最佳闸门所包含的连续N个伪心率的平均值以及偏差CV输出至稳定化区域确定部28。
参照图7A、图7B、图7C以及图7D，对步骤S25、S26的处理的具体例进行说明。作为一个例子，以根据图3所示的相关值波形A求得的伪心率R1～R9为例进行说明。例如若N＝5，则如图7A所示，稳定化波形部分确定部26针对按时间顺序排列的伪心率R1～R5设定闸门，运算伪心率R1～R5的平均值以及偏差CV。接着，如图7B所示，稳定化波形部分确定部26使闸门错开而对伪心率R2～R6设定闸门，运算伪心率R2～R6的平均值以及偏差CV。另外，稳定化波形部分确定部26如图7C所示，运算伪心率R3～R7的平均值以及偏差CV，如图7D所示，运算伪心率R4～R8的平均值以及偏差CV。以下也相同，稳定化波形部分确定部26求得伪心率列的移动平均以及偏差CV。然后，稳定化波形部分确定部26在相关值波形A中确定在多个闸门中偏差CV成为最小的最佳闸门，将最佳闸门所包含的连续N个伪心率的平均值以及偏差CV输出至稳定化区域确定部28。例如，在相关值波形A中，在图7A所示的伪心率R1～R5的偏差CV成为最小的情况下，稳定化波形部分确定部26将伪心率R1～R5的平均值以及偏差CV输出至稳定化区域确定部28。
稳定化波形部分确定部26针对相关值波形A～F的每一个确定偏差CV成为最小的最佳闸门，将最佳闸门所包含的连续N个伪心率的平均值以及偏差CV输出至稳定化区域确定部28。此外，在图7A～7D所示的例子中，N＝5，但也可以使用除此以外的值。
然后，稳定化区域确定部28在相关值波形A～F中，确定最佳闸门所包含的连续N个伪心率的偏差CV成为最小的相关值波形，确定与该相关值波形对应的候补区域(稳定化区域)(S27)。即，稳定化区域确定部28使用伪心率列中的连续N个的偏差CV，评价相关值波形A～F的稳定化度，确定稳定化度最高的(偏差CV成为最小的)相关值波形。作为一个例子，在相关值波形A的偏差CV在相关值波形A～F中成为最小的情况下，稳定化区域确定部28将与相关值波形A对应的候补区域50A确定为稳定化区域。
此外，稳定化波形部分确定部26也可以针对各相关值波形的各闸门运算标准偏差σ，确定与标准偏差σ成为最小的闸门对应的相关值波形，将与该相关值波形对应的候补区域确定为稳定化区域。
若如以上那样确定稳定化区域，则心率运算部30运算稳定化区域的心率(S28)。例如，心率运算部30将稳定化区域的伪心率列(全部伪心率)的平均值运算为胎儿的心率。心率运算部30也可以将稳定化区域的伪心率列中的最佳闸门所包含的连续N个的平均值运算为胎儿的心率。胎儿的心率例如被输出到显示部34来显示。作为一个例子，在候补区域50A被确定为稳定化区域的情况下，心率运算部30将相关值波形A中的最佳闸门所包含的连续N个的平均值(例如图7A所示的伪心率R1～R5的平均值)运算为胎儿的心率。
在继续处理的情况下(S29，是)，更新作为处理对象的显示帧列(S30)，以更新后的显示帧列为对象进行步骤S23～S28的处理。在不继续处理的情况下(S29，否)，心率的计测结束。
如以上那样，在实施例2中，针对每个候补区域确定伪心率列中的连续N个的偏差CV成为最小的最佳闸门。而且，基于每个候补区域的最佳闸门所包含的伪心率的偏差CV，评价各候补区域的相关值波形。由此，能够从相关值波形除去不稳定化波形部分来评价相关值波形。其结果，能够不受不稳定化波形部分的影响地确定稳定化区域。
对该点详细地进行说明。最佳闸门所包含的连续N个伪心率的偏差CV比其他闸门中的偏差CV小。换句话说，在最佳闸门中，与其他闸门相比不包含过于极端的伪心率。因此，与最佳闸门的伪心率对应的波形部分与同其他闸门的伪心率对应的波形部分相比稳定，相当于相关值波形的稳定化波形部分。因此，通过使用最佳闸门所包含的连续N个伪心率的偏差CV，能够在除去异常值的状态下评价相关值波形来确定稳定化区域。
而且，基于从稳定化区域获得的相关值波形运算心率，从而能够提高心率的测定精度。另外，最佳闸门所包含的连续N个伪心率与稳定化波形部分对应。根据该稳定化波形部分对心率进行运算，从而心率的测定精度进一步提高。
另外，根据本实施方式，在不同的位置设定多个候补区域，从而能够从候补区域组中确定适于运算心率的候补区域(稳定地进行周期运动的候补区域)的位置。另外，设定尺寸不同的多个候补区域，从而能够从候补区域组中确定适于运算心率的候补区域的尺寸。
另外，通过稳定化区域确定部28确定稳定化区域，因此能够消除用户进行的稳定化区域的指定的繁琐。
此外，也可以将实施例1、2组合。例如，也可以通过实施例1的处理确定稳定化区域，通过实施例2的处理运算心率。具体而言，与实施例1相同地，稳定化波形部分确定部26针对每个候补区域，按值从大到小的顺序或者从小到大的顺序排序伪心率列，运算伪心率列中的中央N个的偏差。稳定化区域确定部28将该偏差成为最小的候补区域确定为稳定化区域。心率运算部30针对稳定化区域的伪心率列设定闸门，边使闸门错开边运算各闸门所包含的N个伪心率的平均值以及偏差。然后，心率运算部30确定在多个闸门中偏差成为最小的最佳闸门，将最佳闸门所包含的N个伪心率的平均值运算为胎儿的心率。
另外，也可以通过实施例2的处理确定稳定化区域，通过实施例1的处理运算心率。具体而言，与实施例2相同地，稳定化波形部分确定部26针对每个候补区域确定最佳闸门。稳定化区域确定部28基于每个候补区域的最佳闸门所包含的伪心率的偏差确定稳定化区域。心率运算部30按值从大到小的顺序或者从小到大的顺序排序稳定化区域的伪心率列，将排序后的伪心率列中的中央N个的平均值运算为胎儿的心率。
如以上那样，即便在将实施例1、2组合的情况下，也基于从稳定化区域获得的相关值波形运算心率，因此心率的测定精度提高。
接下来，对变形例的候补区域组的设定例进行说明。图8A～图8I表示变形例1的候补区域组的设定例。在变形例1中，如图8A～图8I所示，在设定于处理对象的显示帧列的关心区域60内设定了形状以及尺寸相同的九个候补区域(矩形状的候补区域61～69)。候补区域61～69相互局部重叠地设定。另外，图9A～图9F表示变形例2的候补区域组的设定例。在变形例2中，如图9A～图9F所示，在设定于处理对象的显示帧列的关心区域70内设定形状相同的六个候补区域(矩形状的候补区域71～76)。候补区域71～75的大小相同。候补区域76被设定为比候补区域71～75大，包含关心区域70的整体。此外，候补区域的形状、尺寸、个数以及设定位置是任意的。它们不限定于图8A～图8I以及图9A～图9F所示的例子。
接下来，参照图10以及图11A～图11L，对变形例3的候补区域组的设定例进行说明。例如如图10所示，候补区域组设定部22以基准帧为对象进行边界自动提取处理、学习功能等的图像分析处理，从而自动地确定胎儿的心脏44的区域、心脏44内的组织(例如左室等)。作为一个例子，候补区域组设定部22在左室的区域中设定关心区域46，在关心区域46内设定候补区域组。例如，在设定矩形状的候补区域的情况下，如图11A～图11L所示，候补区域组设定部22在与关心区域46内接的区域80内设定候补区域组(候补区域81～92)。候补区域81～92的形状、尺寸以及设定位置是任意的。这样，通过自动地确定对象物并自动地设定候补区域，省去用户进行的候补区域的设定的麻烦。
在本实施方式中，利用伪心率确定稳定化波形部分以及稳定化区域，但也可以利用从相关值波形获得的多个伪的一次心跳时间，确定稳定化波形部分以及稳定化区域。
另外，在本实施方式中，利用数字扫描转换后的显示帧列确定稳定化波形部分以及稳定化区域，但也可以利用数字扫描转换前的接收帧列确定稳定化波形部分以及稳定化区域。在该情况下，将从收发部12输出的接收帧列存储于帧列存储部18。图像处理部16以接收帧列为对象执行处理，从而确定稳定化波形部分以及稳定化区域，运算胎儿的心率。
符号说明
10 探头，
12 收发部，
14 图像形成部，
16 图像处理部，
18 帧列存储部，
20 基准帧选择部，
22 候补区域组设定部，
24 相关值运算部，
26 稳定化波形部分确定部，
28 稳定化区域确定部，
30 心率运算部，
32 操作部，
34 显示部。

Claims (12)

1.一种超声波诊断装置，其特征在于，具有：
帧列生成部，其基于通过对周期性地进行运动的脏器发送接收超声波而获得的信号生成帧列；
候补区域组设定部，其针对所述帧列的每一个帧设定候补区域组；
相关值运算部，其针对每个所述候补区域，在所述帧列中的基准帧与除此以外的各帧之间依次运算相关值，由此针对每个所述候补区域生成表示所述相关值的时间变化的相关值波形；
稳定化波形部分确定部，其在每个所述候补区域的相关值波形中确定稳定化波形部分；
最佳稳定化波形部分确定部，其从通过所述相关值运算部生成的多个相关值波形中被确定的多个稳定化波形部分中确定最佳稳定化波形部分；以及
周期信息运算部，其基于从与所述最佳稳定化波形部分对应的候补区域获得的相关值波形，运算所述脏器的运动的周期信息。
2.根据权利要求1所述的超声波诊断装置，其特征在于，
所述候补区域组由在帧区域的整体内或者一部分内相互具有非一致的关系而设定的多个候补区域构成。
3.根据权利要求2所述的超声波诊断装置，其特征在于，
所述候补区域组包含在所述帧区域的整体内或者一部分内至少设定在不同的位置的多个候补区域。
4.根据权利要求2所述的超声波诊断装置，其特征在于，
所述候补区域组包含在所述帧区域的整体内或者一部分内至少具有不同的尺寸的多个候补区域。
5.根据权利要求1所述的超声波诊断装置，其特征在于，
所述稳定化波形部分确定部通过所述相关值波形的波形分析来确定所述稳定化波形部分。
6.根据权利要求5所述的超声波诊断装置，其特征在于，
所述稳定化波形部分确定部包含：
生成部，其在所述相关值波形中针对每个邻接峰值间隔运算伪周期信息，由此生成伪周期信息列；以及
判定部，其在所述伪周期信息列中判定满足稳定化条件的多个伪周期信息，由此确定所述稳定化波形部分。
7.根据权利要求6所述的超声波诊断装置，其特征在于，
所述判定部包含：
排序部，其根据排序条件对所述伪周期信息列进行排序；以及
确定部，其从所述排序后的伪周期信息列中确定以排序方向排列的基准数的伪周期信息，作为所述多个伪周期信息。
8.根据权利要求7所述的超声波诊断装置，其特征在于，
所述排序部按值从大到小的顺序或者从小到大的顺序对所述伪周期信息列进行排序，
所述确定部将所述排序后的伪周期信息列的中间部分确定为所述基准数的伪周期信息。
9.根据权利要求6所述的超声波诊断装置，其特征在于，
所述判定部包含：
针对所述伪周期信息列设定多个偏差参照窗并运算多个偏差的功能；以及
从所述多个偏差中确定最小的偏差，由此确定所述相关值波形中的所述稳定化波形部分的功能。
10.根据权利要求1所述的超声波诊断装置，其特征在于，
所述最佳稳定化波形部分确定部将在所述多个稳定化波形部分中偏差为最小的稳定化波形部分确定为所述最佳稳定化波形部分。
11.根据权利要求1所述的超声波诊断装置，其特征在于，
所述周期信息运算部根据所述最佳稳定化波形部分运算所述周期信息。
12.一种超声波图像处理方法，其特征在于，包含：
接受基于通过对周期性地进行运动的脏器发送接收超声波而获得的信号而生成的帧列，并针对所述帧列的每一个帧设定候补区域组的工序；
针对每个所述候补区域，在所述帧列中的基准帧与除此以外的各帧之间依次运算相关值，由此针对每个所述候补区域生成表示所述相关值的时间变化的相关值波形的工序；
针对每个所述候补区域，在相关值波形中确定稳定化波形部分的工序；
从在生成的多个所述相关值波形中被确定的多个稳定化波形部分中确定最佳稳定化波形部分的工序；以及
基于从与所述最佳稳定化波形部分对应的候补区域获得的相关值波形运算所述脏器的运动的周期信息的工序。
CN201480076573.0A 2014-02-27 2014-10-08 超声波诊断装置以及超声波图像处理方法 Pending CN106061397A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2014-037116 2014-02-27
JP2014037116A JP5651258B1 (ja) 2014-02-27 2014-02-27 超音波診断装置及びプログラム
PCT/JP2014/076940 WO2015129090A1 (ja) 2014-02-27 2014-10-08 超音波診断装置及び超音波画像処理方法

Publications (1)

Publication Number Publication Date
CN106061397A true CN106061397A (zh) 2016-10-26

Family

ID=52344877

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480076573.0A Pending CN106061397A (zh) 2014-02-27 2014-10-08 超声波诊断装置以及超声波图像处理方法

Country Status (4)

Country Link
US (1) US20170065257A1 (zh)
JP (1) JP5651258B1 (zh)
CN (1) CN106061397A (zh)
WO (1) WO2015129090A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111374706A (zh) * 2018-12-28 2020-07-07 深圳迈瑞生物医疗电子股份有限公司 一种胎儿心率显示方法及超声成像装置、存储介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3508132A1 (en) * 2018-01-04 2019-07-10 Koninklijke Philips N.V. Ultrasound system and method for correcting motion-induced misalignment in image fusion

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5944241A (ja) * 1982-09-07 1984-03-12 穂垣　正暢 胎児心拍測定装置
JP2000225115A (ja) * 1999-02-05 2000-08-15 Shimadzu Corp 超音波診断装置
CN1988850A (zh) * 2004-07-30 2007-06-27 株式会社日立医药 医用图像诊断辅助方法、装置和图像处理程序
JP2010233966A (ja) * 2009-03-31 2010-10-21 Toshiba Corp 超音波診断装置及び擬似心拍音出力用制御プログラム
CN102596050A (zh) * 2009-10-27 2012-07-18 株式会社日立医疗器械 超声波成像装置、超声波成像方法及超声波成像用程序
CN102939050A (zh) * 2010-06-04 2013-02-20 株式会社日立医疗器械 超声波诊断装置以及超声波收发方法
CN103281963A (zh) * 2010-12-27 2013-09-04 株式会社日立医疗器械 超声波诊断装置以及图像处理方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8298143B2 (en) * 2005-05-30 2012-10-30 Panasonic Corporation Ultrasonograph that determines tissue properties utilizing a reference waveform
RU2012150806A (ru) * 2010-04-28 2014-06-10 Конинклейке Филипс Электроникс Н.В. Визуализация размера инфаркта миокарда в экг-диагностике
JP5597492B2 (ja) * 2010-09-08 2014-10-01 株式会社東芝 超音波診断装置、画像処理装置およびプログラム
KR101511084B1 (ko) * 2012-10-11 2015-04-10 삼성메디슨 주식회사 의료 화상 표시 방법, 장치 및 사용자 인터페이스 화면 생성 방법
CN104066382B (zh) * 2012-01-12 2016-11-16 株式会社日立制作所 图像诊断装置以及图像显示方法
JP5386001B2 (ja) * 2012-03-26 2014-01-15 雅彦 中田 超音波診断装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5944241A (ja) * 1982-09-07 1984-03-12 穂垣　正暢 胎児心拍測定装置
JP2000225115A (ja) * 1999-02-05 2000-08-15 Shimadzu Corp 超音波診断装置
CN1988850A (zh) * 2004-07-30 2007-06-27 株式会社日立医药 医用图像诊断辅助方法、装置和图像处理程序
JP2010233966A (ja) * 2009-03-31 2010-10-21 Toshiba Corp 超音波診断装置及び擬似心拍音出力用制御プログラム
CN102596050A (zh) * 2009-10-27 2012-07-18 株式会社日立医疗器械 超声波成像装置、超声波成像方法及超声波成像用程序
CN102939050A (zh) * 2010-06-04 2013-02-20 株式会社日立医疗器械 超声波诊断装置以及超声波收发方法
CN103281963A (zh) * 2010-12-27 2013-09-04 株式会社日立医疗器械 超声波诊断装置以及图像处理方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111374706A (zh) * 2018-12-28 2020-07-07 深圳迈瑞生物医疗电子股份有限公司 一种胎儿心率显示方法及超声成像装置、存储介质

Also Published As

Publication number Publication date
JP2015159980A (ja) 2015-09-07
US20170065257A1 (en) 2017-03-09
JP5651258B1 (ja) 2015-01-07
WO2015129090A1 (ja) 2015-09-03

Similar Documents

Publication Publication Date Title
US10682118B2 (en) Ultrasound system and method for analyzing cardiac periodicity
CN104114102B (zh) 超声波诊断装置、图像处理装置以及图像处理方法
JP4966108B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
US20140031690A1 (en) Ultrasound diagnostic apparatus and method for identifying blood vessel
CN1915178B (zh) 超声波诊断装置和超声波图像处理方法
US9675320B2 (en) Diagnostic ultrasound apparatus
CN107106120A (zh) 用于通过超声换能器的持续振动进行超声弹性成像的方法
CN101317773A (zh) 超声波图像处理装置
CN104080408A (zh) 超声波诊断装置、图像处理装置以及图像处理方法
JP5920786B2 (ja) 超音波診断装置と心機能検査区間の検索表示方法
US11883242B2 (en) System and method for scanning for a second object within a first object using an adaptive scheduler
US20160174937A1 (en) Wireless ultrasound probe
US9192360B2 (en) Ultrasound system and method of obtaining ultrasound image
US9629608B2 (en) Ultrasound diagnosis apparatus and control method
US11944485B2 (en) Ultrasound device, systems, and methods for lung pulse detection by plueral line movement
CN106061397A (zh) 超声波诊断装置以及超声波图像处理方法
CN108135573A (zh) 超声波诊断***及超声波诊断方法
EP2659839A1 (en) Ultrasonic diagnosis device and image processing method
JPH11221210A (ja) 超音波診断装置
CN106659470B (zh) 超声波诊断装置
JP5443781B2 (ja) 超音波診断装置
JP2017012607A (ja) 医用画像診断装置、医用画像処理装置及びプログラム
US20140107485A1 (en) Method and apparatus for generating volume image
JP2008061835A (ja) 超音波診断装置及び超音波診断装置の制御プログラム
CN110037734A (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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20161026

WD01 Invention patent application deemed withdrawn after publication