CN105916448B - 超声波观测装置以及超声波观测装置的工作方法 - Google Patents
超声波观测装置以及超声波观测装置的工作方法 Download PDFInfo
- Publication number
- CN105916448B CN105916448B CN201580004506.2A CN201580004506A CN105916448B CN 105916448 B CN105916448 B CN 105916448B CN 201580004506 A CN201580004506 A CN 201580004506A CN 105916448 B CN105916448 B CN 105916448B
- Authority
- CN
- China
- Prior art keywords
- attenuation rate
- characteristic quantity
- frequency
- correction
- attenuation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/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/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- 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/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- 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
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
-
- 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
Abstract
本发明提供超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序。本发明所涉及的超声波观测装置具备:特征量计算部,其分别计算通过对超声波信号的频率进行分析而得到的多个频谱的特征量;衰减率设定部,其针对将超声波图像分割为多个区域得到的各分割区域,使用赋予互不相同的衰减特性的、多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除超声波的影响的衰减校正,由此针对每个衰减率候选值计算的各频谱的预备校正特征量,设定多个衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及特征量校正部,其使用超声波振子的表面与采样点之间存在的分割区域的最佳的衰减率,来计算采样点处的每单位频率的累积衰减率,使用该累积衰减率来进行特征量的衰减校正,由此计算校正特征量。
Description
技术领域
本发明涉及一种使用超声波对观测对象的组织进行观测的超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序。
背景技术
以往,已知如下技术:在使用超声波对观测对象的组织进行观察的超声波观测装置中,为了对接收信号补偿生物体组织的频率依赖衰减而进行校正(例如,参照专利文献1)。在该技术中,使用对来自被检体的反射波依次进行与接收点的深度相应的动态校正处理和脉冲压缩处理得到的接收信号来形成超声波图像。
专利文献1:日本特开2010-246640号公报
发明内容
发明要解决的问题
然而,在上述的专利文献1所记载的技术中,在对远离发送和接收超声波的探头表面的关心区域的接收信号进行校正时,没有考虑超声波到达该关心区域之前的衰减率的变化。因此,只能应用于观测对象的衰减率均一的情况,从而难以高精度地辨别衰减率不均的观测对象的组织性状。
本发明是鉴于上述情形而完成的,其目的在于提供一种能够高精度地辨别衰减率不均的观测对象的组织性状的超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序。
用于解决问题的方案
为了解决上述的课题并达成目的,本发明所涉及的超声波观测装置基于由超声波探头获取到的超声波信号来生成超声波图像,该超声波探头具备向观测对象发送超声波并接收由该观测对象反射的超声波的超声波振子,该超声波观测装置的特征在于,具备:频率分析部,其通过对所述超声波信号的频率进行分析,来计算与所述超声波信号的接收深度及接收方向相应的多个频谱;特征量计算部,其分别计算所述多个频谱的特征量;衰减率设定部,其针对将所述超声波图像分割为多个区域得到的各分割区域,使用在所述超声波在所述观测对象中传播时赋予互不相同的衰减特性的、多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除所述超声波的影响的衰减校正,由此针对每个所述衰减率候选值计算各所述频谱的预备校正特征量,基于该计算的结果来设定多个所述衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及特征量校正部,其使用由所述衰减率设定部针对每个所述分割区域设定的最佳的衰减率中的、所述超声波振子的表面与采样点之间存在的所述分割区域的最佳的衰减率,来计算所述采样点处的每单位频率的累积衰减率,使用该累积衰减率进行所述特征量的衰减校正,由此计算校正特征量。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,所述特征量校正部将利用所述超声波振子的表面位置与所述采样点之间存在的各个分割区域的在深度方向上的往返距离对该各分割区域的所述最佳的衰减率进行加权得到的结果累加,由此计算所述采样点处的所述累积衰减率。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,关于在深度方向上相邻的两个所述分割区域,离所述超声波振子较远一方的所述分割区域的所述深度方向上的长度为离所述超声波振子较近一方的所述分割区域的所述深度方向上的长度以上。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,所述衰减率设定部针对每个所述衰减率候选值计算所述预备校正特征量的统计性偏差,将该统计性偏差最小的衰减率候选值设为所述最佳的衰减率。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,还具备特征量图像数据生成部,该特征量图像数据生成部生成将关于所述校正特征量的信息与所述超声波图像一并显示的特征量图像数据。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,所述特征量计算部通过进行利用n次式对各所述频谱进行近似的处理来计算所述特征量,其中,n为正整数。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,所述特征量计算部利用一次式对所述频谱中的规定的频带进行近似,计算所述一次式的截距、斜率、以及所述频带的中间频率下的所述一次式的值即频带中心对应强度中的一个或多个作为所述特征量,且所述特征量包含所述斜率和所述频带中心对应强度中的任一个,所述衰减率设定部基于所述斜率和所述频带中心对应强度中的任一个来设定所述最佳的衰减率。
本发明所涉及的超声波观测装置的特征在于,在上述发明中,所述衰减率设定部在将所述斜率设为所述特征量的情况下,基于所述斜率来设定所述最佳的衰减率,在将所述频带中心对应强度设为所述特征量的情况下,基于所述频带中心对应强度来设定所述最佳的衰减率。
本发明所涉及的超声波观测装置的工作方法是如下的超声波观测装置的工作方法,该超声波观测装置基于由超声波探头获取到的超声波信号来生成超声波图像,该超声波探头具备向观测对象发送超声波并接收由该观测对象反射的超声波的超声波振子,该超声波观测装置的工作方法的特征在于,包括:频率分析步骤,频率分析部通过对所述超声波信号的频率进行分析,来计算与所述超声波信号的接收深度及接收方向相应的多个频谱;特征量计算步骤,特征量计算部分别计算所述多个频谱的特征量;衰减率设定步骤,衰减率设定部针对将所述超声波图像分割为多个区域得到的各分割区域,使用在所述超声波在所述观测对象中传播时赋予互不相同的衰减特性的、多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除所述超声波的影响的衰减校正,由此针对每个所述衰减率候选值计算各所述频谱的预备校正特征量,基于该计算的结果来设定多个所述衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及特征量校正步骤,特征量校正部使用针对每个所述分割区域设定的最佳的衰减率中的、所述超声波振子的表面与采样点之间存在的所述分割区域的最佳的衰减率,来计算所述采样点处的每单位频率的累积衰减率,使用该累积衰减率来进行所述特征量的衰减校正,由此计算校正特征量。
本发明所涉及的超声波观测装置的工作程序的特征在于,使基于由具备向观测对象发送超声波并接收由该观测对象反射的超声波的超声波振子的超声波探头获取到的超声波信号来生成超声波图像的超声波观测装置执行以下步骤:频率分析步骤,频率分析部通过对所述超声波信号的频率进行分析,来计算与所述超声波信号的接收深度及接收方向相应的多个频谱;特征量计算步骤,特征量计算部分别计算所述多个频谱的特征量;衰减率设定步骤,衰减率设定部针对将所述超声波图像分割为多个区域得到的各分割区域,使用在所述超声波在所述观测对象中传播时赋予互不相同的衰减特性的、每单位长度和每单位频率的多个衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除所述超声波的影响的衰减校正,由此针对每个所述衰减率候选值计算各所述频谱的预备校正特征量,基于该计算的结果来设定所述多个衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及特征量校正步骤,特征量校正部使用针对每个所述分割区域设定的最佳的衰减率中的、所述超声波振子的表面与采样点之间存在的所述分割区域的最佳的衰减率,来计算所述采样点处的每单位频率的累积衰减率,使用该累积衰减率进行所述特征量的衰减校正,由此计算校正特征量。
发明的效果
根据本发明,在将超声波图像分割为多个区域得到的各分割区域,计算与赋予互不相同的衰减特性的、每单位长度和每单位频率的多个衰减率候选值相应的频谱的预备校正特征量,基于该计算结果来从多个衰减率候选值中对每个分割区域设定对于观测对象最佳的衰减率,使用超声波振子的表面与采样点之间存在的分割区域的最佳的衰减率,来计算采样点处的每单位频率的累积衰减率,使用该累积衰减率进行特征量的衰减校正,由此计算校正特征量,因此即使是衰减率不均一的观测对象,也能够计算考虑到不均一性的校正特征量。因而,根据本发明,能够高精度地辨别衰减率不均一的观测对象的组织性状。
附图说明
图1是表示具备本发明的一个实施方式所涉及的超声波观测装置的超声波诊断系统的功能结构的框图。
图2是表示本发明的一个实施方式所涉及的超声波观测装置的信号放大部进行的放大处理中的接收深度与放大率之间的关系的图。
图3是表示本发明的一个实施方式所涉及的超声波观测装置的放大校正部进行的放大校正处理中的接收深度与放大率之间的关系的图。
图4是示意性地表示超声波信号的一个声线中的数据排列的图。
图5是表示本发明的一个实施方式所涉及的超声波观测装置的频率分析部计算出的频谱的例子的图。
图6是示意性地表示超声波图像中的分割区域的设定例的图。
图7是表示由具有本发明的一个实施方式所涉及的超声波观测装置的衰减率设定部进行校正得到的预备校正特征量来作为参数的直线的图。
图8是示意性地表示基于两个不同的衰减率候选值对相同的观测对象分别进行衰减校正得到的预备校正特征量的分布例的图。
图9是表示本发明的一个实施方式所涉及的超声波观测装置进行的处理的概要的流程图。
图10是表示本发明的一个实施方式所涉及的超声波观测装置的频率分析部执行的处理的概要的流程图。
图11是表示本发明的一个实施方式所涉及的超声波观测装置的衰减率设定部进行的处理的概要的图。
图12是表示超声波图像中的分割区域的其它设定例的图。
图13是表示本发明的一个实施方式所涉及的超声波观测装置的衰减率设定部进行的最佳的衰减率的其它设定方法的概要的图。
具体实施方式
下面,参照附图来说明用于实施本发明的方式(以下称为“实施方式”)。
图1是表示具备本发明的一个实施方式所涉及的超声波观测装置的超声波诊断系统的功能结构的框图。该图所示的超声波诊断系统1具备:超声波内窥镜2,其向作为观测对象的被检体发送超声波,并接收由该被检体反射的超声波;超声波观测装置3,其基于由超声波内窥镜2获取到的超声波信号来生成超声波图像;以及显示装置4,其显示由超声波观测装置3生成的超声波图像。
超声波内窥镜2在其前端部具有超声波振子21,该超声波振子21将从超声波观测装置3接收到的电脉冲信号转换为超声波脉冲(声脉冲)来向被检体照射,并且将由被检体反射的超声波回波转换为以电压变化表现的电回波信号来输出。超声波振子21可以是凸起型振子、线性振子以及径向型振子中的任一个。超声波内窥镜2既可以是使超声波振子21进行机械式扫描的结构,也可以是如下结构:将多个元件阵列状地设置成超声波振子21,对与发送和接收有关的元件以电子方式进行切换、或者使各元件的发送和接收延迟,由此使超声波振子21进行电子式扫描。下面,在本实施方式中,为了便于说明,设超声波振子21为线性振子。
超声波内窥镜2通常具有摄像光学系统和摄像元件,能够被插入到被检体的消化管(食道、胃、十二指肠、大肠)或呼吸器官(气管、支气管)来对消化管、呼吸器官、其周围脏器(胰脏、胆囊、胆管、胆道、淋巴结、纵隔脏器、血管等)进行拍摄。另外,超声波内窥镜2具有在摄像时引导向被检体照射的照明光的光导件。该光导件的前端部到达超声波内窥镜2的向被检体插入的插入部的前端,另一方面,该光导件的基端部与产生照明光的光源装置连接。
超声波观测装置3与超声波内窥镜2电连接,该超声波观测装置3具备:发送和接收部31,其基于规定的波形和发送定时来向超声波振子21发送由高电压脉冲构成的发送信号(脉冲信号),并且从超声波振子21接收作为电接收信号的回波信号来生成数字的高频(RF:Radio Frequency(无线电频率))信号的数据(以下称为RF数据)并输出该数据;信号处理部32,其基于从发送和接收部31接收到的RF数据来生成数字的B模式用接收数据;运算部33,其对从发送和接收部31接收到的RF数据实施规定的运算;图像处理部34,其生成各种图像数据;输入部35,其使用键盘、鼠标、触摸面板等用户接口来实现,接受各种信息的输入;控制部36,其对超声波诊断系统1整体进行控制;以及存储部37,其存储超声波观测装置3进行动作所需要的各种信息。
发送和接收部31具有将回波信号放大的信号放大部311。信号放大部311进行STC(Sensitivity Time Control:灵敏度时间控制)校正,该STC校正为回波信号的接收深度越大则以越高的放大率进行放大的校正。图2是表示信号放大部311进行的放大处理中的接收深度与放大率之间的关系的图。图2所示的接收深度z是基于从超声波的接收开始时刻起的经过时间计算出的量。如图2所示,在接收深度z小于阈值zth的情况下,放大率β(dB)随着接收深度z的增加而从β0向βth(>β0)线性地增加。另外,在接收深度z为阈值zth以上的情况下,放大率β(dB)取固定值βth。阈值zth的值为从观测对象接收的超声波信号几乎全部衰减而噪声处于支配地位那样的值。更一般地,只要在接收深度z小于阈值zth的情况下放大率β随着接收深度z的增加而单调增加即可。图2所示的关系被预先存储于存储部37。
发送和接收部31在对通过信号放大部311而被放大了的回波信号实施滤波等处理之后,通过进行A/D转换来生成时域的RF数据并向信号处理部32和运算部33输出该RF数据。在超声波探头2具有使将多个元件阵列状地设置而成的超声波振子21进行电子式扫描的结构的情况下,发送和接收部31具有与多个元件对应的光束合成用的多通道电路。
发送和接收部31发送的脉冲信号的频带设为大致覆盖超声波振子21中的脉冲信号向超声波脉冲进行电声转换的线性响应频带的宽频带即可。另外,信号放大部311中的回波信号的各种处理频带设为大致覆盖由超声波振子21进行的超声波回波向回波信号的声电转换的线性响应频带的宽频带即可。由此,在执行后述的频谱的近似处理时,能够进行高精度的近似。
发送和接收部31还具有如下功能:对超声波内窥镜2发送由控制部36输出的各种控制信号,并且从超声波内窥镜2接收包含识别用的ID的各种信息并向控制部36发送上述各种信息。
信号处理部32对RF数据实施带通滤波、包络线检波、对数转换等公知的处理,来生成数字的B模式用接收数据。在对数转换中,取对RF数据除以基准电压得到的量的常用对数,以分贝值表示。信号处理部32向图像处理部34输出所生成的B模式用接收数据。信号处理部32使用CPU(Central Proccesing Unit:中央处理器)、各种运算电路等来实现。
运算部33具有:放大校正部331,其对由发送和接收部31输出的RF数据以使放大率不依赖于收深度而为固定的方式进行放大校正;频率分析部332,其对进行了放大校正的RF数据实施高速傅立叶变换(FFT:Fast Fourier Transform)来进行频率分析,由此计算与超声波信号的接收深度及接收方向相应的多个频谱;特征量计算部333,其计算各频谱的特征量;衰减率设定部334,其针对将超声波图像分割为多个区域时的各分割区域,使用在超声波在观测对象中传播时赋予互不相同的衰减特性的多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除超声波的影响的衰减校正,由此针对每个衰减率候选值计算各频谱的预备校正特征量,基于该计算的结果来设定多个衰减率候选值中的对于观测对象而言最佳的衰减率;以及特征量校正部335,其使用由衰减率设定部334针对每个分割区域设定的最佳的衰减率中的、超声波振子的表面与采样点之间存在的分割区域的最佳的衰减率,来计算采样点处的每单位频率的累积衰减率,使用该累积衰减率来进行特征量的衰减校正,由此计算校正特征量。运算部33使用CPU、各种运算电路等来实现。
图3是表示放大校正部331进行的放大校正处理中的接收深度与放大率之间的关系的图。如图3所示,关于放大校正部331进行的放大处理中的放大率β(dB),在接收深度z为零时取最大值βth-β0,在接收深度z从零起至达到阈值zth为止呈线性地减少,在接收深度z为阈值zth以上时为零。此外,图3所示的关系被预先存储于存储部37。放大校正部331基于图3所示的关系对数字RF信号进行放大校正,由此能够抵消信号放大部311中的STC校正的影响,从而输出放大率βth固定的信号。此外,放大校正部331进行的放大校正处理中的接收深度z与放大率β之间的关系根据信号放大部311进行的放大处理中的接收深度与放大率之间的关系不同而不同,这是不言而喻的。
说明进行这样的放大校正的理由。STC校正为如下的校正处理:使模拟信号波形的振幅在整个频带中均一地且以相对于深度而单调增加的放大率放大,由此从模拟信号波形的振幅排除衰减的影响。因此,在生成将回波信号的振幅转换为亮度来进行显示的B模式图像的情况下、且扫描一样的组织的情况下,通过进行STC校正来使亮度值不依赖于深度而为固定。即,能够获得从B模式图像的亮度值排除了衰减的影响的效果。
另一方面,在如本实施方式那样利用对超声波的频谱进行计算并进行分析得到的结果的情况下,存在即使进行STC校正也未必能够准确地排除伴随超声波的传播而发生的衰减的影响这样的问题。原因在于,一般地,衰减量根据频率不同而不同(参照后述的式(1)),但是STC校正的放大率仅根据距离发生变化,不具有频率依赖性。
为了解决上述的问题,考虑在生成B模式图像时输出实施了STC校正的接收信号,另一方面,在生成基于频谱的图像时,进行与用于生成B模式图像的发送不同的新的发送,输出未实施STC校正的接收信号。可是,在该情况下,存在基于接收信号生成的图像数据的帧频下降的问题。
因此,在本实施方式中,利用放大校正部331进行放大率的校正,以维持所生成的图像数据的帧频并且针对为了用于B模式图像而被实施了STC校正的信号排除STC校正的影响。
频率分析部332以规定的时间间隔对由放大校正部331进行了放大校正的各声线的RF数据(线数据)进行采样,来生成采样数据。频率分析部332通过对采样数据群实施FFT处理来计算RF数据上的多个位置(数据位置)处的频谱。
图4是示意性地表示超声波信号的一个声线中的数据排列的图。在该图所示的声线SRk中,白色或黑色的长方形意味着一个采样点处的数据。另外,在声线SRk中,数据的位置越靠右侧,则该数据是在沿声线SRk测量的情况下的从超声波振子21起的深度越深的位置处的采样数据(参照图4的箭头)。声线SRk以与发送和接收部31所进行的A/D转换中的采样频率(例如50MHz)对应的时间间隔被离散化。在图4中,示出了将编号k的声线SRk的第8个数据位置设定为接收深度z的方向上的初始值Z(k) 0的情况,但是初始值的位置能够任意地设定。关于频率分析部332的计算结果,能够以复数形式获得,并保存于存储部37。
图4所示的数据群Fj(j=1、2、…、K)是成为FFT处理的对象的采样数据群。一般地,为了进行FFT处理,采样数据群需要具有2的乘方的数据数。在该意义上,采样数据群Fj(j=2、…、K-1)的数据数为16(=24),是正常的数据群,另一方面,采样数据群FK的数据数为12,因此是异常的数据群。在对异常的数据群进行FFT处理时,进行通过对不足的部分插入零数据来生成正常的采样数据群的处理。关于这一点,在说明频率分析部332的处理时详细记述(参照图10)。
图5是表示频率分析部332计算出的频谱的例子的图。在此所说的“频谱”意味着通过对采样数据群实施FFT处理而得到的“某一接收深度z的强度的频率分布”。另外,在此所说的“强度”例如是指回波信号的电压、回波信号的电力、超声波回波的声压、超声波回波的声能等参数、这些参数的振幅、时间积分值、其组合中的任一个。
在图5中,取频率f为横轴。另外,在图5中,取对强度I0除以基准强度Ic(常数)得到的量的常用对数(用分贝表示)I=10log10(I0/Ic)为纵轴。在图5中,接收深度z是固定的。关于图5所示的直线L10,在后面记述。此外,在本实施方式中,曲线和直线由离散的点的集合构成。
在图5所示的频谱C1上,在以后的运算中使用的频带的下限频率fL和上限频率fH是基于超声波振子21的频带、发送和接收部31发送的脉冲信号的频带等决定的参数,例如,fL=3MHz,fH=10MHz。下面,在图5中,将由下限频率fL和上限频率fH确定的频带称为“频带U”。
一般地,在观测对象为生物体组织的情况下,频谱根据超声波所扫描的生物体组织的性状不同而表示出不同的倾向。这是因为,频谱与使超声波散射的散射体的大小、数密度、声音阻抗等具有相关性。在此所说的“生物体组织的性状”例如是指恶性肿瘤(癌)、良性肿瘤、内分泌肿瘤、粘液性肿瘤、正常组织、囊肿、血管等。
特征量计算部333通过对规定频带中的频谱进行回归分析来用一次式对该频谱进行近似,由此计算使该近似的一次式具有特征的特征量。例如,在图5所示的频谱C1的情况下,特征量计算部333通过针对频带U进行回归分析而获得近似直线L10。下面,当以频率f的一次式I=a0f+b0来表示近似直线L10时,特征量计算部333计算斜率a0、截距b0以及频带U的中心频率fM=(fs+fe)/2处的强度I的值即频带中心对应强度(Mid-band fit,日语:ミッドバンドフィット)c0=a0fM+b0来作为与直线L10对应的特征量。此外,特征量计算部333也可以通过回归分析来用二次以上的多项式对频谱进行近似。
三个校正前特征量中的斜率a0与超声波的散射体的大小具有相关性,一般认为散射体越大则斜率具有越小的值。另外,截距b0与散射体的大小、声音阻抗的差、散射体的数密度(浓度)等具有相关性。具体地说,认为散射体越大则截距b0具有越大的值,声音阻抗的差越大则截距b0具有越大的值,散射体的数密度越大则截距b0具有越大的值。频带中心对应强度c0是基于斜率a0和截距b0导出的间接的参数,被赋予有效的频带内的中心处的频谱的强度。因此,认为频带中心对应强度c0除了与散射体的大小、声音阻抗的差、散射体的数密度具有相关性之外,还与B模式图像的亮度具有某种程度的相关性。
衰减率设定部334针对将超声波图像分割为多个区域时的各分割区域,使用赋予超声波的每单位长度和每单位频率的衰减量的衰减率来进行特征量的衰减校正。图6是示意性地表示超声波图像中的分割区域的设定例的图。在图6中,示出了将呈矩形状的超声波图像101分割为30个分割区域Pij(i=1~5,j=1~6)的情况。此外,图6的虚线表示声线。在图6所示的情况下,分割区域Pij全部具有相同的面积。因此,分割区域Pij的在相对于超声波振子21的表面位置102的深度的方向上的长度H(以下称为分割区域Pij的高度H)是固定的。与分割区域Pij有关的信息被保存于存储部37所具有的分割区域信息存储部371。此外,声线的数量只是一例,并不限定于图6所示的情况。另外,图6是超声波振子21为线性振子的情况下的超声波图像,因此呈矩形状,但是在其它类型的超声波振子21的情况下,超声波图像的外形形状是不同的,这是不言而喻的。
一般地,超声波的衰减量A(f,z)是超声波在往返于接收深度0与接收深度z之间的期间发生的衰减,被定义为往返前后的强度变化(用分贝表示的差)。根据经验可知,在一样的组织内,衰减量A(f,z)与频率成比例,用以下的式(1)表示。
A(f,z)=2αzf…(1)
在此,比例常数α是被称为衰减率的量,赋予每单位长度和每单位频率的超声波的衰减量。另外,z是超声波的接收深度,f是频率。在观测对象为生物体的情况下,衰减率α的具体的值根据生物体的部位来确定。衰减率α的单位例如是dB/cm/MHz。
衰减率设定部334设定多个衰减率候选值中的最佳的衰减率。此时,衰减率设定部334使用衰减率候选值α,按照以下所示的式(2)~(4)来对由特征量计算部333所计算出的特征量(斜率a0、截距b0、频带中心对应强度c0)进行衰减校正,由此计算预备校正特征量a、b、c。
a=a0+2αz…(2)
b=b0…(3)
c=c0+A(fM,z)=c0+2αzfM(=afM+b)…(4)
从式(2)、(4)也显而易见地,衰减率设定部334进行超声波的接收深度z越大则校正量越大的校正。另外,根据式(3),关于截距的校正是恒等变换。这是因为,截距是与频率0(Hz)对应的频率成分,不受衰减的影响。
图7是表示具有由衰减率设定部334进行校正得到的预备校正特征量a、b、c来作为参数的直线的图。直线L1的式子表示如下。
I=af+b=(a0+2αz)f+b0…(5)
从该式(5)也显而易见地,与衰减校正前的直线L10相比,直线L1的斜率较大(a>a0)且截距相同(b=b0)。
衰减率设定部334针对各分割区域,设定按每个衰减率候选值计算出的预备校正特征量的统计性偏差最小的衰减率候选值来作为最佳的衰减率。在本实施方式中,应用方差来作为表示统计性偏差的量。在该情况下,衰减率设定部334将方差最小的衰减率候选值设定为最佳的衰减率。上述的三个预备校正特征量a、b、c中的独立的预备校正特征量为两个。此外,预备校正特征量b不依赖于衰减率。因而,在针对预备校正特征量a、c来设定最佳的衰减率的情况下,衰减率设定部334只要计算预备校正特征量a和c中的某一方的方差即可。
但是,更优选的是,衰减率设定部334在使用预备校正特征量a来设定最佳的衰减率的情况下,应用预备校正特征量a的方差,在使用预备校正特征量c来设定最佳的衰减率的情况下,应用预备校正特征量c的方差。原因在于,赋予衰减量A(f,z)的式(1)只不过是理想的数式,在现实中以下的式(6)更适当。
A(f,z)=2αzf+2α1z…(6)
式(6)的右边第二项的α1是表示信号强度与超声波的接收深度z成比例地变化的大小的系数,α1是表示由于观测对象的组织不均、光束合成时的通道数的变更等而发生的信号强度的变化的系数。由于存在式(6)的右边第二项,因此在使用预备校正特征量c来设定最佳的衰减率的情况下,应用预备校正特征量c的方差能够更准确地对衰减进行校正(参照式(4))。另一方面,在使用与频率f成比例的系数即预备校正特征量a来设定最佳的衰减率的情况下,应用预备校正特征量a的方差能够排除式(6)的右边第二项的影响来更准确地对衰减进行校正。
在此,说明能够基于统计性偏差来设定最佳的衰减率的理由。认为在应用了对于观测对象而言最佳的衰减率的情况下,同观测对象与超声波振子21之间的距离无关地,特征量向观测对象所固有的值收敛,统计性偏差变小。另一方面,认为在将不适合于观测对象的衰减率候选值设为最佳的衰减率的情况下,衰减校正过量或不足,因此根据观测对象与超声波振子21之间的距离而在特征量中产生偏离,特征量的统计性偏差变大。因而,可以说统计性偏差最小的衰减率候选值是对于观测对象而言最佳的衰减率。
图8是示意性地表示基于两个不同的衰减率候选值分别对相同的观测对象进行衰减校正而得到的预备校正特征量的分布例的图。在图8中,将横轴设为预备校正特征量,将纵轴设为频度。关于图8所示的两条分布曲线N1、N2,频度的总和是相同的。在图8所示的情况下,与分布曲线N2相比,分布曲线N1的特征量的统计性偏差小(方差小),形成波形陡峭的形状。因而,衰减率设定部334在设定与这两条分布曲线N1、N2对应的两个衰减率候选值中的最佳的衰减率的情况下,将与分布曲线N1对应的衰减率候选值设定为最佳的衰减率。
特征量校正部335使用每个分割区域的最佳的衰减率来计算采样点处的每单位频率的累积衰减率(以下,也仅称为累积衰减率),使用该累积衰减率来进行特征量的衰减校正。任意的采样点处的累积衰减率是使用该采样点距超声波振子21的表面的距离以及超声波振子21的表面与采样点之间存在的分割区域的最佳的衰减率来计算的。在分割区域PIJ(1≤I≤imax、1≤J≤jmax;I、J为常数)中距靠近超声波振子21的一侧的边界的距离为h(0<h≤H)那样的采样点SIJ(h)处的累积衰减率γIJ(h)表示如下。
γIJ(h)=[∑j=1,2,…,J-1(2H·α(PIj))]+2h·α(PIJ)…(7)
在此,Σj=1,2,…,J-1意味着取j=1~J-1的和,α(PIj)表示分割区域PIj的最佳的衰减率。式(7)的右边第一项是对超声波的在每个分割区域的往返距离2H乘以最佳的衰减率α(PIj)得出的衰减率的和。另外,式(7)的右边第二项是对超声波到分割区域PIJ中的采样点SIJ(h)为止的往返距离2h乘以分割区域PIJ的最佳的衰减率α(PIJ)得出的衰减率。通过这样,特征量校正部335将从超声波振子21的表面开始的衰减率累积起来,由此计算累积衰减率γIJ(h)。在将最佳的衰减率的单位设为dB/cm/MHz的情况下,累积衰减率的单位为dB/MHz。
例如,在图6中例示了在式(7)中I=2、J=3(其中,imax=5、jmax=6)的情况。在该情况下,根据式(7),采样点S23(h)处的累积衰减率γ23(h)如下。
γ23(h)=2H·α(P21)+2H·α(P22)+2h·α(P23)…(8)
特征量校正部335使用累积衰减率γIJ(h),如以下那样对采样点SIJ(h)处的特征量进行衰减校正。
aIJ(h)=a0+2γIJ(h)…(9)
bIJ(h)=b0…(10)
cIJ(h)=c0+2fMγIJ(h)…(11)
图像处理部34具有:B模式图像数据生成部341,其生成将回波信号的振幅转换为亮度来进行显示的超声波图像即B模式图像数据;以及特征量图像数据生成部342,其生成显示与由特征量计算部333计算出的特征量有关的信息的特征量图像数据。
B模式图像数据生成部341对从信号处理部32接收到的B模式用接收数据进行增益处理、对比度处理等使用公知技术的信号处理,并且进行与根据显示装置4中的图像的显示范围确定的数据步长相应的数据的间除等,由此生成B模式图像数据。B模式图像是使采用RGB颜色系统来作为颜色空间的情况下的变量即R(红色)、G(绿色)、B(蓝色)的值一致后的灰度图像。
B模式图像数据生成部341对B模式用接收数据实施重新排列的坐标转换使得能够在空间上正确地表现扫描范围,之后实施B模式用接收数据间的插值处理,由此填充B模式用接收数据间的空隙,从而生成B模式图像数据。B模式图像数据生成部341将所生成的B模式图像数据输出到特征量图像数据生成部342。
特征量图像数据生成部342通过将对B模式图像数据中的图像的各像素叠加与由特征量计算部333计算出的特征量相关联的视觉信息,来生成特征量图像数据。特征量图像数据生成部342例如对与图4所示的一个振幅数据群Fj(j=1、2、…、K)的数据量对应的像素区域分配与根据该振幅数据群Fj计算出的频谱的特征量对应的视觉信息。特征量图像数据生成部342例如通过将作为视觉信息的色调与上述的斜率、截距、频带中心对应强度中的任一个对应起来,来生成特征量图像数据。此外,特征量图像数据生成部342也可以将色调与从斜率、截距、频带中心对应强度中选择的两个特征量中的一个对应起来,并且将亮暗与另一个对应起来,由此生成特征量图像数据。作为与特征量相关联的视觉信息,除了有色调、亮暗(明度)之外,例如还能够列举饱和度、亮度值、R(红色)、G(绿色)、B(蓝色)等构成规定的颜色系统的颜色空间的变量。
控制部36使用具有运算和控制功能的CPU(Central Proccesing Unit:中央处理单器)、各种运算电路等来实现。控制部36从存储部37读出存储部37所存储、保存的信息,执行与超声波观测装置3的工作方法相关联的各种运算处理,由此对超声波观测装置3进行综合控制。此外,也能够使用与信号处理部32及运算部33共用的CPU等来构成控制部36
存储部37具有:分割区域信息存储部371,其存储与分割区域有关的信息;频谱信息存储部372,其将由频率分析部332计算出的频谱的信息与接收深度及接收方向一同存储;特征量信息存储部373,其存储与由特征量计算部333计算出的特征量及由特征量校正部335进行校正得到的校正特征量有关的信息;以及衰减率信息存储部374,其存储与由衰减率设定部334针对每个分割区域设定的最佳的衰减率及由特征量校正部335计算出的每个采样点的累积衰减率有关的信息。
存储部37除了存储上述信息以外,例如还存储进行放大处理所需要的信息(图2所示的放大率与接收深度之间的关系)、进行放大校正处理所需要的信息(图3所示的放大率与接收深度之间的关系)、进行衰减校正处理所需要的信息(参照式(1))、进行频率分析处理所需要的窗函数(Hamming(汉明窗)、Hanning(汉宁窗)、Blackman(布莱克曼窗)等)的信息等。
另外,存储部37还存储包含用于执行超声波观测装置3的工作方法的工作程序的各种程序。工作程序也能够记录于硬盘、快闪存储器、CD-ROM、DVD-ROM、软盘等计算机可读取的记录介质来广泛地流通。此外,上述的各种程序也能够通过经由通信网络下载来获取。在此所说的通信网络例如通过已有的公共线路网、LAN(Local Area Network:局域网)、WAN(Wide Area Network:广域网)等来实现,不论有线、无线方式都可以。
具有以上的结构的存储部37使用预先安装有各种程序等的ROM(Read OnlyMemory:只读存储器)以及用于存储各处理的运算参数、数据等的RAM(Random AccessMemory:随机存取存储器)等来实现。
图9是表示具有以上结构的超声波观测装置3进行的处理的概要的流程图。具体地说,是表示超声波观测装置3从超声波内窥镜2接收回波信号以后的处理的概要的流程图。下面,参照图9来说明超声波观测装置3进行的处理。首先,超声波观测装置3从超声波内窥镜2接收作为利用超声波振子21得到的观测对象的测量结果的回波信号(步骤S1)。
从超声波振子21接收到回波信号的信号放大部311进行该回波信号的放大(步骤S2)。在此,信号放大部311例如基于图2所示的放大率与接收深度之间的关系来进行回波信号的放大(STC校正)。
接着,B模式图像数据生成部341使用由信号放大部311放大后的回波信号来生成B模式图像数据,并向显示装置4输出该B模式图像数据(步骤S3)。接收到B模式图像数据的显示装置4显示与该B模式图像数据对应的B模式图像。
放大校正部331对从发送和接收部31输出的RF数据以使放大率与接收深度无关地而为固定的方式进行放大校正(步骤S4)。在此,放大校正部331例如基于图3所示的放大率与接收深度之间的关系来进行放大校正。
之后,频率分析部332通过对放大校正后的各声线的RF数据进行基于FFT的频率分析来计算针对所有的采样数据群的频谱,并将该频谱保存到频谱信息存储部372(步骤S5)。图10是表示在步骤S5中频率分析部332进行的处理的概要的流程图。下面,参照图10所示的流程图来详细地说明频率分析处理。
首先,频率分析部332将用于识别分析对象的声线的计数值k设为k0(步骤S21)。
接着,频率分析部332对代表为了使用于FFT运算而生成的一系列的数据群(采样数据群)的数据位置(相当于接收深度)Z(k)的初始值Z(k) 0进行设定(步骤S22)。例如,在图4中,如上述那样示出了将声线SRk的第8个数据位置设定为初始值Z(k) 0的情况。
之后,频率分析部332获取采样数据群(步骤S23),对获取到的采样数据群作用存储部37所存储的窗函数(步骤S24)。通过像这样对采样数据群作用窗函数,能够避免采样数据群在边界处变得不连续,能够防止产生伪像。
接着,频率分析部332判定数据位置Z(k)的采样数据群是否为正常的数据群(步骤S25)。如参照图4时所说明的那样,采样数据群需要具有2的乘方的数据数。下面,将正常的采样数据群的数据数设为2n(n为正整数)。在本实施方式中,尽可能将数据位置Z(k)设定为Z(k)所属的采样数据群的中心。具体地说,由于采样数据群的数据数为2n,因此将Z(k)设定为接近该采样数据群的中心的第2n/2(=2n-1)个位置。在该情况下,采样数据群正常意味着在比数据位置Z(k)浅的部分侧存在2n-1-1(设为2n-1-1=N)个数据,在比数据位置Z(k)深的部分侧存在2n-1(设为2n-1=M)个数据。在图4所示的情况下,采样数据群Fj(j=1、2、…、K-1)是正常的。此外,在图4中例示了n=4(N=7、M=8)的情况。
在步骤S25中的判定的结果为数据位置Z(k)的采样数据群正常的情况下(步骤S25:“是(Yes)”),频率分析部332转移到后述的步骤S27。
在步骤S25中的判定的结果为数据位置Z(k)的采样数据群不正常的情况下(步骤S25:“否(No)”),频率分析部332通过对不足的部分插入零数据来生成正常的采样数据群(步骤S26)。关于在步骤S25中被判定为不正常的采样数据群(例如图4的采样数据群FK),在追加零数据之前作用窗函数。因此,即使在采样数据群中插入零数据也不产生数据的不连续。在步骤S26之后,频率分析部332转移到后述的步骤S27。
在步骤S27中,频率分析部332使用采样数据群来进行FFT运算,由此得到振幅的频率分布、即频谱(步骤S27)。
接着,频率分析部332使数据位置Z(k)以步长D变化(步骤S28)。设步长D预先存储于存储部37。在图4中,例示了D=15的情况。期望的是步长D与在B模式图像数据生成部341生成B模式图像数据时利用的数据步长一致,但是在想要削减频率分析部332中的运算量的情况下,也可以设定比数据步长大的值来作为步长D。
之后,频率分析部332判定数据位置Z(k)是否大于声线SRk中的最大值Z(k) max(步骤S29)。在数据位置Z(k)大于最大值Z(k) max的情况下(步骤S29:“是”),频率分析部332使计数值k增加1(步骤S30)。这意味着将处理移向旁边的声线。另一方面,在数据位置Z(k)为最大值Z(k) max以下的情况下(步骤S29:“否”),频率分析部332返回步骤S23。
在步骤S30之后,频率分析部332判定计数值k是否大于最大值kmax(步骤S31)。在计数值k大于kmax的情况下(步骤S31:“是”),频率分析部332结束一系列的频率分析处理。另一方面,在计数值k为kmax以下的情况下(步骤S31:“否”),频率分析部332返回步骤S22。该最大值kmax设为手术操作者等用户通过输入部35任意地指示输入的值、或存储部37中预先设定的值。
通过这样,频率分析部332对分析对象区域内的(kmax-k0+1)个声线分别进行多次的FFT运算。作为FFT运算的结果而得到的频谱与接收深度、接收方向一同被保存于频谱信息存储部372。
此外,在以上的说明中,设为频率分析部332对接收到超声波信号的全部区域进行频率分析处理,但是也能够设为输入部35构成为能够接受按特定的深度幅度和声线宽度划分出的部分区域的设定输入,仅在所设定的部分区域内进行频率分析处理。
继以上所说明的步骤S5的频率分析处理之后,特征量计算部333针对各分割区域中所包含的采样点计算频谱的特征量(步骤S6)。具体地说,特征量计算部333通过对规定的频带的频谱进行回归分析来用一次式I=a0f+b0进行近似,计算斜率a0、截距b0、频带中心对应强度c0来作为特征量。例如,图5所示的直线L10是特征量计算部333通过对频带U的频谱C1进行回归分析而得到近似的回归直线。
之后,衰减率设定部334将在进行后述的衰减校正时应用的衰减率候选值α的值设定为规定的初始值α0(步骤S7)。该初始值α0的值预先被存储于衰减率信息存储部374即可。
接着,衰减率设定部334将衰减率候选值设为α,来对由特征量计算部333对各频谱近似而得到的特征量进行衰减校正,由此计算预备校正特征量,并将计算出的预备校正特征量与衰减率候选值α一同保存于特征量信息存储部373(步骤S8)。图7所示的直线L1是衰减率设定部334进行衰减校正处理而得到的直线的例子。
在步骤S8中,衰减率设定部334是通过将使用超声波信号的声线的数据排列而得到的数据位置Z=(fsp/2vs)Dn代入到上述的式(2)、(4)中的接收深度z来计算的。在此,fsp是数据的采样频率,vs是声速,D是数据步长,n是从声线的第一个数据起到处理对象的振幅数据群的数据位置为止的数据步数。例如,如果将数据的采样频率fsp设为50MHz、将声速vs设为1530m/sec,采用图4所示的数据排列而将步长D设为15,则z=0.2295n(mm)。
衰减率设定部334计算由衰减率设定部334对各频谱进行衰减校正而得到的多个预备校正特征量中的作为代表的预备校正特征量的方差,并将该方差与衰减率候选值α以相对应的方式保存到特征量信息存储部373(步骤S9)。例如,在预备校正特征量为斜率a、频带中心对应强度c的情况下,如上述那样,衰减率设定部334计算预备校正特征量a和c中的任一方的方差。在该步骤S9中,优选的是,在使用预备校正特征量a来生成特征量图像的情况下,应用各分割区域的预备校正特征量a的方差,在使用预备校正特征量c来生成特征量图像的情况下,应用各分割区域的预备校正特征量c的方差。
之后,衰减率设定部334使衰减率候选值α的值增加Δα(步骤S10),对增加后的衰减率候选值α与规定的最大值αmax的大小进行比较(步骤S11)。在步骤S11中的比较的结果为衰减率候选值α大于最大值αmax的情况下(步骤S11:“是”),超声波观测装置3转移到步骤S12。另一方面,在步骤S11中的比较的结果为衰减率候选值α为最大值αmax以下的情况下(步骤S11:“否”),超声波观测装置3返回步骤S8。
在步骤S12中,衰减率设定部334针对每个分割区域,参照特征量信息存储部373所存储的每个衰减率候选值的预备校正特征量的方差,将该方差最小的衰减率候选值设定为分割区域的最佳的衰减率(步骤S12)。
图11是表示衰减率设定部334进行的处理的概要的图。是表示设α0=0(dB/cm/MHz)、αmax=1.0(dB/cm/MHz)、Δα=0.2(dB/cm/MHz)的情况下的衰减率候选值α与方差S(α)之间的关系的例子的图。在图11所示的情况下,在衰减率候选值α为0.2(dB/cm/MHz)时,方差取最小值S(α)min。因而,在图11所示的情况下,衰减率设定部334将α=0.2(dB/cm/MHz)设定为最佳的衰减率。
之后,特征量校正部335使用由衰减率设定部334针对各分割区域设定的最佳的衰减率,来计算采样点处的累积衰减率(步骤S13)。例如,图6所示的分割区域PIJ的采样点SIJ(h)处的累积衰减率γIJ(h)由式(7)来表示。
接着,特征量校正部335使用累积衰减率来对特征量进行衰减校正,由此计算校正特征量(步骤S14)。例如,关于图6所示的分割区域PIJ的采样点SIJ(h)处的斜率a0、截距b0、频带中心对应强度c0的校正特征量aIJ(h)、bIJ(h)、cIJ(h),分别使用式(9)~(11)来计算。
特征量图像数据生成部342通过对由B模式图像数据生成部341生成的B模式图像数据中的各像素叠加与在步骤S14中计算出的校正特征量相关联的视觉信息(例如色调),来生成特征量图像数据(步骤S15)。特征量图像数据生成部342向显示装置4发送所生成的特征量图像数据。接收到特征量图像数据的显示装置4显示与接收到的特征量图像数据对应的特征量图像。
在步骤S15之后,超声波观测装置3结束一系列的处理。此外,超声波观测装置3周期性地重复执行步骤S1~S15的处理。
根据以上说明的本发明的一个实施方式,针对将超声波图像分割为多个区域得到的各分割区域,计算与赋予互不相同的衰减特性的多个每单位长度和每单位频率的衰减率候选值相应的频谱的预备校正特征量,基于该计算结果针对每个分割区域设定多个衰减率候选值中的对于观测对象而言最佳的衰减率,使用超声波振子的表面与采样点之间存在的分割区域的最佳的衰减率来计算采样点处的每单位频率的累积衰减率,使用该累积衰减率进行特征量的衰减校正,由此计算校正特征量,因此即使是衰减率不均的观测对象,也能够计算考虑到不均性的校正特征量。因而,根据本发明,能够高精度地辨别衰减率不均的观测对象的组织性状。
另外,根据本实施方式,将利用超声波振子的表面位置与采样点之间存在的各分割区域的在深度方向上的往返距离对该各分割区域的最佳的衰减率进行加权得到的结果累加,由此计算采样点处的累积衰减率,因此能够适当地设定从超声波振子的表面到采样点为止的中途的衰减率。
目前为止说明了用于实施本发明的方式,但是本发明不应仅限定于上述的实施方式。例如,超声波图像的分割区域的设定方法并不限于图6所示的方法。图12是表示超声波图像中的分割区域的其它设定例的图。在图12中,示出了将超声波图像201分割为30个分割区域Qij(i=1~5,j=1~6)的情况。关于分割区域Qij,深度相同的区域、即j的值相同的区域全部具有相同的面积。另一方面,关于分割区域Qij,在深度不同的区域中,距超声波振子21的表面位置202的深度越大,则分割区域Qij的高度Hj(j=1~6)越大(H1<H2<…<H6)。一般地,深度越大则超声波的衰减越大,因此通过使离超声波振子21较远的分割区域的面积变大,能够提高远处的S/N比。
在该情况下,在分割区域QIJ(1≤I≤imax、1≤J≤jmax,I、J是常数)中距靠近超声波振子21的一侧的边界的距离为h(0<h≤HJ)那样的采样点SIJ(h)处的累积衰减率γIJ(h)表示如下。
γIJ(h)=[∑j=1,2,…,J-1(2Hj·α(QIj))]+2h·α(QIJ)…(12)
该式(12)式(7)相比,不同点在于右边第一项的分割区域QIj的深度方向上的长度Hj依赖于区域位置。
另外,衰减率设定部334设定最佳的衰减率的设定方法并不限定于上述的方法。图13是表示衰减率设定部334设定最佳的衰减率的其它设定方法的概要的图。在图13中,示出了设α0=0(dB/cm/MHz)、αmax=1.0(dB/cm/MHz)、Δα=0.2(dB/cm/MHz)的情况下的衰减率候选值α与方差S(α)之间的关系的例子,衰减率候选值α=0、0.2、0.4、0.6、0.8、1.0(单位均为dB/cm/MHz)的方差S(α)的值分别与图11相同。在该情况下,在衰减率设定部334设定最佳的衰减率之前,特征量计算部333通过进行回归分析来计算对衰减率候选值α的方差S(α)的值进行插值得到的曲线R。之后,衰减率设定部334针对该曲线R计算0(dB/cm/MHz)≤α≤1.0(dB/cm/MHz)中的最小值S’(α)min,将此时的衰减率候选值的值α’设定为最佳的衰减率。因而,在图13所示的情况下,最佳的衰减率α’为0(dB/cm/MHz)与0.2(dB/cm/MHz)之间的值。
另外,也可以是,衰减率设定部334针对超声波图像的全部帧分别计算与最佳的衰减率相当的最佳衰减率相当值,将包含最新的帧的最佳衰减率相当值在内的规定数量的最佳衰减率相当值的平均值、中间值或最频值设定为最佳的衰减率。在该情况下,与针对各帧设定最佳的衰减率的情况相比,最佳的衰减率的变化减少,从而能够使其值稳定。
另外,也可以是,衰减率设定部334以超声波图像的规定的帧为间隔设定最佳的衰减率。由此,能够大幅地削减计算量。在该情况下,在下一次设定最佳的衰减率之前的期间,使用最后设定的最佳的衰减率的值即可。
另外,也可以设为输入部35能够接受衰减率候选值的初始值α0的设定变更的输入的结构。
另外,作为赋予统计性偏差的量,例如也能够应用标准偏差、总体中的特征量的最大值与最小值之差、特征量的分布的半值宽度中的任一个。此外,还考虑应用方差的倒数来作为赋予统计性偏差的量的情况,在该情况下,值最大的衰减率候选值为最佳的衰减率,这是不言而喻的。
另外,衰减率设定部334也能够分别计算多个种类的预备校正特征量的统计性偏差,将统计性偏差最小的情况下的衰减率候选值设定为最佳的衰减率。
另外,衰减率设定部334也可以使用多个衰减率候选值来对频谱进行衰减校正,通过对衰减校正后的频谱进行回归分析来计算预备校正特征量。
另外,对超声波内窥镜以外的超声波探头也能够应用。作为超声波探头,例如也可以应用没有光学系统的细径的超声波微型探头。超声波微型探头通常在被插入到胆道、胆管、胰管、气管、支气管、尿道、尿管来观察其周围脏器(胰脏、肺、前列腺、膀胱、淋巴结等)时使用。另外,作为超声波探头,还可以应用从被检体的体表照射超声波的体外式超声波探头。体外式超声波探头通常在观察腹部脏器(肝脏、胆囊、膀胱)、乳房(特别是乳腺)、甲状腺时使用。
这样,本发明能够在不脱离权利要求书所记载的技术思想的范围内包含各种实施方式。
产业上的可利用性
如以上那样,本发明所涉及的超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序对于高精度地辨别衰减率不均的观测对象的组织性状是有用的。
附图标记说明
1:超声波诊断系统;2:超声波内窥镜;3:超声波观测装置;4:显示装置;21:超声波振子;31:发送和接收部;32:信号处理部;33:运算部;34:图像处理部;35:输入部;36:控制部;37:存储部;101、201:超声波图像;102、202:表面位置;311:信号放大部;331:放大校正部;332:频率分析部;333:特征量计算部;334:衰减率设定部;335:特征量校正部;341:B模式图像数据生成部;342:特征量图像数据生成部;371:分割区域信息存储部;372:频谱信息存储部;373:特征量信息存储部;374:衰减率信息存储部。
Claims (9)
1.一种超声波观测装置,基于由超声波探头获取到的超声波信号来生成超声波图像,该超声波探头具备向观测对象发送超声波并接收由该观测对象反射的超声波的超声波振子,该超声波观测装置的特征在于,具备:
频率分析部,其通过对所述超声波信号的频率进行分析,来计算与所述超声波信号的接收深度及接收方向相应的多个频谱;
特征量计算部,其分别计算所述多个频谱的特征量;
衰减率设定部,其针对将所述超声波图像分割为多个区域得到的各分割区域,使用在所述超声波在所述观测对象中传播时赋予互不相同的衰减特性的、多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除所述超声波的影响的衰减校正,由此针对每个所述衰减率候选值计算各所述频谱的预备校正特征量,基于该计算的结果来设定多个所述衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及
特征量校正部,其使用由所述衰减率设定部针对每个所述分割区域设定的最佳的衰减率中的、所述超声波振子的表面与采样点之间存在的所述分割区域的最佳的衰减率,来计算所述采样点处的每单位频率的累积衰减率,使用该累积衰减率来进行所述特征量的衰减校正,由此计算校正特征量,
其中,所述衰减率设定部针对每个所述衰减率候选值计算所述预备校正特征量的统计性偏差,将该统计性偏差最小的衰减率候选值设为所述最佳的衰减率。
2.根据权利要求1所述的超声波观测装置,其特征在于,
所述特征量校正部将利用所述超声波振子的表面位置与所述采样点之间存在的各分割区域的在深度方向上的往返距离对该各分割区域的所述最佳的衰减率进行加权得到的结果累加,由此计算所述采样点处的所述累积衰减率。
3.根据权利要求1所述的超声波观测装置,其特征在于,
关于在深度方向上相邻的两个所述分割区域,离所述超声波振子较远的一方的所述分割区域的所述深度方向上的长度为离所述超声波振子较近的一方的所述分割区域的所述深度方向上的长度以上。
4.根据权利要求1所述的超声波观测装置,其特征在于,
还具备特征量图像数据生成部,该特征量图像数据生成部生成将关于所述校正特征量的信息与所述超声波图像一并显示的特征量图像数据。
5.根据权利要求1所述的超声波观测装置,其特征在于,
所述特征量计算部通过进行利用n次式对各所述频谱进行近似的处理来计算所述特征量,其中,n为正整数。
6.根据权利要求5所述的超声波观测装置,其特征在于,
所述特征量计算部利用一次式对所述频谱中的规定的频带进行近似,计算所述一次式的截距、斜率、以及所述频带的中间频率下的所述一次式的值即频带中心对应强度中的一个或多个作为所述特征量,且所述特征量包含所述斜率和所述频带中心对应强度中的任一个,
所述衰减率设定部基于所述斜率和所述频带中心对应强度中的任一个来设定所述最佳的衰减率。
7.根据权利要求6所述的超声波观测装置,其特征在于,
所述衰减率设定部在将所述斜率设为所述特征量的情况下,基于所述斜率来设定所述最佳的衰减率,在将所述频带中心对应强度设为所述特征量的情况下,基于所述频带中心对应强度来设定所述最佳的衰减率。
8.一种超声波观测装置的工作方法,该超声波观测装置基于由超声波探头获取到的超声波信号来生成超声波图像,该超声波探头具备向观测对象发送超声波并接收由该观测对象反射的超声波的超声波振子,该超声波观测装置的工作方法的特征在于,包括:
频率分析步骤,频率分析部通过对所述超声波信号的频率进行分析,来计算与所述超声波信号的接收深度及接收方向相应的多个频谱;
特征量计算步骤,特征量计算部分别计算所述多个频谱的特征量;
衰减率设定步骤,衰减率设定部针对将所述超声波图像分割为多个区域得到的各分割区域,使用在所述超声波在所述观测对象中传播时赋予互不相同的衰减特性的、多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施用于排除所述超声波的影响的衰减校正,由此针对每个所述衰减率候选值计算各所述频谱的预备校正特征量,基于该计算的结果来设定多个所述衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及
特征量校正步骤,特征量校正部使用针对每个所述分割区域设定的最佳的衰减率中的、所述超声波振子的表面与采样点之间存在的所述分割区域的最佳的衰减率,来计算所述采样点处的每单位频率的累积衰减率,使用该累积衰减率来进行所述特征量的衰减校正,由此计算校正特征量,
其中,所述衰减率设定部针对每个所述衰减率候选值计算所述预备校正特征量的统计性偏差,将该统计性偏差最小的衰减率候选值设为所述最佳的衰减率。
9.一种计算机可读存储介质,其上存储有超声波观测装置的工作程序,该超声波观测装置基于由超声波探头获取到的超声波信号来生成超声波图像,该超声波探头具备向观测对象发送超声波并接收由该观测对象反射的超声波的超声波振子,其特征在于,该超声波观测装置的工作程序被所述计算机执行时使超声波观测装置执行以下步骤:
频率分析步骤,频率分析部通过对所述超声波信号的频率进行分析,来计算与所述超声波信号的接收深度及接收方向相应的多个频谱;
特征量计算步骤,特征量计算部分别计算所述多个频谱的特征量;
衰减率设定步骤,衰减率设定部针对将所述超声波图像分割为多个区域得到的各分割区域,使用在所述超声波在所述观测对象中传播时赋予互不相同的衰减特性的、多个每单位长度和每单位频率的衰减率候选值中的每个衰减率候选值,来对各频谱的特征量实施排除所述超声波的影响的衰减校正,由此针对每个所述衰减率候选值计算所述各频谱的预备校正特征量,基于该计算的结果来设定多个所述衰减率候选值中的对于所述观测对象而言最佳的衰减率;以及
特征量校正步骤,特征量校正部使用针对每个所述分割区域设定的最佳的衰减率中的、所述超声波振子的表面与采样点之间所存在的所述分割区域的最佳的衰减率,来计算所述采样点处的每单位频率的累积衰减率,使用该累积衰减率来进行所述特征量的衰减校正,由此计算校正特征量,
其中,所述衰减率设定部针对每个所述衰减率候选值计算所述预备校正特征量的统计性偏差,将该统计性偏差最小的衰减率候选值设为所述最佳的衰减率。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014259474 | 2014-12-22 | ||
JP2014-259474 | 2014-12-22 | ||
PCT/JP2015/078548 WO2016103849A1 (ja) | 2014-12-22 | 2015-10-07 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105916448A CN105916448A (zh) | 2016-08-31 |
CN105916448B true CN105916448B (zh) | 2019-03-01 |
Family
ID=56149886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580004506.2A Active CN105916448B (zh) | 2014-12-22 | 2015-10-07 | 超声波观测装置以及超声波观测装置的工作方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20160317126A1 (zh) |
EP (1) | EP3238632B1 (zh) |
CN (1) | CN105916448B (zh) |
WO (1) | WO2016103849A1 (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102802536A (zh) * | 2010-11-11 | 2012-11-28 | 奥林巴斯医疗株式会社 | 超声波诊断装置、超声波诊断装置的动作方法以及超声波诊断装置的动作程序 |
CN102905624A (zh) * | 2010-11-11 | 2013-01-30 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN103945773A (zh) * | 2012-07-27 | 2014-07-23 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN104125804A (zh) * | 2012-10-01 | 2014-10-29 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
EP2853918A1 (en) * | 2013-09-24 | 2015-04-01 | Siemens Medical Solutions USA, Inc. | Shear wave detection in medical ultrasound imaging |
EP3005945A1 (en) * | 2013-05-29 | 2016-04-13 | Olympus Corporation | Ultrasonic observation device, operation method for ultrasonic observation device, and operation program for ultrasonic observation device |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040158154A1 (en) * | 2003-02-06 | 2004-08-12 | Siemens Medical Solutions Usa, Inc. | Portable three dimensional diagnostic ultrasound imaging methods and systems |
JP4585326B2 (ja) * | 2005-02-08 | 2010-11-24 | 富士フイルム株式会社 | 超音波撮像装置及び超音波撮像方法 |
WO2007003058A1 (en) * | 2005-07-06 | 2007-01-11 | National Research Council Of Canada | Method and system for determining material properties using ultrasonic attenuation |
US8187191B2 (en) * | 2009-01-08 | 2012-05-29 | Volcano Corporation | System and method for equalizing received intravascular ultrasound echo signals |
JP5355194B2 (ja) * | 2009-04-13 | 2013-11-27 | 日立アロカメディカル株式会社 | 超音波診断装置 |
US20110288410A1 (en) * | 2010-02-22 | 2011-11-24 | Speyer Gavriel A | Methods and systems for diagnostic ultrasound based monitoring of high intensity focused ultrasound therapy |
CN103200876B (zh) * | 2010-11-11 | 2015-09-09 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法 |
US8663113B2 (en) * | 2011-08-11 | 2014-03-04 | Delphinus Medical Technologies, Inc. | Method for imaging a volume of tissue |
US10512450B2 (en) * | 2013-09-24 | 2019-12-24 | Siemens Medical Solutions Usa, Inc. | Shear wave estimation from analytic data |
-
2015
- 2015-10-07 WO PCT/JP2015/078548 patent/WO2016103849A1/ja active Application Filing
- 2015-10-07 EP EP15872418.7A patent/EP3238632B1/en active Active
- 2015-10-07 CN CN201580004506.2A patent/CN105916448B/zh active Active
-
2016
- 2016-07-12 US US15/207,839 patent/US20160317126A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102802536A (zh) * | 2010-11-11 | 2012-11-28 | 奥林巴斯医疗株式会社 | 超声波诊断装置、超声波诊断装置的动作方法以及超声波诊断装置的动作程序 |
CN102905624A (zh) * | 2010-11-11 | 2013-01-30 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN103945773A (zh) * | 2012-07-27 | 2014-07-23 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
CN104125804A (zh) * | 2012-10-01 | 2014-10-29 | 奥林巴斯医疗株式会社 | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 |
EP3005945A1 (en) * | 2013-05-29 | 2016-04-13 | Olympus Corporation | Ultrasonic observation device, operation method for ultrasonic observation device, and operation program for ultrasonic observation device |
EP2853918A1 (en) * | 2013-09-24 | 2015-04-01 | Siemens Medical Solutions USA, Inc. | Shear wave detection in medical ultrasound imaging |
Also Published As
Publication number | Publication date |
---|---|
EP3238632A1 (en) | 2017-11-01 |
EP3238632B1 (en) | 2019-08-07 |
CN105916448A (zh) | 2016-08-31 |
US20160317126A1 (en) | 2016-11-03 |
WO2016103849A1 (ja) | 2016-06-30 |
EP3238632A4 (en) | 2018-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106170253B (zh) | 超声波观测装置以及超声波观测装置的工作方法 | |
JP5974210B2 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
CN105939674B (zh) | 超声波诊断装置以及超声波诊断装置的工作方法 | |
WO2018142937A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
US10201329B2 (en) | Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and computer-readable recording medium | |
CN108472018A (zh) | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 | |
CN108366784A (zh) | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 | |
JP6892320B2 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP6253869B2 (ja) | 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム | |
CN105939672B (zh) | 超声波观测装置以及超声波观测装置的工作方法 | |
CN105916448B (zh) | 超声波观测装置以及超声波观测装置的工作方法 | |
WO2016181869A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP6157790B1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JPWO2015198712A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP5981072B1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP5927367B1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP6138402B2 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP5953457B1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
CN107530057A (zh) | 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序 | |
JP6010274B1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JPWO2015198713A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
JP2017217313A (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |