CN102834059A - 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 - Google Patents
超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 Download PDFInfo
- Publication number
- CN102834059A CN102834059A CN2011800175532A CN201180017553A CN102834059A CN 102834059 A CN102834059 A CN 102834059A CN 2011800175532 A CN2011800175532 A CN 2011800175532A CN 201180017553 A CN201180017553 A CN 201180017553A CN 102834059 A CN102834059 A CN 102834059A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- characteristic
- characteristic quantity
- point
- characteristic point
- 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/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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- 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/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- 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/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- 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
-
- 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
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Radar, Positioning & Navigation (AREA)
- Medical Informatics (AREA)
- Remote Sensing (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Vascular Medicine (AREA)
- Data Mining & Analysis (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
超声波观测装置在能够设定以通过对接收到的超声波的频谱进行近似而提取出的特征量的至少一部分为坐标成分的坐标系的特征量空间内,在规定的数据位置处的频谱的特征点在显示部中在第i帧(i是自然数)中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于第一区域的第二区域接近的情况下,在与最新特征点的位置相比远离第二区域且处于第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与规定的数据位置的在第(i+1)帧的最新特征点对应的视觉信息,由此生成特征量图像数据,进行显示。
Description
技术领域
本发明涉及一种使用超声波来观测检体的组织的超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序。
背景技术
以往,作为使用超声波进行的乳腺癌等的检查技术,已知一种超声波弹性成像技术(例如参照专利文献1)。超声波弹性成像是利用生物体内的癌、肿瘤组织的硬度由于疾病的发展状况、生物体不同而不同的技术。在该技术中,在从外部挤压检查部位的状态下,使用超声波对该检查部位的生物体组织的应变量、弹性率进行测量,将该测量结果转换为断层图像来进行图像显示。
专利文献1:国际公开第2005/122906号
发明内容
发明要解决的问题
通常,使用超声波观测检体的组织的超声波观测装置以在检体内具有一定厚度的切片平面为观测对象。另外,在进行观测时存在以下情况:由于超声波观测装置所具备的超声波探头(内窥镜)和检体中的至少一个移动而两者的相对位置关系发生变化。在该情况下,超声波探头观测的切片平面发生变化。当超声波观测装置观测的切片平面发生变化时,有时在最初的帧中观测到的检体内的组织在其后的帧中观测不到。认为产生这种情形是由于在最初的帧中观测到的组织微小或者位于远离超声波探头的位置。
如上所述,在以往的超声波观测装置中存在以下问题:与观测对象之间的相对位置关系随着时间经过而发生变化,由此无法长时间观测期望的组织。
本发明是鉴于上述情况而完成的,目的在于提供一种即使与观测对象之间的相对位置关系随着时间经过发生变化也能够长时间观测期望的组织的超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序。
用于解决问题的方案
为了解决上述问题而达到目的,本发明所涉及的超声波观测装置对检体发送超声波,并且接收由上述检体反射的超声波,该超声波观测装置的特征在于,具备:频率分析部,其对接收到的超声波算出决定的多个数据位置处的频谱;特征量提取部,其通过对上述频率分析部算出的频谱进行近似来提取针对上述频谱的一个或者多个特征量;图像处理部,其依次生成包含与上述特征量提取部提取出的特征量有关的信息的特征量图像数据;以及显示部,其依次显示与上述图像处理部依次生成的特征量图像数据对应的图像,其中,在能够设定以上述一个或者多个特征量的至少一部分为坐标成分的坐标系的特征量空间内,在规定的数据位置处的频谱的特征点在上述显示部上的第i帧中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于上述第一区域的第二区域接近的情况下,其中,i是自然数,上述图像处理部在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与上述规定的数据位置的在第(i+1)帧中的上述最新特征点对应的视觉信息,由此生成上述特征量图像数据。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述图像处理部将同一数据位置处的频谱的特征点,即生成上述第(i+1)帧时算出的第一特征点以及生成上述第i帧时算出的第二特征点投影到分别通过上述第一区域的代表点和上述第二区域的代表点的基准轴,根据投影点与上述第二区域的代表点之间的距离来判断上述第一特征点和上述第二特征点与上述第二区域的位置关系。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述第一特征点是根据上述特征量提取部提取出的特征量而决定的特征点,上述第二特征点是上述虚拟的特征点。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述图像处理部以与分别通过上述第一区域的代表点和上述第二区域的代表点的基准轴正交且通过生成上述第i帧时作为同一数据位置处的特征点而分别算出的第三特征点和第四特征点中的某一个的直线或者平面为基准边界,来判断生成上述第(i+1)帧时算出的特征点与上述第二区域的位置关系,其中,上述第三特征点是根据上述特征量提取部提取出的特征量而决定的特征点,上述第四特征点是上述虚拟的特征点。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述图像处理部使用上述第i帧中的上述第三特征点和上述第四特征点中的包含于上述基准边界的特征点的位置以及上述基准轴,来决定上述第(i+1)帧中的上述虚拟的特征点的位置。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述图像处理部使用上述第i帧中的上述第三特征点和上述第四特征点中的包含上述基准边界的特征点的位置以及上述第二区域的代表点的位置,来决定上述第(i+1)帧中的上述虚拟的特征点的位置。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述特征量提取部对上述频率分析部算出的频谱进行衰减校正处理以及近似处理来提取上述特征量,该衰减校正处理用于削减超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述特征量提取部具有:近似部,其对由上述频率分析部算出的频谱进行上述近似处理,由此提取进行上述衰减校正处理前的校正前特征量;以及衰减校正部,其对上述近似部提取出的校正前特征量进行上述衰减校正处理,由此提取上述频谱的特征量。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述特征量提取部具有:衰减校正部,其对上述频谱进行上述衰减校正处理;以及近似部,其对由上述衰减校正部校正后的频谱进行上述近似处理,由此提取上述频谱的特征量。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,超声波的接收深度越深,上述衰减校正部进行的校正越大。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述近似部通过回归分析以多项式来对上述频谱进行近似。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述视觉信息是构成颜色空间的变量。
另外,本发明涉及一种超声波观测装置的动作方法,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作方法的特征在于,具有以下步骤:频率分析步骤,由频率分析部通过对接收到的超声波的频率进行分析来算出频谱;特征量提取步骤,对通过上述频率分析步骤算出的频谱进行近似来提取针对上述频谱的一个或者多个特征量;图像处理步骤,通过图像处理部生成包含与通过上述特征量提取步骤提取出的特征量有关的信息的特征量图像数据;以及显示步骤,通过显示部依次显示与通过上述图像处理步骤依次生成的特征量图像数据对应的图像,反复进行上述频率分析步骤至上述显示步骤,在能够设定以上述一个或者多个特征量的至少一部分为坐标成分的坐标系的特征量空间上,在规定的数据位置处的频谱的特征点在上述显示部上的第i帧中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于上述第一区域的第二区域接近的情况下,其中,i是自然数,在上述图像处理步骤中在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与上述规定的数据位置的在第(i+1)帧中的上述最新特征点对应的视觉信息,由此生成上述特征量图像数据。
另外,本发明所涉及一种超声波观测装置的动作程序,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作程序的特征在于,使该超声波观测装置反复执行以下步骤:频率分析步骤,由频率分析部通过对接收到的超声波的频率进行分析来算出频谱;特征量提取步骤,对通过上述频率分析步骤算出的频谱进行近似来提取针对上述频谱的一个或者多个特征量;图像处理步骤,通过图像处理部生成包含与通过上述特征量提取步骤提取出的特征量有关的信息的特征量图像数据;以及显示步骤,通过显示部依次显示与通过上述图像处理步骤依次生成的特征量图像数据对应的图像,其中,在能够设定以上述一个或者多个特征量的至少一部分为坐标成分的坐标系的特征量空间上,在规定的数据位置处的频谱的特征点在上述显示部上的第i帧中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于上述第一区域的第二区域接近的情况下,其中,i是自然数,在上述图像处理步骤中在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与上述规定的数据位置的在第(i+1)帧中的上述最新特征点对应的视觉信息,由此生成上述特征量图像数据。
发明的效果
根据本发明,在特征量空间上,在规定的数据位置处的频谱的特征点在第i帧(i是自然数)中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于第一区域的第二区域接近的情况下,在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与规定的数据位置的在第(i+1)帧中的最新特征点对应的视觉信息,由此生成包含与特征量有关的信息的特征量图像数据,依次显示与所生成该特征量图像数据对应的图像,因此能够尽可能地延长显示具有与图像显示优先级高的区域附近的特征点对应的像素值的图像。因而,即使与观测对象之间的相对位置关系随着时间经过发生变化,也能够长时间观测期望的组织。
附图说明
图1是表示本发明的实施方式1所涉及的超声波观测装置的结构的框图。
图2是示意性地表示本发明的实施方式1所涉及的超声波观测装置所存储的特征量空间的结构的图。
图3是表示本发明的实施方式1所涉及的超声波观测装置的处理的概要的流程图。
图4是表示本发明的实施方式1所涉及的超声波观测装置的显示部中的B模式图像的显示例的图。
图5是表示本发明的实施方式1所涉及的超声波观测装置的频率分析部所进行的处理的概要的流程图。
图6是示意性地表示一个声线的数据排列的图。
图7是表示本发明的实施方式1所涉及的超声波观测装置的频率分析部算出的频谱的例子(第一例)的图。
图8是表示本发明的实施方式1所涉及的超声波观测装置的频率分析部算出的频谱的例子(第二例)的图。
图9是表示根据对与图7示出的直线有关的特征量进行衰减校正后的特征量而确定的新的直线的图。
图10是表示本发明的实施方式1所涉及的超声波观测装置的特征量图像数据生成部所进行的特征量图像数据生成处理的概要的图。
图11是表示本发明的实施方式1所涉及的超声波观测装置的特征量图像数据生成部所进行的余像条件判定处理的概要的流程图。
图12是示意性地表示在本发明的实施方式1所涉及的超声波观测装置所进行的余像条件判定处理中判断为满足余像条件的情况的图。
图13是表示在本发明的实施方式1所涉及的超声波观测装置所进行的余像条件判定处理中判断为没有满足余像条件的情况的例子(第一例)的图。
图14是表示在本发明的实施方式1所涉及的超声波观测装置所进行的余像条件判定处理中判断为没有满足余像条件的情况的例子(第二例)的图。
图15是表示本实施方式1所涉及的超声波观测装置所进行的余像用特征点的计算方法的概要的图。
图16是表示本发明的实施方式1所涉及的超声波观测装置的显示部所显示的特征量图像的显示例的图。
图17是以黑白方式示意性地表示图16示出的图像的图。
图18是示意性地表示本发明的实施方式1所涉及的超声波观测装置进行观测的状况的图。
图19是示意性地表示在以往的超声波观测装置中在超声波探头与检体的相对位置关系随着时间经过而发生变化的状况下显示部所显示的特征量图像的显示例(第一例)的图。
图20是示意性地表示在以往的超声波观测装置中在超声波探头与检体的相对位置关系随着时间经过而发生变化的状况下显示部所显示的特征量图像的显示例(第二例)的图。
图21是示意性地表示在本发明的实施方式1所涉及的超声波观测装置中在与图19相同的状况下显示部所显示的特征量图像的显示例的图。
图22是示意性地表示在本发明的实施方式1所涉及的超声波观测装置中在与图20相同的状况下显示部所显示的特征量图像的显示例的图。
图23是说明本发明的实施方式1所涉及的超声波观测装置所进行的衰减校正处置的效果的图。
图24是表示本发明的实施方式2所涉及的超声波观测装置的特征量图像数据生成部所进行的特征量图像数据生成处理的概要的流程图。
图25是表示本发明的实施方式2所涉及的超声波观测装置所进行的余像条件判定处理的概要的流程图。
图26是表示在本发明的实施方式2所涉及的超声波观测装置所进行的余像条件判定处理中判断为满足余像条件的情况下的特征量空间的结构的图。
图27是表示在本发明的实施方式2所涉及的超声波观测装置所进行的余像条件判定处理中判断为没有满足余像条件的情况下的特征量空间的结构的图。
图28是表示本发明的实施方式2所涉及的超声波观测装置所进行的余像用特征点算出处理(第一例)的概要的图。
图29是表示本发明的实施方式2所涉及的超声波观测装置所进行的余像用特征点算出处理(第二例)的概要的图。
图30是表示本发明的实施方式3所涉及的超声波观测装置的处理的概要的流程图。
图31是表示本发明的实施方式3所涉及的超声波观测装置所进行的衰减校正处理的概要的图。
具体实施方式
下面,参照附图说明用于实施本发明的方式(以下称为“实施方式”)。
(实施方式1)
图1是表示本发明的实施方式1所涉及的超声波观测装置的结构的框图。该图示出的超声波观测装置1是使用超声波来判断作为观测对象的检体的组织性状的装置。
超声波观测装置1具备:超声波探头2,其向外部输出超声波脉冲,并且接收在外部反射的超声波回波;发送和接收部3,其与超声波探头2之间进行电信号的发送和接收;运算部4,其对变换超声波回波而得到的电回波信号实施规定的运算;图像处理部5,其使用变换超声波回波而得到的电回波信号来生成各种图像数据;输入部6,其使用键盘、鼠标、触摸面板等接口来实现,接受各种信息的输入;显示部7,其使用由液晶或者有机EL等构成的显示面板来实现,显示包括由图像处理部5生成的图像的各种信息;存储部8,其存储包括与已知检体的组织性状有关的信息的各种信息;以及控制部9,其对超声波观测装置1进行动作控制。
超声波探头2具有信号变换部21,该信号变换部21将从发送和接收部3接收到的电脉冲信号变换为超声波脉冲(声脉冲信号),并且将由外部的检体反射而得到的超声波回波变换为电回波信号。超声波探头2可以是以机械方式使超声波振子进行扫描的探头,也可以是以电子方式使多个超声波振子进行扫描的探头。
发送和接收部3与超声波探头2电连接,将脉冲信号发送给超声波探头2,并且从超声波探头2接收作为接收信号的回波信号。具体地说,发送和接收部3根据预先设定的波形和发送定时来生成脉冲信号,将所生成的该脉冲信号发送给超声波探头2。
发送和接收部3与超声波探头2电连接,将脉冲信号发送给超声波探头2,并且从超声波探头2接收回波信号。具体地说,发送和接收部3根据预先设定的波形和发送定时来生成脉冲信号,将所生成的该脉冲信号发送给超声波探头2。另外,发送和接收部3在对接收到的回波信号实施放大、滤波等处理之后,通过进行A/D变换来生成数字RF信号而输出。此外,在超声波探头2是以电子方式使多个超声波振子进行扫描的探头的情况下,发送和接收部3具有与多个超声波振子对应的束合成用的多通道电路。
运算部4具有:频率分析部41,其通过对由发送和接收部3输出的数字RF信号实施高速傅里叶变换(FFT)来对回波信号进行频率分析;以及特征量提取部42,其对由频率分析部41算出的频谱进行衰减校正处理以及近似处理,由此提取检体的特征量,其中,在该衰减校正处理中削减超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献。
频率分析部41针对各声线(行数据),对由规定的数据量构成的FFT数据群进行高速傅里叶变换,由此算出频谱。频谱示出根据检体的组织性状不同而不同的趋势。这是由于,频谱与作为使超声波散射的散射体的检体的大小、密度、声阻抗等具有相关性。在此所指的“组织性状”例如是指癌症、内分泌肿瘤、粘液性肿瘤、正常组织、脉管等中的任一个,在检体为胰腺的情况下,作为组织性状还包含慢性胰腺炎、自身免疫性胰腺炎等。
特征量提取部42具有:近似部421,其对由频率分析部41算出的频谱进行近似处理,由此算出进行衰减校正处理前的校正前特征量;以及衰减校正部422,其对由近似部421进行近似而得到的校正前特征量进行衰减校正处理,由此提取特征量。
近似部421通过回归分析以一次式对频谱进行近似,由此提取使该近似的一次式具有特征的校正前特征量。具体地说,近似部421通过回归分析来算出一次式的斜率a0和截距b0,并且算出频谱中的频带内的特定频率处的强度作为校正前特征量。在本实施方式1中,近似部421算出中心频率fMID=(fLOW+fHIGH)/2处的强度(Mid-band fit)c0=a0fMID+b0,但是这严格说来仅是一例。在此所指的“强度”是指电压、电力、声压、声能等参数中的任一个。
三个特征量中的斜率a0与超声波散射体的大小具有相关性,通常认为散射体越大则斜率具有越小的值。另外,截距b0与散射体的大小、声阻抗的差、散射体的密度(浓度)等具有相关性。具体地说,认为对于截距b0,散射体越大则具有越大的值,声阻抗越大则具有越大的值,散射体的密度(浓度)越大则具有越大的值。中心频率fMID处的强度(以下简单称为“强度”)c0是根据斜率a0和截距b0导出的间接的参数,提供有效频带内的中心处的频谱强度。因此,认为强度c0除了与散射体的大小、声阻抗的差、散射体的密度相关以外,还与B模式图像的亮度具有某种程度的相关性。此外,由特征量提取部42算出的近似多项式并不限定于一次式,还能够使用二次以上的近似多项式。
说明由衰减校正部422进行的校正。超声波的衰减量A能够表示为如下。
A=2αzf …(1)
在此,α是衰减率,z是超声波的接收深度,f是频率。根据式(1)可知,衰减量A与频率f成比例。对于衰减率α的具体的值,在生物体的情况下为0~1.0(dB/cm/MHz),更优选为0.3~0.7(dB/cm/MHz),根据作为观察对象的脏器的种类来决定。例如,在作为观察对象的脏器为胰腺的情况下,决定为α=0.6(dB/cm/MHz)。此外,在本实施方式1中,还能够设为以下结构:能够根据来自输入部6的输入来变更衰减率α的值。
衰减校正部422如下那样对由近似部421提取出的校正前特征量(斜率a0、截距b0、强度c0)进行校正。
a=a0+2αz …(2)
b=b0 …(3)
c=c0+2αzfMID(=afMID+b) …(4)
根据式(2)、(4)可知,衰减校正部422进行如下校正:超声波的接收深度z越大则校正量越大。另外,根据式(3),与截距相关的校正是恒等变换。这是由于,截距是与频率0(Hz)对应的频率成分而没有被衰减。
图像处理部5具有:B模式图像数据生成部51,其根据回波信号来生成B模式图像数据;以及特征量图像数据生成部52,其生成包含与特征量有关的信息的特征量图像数据。
B模式图像数据生成部51对数字信号进行带通滤波、对数变换、增益处理、对比度处理等使用公知技术的信号处理,并且进行数据的间隔剔除等,由此生成B模式图像数据,其中,与根据显示部7中的图像的显示范围而决定的数据步长相应地进行数据的间隔剔除。
特征量图像数据生成部52使用由B模式图像数据生成部51生成的B模式图像数据以及由特征量提取部42提取出的检体的特征量,在时间上连续地生成特征量图像数据。具体地说,特征量图像数据生成部52在能够设定以特征量提取部42提取出的特征量的至少一部分为坐标成分的坐标系的特征量空间上,在同一数据位置处的频谱的特征点随着时间经过而从规定的第一区域内向图像显示的优先级低于第一区域的第二区域接近的情况下,在与特征点的最新位置相比远离第二区域且处于第一区域内部或者周边的位置处设置虚拟的特征点,对上述同一数据位置附加与该虚拟的特征点对应的视觉信息,由此连续地生成包含与特征量有关的信息的特征量图像数据。在此,优选将第一区域设定为与例如癌症那样要通过观测发现的重要度高的组织性状对应的区域。
在本实施方式1中,对第一区域设置如上所述那样用于显示余像的虚拟的特征点,但是对于第二区域,即使同一数据位置处的频谱的特征点随着时间经过偏离第二区域,也不设置虚拟的特征点。换言之,特征量图像数据生成部52仅对存在于特征量空间上的第一区域内的特征点进行余像处理。
存储部8具有:已知检体信息存储部81,其存储包含已知检体的特征量的已知检体信息;窗函数存储部82,其存储频率分析部41进行频率分析处理时使用的窗函数;校正信息存储部83,其存储衰减校正部422进行衰减校正处理时参照的校正信息;特征量空间信息存储部84,其存储与特征量空间有关信息,该特征量空间是根据已知检体信息存储部81所存储的已知检体的特征量而设定的;以及特征量信息存储部85,其存储与特征量有关的信息,该特征量被算出为特征量空间信息存储部84所存储的特征量空间内的点的坐标值。
已知检体信息存储部81存储对已知检体进行频率分析而提取出的频谱的特征量。另外,已知检体信息存储部81存储对与已知检体相关联的频谱的特征量按根据已知检体的组织性状的种类分类的每个组而算出的平均和标准偏差。在本实施方式1中,超声波接收信号的频谱的特征量的平均和标准偏差反应检体中的核肿大、异形等细胞水平的变化、间质中的纤维增生、实质组织被纤维替换等组织变化,与组织性状相应地呈现特有的值,鉴于此使用已知检体的频谱特征量的平均和标准偏差来进行组织性状的分类。
窗函数存储部82存储Hamming、Hanning、Blackman等窗函数中的至少一个窗函数。校正信息存储部83存储与式(2)~(4)的变换有关的信息。
特征量空间信息存储部84存储与根据由已知检体信息存储部81存储的已知检体信息而设定的特征量空间有关的信息,作为该信息,存储根据特征量对多个已知检体进行分类而得到的多个组和各组的代表点。对于该代表点,例如可以设为组内的特征量的平均,也可以设为组内的特征量的重心。
图2是示意性地表示特征量空间信息存储部84所存储的特征量空间的结构的图。在图2示出的特征量空间内,横轴为截距b、纵轴为强度c(参照式(3)、(4))。另外,图2示出的区域Gμ、Gν表示已知检体信息存储部81所存储的已知检体的组织性状分别为μ、ν的组。区域Gμ内的点μ0和区域Gν内的点ν0是构成各区域的要素的特征点的代表点。在本实施方式1中,设区域Gμ为第一区域(以下称为余像区域),区域Gν为第二区域。
特征量信息存储部85存储算出作为视觉信息的像素值时使用的特征点。作为这种特征点,除了由特征量提取部42提取出的特征点以外,还存在根据特征点的位置而决定的虚拟的特征点即余像用特征点。另外,特征量信息存储部85在规定的存储区域存储余像区域标志,根据算出像素值时使用的特征点是否为余像区域内的点,该余像区域标志取不同的值。例如,特征量信息存储部85将使用余像区域内的特征点或者余像用特征点来算出像素值的情况下的余像区域标志设为1进行存储,另一方面将使用余像区域外的特征点或者余像用特征点算出像素值的情况下的余像区域标志设为0进行存储。
在图2中,还记载了通过代表点μ0和ν0并以μ0→ν0的朝向为正向的h轴(基准轴)。关于该h轴,在特征量图像数据生成部52决定帧中像素的像素值时、判断使用最新特征点的像素值还是使用已经存储的余像用特征点的像素值时使用。
在图2示出的情况下,两个区域Gμ、Gν在特征量空间上存在于相互不交叉的区域内。这样,在本实施方式1中,以通过频率分析得到的频谱的特征量为指标进行分组,因此能够严加区别相互不同的组。特别是,在本实施方式1中,对超声波的回波信号进行衰减校正,因此与不进行衰减校正的情况相比,还能够以更清楚地分离的状态得到特征量空间上的各组的区域。此外,在特征量空间上的b轴成分与c轴成分的刻度大不相同的情况下,期望适当地进行加权以使各距离的贡献大致均等。
另外,特征量空间信息存储部84存储特征量空间上的点与像素值之间的关系作为按每个像素决定的视觉信息即像素值的信息。例如,特征量空间信息存储部84存储对截距b与强度c分配的构成颜色空间的变量的值。在此所指的颜色空间例如是指RGB表色系、补色系的变量、表示光的三个属性(色调、明度、彩度)的表色系等。
通过预先存储了本实施方式1所涉及的超声波观测装置的动作程序、启动规定的OS的程序等的ROM以及存储各处理的运算参数、数据等的RAM等来实现存储部8。
具有上述功能结构的超声波观测装置1的超声波探头2以外的结构要素通过使用具备CPU的计算机来实现,该CPU具有运算和控制功能。超声波观测装置1所具备的CPU从存储部8读出存储部8存储、保存的信息以及包含上述超声波观测装置的动作程序的各种程序,由此执行与本实施方式1所涉及的超声波观测装置的动作方法相关联的运算处理。
此外,也能够将本实施方式1所涉及的超声波观测装置的动作程序记录到硬盘、快闪存储器、CD-ROM、DVD-ROM、软盘等计算机可读取的记录介质而广泛流通。
图3是表示具有上述结构的超声波观测装置1的处理的概要的流程图。图3示出的流程图示出对一帧进行的处理的概要。在以下说明中,将用于识别帧的变量设为i(=1、2、…)。
在图3中,首先,超声波观测装置1将帧识别用的变量i设定为1(步骤S1)。接着,超声波观测装置1首先通过超声波探头2测量新检体(步骤S2)。
之后,B模式图像数据生成部51使用从发送和接收部3输出的B模式图像用回波信号来生成B模式图像数据(步骤S3)。
接着,控制部9进行控制,使显示部7显示与由B模式图像数据生成部51生成的B模式图像数据对应的B模式图像(步骤S4)。图4是表示显示部7中的B模式图像的显示例的图。该图示出的B模式图像100是使采用RGB表色系作为颜色空间的情况下的变量即R(红)、G(绿)、B(蓝)的值一致的灰度等级图像。
之后,频率分析部41通过FFT运算来进行频率分析,由此算出频谱(步骤S5)。在此,参照图5示出的流程图来详细说明由频率分析部41进行的处理(步骤S5)。首先,频率分析部41将最初成为分析对象的声线的声线编号L设为初始值L0(步骤S21)。例如可以对发送和接收部3最初接收到的声线附加初始值L0,也可以对与通过输入部6设定的关心区域的左右一侧的边界位置对应的声线附加初始值L0。
接着,频率分析部41算出设定于一个声线上的多个数据位置的全部位置的频谱。首先,频率分析部41设定代表用于进行FFT运算而获取的一系列数据群(FFT数据群)的数据位置Z(相当于接收深度)的初始值Z0(步骤S22)。图6是示意性地表示一个声线的数据排列的图。在该图示出的声线LD中,白色或者黑色的长方形意味着一个数据。以与发送和接收部3进行的A/D变换中的采样频率(例如50MHz)对应的时间间隔使声线LD离散化。在图6中,示出将声线LD的第一个数据设定为数据位置Z的初始值Z0的情况。此外,图6严格说来仅是一例,能够任意地设定初始值Z0的位置。例如,也可以将与关心区域的上端位置对应的数据位置Z设定为初始值Z0。
之后,频率分析部41获取数据位置Z的FFT数据群(步骤S23),使由窗函数存储部82存储的窗函数作用于获取到的FFT数据群(步骤S24)。这样使窗函数作用于FFT数据群,由此避免FFT数据群在边界处不连续,从而能够防止产生伪像。
接着,频率分析部41判断数据位置Z的FFT数据群是否为正常的数据群(步骤S25)。在此,FFT数据群需要具有2的乘方个数据数。以下,将FFT数据群的数据数设为2n(n为正整数)。FFT数据群为正常是指数据位置Z在FFT数据群中位于从前方起第2n-1的位置。换言之,FFT数据群正常是指在数据位置Z的前方存在2n-1-1(设为=N)个数据而在数据位置Z的后方存在2n-1(设为=M)个数据。在图6示出的情况下,FFT数据群F2、F3、FK-1正常,另一方面FFT数据群F1、FK异常。其中,在图6中,设为n=4(N=7、M=8)。
在步骤S25中的判断的结果是数据位置Z的FFT数据群正常的情况下(步骤S25:“是”),频率分析部41转移到后述的步骤S27。
在步骤S25中的判断的结果是数据位置Z的FFT数据群不正常的情况下(步骤S25:“否”),频率分析部41对不足部分插入零数据,由此来生成正常的FFT数据群(步骤S26)。在追加零数据之前使窗函数作用于在步骤S25中判断为不正常的FFT数据群。因此,即使对FFT数据群插入零数据,也不产生数据的不连续。在步骤S26之后,频率分析部41转移到后述的步骤S27。
在步骤S27中,频率分析部41使用FFT数据群来进行FFT运算,由此得到频谱(步骤S27)。图7和图8是表示由频率分析部41算出的频谱例的图。在图7和图8中,横轴f是频率,纵轴I为强度。在图7和图8分别示出的频谱曲线C1和C2中,频谱的下限频率fLOW和上限频率fHIGH是根据超声波探头2的频带、由发送和接收部3发送的脉冲信号的频带等来决定的参数,例如是fLOW=3MHz、fHIGH=10MHz。此外,通过后述的特征量提取处理来说明图7示出的直线L1和图8示出的直线L2。在本实施方式1中,曲线和直线由离散的点的集合构成。这一点在后述的实施方式中也相同。
接着,频率分析部41对数据位置Z加上规定的数据步长D而算出下一个分析对象的FFT数据群的数据位置Z(步骤S28)。期望此处的数据步长D与由B模式图像数据生成部51生成B模式图像数据时利用的数据步长一致,但是在想要削减频率分析部41中的运算量的情况下,也可以设定比B模式图像数据生成部51所利用的数据步长大的值。在图6中,示出D=15的情况。
之后,频率分析部41判断数据位置Z是否大于最终数据位置Zmax(步骤S29)。在此,最终数据位置Zmax可以是声线LD的数据长度,也可以是与关心区域的下端对应的数据位置。在判断的结果是数据位置Z大于最终数据位置Zmax的情况下(步骤S29:“是”),频率分析部41使声线编号L增加1(步骤S30)。另一方面,在数据位置Z为最终数据位置Zmax以下的情况下(步骤S29:“否”),频率分析部41返回到步骤S23。这样,频率分析部41针对一个声线LD,对[{(Zmax-Z0)/D}+1](=K)个FFT数据群进行FFT运算。在此,[X]表示不超过X的最大整数。
在通过步骤S30增加之后的声线编号L大于最终声线编号Lmax的情况下(步骤S31:“是”),频率分析部41返回到图2示出的主例程。另一方面,在通过步骤S30增加之后的声线编号L为最终声线编号Lmax以下的情况下(步骤S31:“否”),频率分析部41返回到步骤S22。
这样,频率分析部41对(Lmax-L0+1)个声线各自进行K次的FFT运算。此外,例如可以对由发送和接收部3接收到的最终的声线附加最终声线编号Lmax,也可以对与关心区域左右的任一个边界对应的声线附加最终声线编号Lmax。以下,将频率分析部41对全部声线进行的FFT运算的总数(Lmax-L0+1)×K设为P。
接在上述说明的步骤S5的频率分析处理之后,作为近似处理,近似部421对频率分析部41算出的P个频谱进行回归分析,由此提取校正前特征量(步骤S6)。具体地说,近似部421通过回归分析来算出对频带fLOW<f<fHIGH的频谱进行近似的一次式,由此作为校正前特征量而提取出使该一次式具有特征的斜率a0、截距b0、强度c0。图7示出的直线L1和图8示出的直线L2是在该步骤S6中对频谱曲线C1和C2分别进行回归分析而得到的回归直线。
之后,衰减校正部422对由近似部421提取出的校正前特征量进行衰减校正处理(步骤S7)。例如在数据的采样频率为50MHz的情况下,数据采样的时间间隔为20(nsec)。在此,当将声速设为1530(m/sec)时,数据采样距离间隔成为1530(m/sec)×20(nsec)/2=0.0153(mm)。当将从声线LD的第一个数据起到处理对象的FFT数据群的数据位置为止的数据步长数为k时,该数据位置Z成为0.0153k(mm)。衰减校正部422将这样求得的数据位置Z的值代入到上述式(2)~(4)的接收深度z,由此算出作为频谱的特征量的斜率a、截距b、强度c。图9是表示根据对与图7示出的直线L1相关联的特征量进行衰减校正之后的特征量而决定的直线的图。表示图9示出的直线L1’的式如下。
I=af+b=(a0+2αZ)f+b0 …(5)
根据该式(5)也可知,与直线L1相比,直线L1’的斜率大且截距的值相同。
之后,特征量图像数据生成部52生成特征量图像数据(步骤S8)。图10是表示特征量图像数据生成部52所进行的特征量图像数据生成处理的概要的流程图。在以下说明中,将用于识别帧所包含的像素的变量设为j(=1,2,…,jmax)。另外,将与第i个帧(第i帧)的像素j对应的特征量即截距bi(j)和强度ci(j)的组(bi(j),ci(j))称为特征点,记载为Si(j)。
在图10中,首先,特征量图像数据生成部52将像素识别用的变量j设定为1(步骤S41)。
之后,特征量图像数据生成部52获取由特征量提取部42提取出的特征点(第一特征点)Si(j)(步骤S42),从特征量空间信息存储部84读出获取到的特征点Si(j)的余像区域标志,判断余像区域标志的值(步骤S43)。在余像区域标志为0的情况下(步骤S43:“0”),特征量图像数据生成部52使用特征点Si(j)来算出像素j的像素值(步骤S44)。
与此相对,在获取到的特征点Si(j)的余像区域标志为1的情况下(步骤S43:“1”),特征量图像数据生成部52判断是否满足余像条件(步骤S45)。
图11是表示余像条件判定处理的概要的流程图。在图11中,特征量图像数据生成部52分别算出将特征点Si(j)和余像用特征点(第二特征点)Ti-1(j)投影到特征量空间内的h轴的投影点S’i(j)和T’i-1(j)的坐标(步骤S61)。
接着,特征量图像数据生成部52将代表点μ0设为起点,判断分别以投影点S’i(j)和T’i-1(j)为终点的两个矢量μ0S’、μ0T’的朝向是否相同(步骤S62)。在判定的结果是两个矢量μ0S’、μ0T’的朝向相同的情况下(步骤S62:“是”),特征量图像数据生成部52算出投影点S’i(j)与代表点μ0之间的距离dS以及投影点T’i-1(j)与代表点μ0之间的距离dT(步骤S63)。
之后,特征量图像数据生成部52对两个距离dS与dT进行比较(步骤S64)。如果比较的结果是dS≥dT(步骤S64:“是”),则特征量图像数据生成部52判断为满足余像条件(步骤S65),转移到步骤S46。图12是表示两个矢量μ0S’、μ0T’具有相同朝向且dS≥dT的特征量空间的状况的图。在该情况下,投影点S’i(j)比投影点T’i-1(j)更接近区域Gν的代表点ν0。因而,基于余像用特征点Ti-1(j)来决定像素值,能够以与特征量处于区域Gμ内时的显示方式接近的显示方式显示图像。换言之,将与特征点Si(j)对应的像素值替换为与虚拟的特征点即余像用特征点Ti-1(j)对应的像素值来进行显示,能够以与特征量处于区域Gμ内时的显示方式接近的显示方式显示图像。
如果步骤S64中的比较结果是dS<dT(步骤S64:“否”),则特征量图像数据生成部52判断为不满足余像条件(步骤S66),转移到步骤S46。图13是表示两个矢量μ0S’、μ0T’具有相同朝向且dS<dT的情况下的特征量空间的状况的图。在该情况下,与投影点T’i-1(j)相比投影点S’i(j)远离区域Gν的代表点ν0。因而,基于区域Gμ内的特征点Si(j)来决定像素值能够显示具有根据区域Gμ内的更正确的位置而附加的像素值的图像。
在步骤S62中,两个矢量μ0S’、μ0T’的朝向不同的情况下(步骤S62:“否”),特征量图像数据生成部52转移到步骤S66。图14是表示两个矢量μ0S’、μ0T’具有不同朝向的情况下的特征量空间的状况的图。在该情况下,与投影点T’i-1(j)相比投影点S’i(j)远离区域Gν的代表点ν0。因而,基于特征点Si(j)来决定像素值能够显示具有根据区域Gμ内的更正确的位置而附加的像素值的图像。
再次参照图10进行说明。在步骤S45的余像条件判定处理的结果是满足余像条件的情况下(步骤S46:“是”),特征量图像数据生成部52使用余像用特征点Ti-1(j)来算出像素值(步骤S47)。另一方面,在没有满足余像条件的情况下(步骤S46:“否”),特征量图像数据生成部52转移到步骤S44。
在步骤S44或者步骤S47之后,特征量图像数据生成部52判断使用于像素值的计算的特征点Si(j)或者余像用特征点Ti-1(j)是否存在于余像区域内(步骤S48)。在特征点Si(j)或者余像用特征点Ti-1(j)不存在于余像区域内的情况下(步骤S48:“否”),特征量图像数据生成部52进行处理将余像区域标志设为0(步骤S49)。
之后,如果像素值的变量j为最大值jmax(步骤S50:“是”),则特征量图像数据生成部52结束特征量图像数据生成处理。另一方面,如果像素值的变量j小于最大值jmax(步骤S50:“否”),则特征量图像数据生成部52使变量j增加1(步骤S51),返回到步骤S42。此外,在此的变量j的增加只不过意味着处理对象的像素的变更。因而,对各个像素附加的特征点、余像用特征点的变量不发生变化。
说明在步骤S48中特征点Si(j)或者余像用特征点Ti-1(j)存在于余像区域内的情况(步骤S48:“是”)。在该情况下,特征量图像数据生成部52进行处理将余像区域标志设为1(步骤S52),使用特征点Si(j)算出余像用特征点Ti(j),记录到特征量空间信息存储部84(步骤S53),转移到步骤S50。图15是表示余像用特征点Ti(j)的计算方法的概要的图。如图15所示,余像用特征点Ti(j)是使特征点Si(j)的各成分偏移规定量而得到的。即通过(bi(j)+β、ci(j)+γ)来提供余像用特征点Ti(j)的坐标。在此,β、γ均为正的常数,被设定为使余像用特征点Ti(j)比特征点Si(j)更接近区域Gν的代表点ν0的值。此外,在图15中,图示要偏移的方向与h轴平行的情况,但是这严格说来仅是一例。
接着,显示部7显示由特征量图像数据生成部52生成的特征量图像(步骤S9)。图16是表示由显示部7显示的特征量图像的显示例的图。图17是以黑白方式示意性地表示图16示出的图像的图。这些图示的特征量图像200与B模式图像100相比,图像彩色化,与组相应的颜色差异变得清晰。特征量图像200大致由绿色系区域200g和红色系区域200r构成,用黄色系的颜色显示两个区域的边界部(在图17中未图示)。如图16所示,各区域并非由单一颜色构成。例如,绿色系区域200g是由接近绿色的颜色构成的像素聚集的区域。同样地,红色系区域200r是由接近红色的颜色构成的像素聚集的区域。观察这种特征量图像200的观测者能够明确地识别组的差异即组织性状的差异。此外,控制部9也可以使显示部7排列显示特征量图像和B模式图像。由此,能够在一个画面上识别两个图像的差异。
之后,在通过输入部6输入了结束处理的指示信号的情况下(步骤S10:“是”),超声波观测装置1结束一系列处理。与此相对,在没有通过输入部6输入结束处理的指示信号的情况下(步骤S10:”否”),超声波观测装置1使帧识别用的变量i增加1(步骤S11),返回到步骤S2。
在此,说明本实施方式1的效果。通常,超声波观测装置1以在检体内具有一定厚度的切片平面为观测对象。图18是示意性地表示超声波观测装置1进行观测的状况的图。如图18所示,在超声波观测装置1的前端部设置有超声波探头2。超声波探头2具有由多个振子构成的振子群22。构成振子群22的多个振子分别在规定的定时依次发送超声波。所送信的超声波到达在检体内具有一定厚度的切片平面PL。此外,在图18中,将切片平面PL记载为长方体状,但是这是为了方便起见而设置的,切片平面PL的实际形状根据检体的观测部分不同而发生变化。
在图18示出的状况下,在以往的超声波观测装置中,存在超声波探头2和检体中的至少一个移动而两者的相对位置关系发生变化的情况。在该情况下,作为观测对象的切片平面PL发生变化。当切片平面PL发生变化时,有时会产生在最初的帧中观测到的检体内的组织Or在其后的帧中观测不到的情形。认为产生这种情形是由于在最初的帧中观测到的组织Or微小或者位于远离超声波探头的位置。
图19是示意性地表示在以往的超声波观测装置中在超声波探头与检体的相对位置关系随着时间经过而发生变化的状况下显示部所显示的特征量图像的显示例(第一例)的图。在图19示出的情况下,在最初的特征量图像301中显示的组织Or1在其后连续显示的特征量图像302、303、304中完全消失。
图20是示意性地表示在以往的超声波观测装置中在超声波探头与检体的相对位置关系随着时间经过而发生变化的状况下显示部所显示的特征量图像的显示例(第二例)的图。在图20示出的情况下,在最初的特征量图像401中显示的组织Or2在其后显示的特征量图像402中消失。另外,在特征量图像402之后显示的特征量图像403中,显示形状与组织Or2不同的组织Or3。该组织Or3在其后的特征量图像404中消失。
如图19以及图20所示,在以往的超声波观测装置中,存在以下问题:由于超声波探头与检体的相对位置关系随着时间经过发生变化而无法长时间观测期望的组织。
图21以及图22是说明本实施方式1所涉及的超声波观测装置1所进行的处理的效果的图。其中,图21是示意性地表示与图19相同的状况下的显示部7的显示例的图。如图21示出的情况下,与图19同样地,组织Or1仅在拍摄最初显示的特征量图像501的时刻被观测到。然而,通过进行上述处理,在后续的特征量图像502、503、504中,在与组织Or1相同的位置处也分别显示相同形状的余像Or11、Or12、Or13。其中,这些余像Or11、Or12、Or13的显示颜色与组织Or1的显示颜色相似,但是有少许不同。这是由于,在特征量空间内,与构成各余像的区域内的像素对应的特征点和与组织Or1内的像素对应的特征点不同。
图22是示意性地表示与图20相同状况下的显示部7的显示例的图。在图22示出的情况下,组织Or2仅在拍摄最初显示的特征量图像601的时刻被观测到而显示。在图22中,在与图20的特征量图像402对应的特征量图像602中,显示组织Or2的余像Or21。另外,在特征量图像602后续的特征量图像603中,与图20的特征量图像403同样地显示组织Or3。在特征量图像603后续的特征量图像604中,显示组织Or3的余像Or31。此外,在图22示出的情况下,组织Or2的显示颜色与余像Or21的显示颜色不同,组织Or3的显示颜色与余像Or31的显示颜色不同。
这样,在本实施方式1中,在显示特征量图像时,即使在以往没有显示组织的帧中,也基于前一帧显示的组织或者余像来描绘余像,因此即使在某一帧中无法观测显示优先级高的组织,通过使用其前一帧的图像来进行余像显示,也能够长时间显示期望的组织。其结果是用户能够连续地观察想看的图像。
图23是说明本实施方式1中的衰减校正的效果的图。图23示出的图像700是没有进行衰减校正的情况下的特征量图像。该情况下的特征量图像是针对由B模式图像数据生成部51生成的B模式图像将特征量b均匀地分配给R(红)、G(绿)、B(蓝)的灰度图像。特征量图像700在接收深度大的区域(图的下方区域)内由于衰减的影响而信号强度下降,图像变暗。与此相对,可知在使用同一B模式图像来进行衰减校正后的特征量图像800中,得到画面整体亮度均匀的图像。
根据上述说明的本发明的实施方式1,在特征量空间内,在规定的数据位置处的频谱的特征点在第i帧(i是自然数)中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于第一区域的第二区域接近的情况下,在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息,替换与规定的数据位置的在第(i+1)帧中的最新特征点对应的视觉信息,由此生成包含与特征量有关的信息的特征量图像数据,依次显示与所生成的该特征量图像数据对应的图像,因此尽可能地延长显示具有与图像显示优先级高的区域附近的特征点对应的像素值的图像。因而,即使与观测对象之间的相对位置关系随着时间经过发生变化,也能够长时间观测期望的组织。
另外,根据本实施方式1,通过对分析接收到的超声波的频率而得到的频谱的校正前特征量进行衰减校正来提取特征量,使用提取出的该特征量并且使用根据由多个已知检体反射的超声波而提取出的频谱的特征量,因此不使用生物体组织的应变量、弹性率就能够明确地严加区别组织的差异。因而,通过高精度地观测检体,能够提高观测结果的可靠性。
另外,根据本实施方式1,在求出已知检体的特征量时,以对通过频率分析得到的频谱的校正前特征量进行衰减校正而得到的特征量为指标进行组织性状的分类、判断,因此能够严加区别相互不同的组织性状。特别是,在本实施方式1中,使用进行了衰减校正的特征量,因此与使用不进行衰减校正而提取出的特征量的情况相比,能够以更清楚地分离的状态得到特征量空间上的各组的区域。
在超声波弹性成像中,存在按压的压力难以传递到血管、淋巴管等脉管下部的问题。因此,在脉管附近形成肿瘤的情况下,肿瘤的边界不清楚,也难以鉴别肿瘤向脉管内的浸润,有时无法高精度地进行检体的观测。另外,在超声波弹性成像中,在检查者挤压检查部位时的压力、挤压速度上容易产生个体差异,因此还存在测量结果的可靠性低这种问题。根据本实施方式1,能够高精度地观测检体来提高观测结果的可靠性,因此能够提供一种适合于解决超声波弹性成像特有的问题的技术。
(实施方式2)
在本发明的实施方式2中,由特征图像数据生成部进行的特征图像数据生成处理与实施方式1不同。本实施方式2所涉及的超声波观测装置的结构与在实施方式1中说明的超声波观测装置1的结构相同。因此,在以下说明中,对与超声波观测装置1的结构要素对应的结构要素附加相同的附图标记。
在本实施方式2中,特征量信息存储部85还存储了边界用特征点Ui(j)。该边界用特征点Ui(j)是使用于像素值的计算的特征点Si(j)或者余像用特征点Ti(j),用于在下一帧中设定在特征量空间中判断余像条件时所需的边界。
图24是表示本实施方式2所涉及的超声波观测装置1的特征量图像数据生成部52所进行的特征量图像数据生成处理的概要的流程图。在图24中,特征量图像数据生成部52将像素识别用的变量j设定为1(步骤S71)。
之后,特征量图像数据生成部52获取由特征量提取部42提取出的特征点(第三特征点)Si(j)(步骤S72),从特征量空间信息存储部84读出获取到的特征点Si(j)的余像区域标志,判断余像区域标志的值(步骤S73)。在余像区域标志为0的情况下(步骤S73:“0”),特征量图像数据生成部52使用特征点Si(j)来算出像素j的像素值(步骤S74)。
与此相对,在获取到的特征点Si(j)的余像区域标志为1的情况下(步骤S73:“1”),特征量图像数据生成部52判断是否满足余像条件(步骤S75)。
图25是表示余像条件判定处理的概要的流程图。在图25中,首先,特征量图像数据生成部52在特征量空间内设定通过边界用特征点Ui-1(j)并且与基准轴(h轴)正交的基准边界(步骤S91)。
接着,特征量图像数据生成部52判断相对于基准边界特征点Si(j)与第二区域是否处于相同一侧(步骤S92)。在相对于基准边界特征点Si(j)与第二区域处于相同一侧的情况下(步骤S92:“是”),特征量图像数据生成部52判断为满足余像条件(步骤S93),转移到步骤S76。图26是表示在由超声波观测装置1所进行的余像条件判定处理中判断为满足余像条件的情况下的特征量空间的结构的图。在图26中,相对于由通过边界用特征点Ui-1(j)并且与h轴正交的直线构成的基准边界B,特征点Si(j)与第二区域即区域Gν处于相同一侧。此外,在特征量空间为三维空间的情况下,基准边界为平面。
另一方面,在相对于基准边界特征点Si(j)与第二区域不处于相同一侧的情况下(步骤S92:“否”),特征量图像数据生成部52判断为没有满足余像条件(步骤S94),转移到步骤S76。图27是表示在超声波观测装置1所进行的余像条件判定处理中判断为没有满足余像条件的情况下的特征量空间的结构的图。在图27中,相对于基准边界B特征点Si(j)与区域Gν位于不同侧。
再次参照图24进行说明。在步骤S75的余像条件判定处理的结果是满足余像条件的情况下(步骤S76:“是”),特征量图像数据生成部52算出余像用特征点(第四特征点)Ti(j)(步骤S77)。
图28是表示步骤S77的余像用特征点算出处理(第一例)的概要的图。特征量图像数据生成部52首先求出通过特征点Si(j)且与基准边界B正交的直线(与h轴平行的直线),求出该直线与基准边界B的交点U’i-1(j),算出该交点U’i-1(j)与特征点Si(j)之间的距离dS1。接着,特征量图像数据生成部52将距边界用特征点Ui(j)距离为dT1=dS1/n(n>1,n是常数)、比边界用特征点Ui(j)更接近区域Gν且位于通过边界用特征点Ui-1(j)而与基准边界B正交的直线上的点设为余像用特征点Ti(j)。此外,还能够将距离dT1定义为固定值。
之后,特征量图像数据生成部52使用余像用特征点Ti(j)来算出像素j的像素值(步骤S78)。
在步骤S76中没有满足余像条件的情况下(步骤S76:“否”),特征量图像数据生成部52转移到步骤S74。
在步骤S74或者步骤S78之后,特征量图像数据生成部52判断使用于像素值的计算的特征点Si(j)或者余像用特征点Ti(j)是否存在于余像区域内(步骤S79)。在特征点Si(j)或者余像用特征点Ti(j)不存在于余像区域内的情况下(步骤S79:“否”),特征量图像数据生成部52进行处理将余像区域标志设为0(步骤S80)。
之后,如果像素值的变量j为最大值jmax(步骤S81:“是”),则特征量图像数据生成部52结束特征量图像数据生成处理。另一方面,如果像素值的变量j小于最大值jmax(步骤S81:“否”),则特征量图像数据生成部52使变量j增加1(步骤S82),返回到步骤S72。此外,在此,变量j的增加只不过意味着处理对象的像素的变更。因而,对各个像素附加的特征点、余像用特征点的变量不会发生变化。
说明在步骤S79中特征点Si(j)或者余像用特征点Ti(j)存在于余像区域内的情况下(步骤S79:“是”)。在该情况下,特征量图像数据生成部52进行处理将余像区域标志设为1(步骤S83),将使用于像素值的计算的特征点Si(j)或者余像用特征点Ti(j)作为边界用特征点Ui(j)而记录到特征量信息存储部85(步骤S84),转移到步骤S81。
根据上述说明的本发明的实施方式2,在特征量空间中,在规定的数据位置处的频谱的特征点在第i帧(i是自然数)中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于第一区域的第二区域接近的情况下,在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与规定的数据位置的在第(i+1)帧中的最新特征点对应的视觉信息,由此生成包含与特征量有关的信息的特征量图像数据,依次显示与所生成的该特征量图像数据对应的图像,因此能够尽可能地延长显示具有与图像显示优先级高的区域附近的特征点对应的像素值的图像。因而,即使与观测对象之间的相对位置关系随着时间经过发生变化,也能够长时间观测期望的组织。
另外,根据本实施方式2,通过对分析接收到的超声波的频率而得到的频谱的校正前特征量进行衰减校正来提取特征量,使用提取出的该特征量并且使用根据由多个已知检体反射的超声波而提取出的频谱的特征量,因此不使用生物体组织的应变量、弹性率就能够清楚地严加区别组织的差异。因而,能够高精度地观测检体,从而提高观测结果的可靠性。其结果是能够提供适合于超声波弹性成像的技术。
图29是表示图24的步骤S77的余像用特征点算出处理的其它例(第二例)的概要的图。特征量图像数据生成部52首先算出边界用特征点Ui-1(j)与特征点Si(j)之间的距离dS2。接着,特征量图像数据生成部52将距边界用特征点Ui-1(j)的距离为dT2=dS2/n(n>1,n是常数)、比边界用特征点Ui-1(j)更接近区域Gν且位于通过边界用特征点Ui-1(j)和区域Gν的代表点ν0的直线上的点设为余像用特征点Ti(j)。
(实施方式3)
在本发明的实施方式3中,由特征量提取部进行的特征量提取处理与实施方式1不同。本实施方式3所涉及的超声波观测装置的结构与在实施方式1中说明的超声波观测装置1的结构相同。因此,在以下说明中,对与超声波观测装置1的结构要素对应的结构要素附加相同的附图标记。
在本实施方式3的特征量提取处理中,首先衰减校正部422对由频率分析部41算出的频谱进行衰减校正处理。之后,近似部421对由衰减校正部422进行衰减校正后的频谱进行近似处理,由此提取频谱的特征量。
图30是表示本实施方式3所涉及的超声波观测装置的处理的概要的流程图。在图30中,步骤S101~S105的处理与图3的步骤S1~S5的处理依次对应。
在步骤S106中,衰减校正部422对由频率分析部41通过FFT演算而算出的全部频谱进行衰减校正(步骤S106)。图31是示意性地表示该步骤S106的处理的概要的图。如图31所示,衰减校正部422针对频谱曲线C3,对所有频率f进行将上述式(1)的衰减量A加到强度I的校正,由此得到新的频谱曲线C3’。由此,能够得到削减了伴随超声波传播的衰减的贡献的频谱。
之后,近似部421对由衰减校正部422进行衰减校正后的所有频谱进行回归分析,由此提取频谱的特征量(步骤S107)。具体地说,近似部421通过回归分析来算出一次式的斜率a、截距b以及中心频率fMID处的强度c。图31示出的直线L3是在该步骤S107中对频谱曲线C3进行特征量提取处理而得到的回归直线(截距b3)。
步骤S108~S111的处理与图3的步骤S8~S11的处理依次对应。
根据上述说明的本发明的实施方式3,在特征量空间内,在规定的数据位置处的频谱的特征点在第i帧(i是自然数)中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于第一区域的第二区域接近的情况下,在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与规定的数据位置的在第(i+1)帧中的最新特征点对应的视觉信息,由此生成包含与特征量有关的信息的特征量图像数据,依次显示与所生成的该特征量图像数据对应的图像,因此能够尽可能地延长显示具有与图像显示优先级高的区域附近的特征点对应的像素值的图像。因而,即使与观测对象之间的相对位置关系随着时间经过发生变化,也能够长时间观测期望的组织。
另外,根据本实施方式3,对分析接收到的超声波的频率而得到的频谱进行衰减校正,根据衰减校正后的频谱提取特征量,使用提取出的该特征量并且使用根据由多个已知检体反射的超声波而提取出的频谱的特征量,因此不使用生物体组织的应变量、弹性率就能够清楚地严加区别组织的差异。因而,能够高精度地观测检体,从而提高观测结果的可靠性。其结果是能够提供也适合于超声波弹性成像的技术。
附图标记说明
1:超声波观测装置;2:超声波探头;3:发送和接收部;4:运算部;5:图像处理部;6:输入部;7:显示部;8:存储部;9:控制部;21:信号变换部;41:频率分析部;42:特征量提取部;51:B模式图像数据生成部;52:特征量图像数据生成部;81:已知检体信息存储部;82:窗函数存储部;83:校正信息存储部;84:特征量空间信息存储部;85:特征量信息存储部;100:B模式图像;200、301~304、401~404、501~504、601~604、700、800:特征量图像;200g:绿色系统区域;200r:红色系统区域;421:近似部;422:衰减校正部。
Claims (14)
1.一种超声波观测装置,对检体发送超声波,并且接收由上述检体反射的超声波,该超声波观测装置的特征在于,具备:
频率分析部,其对接收到的超声波算出决定的多个数据位置处的频谱;
特征量提取部,其通过对上述频率分析部算出的频谱进行近似来提取针对上述频谱的一个或者多个特征量;
图像处理部,其依次生成包含与上述特征量提取部提取出的特征量有关的信息的特征量图像数据;以及
显示部,其依次显示与上述图像处理部依次生成的特征量图像数据对应的图像,
其中,在能够设定以上述一个或者多个特征量的至少一部分为坐标成分的坐标系的特征量空间内,在规定的数据位置处的频谱的特征点在上述显示部上的第i帧中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于上述第一区域的第二区域接近的情况下,
上述图像处理部在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与上述规定的数据位置的在第(i+1)帧中的上述最新特征点对应的视觉信息,由此生成上述特征量图像数据,其中,i是自然数。
2.根据权利要求1所述的超声波观测装置,其特征在于,
上述图像处理部将同一数据位置处的频谱的特征点,即生成上述第(i+1)帧时算出的第一特征点以及生成上述第i帧时算出的第二特征点投影到分别通过上述第一区域的代表点和上述第二区域的代表点的基准轴,根据投影点与上述第二区域的代表点之间的距离来判断上述第一特征点和上述第二特征点与上述第二区域的位置关系。
3.根据权利要求2所述的超声波观测装置,其特征在于,
上述第一特征点是根据上述特征量提取部提取出的特征量而决定的特征点,
上述第二特征点是上述虚拟的特征点。
4.根据权利要求1所述的超声波观测装置,其特征在于,
上述图像处理部以与分别通过上述第一区域的代表点和上述第二区域的代表点的基准轴正交且通过生成上述第i帧时作为同一数据位置处的特征点而分别算出的第三特征点和第四特征点中的某一个的直线或者平面为基准边界,来判断生成上述第(i+1)帧时算出的特征点与上述第二区域的位置关系,
其中,上述第三特征点是根据上述特征量提取部提取出的特征量而决定的特征点,
上述第四特征点是上述虚拟的特征点。
5.根据权利要求4所述的超声波观测装置,其特征在于,
上述图像处理部使用上述第i帧中的上述第三特征点和上述第四特征点中的包含于上述基准边界的特征点的位置以及上述基准轴,来决定上述第(i+1)帧中的上述虚拟的特征点的位置。
6.根据权利要求4所述的超声波观测装置,其特征在于,
上述图像处理部使用上述第i帧中的上述第三特征点和上述第四特征点中的包含上述基准边界的特征点的位置以及上述第二区域的代表点的位置,来决定上述第(i+1)帧中的上述虚拟的特征点的位置。
7.根据权利要求1~6中的任一项所述的超声波观测装置,其特征在于,
上述特征量提取部对上述频率分析部算出的频谱进行衰减校正处理以及近似处理来提取上述特征量,该衰减校正处理用于削减超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献。
8.根据权利要求7所述的超声波观测装置,其特征在于,
上述特征量提取部具有:
近似部,其对由上述频率分析部算出的频谱进行上述近似处理,由此提取进行上述衰减校正处理前的校正前特征量;以及
衰减校正部,其对上述近似部提取出的校正前特征量进行上述衰减校正处理,由此提取上述频谱的特征量。
9.根据权利要求7所述的超声波观测装置,其特征在于,
上述特征量提取部具有:
衰减校正部,其对上述频谱进行上述衰减校正处理;以及
近似部,其对由上述衰减校正部校正后的频谱进行上述近似处理,由此提取上述频谱的特征量。
10.根据权利要求8或者9所述的超声波观测装置,其特征在于,
超声波的接收深度越深,上述衰减校正部的校正量越大。
11.根据权利要求7~10中的任一项所述的超声波观测装置,其特征在于,
上述近似部通过回归分析以多项式来对上述频谱进行近似。
12.根据权利要求1~11中的任一项所述的超声波观测装置,其特征在于,
上述视觉信息是构成颜色空间的变量。
13.一种超声波观测装置的动作方法,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作方法的特征在于,具有以下步骤:
频率分析步骤,由频率分析部通过对接收到的超声波的频率进行分析来算出频谱;
特征量提取步骤,对通过上述频率分析步骤算出的频谱进行近似来提取针对上述频谱的一个或者多个特征量;
图像处理步骤,通过图像处理部生成包含与通过上述特征量提取步骤提取出的特征量有关的信息的特征量图像数据;以及
显示步骤,通过显示部依次显示与通过上述图像处理步骤依次生成的特征量图像数据对应的图像,
反复进行上述频率分析步骤至上述显示步骤,
在能够设定以上述一个或者多个特征量的至少一部分为坐标成分的坐标系的特征量空间上,在规定的数据位置处的频谱的特征点在上述显示部上的第i帧中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于上述第一区域的第二区域接近的情况下,
在上述图像处理步骤中在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与上述规定的数据位置的在第(i+1)帧中的上述最新特征点对应的视觉信息,由此生成上述特征量图像数据,其中,i是自然数。
14.一种超声波观测装置的动作程序,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作程序的特征在于,使该超声波观测装置反复执行以下步骤:
频率分析步骤,由频率分析部通过对接收到的超声波的频率进行分析来算出频谱;
特征量提取步骤,对通过上述频率分析步骤算出的频谱进行近似来提取针对上述频谱的一个或者多个特征量;
图像处理步骤,通过图像处理部生成包含与通过上述特征量提取步骤提取出的特征量有关的信息的特征量图像数据;以及
显示步骤,通过显示部依次显示与通过上述图像处理步骤依次生成的特征量图像数据对应的图像,
其中,在能够设定以上述一个或者多个特征量的至少一部分为坐标成分的坐标系的特征量空间上,在规定的数据位置处的频谱的特征点在上述显示部上的第i帧中存在于规定的第一区域内而在接下来的第(i+1)帧中向图像显示的优先级低于上述第一区域的第二区域接近的情况下,
在上述图像处理步骤中在与最新特征点的位置相比远离上述第二区域且处于上述第一区域的内部或者周边的位置处设置虚拟的特征点,利用与该虚拟的特征点对应的视觉信息来替换与上述规定的数据位置的在第(i+1)帧中的上述最新特征点对应的视觉信息,由此生成上述特征量图像数据,其中,i是自然数。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010253291 | 2010-11-11 | ||
JP2010-253291 | 2010-11-11 | ||
PCT/JP2011/076605 WO2012063977A1 (ja) | 2010-11-11 | 2011-11-11 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102834059A true CN102834059A (zh) | 2012-12-19 |
CN102834059B CN102834059B (zh) | 2013-12-04 |
Family
ID=46051108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011800175532A Active CN102834059B (zh) | 2010-11-11 | 2011-11-11 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8447091B2 (zh) |
EP (1) | EP2526870B1 (zh) |
JP (1) | JP5066306B2 (zh) |
CN (1) | CN102834059B (zh) |
WO (1) | WO2012063977A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104582584A (zh) * | 2013-05-29 | 2015-04-29 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN105828726A (zh) * | 2014-06-11 | 2016-08-03 | 奥林巴斯株式会社 | 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103717146B (zh) * | 2012-05-30 | 2016-01-13 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的动作方法 |
WO2014017255A1 (ja) * | 2012-07-27 | 2014-01-30 | オリンパスメディカルシステムズ株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP5568199B1 (ja) * | 2012-10-01 | 2014-08-06 | オリンパスメディカルシステムズ株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
KR101501479B1 (ko) * | 2013-05-09 | 2015-03-11 | 알피니언메디칼시스템 주식회사 | 초음파 최적화 방법과 그를 위한 초음파 의료 장치 |
EP3395257A4 (en) * | 2015-12-24 | 2019-08-07 | Olympus Corporation | ULTRASONIC OBSERVATION DEVICE, METHOD OF OPERATING THE ULTRASONIC OBSERVATION DEVICE AND A PROGRAM FOR OPERATING THE ULTRASONIC OBSERVATION DEVICE |
US10776989B1 (en) * | 2019-05-13 | 2020-09-15 | Robert Edwin Douglas | Method and apparatus for prioritized volume rendering |
US11275159B2 (en) | 2019-01-08 | 2022-03-15 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Attenuation estimation using ultrasound |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040019276A1 (en) * | 2002-07-23 | 2004-01-29 | Medison Co., Ltd., | Apparatus and method for identifying an organ from an input ultrasound image signal |
US20050203405A1 (en) * | 2004-03-15 | 2005-09-15 | Fuji Photo Film Co., Ltd. | Ultrasonic imaging method and apparatus |
JP2007097671A (ja) * | 2005-09-30 | 2007-04-19 | Toshiba Corp | 画像診断支援装置、及び画像診断支援プログラム |
US20070160275A1 (en) * | 2006-01-11 | 2007-07-12 | Shashidhar Sathyanarayana | Medical image retrieval |
JP2007524431A (ja) * | 2003-02-03 | 2007-08-30 | ザ クリーヴランド クリニック ファウンデーション | 非侵襲性組織特性顕示システムと方法 |
CN101209210A (zh) * | 2006-12-27 | 2008-07-02 | 株式会社东芝 | 超声波图像取得装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0734742B1 (en) * | 1995-03-31 | 2005-05-11 | Kabushiki Kaisha Toshiba | Ultrasound therapeutic apparatus |
WO2002043564A2 (en) * | 2000-11-28 | 2002-06-06 | Allez Physionix Limited | Systems and methods for making non-invasive physiological assessments |
WO2005122906A1 (ja) | 2004-06-18 | 2005-12-29 | Hitachi Medical Corporation | 超音波診断装置 |
-
2011
- 2011-11-11 WO PCT/JP2011/076605 patent/WO2012063977A1/ja active Application Filing
- 2011-11-11 JP JP2012526216A patent/JP5066306B2/ja active Active
- 2011-11-11 CN CN2011800175532A patent/CN102834059B/zh active Active
- 2011-11-11 EP EP11840634.7A patent/EP2526870B1/en not_active Not-in-force
-
2012
- 2012-06-08 US US13/491,782 patent/US8447091B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040019276A1 (en) * | 2002-07-23 | 2004-01-29 | Medison Co., Ltd., | Apparatus and method for identifying an organ from an input ultrasound image signal |
JP2007524431A (ja) * | 2003-02-03 | 2007-08-30 | ザ クリーヴランド クリニック ファウンデーション | 非侵襲性組織特性顕示システムと方法 |
US20050203405A1 (en) * | 2004-03-15 | 2005-09-15 | Fuji Photo Film Co., Ltd. | Ultrasonic imaging method and apparatus |
JP2007097671A (ja) * | 2005-09-30 | 2007-04-19 | Toshiba Corp | 画像診断支援装置、及び画像診断支援プログラム |
US20070160275A1 (en) * | 2006-01-11 | 2007-07-12 | Shashidhar Sathyanarayana | Medical image retrieval |
CN101209210A (zh) * | 2006-12-27 | 2008-07-02 | 株式会社东芝 | 超声波图像取得装置 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104582584A (zh) * | 2013-05-29 | 2015-04-29 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN104582584B (zh) * | 2013-05-29 | 2016-09-14 | 奥林巴斯株式会社 | 超声波观测装置以及超声波观测装置的动作方法 |
CN105828726A (zh) * | 2014-06-11 | 2016-08-03 | 奥林巴斯株式会社 | 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序 |
CN105828726B (zh) * | 2014-06-11 | 2019-06-18 | 奥林巴斯株式会社 | 医用诊断装置以及医用诊断装置的工作方法 |
Also Published As
Publication number | Publication date |
---|---|
US8447091B2 (en) | 2013-05-21 |
EP2526870A4 (en) | 2013-03-06 |
JP5066306B2 (ja) | 2012-11-07 |
EP2526870B1 (en) | 2014-01-01 |
JPWO2012063977A1 (ja) | 2014-05-12 |
CN102834059B (zh) | 2013-12-04 |
WO2012063977A1 (ja) | 2012-05-18 |
US20130011038A1 (en) | 2013-01-10 |
EP2526870A1 (en) | 2012-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102834059B (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN103200876B (zh) | 超声波观测装置、超声波观测装置的动作方法 | |
CN102802536A (zh) | 超声波诊断装置、超声波诊断装置的动作方法以及超声波诊断装置的动作程序 | |
CN103153195B (zh) | 超声波观测装置、超声波观测装置的动作方法 | |
CN102905624A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN102858251B (zh) | 超声波观测装置及超声波观测方法 | |
CN102958447A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN104125804B (zh) | 超声波观测装置与超声波观测装置的动作方法 | |
CN104582584B (zh) | 超声波观测装置以及超声波观测装置的动作方法 | |
CN105828726A (zh) | 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序 | |
US20130096429A1 (en) | Ultrasonic diagnosis apparatus, operation method of the same, and computer readable recording medium | |
CN103717146B (zh) | 超声波观测装置、超声波观测装置的动作方法 | |
CN104411250A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN105392427A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
JPWO2012011414A1 (ja) | 超音波診断装置 | |
CN103945773A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN113749690A (zh) | 血管的血流量测量方法、装置及存储介质 | |
CN108366784A (zh) | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20151118 Address after: Tokyo, Japan, Japan Patentee after: Olympus Corporation Address before: Tokyo, Japan, Japan Patentee before: Olympus Medical Systems Corp. |