CN102596048B - 超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置 - Google Patents

超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置 Download PDF

Info

Publication number
CN102596048B
CN102596048B CN201180002843.XA CN201180002843A CN102596048B CN 102596048 B CN102596048 B CN 102596048B CN 201180002843 A CN201180002843 A CN 201180002843A CN 102596048 B CN102596048 B CN 102596048B
Authority
CN
China
Prior art keywords
mentioned
correcting unit
gray correction
histogrammic
correction function
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
Application number
CN201180002843.XA
Other languages
English (en)
Other versions
CN102596048A (zh
Inventor
佐佐木琢也
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN102596048A publication Critical patent/CN102596048A/zh
Application granted granted Critical
Publication of CN102596048B publication Critical patent/CN102596048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration by the use of histogram techniques
    • G06T5/92
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/483Diagnostic techniques involving the acquisition of a 3D volume of data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20004Adaptive image processing
    • G06T2207/20008Globally adaptive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30044Fetus; Embryo

Abstract

本发明涉及的超声波诊断装置,具备:图像数据取得单元,对被检体的包含诊断对象的规定区域内发送超声波,接收来自规定区域内的反射波,并根据反射波取得超声波图像数据;校正单元,执行超声波图像数据的灰度校正,校正单元在灰度校正中,算出与图像数据的亮度相关的直方图,使用直方图来算出与诊断对象对应的亮度分布范围与灰度校正函数。

Description

超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置
技术领域
本发明涉及根据使用超声波探头进行超声波扫描而得到的信号生成并显示图像时,以进行适当的灰度校正为目的而被使用的超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置。
背景技术
超声波诊断装置是对被检体内放射从设置于超声波探头的振动元件发生的超声波脉冲,由上述振动元件接收由于被检体组织的声阻抗的差异而生成的超声波反射波并收集生物体信息的装置。能够以只使超声波探头与体表接触的简单的操作来实时显示图像数据,例如,由于可以观察心脏等存在活动的对象物,因此,广泛用于循环器官区域、各种脏器的形态诊断或功能诊断中。为了易于观察这样的由超声波诊断装置取得的超声波图像,需要根据对象物适当设定、调整亮度的浓淡(灰度)等。
发明内容
但是,由于生物体内的构造、部位(脏器),在诊断中不需要的噪音等变多,存在难以进行用于使超声波图像的对比度适当的调整,或在用户自身可调整的范围内不能够进行适当地调整的情况等。然而,现有的超声波诊断装置操作性较低,例如,在上述各情况下,不能适当地帮助图像灰度的设定、调整。
本发明是鉴于上述情形而完成的,其目的在于,提供一种与现有的超声波诊断装置相比较,能够适当设定、调整图像的灰度的超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置。
一实施方式涉及的超声波诊断装置,其特征在于,具备:图像数据取得单元,对被检体的包含诊断对象的规定区域内发送超声波,接收来自上述规定区域内的反射波,并根据上述反射波取得超声波图像数据;以及校正单元,执行上述超声波图像数据的灰度校正,上述校正单元在上述灰度校正中,算出与上述图像数据的亮度相关的直方图,使用上述直方图来算出与上述诊断对象对应的亮度分布范围与灰度校正函数。
根据上述本发明,可以实现一种与现有的超声波诊断装置相比较,能够适当地设定、调整图像的灰度的超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置。
附图说明
图1是表示实施方式涉及的超声波诊断装置1的方框结构图。
图2是表示在依照本灰度优化功能的处理中执行的各处理的流程的流程图。
图3是表示在依照本灰度优化功能的处理中执行的各处理的流程的流程图。
图4是表示在灰度优化处理中取得的直方图及直方图的累积和的一个例子的图。
图5是表示在灰度优化处理中取得的直方图的微分的一个例子的图。
图6是表示在灰度优化处理中取得的直方图的二阶微分的一个例子的图。
图7是用于说明使用了特征点的亮度分布范围的设定处理的图。
图8是表示未实施本灰度优化处理的超声波图像的一个例子的图。
图9是表示对于图8的图像实施本灰度优化处理的超声波图像的一个例子的图。
图10是表示采用使第2零点P2中的平滑化累积值的斜率标准化而算出的A2的斜率的变形例的图。
图11是用于说明本实施方式的变形例2的图。
图12是用于说明本实施方式的变形例3的图。
图13是用于说明本实施方式的变形例3的图。
图14是用于说明本实施方式的变形例3的图。
图15是使本灰度优化处理应用于运动图像显示时(变形例5涉及的情况)的流程图的一个例子。
图16是表示与规定的超声波图像的增益调整前的亮度值相关的直方图的一个例子的图。
图17是表示图16所示的直方图的一阶导数的曲线的图。
图18是表示图16所示的直方图的二阶导数的曲线的图。
图19是表示图16所示的直方图的三阶导数的曲线的图。
图20是与提高了增益时的超声波图像的亮度值相关的直方图的一个例子的图。
图21是表示图20所示的直方图的一阶导数的曲线的图。
图22是表示图20所示的直方图的二阶导数的曲线的图。
图23是表示图20所示的直方图的三阶导数的曲线的图。
图24是表示与降低了增益时的超声波图像的亮度值相关的直方图的一个例子的图。
图25是表示图24所示的直方图的一阶导数的曲线的图。
图26是表示图24所示的直方图的二阶导数的曲线的图。
图27是表示图24所示的直方图的三阶导数的曲线的图。
符号说明
10...超声波诊断装置  12...超声波探头  13...输入装置  14...监视器  21...超声波发送单元  22...超声波接收单元  23...B模式处理单元  24...血流检测单元  25...RAW数据存储器  26...体数据生成单元  27...管腔附近血流描绘单元  28...图像处理单元  29...控制处理器  30...显示处理单元  31...存储单元  32...接口单元
具体实施方式
下面,参照附图对本实施方式进行说明。另外,在以下的说明中,对于具有大致相同功能及结构的构成要素,添加同一符号,只在必要时进行重复说明。
(第1实施方式)
图1是表示本实施方式涉及的超声波诊断装置1的方框结构图。如该图所示,本超声波诊断装置1具备超声波探头12、输入装置13、监视器14、超声波发送单元21、超声波接收单元22、B模式处理单元23、血流检测单元24、RAW数据存储器25、图像处理单元26、具有灰度优化单元280的显示处理单元28、控制处理器(CPU)29、存储单元30、接口单元31。以下。针对各个构成要素的功能进行说明。
超声波探头12是对于被检体发送超声波,并接收基于该发送的超声波的来自被检体的反射波的装置(探针),具有多个在其前端排列的压电振子、匹配层背衬材料等。超声波探头12根据来自超声波发送单元21的驱动信号对扫描区域内的所希望的方向发送超声波,压电振子将来自该被检体的反射波转换成电信号。匹配层设置于该压电振子上,是用于使超声波能量有效地传播的中间层。背衬材料防止超声波从该压电振子向后方传播。如果从该超声波探头12对被检体P发送超声波,则该发送超声波在体内组织的声阻抗的不连续面上依次被反射,并作为回波信号被超声波探头12接收。该回波信号的振幅依存于反射时的不连续面上的声阻抗的差。另外,被发送的超声波脉冲因移动的血流而被反射时的回波由于多普勒效应依存于移动体的超声波发送方向的速度分量,而受到频移。
另外,在取得体数据时,作为超声波探头12,例如可以采用二维阵列探头(多个超声波振子二维矩阵状排列的探头)、或机械4D探头(能够一边使超声波振子列机械地在正交于其排列方向的方向上摇动一边执行超声波扫描的探头)。但是,并不拘泥于该例子,例如也可以将一维阵列探头作为超声波探头12来使用,即使将其一边手动地进行摆动,一边进行超声波扫描,也能够取得体数据。
输入装置13与装置主体11连接,具有用于将来自操作者的各种指示、条件、关心区域(ROI)的设定指示、各种画质条件设定指示等取入装置主体11的各种开关、按钮、轨迹球、鼠标、键盘等。另外,输入装置13在后述的灰度优化功能中,具有其开始指示或作为对象的超声波图像的选择指示等。
监视器14根据来自图像处理单元28的视频信号,将生物体内的形态学信息或血流信息作为图像来显示。
超声波发送单元21具有未图示的触发发生电路、延迟电路及脉冲发生器电路等。在触发发生电路中,以规定的额定频率fr Hz(周期;1/fr秒)反复发生用于形成发送超声波的触发脉冲。另外,在延迟电路中,针对每个通道将超声波会集成束状并决定发送指向性所需的延迟时间提供给各触发脉冲。脉冲发生器电路以基于该触发脉冲的定时来对探头12施加驱动脉冲。
另外,超声波发送单元21为了按照控制处理器28的指示执行规定的扫描序列,而具有能够瞬时变更发送频率、发送驱动电压等功能。特别对于发送驱动电压的变更,通过能够瞬时切换其值的线性放大器型的发送电路、或电切换多个电源单元的机构来实现。
超声波接收单元22具有未图示的放大电路、A/D转换器、加法器等。在放大电路中,针对每个通道将经由探头12取入的回波信号放大。在A/D转换器中,对于被放大的回波信号决定接收指向性,并提供进行接收动态聚焦所需的延迟时间,之后在加法器中进行相加处理。通过该相加,强调来自与回波信号的接收指向性相应的方向的反射分量,并根据接收指向性与发送指向性形成超声波发送接收的综合性的波束。
B模式处理单元23从接收单元22来接收回波信号,并实施对数放大、包络线检波处理等,并生成由亮度的明暗来表现信号强度的数据。
血流检测单元24根据从接收单元22接收到的回波信号来检测血流信号,并生成血流数据。血流信号的检测通常由CFM(Color Flowmapping:彩色血流成像)来进行。此时,解析血流信号,作为血流数据就多点求得平均速度、分散、能量等血流信息。
RAW数据存储器25使用从B模式处理单元23接收到的B模式数据、从血流检测单元24接收到的血流数据,分别生成每一帧的B模式RAW数据、血流RAW数据。另外,RAW数据存储器25根据需要通过执行RAW体素转换,从而从RAW数据生成体数据。
图像处理单元28对于从RAW数据存储器25接收到的RAW数据,执行扫描转换处理等。另外,图像处理单元28对于从RAW数据存储器25接收到的体数据,进行体绘制、多剖面转换显示(MPR:multi planar reconstruction:多平面重建)、最大值投影显示(MIP:maximum intensity Projection)等规定的图像处理。另外,为了使噪音降低或图像的连接优良,可以在图像处理单元28之后插入二维滤波器,进行空间性的平滑处理。
显示处理单元28对于在图像处理单元28中生成、处理的各种图像数据,执行各种动态范围、亮度(brightness)、对比度、γ曲线校正、RGB转换等。并且,显示处理单元28所具有的灰度优化单元280按照来自控制处理器29的控制来执行后述的依照灰度优化功能的处理。
控制处理器29具有作为信息处理装置(计算机)的功能,控制本超声波诊断装置主体的动作。控制处理器29从存储单元30读出用于实现后述的灰度优化功能的专用程序,并在自身所具有的存储器上展开,执行与各种处理相关的运算、控制等。
存储单元30保管有用于实现后述的灰度优化功能的专用程序、诊断信息(患者ID、医师的评论等)、诊断协议、发送接收条件、用于实现散斑除去功能的程序、体标生成程序以及其他数据组。另外,根据需要,也可以用于保管RAW数据存储器中的图像等。存储单元30的数据能够经由接口单元31向外部周边装置传送。
接口单元31是与输入装置13、网络、新的外部存储装置(未图示)相关的接口。由该装置得到的超声波图像等数据和解析结果等能够通过接口单元31,经由网络传送给其他装置。
(灰度优化功能)
然后,针对本超声波诊断装置1所具有的灰度优化功能进行说明。该功能根据使用所取得的各个超声波图像数据而生成的直方图、该直方图的一阶微分、二阶微分来算出特征点,并设定与诊断对象对应的亮度分布范围。另外,使用算出的特征点来设定控制条件,并使用设定的控制条件算出每个图像的灰度校正函数。
图2、图3是表示在依照本灰度优化功能的处理(灰度优化处理)中执行的各处理的流程的流程图。以下,针对各步骤中的处理的内容进行说明。
另外,在以下的说明中,为了使说明具体化,以对于通过B模式对包含心肌的二维区域进行超声波扫描所得到的超声波图像进行灰度优化处理时为例进行说明。另外,在超声波图像处理装置中实现本灰度优化处理时,例如,预先存储在步骤S2中生成的超声波图像数据,然后执行步骤S3以后的处理。此时的超声波图像处理装置具备图1的虚线内的结构。
[输入接受患者信息、发送接收条件:步骤S1]
经由输入装置13执行患者信息的输入,发送接收条件(被扫描区域的大小、焦点位置、发送电压等)、用于将被检体的规定区域进行超声波扫描的摄影模式、扫描序列等的选择(步骤S1)。输入、选择的各种信息、条件等被自动地存储在存储单元30中。
[超声波扫描、超声波图像的生成:步骤S2、S3]
超声波探头12抵接在被检体表面的希望位置上,将包含诊断部位(此时,为心脏)的二维区域作为被扫描区域,执行基于B模式的超声波扫描(步骤S2)。由基于B模式的超声波扫描而取得的回波信号依次经由超声波接收单元22被送至B模式处理单元23。B模式处理单元23执行对数放大处理、包络线检波处理等并生成多个B模式数据。RAW数据存储器25使用多个从B模式处理单元23接收到的B模式数据来生成B模式RAW数据。图像处理单元26通过对于生成的B模式RAW数据执行扫描转换来生成每一帧的超声波图像(步骤S3)。
[灰度优化处理:步骤S4]
然后,灰度优化单元280对于在图像处理单元26中生成的超声波图像执行图3所示那样的灰度优化处理。即,灰度优化单元280首先生成与从图像处理单元26接收到的超声波图像的亮度相关的直方图,并对于得到的直方图执行平滑化处理(步骤S41)。之后,灰度优化单元280使用得到的直方图,如图4所示那样算出直方图的累积和,并执行平滑化处理(步骤S42)。
然后,灰度优化单元280算出直方图的微分(或差分)之后,执行平滑化处理,得到图5所示的那样的结果(步骤S43)。并且,灰度优化单元280算出直方图的二阶微分(或二阶差分)之后,执行平滑化处理,得到图6所示的那样的结果(步骤S44)。
然后,灰度优化处理280使用算出的直方图的二阶微分,算出直方图的特征点(步骤S45)。在本实施方式中,作为直方图的特征点的例子,假设使用直方图二阶微分的负的最低(极大)值(相当于直方图的正的峰值位置(最大值或极大值))、直方图二阶微分的第2零点(直方图的底部区域中的拐点:相当于直方图微分的负的峰值位置)、直方图的底部区域中的拐点、直方图二阶微分的第3零点(例如相当于来自心肌区域的信号与来自其他区域的信号的分支点)这4个点。
然后,灰度优化处理280使用算出的特征点来设定亮度分布范围(步骤S46)。例如,如图7所示,灰度优化单元280设定将直方图二阶微分的负的最低值(极大值)P1作为下限,将直方图二阶微分的第3零点P4作为上限的亮度分布范围L。
然后,灰度优化单元280使用算出的特征点、设定的控制条件及亮度分布范围L来算出灰度特性函数(步骤S47)。如果列举出一个例子,灰度优化单元280根据以下的控制条件(1)~(5)来算出灰度特性函数。
(1)使直方图二阶微分的负的最低值(极大值)P1中的平滑化累积值的斜率标准化,并作为灰度特性函数的斜率。
(2)使输入最小值与输出最小值对应(图7下方的点B1)。
(3)噪音水平上限作为拐点(二阶微分的第2零点),与所希望的输出值对应。(图7下方的点B2)
(4)对象物区域(此时为心肌区域)根据噪音水平条件(拐点),作为第3零点,第3零点的输出值与由所希望的斜率求得的值对应(图7下方的点B3)。
(5)输入最大值与输出最大值对应(图7下方的点B4)。
然后,灰度优化单元280将在步骤S47中得到的灰度特性函数进行插补(例如,样条插补),并算出例如图7下方的灰度校正函数C1(步骤S48),使用算出的灰度校正函数C1,对于超声波图像执行灰度优化处理(步骤S49)。
[超声波图像显示:步骤S5]
然后,监视器14将被实施了灰度优化处理的超声波图像以规定的形态进行显示(步骤S5)。
图8是表示未实施本灰度优化处理的超声波图像的一个例子的图,图9是表示对于图8的图像实施了本灰度优化处理了的超声波图像的一个例子的图。比较图8与图9得知,通过本灰度优化处理进行灰度调整的超声波图像变成对比度清晰、易于观察的图像。
另外,灰度优化处理的内容并不拘泥于上述例子。以下,针对本灰度优化处理的变形例进行说明。
(变形例1)
例如,如图10所示,也能够采用使第2零点P2的平滑化累积值的斜率标准化而算出的A2的斜率,算出灰度校正函数C2。另外,也可以采用使与对象物区域对应的范围(例如亮度分布范围P2-P3)的斜率平均而算出的斜率。
(变形例2)
对于上述(1)~(5),也可以追加进一步的控制条件并算出灰度特性函数。例如,可以增加以抑制高亮度区域的斜率为目的的以下控制条件。
(6)选择存在于从第3零点到输入最大值(例如灰度为256)之间的输入值,使该输入值与所希望的输出值对应(例如,图11下方的点B5)。
另外,也可以将这样的控制条件(6)的追加,在算出灰度特性函数中相对于控制条件(1)~(5)并列地追加。另外,也可以基于使用控制条件(1)~(5)算出的灰度特性函数的结果,根据需要追加控制条件(6),并重新算出灰度特性函数。另外,例如当与追加的控制条件对应的位置的灰度特性低于255时,为了抑制特性的混乱,优选校正为灰度特性函数变得单调增加。
(变形例3)
也可以根据需要进行灰度扩展处理。作为是否进行灰度扩展处理的判定例子,例如,当输入图像的直方图累积和成为规定量以上的位置达到比进行平滑处理的端部处理区域还往前的情况下,设为图像取得时的动态范围(DR)或灰度特性不恰当,例如在步骤S47的算出中,进行图12至图13所示的那样的灰度特性的扩展。
(变形例4)
亮度分布范围可以任意地进行变更。例如,可以将上述图7的亮度分布范围L(下限:直方图二阶微分的负的最低值。上限:第3零点)在任意的定时变更成图14所示的亮度分布范围L1(下限:第2零点。上限:第3零点)。由此,有时可以更好地提取出对象物。
(变形例5)
本灰度优化处理也能够适用于取得并显示运动图像的情况。
图15是使本灰度优化处理应用于运动图像显示时的流程图的一个例子。与图3所示的流程图相比较,在追加了步骤S40“直方图解析图像的选择处理”这一点不同。
在直方图解析图像的选择处理(步骤S40)中,在取得的多个超声波图像(一次心跳或多次心跳时的心脏的超声波图像)中,选择任一图像并执行步骤S41以后的灰度优化处理。直方图解析图像例如可以选择与预先设定的时相或定时对应的图像、与规定期间的最初、中间或最后的时相对应的图像、图像的亮度平均最大的图像、图像的灰度平均为多个图像中的平均数的图像等。虽然将哪一图像作为直方图解析图像是自由的,但当将诊断对象设为心脏时,由于直方图累积和配合收缩、扩张而周期性地进行变动等原因,应该不要选择较暗的图像及为将对象物适当进行成像的图像等。
(变形例6)
本灰度优化处理例如可以每当取得超声波图像数据时执行。另外,不需要对于所有的图像总是执行时,例如可以应答从用户经由输入装置13输入的指示,在任意的定时来执行。
(应用例)
在以上所述的实施方式中,示出了将与超声波图像的亮度值相关的直方图的二阶导数的第1极值、第2、第3零点等作为特征点而算出的例子。但是,本灰度优化功能不拘泥于该例,例如,也可以将n阶导数的第i极值、第j零点(其中,n,i,j分别为任意的自然数)作为特征点来算出,并使用得到的特征点算出亮度分布范围、灰度校正曲线。
根据以上所述的结构,可以得到以下的效果。
根据本超声波诊断装置,根据使用取得超声波图像数据而生成的直方图、该直方图的一阶微分、二阶微分等算出特征点,并设定与诊断对象对应的亮度分布范围。另外,使用算出的特征点设定控制条件,并使用设定的控制条件算出灰度校正函数。从而,可以将与各个图像的灰度特性相应的适当的灰度分布范围、灰度校正函数进行定义。其结果可以根据诊断对象或个体差、摄影状况,迅速且简单地将灰度校正进行优化。
(第2实施方式)
然后,针对第2实施方式涉及的超声波诊断装置进行说明。在进行成为本灰度优化处理的对象的超声波图像的增益调整时,与增益调整后的亮度值相关的直方图的形状发生变化,其结果直方图的累积和、n阶导数(其中,n为自然数)的形状也发生变化。在本实施方式中,针对这样直方图的形状由于增益调整而发生变化时的灰度优化处理进行说明。
假设如下情况:取得规定的超声波图像,得到如图16那样与该图像的增益调整前的亮度值相关的直方图。此时,直方图的一阶导数、二阶导数、三阶导数分别如图17、图18、图19所示的那样。
对于现在的超声波图像(与图16的直方图对应的超声波图像)实施提高增益的处理时,与提高增益后的亮度值相关的直方图、直方图的一阶导数、二阶导数、三阶导数分别变成图20、图21、图22、图23所示的那样。在此,例如,如果将图16与图20相比较,则得知提高了增益的结果使直方图的峰值位置发生移动(向右方向偏移)。另外,如果分别比较图17与图21、图18与图22、图19与图23,则得知提高了增益的结果,在各阶导数中各极值与各零点的位置也发生移动(向右方向偏移)。特别是如果观察图19与图23,则在图19中最初出现的极值的极性为负,与此相对,在图23中最初出现的极值的极性变成正。因而,例如,增益调整前“采用图19所示的三阶导数的第1、第2极值作为特征点”的条件与增益调整后“采用图23所示的增益调整前的三阶导数的第2、第3极值作为特征点”的条件对应。
另外,对于现在的超声波图像实施降低增益的处理时,与降低了增益后的亮度值相关的直方图、直方图的一阶导数、二阶导数、三阶导数分别变成图24、图25、图26、图27所示的那样。在此,例如,如果将图16与图24相比较,则得知降低了增益的结果使直方图的峰值位置向左方向发生偏移。另外,如果分别比较图17与图25、图18与图26、图19与图27,则得知降低了增益的结果在各次(各阶)导数中,各极值与各零点的位置也向左方向发生偏移。特别地,如果在增益上升前后分别比较一阶导数、二阶导数、三阶导数,则得知相同顺序的极值的极性分别发生逆转(例如,在一阶导数中,增益上升前的最初的极值的极性为正,与此相对,增益上升后的最初的极值的极性变为负)。因而,例如,增益调整前“采用图18所示的二阶导数的第2、第3零点作为特征点”的条件与增益调整后“采用图26所示的增益调整前的二阶导数的第1、第2零点作为特征点”的条件对应。
即,直方图的峰值位置、直方图的n阶导数的各极值的极性、各零点的位置根据增益调整而发生变化。因而,灰度优化处理中的算出特征点中所使用的第i个极值、第j个极值的零点(i,j为自然数)最好根据增益校正后的直方图、n阶导数的位置、形状而单独地适当地进行选择。
另外,本发明并不限定于上述实施方式本身,在实施阶段在不脱离本发明的要旨范围内,可以将构成要素进行变形并进行具体化。作为具体的变形例,例如存在以下形式。
(1)本实施方式涉及的各功能可以通过将执行该处理的程序安装在工作站等计算机上,并将他们在存储器上展开来实现。此时,可以使计算机执行该方法的程序能够存储并分布在磁盘(软盘(注册商标)、硬盘等)、光盘(CD-ROM、DVD等)、半导体存储器等存储介质中。
(2)在上述实施方式中,以对于利用B模式将包含心肌的二维区域进行超声波扫描而得到的超声波图像进行灰度优化处理时为例进行说明。但是,本灰度优化功能并不拘泥于该例子。例如,也可以适用于以心脏以外的部位作为诊断对象的情况、不利用B模式而利用血流检测模式(CFM模式等)进行摄影的情况、超声波扫描三维区域的情况等情况。特别地,在超声波扫描三维区域时,能够对于构成体数据的各二维图像、MPR图像等,分别应用上述灰度优化处理,从而取得最优的超声波图像。
(3)在上述各实施方式中,以对于由超声波诊断装置取得的超声波图像数据执行灰度优化处理的情况为例进行了说明。但是,在本灰度优化处理中对象并不限定于超声波图像数据,也能够适用于通过X射线诊断装置、X射线计算机断层摄影装置、磁共振成像装置等医用图像诊断装置取得的各医用图像。另外,也可以将由各种医用图像诊断装置取得的图像使用医用图像处理装置,事后进行灰度优化处理。
另外,可以根据上述实施方式中公开的多个的构成要素的适宜的组合,形成各种的发明。例如:既可以从实施方式中显示的全部构成要素中削除几个构成要素。进而还可以适当地组合不同实施方式内的构成要素。

Claims (30)

1.一种超声波诊断装置,其特征在于,具备:
图像数据取得单元,对被检体的包含诊断对象的规定区域内发送超声波,接收来自上述规定区域内的反射波,并根据上述反射波取得超声波图像数据;以及
校正单元,执行上述超声波图像数据的灰度校正,
上述校正单元在上述灰度校正中,
算出与上述超声波图像数据的亮度相关的直方图,
使用上述直方图来算出与上述诊断对象对应的亮度分布范围与灰度校正函数。
2.根据权利要求1所述的超声波诊断装置,其特征在于,
上述校正单元算出上述直方图的累积和,利用上述亮度分布范围中的上述累积和的斜率来算出上述灰度校正函数的斜率。
3.根据权利要求1所述的超声波诊断装置,其特征在于,
上述校正单元使用从上述直方图的二阶微分得到的特征点来算出上述亮度分布范围与上述灰度校正函数。
4.根据权利要求3所述的超声波诊断装置,其特征在于,
上述校正单元将上述直方图的拐点作为上述特征点,算出上述亮度分布范围与上述灰度校正函数。
5.根据权利要求3所述的超声波诊断装置,其特征在于,
上述校正单元将上述直方图的二阶微分的最低值、最大值、第2零点、第3零点中的至少一个作为上述特征点,算出上述灰度校正函数。
6.根据权利要求3所述的超声波诊断装置,其特征在于,
上述校正单元将经由输入装置指定的任意点作为上述特征点,算出上述灰度校正函数。
7.根据权利要求1所述的超声波诊断装置,其特征在于,
上述校正单元使用从上述直方图的n阶微分得到的特征点,算出上述亮度分布范围与上述灰度校正函数,其中,n为自然数。
8.根据权利要求7所述的超声波诊断装置,其特征在于,
上述校正单元将上述直方图的拐点作为上述特征点,算出上述亮度分布范围与上述灰度校正函数。
9.根据权利要求7所述的超声波诊断装置,其特征在于,
上述校正单元将上述直方图的n阶微分的第i极值、第j零点中的至少一个作为上述特征点,算出上述灰度校正函数,其中,n,i,j分别为自然数。
10.根据权利要求7所述的超声波诊断装置,其特征在于,
上述校正单元将经由输入装置指定的任意点作为上述特征点,算出上述灰度校正函数。
11.根据权利要求1所述的超声波诊断装置,其特征在于,
在增益调整被进行了的情况下,上述校正单元算出与增益调整后的上述超声波图像数据的亮度相关的上述直方图。
12.根据权利要求1所述的超声波诊断装置,其特征在于,
上述诊断对象为心脏,
上述校正单元将与上述直方图的心肌对应的区域决定为上述亮度分布范围。
13.根据权利要求1所述的超声波诊断装置,其特征在于,
还具备输入单元,用于输入上述灰度校正的开始指示,
上述校正单元应答经由上述输入单元输入的上述开始指示,执行上述灰度校正。
14.根据权利要求1所述的超声波诊断装置,其特征在于,
每当上述图像数据取得单元取得超声波图像数据,上述校正单元就会执行上述灰度校正。
15.一种超声波图像处理装置,其特征在于,具备:
存储单元,对被检体的包含诊断对象的规定区域内发送超声波,接收来自上述规定区域内的反射波,并存储根据上述反射波取得的超声波图像数据;以及
校正单元,执行上述超声波图像数据的灰度校正,
上述校正单元在上述灰度校正中,
算出与上述超声波图像数据的亮度相关的直方图;
使用上述直方图算出与上述诊断对象对应的亮度分布范围与灰度校正函数。
16.根据权利要求15所述的超声波图像处理装置,其特征在于,
上述校正单元算出上述直方图的累积和,利用上述亮度分布范围中的上述累积和的斜率算出上述灰度校正函数的斜率。
17.根据权利要求15所述的超声波图像处理装置,其特征在于,
上述校正单元使用从上述直方图的二阶微分得到的特征点,算出上述亮度分布范围与上述灰度校正函数。
18.根据权利要求17所述的超声波图像处理装置,其特征在于,
上述校正单元以上述直方图的拐点作为上述特征点,算出上述亮度分布范围与上述灰度校正函数。
19.根据权利要求17所述的超声波图像处理装置,其特征在于,
上述校正单元将上述直方图的二阶微分的最低值、最大值、第2零点、第3零点中的至少一个作为上述特征点,算出上述灰度校正函数。
20.根据权利要求17所述的超声波图像处理装置,其特征在于,
上述校正单元将经由输入装置指定的任意点作为上述特征点,算出上述灰度校正函数。
21.根据权利要求15所述的超声波图像处理装置,其特征在于,
上述校正单元使用从上述直方图的n阶微分得到的特征点,算出上述亮度分布范围与上述灰度校正函数,其中,n为自然数。
22.根据权利要求21所述的超声波图像处理装置,其特征在于,
上述校正单元将上述直方图的拐点作为上述特征点,算出上述亮度分布范围与上述灰度校正函数。
23.根据权利要求21所述的超声波图像处理装置,其特征在于,
上述校正单元将上述直方图的n阶微分的第i极值、第j零点中的至少一个作为上述特征点,算出上述灰度校正函数,其中,n,i,j分别为自然数。
24.根据权利要求21所述的超声波图像处理装置,其特征在于,
上述校正单元将经由输入装置指定的任意点作为上述特征点,算出上述灰度校正函数。
25.根据权利要求15所述的超声波图像处理装置,其特征在于,
上述诊断对象为心脏,
上述校正单元将与上述直方图的心肌对应的区域决定为上述亮度分布范围。
26.根据权利要求15所述的超声波图像处理装置,其特征在于,
还具备输入单元,用于输入上述灰度校正的开始指示,
上述校正单元应答经由上述输入单元输入的上述开始指示,执行上述灰度校正。
27.根据权利要求15所述的超声波图像处理装置,其特征在于,
每次上述图像数据取得单元取得超声波图像数据,上述校正单元执行上述灰度校正。
28.根据权利要求15所述的超声波图像处理装置,其特征在于,
在增益调整被进行了的情况下,上述校正单元算出与增益调整后的上述超声波图像数据的亮度相关的上述直方图。
29.一种医用图像诊断装置,其特征在于,具备:
图像数据取得单元,取得与被检体的包含诊断对象的规定区域相关的图像数据;以及
校正单元,执行上述图像数据的灰度校正,
上述校正单元在上述灰度校正中,
算出与上述图像数据的亮度相关的直方图,
使用上述直方图算出与上述诊断对象对应的亮度分布范围与灰度校正函数。
30.一种医用图像处理装置,其特征在于,具备:
存储单元,存储与被检体的包含诊断对象的规定区域相关的图像数据;以及
校正单元,执行上述图像数据的灰度校正,
上述校正单元在上述灰度校正中,
算出与上述图像数据的亮度相关的直方图,
使用上述直方图算出与上述诊断对象对应的亮度分布范围与灰度校正函数。
CN201180002843.XA 2010-10-26 2011-10-26 超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置 Active CN102596048B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2010240112 2010-10-26
JP2010-240112 2010-10-26
PCT/JP2011/074684 WO2012057211A1 (ja) 2010-10-26 2011-10-26 超音波診断装置、超音波画像処理装置、医用画像診断装置及び医用画像処理装置

Publications (2)

Publication Number Publication Date
CN102596048A CN102596048A (zh) 2012-07-18
CN102596048B true CN102596048B (zh) 2014-12-10

Family

ID=45993917

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180002843.XA Active CN102596048B (zh) 2010-10-26 2011-10-26 超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置

Country Status (5)

Country Link
US (1) US9743910B2 (zh)
EP (1) EP2633818B1 (zh)
JP (1) JP5943373B2 (zh)
CN (1) CN102596048B (zh)
WO (1) WO2012057211A1 (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5943373B2 (ja) * 2010-10-26 2016-07-05 東芝メディカルシステムズ株式会社 超音波診断装置、超音波画像処理装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム
JP5788624B2 (ja) * 2013-07-18 2015-10-07 オリンパス株式会社 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP6274495B2 (ja) 2013-10-25 2018-02-07 キヤノンメディカルシステムズ株式会社 画像処理装置および超音波診断装置
CN104918029A (zh) * 2014-03-11 2015-09-16 南京巨鲨显示科技有限公司 增强显示器图像景深的方法
US9460499B2 (en) * 2014-05-30 2016-10-04 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Systems and methods for selective enhancement of a region of interest in an image
EP2989986B1 (en) * 2014-09-01 2019-12-18 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
CN104490418A (zh) * 2014-09-25 2015-04-08 深圳市恩普电子技术有限公司 一种基于信号统计分析的超声图像自动优化方法
CN107894204B (zh) * 2016-10-04 2020-02-21 财团法人工业技术研究院 干涉仪及其成像方法
DE102017107151A1 (de) * 2017-04-03 2018-10-04 Herrmann Ultraschalltechnik Gmbh & Co. Kg Ultraschallbearbeitungsmaschine mit zwei Sonotroden und Verfahren zum Betreiben einer solchen
CN107582099B (zh) * 2017-09-22 2019-12-27 杭州创影健康管理有限公司 回声强度处理方法、装置及电子设备
CN110200640B (zh) * 2019-05-14 2022-02-18 南京理工大学 基于双模态传感器的非接触式情绪识别方法
KR102355438B1 (ko) * 2019-11-28 2022-01-26 주식회사 씨자인 초음파 센서를 이용한 사물 인식 시스템 및 그 인식 방법

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101467892A (zh) * 2007-12-26 2009-07-01 株式会社东芝 超声波诊断装置、超声波图像处理装置及方法

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4662380A (en) * 1985-10-31 1987-05-05 General Electric Company Adaptive time gain compensation system for ultrasound imaging
US4868651A (en) * 1988-05-17 1989-09-19 S&S Inficon, Inc. Digital radiography with image brightness and contrast normalization
US5313948A (en) * 1991-11-28 1994-05-24 Aloka Co., Ltd. Ultrasonic diagnostic apparatus
JP3696616B2 (ja) * 1994-05-27 2005-09-21 フクダ電子株式会社 超音波診断装置
JPH08206117A (ja) 1994-05-27 1996-08-13 Fujitsu Ltd 超音波診断装置
US5579768A (en) * 1995-03-21 1996-12-03 Acuson Corporation Automatic gain compensation in an ultrasound imaging system
JP4054454B2 (ja) * 1998-09-25 2008-02-27 キヤノン株式会社 画像処理装置、画像処理方法及びコンピュータ読み取り可能な記憶媒体
JP2002209891A (ja) * 2001-01-18 2002-07-30 Matsushita Electric Ind Co Ltd 超音波装置、超音波信号変換方法、およびプログラム
KR100438903B1 (ko) * 2002-01-31 2004-07-02 한국전자통신연구원 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비
JP4626138B2 (ja) * 2003-09-30 2011-02-02 コニカミノルタエムジー株式会社 画像処理装置、画像処理方法およびプログラム
WO2005065028A2 (en) * 2004-01-07 2005-07-21 Ramot At Tel Aviv University Ltd. Methods and apparatus for analysing ultrasound images
US7248725B2 (en) * 2004-01-07 2007-07-24 Ramot At Tel Avia University Ltd. Methods and apparatus for analyzing ultrasound images
JP4575738B2 (ja) * 2004-09-29 2010-11-04 富士フイルム株式会社 超音波画像境界抽出方法及び超音波画像境界抽出装置、並びに、超音波撮像装置
US7586653B2 (en) * 2005-04-22 2009-09-08 Lexmark International, Inc. Method and system for enhancing an image using luminance scaling
JP4807824B2 (ja) * 2005-07-07 2011-11-02 株式会社日立メディコ 医用画像診断システム
US20090279764A1 (en) 2005-10-27 2009-11-12 Konica Minolta Medical & Graphic, Inc. Small-scale diagnosis system
JP2007117351A (ja) * 2005-10-27 2007-05-17 Konica Minolta Medical & Graphic Inc 小規模診断支援システム
JP4525719B2 (ja) * 2007-08-31 2010-08-18 カシオ計算機株式会社 階調補正装置、階調補正方法、及び、プログラム
JP2010154132A (ja) * 2008-12-24 2010-07-08 Toshiba Corp ガンマ補正回路およびガンマ補正方法
JP2011188391A (ja) * 2010-03-10 2011-09-22 Sony Corp 画像処理装置および方法、並びにプログラム
JP5943373B2 (ja) * 2010-10-26 2016-07-05 東芝メディカルシステムズ株式会社 超音波診断装置、超音波画像処理装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム
JP6476898B2 (ja) * 2014-03-07 2019-03-06 株式会社リコー 画像処理装置、画像処理方法、プログラム及び記憶媒体

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101467892A (zh) * 2007-12-26 2009-07-01 株式会社东芝 超声波诊断装置、超声波图像处理装置及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JP特开2002-209891A 2002.07.30 *
JP特开2005-169155A 2005.06.30 *
JP特开2007-117351A 2007.05.17 *

Also Published As

Publication number Publication date
JP2012105973A (ja) 2012-06-07
JP5943373B2 (ja) 2016-07-05
US20130237825A1 (en) 2013-09-12
EP2633818B1 (en) 2020-01-15
US9743910B2 (en) 2017-08-29
EP2633818A4 (en) 2015-09-02
EP2633818A1 (en) 2013-09-04
WO2012057211A1 (ja) 2012-05-03
CN102596048A (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN102596048B (zh) 超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置
CN106037797B (zh) 超声成像中感兴趣的三维容积
KR100908252B1 (ko) 영상 처리 시스템 및 방법
CN102247171B (zh) 超声波诊断装置、超声波图像处理装置以及医用图像诊断装置
KR102223048B1 (ko) 정량적 초음파 이미징을 위한 관심 구역 배치
CN102415902B (zh) 超声波诊断装置以及超声波图像处理装置
KR102011545B1 (ko) 초음파 진단 장치, 화상 처리 장치 및 화상 처리 방법
CN103889337A (zh) 超声波诊断装置以及超声波诊断装置控制方法
CN103169500A (zh) 超声波诊断装置、医用图像诊断装置以及医用图像处理方法
JP6309340B2 (ja) 超音波診断装置及び超音波イメージングプログラム
US11308609B2 (en) System and methods for sequential scan parameter selection
JP6687336B2 (ja) 超音波診断装置及び制御プログラム
KR20140131808A (ko) 초음파 영상 장치 및 그 제어 방법
JP2016067720A (ja) 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
US20190117196A1 (en) Ultrasound diagnostic apparatus and ultrasound probe
JP2001017428A (ja) オパシティ設定方法、3次元像形成方法および装置並びに超音波撮像装置
JP2013143978A (ja) 超音波診断装置
JP6169361B2 (ja) 超音波診断装置および輝度補正方法
JP4530834B2 (ja) 超音波画像処理方法及び超音波画像処理装置、並びに、超音波画像処理プログラム
JP6727363B2 (ja) 医用診断装置、医用画像処理装置及び医用画像処理方法
JP2022158712A (ja) 超音波診断装置、画像処理装置、及び画像処理プログラム
JP5337446B2 (ja) 超音波画像診断装置、画像処理装置及び超音波画像診断支援プログラム
CN113301854A (zh) 图像解析装置
JP7469877B2 (ja) 超音波診断装置、医用画像処理装置、および医用画像処理プログラム
JP2013244138A (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
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: 20160718

Address after: Japan Tochigi

Patentee after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Toshiba Corp

Patentee before: Toshiba Medical System Co., Ltd.