CN101209210B - 超声波图像取得装置 - Google Patents

超声波图像取得装置 Download PDF

Info

Publication number
CN101209210B
CN101209210B CN200710167811.0A CN200710167811A CN101209210B CN 101209210 B CN101209210 B CN 101209210B CN 200710167811 A CN200710167811 A CN 200710167811A CN 101209210 B CN101209210 B CN 101209210B
Authority
CN
China
Prior art keywords
mentioned
interpolation
data
coordinate system
coordinate
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
CN200710167811.0A
Other languages
English (en)
Other versions
CN101209210A (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 CN101209210A publication Critical patent/CN101209210A/zh
Application granted granted Critical
Publication of CN101209210B publication Critical patent/CN101209210B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • 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/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52044Scan converters

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Acoustics & Sound (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

在本发明的超声波图像取得装置中,发送接收部件经由超声波探头用超声波扫描被检体,取得由沿着扫描线的第一座标系中的超声波数据构成的第一超声波数据,包含在图像生成部件中的座标变换部件使用GPU进行计算,将上述第一超声波数据从上述第一座标系变换为用于图像显示的第二座标系。

Description

超声波图像取得装置
技术领域
本发明涉及在根据用超声波进行扫描而得到的信号生成图像时,从扫描时的座标向显示的座标进行座标变换进行显示的超声波图像取得装置。
背景技术
超声波图像取得装置是通过重复进行从内置于超声波探头中的压电振子向被检体内的超声波的发送接收,并进行各种处理,从而根据被检体内的超声波图像得到身体信息的装置。在超声波图像取得装置中,能够进行基于超声波断层法的诊断、基于超声波多普勒(Doppler)法的诊断。
基于超声波断层法的诊断通过用超声波对被检体内的断面进行扫描,将反射波信号的振幅变换为亮度,来得到被检体的二维断层像。通过超声波断层法得到的二维断层像被称为B模式(mode)图像。将作为作成该B模式图像的基础的数据(data)称为B模式数据。
另外,基于超声波多普勒法的诊断是利用超声波的多普勒效果测量被检体内的血流和组织的速度,或彩色地将血流显示为二维图像(以下也有称为彩色多普勒(color Doppler)像的时候)的方法。在超声波多普勒法中,特别将检测血流速度而彩色地显示血流的方法称为彩色多普勒法。
在此,在接收到的超声波信号以设定了的速度以上通过了抽出其中运动的部分的MTI(Moving Target Indicator)滤波器后,如果设所得到的复数多普勒信号为x1、x2、x3、……xN,则如下:
能量(power)P=∑|xi|2
自相关函数C1=∑xi*xi+1(在此,x*为共轭复数)
速度V=tan-1C1
分布T=1-|C1|/P。
进而,将通过彩色多普勒法得到的彩色多普勒像重叠在通过超声波断层法得到的B模式像上进行显示的方法被称为彩色多普勒断层法。
这样,在超声波图像取得装置中,通过压电振子使超声波波束(beam)大致从固定点向规定的方向偏向而进行发送和接收,使该方向顺序地变化而进行超声波扫描,由此得到二维断面内或三维空间内的数据。然后,为了根据基于取得数据时的座标系(该座标系是(扫描方向,距离)的座标系。以下,将该座标系称为“处理座标系”(scanning coordinate system),将座标称为“处理座标”(scanningcoordinate))进行超声波扫描而得到的采样(sample)点(参考图1A)而显示在图像显示用的监视器上,必须变换为作为显示监视器的栅格的正交座标系(参考图1B)。在此,图1A是用处理座标表示的超声波信号的图。图1B是作为变换的目标的正交座标的图。另外,图1C是将处理座标和正交座标重叠了的座标变换的概念图。另外,如在图1C中重叠表示的那样,必须通过座标变换将处理座标的点表示为正交座标的点。
这些超声波图像取得装置中的座标变换和内插的处理在以前通过被称为DSC(Digital Scan Converter)的专用硬件进行(特开平11-9603号公报)。以下,说明通过DSC进行的内插处理。
在对B模式数据、或作为血流信息的血流的速度、分布以及能量的各数据(以下在不对这些数据进行区别的情况下,也有将任意一个或多个数据称为超声波数据的情况)进行内插的情况下,进行双线性(bi-linear)内插。例如是以下这样的方法:在进行核心大小(kernelsize)(核心大小是指施加用于进行内插的滤波器的大小。在此,滤波器表示所谓的平均化的程度,其平均化的程度依存于内插点的个数。)为2×2的双线性内插的情况下,将4点分为X或Y方向的2组,以通常的比例进行内插(线性内插),利用通过该内插求出的2点进而进行线性内插,由此进行原来的4点间的内插。换一种说法,是以下这样的方法:在将相邻的4点数据设为A、B、C、D的情况下,根据(1-β){(1-α)A+αB}+β{(1-α)C+αD}这样的公式计算内插值。在此,α、β是处理座标中的距离方向和方位方向(X和Y)的从内插中心的偏离的比例。
进而,根据生成B模式图像的情况说明进行核心大小2×4的内插的情况下的内插处理。图2是用于说明B模式数据的内插的图,该图被称为构造(texture)。图2A是表示处理座标中的取得的B模式数据的图。图2B是表示变换的目标的正交座标的图。图2C是放大了处理座标系中的取得的B模式数据的图。参考图2A所示的取得的B模式数据的处理座标、图2B所示的变换目标的正交座标的信息,为了从处理座标变换到正交座标,而将对应的顶点存储在暂时存储部件中。即,对应的顶点是指(0,0)与(x0,y0)的对应、(0,1)与(x0,y0)的对应、(1,0)与(x1,y1)的对应、(1,1)与(x2,y2)的对应这样的组合。在此,如图2A所示那样,用将扫描方向(v)作为纵轴,将距离(u)作为横轴的图来表示B模式数据的处理座标。另外,如图2B所示那样,用将横轴设为X座标、将纵轴设为Y座标的图,来表示变换后的正交座标。
在此,考虑抽出正交座标中的(x,y)作为内插的点的情况。说明在该情况下,计算出与正交座标系的点(x,y)对应的处理座标系的点(u,v)的方法。
接收点(u,v),根据核心大小2×4,舍去为了求出点(u,v)所需要的点以外的点,由此得到与点(u,v)对应的图2A所示的原始的采样点、A1、A2、B1、B2、C1、C2、D1、D2。如图2C所示,各采样点成为2×4的列。在此,图2是用于说明本发明的超声波图像取得装置的内插的图。如图2C所示,在设B1与C1的距离为1时,点(u,v)到将B1和B2连接起来的线为止的比例为dv,在设A1与A2的距离为1时,点(u,v)到将A1和D1连接起来的线为止的比例为du。
接着,进行A1与A2、B1与B2、C1与C2、D1与D2之间的位于比例du的位置的点的核心大小2×2的双线性内插。在此,在各内插点中,v轴方向的值分别与A1、B1、C1、D1相同,因此该2×2的内插与2点间的内插等价。即,求出(1-du)A1+duA2、(1-du)B1+duB2、(1-du)C1+duC2、(1-du)D1+duD2作为内插值。如图2C所示,假设该得到的内插值为A3、B3、C3、D3。
接着,说明根据内插值A3、B3、C3、D3以及内插函数f(x)计算出内插系数。图3是用于说明根据内插函数计算出内插系数的图。内插函数f(x)是表示从粗到密的精细度的函数。在此,假设内插函数f(x)是用于得到精细的图像的内插函数。因此,如图3所示,f(x)是在将通过上述的2×2的双线性内插计算出的4点的内插值的中间点的x座标设为0时,y座标为0的函数。图3所示的点501、点502、点503、以及点504例如表示内插点与采样点一致的情况。另外,在内插点与采样点不一致的情况下,从实际希望内插的点(u,v)到B3的比例是dv,因此在假设点(u,v)的x座标为0时,即,如点509为点(u,v)的情况那样,求出将点501移动到点505,将点502移动到点506,将点503移动到点507,将点504移动到点508,而移动了距离511后的4点的y座标的值。如图3所示,该y座标的值分别被计算为a、b、c、d,其值为内插系数。
根据所求出的内插系数,基于A1、A2、B1、B2、C1、C2、D1、D2计算出(u,v)的内插值。根据公式1求出该内插值。
公式1
a((1-du)A1+duA2)+b((1-du)B1+duB2)+c((1-du)C1+duC2)+d((1-du)D1+duD2)
du:距离方向的从内插中心的偏离
a,b,c,d:方位方向的4点的内插系数
为了通过DSC进行内插处理和座标变换的处理,必须作成专用的硬件,不但花费开发和生产的成本,在硬件的处理中只能进行开发时嵌入了的处理,在以后变更处理内容或对应新的处理的情况下,必须再次作成专用的硬件。
对于该点,近年来,也考虑了提高CPU(中央处理单元)的处理速度,利用软件通过CPU进行与DSC同等的全部处理的方法。但是,在通过CPU对帧速率为每秒60帧的图像等进行变换而显示的处理中,对CPU的负荷增大,产生显示以外的处理的响应性变慢、或丢失显示图像的问题。
另外,近年来,为了在计算机等中进行图像显示,使用GPU(图形处理单元)的情况变多。在此,GPU是指具有图像显示功能、描影(shading)功能、内插功能这样的计算功能的图像处理芯片。近年的GPU还支持座标变换功能,通过利用该功能,能够减轻座标变换对CPU的负荷。但是,由GPU支持的内插功能最大只是核心大小为2×2的双线性内插,处理能力低下。
在超声波图像取得装置中,要求高画质的图像,因此必须通过高次的核心大小进行内插,最低也需要2×4的核心大小。对于该点,在超声波图像取得装置中,在2×2的核心大小的内插和2×4的核心大小的内插中,画质差是很明显的。因此,难以只通过2×2核心大小的双线性内插的处理来作成所需要的图像。因此,在现有的超声波图像取得装置中,难以安装GPU。
进而,在对血流信号(彩色多普勒像)的速度信号进行内插的情况下,需要考虑折返(杂波:aliasing)。该折返是指表示速度信号的频谱成分超过用于取得数据的采样脉冲的循环频率的一半,即超过了Nyquist频率的情况下,将频谱成分观测为负方向的流的现象。图4A和图4B是用于说明折返的图。图4A是表示血流的流动和超声波信号的图。图4B是用于表示超声波信号和颜色的关系的图。另外,图4B是纵轴表示速度和与该速度对应的颜色的图,设中间值310为0,设最高速度(值309)为127,设最低速度(值311)为-128。例如如图4所示,考虑血液301在箭头300方向上流动的情况。通过用速度100的红色表示的点304来表示血液301的点302。另外,通过用速度-100的青色表示的点305来表示血液301的点303。在此,由于折返,将点305表示为反方向的血流。如果在此进行通常的内插,则进行点304和点305的内插,成为用速度0的黑色表示的点306。但是,实际上由于折返现象,速度从100(红色)变化到-100(青色),希望的内插值是-128(青色)的点307。难以通过处理能力低的处理器实现这样的特殊的内插。对于该点,也难以在现有的超声波图像取得装置中安装GPU。
另外,在超声波图像取得装置的图像处理中,还需要按照规定的理论对组织像(B模式像)和血流像(彩色多普勒像)进行合成而处理成为一张图像。但是,该合成理论有根据B模式像和彩色多普勒像两者的值而变化的情况,通过alpha blending这样的单纯重叠图像的方法是难以实现的。对于该点,难以将GPU安装在超声波图像取得装置中。
另外,在初始的GPU构造中,只能进行在GPU开发时嵌入的图形处理,但在近年的可编程GPU构造中,通过更新系统,能够即时地适用新开发的技术。
发明内容
本发明就是鉴于这样的问题而提出的,其目的在于:提供一种能够组合硬件和软件的双方进行内插,显示高画质的超声波图像的使用了超声波图像处理用的可编程GPU的超声波图像取得装置。
进而,其目的还在于:提供一种通过将速度、分布、能量作为复数进行处理和计算,能够容易地对血流速度信息中的折返进行处理的超声波图像取得装置。
本发明的第一形式的特征在于包括:经由超声波探头通过超声波对被检体进行扫描,取得由沿着扫描线的第一座标系中的超声波数据构成的第一超声波数据的发送接收部件;使用GPU进行计算,将上述第一超声波数据从上述第一座标系变换为用于图像显示的第二座标系的座标变换部件。该技术可以适用于超声波图像取得装置中。
本发明的第二形式的特征在于包括:经由超声波探头通过超声波对被检体进行扫描,取得由沿着扫描线的第一座标系中的超声波数据构成的第一超声波数据的发送接收部件;使用能够执行并行处理并且可编程的处理器进行计算,将上述第一超声波数据从上述第一座标系变换为用于图像显示的第二座标系的座标变换部件。该技术可以适用于超声波图像取得装置中。
根据以上形式,通过用硬件进行用于在正交座标系中显示超声波信号的内插的基础的一部分,使得进行高速的处理,进而通过程序进行一部分处理,因此通过对其程序进行编程,能够根据各种内插函数进行内插。由此,不需要开发用于进行内插处理的专用的硬件,能够削减成本,进而容易地变更处理内容和对应新的处理。
附图说明
图1A是用于说明处理座标中的超声波数据的图。
图1B是用于说明变换的目标的正交座标的图。
图1C是用于说明从处理座标变换到正交座标的图。
图2A是表示处理座标中的取得的B模式数据的图。
图2B是表示变换的目标的正交座标的图。
图2C是放大了处理座标中的取得的B模式数据的图。
图3是用于说明根据内插函数计算内插系数的图。
图4A是表示血流的流动和超声波信号的图。
图4B是用于表示超声波信号和颜色的关系的图。
图5是本发明的一个实施例的超声波诊断装置的框图。
图6是用于说明使用了复数的血流速度的内插方法的图。
图7是实施例1的超声波诊断装置的图像形成的流程图。
图8是实施例2的超声波诊断装置的图像形成的流程图。
具体实施方式
实施例1
以下,说明本发明的实施例1的超声波图像取得装置。图5是表示本发明的一个实施例的超声波图像取得装置的功能的框图。在此,图5中的B模式处理部件003、血流信息处理部件004、执行控制部件005由CPU构成。另外,图像生成部件006由GPU100构成,该GPU100具有进行用于从处理座标变换为正交座标的内插的内插功能部件110、图像显示控制功能部件120。在此,处理座标系与本发明的“第一座标系”对应。另外,正交座标系与本发明的“第二座标系”对应。该GPU100是指能够进行并行处理的可编程的处理器。另外,GPU100具有图形的描绘这样的图像显示功能、描影功能、内插功能这样的计算功能。进而,内插功能部件110具有低维内插功能部件111和高维内插功能部件112的2个功能。在此,低维内插功能部件111通过GPU100的硬件或GPU100的微代码(micro cord)进行超声波数据的内插的基础部分的计算。另外,高维内插功能部件112通过用于进行预先编程了的描影处理的程序,使用通过硬件计算出的值进行高次的内插。该描影处理一般是指图像的变形、移动以及颜色的效果的处理。另外,本发明的描影处理是指主要进行内插和座标变换的处理。在本发明的超声波图像取得装置中,实际由执行控制部件005进行全体的控制,但在以下,为了说明的方便,说明为各部件直接进行信息的交换。这样,GPU100通过硬件或微代码对简单的内插进行处理,通过程序对高度的内插进行处理。
发送接收部件002将经由超声波探头001接收到的超声波信号发送到B模式处理部件003和血流信息处理部件004。
B模式处理部件003将取得的超声波信号变换为B模式数据。在此,作为B模式数据得到的数据具有用扫描方向和距离的坐标表示的位置处的强度的信息。
接着,B模式处理部件003将得到的B模式数据发送到图像生成部件006。
血流信息处理部件004将取得的超声波信号变换为作为血流信息的血流的速度、分布、能量(以下简称为“速度、分布和能量”)的数据。在此,如背景技术所述那样,根据从检测出的超声波信号得到的复数多普勒信号,计算出作为血流信息得到的数据。
接着,血流信息处理部件004将得到的血流信息的速度、分布和能量的数据发送到图像生成部件006。
以下,将从B模式处理部件003和血流信息处理部件004发送到图像生成部件006的没有进行图像处理的数据称为原始数据。该原始数据与本发明中的“第一超声波数据”对应。另外,表示第一超声波数据的扫描方向和距离的座标是第一座标系。
执行控制部件005将由操作者输入的核心大小和内插函数发送到图像生成部件006。在此,设核心大小为2×4,内插函数为f(x)。
图像生成部件006将接收到的B模式数据、或作为血流信息的速度、分布和能量的数据存储在暂时存储部件007中。另外,根据存储在暂时存储部件007中的数据,图像生成部件006对B模式数据、或作为血流信息的速度、分布和能量的数据进行为了显示在显示部件008上的座标变换和内插。以下,分为B模式数据的情况和血流信息的情况,详细说明图像生成部件006的座标变换、内插、图像显示。
B模式数据的情况
GPU100参考图2A所示的接收到的B模式数据的处理座标处的位置信息、图2B所示的希望进行变换的正交座标的信息,为了从B模式数据位置信息变换到正交座标,将对应的顶点存储在暂时存储部件007中。在此,上述的为了从B模式数据位置信息变换为正交座标而使各顶点对应的数据是使本发明的“第一座标系”与“第二座标系”对应起来的“几何变换数据”。
GPU100抽出内插的点。在此,考虑抽出正交座标系中的(x,y)作为内插的点的情况。接着,GPU100计算出与正交座标系的点(x,y)对应的处理座标系的点(u,v)。
作为用于执行2×4的核心大小的内插的第一步骤,执行控制部件005对低维内插功能部件111进行以下的控制。
低维内插功能部件111接收点(u,v),通过舍去求解所需要的点以外的点,而得到与(u,v)对应的图2A所示的原始数据中的采样点、A1、A2、B1、B2、C1、C3、D1、D2。如图2C所示,设在B1与C1的距离为1时的(u,v)到将B1、B2连接起来的线为止的比例为dv,设在A1与A2的距离为1时的(u,v)到将A1、D1连接起来的线为止的比例为du。
低维内插功能部件111设dv=0,对A1与A2、B1与B2、C1与C2、D1与D2之间的位于比例du的位置处的点进行核心大小2×2的双线性内插。在此,对于各内插点,由于v轴方向的值分别与A1、B1、C1、D1相同,所以该2×2的内插与2点之间的内插等价。即,低维内插功能部件111计算出A3=(1-du)A1+duA2、B3 =(1-du)B1+duB2、C3=(1-du)C1+duC2、D3 =(1-du)D1+duD2作为内插值。
在本实施例的超声波图像取得装置中,由于使用GPU进行内插的计算,所以可以在求出内插点的计算中使用浮点小数。因此,可以以浮点小数的位数精度求出用于求出内插点的比例du,与现有的使用了固定小数点的内插相比,能够提高浮点小数的位数的内插精度。
低维内插功能部件111将计算出的内插值A3、B3、C3、D3发送到高维内插功能部件112。
作为用于执行2×4的核心大小的内插的步骤,执行控制部件005对高维内插功能部件112进行以下的控制。
高维内插功能部件112从执行控制部件005接收内插函数f,求出由低维内插功能部件11 1计算出的4点的内插值各自的内插系数a、b、c、d。在此,内插函数f是dv的函数,表示为(a、b、c、d)=f(dv)。作为处理内容,到此为止的处理是与现有技术一样的处理,但在本实施例中,通过用GPU100的硬件来进行该处理,能够高速地执行处理。
高维内插功能部件112根据求出的内插系数,计算出内插值。作为aA3+bB3+cC3+dD3而求出该内插值。
高维内插功能部件112将计算出的内插值发送到图像显示控制功能部件120。
图像显示控制功能部件120使用LUT(查找表:Look Up Table),根据内插值的数据,变换为所显示的颜色和灰度等级的RGB数据,在显示部件008的正交座标系的点(x,y)处显示该RGB数据。
图像生成部件006以必要的点个数进行上述内插作业,将基于取得的B模式数据的B模式图像显示在显示部件008上。
血流信息的情况
GPU100利用LUT,根据作为血流信息的速度(V)、分布(T)、能量(P),使用下式,求出具有实数项Re、虚数项Im的复数数据。
Re=((255-T)/2)cos(πV/128)
Im=((255-T)/2)sin(πV/128)
在此,由于使用8比特的信号进行处理,所以系数使用上述公式中的系数的值,但该系数根据所使用的信号而变化,并不只限于该值。
GPU100通过内插功能部件110,对由该Re、Im以及P的3个值构成的数据,进行用于从用X轴表示Re、Y轴表示Im的处理座标变换为正交座标的内插,求出内插值。该内插方法是与在B模式数据的情况下说明的方法一样的方法,分别对Re、Im、P进行。
如果取得与上述内插了的Re和Im的相同位置的值,考虑复数Re+jIm,则成为对其偏角进行了内插的速度。如图6所示,折返速度附近的2个向量a和向量b的内插向量为向量c,得到与图4所示的内插大致同等的结果。在此,图6是用于说明使用了复数的血流速度的内插方法的图。进而,如果参考上述Re和Im,则在该内插中,用1-T=|C1|/P这样的函数进行内插。另外,由于如背景技术所说明的那样,表示为P=∑|xi|2、C1=∑xi*xi+1,所以如果C1变大,则P也变大。即,由于以某程度的大小抑制1-T,所以能够抑制动态范围,能够容易地进行内插处理。
另外,GPU100能够通过内插功能部件110进行帧的内插。帧的内插是指:为了使帧速率变慢,平滑地显示彩色多普勒像的时间过程,实际在扫描生成的彩色多普勒像之间插入内插帧进行显示。该内插处理也利用Re、Im、P的值进行。通过2帧之间的线性内插进行该内插处理。进而,内插功能部件110根据求出的Re、Im的内插值,求出速度(V)、分布(T)、能量(P)的内插值。
接着,GPU100进行原始数据和内插了的数据的着色处理。该处理是使求出的数据可视化的处理。该处理根据在彩色多普勒像中显示的信息(以下称为“显示模式”)是速度(V)、分布(T)、能量(P)中的哪一个信息而变化。
在彩色多普勒像的显示模式是速度(V)和分布(T)的情况下,GPU100根据以下的公式从Re、Im计算出V、T。
V=(128/π)atan2(Im,Re)
T=255-2√(Re2+Im2)
P=(1-|C1|)/T
在此,atan2是在-π~π的范围内求出角度的arctangent函数。
通过以上那样的内插处理得到的内插速度在基于现有的折返对应内插得到的速度和图像上近似到难以识别的程度,因此可以说该内插处理具有充分的精度。
接着,GPU100通过图像显示控制功能部件120,使用LUT将所得到的速度(V)和分布(T)变换为颜色和灰度等级的RGB,并显示在显示部件008上。
在能量(P)的情况下,根据原始的能量值或内插了的能量(P)值而决定彩色多普勒像的显示模式。另外,GPU100通过图像显示控制功能部件120,使用LUT,将能量(P)变换为颜色和灰度等级的RGB,并显示在显示部件008上。
如上所述,在作成彩色多普勒像的情况下,内插处理需要进行速度(V)、分布(T)、以及能量(P)这样的3个数据的内插处理。然后,GPU100的图像显示控制功能部件120具备:在图像显示的处理中,能够并行地对红、青、绿、阿尔法(alpha)这样的4个数据进行处理的并行电路。因此,也可以构成为在对红、青、绿的颜色进行处理的3个并行处理电路中分别对速度(V)、分布(T)、能量(P)这样的3个数据进行处理,由此并行地进行内插处理。
在此,本发明中的“座标变换部件”包含内插功能部件110,进行以上说明了的GPU100的内插和座标变换。另外,座标变换部件是图像形成部件006的构成要素之一。
接着,参考图7说明本实施例的超声波图像的形成流程。在此,图7是本实施例的超声波图像的形成的流程图。
步骤S001:操作者输入核心大小和内插函数。
步骤S002:发送接收部件002将经由超声波探头001接收到的超声波信号发送到B模式处理部件003和血流信息处理部件004。
B模式处理
步骤S003:低维内插功能部件111将发送来的B模式数据存储在暂时存储部件007中,参考输入的核心大小,基于暂时存储部件007的B模式数据进行2×2的内插,并将该内插值发送到高维内插功能部件112。
步骤S004:高维内插功能部件112参考内插函数计算出内插系数,使用该内插系数进行比2×2更高维的内插,内插部件110将该内插值和原样使用的原始数据发送到图像显示控制功能部件120。
血流信息处理
步骤S005:GPU100将作为发送来的血流信息的速度(V)、分布(T)、能量(P)的数据存储在暂时存储部件007中,基于该数据求出Re、Im。
步骤S006:低维内插功能部件111参考输入的核心大小,基于求出的Re、Im进行2×2的内插,并将该内插值发送到高维内插功能部件112。
步骤S007:高维内插功能部件112参考内插函数计算出内插系数,使用该内插系数,进行比2×2更高维的内插,例如2×4的内插。
步骤S008:内插功能部件110进行帧内插。
步骤S009:内插功能部件110基于Re、Im的内插值,求出速度(V)、分布(T)、能量(P)的内插值,将该内插值的数据和原样使用的原始数据发送到图像显示控制功能部件120。
(B模式处理和血流信息处理)
步骤S010:图像显示控制功能部件120使用LUT,将发送来的数据变换为RGB。
步骤S011:图像显示控制功能部件120基于变换了的RGB的值,将图像显示在显示部件008上。
如上所述,在本发明的超声波图像取得装置中,可以使用通常的GPU,使用低维内插功能部件111,通过硬件或GPU的微代码,高速地进行核心大小直到2×2为止的内插,基于该内插值,使用高维内插功能部件112通过程序进行核心大小2×4的内插。由此,在能够削减成本的同时,能够自由地对内插处理进行编程,因此容易地更新处理内容和对应新的处理。另外,由于能够使用通用的GPU,所以能够抑制超声波图像取得装置的制造费用,削减成本。进而,还能够降低超声波图像取得装置的开发费用。另外,由于通过GPU进行内插计算,所以能够基于GPU的处理能力的浮点小数计算执行内插处理,由此,能够提高与GPU的处理能力对应的浮点小数位数的内插精度,能够提高所生成的超声波断层像的画质。另外,由于能够通过GPU的几何变换执行座标变换,所以能够进行从极座标到正交座标的变换所不能进行的座标变换(无法用公式表示的座标变换)。因此,能够容易地高精度地进行需要无法通过从极座标向正交座标的变换来进行的座标变换的特殊检测和使用了2维阵列探头的扫描。进而,通过将血流的速度、分布和能量的处理分别分配给作为并行处理系统而由GPU所具有的R、G、B的处理系统等,能够同时进行计算,可以不必维持3个系统的相同电路而采用单纯的结构,还能够抑制成本。进而,还能够通过并行处理而提高处理速度。
实施例2
以下,说明本发明的实施例2的超声波图像取得装置。本实施例的超声波图像取得装置是在实施例1中,进行所生成的B模式图像与彩色多普勒像的合成处理,而重叠进行显示。因此,说明B模式图像和彩色多普勒像的合成处理。在此,本实施例的超声波图像取得装置也具有图5所示的功能模块。另外,由GPU100进行合成处理。
GPU100的图像显示控制功能部件120具有红、青、绿、阿尔法(alpha)这样的4个并行处理电路。另外,作为通常的处理电路的使用方法,GPU100通过前面的3个处理电路对颜色进行处理,通过阿尔法的处理电路对不透明度进行处理。在本发明中,在GPU100中,使用处理红、青、绿的处理电路作为在作成B模式图像和作成彩色多普勒像时使用的处理电路,因此不使用处理阿尔法数据的电路。因此,在本实施例的超声波图像取得装置中,预先将符合图像条件的阈值信息,存储在GPU100中的图像显示控制功能部件120的上述红、青、绿、阿尔法这样的4个并行处理电路中的处理阿尔法数据的电路中。
该阈值信息是指用于进行以下处理等的一定值,即在B模式数据为一定值以下或一定值以上时不显示B模式图像的处理、在血流信息的值为某一定值以下或一定值以上时不显示彩色多普勒像的处理。这是例如在生成心脏的超声波图像的情况下,在B模式数据的值的亮度高的情况下,此处应该有心脏的壁而没有血流,因此将该处的部分的血流信息认为是噪声。因此,必须消去该处的部分的彩色多普勒图像。另外,在生成腹部的超声波图像的情况下,血流的信息是重要的,如果血流为某一定以上,则需要无条件地显示血流。因此,如果血流为某一定以下,则必须作为噪声消去。进而,在能量值是低的值的情况下,是噪声的可能性高,因此,必须不显示彩色多普勒像。这样进行重叠时,需要根据各种条件进行合成处理,因此必须存储该图像条件和阈值信息。
GPU100的图像显示控制功能部件120从高维内插功能部件112接收B模式数据和作为血流信息的速度(V)、分布(T)和能量(P)的数据,参考存储在对阿尔法数据进行处理的电路中的图像条件和阈值信息,决定在各点显示哪个数据。
图像显示控制功能部件120基于接收到的B模式数据和作为血流信息的速度(V)、分布(T)和能量(P)的数据、以及在各点显示哪个数据的决定,使用LUT,变换为所显示的颜色和灰度等级的RGB数据,将对B模式图像和彩色多普勒图像进行了合成处理后的图像显示在显示部件008上。
以上,说明了优先显示B模式图像和彩色多普勒像的任意一个的情况,但也可以构成为在一定的条件下显示双方。例如在根据B模式数据和血流信息中的一个数据(在此,考虑速度(V))而显示为与速度(V)的数据重叠地半透明地可以看到B模式数据时,图像显示控制功能部件120将B模式数据和速度(V)的数据的双方发送到阿尔法数据的处理电路,在计算出将B模式数据设定为怎样的透明度并进行了合成处理后,使用LUT变换为RGB数据,并显示在显示部件008上。
接着,参考图8说明本实施例的超声波图像的形成流程。在此,图8是本实施例的超声波图像的形成的流程图。
步骤S101:操作者输入核心大小和内插函数。
步骤S102:发送接收部件002将经由超声波探头001接收到的超声波信号发送到B模式处理部件003和血流信息处理部件004。
B模式处理
步骤S103:低维内插功能部件111将发送来的B模式数据存储在暂时存储部件007中,参考输入的核心大小,基于暂时存储部件007的B模式数据进行2×2的内插,并将该内插值发送到高维内插功能部件112。
步骤S104:高次内插功能部件112参考内插函数计算出内插系数,使用该内插系数进行比2×2更高维的内插,内插部件110将该内插值和原样使用的原始数据发送到图像显示控制功能部件120。
血流信息处理
步骤S105:GPU100将发送来的作为血流信息的速度(V)、分布(T)、能量(P)的数据存储在暂时存储部件007中,基于该数据求出Re、Im。
步骤S106:低维内插部件111参考输入的核心大小,基于求出的Re、Im进行2×2的内插,并将该内插值发送到高维内插功能部件112。
步骤S107:高维内插功能部件112参考内插函数计算出内插系数,使用该内插系数,进行比2×2更高维的内插。
步骤S108:内插功能部件110进行帧内插。
步骤S109:内插功能部件110基于Re、Im的内插值,求出速度(V)、分布(T)、能量(P)的内插值,将该内插值的数据和原样使用的原始数据发送到图像显示控制功能部件120。
(B模式处理和血流信息处理)
步骤S110:图像显示控制功能部件120基于存储在阿尔法数据的处理电路中的图像条件和阈值信息,决定显示的数据。
步骤S111:图像显示控制功能部件120使用LUT,将决定了的数据变换为RGB。
步骤S112:图像显示控制功能部件120基于变换了的RGB的值,将图像显示在显示部件008上。
如上所述,在将B模式图像和彩色多普勒像重叠而进行合成处理的情况下,通过使用GPU的处理阿尔法数据的电路,还能够与图像的作成并行地进行合成处理,因此能够高速地处理合成图像那样的复杂图像并进行显示,能够对超声波诊断的效果提高起作用。

Claims (7)

1.一种超声波图像取得装置,其特征在于包括:
经由超声波探头用超声波对被检体进行扫描,取得由沿着扫描线的第一座标系中的超声波数据构成的第一超声波数据的发送接收部件;和
使用GPU进行计算,将上述第一超声波数据从上述第一座标系变换为用于图像显示的第二座标系的座标变换部件,
上述座标变换部件根据上述第一座标系中的多个点,利用上述GPU进行用于内插上述第二座标系中的1点的内插计算,
在上述内插计算中,接受施加用于进行内插的滤波器的大小即核心大小的输入,通过GPU的硬件或GPU的微代码,进行2×2的双线性内插,求出内插值,通过GPU的程序,基于该内插值和规定的内插函数,求出输入的核心大小的内插值。
2.根据权利要求1所述的超声波图像取得装置,其特征在于:
上述座标变换部件取得使上述第一座标系和上述第二座标系对应起来的几何变换数据,上述GPU基于上述几何变换数据进行几何变换,而进行上述座标变换。
3.根据权利要求1所述的超声波图像取得装置,其特征在于:
上述第一超声波数据包含多种血流信息,
上述座标变换部件同时对上述多种血流信息中的多个座标变换进行处理。
4.根据权利要求1所述的超声波图像取得装置,其特征在于:
上述座标变换部件求出以包含在上述第一超声波数据中的血流的速度、分布和能量为要素的复数,在计算出该复数数据的内插值后,基于该内插值分别变换为上述血流的速度、分布和能量,由此求出上述血流的速度、分布和能量的内插值。
5.根据权利要求3所述的超声波图像取得装置,其特征在于:
上述多个血流信息包含与血流的速度、分布和能量有关的信息,
上述座标变换部件使速度、分布和能量的各个座标变换和内插处理与GPU所具有的用于处理红、青或绿的各个颜色的并行处理电路对应地同时进行处理。
6.一种超声波图像取得装置,其特征在于包括:
经由超声波探头用超声波对被检体进行扫描,取得由沿着扫描线的第一座标系中的超声波数据构成的第一超声波数据的发送接收部件;和
使用能够执行并行处理并且可编程的处理器进行计算,将上述第一超声波数据从上述第一座标系变换为用于图像显示的第二座标系的座标变换部件,
上述座标变换部件根据上述第一座标系中的多个点,利用上述处理器进行用于内插上述第二座标系中的1点的内插计算,
在上述内插计算中,接受施加用于进行内插的滤波器的大小即核心大小的输入,通过上述处理器的硬件或上述处理器的微代码,进行2×2的双线性内插,求出内插值,通过上述处理器的程序,基于该内插值和规定的内插函数,求出输入的核心大小的内插值。
7.根据权利要求6所述的超声波图像取得装置,其特征在于:
上述处理器代替CPU进行图形的描绘。
CN200710167811.0A 2006-12-27 2007-10-26 超声波图像取得装置 Active CN101209210B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2006351833 2006-12-27
JP2006351833A JP5022700B2 (ja) 2006-12-27 2006-12-27 超音波診断装置
JP2006-351833 2006-12-27

Publications (2)

Publication Number Publication Date
CN101209210A CN101209210A (zh) 2008-07-02
CN101209210B true CN101209210B (zh) 2010-07-07

Family

ID=39584985

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200710167811.0A Active CN101209210B (zh) 2006-12-27 2007-10-26 超声波图像取得装置

Country Status (3)

Country Link
US (1) US8708911B2 (zh)
JP (1) JP5022700B2 (zh)
CN (1) CN101209210B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2143384A1 (en) * 2008-07-09 2010-01-13 Medison Co., Ltd. Enhanced ultrasound data processing in an ultrasound system
EP2344043B1 (en) * 2008-08-01 2014-05-07 Esaote Europe B.V. Portable ultrasound system
JP2010103642A (ja) * 2008-10-21 2010-05-06 Toshiba Corp シェーディング補正装置
US20100135580A1 (en) * 2008-11-28 2010-06-03 Inventec Appliances Corp. Method for adjusting video frame
EP2199971A1 (en) 2008-12-15 2010-06-23 Medison Co., Ltd. Image processing system and method
JP5449852B2 (ja) * 2009-05-08 2014-03-19 株式会社東芝 超音波診断装置
WO2012063977A1 (ja) * 2010-11-11 2012-05-18 オリンパスメディカルシステムズ株式会社 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5936857B2 (ja) * 2011-01-27 2016-06-22 株式会社東芝 超音波診断装置及びその画像処理プログラム
KR101310932B1 (ko) 2012-01-11 2013-09-25 알피니언메디칼시스템 주식회사 방향성의 왜곡없이 실시간 보간이 가능한 주사 변환기, 초음파 진단장치 및 방법
CN102631222B (zh) * 2012-04-26 2014-01-01 珠海医凯电子科技有限公司 基于gpu实现的超声成像扫描变换方法
CN103767737B (zh) * 2012-10-23 2015-09-23 深圳市蓝韵实业有限公司 提高彩色血流检测灵敏度的方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5682896A (en) * 1996-03-28 1997-11-04 Diasonics Ultrasound, Inc. Method and apparatus for generating volume flow measurement
US6685641B2 (en) * 2002-02-01 2004-02-03 Siemens Medical Solutions Usa, Inc. Plane wave scanning reception and receiver

Family Cites Families (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2674157A (en) * 1949-08-15 1954-04-06 Leitz Ernst Gmbh Phase microscope
DE3442218A1 (de) * 1984-11-19 1986-05-28 Fa. Carl Zeiss, 7920 Heidenheim Auflichtbeleuchtungsapparat fuer mikroskope
US4974094A (en) * 1989-12-04 1990-11-27 Yuhkoh Morito Direct lighting/illuminating system for miniature CCD camera
US5429137A (en) * 1994-06-03 1995-07-04 Siemens Medical Systems, Inc. Acoustic scan conversion method and apparatus for velocity flow
US5820250A (en) * 1995-10-24 1998-10-13 Dolan-Jenner Industries, Inc. Dark field illuminator ringlight adaptor
US5720291A (en) * 1996-03-22 1998-02-24 Advanced Technology Laboratories, Inc. Three dimensional medical ultrasonic diagnostic image of tissue texture and vasculature
US5860924A (en) * 1996-11-26 1999-01-19 Advanced Technology Laboratories, Inc. Three dimensional ultrasonic diagnostic image rendering from tissue and flow images
JP3746122B2 (ja) * 1996-12-06 2006-02-15 株式会社東芝 2次元ドプラ形超音波診断装置
JP3806229B2 (ja) * 1997-06-27 2006-08-09 アロカ株式会社 超音波診断装置
JP2000175914A (ja) * 1998-12-17 2000-06-27 Fukuda Denshi Co Ltd 超音波診断装置
US7096079B2 (en) * 1999-10-14 2006-08-22 Sony Computer Entertainment Inc. Audio processing and image generating apparatus, audio processing and image generating method, recording medium and program
WO2001041891A1 (fr) * 1999-12-10 2001-06-14 Namco, Ltd. Machine de jeu
US20010033287A1 (en) * 2000-01-11 2001-10-25 Sun Microsystems, Inc. Graphics system having a super-sampled sample buffer which utilizes a window ID to specify pixel characteristics
JP2001224850A (ja) * 2000-02-16 2001-08-21 Konami Co Ltd ゲーム装置、ゲーム装置の制御方法、情報記憶媒体、ゲーム配信装置及びゲーム配信方法
JP2002014672A (ja) * 2000-06-13 2002-01-18 Doogi Doogi Drm Co Ltd ドラム教育兼用娯楽装置
US6541692B2 (en) * 2000-07-07 2003-04-01 Allan Miller Dynamically adjustable network enabled method for playing along with music
JP2002066128A (ja) * 2000-08-31 2002-03-05 Konami Co Ltd ゲーム装置、ゲーム処理方法及び情報記憶媒体
JP2002078970A (ja) * 2000-09-08 2002-03-19 Alps Electric Co Ltd ゲーム用入力装置
JP2002258842A (ja) * 2000-12-27 2002-09-11 Sony Computer Entertainment Inc 音声制御装置、音声制御方法、音声制御プログラム、音声制御プログラムを記憶したコンピュータ読み取り可能な記憶媒体、音声制御プログラムを実行するプログラム実行装置
US7191023B2 (en) * 2001-01-08 2007-03-13 Cybermusicmix.Com, Inc. Method and apparatus for sound and music mixing on a network
US6898637B2 (en) * 2001-01-10 2005-05-24 Agere Systems, Inc. Distributed audio collaboration method and apparatus
JP3685731B2 (ja) * 2001-03-28 2005-08-24 任天堂株式会社 ゲーム装置及びそのプログラム
JP3576994B2 (ja) * 2001-04-27 2004-10-13 株式会社コナミコンピュータエンタテインメントスタジオ ゲームサーバ、ネットゲーム進行制御プログラム及びネットゲーム進行制御方法
US6679843B2 (en) * 2002-06-25 2004-01-20 Siemens Medical Solutions Usa , Inc. Adaptive ultrasound image fusion
US20040138560A1 (en) * 2002-12-02 2004-07-15 Gianluca Paladini Real-time scan conversion and rendering of ultrasound data
US7037263B2 (en) * 2003-08-20 2006-05-02 Siemens Medical Solutions Usa, Inc. Computing spatial derivatives for medical diagnostic imaging methods and systems
US7119810B2 (en) * 2003-12-05 2006-10-10 Siemens Medical Solutions Usa, Inc. Graphics processing unit for simulation or medical diagnostic imaging
US8021300B2 (en) * 2004-06-16 2011-09-20 Siemens Medical Solutions Usa, Inc. Three-dimensional fly-through systems and methods using ultrasound data
JP3943563B2 (ja) * 2004-08-11 2007-07-11 ザイオソフト株式会社 画像表示方法及び画像表示プログラム
US7744538B2 (en) * 2004-11-01 2010-06-29 Siemens Medical Solutions Usa, Inc. Minimum arc velocity interpolation for three-dimensional ultrasound imaging
US20060114254A1 (en) * 2004-11-26 2006-06-01 Tim Day Volume rendering apparatus and method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5682896A (en) * 1996-03-28 1997-11-04 Diasonics Ultrasound, Inc. Method and apparatus for generating volume flow measurement
US6685641B2 (en) * 2002-02-01 2004-02-03 Siemens Medical Solutions Usa, Inc. Plane wave scanning reception and receiver

Also Published As

Publication number Publication date
JP5022700B2 (ja) 2012-09-12
US20080161690A1 (en) 2008-07-03
CN101209210A (zh) 2008-07-02
US8708911B2 (en) 2014-04-29
JP2008161289A (ja) 2008-07-17

Similar Documents

Publication Publication Date Title
CN101209210B (zh) 超声波图像取得装置
CN100518656C (zh) 超声波摄像装置及其摄像数据处理方法以及处理程序
US6413219B1 (en) Three-dimensional ultrasound data display using multiple cut planes
JP4297555B2 (ja) 超音波カラー・ドップラー速度/方向撮像
JP2934402B2 (ja) 3次元超音波画像作成方法および画像処理装置
US8187192B2 (en) Method and apparatus for scan conversion and interpolation of ultrasonic linear array steering imaging
US20050228280A1 (en) Acquisition and display methods and systems for three-dimensional ultrasound imaging
US10675006B2 (en) Registration for multi-modality medical imaging fusion with narrow field of view
CN103117010B (zh) 一种超声成像模拟系统
JP5525930B2 (ja) 三次元超音波画像を生成表示する超音波診断装置
CN102487603A (zh) 超声波诊断装置
US8070683B2 (en) Scan conversion for ultrasonic imaging and apparatus using the same
WO2013153857A1 (ja) 超音波診断装置及び軌跡表示方法
JP7262298B2 (ja) 水中探知装置、および、水中探知方法
US11364012B2 (en) 3-D imaging via free-hand scanning with a multiplane US transducer
CN108022209A (zh) 颜色值的获取方法及装置、空间复合成像的方法及设备
JPH0428354A (ja) 超音波診断装置
US6458082B1 (en) System and method for the display of ultrasound data
CN109580786B (zh) 一种超声探头校准方法
US6306092B1 (en) Method and apparatus for calibrating rotational offsets in ultrasound transducer scans
US8911373B2 (en) Vector flow ultrasound imaging
US8002702B2 (en) Ultrasonograph and method for controlling movement of display body of ultrasonograph
JPH02211134A (ja) カラー超音波診断装置
CN102958446B (zh) 超声波诊断装置以及医用图像处理装置
JP4703144B2 (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: 20160720

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.