CN103153195A - 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 - Google Patents
超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 Download PDFInfo
- Publication number
- CN103153195A CN103153195A CN2011800487962A CN201180048796A CN103153195A CN 103153195 A CN103153195 A CN 103153195A CN 2011800487962 A CN2011800487962 A CN 2011800487962A CN 201180048796 A CN201180048796 A CN 201180048796A CN 103153195 A CN103153195 A CN 103153195A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- characteristic quantity
- correction
- corpse
- amplification
- 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
-
- 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/52098—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 related to workflow protocols
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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/5292—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves using additional data, e.g. patient information, image labeling, acquisition parameters
-
- 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/52033—Gain control of receivers
-
- 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
- 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/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/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physiology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
超声波观测装置以与接收深度相应的放大率对从检体接收到的超声波信号进行放大,生成将放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据,对放大后的超声波信号进行放大校正以使不管接收深度如何放大率均固定,通过对放大校正后的超声波信号的频率进行分析来算出频谱,通过对算出的频谱进行近似来提取频谱的特征量,对频谱和特征量中的任一个进行校正来消减随着超声波传播而产生的强度衰减的贡献,生成用于显示与特征量对应的视觉信息的特征量图像数据。
Description
技术领域
本发明涉及一种使用超声波来观测检体的组织的超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序。
背景技术
以往,作为使用超声波进行的乳腺癌等的检查技术,已知一种超声波弹性成像技术(例如参照专利文献1)。超声波弹性成像是利用生物体内的癌、肿瘤组织的硬度由于疾病的发展状况、生物体不同而不同的技术。在该技术中,在从外部挤压检查部位的状态下,使用超声波对该检查部位的生物体组织的应变量、弹性率进行测量,将该测量结果变换为断层图像来进行图像显示。
在观测超声波时,为了对随着超声波传播而产生的衰减进行校正而以最佳范围进行A/D变换,通常对B模式图像生成用的接收数据(模拟信号)进行基于深位置的接收数据提供高放大率的STC(Sensitivity Time Control:灵敏度时间控制)校正(例如,参照专利文献2)。
专利文献1:国际公开第2005/122906号
专利文献2:日本特开平10-216143号公报
发明内容
发明要解决的问题
然而,在上述STC校正中,仅使模拟信号波形的振幅在全频带中均匀地放大,没有起到对依赖于频率的衰减进行校正的效果。因此,在生成利用超声波振幅的B模式图像时,通过进行STC校正能够得到充分的效果,但是在根据频谱来鉴别检体的组织性状时,存在无法正确地排除随着超声波传播而产生的衰减的影响而鉴别精度降低的问题。
为了解决该问题,考虑在生成B模式图像时输出实施了STC校正的接收信号,另一方面在生成基于频谱的图像时,进行与用于生成B模式图像的发送不同的新发送,输出没有实施STC校正的接收信号。然而,在该情况下存在根据接收信号生成的图像数据的帧频降低这种问题。
本发明是鉴于上述情况而完成的,目的在于提供一种能够正确地排除随着超声波传播而产生的衰减的影响并且防止根据接收到的超声波生成的图像数据的帧频降低的超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序。
用于解决问题的方案
为了解决上述问题而达到目的,本发明所涉及的超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该超声波观测装置的特征在于,具备:信号放大部,其以与接收深度相应的放大率来对从上述检体接收到的超声波信号进行放大;B模式图像数据生成部,其生成将上述信号放大部放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据;放大校正部,其对上述信号放大部放大后的上述超声波信号进行放大校正以使不管接收深度如何放大率均固定;频率分析部,其通过对上述放大校正部放大校正后的上述超声波信号的频率进行分析来算出频谱;特征量提取部,其对上述频率分析部算出的频谱进行衰减校正处理和近似处理来提取上述检体的特征量,其中,上述衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献;以及特征量图像数据生成部,其生成特征量图像数据,该特征量图像数据用于显示与上述特征量提取部提取出的特征量对应的视觉信息。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述信号放大部进行放大时的放大率单调增加直到规定的接收深度为止。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述特征量提取部具有:近似部,其对上述频率分析部算出的频谱进行上述近似处理,由此提取进行上述衰减校正处理前的校正前特征量;以及衰减校正部,其对上述近似部提取出的校正前特征量进行上述衰减校正处理,由此提取上述频谱的特征量。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述特征量提取部具有:衰减校正部,其对上述频谱进行上述衰减校正处理;以及近似部,其对上述衰减校正部校正后的频谱进行上述近似处理,由此提取上述频谱的特征量。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,超声波的接收深度越大,上述衰减校正部进行越大的校正。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,还具备一并进行上述放大校正部和上述衰减校正部的校正的控制部。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述近似部通过回归分析以多项式来对上述频谱进行近似。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述近似部以一次式对上述频谱进行近似并提取多个特征量,上述多个特征量包含上述一次式的斜率、上述一次式的截距以及强度中的至少两个,该强度是使用上述斜率、上述截距以及包含在上述频谱的频带内的特定的频率而决定的。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,还具备:存储部,其将根据由多个已知检体分别反射的超声波而提取出的频谱的特征量与上述多个已知检体的组织性状相关联地存储;以及组织性状判断部,其通过使用上述存储部与上述多个已知检体相关联地存储的特征量以及上述特征量提取部提取出的特征量,来判断上述检体的规定区域的组织性状。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述特征量提取部提取多个特征量,上述存储部存储按每个组织性状对多个已知检体进行分类而得到的组中的各特征量的平均,上述组织性状判断部设定以上述多个特征量中的至少一个为成分的特征量空间,根据检体点与已知检体平均点在上述特征量空间上的距离,来判断上述检体的组织性状,其中,该检体点具有上述检体的频谱的特征量中的形成上述特征量空间成分的特征量作为上述特征量空间的坐标,该已知检体平均点具有上述多个已知检体的上述组中的各特征量中的形成上述特征量空间成分的特征量的平均作为上述特征量空间的坐标。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述组织性状判断部算出将上述多个已知检体中的按每个组织性状来分类的组加上上述检体的特征量而得到的母集团中的特征量的标准偏差,将具有该标准偏差与上述组中的特征量的标准偏差的差最小的特征量的组所对应的组织性状设为上述检体的组织性状。
另外,本发明所涉及的超声波观测装置的特征在于,在上述发明中,上述视觉信息为构成颜色空间的变量。
另外,本发明所涉及的超声波观测装置的动作方法,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作方法的特征在于,具有以下步骤:信号放大步骤,由信号放大部按每个接收深度以从存储部读出的放大率来对从上述检体接收到的超声波信号进行放大,其中,上述存储部存储与接收深度相应的放大率信息;B模式图像数据生成步骤,由B模式图像数据生成部生成将通过上述信号放大步骤放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据;放大校正步骤,由放大校正部对通过上述信号放大步骤放大后的上述超声波信号进行放大校正以使不管接收深度如何放大率均固定;频率分析步骤,由频率分析部对通过上述放大校正步骤放大校正后的上述超声波信号的频率进行分析来算出频谱;特征量提取步骤,由特征量提取部对上述频率分析部算出的频谱进行衰减校正处理和近似处理来提取上述检体的特征量,其中,上述衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献;以及特征量图像数据生成步骤,由特征量图像数据生成部生成特征量图像数据,该特征量图像数据用于显示与通过上述特征量提取步骤提取出的特征量对应的视觉信息。
另外,本发明所涉及的超声波观测装置的动作程序,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作程序的特征在于,使该超声波观测装置执行以下步骤:信号放大步骤,由信号放大部按每个接收深度以从存储部读出的放大率来对从上述检体接收到的超声波信号进行放大,其中,上述存储部存储与接收深度相应的放大率信息;B模式图像数据生成步骤,由B模式图像数据生成部生成将通过上述信号放大步骤放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据;放大校正步骤,由放大校正部对通过上述信号放大步骤放大后的上述超声波信号进行放大校正以使不管接收深度如何放大率均固定;频率分析步骤,由频率分析部对通过上述放大校正步骤放大校正后的上述超声波信号的频率进行分析来算出频谱;特征量提取步骤,由特征量提取部对上述频率分析部算出的频谱进行衰减校正处理和近似处理来提取上述检体的特征量,其中,上述衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献;以及特征量图像数据生成步骤,由特征量图像数据生成部生成特征量图像数据,该特征量图像数据用于显示与通过上述特征量提取步骤提取出的特征量对应的视觉信息。
发明的效果
根据本发明,根据施加了以与接收深度相应的放大率进行放大的STC校正的信号来生成B模式图像数据,另一方面,抵消STC校正的影响而进行将放大率设为谱固定的放大校正之后对频谱进行衰减校正处理和近似处理,由此提取检体的特征量,生成用于显示与提取的特征量对应的视觉信息的特征量图像数据,因此不仅在特征量图像数据中排除随着超声波传播而产生的衰减的影响,还不需要分开发送B模式图像用的信号和特征量图像用的信号,其中,衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献。因而,正确地排除随着超声波传播而产生的衰减的影响,并且能够防止根据接收到的超声波生成的图像数据的帧频降低。
附图说明
图1是表示本发明的实施方式1所涉及的超声波观测装置的结构的框图。
图2是表示本发明的实施方式1所涉及的超声波观测装置的信号放大部所进行的放大处理中的接收深度与放大率的关系的图。
图3是表示本发明的实施方式1所涉及的超声波观测装置的放大校正部所进行的放大处理中的接收深度与放大率的关系的图。
图4是表示本发明的实施方式1所涉及的超声波观测装置的处理的概要的流程图。
图5是表示本发明的实施方式1所涉及的超声波观测装置的显示部中的B模式图像的显示例的图。
图6是表示本发明的实施方式1所涉及的超声波观测装置的频率分析部所进行的处理的概要的流程图。
图7是示意性地表示一个声线的数据排列的图。
图8是表示本发明的实施方式1所涉及的超声波观测装置的频率分析部所算出的频谱的例子(第一例)的图。
图9是表示本发明的实施方式1所涉及的超声波观测装置的频率分析部所算出的频谱的例子(第二例)的图。
图10是表示根据对与图8示出的直线有关的特征量进行衰减校正之后的特征量而决定的新的直线的图。
图11是表示本发明的实施方式1所涉及的超声波观测装置的组织性状判断部所进行的处理的概要的流程图。
图12是表示本发明的实施方式1所涉及的超声波观测装置的组织性状判断部所设定的特征量空间的一例的图。
图13是表示本发明的实施方式1所涉及的超声波观测装置的显示部所显示的判断结果显示图像的显示例的图。
图14是说明本发明的实施方式1所涉及的超声波观测装置所进行的衰减校正处理的效果的图。
图15是表示本发明的实施方式2所涉及的超声波观测装置的处理的概要的流程图。
图16是示意性地表示本发明的实施方式2所涉及的超声波观测装置所进行的衰减校正处理的概要的图。
具体实施方式
下面,参照附图说明用于实施本发明的方式(以下称为“实施方式”)。
(实施方式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具有信号放大部31,该信号放大部31放大回波信号。具体地说,信号放大部31进行STC校正,接收深度越大的回波信号以越高放大率进行放大。图2是表示回波信号的接收深度与放大率的关系的图。图2示出的接收深度z是根据从超声波的接收开始时间点起的经过时间而算出的量。如图2所示,在接收深度z小于阈值zth的情况下,放大率β随着接收深度z的增加而从β0向βth(>β0)线性地增加。另外,在接收深度z为阈值zth以上的情况下,放大率β取固定值βth。阈值zth的值是从检体接收的超声波信号几乎衰减而噪声处于支配地位的值。此外,更一般地,只要在接收深度z小于阈值zth的情况下放大率β随着接收深度z的增加而单调增加即可。
发送和接收部3在对由信号放大部31放大的回波信号实施过滤等处理之后,通过进行A/D变换而生成数字RF信号并输出。此外,在超声波探头2是使多个超声波振子电扫描的超声波探头的情况下,发送和接收部3具有与多个超声波振子对应的束合成用的多通道电路。
运算部4具有:放大校正部41,其对由发送和接收部3输出的数字RF信号进行放大校正以使不管接收深度如何放大率均固定;频率分析部42,其通过对进行了放大校正的数字RF信号实施高速傅里叶变换(FFT)来进行频率分析;特征量提取部43,其对由频率分析部42算出的频谱(功率谱)进行近似处理以及校正处理,由此提取频谱的特征量,其中,在该校正处理中削减依赖于超声波的接收深度和频率的超声波的衰减的贡献;以及组织性状判断部44,其使用由特征量提取部43提取的特征量来判断检体的规定区域的组织性状。
图3是表示放大校正部41所进行的放大处理中的接收深度与放大率的关系的图。如图3所示,关于放大校正部41所进行的放大处理的放大率,在接收深度为零时是最大值βth-β0,到接收深度zth为止线性地减少,在接收深度zth以上为零。通过利用这种放大率对数字RF信号进行放大处理,由此抵消信号放大部31中的STC校正的影响,能够输出固定放大率βth的信号。此外,放大校正部41所进行的接收深度z与放大率β的关系当然根据信号放大部31中的接收深度与放大率的关系不同而不同。
频率分析部42针对各声线(行数据),对由规定的数据量构成的FFT数据群进行高速傅里叶变换,由此算出频谱。频谱示出根据检体的组织性状不同而不同的趋势。这是由于,频谱与作为使超声波散射的散射体的检体的大小、密度、声阻抗等具有相关性。
特征量提取部43具有:近似部431,其对由频率分析部42算出的频谱进行近似处理,由此算出进行衰减校正处理前的校正前特征量;以及衰减校正部432,其对由近似部431进行近似而得到的校正前特征量进行衰减校正处理,由此提取特征量。
近似部431通过回归分析以一次式对频谱进行近似,由此提取使该近似的一次式具有特征的校正前特征量。具体地说,近似部431通过回归分析来算出一次式的斜率a0和截距b0,并且算出频谱中的频带内的特定频率处的强度作为校正前特征量。在本实施方式1中,近似部431算出中心频率fMID=(fLOW+fHIGH)/2处的强度(Mid-band fit)c0=a0fMID+b0,但是这严格说来仅是一例。在此所指的“强度”是指电压、电力、声压、声能等参数中的任一个。
三个特征量中的斜率a0与超声波散射体的大小具有相关性,通常认为散射体越大则斜率具有越小的值。另外,截距b0与散射体的大小、声阻抗的差、散射体的密度(浓度)等具有相关性。具体地说,认为对于截距b0,散射体越大则具有越大的值,声阻抗越大则具有越大的值,散射体的密度(浓度)越大则具有越大的值。中心频率fMID处的强度(以下简单称为“强度”)c0是根据斜率a0和截距b0导出的间接的参数,提供有效频带内的中心处的频谱强度。因此,认为强度c0除了与散射体的大小、声阻抗的差、散射体的密度相关以外,还与B模式图像的亮度具有某种程度的相关性。此外,由特征量提取部43算出的近似多项式并不限定于一次式,还能够使用二次以上的近似多项式。
说明由衰减校正部432进行的校正。超声波的衰减量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的输入来变更衰减率α的值。
衰减校正部432如下那样对由近似部431提取出的校正前特征量(斜率a0、截距b0、强度c0)进行校正。
a=a0+2αz···(2)
b=b0···(3)
c=c0+2αzfMID(=afMID+b)···(4)
根据式(2)、(4)可知,衰减校正部432进行如下校正:超声波的接收深度z越大则校正量越大。另外,根据式(3),与截距相关的校正是恒等变换。这是由于,截距是与频率0(Hz)对应的频率成分而没有被衰减。
组织性状判断部44按照每个特征量来算出由特征量提取部43提取出并由衰减校正部432校正后的频谱的特征量的平均和标准偏差。组织性状判断部44使用算出的平均和标准偏差以及由存储部8存储的已知检体的频谱的特征量的平均和标准偏差来判断检体规定区域的组织性状。在此所指的“规定区域”是观察到由图像处理部5生成的图像的超声波观测装置1的操作者通过输入部6指定的图像中的区域(以下称为“关心区域”)。另外,在此所指的“组织性状”例如是指癌症、内分泌肿瘤、粘液性肿瘤、正常组织、脉管等中的任一个。此外,在检体为胰腺的情况下,作为组织性状还包含慢性胰腺炎、自身免疫性胰腺炎等。
由组织性状判断部44算出的特征量的平均和标准偏差反应了核肿大、异形等细胞水平的变化、间质中的纤维增生、实质组织被纤维替换等组织变化,与组织性状相应地呈现特有的值。因而,通过使用这种特征量的平均和标准偏差,能够正确地判断检体的规定区域的组织性状。
图像处理部5具有:B模式图像数据生成部51,其根据回波信号来生成B模式图像数据;特征量图像数据生成部52,其生成特征量图像数据,该特征量图像数据显示与检体的特征量对应的视觉信息;以及判断结果显示图像数据生成部53,其使用由B模式图像数据生成部51和运算部4分别输出的数据来生成用于显示关心区域的组织性状的判断结果和与该判断结果有关的信息的判断结果显示图像数据。
B模式图像数据生成部51对数字信号进行带通滤波、对数变换、增益处理、对比度处理等使用公知技术的信号处理,并且进行数据的间隔剔除等,由此生成B模式图像数据,其中,与根据显示部7中的图像的显示范围而决定的数据步长相应地进行数据的间隔剔除。
特征量图像数据生成部52生成特征量图像数据时使用的视觉信息例如是构成以R(红)、G(绿)、B(蓝)为变量的RGB表色系等颜色空间的变量。
判断结果显示图像数据生成部53生成判断结果显示图像数据,该判断结果显示图像数据包括由B模式图像数据生成部51生成的B模式图像数据、由特征量图像数据生成部52生成的特征量图像数据、由特征量提取部43提取的特征量以及由组织性状判断部44判断的判断结果。
存储部8具有:已知检体信息存储部81,其存储已知检体的信息;放大率信息存储部82,其存储信号放大部31和放大校正部41进行放大处理时参照的放大率的信息;窗函数存储部83,其存储频率分析部42进行频率分析处理时使用的窗函数;以及校正信息存储部84,其存储衰减校正部432进行处理时参照的校正信息。
已知检体信息存储部81将针对已知检体提取出的频谱的特征量与已知检体的组织性状相关联地进行存储。另外,已知检体信息存储部81针对与已知检体相关联的频谱的特征量,将按根据已知检体的组织性状分类的每个组算出的平均和标准偏差与已知检体的特征量的全部数据一起进行存储。在此,已知检体的特征量是通过与本实施方式1相同的处理提取出的。其中,不需要利用超声波观测装置1进行已知检体的特征量提取处理。由已知检体信息存储部81存储的已知检体的信息期望是与组织性状有关的可靠性高的信息。放大率信息存储部82存储图2和图3示出的放大率与接收深度的关系。窗函数存储部83存储Hamming、Hanning、Blackman等窗函数中的至少一个窗函数。校正信息存储部84存储与式(2)~(4)的变换有关的信息。
存储部8是使用预先存储了本实施方式1所涉及的超声波观测装置的动作程序、启动规定的OS的程序等的ROM以及存储各处理的运算参数、数据等的RAM等来实现的。
具有上述功能结构的超声波观测装置1的除了超声波探头2以外的结构要素是使用具备具有运算和控制功能的CPU的计算机来实现的。超声波观测装置1所具备的CPU从存储部8读出由存储部8存储、保存的信息以及包括上述超声波观测装置的动作程序的各种程序,由此执行与本实施方式1所涉及的超声波观测装置的动作方法有关的运算处理。
此外,本实施方式1所涉及的超声波观测装置的动作程序记录到硬盘、快闪存储器、CD-ROM、DVD-ROM、软盘等计算机可读取的记录介质而广泛流通。
图4是表示具有以上结构的超声波观测装置1的处理的概要的流程图。在图4中,超声波观测装置1首先通过超声波探头2进行新检体的测量(步骤S1)。
接着,从超声波探头2接收到回波信号的信号放大部31进行该回波信号的放大(步骤S2)。在此,信号放大部31根据图2示出的放大率与接收深度的关系来进行放大。
之后,B模式图像数据生成部51使用从发送和接收部3输出的B模式图像用回波信号来生成B模式图像数据(步骤S3)。
接着,控制部9进行控制,使显示部7显示与由B模式图像数据生成部51生成的B模式图像数据对应的B模式图像(步骤S4)。图5是表示显示部7中的B模式图像的显示例的图。该图示出的B模式图像100是使采用RGB表色系作为颜色空间的情况下的变量即R(红)、G(绿)、B(蓝)的值一致的灰度等级图像。
之后,在通过输入部6设定了关心区域的情况下(步骤S5:“是”),放大校正部41对从发送和接收部3输出的信号进行校正以使不管接收深度如何放大率均固定(步骤S6)。在此,放大校正部41根据图3示出的放大率与接收深度的关系来进行放大处理。
另一方面,在没有设定关心区域的情况下(步骤S5:“否”),在由输入部6输入结束处理的指示时(步骤S7:“是”),超声波观测装置1结束处理。与此相对,在没有设定关心区域的情况下(步骤S5:“否”),在没有由输入部6输入结束处理的指示时(步骤S7:“否”),超声波观测装置1返回到步骤S5。
在步骤S6之后,频率分析部42通过FFT运算进行频率分析,由此算出频谱(步骤S8)。在该步骤S8中,还能够将图像的全部区域设定为关心区域。
在此,参照图6示出的流程图来详细说明频率分析部42进行的处理(步骤S8)。首先,频率分析部42将最初成为分析对象的声线的声线编号L设为初始值L0(步骤S21)。例如可以对发送和接收部3最初接收到的声线附加初始值L0,也可以对与通过输入部6设定的关心区域的左右一侧的边界位置对应的声线附加初始值L0。
接着,频率分析部42算出设定于一个声线上的多个数据位置的全部位置的频谱。首先,频率分析部42设定代表用于进行FFT运算而获取的一系列数据群(FFT数据群)的数据位置Z(相当于接收深度)的初始值Z0(步骤S22)。图7是示意性地表示一个声线的数据排列的图。在该图示出的声线LD中,白色或者黑色的长方形意味着一个数据。以与发送和接收部3进行的A/D变换中的采样频率(例如50MHz)对应的时间间隔使声线LD离散化。在图7中,示出将声线LD的第一个数据设定为数据位置Z的初始值Z0的情况。此外,图7严格说来仅是一例,能够任意地设定初始值Z0的位置。例如,也可以将与关心区域的上端位置对应的数据位置Z设定为初始值Z0。
之后,频率分析部42获取数据位置Z的FFT数据群(步骤S23),使由窗函数存储部83存储的窗函数作用于获取到的FFT数据群(步骤S24)。这样使窗函数作用于FFT数据群,由此避免FFT数据群在边界处不连续,从而能够防止产生伪像。
接着,频率分析部42判断数据位置Z的FFT数据群是否为正常的数据群(步骤S25)。在此,FFT数据群需要具有2的乘方个数据数。以下,将FFT数据群的数据数设为2n(n为正整数)。FFT数据群为正常是指数据位置Z在FFT数据群中位于从前方起第2n-1的位置。换言之,FFT数据群正常是指在数据位置Z的前方存在2n-1-1(设为=N)个数据而在数据位置Z的后方存在2n-1(设为=M)个数据。在图7示出的情况下,FFT数据群F2、F3、FK-1正常,另一方面FFT数据群F1、FK异常。其中,在图7中,设为n=4(N=7、M=8)。
在步骤S25中的判断的结果是数据位置Z的FFT数据群正常的情况下(步骤S25:“是”),频率分析部42转移到后述的步骤S27。
在步骤S25中的判断的结果是数据位置Z的FFT数据群不正常的情况下(步骤S25:“否”),频率分析部42对不足部分插入零数据,由此来生成正常的FFT数据群(步骤S26)。在追加零数据之前使窗函数作用于在步骤S25中判断为不正常的FFT数据群。因此,即使对FFT数据群插入零数据,也不产生数据的不连续。在步骤S26之后,频率分析部42转移到后述的步骤S27。
在步骤S27中,频率分析部42使用FFT数据群来进行FFT运算,由此得到频谱(步骤S27)。图8和图9是表示由频率分析部42算出的频谱例的图。在图8和图9中,横轴f是频率,纵轴I为强度。在图8和图9分别示出的频谱曲线C1和C2中,频谱的下限频率fLOW和上限频率fHIGH是根据超声波探头2的频带、由发送和接收部3发送的脉冲信号的频带等而决定的参数,例如是fLOW=3MHz、fHIGH=10MHz。此外,通过后述的特征量提取处理来说明图8示出的直线L1和图9示出的直线L2。另外,在本实施方式1中,曲线和直线由离散的点的集合构成。这一点在后述的实施方式中也相同。
接着,频率分析部42对数据位置Z加上规定的数据步长D而算出下一个分析对象的FFT数据群的数据位置Z(步骤S28)。期望此处的数据步长D与由B模式图像数据生成部51生成B模式图像数据时利用的数据步长一致,但是在想要削减频率分析部42中的运算量的情况下,也可以设定比B模式图像数据生成部51所利用的数据步长大的值。在图7中,示出D=15的情况。
之后,频率分析部42判断数据位置Z是否大于最终数据位置Zmax(步骤S29)。在此,最终数据位置Zmax可以是声线LD的数据长度,也可以是与关心区域的下端对应的数据位置。在判断的结果是数据位置Z大于最终数据位置Zmax的情况下(步骤S29:“是”),频率分析部42使声线编号L增加1(步骤S30)。另一方面,在数据位置Z为最终数据位置Zmax以下的情况下(步骤S29:“否”),频率分析部42返回到步骤S23。这样,频率分析部42针对一个声线LD,对[{(Zmax-Z0)/D}+1](=K)个FFT数据群进行FFT运算。在此,[X]表示不超过X的最大整数。
在通过步骤S30增加之后的声线编号L大于最终声线编号Lmax的情况下(步骤S31:“是”),频率分析部42返回到图4示出的主例程。另一方面,在通过步骤S30增加之后的声线编号L为最终声线编号Lmax以下的情况下(步骤S31:“否”),频率分析部42返回到步骤S22。
这样,频率分析部42对(Lmax-L0+1)个声线各自进行K次的FFT运算。此外,例如可以对由发送和接收部3接收到的最终的声线附加最终声线编号Lmax,也可以对与关心区域左右的任一个边界对应的声线附加最终声线编号Lmax。以下,将频率分析部42对全部声线进行的FFT运算的总数(Lmax-L0+1)×K设为P。
接在上述说明的步骤S8的频率分析处理之后,作为近似处理,近似部431对频率分析部42算出的P个频谱进行回归分析,由此提取校正前特征量(步骤S9)。具体地说,近似部431通过回归分析来算出对频带fLOW<f<fHIGH的频谱进行近似的一次式,由此作为校正前特征量而提取出使该一次式具有特征的斜率a0、截距b0、强度c0。图8示出的直线L1和图9示出的直线L2是在该步骤S9中对频谱曲线C1和C2分别进行特征量提取处理而得到的回归直线。
之后,衰减校正部432对由近似部431提取出的校正前特征量进行衰减校正处理(步骤S10)。例如在数据的采样频率为50MHz的情况下,数据采样的时间间隔为20(nsec)。在此,当将声速设为1530(m/sec)时,数据采样距离间隔成为1530(m/sec)×20(nsec)/2=0.0153(mm)。当将从声线LD的第一个数据起到处理对象的FFT数据群的数据位置为止的数据步长数设为k时,该数据位置Z成为0.0153k(mm)。衰减校正部432将这样求得的数据位置Z的值代入到上述式(2)~(4)的接收深度z,由此算出作为频谱的特征量的斜率a、截距b、强度c。图10是表示根据对与图8示出的直线L1相关联的特征量进行衰减校正之后的特征量而决定的直线的图。表示图10示出的直线L1’的式如下。
I=af+b=(a0+2αZ)f+b0···(5)
根据该式(5)也可知,与直线L1相比,直线L1’的斜率大且截距的值相同。
之后,组织性状判断部44根据由特征量提取部43提取出的、由衰减校正部432校正的特征量以及由已知检体信息存储部81存储的已知检体信息来判断检体的关心区域中的组织性状(步骤S11)。
在此,参照图11示出的流程图来详细说明由组织性状判断部44进行的处理(步骤S11)。首先,组织性状判断部44设定判断组织性状时使用的特征量空间(步骤S41)。在本实施方式1中,作为三个特征量的斜率a、截距b、强度c中的独立的参数为两个。因而,能够将以三个特征量中的任意两个特征量为成分的二维空间设定为特征量空间。另外,还能够将以三个特征量中的任意一个特征量为成分的一维空间设定为特征量空间。在该步骤S41中,设为要设定的特征量空间为预先决定的特征量空间,但是也可以由操作者通过输入部6来选择期望的特征量空间。
图12是表示组织性状判断部44所设定的特征量空间的一例的图。在图12示出的特征量空间中,横轴为截距b,纵轴为强度c。图12示出的点Sp示出对判断对象的检体算出的具有截距b和强度c作为特征量空间的坐标的点(以下,将该点称为“检体点”)。另外,图12示出的区域Gμ、Gν、Gρ表示由已知检体信息存储部81存储的已知检体的组织性状分别为μ、ν、ρ的组。在图12示出的情况下,三个组Gμ、Gν、Gρ在特征量空间上存在于与其它组相互不交叉的区域内。
在本实施方式1中,在求已知检体的特征量时也将对通过频率分析得到的频谱的校正前特征量进行衰减校正而得到的特征量作为指标而进行组织性状的分类、判断,因此能够鉴别相互不同的组织性状。特别是,在本实施方式1中,使用进行了衰减校正的特征量,因此与不进行衰减校正而使用提取出的特征量的情况相比,能够在更明确地分离的状态下得到特征量空间内的各组织性状的区域。
在步骤S41之后,组织性状判断部44分别算出检体点Sp与具有分别包含在组Gμ、Gν、Gρ中的FFT数据群的频谱的截距b和强度c的各平均作为特征量空间的坐标的点μ0、ν0、ρ0(以下,将这些点称为“已知检体平均点”)之间的特征量空间上的距离dμ、dν、dρ(步骤S42)。在此,在特征量空间内的b轴成分与c轴成分的刻度大不相同的情况下,期望适当地进行加权以使各距离的贡献大致均等。
接着,组织性状判断部44根据通过步骤S42算出的距离来判断包括检体点Sp的全部检体点的组织性状(步骤S43)。例如在图12示出的情况下,由于距离dμ为最小,因此组织性状判断部44判断为检体的组织性状为μ。此外,在检体点Sp与已知检体平均点μ0、ν0、ρ0极端分离的情况下,即使求出距离dμ、dν、dρ的最小值,组织性状的判断结果的可靠性也低。因此,在dμ、dν、dρ大于规定阈值的情况下,组织性状判断部44也可以输出错误信号。另外,在dμ、dν、dρ中的最小值产生两个以上的情况下,组织性状判断部44可以选择与最小值对应的所有组织性状作为候选,也可以按照规定的规则来选择任一个组织性状。在后者的情况下,例如能够举出将癌症等恶性高的组织性状的优先级设定得高的方法。另外,在dμ、dν、dρ中最小值产生两个以上的情况下,组织性状判断部44也可以输出错误信号。
之后,组织性状判断部44输出步骤S42中的距离计算结果以及步骤S43中的判断结果(步骤S44)。由此,结束步骤S11的组织性状判断处理。
在以上说明的步骤S11之后,特征量图像数据生成部52生成特征量图像数据,该特征量图像数据用于显示与由特征量提取部43提取出的特征量对应的视觉信息(步骤S12)。
接着,判断结果显示图像数据生成部53通过使用由B模式图像数据生成部51生成的B模式图像数据、由特征量提取部43算出的特征量、由特征量图像数据生成部52生成的特征量图像数据以及由组织性状判断部44判断的判断结果,来生成判断结果显示图像数据(步骤S13)。
之后,显示部7显示由判断结果显示图像数据生成部53生成的判断结果显示图像(步骤S14)。图13是表示显示部7所显示的判断结果显示图像的显示例的图。该图示出的判断结果显示图像200具有信息显示部201和图像显示部202,该信息显示部201显示包括组织性状的判断结果的各种关联信息,该图像显示部202显示根据B模式图像来显示特征量的特征量图像。
在信息显示部201中例如显示检体的识别信息(ID编号、姓名、性别等)、由组织性状判断部44算出的组织性状判断结果、与进行组织性状判断时的特征量有关的信息、增益、对比度等超声波画质信息。在此,作为与特征量有关的信息,能够利用位于关心区域内部的Q组FFT数据群的频谱的特征量的平均、标准偏差进行显示。具体地说,在信息显示部201中,例如能够显示为斜率=1.5±0.3(dB/MHz)、截距=-60±2(dB)、强度=-50±1.5(dB)。
在图像显示部202中显示的特征量图像300是针对图5示出的B模式图像100将截距b均等地分配给R(红)、G(绿)、B(蓝)的灰度等级图像。
显示部7显示具有以上结构的判断结果显示图像200,由此操作者能够更正确地掌握关心区域的组织性状。此外,判断结果显示图像并不限定于上述结构。例如,作为判断结果显示图像,也可以排列显示特征量图像和B模式图像。由此,能够在一个画面上识别两个图像的差异。
图14是说明超声波观测装置1所进行的衰减校正处理的效果的图。图14示出的图像400是不进行衰减校正的情况下的特征量图像。对于特征量图像400,在接收深度大的区域(图的下方区域)由于衰减的影响而信号强度下降,图像变暗。与此相对,可知在进行了衰减校正的特征量图像300中,得到画面整体明度均匀的图像。
此外,图13、14示出的特征量图像300严格说来仅是一例。除此以外,例如还能够将三个特征量a’、b’、c’分别分配至R、G、B,通过彩色图像来显示特征量图像。在该情况下,以固有的颜色对组织性状进行表现,因此操作者能够根据图像的颜色分布来掌握关心区域的组织性状。另外,代替以RGB表色系构成颜色空间,也可以通过青色、品红、黄色那样的补色系的变量来构成颜色空间,对各变量分配特征量。另外,也可以按照规定比率来混合B模式图像数据与彩色图像数据,由此生成特征量图像数据。另外,也可以仅将关心区域替换为彩色图像数据来生成特征量图像数据。
根据以上说明的本发明的实施方式1,根据施加了以与接收深度相应的放大率进行放大的STC校正的信号来生成B模式图像数据,另一方面,在抵消STC校正的影响而进行使放大率谱固定的放大校正之后算出频谱并提取校正前特征量,对提取出的校正前特征量实施衰减校正,由此提取检体的特征量,生成用于显示与提取出的特征量对应的视觉信息的特征量图像数据,因此不仅在特征量图像数据中排除随着超声波传播而产生的衰减的影响,也不需要分开发送B模式图像用的信号和特征量图像用的信号。因而,正确地排除随着超声波传播而产生的衰减的影响,并且能够防止根据接收到的超声波生成的图像数据的帧频降低。
另外,根据本实施方式1,使用正确地进行衰减校正后的频谱的特征量来判断检体的规定区域的组织性状,因此不使用生体组织的应变量、弹性率而能够明确地鉴别组织的差异。因而,能够高精度地鉴别组织性状,并且能够提高观测结果的可靠性。
此外,作为本实施方式1的变形例,控制部9也可以一并进行放大校正部41的放大校正处理以及衰减校正部432的衰减校正处理。该处理相当于删除图4的步骤S6中的放大处理并且如以下式(6)那样变更图4的步骤S10中的衰减校正处理的衰减量。
A’=2αzf+γ(z)···(6)
在此,右边的γ(z)为接收深度z处的放大率β与β0的差,表示为如下。
γ(z)=-{(βth-β0)/zth}z+βth-β0(z≤zth)···(7)
γ(z)=0(z>zth)···(8)
(实施方式2)
在本发明的实施方式2中,特征量提取部进行的特征量提取与实施方式1不同。本实施方式2所涉及的超声波观测装置的结构与在实施方式1中说明的超声波观测装置1的结构相同。因此,在以下说明中,对与超声波观测装置1的结构要素对应的结构要素附加相同的附图标记。
在本实施方式2的特征量提取处理中,首先衰减校正部432对由频率分析部42算出的频谱进行衰减校正处理。之后,近似部431对由衰减校正部432进行衰减校正后的频谱进行近似处理,由此提取频谱的特征量。
图15是表示本实施方式2所涉及的超声波观测装置的处理的概要的流程图。在图15中,步骤S51~S58的处理与图4的步骤S1~S8的处理依次对应。
在步骤S59中,衰减校正部432对频率分析部42通过FFT运算而算出的频谱进行衰减校正(步骤S59)。图16是示意性地表示该步骤S59的衰减校正处理的概要的图。如图16所示,衰减校正部432针对频谱曲线C3,对全部频率f进行将上述式(1)的衰减量A加到强度I的校正,由此得到新频谱曲线C3’。由此,能够得到削减了随着超声波传播而产生的衰减的贡献的频谱。
之后,近似部431对由衰减校正部432进行衰减校正后的全部频谱进行回归分析,由此提取频谱的特征量(步骤S60)。具体地说,近似部431通过回归分析来算出一次式的斜率a、截距b以及中心频率fMID处的强度c。图16示出的直线L3是在该步骤S60中对频谱曲线C3进行特征量提取处理而得到的回归直线(截距b3)。
步骤S61~S64的处理与图4的步骤S11~S14的处理依次对应。
根据上述说明的本发明的实施方式2,根据施加了以与接收深度相应的放大率进行放大的STC校正的信号来生成B模式图像数据,另一方面,在抵消STC校正的影响而进行使放大率谱固定的放大校正之后算出频谱,在对该频谱实施衰减校正之后提取特征量,生成用于显示与提取出的特征量对应的视觉信息的特征量图像数据,因此不仅在特征量图像数据中排除随着超声波传播而产生的衰减的影响,也不需要分开发送B模式图像用的信号和特征量图像用的信号。因而,与上述实施方式1同样地,正确地排除随着超声波传播而产生的衰减的影响,并且能够防止根据接收到的超声波生成的图像数据的帧频降低。
另外,根据本实施方式2,使用正确地进行衰减校正后的频谱的特征量来判断检体的规定区域的组织性状,因此不使用生体组织的应变量、弹性率而能够明确地鉴别组织的差异。因而,能够高精度地鉴别组织性状,并且能够提高测量结果的可靠性。
此外,在本实施方式2中也能够删除图15的步骤S56的放大校正处理,在图15的步骤S59中将进行频谱的衰减校正时的衰减量设为式(6)的A’来进行处理。
(实施方式3)
本发明的实施方式3的组织性状判断部中的组织性状判断处理与实施方式1不同。本实施方式3所涉及的超声波观测装置的结构与实施方式1所说明的超声波观测装置1的结构相同。因此,在以下说明中,对与超声波观测装置1的结构要素对应的结构要素附加相同的附图标记。
组织性状判断部44在将特征量(a、b、c)分别加到构成组织性状μ、ν、ρ的组Gμ、Gν、Gρ(参照图12)而构成新的母集团之后,求出构成各组织性状的每个数据的特征量的标准偏差。
之后,组织性状判断部44算出仅由已知检体构成的原来的母集团中的组Gμ、Gν、Gρ的各特征量的标准偏差与分别加上新检体的新的母集团中的组Gμ、Gν、Gρ的各特征量的标准偏差的差(以下,简单称为“标准偏差的差”),将包括该标准偏差的差最小的特征量的组所对应的组织性状判断为检体的组织性状。
在此,组织性状判断部44还可以仅对从多个特征量中预先选择的特征量的标准偏差算出标准偏差的差。该情况下的特征量的选择可以由操作者任意地进行,也可以由超声波观测装置1自动地进行。
另外,组织性状判断部44也可以算出按每组对全部特征量的标准偏差的差适当地进行加权后相加的值,将该值最小的组所对应的组织性状判断为检体的组织性状。在该情况下,例如在特征量为斜率a、截距b、强度c时,组织性状判断部44将与斜率a、截距b、强度c分别对应的权重设为wa、wb、wc而算出wa×(a的标准偏差的差)+wb×(b的标准偏差的差)+wc×(c的标准偏差的差),根据该算出的值来判断检体的组织性状。此外,权重wa、wb、wc的值可以由操作者任意地设定,也可以由超声波观测装置1自动地设定。
另外,组织性状判断部44也可以算出按每组对全部特征量的标准偏差的差的乘方适当地进行加权后相加的值的平方根,将该平方根最小的组所对应的组织性状判断为检体的组织性状。在该情况下,例如在特征量为斜率a、截距b、强度c时,组织性状判断部44将与斜率a、截距b、强度c分别对应的权重设为w’a、w’b、w’c而算出{w’a×(a的标准偏差的差)2+w’b×(b的标准偏差的差)2+w’c×(c的标准偏差的差)2}1/2,根据算出的该值来判断组织性状。此外,在该情况下也同样,权重w’a、w’b、w’c的值可以由操作者任意地设定,也可以由超声波观测装置1自动地设定。
根据上述说明的本发明的实施方式3,与上述实施方式1同样地,正确地排除随着超声波传播而产生的衰减的影响,并且能够防止根据接收到的超声波生成的图像数据的帧频降低。
另外,根据本实施方式3,能够高精度地鉴别组织性状,并且能够提高测量结果的可靠性,通过消除随着超声波传播而产生的衰减的影响能够进行更加高精度的组织性状判断。
此外,在本实施方式3中,组织性状判断部44根据原来的母集团与加上新检体的母集团之间的各特征量的标准偏差的变化来进行组织性状的判断,但是这仅是一例。例如,组织性状判断部44也可以根据原来的母集团与加上新检体的母集团之间的各特征量的平均的变化来进行组织性状的判断。
到此为止说明了用于实施本发明的方式,但是本发明并不仅限定于上述实施方式1~3。即,本发明在不脱离记载于权利要求范围内的技术思想的范围内可包括各种实施方式。
附图标记说明
1:超声波观测装置;2:超声波探头;3:发送和接收部;4:运算部;5:图像处理部;6:输入部;7:显示部;8:存储部;9:控制部;21:信号变换部;31:信号放大部;41:放大校正部;42:频率分析部;43:特征量提取部;44:组织性状判断部;51:B模式图像数据生成部;52:特征量图像数据生成部;81:已知检体信息存储部;82:放大率信息存储部;83:窗函数存储部;84:校正信息存储部;100:B模式图像;200:判断结果显示图像;201:信息显示部;202:图像显示部;300、400:特征量图像;431:近似部;432:衰减校正部。
Claims (14)
1.一种超声波观测装置,对检体发送超声波并且接收由上述检体反射的超声波,该超声波观测装置的特征在于,具备:
信号放大部,其以与接收深度相应的放大率来对从上述检体接收到的超声波信号进行放大;
B模式图像数据生成部,其生成将上述信号放大部放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据;
放大校正部,其对上述信号放大部放大后的上述超声波信号进行放大校正以使不管接收深度如何放大率均固定;
频率分析部,其通过对上述放大校正部放大校正后的上述超声波信号的频率进行分析来算出频谱;
特征量提取部,其对上述频率分析部算出的频谱进行衰减校正处理和近似处理来提取上述检体的特征量,其中,上述衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献;以及
特征量图像数据生成部,其生成特征量图像数据,该特征量图像数据用于显示与上述特征量提取部提取出的特征量对应的视觉信息。
2.根据权利要求1所述的超声波观测装置,其特征在于,
上述信号放大部进行放大时的放大率单调增加直到规定的接收深度为止。
3.根据权利要求1或者2所述的超声波观测装置,其特征在于,
上述特征量提取部具有:
近似部,其对上述频率分析部算出的频谱进行上述近似处理,由此提取进行上述衰减校正处理前的校正前特征量;以及
衰减校正部,其对上述近似部提取出的校正前特征量进行上述衰减校正处理,由此提取上述频谱的特征量。
4.根据权利要求1或者2所述的超声波观测装置,其特征在于,
上述特征量提取部具有:
衰减校正部,其对上述频谱进行上述衰减校正处理;以及
近似部,其对上述衰减校正部校正后的频谱进行上述近似处理,由此提取上述频谱的特征量。
5.根据权利要求3或者4所述的超声波观测装置,其特征在于,
超声波的接收深度越大,上述衰减校正部进行越大的校正。
6.根据权利要求3~5中的任一项所述的超声波观测装置,其特征在于,
还具备一并进行上述放大校正部和上述衰减校正部的校正的控制部。
7.根据权利要求3~6中的任一项所述的超声波观测装置,其特征在于,
上述近似部通过回归分析以多项式来对上述频谱进行近似。
8.根据权利要求7所述的超声波观测装置,其特征在于,
上述近似部以一次式对上述频谱进行近似并提取多个特征量,上述多个特征量包含上述一次式的斜率、上述一次式的截距以及强度中的至少两个,该强度是使用上述斜率、上述截距以及包含在上述频谱的频带内的特定的频率而决定的。
9.根据权利要求1~8中的任一项所述的超声波观测装置,其特征在于,还具备:
存储部,其将根据由多个已知检体分别反射的超声波而提取出的频谱的特征量与上述多个已知检体的组织性状相关联地存储;以及
组织性状判断部,其通过使用上述存储部与上述多个已知检体相关联地存储的特征量以及上述特征量提取部提取出的特征量,来判断上述检体的规定区域的组织性状。
10.根据权利要求9所述的超声波观测装置,其特征在于,
上述特征量提取部提取多个特征量,
上述存储部存储按每个组织性状对多个已知检体进行分类而得到的组中的各特征量的平均,
上述组织性状判断部设定以上述多个特征量中的至少一个为成分的特征量空间,根据检体点与已知检体平均点在上述特征量空间上的距离,来判断上述检体的组织性状,其中,该检体点具有上述检体的频谱的特征量中的形成上述特征量空间成分的特征量作为上述特征量空间的坐标,该已知检体平均点具有上述多个已知检体的上述组中的各特征量中的形成上述特征量空间成分的特征量的平均作为上述特征量空间的坐标。
11.根据权利要求9所述的超声波观测装置,其特征在于,
上述组织性状判断部算出将上述多个已知检体中的按每个组织性状来分类的组加上上述检体的特征量而得到的母集团中的特征量的标准偏差,将具有该标准偏差与上述组中的特征量的标准偏差的差最小的特征量的组所对应的组织性状设为上述检体的组织性状。
12.根据权利要求1~11中的任一项所述的超声波观测装置,其特征在于,
上述视觉信息为构成颜色空间的变量。
13.一种超声波观测装置的动作方法,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作方法的特征在于,具有以下步骤:
信号放大步骤,由信号放大部按每个接收深度以从存储部读出的放大率来对从上述检体接收到的超声波信号进行放大,其中,上述存储部存储与接收深度相应的放大率信息;
B模式图像数据生成步骤,由B模式图像数据生成部生成将通过上述信号放大步骤放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据;
放大校正步骤,由放大校正部对通过上述信号放大步骤放大后的上述超声波信号进行放大校正以使不管接收深度如何放大率均固定;
频率分析步骤,由频率分析部对通过上述放大校正步骤放大校正后的上述超声波信号的频率进行分析来算出频谱;
特征量提取步骤,由特征量提取部对上述频率分析部算出的频谱进行衰减校正处理和近似处理来提取上述检体的特征量,其中,上述衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献;以及
特征量图像数据生成步骤,由特征量图像数据生成部生成特征量图像数据,该特征量图像数据用于显示与通过上述特征量提取步骤提取出的特征量对应的视觉信息。
14.一种超声波观测装置的动作程序,该超声波观测装置对检体发送超声波并且接收由上述检体反射的超声波,该动作程序的特征在于,使该超声波观测装置执行以下步骤:
信号放大步骤,由信号放大部按每个接收深度以从存储部读出的放大率来对从上述检体接收到的超声波信号进行放大,其中,上述存储部存储与接收深度相应的放大率信息;
B模式图像数据生成步骤,由B模式图像数据生成部生成将通过上述信号放大步骤放大后的超声波信号的振幅变换为亮度来显示的B模式图像数据;
放大校正步骤,由放大校正部对通过上述信号放大步骤放大后的上述超声波信号进行放大校正以使不管接收深度如何放大率均固定;
频率分析步骤,由频率分析部对通过上述放大校正步骤放大校正后的上述超声波信号的频率进行分析来算出频谱;
特征量提取步骤,由特征量提取部对上述频率分析部算出的频谱进行衰减校正处理和近似处理来提取上述检体的特征量,其中,上述衰减校正处理消减在超声波传播时与该超声波的接收深度和频率相应地产生的衰减的贡献;以及
特征量图像数据生成步骤,由特征量图像数据生成部生成特征量图像数据,该特征量图像数据用于显示与通过上述特征量提取步骤提取出的特征量对应的视觉信息。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010253290 | 2010-11-11 | ||
JP2010-253290 | 2010-11-11 | ||
PCT/JP2011/076027 WO2012063929A1 (ja) | 2010-11-11 | 2011-11-11 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103153195A true CN103153195A (zh) | 2013-06-12 |
CN103153195B CN103153195B (zh) | 2016-08-03 |
Family
ID=46051066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180048796.2A Active CN103153195B (zh) | 2010-11-11 | 2011-11-11 | 超声波观测装置、超声波观测装置的动作方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9028414B2 (zh) |
EP (1) | EP2599441B1 (zh) |
JP (1) | JP5307939B2 (zh) |
CN (1) | CN103153195B (zh) |
WO (1) | WO2012063929A1 (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104411250A (zh) * | 2013-01-23 | 2015-03-11 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN105828726A (zh) * | 2014-06-11 | 2016-08-03 | 奥林巴斯株式会社 | 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序 |
CN105939672A (zh) * | 2014-12-22 | 2016-09-14 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN106068099A (zh) * | 2014-12-22 | 2016-11-02 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN106659472A (zh) * | 2014-07-11 | 2017-05-10 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN106659478A (zh) * | 2015-05-13 | 2017-05-10 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN107427284A (zh) * | 2015-10-23 | 2017-12-01 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN107530057A (zh) * | 2015-05-13 | 2018-01-02 | 奥林巴斯株式会社 | 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序 |
CN108366782A (zh) * | 2015-12-08 | 2018-08-03 | 奥林巴斯株式会社 | 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序 |
CN108472018A (zh) * | 2015-12-24 | 2018-08-31 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10226234B2 (en) | 2011-12-01 | 2019-03-12 | Maui Imaging, Inc. | Motion detection using ping-based and multiple aperture doppler ultrasound |
US9668714B2 (en) * | 2010-04-14 | 2017-06-06 | Maui Imaging, Inc. | Systems and methods for improving ultrasound image quality by applying weighting factors |
EP2719337A4 (en) * | 2012-05-30 | 2015-04-08 | Olympus Medical Systems Corp | ULTRASONIC OBSERVATION DEVICE, OPERATING METHOD FOR A ULTRASONIC OBSERVATION DEVICE AND OPERATING PROGRAM FOR AN ULTRASONIC OBSERVATION DEVICE |
CA2896505A1 (en) * | 2013-01-08 | 2014-07-17 | Volcano Corporation | Method for focused acoustic computed tomography (fact) |
KR102218308B1 (ko) | 2013-05-21 | 2021-02-23 | 삼성전자주식회사 | 초음파 영상 처리 장치 및 방법 |
WO2014192954A1 (ja) * | 2013-05-29 | 2014-12-04 | オリンパスメディカルシステムズ株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP5788624B2 (ja) | 2013-07-18 | 2015-10-07 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
EP3078330A4 (en) * | 2013-12-05 | 2017-10-18 | Olympus Corporation | Ultrasonic observation device, ultrasonic observation device operation method, and ultrasonic observation device operation program |
WO2016038944A1 (ja) * | 2014-09-10 | 2016-03-17 | 富士フイルム株式会社 | 音響波画像生成装置およびその制御方法 |
CN106572843B (zh) * | 2014-12-25 | 2020-03-10 | 奥林巴斯株式会社 | 超声波观测装置以及超声波观测装置的工作方法 |
JP5953457B1 (ja) * | 2015-03-23 | 2016-07-20 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
EP3275376A4 (en) * | 2015-03-23 | 2019-01-16 | Olympus Corporation | ULTRASONIC OBSERVATION DEVICE, ULTRASONIC OBSERVATION DEVICE OPERATING PROCEDURE AND ULTRASONIC OBSERVATION DEVICE OPERATING PROGRAM |
JP6138402B2 (ja) * | 2015-03-31 | 2017-05-31 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP6010274B1 (ja) * | 2015-05-13 | 2016-10-19 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
KR20200093743A (ko) | 2019-01-28 | 2020-08-06 | 삼성전자주식회사 | 정보 저장 패턴을 포함하는 반도체 소자 및 그 제조 방법 |
WO2020157870A1 (ja) * | 2019-01-30 | 2020-08-06 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4092867A (en) * | 1977-02-10 | 1978-06-06 | Terrance Matzuk | Ultrasonic scanning apparatus |
US20040243001A1 (en) * | 2003-04-22 | 2004-12-02 | Zagzebski James A. | Parametric ultrasound imaging using angular compounding |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4016750B1 (en) * | 1975-11-06 | 1994-04-05 | Stanford Research Inst | Ultrasonic imaging method and apparatus |
FR2554238B1 (fr) | 1983-10-28 | 1986-02-28 | Labo Electronique Physique | Appareil d'exploration de milieux par echographie ultrasonore |
JP3779410B2 (ja) | 1997-02-03 | 2006-05-31 | 株式会社東芝 | 超音波照射装置 |
JP4336516B2 (ja) * | 2003-04-10 | 2009-09-30 | 株式会社リコー | 特徴量次元圧縮装置、マッチング装置、プログラム及び記憶媒体 |
US7601122B2 (en) | 2003-04-22 | 2009-10-13 | Wisconsin Alumni Research Foundation | Ultrasonic elastography with angular compounding |
JP4431354B2 (ja) * | 2003-10-06 | 2010-03-10 | アロカ株式会社 | 超音波診断装置 |
JP5203605B2 (ja) | 2004-06-18 | 2013-06-05 | 株式会社日立メディコ | 超音波診断装置 |
DE602005027617D1 (de) * | 2004-09-20 | 2011-06-01 | Resonant Medical Inc | Radiotherapiebehandlungsüberwachung unter einsatz von ultraschall |
JP2009082624A (ja) * | 2007-10-02 | 2009-04-23 | Olympus Medical Systems Corp | 超音波観測装置及びこの超音波観測装置を用いた超音波診断装置 |
US20100130860A1 (en) * | 2008-11-21 | 2010-05-27 | Kabushiki Kaisha Toshiba | Medical image-processing device, medical image-processing method, medical image-processing system, and medical image-acquiring device |
-
2011
- 2011-11-11 WO PCT/JP2011/076027 patent/WO2012063929A1/ja active Application Filing
- 2011-11-11 EP EP11840305.4A patent/EP2599441B1/en active Active
- 2011-11-11 CN CN201180048796.2A patent/CN103153195B/zh active Active
- 2011-11-11 JP JP2012533410A patent/JP5307939B2/ja active Active
-
2012
- 2012-07-30 US US13/561,443 patent/US9028414B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4092867A (en) * | 1977-02-10 | 1978-06-06 | Terrance Matzuk | Ultrasonic scanning apparatus |
US20040243001A1 (en) * | 2003-04-22 | 2004-12-02 | Zagzebski James A. | Parametric ultrasound imaging using angular compounding |
Non-Patent Citations (1)
Title |
---|
AMIN KATOUZIAN ET AL: "challenges in atherosclerotic plaque characterization with intravascular ultrasound(IVUS):from data collection to classification", 《IEEE TRANSACTIONS ON INFORMATION TECHNOLOGY IN BIOMEDICINE》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104411250A (zh) * | 2013-01-23 | 2015-03-11 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN105828726B (zh) * | 2014-06-11 | 2019-06-18 | 奥林巴斯株式会社 | 医用诊断装置以及医用诊断装置的工作方法 |
CN105828726A (zh) * | 2014-06-11 | 2016-08-03 | 奥林巴斯株式会社 | 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序 |
CN106659472A (zh) * | 2014-07-11 | 2017-05-10 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN106659472B (zh) * | 2014-07-11 | 2020-05-15 | 奥林巴斯株式会社 | 超声波观测装置以及超声波观测装置的工作方法 |
CN105939672A (zh) * | 2014-12-22 | 2016-09-14 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN106068099A (zh) * | 2014-12-22 | 2016-11-02 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN106068099B (zh) * | 2014-12-22 | 2019-08-02 | 奥林巴斯株式会社 | 超声波观测装置以及超声波观测装置的工作方法 |
CN106659478B (zh) * | 2015-05-13 | 2020-01-10 | 奥林巴斯株式会社 | 超声波观测装置以及超声波观测装置的工作方法 |
CN107530057A (zh) * | 2015-05-13 | 2018-01-02 | 奥林巴斯株式会社 | 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序 |
CN106659478A (zh) * | 2015-05-13 | 2017-05-10 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN107427284A (zh) * | 2015-10-23 | 2017-12-01 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN107427284B (zh) * | 2015-10-23 | 2021-04-13 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法及存储介质 |
CN108366782A (zh) * | 2015-12-08 | 2018-08-03 | 奥林巴斯株式会社 | 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序 |
CN108472018A (zh) * | 2015-12-24 | 2018-08-31 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN108472018B (zh) * | 2015-12-24 | 2021-01-05 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
Also Published As
Publication number | Publication date |
---|---|
WO2012063929A1 (ja) | 2012-05-18 |
JPWO2012063929A1 (ja) | 2014-05-12 |
JP5307939B2 (ja) | 2013-10-02 |
US20130030296A1 (en) | 2013-01-31 |
EP2599441A1 (en) | 2013-06-05 |
US9028414B2 (en) | 2015-05-12 |
CN103153195B (zh) | 2016-08-03 |
EP2599441B1 (en) | 2019-04-17 |
EP2599441A4 (en) | 2018-01-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103153195A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN102802536B (zh) | 超声波诊断装置、超声波诊断装置的动作方法以及超声波诊断装置的动作程序 | |
CN103200876B (zh) | 超声波观测装置、超声波观测装置的动作方法 | |
CN104125804B (zh) | 超声波观测装置与超声波观测装置的动作方法 | |
CN102958447B (zh) | 超声波观测装置、超声波观测装置的动作方法 | |
CN102834059B (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
JP5897227B1 (ja) | 医用診断装置、医用診断装置の作動方法および医用診断装置の作動プログラム | |
CN102858251B (zh) | 超声波观测装置及超声波观测方法 | |
CN102905624A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN103717146B (zh) | 超声波观测装置、超声波观测装置的动作方法 | |
CN104582584B (zh) | 超声波观测装置以及超声波观测装置的动作方法 | |
CN104411250A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
CN105392427B (zh) | 超声波观测装置以及超声波观测装置的动作方法 | |
WO2012063930A1 (ja) | 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム | |
CN103945773A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20151111 Address after: Tokyo, Japan, Japan Applicant after: Olympus Corporation Address before: Tokyo, Japan, Japan Applicant before: Olympus Medical Systems Corp. |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |