CN101822547A - 超声波诊断装置及其控制方法、图像处理装置和方法 - Google Patents
超声波诊断装置及其控制方法、图像处理装置和方法 Download PDFInfo
- Publication number
- CN101822547A CN101822547A CN201010128835A CN201010128835A CN101822547A CN 101822547 A CN101822547 A CN 101822547A CN 201010128835 A CN201010128835 A CN 201010128835A CN 201010128835 A CN201010128835 A CN 201010128835A CN 101822547 A CN101822547 A CN 101822547A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- volume data
- voxel
- high frequency
- low frequency
- 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
- 238000002604 ultrasonography Methods 0.000 title claims abstract description 177
- 238000000034 method Methods 0.000 title claims abstract description 57
- 239000000523 sample Substances 0.000 claims abstract description 50
- 238000001914 filtration Methods 0.000 claims abstract description 28
- 238000003672 processing method Methods 0.000 claims abstract description 17
- 238000009792 diffusion process Methods 0.000 claims description 129
- 230000009466 transformation Effects 0.000 claims description 71
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000005540 biological transmission Effects 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 description 23
- 241001269238 Data Species 0.000 description 19
- 150000001875 compounds Chemical class 0.000 description 13
- 230000004069 differentiation Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000009206 nuclear medicine Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000003696 structure analysis method Methods 0.000 description 1
- 238000010189 synthetic method Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/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/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- 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/52077—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 with means for elimination of unwanted signals, e.g. noise or interference
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
- G06T2207/10136—3D ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Gynecology & Obstetrics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明提供一种超声波诊断装置、图像处理装置、超声波诊断装置的控制方法以及图像处理方法。发送接收部经由超声波探头向被检体发送超声波,接收由上述被检体反射的超声波,输出与上述接收到的超声波对应的回波信号。B模式体数据产生部或多普勒体数据产生部根据回波信号产生第1体数据。散斑/噪声减少处理部对第1体数据中包含的各体素,按照与体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据。三维图像处理部对第2体数据进行三维图像处理,产生二维的超声波图像的数据。
Description
相关申请的交叉引用
本申请基于2009年3月4日提交的在先的日本专利申请No.2009-051145并要求其为优先权,其日本专利申请的全部内容在本申请中被援引。
技术领域
本发明涉及减少在体数据中包含的散斑、噪声的超声波诊断装置、图像处理装置、超声波诊断装置的控制方法以及图像处理方法。
背景技术
超声波诊断装置向被检体内放射由内置在超声波探头内的振子所产生的超声波脉冲,并经由振子接收由被检体组织反射出的超声波。并且,超声波诊断装置产生与接收到的超声波相应的回波信号,根据产生的回波产生超声波图像数据,并显示超声波图像。一般而言,超声波诊断装置通过使用具有排列成一维的多个振子的超声波探头,扫描被检体内的断面来收集二维超声波图像。另外,近年来,超声波诊断装置通过使用具有排列成二维的多个振子的超声波探头等,也能够扫描被检体内的体积,收集三维超声波图像(体数据)。
但是,由多个接近的被检体组织反射的超声波因各自的相位而发生干涉。由此干涉,产生与只合成振幅的情况所能看到的内容不同的图像图案即散斑(speckle)。散斑会妨碍被检体组织的边界的位置、形状的正确观测。因此,提出了用于减少(reduce)散斑的各种处理方法。
例如,如第1文献(日本特开2006-116307号公报)所记载的那样,提出了通过小波(wavelet)变换/逆变换等对超声波图像进行多分辨率分解,并检测出在各等级下分解出的图像的边缘,对每个象素算出边缘的朝向,实施在边缘的切线方向平滑化、在法线方向尖锐化的滤波的方法。但是,该第1文献的应用被限定于二维超声波图像。
另一方面,如第2文献(K.Z.Abd-Elmoniem,A.M.Youssef,andY.M.Kadah,″Real-Time Speckle Reduction and CoherenceEnhancement in Ultrasound Imaging via Nonlinear AnisotropicDiffusion″,IEEE transactions on biomedical engineering,vol.49,NO.9,Sep.2002)所记载的那样,作为图像数据中的象素区域的结构分析方法,有使用结构张量(Structure Tensor)的方法。第2文献中记载了基于非线性各向异性扩散滤波器(Nonlinear Anisotropic DiffusionFilter)减少超声波图像的散斑的应用。但是第2文献的应用也被限定于二维超声波图像。
但是,通过在与特定的坐标轴垂直的断面上分割体数据,并在那些各断面上应用二维散斑减少滤波器,可以减少体数据内的散斑。但是,此时与在应用了二维滤波器的面中的散斑的减少精度相比,在与其垂直相交的2个面中的散斑的减少精度差。
三维地扩充非线性各向异性扩散滤波器时,需要根据三维结构张量的3个固有值的大小关系对结构的特征进行分类。在第3文献(Z.Yu,C.Bajaj,″A Structure Tensor Approach for 3D Image Skeletonization:Applications in Protein Secondary Structure Analysis,Image Processing,2006IEEE International Conference on,2006)中记载了其具体例子。但是,在该第3文献中没有记载在三维扩散方程式中结合三维结构张量的各固有值的方法。
发明内容
本发明的目的在于:提供一种实现对体数据中包含的散斑、噪声的减少精度以及减少处理速度进行提高的超声波诊断装置、图像处理装置、超声波诊断装置的控制方法以及图像处理方法。
根据本发明的第1实施方式提供一种超声波诊断装置,其特征在于,包括:超声波探头;发送接收部,经由上述超声波探头向被检体发送超声波,接收由来于由上述被检体反射的超声波的回波信号;第1体数据产生部,根据上述接收到的回波信号产生第1体数据;低频/高频体数据产生部,根据上述产生的第1体数据,产生在空间频率上的低频体数据与高频体数据;滤波器部,计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算边缘信息,对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器;高频体素控制部,根据上述产生的边缘信息,对上述高频体数据实施边缘强调处理;以及第2体数据产生部,根据在上述滤波器部中应用了非线性各向异性扩散滤波器后的低频体数据与在上述高频体素控制部中实施了边缘强调处理后的高频体数据,产生散斑和噪声被减少了的第2体数据。
根据本发明的第2实施方式提供一种超声波诊断装置,其特征在,包括:超声波探头;发送接收部,经由上述超声波探头向被检体发送超声波,接收由上述被检体反射的超声波,输出与上述接收到的超声波对应的回波信号;第1体数据产生部,根据上述回波信号产生第1体数据;第2体数据产生部,对上述第1体数据中包含的各个体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据;以及图像数据产生部,根据上述第2体数据产生二维的超声波图像的数据。
根据本发明的第3实施方式提供一种超声波诊断装置,其特征在于,包括:超声波探头;发送接收部,经由上述超声波探头向被检体发送超声波,接收由上述被检体反射的超声波,输出与上述接收的超声波对应的回波信号;产生部,根据上述回波信号产生第1体数据;减少部,通过减少上述第1体数据中包含的具有三维各向同性结构的分量,强调具有三维各向异性结构的分量,从而根据上述第1体数据产生第2体数据;以及三维图像处理部,对上述第2体数据进行三维图像处理,产生二维的超声波图像的数据。
根据本发明的第4实施方式提供一种图像处理装置,其特征在于,包括:存储部,存储与被检体相关的第1体数据;第2体数据产生部,对上述第1体数据中包含的各体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据;图像数据产生部,根据上述第2体数据产生二维的超声波图像的数据。
根据本发明的第5实施提供一种图像处理装置,其特征在于,包括:存储部,存储与被检体相关的第1体数据;减少部,通过减少上述第1体数据中包含的具有三维各向同性结构的分量,强调具有三维各向异性结构的分量,从而根据上述第1体数据产生第2体数据;三维图像处理部,对上述第2体数据进行三维图像处理,产生二维图像的数据。
根据本发明的第6实施方式提供一种超声波诊断装置的控制方法,该超声波诊断装置具备:超声波探头、经由上述超声波探头发送接收超声波的发送接收部、产生第1体数据的第1体数据产生部、产生低频体数据与高频体数据的低频/高频体数据产生部、应用滤波器的滤波器部、控制高频体素的高频体素控制部、产生第2体数据的第2体数据产生部、***控制部,该超声波诊断装置的控制方法的特征在于:上述***控制部控制上述发送接收部,以使得经由上述超声波探头向被检体发送超声波,经由上述超声波探头接收由上述被检体反射的超声波,输出与上述接收的超声波对应的回波信号,上述***控制部控制上述第1体数据产生部,以使得根据上述输出的回波信号产生第1体数据,上述***控制部控制上述低频/高频体数据产生部,以使得根据上述产生的第1体数据产生在空间频率上的低频体数据与高频体数据,上述***控制部控制上述滤波器部,以使得计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算边缘信息,对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器,上述***控制部控制上述高频体素控制部,以使得根据上述产生的边缘信息对上述高频体数据实施边缘强调处理,上述***控制部控制上述第2体数据产生部,以使得根据应用了上述非线性各向异性扩散滤波器后的低频体数据与实施了上述边缘强调处理后的高频体数据,产生散斑与噪声被减少的第2体数据。
根据本发明的第7实施方式提供一种超声波诊断装置的控制方法,该超声波诊断装置具备:超声波探头、经由上超声波探头发送接收超声波的发送接收部、产生第1体数据的第1体数据产生部、产生第2体数据的第2体数据产生部、实施三维图像处理的三维图像处理部、***控制部,该超声波诊断装置的控制方法的特征在于:上述***控制部控制上述接收部,以使得经由上述超声波探头向被检体发送超声波,经由上述超声波探头接收由上述被检体反射的超声波,输出与上述接收的超声波对应的回波信号,上述***控制部控制上述第1体数据产生部,以使得根据上述输出的回波信号产生第1体数据,上述***控制部控制上述第2体数据产生部,以使得通过对上述产生的第1体数据中包含的各个体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波产生第2体数据,上述***控制部控制上述三维图像处理部,以使得根据上述产生的第2体数据产生二维超声波图像的数据。
根据本发明的第8实施方式提供一种图像处理方法,其特征在于:根据与被检体相关的第1体数据,产生在空间频率上的低频体数据与高频体数据,计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算边缘信息,对上述低频体数据应用了具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器,根据上述产生的边缘信息,对上述高频体数据实施边缘强调处理,根据应用了上述非线性各向异性扩散滤波器后的低频体数据与实施了上述边缘强调处理的高频体数据,产生散斑与噪声被减少了的第2体数据。
根据本发明的第9实施方式提供一种图像处理方法,其特征在于:通过对与被检体相关的第1体数据中包含的各体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据,根据上述第2体数据产生二维图像的数据。
在下面的描述中将提出本发明的其它目的和优点,部分内容可以从说明书的描述中变得明显,或者通过实施本发明可以明确上述内容。通过下文中详细指出的手段和组合可以实现和得到本发明的目的和优点。
附图说明
结合在这里并构成说明书的一部分的附图描述本发明当前优选的实施方式,并且与上述的概要说明以及下面的对优选实施方式的详细描述一同用来说明本发明的原理。
图1为表示与本发明的实施方式相关的超声波诊断装置的结构的图。
图2为表示图1的散斑/噪声减少处理部的结构的图。
图3为表示由图2的三维小波变换部进行三维小波变换的概念的图。
图4为表示在图3的小波变换中的各坐标轴方向应用的滤波器的应用模式的图。
图5为表示根据三维结构张量的3个固有值的大小关系分类的3个结构模式之一的图。
图6为表示根据三维结构张量的3个固有值的大小关系分类的3个结构模式之一的图。
图7为表示根据三维结构张量的3个固有值的大小关系分类的3个结构模式之一的图。
图8为表示图2的非线性各向异性扩散滤波器部进行处理的典型流程的图。
图9为表示与本实施方式的变形例1相关的超声波诊断装置结构的图。
具体实施方式
以下,参照附图说明与本发明的实施方式相关的超声波诊断装置、图像处理装置、超声波诊断装置的控制方法以及图像处理方法。
图1为表示与本实施方式相关的超声波诊断装置的结构的图。如图1所示,本超声波诊断装置具备超声波探头1、发送接收部2、B模式处理部3、B模式体数据产生部4、多普勒处理部5、多普勒体数据产生部6、散斑/噪声减少处理部7、三维图像处理部8以及显示部9。
超声波探头1具有排列成二维状的多个振子。超声波探头1接收来自发送接收部2的驱动信号并向被检体发送超声波。向被检体发送的超声波在体内组织的声阻抗的不连续面上依次被反射。被反射的超声波作为回波信号由超声波探头1接收。该回波信号的振幅取决于被反射的不连续面中的声阻抗的差异。另外,被发送的超声波在移动的血流或心脏壁等移动体的表面被反射时,回波信号因多普勒效应而受到取决于移动体的超声波发送方向的速度分量的频率偏移。另外,超声波探头1如果能够三维扫描,不是必须是二维阵列型。例如,超声波探头1也可以是能够机械地摇动的一维阵列型。
发送接收部2经由超声波探头1使用超声波重复三维扫描作为被检体的扫描对象的体积。通过该三维扫描由发送接收部2输出与扫描体有关的关于多个扫描线的多个回波信号。
更详细而言,发送接收部2具有超声波发送用的未图示的速率脉冲产生电路、发送延迟电路以及驱动脉冲产生电路等。速率脉冲产生电路以规定的速率频率frHz(周期:1/fr秒)按每一通道重复产生速率脉冲。延迟电路将按照每一通道将超声波聚集为波束状且决定发送指向性所需的延迟时间给予各速率脉冲。驱动脉冲产生电路按照基于各延迟的速率脉冲的定时向超声波探头1施加驱动脉冲。
另外,发送接收部2具有超声波接收用的未图示的放大器电路、A/D变换器、接收延迟电路以及加法器等。放大器电路接收来自超声波探头1的回波信号,并将接收到的回波信号按每一通道加以放大。A/D变换器将放大的回波信号按每一通道从模拟信号变换为数字信号。接收延迟电路对变换为数字信号的回波信号给予按每一通道聚集为束状且决定接收指向性所需要的延迟时间。加法器将给予了延迟时间的各回波信号相加。通过此加法处理,强调来自与回波信号的接收指向性相应的方向的反射分量,通过接收指向性与发送指向性形成超声波束。1个超声波束对应于1个超声波扫描线。每个扫描线的回波信号被供给到B模式处理部3与多普勒处理部5。
B模式处理部3通过对数放大来自发送接收部2的回波信号,包络线检波被对数放大后的回波信号,从而产生使用亮度来表现回波信号的强度的B模式信号数据。所产生的B模式信号数据被供给到B模式体数据产生部4。
B模式体数据产生部4根据来自B模式处理部3的B模式信号,产生与被检体相关的体数据(以下,称为B模式体数据)。具体而言,B模式体数据产生部4按照扫描线的位置信息在存储器上三维地配置B模式信号数据,插补扫描线间的B模式信号数据。通过该配置处理与插补处理,产生由多个体素构成的B模式体数据。各体素具有与原来(由来)的B模式信号数据的强度对应的体素值。所产生的B模式体数据被供给到散斑/噪声减少处理部7。
多普勒处理部5频率解析来自发送接收部2的回波信号,提取由多普勒效应所产生的血流、组织、造影剂回波分量,产生使用颜色来表现平均速度、分散、功率等血流信息的强度的多普勒信号数据。所产生的多普勒信号数据被供给到多普勒体数据产生部6。
多普勒体数据产生部6根据来自多普勒处理部5的多普勒信号,产生与被检体相关的体数据(以下,称为多普勒体数据)。具体而言,多普勒体数据产生部6按照扫描线的位置信息在存储器上三维地配置多普勒信号数据,插补扫描线间的多普勒信号数据。通过该配置处理与插补处理,产生由多个体素构成的多普勒体数据。各体素具有与原始(由来)的多普勒信号数据的强度对应的体素值。所产生的多普勒体数据被供给到散斑/噪声减少处理部7。
散斑/噪声减少处理部7对来自B模式体数据产生部4的B模式体数据、来自多普勒体数据产生部6的多普勒体数据执行散斑/噪声减少处理,产生散斑/噪声被减少了的B模式体数据、多普勒体数据。更详细而言,散斑/噪声减少处理部7通过减少B模式体数据、多普勒体数据中包含的具有三维各向同性结构的分量,强调具有三维各向异性结构的分量,从而减少这些体数据中的散斑、噪声。散斑、噪声被减少的B模式体数据、多普勒体数据被供给到三维图像处理部8。关于散斑/噪声减少处理部7的详细情况在后面进行叙述。另外,散斑/噪声减少处理部7并不是必需都对B模式体数据与多普勒体数据两者进行散斑减少处理,也可以只对任何一方进行。
三维图像处理部8对来自散斑/噪声减少处理部7的B模式体数据进行三维图像处理,产生二维B模式图像数据。另外,三维图像处理部8对来自散斑/噪声减少处理部7的多普勒体数据进行三维图像处理,产生二维多普勒图像数据。作为所使用的三维图像处理可以举出MPR(multi planar reconstruction:多平面重建)处理或CPR(curvedplanar reconstruction:曲面重建)处理、SPR(streched CPR)处理、体绘制、面绘制、MIP(maximum intensity projection:最大强度投影)等。所产生的B模式图像数据、多普勒图像数据被供给到显示部9。
显示部9显示来自三维图像处理部8的B模式图像。并且,显示部9显示来自三维图像处理部的多普勒图像。另外,显示部9也可以重叠显示B模式图像与多普勒图像。作为显示部9,例如能够适当使用CRT显示器或液晶显示器、有机EL显示器、等离子显示器等显示装置。
并且,超声波诊断装置具备存储部10。存储部10存储来自B模式体数据产生部4的B模式体数据与来自多普勒体数据产生部6的多普勒体数据。由存储部10存储的B模式体数据或多普勒体数据通过散斑/噪声减少处理部7读出,并被供给到散斑/噪声减少处理中。
并且,超声波诊断装置具备***控制部11。***控制部11综合地控制超声波诊断装置的各部分。例如,***控制部11为了执行本实施方式中特有的散斑/噪声减少处理而控制各部分。
散斑/噪声减少处理部7、三维图像处理部8、显示部9、存储部10以及***控制部11构成图像处理装置12。要是典型地来说,图像处理装置12为内置在与本实施方式相关的超声波诊断装置中的计算机。
其次,针对由散斑/噪声减少处理部7执行的散斑/噪声减少处理详细进行说明。如上所述,散斑/噪声减少处理部7能够分别对B模式体数据与多普勒体数据执行散斑减少处理。因此,以下,为了简单说明不对B模式体数据与多普勒体数据进行区别,仅仅称为体数据。另外,在不特别区别B模式图像与多普勒图像时,都称为超声波图像。
散斑/噪声减少处理部7对体数据中包含的各体素,按照与体素的三维方向性对应的滤波特性应用三维滤波。通过该三维滤波器,散斑/噪声减少处理部7减少体数据中包含的具有三维各向同性结构的分量,强调体数据中包含的具有各向异性结构的分量。另外,散斑/噪声减少处理的处理对象未必必需是三维各向同性结构分量与三维各向异性结构分量这两者。即、散斑/噪声减少处理部7也可以只进行三维各向同性结构分量的减少与三维各向异性结构分量的强调中的任何一方。
图2为表示散斑/噪声减少处理部7的结构的图。如图2所示,散斑/噪声减少处理部7为了进行多分辨率分解/合成而具有由多层构成的多重结构。在本实施方式中,多分辨率分解/合成的最高次数为等级3。然而,本实施方式并不需要限定于此。即、多分辨率分解/合成也可以在从第1等级到第n等级(其中,n为2以上的自然数)的范围进行。另外,在本实施方式中,作为多分辨率分解/合成的一个例子,采用小波变换/逆变换。然而,本实施方式并不需要限定于此。例如,作为多分辨率分解/合成,也可以采用拉普拉斯金字塔(Laplacian pyramid)法或Gabor变换/逆变换等已有的多分辨率分解/合成法。另外,本实施方式中的小波变换/逆变换指的是所谓的离散小波变换/逆变换。
散斑/噪声减少处理部7在每层上具备三维小波变换部71(71-1、71-2、71-3)、非线性各向异性扩散滤波器部73(73-1、73-2、73-3)、高频体素控制部75(75-1、75-2、75-3)以及三维小波逆变换部77(77-1、77-2、77-3)。
三维小波变换部71对输入体数据应用三维小波变换,并将输入体数据分解为1种低频体数据与7种高频体数据。
图3为三维小波变换的多分辨率分解的概念图。图4为表示在三维小波变换中的各坐标轴方向应用的滤波器的应用模式的图。如图3与图4所示,对分解前的体数据IV,在xyz正交坐标的各轴方向上分别应用一维低通(L)滤波器与高通(H)滤波器,并对体数据IV进行小波变换。通过这样的三维小波变换,体数据IV被分解为1种低频体数据DV1与7种高频体数据DV2至DV8。低频体数据DV1包含体数据IV具有的空间频率分量中的低频分量。各高频体数据DV2至DV8包含与体数据IV具有的空间频率分量中的至少一个方向相关的高频分量。例如,通过分别对x轴、y轴以及z轴上应用低通滤波器,从而关于x轴、y轴以及z轴方向产生只具有低频分量的低频体数据DV1。如果再举一个例子,则通过对x轴上应用高通滤波器,分别对y轴以及z轴上应用低通滤波器,从而产生强调x轴方向的高频分量的高频体数据DV2。分解后的各体数据DV1至DV8在每个坐标轴的样本数缩小为分解前的体数据IV在每个坐标轴的样本数的一半。
在三维小波变换部71不属于最高层(等级3)时,所产生的低频体数据被供给到1层上的三维小波变换部71。当三维小波变换部71属于最高层时,所产生的低频体数据被供给到同最高层的非线性各向异型扩散滤波器部73。另外,高频体数据被供给到属于同层的高频体素控制部75。
非线性各向异性扩散滤波器部73算出与所供给的低频体数据中包含的边缘分量的大小相关的信息。以下,将与边缘分量的大小相关的信息称为边缘信息。边缘信息被供给到同层的高频体素控制部75。另外,非线性各向异性扩散滤波器部73对所供给的低频体数据应用三维非线性各向异性扩散滤波器。非线性各向异性扩散滤波器强调低频体数据中包含的边缘分量,使低频体数据中包含的非边缘分量平滑化。非线性各向异性扩散滤波器为强调非线性地各向异性扩散的分量的滤波器。低频体数据中的边缘分量发生非线性地各向异性扩散。另一方面,散斑、噪声所代表的非边缘分量发生非线性地各向同性扩散。因此,当对低频体数据应用非线性各向异性扩散滤波器时,边缘分量的体素值增大,非边缘分量减少。为了使与本实施方式相关的非线性各向异性扩散滤波器能够一边提高散斑/噪声的减少精度一边提高处理速度而下了工夫。对此后面进行叙述。滤波后的低频体数据被供给到同层的三维小波逆变换部77。
在此,关于非线性各向异性扩散滤波器部73的处理详细说明。非线性各向异性扩散滤波器部73对低频体数据应用非线性各向异性扩散滤波器,并输出边缘信息与滤波后的低频体数据,该非线性各向异性扩散滤波器是考虑低频体数据内的结构物(体素区域)的方向性的滤波器。非线性各向异性扩散滤波器通过下面的扩散方程式(偏微分方程式)(1)来表示。
[数1]
★
I为要处理的体数据的体素值、为其梯度向量(gradient vector)、t为与处理有关的时刻。t在实际处理中,表示该扩散方程式的处理次数。在本实施方式中,处理次数t无论是几次都可以,但为了具体地进行说明,设为1次。D为扩散张量(Diffusion Tensor),如以下式(2)所示。
[数2]
★
(2)式中的λD1、λD2以及λD3为扩散张量D的固有值,R为扩散张量D的固有向量。固有向量R通过下面的(3)式来表示。
[数3]
★R=(ω1ω2ω3) …(3)
扩散张量D的固有值λD1为表示固有向量ω1所示的方向的扩散强度。同样,固有值λD2表示固有向量ω2所示的方向的扩散强度,固有值λD3表示固有向量ω3所示的方向的扩散强度。通过对每个体素控制固有值λD1、λD2以及λD3的值,从而控制非线性各向异性扩散滤波器的扩散强度。另外,固有向量R与非线性各向异性扩散滤波器的滤波器方向是相同含义。即,通过适当设定固有向量R,设定非线性各向异性扩散滤波器的所希望的滤波器方向。
另外,扩散张量D的固有向量ω1、ω2以及ω3与低频体数据中的体素的三维结构张量S的固有向量相等。另外,扩散张量D的固有值λD1、λD2以及λD3取决于根据后面所述的三维结构张量S的固有值所算出的边缘的大小。
三维结构张量S如下面的式(4)那样规定。
[数4]
★
Ix为表示沿着处理对象的体数据I的x方向的空间微分,Iy为表示沿着处理对象的体数据I的y方向的空间微分,Iz为表示沿着处理对象的体数据I的z方向的空间微分。Gρ为三维高斯函数,运算符“*”表示卷积。固有值λS1、λS2以及λS3为三维结构张量S的固有值,其大小关系设为λS1≥λS2≥λS3。空间微分Ix、Iy以及Iz的算出方法也可以不严格按照上述算出方法。例如,也可以用应用索贝尔滤波器(sobel filter)或多分辨率分解的高频分量来代替算出Ix、Iy以及Iz。
三维结构张量S的各要素s被求出时的固有值λS1、λS2以及λS3与固有向量R以及RT可以使用在线形代数学上熟知的方法来算出。即、式(4)中的3×3矩阵的固有值λS1、λS2以及λS3可以使用例如卡尔达诺(Cardano)方法等,通过实际地解答三维方程式来算出。另外,三维结构张量S由于为实对称矩阵,因此,固有值λS1、λS2以及λS3为实数、固有向量R以及RT为实向量。实相量R以及RT相互正交。
体数据中包含的三维结构物的代表性结构模式根据三维结构张量S的固有值λS1、λS2以及λS3的大小关系来分类。图5、图6以及图7为分别表示3个结构模式中的1个的示意图。如图5所示,时,结构物为各向异性结构物,具有平板状结构。如图6所示, 时,结构物为各向异性结构物,具有线状结构。如图7所示,λ时,结构物为各向同性结构物,具有球状结构。作为减少对象的散斑、噪声具有球状结构。另外,作为强调对象的边缘分量具有平面状结构或线状结构。
本发明者凭经验发现在对非线形各向异性扩散滤波器应用体数据的三维结构张量时,不需要特别地区别平板状结构与线状结构。那是因为考虑了第2固有值λS2而产生的超声波图像与没有考虑第2固有值λS2而产生的超声波图像相比,画质没有提高。因此,在本实施方式中,根据第1固有值λS1与第3固有值λS3的差分值,判别三维结构物的形状。即、该差分值DI>>0时,结构物设为具有各向异性结构,当时,设为具有各向同性结构。第1固有值λS1与第3固有值λS3的差分被用于算出边缘信息。边缘信息Pe如下面的式(5)所示。
[数5]
★
边缘信息Pe是通过使边缘的大小归一化为0至1的范围内而所算出的参数。边缘信息Pe的值越接近1各向异性越强,越接近0各向同性越强。即、边缘信息Pe的值越接近1结构物越接近边缘分量,越接近0结构物越接近非边缘分量。像这样边缘信息Pe可以说表示结构物的各向异性或各向同性(边缘分量或非边缘分量)的程度。
参数k表示边缘分量的提取程度。参数k能够由用户经由未图示的操作部任意设定。例如,当使参数k变小时,边缘分量变得容易被提取出。
使用边缘信息Pe算出非线性各向异性扩散滤波器的扩散强度,即、三维扩散张量D的固有值λD1、λD2以及λD3。固有值λD1通过下面的式(6)来表示、固有值λD2通过下面的式(7)来表示,λD3通过下面的式(8)来表示。
[数6]
λD1=β1(1-Pe)+β2·Pe …(6)
★λD2=β3(1-Pe)+β4·Pe …(7)
λD3=β5(1-Pe)+β6·Pe …(8)
式(6)、(7)以及(8)中的β(1-Pe)的部分表示非边缘分量,β·Pe部分表示边缘分量。对于作为减少对象的散斑或噪声所属的非边缘分量,需要使其不具有其方向性。即、必须使该非边缘分量不依赖方向进行扩散。因此,设定为β1=β3=β5>0。另一方面,对于作为强调对象的边缘分量必须更强调其方向性。即。需要使该边缘分量的垂直方向(固有值λD1的固有向量ω1所指的方向)尖锐化,使其以外的方向(固有值λD 2的固有向量ω2所指的方向与固有值λD3的固有向量ω3所指的方向)扩散。因此,β2设定为接近0的值,β4与β6设定为比β2大的规定值。这些参数β1、β2、β3、β4、β5以及β6分别能够经由未图示的操作部任意设定。
另外,如上所述,扩散张量D的固有向量ω1、ω2以及ω3与三维结构张量S的固有向量ω1、ω2以及ω3相等。通过以上可以算出扩散张量D的固有值λD1、λD2以及λD3与固有向量ω1、ω2以及ω3。因此,可以算出式(2)中的扩散张量D的各要素d,可以取得非线性各向异性扩散滤波器(扩散方程式(1))。与本实施方式相关的非线性各向异性扩散滤波器的固有值(扩散强度)如上所述,与表示各向异性或各向同性的程度的边缘信息Pe的大小对应地发生变化。即、通过对低频体数据应用与本实施方式相关的非线性各向异性扩散滤波器,强调具有低频体数据中的各向异性结构的分量,减少具有各向同性结构的分量。
图8为表示基于非线性各向异性扩散滤波器部73的滤波器处理的典型流程的图。另外,步骤S1至步骤S6的处理针对构成处理对象的体数据的各体素进行。
如图8所示,首先,非线性各向异性扩散滤波器部73输入作为低频体数据的处理对象的输入体素的体素值。并且,非线性各向异性扩散滤波器部73根据输入体素的体素值与其近旁的体素的体素值,沿着各坐标轴空间微分输入体素的体素值,算出微分值I(Ix,Iy,Iz)(步骤S1)。当算出微分值I时,非线性各向异性扩散滤波器部73如式(4)所示,将算出的微分值I与三维高斯函数Gρ进行卷积运算,算出结构张量S的各要素s(步骤S2)。另外,在步骤S2的计算中也包含三维高斯函数Gρ的计算。
当计算出三维结构张量S的各要素s时,非线性各向异性扩散滤波器部73如式(4)所示,对所算出的各要素s进行线性代数运算,从而算出三维结构张量S的固有值λs与固有向量ω(步骤S3)。当算出固有值λs时,非线性各向异性扩散滤波器部73如式(5)所示,利用算出的固有值λs与参数k算出边缘信息Pe(步骤S4)。边缘信息Pe被供给到高频体素控制部75。
当算出边缘信息Pe时,非线性各向异性扩散滤波器部73如式(6)、(7)以及(8)所示,利用边缘信息Pe与参数β,算出扩散强度、即扩散张量D的固有值λD1、λD2以及λD3(步骤S5)。当算出扩散张量D的固有值λD1、λD2以及λD3与扩散张量D(结构张量S)的固有向量ω1、ω2、ω3时,非线性各向异性扩散滤波器部73如式(2)所示,对所算出的固有值λD1、λD2以及λD3以及固有向量ω1、ω2以及ω3进行线性代数运算,从而算出扩散张量D的各要素d(步骤S6)。当算出扩散张量D的要素d时,非线性各向异性扩散滤波器部73根据算出的要素d与微分值I,对式(1)的扩散方程式进行数值解析运算,算出输出体素值(步骤S7)。具体而言,在时刻t,根据某体素与其近旁体素中的体素值与扩散张量的各要素值,求出时刻t+Δt中的其体素的新的体素值,其次将t+Δt作为新的t,重复1次到数次同样的计算。所算出的输出体素值被供给到三维小波逆变换部77。
当进行步骤S7时,改变下一个处理对象的输入体素,再次进行步骤S1至步骤S7。由此,当对构成处理对象的输入体数据的全部的体素进行步骤S1至步骤S7时,基于非线性各向异性扩散滤波器部73的处理结束。
即、非线性各向异性扩散滤波器部73对低频体数据沿着与三维结构张量的固有向量对应的滤波方向(扩散张量的固有向量所指的方向)应用具有与三维结构张量的固有值对应的扩散强度的三维非线性各向异性扩散滤波器,产生滤波后的低频体数据。
高频体素控制部75根据来自非线性各向异性扩散滤波器部73的边缘信息,控制来自三维小波变换部71的高频体数据中包含的体素的体素值。更详细而言,高频体素控制部75使高频体数据中包含的具有各向异性结构的分量的体素值增大,减少具有各向同性结构的分量的体素值。即、通过高频体素控制部75能够减少高频体数据中包含的散斑分量、噪声分量。
更详细而言,高频体素控制部75对滤波器应用输入体素值IH,输出输出体素值IH′。该滤波器通过下面的式(9)来表示。
[数7]
★I′H=A·Pe·IH+B·(1-Pe)·IH …(9)
在此,A为针对边缘分量的控制系数,B为针对散斑、噪声所代表的非边缘分量的控制系数。例如,控制系数A设定在1以上,控制系数B设定在1以下。即,边缘分量通过IH与A与Pe的积来强调。另外,非边缘分量通过IH与B与(1-Pe)的积来减少。通过该滤波器,强调了边缘分量并减少了非边缘分量的高频体数据被供给到三维小波逆变换部77。
另外,由于设定为A=0且B≠0,因此高频体素控制部75不强调边缘分量(各向异性结构分量),减少非边缘分量(各向同性结构分量)。通过这样减少非边缘分量,可以相对强调边缘分量。相反,由于设定为A≠0且B=0,因此高频体素控制部75强调边缘分量(各向异性结构分量),不减少非边缘分量(各向同性结构分量)。通过这样强调边缘分量,可以相对减少非边缘分量。
作为基于高频体素控制部75进行的边缘分量的其他强调方法,有以下方法。首先,将规定的边缘信息的大小设定为阈值,将具有阈值以上的边缘信息的体素设定为边缘分量,将具有阈值以下的边缘信息的体素设定为非边缘分量。并且,将属于边缘分量的体素的输出体素值改变成其输入体素值与控制系数A的积。换而言之,属于边缘分量的体素的输入体素值被置换为输入体素值与控制系数A的积。另外,将属于非边缘分量的体素的输出体素值改变成其输入体素值与控制系数B的积。属于非边缘分量的体素的输入体素值被置换为输入体素值与控制系数B的积。
三维小波逆变换部77将来自非线性各向异性扩散滤波器部73的低频体数据与来自高频体素控制部75的高频体数据进行多分辨率合成(作为本实施方式的典型例子是三维小波逆变换),将低频体数据与高频体数据合成,产生合成体数据。所产生的合成体数据与输入体数据相比,强调具有各向异性结构的分量,减少具有各向同性结构的分量。所产生的合成体数据在每个坐标轴的样本数被扩大为所输入的体数据在每个坐标轴的样本数的2倍。在三维小波逆变换部77不属于最下层时,所产生的合成体数据被供给到1层下的非线性各向异性扩散滤波器部73。在三维小波逆变换部77属于最下层时,所产生的合成体数据被供给到三维图像处理部8。
其次,说明各层中的三维小波变换部71、非线性各向异性扩散滤波器部73、高频体素控制部75以及三维小波逆变换部77的处理的流程。
如图2所示,等级1的三维小波变换部71-1对输入体数据进行三维小波变换,输出1种低频体数据与7种高频体数据。低频体数据被供给到等级2的三维小波变换部71-2,高频体数据被供给到等级1的高频体素控制部75-1。
等级2的三维小波变换部71-2对来自等级1的三维小波变换部71-1的低频体数据进行三维小波变换,输出1种低频体数据与7种高频体数据。低频体数据被供给到等级3的三维小波变换部71-3,高频体数据被供给到等级2的高频体素控制部75-2。
等级3的三维小波变换部71-3对来自等级2的三维小波变换部71-2的低频体数据进行三维小波变换,输出1种低频体数据与7种高频体数据。低频体数据被供给到等级3的非线性各向异性扩散滤波器部73-3,高频体数据被供给到等级3的高频体素控制部75-3。
其次,非线性各向异性扩散滤波器部73-3算出与来自三维小波变换部71-3的低频体数据相关的边缘信息。边缘信息被供给到等级3的高频体素控制部75-3。另外,非线性各向异性扩散滤波器部73-3对低频体数据应用非线性各向异性扩散滤波器。滤波后的低频体数据被供给到等级3的三维小波逆变换部77-3。
高频体素控制部75-3输入来自三维小波变换部71-3的7种高频体数据与来自非线性各向异性扩散滤波器部73-3的边缘信息。并且,高频体素控制部75-3根据输入的边缘信息,控制7种高频体数据的各个中包含的体素的体素值,输出控制体素值的7种高频体数据。7种高频体数据被供给到等级3的三维小波逆变换部77-3。
三维小波逆变换部77-3对来自非线性各向异性扩散滤波器部73-3的低频体数据与来自高频体素控制部75-3的高频体数据进行三维小波逆变换,输出单一的合成体数据。合成体数据作为低频体数据被供给到等级2的非线性各向异性扩散滤波器部73-2。
非线性各向异性扩散滤波器部73-2对来自三维小波逆变换部77-3的低频体数据进行与等级3同样的滤波器处理,将滤波器处理后的低频体数据供给到等级2的小波逆变换部77-2。另一方面,等级2的高频体素控制部75-2与等级3同样控制来自等级2的三维小波变换部71-2的7种高频体数据的体素值。并且,高频体素控制部75-2将控制体素值的7种高频体数据供给到等级2的小波逆变换部77-2。小波逆变换部77-2对来自非线性各向异性扩散滤波器部73-2的低频体数据与来自高频体素控制部75-2的7种高频体数据,使用与等级3同样的方法进行三维小波逆变换,输出单一的合成体数据。合成体数据作为低频体数据被供给到等级1的非线性各向异性扩散滤波器部73-1。
被供给到非线性各向异性扩散滤波器部73-1的低频体数据在进行与等级2、3同样的滤波器处理后,被供给到等级1的小波逆变换部77-1。另一方面,从等级1的三维小波变换部71-1输出的7种高频体数据在等级1的高频体素控制部75-1中与等级2或等级3同样进行体素值控制,并被供给到等级1的小波逆变换部77-1。并且,小波逆变换部77-1与等级2、3同样对低频体数据与7种高频体数据进行三维小波逆变换,输出单一的合成体数据。输出的合成体数据在每个坐标轴的样本数与被输入到散斑/噪声减少处理部7的体数据在每个坐标轴的样本数相等。
由等级1的三维小波逆变换部77-1产生的合成体数据被供给到三维图像处理部8。三维图像处理部8对合成体数据进行三维图像处理,产生超声波图像(二维图像)数据。所产生的超声波图像上的各向同性结构减少,各向异性结构被强调。即、超声波图像上的散斑、噪声减少。所产生的超声波图像通过显示部9予以显示。
根据以上所述结构,可以取得以下效果。
如上所述那样与本实施方式相关的超声波诊断装置对输入体数据的各体素根据三维结构张量算出边缘信息。超声波诊断装置根据边缘信息对每个体素算出三维滤波器(扩散方程式(1))。三维滤波器具有与体素的三维方向性对应的滤波器强度与滤波器方向。超声波诊断装置对输入体数据应用所算出的三维滤波器,并沿着输入体数据内的结构物的边缘方向(处理对象体素的接触面内的全部方向)进行平滑化,沿着与边缘方向垂直的方向(处理对象体素的接触面的法线方向)进行尖锐化。例如,结构物具有平面形状时,沿着结构物的平面(图5的ω2ω3平面)平滑化,沿着该平面的法线方向(图5的ω1方向)尖锐化。像这样超声波诊断装置对体数据可以应用与各结构物的形状对应的最佳的三维滤波器。
更详细而言,根据与本实施方式相关的超声波诊断装置,通过对体数据,沿着各坐标轴分别应用低通滤波器以及高通滤波器,执行三维多分辨率分解,将体数据分解为低频体数据与高频体数据。并且,超声波诊断装置对低频体数据应用非线性各向异性扩散滤波器。通过非线性各向异性扩散滤波器的应用,产生低频体数据中包含的边缘分量增大且散斑、噪声等非边缘分量减少的低频体数据。在应用非线性各向异性扩散滤波器时,需要判别与三维结构张量的固有值的大小关系对应的结构物的结构模式。与本实施方式相关的超声波诊断装置只应用3个固有值中的最大值与最小值来判别三维结构模式。由此,与使用3个固有值的情况相比,处理量下降,处理速度提高。这从作为超声波诊断装置的优点的图像显示实时性的观点来看也合适。
另外,与本发明相关的超声波诊断装置利用通过非线性各向异性扩散滤波器取得的边缘信息,增大高频体数据内的边缘分量,减少非边缘分量。即、在强调相对高频体数据的边缘分量中,通过利用预先取得的边缘信息,从而不需要重新确定边缘分量的位置信息等。由此与以往相比也可以提高处理速度。
并且,将强调了边缘的低频体数据与高频体数据继续进行多分辨率合成,产生非边缘分量被减少且边缘分量被强调的输出体数据。该输出体数据是利用三维多分辨率分解/合成、三维扩散滤波器而产生的。因此,与本实施方式相关的散斑/噪声减少处理与像以往那样对体数据应用二维滤波器的情况相比,在坐标轴方向上不受影响,与以往相比可以三维地均衡、均匀地减少散斑、噪声。
这样,与本实施方式相关的超声波诊断装置、图像处理装置、超声波诊断装置的控制方法以及图像处理方法实现提高体数据中包含的散斑、噪声的减少精度以及减少处理速度。
(变形例1)
以下,说明与本实施方式的变形例1相关的超声波诊断装置、图像处理装置12、超声波诊断装置的控制方法以及图像处理方法。另外,在以下的说明中,对具有与本实施方式大致相同的功能的构成要素,附加相同符号,只在需要的情况下重复说明。
与变形例1相关的超声波诊断装置对体数据进行三维图像处理后执行散斑/噪声减少处理。图9为表示与变形例1相关的超声波诊断装置的结构的图。B模式体数据产生部4将产生的B模式体数据供给到三维图像处理部8。多普勒体数据产生部6将产生的多普勒体数据供给到三维图像处理部8。三维图像处理部8对供给的B模式体数据或多普勒体数据等体数据进行三维图像处理,产生具有厚度的断面图像的数据。所产生的具有厚度的断面图像的数据被供给到散斑/噪声减少处理部7。
散斑/噪声减少处理部7对于供给的具有厚度的断面图像执行本实施方式特有的散斑/噪声减少处理,减少散斑、噪声。如果将具有厚度的断面图像的各断面作为xy平面,将厚度方向作为z轴方向,其处理方法与本实施方式中的三维图像处理方法相同。但是,在变形例1中,与x、y轴方向相比z轴方向的样本数相当少。将处理结果作为平面显示时,处理后的z轴方向的样本数为1。另外,在减少处理后进行MIP等其他处理后构成平面时,减少处理后也为具有厚度的断面图像,z轴方向的样本数为多个。散斑、噪声被减少的具有厚度的断面图像的数据再次被供给到三维图像处理部8。
三维图像处理部8对散斑、噪声被减少了的具有厚度的断面图像的数据进行三维图像处理,产生二维的超声波图像的数据。所产生的超声波图像的数据被供给到显示部9。
显示部9显示从三维图像处理部8供给的超声波图像。
这样,根据变形例1,能够提供实现提高体数据中包含的散斑、噪声的减少精度以及减少处理速度的超声波诊断装置、图像处理装置、超声波诊断装置的控制方法以及图像处理方法。
(变形例2)
与上述实施方式相关的图像处理装置12内置在超声波诊断装置中。然而,本实施方式并不需要限定于此。例如,图像处理装置12也可以内置在X射线计算机断层摄影装置或磁共振成像装置、核医学诊断装置中。另外,图像处理装置12也可以是经由网络与超声波诊断装置、X射线计算机断层摄影装置、磁共振成像装置以及核医学诊断装置等图像诊断装置连接的工作站。即、与本实施方式相关的处理对象的体数据并不取决于图像诊断装置的种类。与变形例2相关的图像处理装置12能够对超声波诊断装置或X射线计算机断层摄影装置、磁共振成像装置、核医学诊断装置等所有图像诊断装置产生的体数据执行上述散斑/噪声减少处理。
这样,根据变形例2,能够提供实现提高体数据中包含的散斑、噪声的减少精度以及减少处理速度的图像处理装置以及图像处理方法。
本领域技术人员容易想到其它优点和变更方式。因此,本发明就其更宽的方面而言不限于这里示出和说明的具体细节和代表性的实施方式。因此,在不背离由所附的权利要求书以及其等同物限定的一般发明概念的精神和范围的情况下,可以进行各种修改。
Claims (31)
1.一种超声波诊断装置,其特征在于,包括:
超声波探头;
发送接收部,经由上述超声波探头向被检体发送超声波,接收由来于由上述被检体反射的超声波的回波信号;
第1体数据产生部,根据上述接收到的回波信号产生第1体数据;
低频/高频体数据产生部,根据上述产生的第1体数据,产生在空间频率上的低频体数据与高频体数据;
滤波器部,计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算边缘信息,对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器;
高频体素控制部,根据上述产生的边缘信息,对上述高频体数据实施边缘强调处理;以及
第2体数据产生部,根据在上述滤波器部中应用了非线性各向异性扩散滤波器后的低频体数据与在上述高频体素控制部中实施了边缘强调处理后的高频体数据,产生散斑和噪声被减少了的第2体数据。
2.一种超声波诊断装置,其特征在于,包括:
超声波探头;
发送接收部,经由上述超声波探头向被检体发送超声波,接收由上述被检体反射的超声波,输出与上述接收到的超声波对应的回波信号;
第1体数据产生部,根据上述回波信号产生第1体数据;
第2体数据产生部,对上述第1体数据中包含的各个体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据;以及
图像数据产生部,根据上述第2体数据产生二维的超声波图像的数据。
3.根据权利要求2所述的超声波诊断装置,其特征在于:
上述第2体数据产生部根据上述体素的体素值与位于上述体素的近旁的多个近旁体素的体素值,算出上述体素的三维扩散强度与方向,根据上述算出的三维扩散强度与方向决定上述滤波特性。
4.根据权利要求2所述的超声波诊断装置,其特征在于:
上述第2体数据产生部,包括:
分解部,根据上述产生的第1体数据,产生在空间频率上的低频体数据与高频体数据;
滤波器部,计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算出边缘信息,对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器;
高频体素控制部,根据上述产生的边缘信息,对上述高频体数据实施边缘强调处理;以及
合成部,根据在上述滤波器部中应用了非线性各向异性扩散滤波器后的低频体数据与在上述高频体素控制部中实施了边缘强调处理后的高频体数据,产生散斑与噪声被减少了的第2体数据。
5.根据权利要求4所述的超声波诊断装置,其特征在于:
上述分解部对上述产生的第1体数据实施多分辨率分解,产生上述低频体数据与上述高频体数据,
上述合成部通过对在上述滤波器部中应用了非线性各向异性扩散滤波器后的低频体数据与在上述高频体素控制部中实施了边缘强调处理后的高频体数据实施多分辨率合成,从而产生上述第2体数据。
6.根据权利要求5所述的超声波诊断装置,其特征在于:
上述多分辨率分解与上述多分辨率合成是小波变换/逆变换、拉普拉斯金字塔法、或Gabor变换/逆变换。
7.根据权利要求2所述的超声波诊断装置,其特征在于,还包括:
显示部,显示上述超声波图像。
8.一种超声波诊断装置,其特征在于,包括:
超声波探头;
发送接收部,经由上述超声波探头向被检体发送超声波,接收由上述被检体反射的超声波,输出与上述接收的超声波对应的回波信号;
产生部,根据上述回波信号产生第1体数据;
减少部,通过减少上述第1体数据中包含的具有三维各向同性结构的分量,强调具有三维各向异性结构的分量,从而根据上述第1体数据产生第2体数据;以及
三维图像处理部,对上述第2体数据进行三维图像处理,产生二维的超声波图像的数据。
9.根据权利要求8所述的超声波诊断装置,其特征在于:
上述减少部,包括:
分解部,对上述第1体数据进行多分辨率分解,产生第1低频体数据与第1高频体数据;
算出部,根据上述第1低频体数据算出三维结构张量的固有值与固有向量;
滤波器部,对上述第1低频体数据沿着与上述固有向量对应的方向应用具有与上述固有值对应的扩散强度的滤波器,产生第2低频体数据;
控制部,根据上述固有值控制上述第1高频体数据中包含的体素的体素值,产生第2高频体数据;
合成部,将上述第2低频体数据与上述第2高频体数据进行多分辨率合成,产生上述第2体数据。
10.根据权利要求9所述的超声波诊断装置,其特征在于:
上述算出部算出上述三维结构张量的3个固有值;
上述滤波器部使用上述3个固有值中的最大固有值与最小固有值算出与上述第1低频体数据中包含的边缘分量的大小相关的边缘信息,对上述第1低频体数据应用具有与上述边缘信息的大小对应的上述扩散强度的非线性各向异性扩散滤波器。
11.根据权利要求9所述的超声波诊断装置,其特征在于:
上述多分辨率分解与上述多分辨率合成是小波变换/逆变换、拉普拉斯金字塔法、或Gabor变换/逆变换。
12.根据权利要求8所述的超声波诊断装置,其特征在于,还包括:
显示部,显示上述超声波图像。
13.一种图像处理装置,其特征在于,包括:
存储部,存储与被检体相关的第1体数据;
第2体数据产生部,对上述第1体数据中包含的各体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据;
图像数据产生部,根据上述第2体数据产生二维的超声波图像的数据。
14.根据权利要求13所述的图像处理装置,其特征在于:
上述第2体数据产生部根据上述体素的体素值与位于上述体素的近旁的多个近旁体素的体素值,算出上述体素的三维的扩散强度与方向,根据上述算出的三维的扩散强度与方向决定上述滤波特性。
15.根据权利要求13所述的图像处理装置,其特征在于:
上述第2体数据产生部,包括:
分解部,根据上述产生的第1体数据产生在空间频率上的低频体数据与高频体数据;
滤波器部,计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算出边缘信息,对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器;
高频体素控制部,根据上述产生的边缘信息对上述高频体数据实施边缘强调处理;
合成部,根据在上述滤波器部中应用了非线性各向异性扩散滤波器后的低频体数据与在上述高频体素控制部中实施了边缘强调处理后的高频体数据,产生散斑与噪声被减少的第2体数据。
16.根据权利要求15所述的图像处理装置,其特征在于:
上述分解部对上述产生的第1体数据实施多分辨率分解,产生上述低频体数据与上述高频体数据;
上述合成部对在上述滤波器部中应用了非线性各向异性扩散滤波器后的低频体数据与在上述高频体素控制部中实施了边缘强调处理后的高频体数据实施多分辨率合成,产生上述第2体数据。
17.根据权利要求16所述的图像处理装置,其特征在于:
上述多分辨率分解与上述多分辨率合成是小波变换/逆变换、拉普拉斯金字塔法、或Gabor变换/逆变换。
18.根据权利要求13所述的图像处理装置,其特征在于,还包括:
显示部,显示上述超声波图像。
19.一种图像处理装置,其特征在于,包括:
存储部,存储与被检体相关的第1体数据;
减少部,通过减少上述第1体数据中包含的具有三维各向同性结构的分量,强调具有三维各向异性结构的分量,从而根据上述第1体数据产生第2体数据;
三维图像处理部,对上述第2体数据进行三维图像处理,产生二维图像的数据。
20.根据权利要求19所述的图像处理装置,其特征在于:
上述减少部,包括:
分解部,对上述第1体数据进行多分辨率分解,产生第1低频体数据与第1高频体数据;
算出部,根据上述第1低频体数据算出三维结构张量的固有值与固有向量;
滤波器部,对上述第1低频体数据沿着与上述固有向量对应的方向应用具有与上述固有值对应的扩散强度的滤波器,产生第2低频体数据;
控制部,根据上述固有值控制上述第1高频体数据中包含的体素的体素值,产生第2高频体数据;
合成部,将上述第2低频体数据与上述第2高频体数据进行多分辨率合成,产生上述第2体数据。
21.根据权利要求20所述的图像处理装置,其特征在于:
上述算出部算出上述三维结构张量的3个固有值,
上述滤波器部使用上述3个固有值中的最大固有值与最小固有值,算出与上述第1低频体数据中包含的边缘分量的大小相关的边缘信息,对上述第1低频体数据应用具有与上述边缘信息的大小对应的上述扩散强度的非线性各向异性扩散滤波器。
22.根据权利要求20所述的图像处理装置,其特征在于:
上述多分辨率分解与上述多分辨率合成是小波变换/逆变换、拉普拉斯金字塔法、或Gabor变换/逆变换。
23.根据权利要求19所述的图像处理装置,其特征在于,还包括:
显示部,显示上述二维图像。
24.一种超声波诊断装置的控制方法,该超声波诊断装置具备:超声波探头、经由上述超声波探头发送接收超声波的发送接收部、产生第1体数据的第1体数据产生部、产生低频体数据与高频体数据的低频/高频体数据产生部、应用滤波器的滤波器部、控制高频体素的高频体素控制部、产生第2体数据的第2体数据产生部、***控制部,该超声波诊断装置的控制方法的特征在于:
上述***控制部控制上述发送接收部,以使得经由上述超声波探头向被检体发送超声波,经由上述超声波探头接收由上述被检体反射的超声波,输出与上述接收的超声波对应的回波信号,
上述***控制部控制上述第1体数据产生部,以使得根据上述输出的回波信号产生第1体数据,
上述***控制部控制上述低频/高频体数据产生部,以使得根据上述产生的第1体数据产生在空间频率上的低频体数据与高频体数据,
上述***控制部控制上述滤波器部,以使得计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,根据上述3个固有值中的2个固有值计算边缘信息,对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器,
上述***控制部控制上述高频体素控制部,以使得根据上述产生的边缘信息对上述高频体数据实施边缘强调处理,
上述***控制部控制上述第2体数据产生部,以使得根据应用了上述非线性各向异性扩散滤波器后的低频体数据与实施了上述边缘强调处理后的高频体数据,产生散斑与噪声被减少的第2体数据。
25.一种超声波诊断装置的控制方法,该超声波诊断装置具备:超声波探头、经由上超声波探头发送接收超声波的发送接收部、产生第1体数据的第1体数据产生部、产生第2体数据的第2体数据产生部、实施三维图像处理的三维图像处理部、***控制部,该超声波诊断装置的控制方法的特征在于:
上述***控制部控制上述接收部,以使得经由上述超声波探头向被检体发送超声波,经由上述超声波探头接收由上述被检体反射的超声波,输出与上述接收的超声波对应的回波信号,
上述***控制部控制上述第1体数据产生部,以使得根据上述输出的回波信号产生第1体数据,
上述***控制部控制上述第2体数据产生部,以使得通过对上述产生的第1体数据中包含的各个体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波产生第2体数据,
上述***控制部控制上述三维图像处理部,以使得根据上述产生的第2体数据产生二维超声波图像的数据。
26.根据权利要求25所述的超声波诊断装置的控制方法,其特征在于:
上述***控制部控制上述第2体数据产生部,以使得根据上述体素的体素值与位于上述体素的近旁的多个近旁体素的体素值算出上述体素的三维的扩散强度与方向,根据上述算出的三维的扩散强度与方向决定上述滤波特性。
27.根据权利要求25所述的超声波诊断装置的控制方法,其特征在于:
上述***控制部控制上述第2体数据产生部,以使得根据上述产生的第1体数据产生在空间频率上的低频体数据与高频体数据,
计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,
根据上述3个固有值中的2个固有值计算出边缘信息,
对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器,
根据上述产生的边缘信息对上述高频体数据实施边缘强调处理,
根据应用了上述非线性各向异性扩散滤波器后的低频体数据与实施了上述边缘强调处理后的高频体数据,产生散斑与噪声被减少了的第2体数据。
28.一种图像处理方法,其特征在于:
根据与被检体相关的第1体数据,产生在空间频率上的低频体数据与高频体数据,
计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,
根据上述3个固有值中的2个固有值计算边缘信息,
对上述低频体数据应用了具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器,
根据上述产生的边缘信息,对上述高频体数据实施边缘强调处理,
根据应用了上述非线性各向异性扩散滤波器后的低频体数据与实施了上述边缘强调处理的高频体数据,产生散斑与噪声被减少了的第2体数据。
29.一种图像处理方法,其特征在于:
通过对与被检体相关的第1体数据中包含的各体素,按照与上述体素的三维方向性对应的滤波特性应用三维滤波,产生第2体数据,
根据上述第2体数据产生二维图像的数据。
30.根据权利要求29所述的图像处理方法,其特征在于:
产生上述第2体数据是通过以下操作实现的,即:根据上述体素的体素值与位于上述体素的近旁的多个近旁体素的体素值,算出上述体素的三维的扩散强度与方向,根据上述算出的三维的扩散强度与方向决定上述滤波特性。
31.根据权利要求29所述的图像处理方法,其特征在于:
产生上述第2体数据是通过以下操作实现的,即:
根据上述产生的第1体数据,产生在空间频率上的低频体数据与高频体数据,
计算用于分类上述产生的低频体数据中包含的三维结构物的三维结构模式的3个固有值,
根据上述3个固有值中的2个固有值计算边缘信息,
对上述低频体数据应用具有与上述计算出的边缘信息对应的滤波器特性的非线性各向异性扩散滤波器,
根据上述产生的边缘信息对上述高频体数据实施边缘强调处理,
根据应用了上述非线性各向异性扩散滤波器后的低频体数据与实施了上述边缘强调处理后的高频体数据,产生散斑与噪声被减少了的第2体数据。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009051145 | 2009-03-04 | ||
JP2009-051145 | 2009-03-04 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101822547A true CN101822547A (zh) | 2010-09-08 |
CN101822547B CN101822547B (zh) | 2012-11-14 |
Family
ID=42678848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101288357A Active CN101822547B (zh) | 2009-03-04 | 2010-03-04 | 超声波诊断装置及其控制方法、图像处理装置和方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9113826B2 (zh) |
JP (1) | JP5525867B2 (zh) |
CN (1) | CN101822547B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102523732A (zh) * | 2010-09-09 | 2012-06-27 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置和医用图像处理方法 |
CN102958445A (zh) * | 2011-06-30 | 2013-03-06 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
CN103118599A (zh) * | 2011-09-20 | 2013-05-22 | 株式会社东芝 | 图像处理装置以及医用图像诊断装置 |
CN103169496A (zh) * | 2011-12-26 | 2013-06-26 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置以及处理方法 |
CN103202714A (zh) * | 2012-01-16 | 2013-07-17 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置以及医用图像处理方法 |
CN103338707A (zh) * | 2011-01-26 | 2013-10-02 | 株式会社日立医疗器械 | 超声波诊断装置以及图像处理方法 |
CN103927715A (zh) * | 2014-03-14 | 2014-07-16 | 中瑞科技(常州)有限公司 | 一种超声图像斑点噪声抑制方法 |
CN112136182A (zh) * | 2017-11-16 | 2020-12-25 | 杨晓东 | 基于Gabor光学相干层析术血流成像的***和方法 |
CN113491536A (zh) * | 2020-04-01 | 2021-10-12 | 株式会社日立制作所 | 超声波摄像装置以及图像处理装置 |
CN114022444A (zh) * | 2021-11-04 | 2022-02-08 | 青岛海信医疗设备股份有限公司 | 超声图像的处理方法及装置 |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9184369B2 (en) | 2008-09-18 | 2015-11-10 | Fujifilm Sonosite, Inc. | Methods for manufacturing ultrasound transducers and other components |
US8306296B2 (en) * | 2009-04-30 | 2012-11-06 | Medison Co., Ltd. | Clutter signal filtering using eigenvectors in an ultrasound system |
JP5832737B2 (ja) * | 2010-11-01 | 2015-12-16 | 株式会社東芝 | 超音波診断装置及び超音波画像処理装置 |
KR101819028B1 (ko) | 2011-07-11 | 2018-01-16 | 삼성전자주식회사 | 초음파 영상 처리 방법 및 장치 |
JP6041127B2 (ja) * | 2011-08-29 | 2016-12-07 | 国立大学法人北海道大学 | 柱状物体抽出方法、柱状物体抽出プログラム、及び柱状物体抽出装置 |
CN102999884B (zh) * | 2011-09-16 | 2015-11-25 | 株式会社东芝 | 图像处理设备和方法 |
KR101881924B1 (ko) | 2011-12-29 | 2018-08-27 | 삼성전자주식회사 | 초음파 영상 처리 방법 및 장치 |
WO2014119412A1 (ja) * | 2013-01-30 | 2014-08-07 | 株式会社 日立メディコ | 医用画像処理装置及び医用画像撮像装置 |
JP6274495B2 (ja) * | 2013-10-25 | 2018-02-07 | キヤノンメディカルシステムズ株式会社 | 画像処理装置および超音波診断装置 |
ES2550532B1 (es) * | 2014-05-08 | 2016-08-19 | Universidad De Valladolid | Método de filtrado y realce de imágenes ultrasónicas 2D, 3D y 3D+t |
KR101612306B1 (ko) | 2014-06-24 | 2016-04-15 | 서울여자대학교 산학협력단 | 3차원 자동 유방 초음파 영상에서의 피부층 자동 검출 시스템 및 그 검출 방법 |
US10542961B2 (en) | 2015-06-15 | 2020-01-28 | The Research Foundation For The State University Of New York | System and method for infrasonic cardiac monitoring |
JP6594075B2 (ja) * | 2015-07-22 | 2019-10-23 | キヤノン株式会社 | 画像処理装置、撮像システム、画像処理方法 |
CA3077311A1 (en) | 2017-09-29 | 2019-04-04 | C. R. Bard, Inc. | Apparatus and method for tracking a medical ultrasonic object |
TWI743411B (zh) * | 2017-11-08 | 2021-10-21 | 美商富士膠片索諾聲公司 | 具有高頻細節的超音波系統 |
JP7079680B2 (ja) * | 2018-07-05 | 2022-06-02 | 富士フイルムヘルスケア株式会社 | 超音波撮像装置、および、画像処理装置 |
US11941784B1 (en) * | 2020-04-27 | 2024-03-26 | Canon U.S.A., Inc. | Image processing apparatus, method, and storage medium for spectrally encoded endoscopy ghost noise reducton |
JP7536557B2 (ja) * | 2020-09-02 | 2024-08-20 | キヤノンメディカルシステムズ株式会社 | 画像処理装置及び超音波診断装置 |
US11908085B2 (en) * | 2021-01-15 | 2024-02-20 | Siemens Medical Solutions Usa, Inc. | Entropy-dependent adaptive image filtering |
CN112749493A (zh) * | 2021-01-25 | 2021-05-04 | 中国自然资源航空物探遥感中心 | 基于全磁力梯度张量特征值的地质体边界检测方法及*** |
CN113191961B (zh) * | 2021-03-19 | 2022-10-25 | 聚融医疗科技(杭州)有限公司 | 一种基于双边滤波和特征识别的超声去噪方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1474358A (zh) * | 2002-08-08 | 2004-02-11 | GEҽ��ϵͳ����������˾ | 三维空间滤波设备和方法 |
CN1493258A (zh) * | 2002-10-28 | 2004-05-05 | ��ʽ���綫֥ | 图像处理设备与超声波诊断设备 |
CN1521695A (zh) * | 2003-02-13 | 2004-08-18 | ��ʽ���綫֥ | 从图像中减少噪声的图像处理设备 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6806705B2 (en) * | 2002-05-15 | 2004-10-19 | Koninklijke Philips Electronics N.V. | Diffusion tensor magnetic resonance imaging including local weighted interpolation |
KR100646715B1 (ko) * | 2004-10-18 | 2006-11-23 | 한국과학기술원 | 후처리를 통한 2차원 초음파 영상의 화질 개선 방법 |
KR100752333B1 (ko) * | 2005-01-24 | 2007-08-28 | 주식회사 메디슨 | 3차원 초음파 도플러 이미지의 화질 개선 방법 |
JP4789854B2 (ja) * | 2007-05-09 | 2011-10-12 | 株式会社日立メディコ | 超音波診断装置および超音波診断装置の画質改善方法 |
JP5049773B2 (ja) * | 2007-12-27 | 2012-10-17 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、超音波画像処理プログラム |
-
2010
- 2010-03-02 JP JP2010045682A patent/JP5525867B2/ja active Active
- 2010-03-03 US US12/716,601 patent/US9113826B2/en active Active
- 2010-03-04 CN CN2010101288357A patent/CN101822547B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1474358A (zh) * | 2002-08-08 | 2004-02-11 | GEҽ��ϵͳ����������˾ | 三维空间滤波设备和方法 |
CN1493258A (zh) * | 2002-10-28 | 2004-05-05 | ��ʽ���綫֥ | 图像处理设备与超声波诊断设备 |
CN1521695A (zh) * | 2003-02-13 | 2004-08-18 | ��ʽ���綫֥ | 从图像中减少噪声的图像处理设备 |
Non-Patent Citations (1)
Title |
---|
KHALED Z. ET AL: "Real-Time Speckle Reduction and Coherence Enhancement in Ultrasound Imaging via Nonlinear Anisotropic Diffusion, Vol 49, No 9, IEEE transactions on biomedical engineering", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9795364B2 (en) | 2010-09-09 | 2017-10-24 | Toshiba Medical Systems Corporation | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method |
CN102523732A (zh) * | 2010-09-09 | 2012-06-27 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置和医用图像处理方法 |
CN102523732B (zh) * | 2010-09-09 | 2014-11-26 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置和医用图像处理方法 |
CN103338707A (zh) * | 2011-01-26 | 2013-10-02 | 株式会社日立医疗器械 | 超声波诊断装置以及图像处理方法 |
CN103338707B (zh) * | 2011-01-26 | 2015-09-30 | 株式会社日立医疗器械 | 超声波诊断装置以及图像处理方法 |
CN102958445A (zh) * | 2011-06-30 | 2013-03-06 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
CN102958445B (zh) * | 2011-06-30 | 2015-05-13 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
CN103118599B (zh) * | 2011-09-20 | 2015-07-22 | 株式会社东芝 | 图像处理装置以及医用图像诊断装置 |
CN103118599A (zh) * | 2011-09-20 | 2013-05-22 | 株式会社东芝 | 图像处理装置以及医用图像诊断装置 |
CN103169496B (zh) * | 2011-12-26 | 2015-10-28 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置以及处理方法 |
US9585636B2 (en) | 2011-12-26 | 2017-03-07 | Toshiba Medical Systems Corporation | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method |
CN103169496A (zh) * | 2011-12-26 | 2013-06-26 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置以及处理方法 |
CN103202714B (zh) * | 2012-01-16 | 2015-06-17 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置以及医用图像处理方法 |
CN103202714A (zh) * | 2012-01-16 | 2013-07-17 | 株式会社东芝 | 超声波诊断装置、医用图像处理装置以及医用图像处理方法 |
CN103927715A (zh) * | 2014-03-14 | 2014-07-16 | 中瑞科技(常州)有限公司 | 一种超声图像斑点噪声抑制方法 |
CN112136182A (zh) * | 2017-11-16 | 2020-12-25 | 杨晓东 | 基于Gabor光学相干层析术血流成像的***和方法 |
CN113491536A (zh) * | 2020-04-01 | 2021-10-12 | 株式会社日立制作所 | 超声波摄像装置以及图像处理装置 |
CN114022444A (zh) * | 2021-11-04 | 2022-02-08 | 青岛海信医疗设备股份有限公司 | 超声图像的处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101822547B (zh) | 2012-11-14 |
JP5525867B2 (ja) | 2014-06-18 |
US9113826B2 (en) | 2015-08-25 |
JP2010227554A (ja) | 2010-10-14 |
US20100228129A1 (en) | 2010-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101822547B (zh) | 超声波诊断装置及其控制方法、图像处理装置和方法 | |
US11373303B2 (en) | Systems and methods for ultrasound imaging | |
CN101879073B (zh) | 超声波诊断装置 | |
CN102370498B (zh) | 超声波诊断装置及超声波图像处理装置 | |
Sakas et al. | Preprocessing and volume rendering of 3D ultrasonic data | |
Shams et al. | Real-time simulation of medical ultrasound from CT images | |
Kutter et al. | Visualization and GPU-accelerated simulation of medical ultrasound from CT images | |
US20220233169A1 (en) | Ultrasonic blood flow parameter displaying method, and ultrasonic imaging system therefor | |
CN101467897B (zh) | 超声波诊断装置、超声波图像处理装置、超声波图像处理方法 | |
US20060184021A1 (en) | Method of improving the quality of a three-dimensional ultrasound doppler image | |
CN102436672B (zh) | 超声波图像处理装置 | |
CN102958445B (zh) | 超声波诊断装置、图像处理装置以及图像处理方法 | |
CN103476345B (zh) | 超声波诊断装置 | |
CN106600550B (zh) | 超声图像处理方法及*** | |
US20110144497A1 (en) | Ultrasonic 3-dimensional image reconstruction method and ultrasonic wave system thereof | |
CN102523732B (zh) | 超声波诊断装置、医用图像处理装置和医用图像处理方法 | |
US6458082B1 (en) | System and method for the display of ultrasound data | |
Kuo et al. | Interactive volume rendering of real-time three-dimensional ultrasound images | |
Starkov et al. | Ultrasound simulation with animated anatomical models and on-the-fly fusion with real images via path-tracing | |
Chen et al. | Neural Implicit Representation for Three-dimensional Ultrasound Carotid Surface Reconstruction using Unsigned Distance Function | |
ÖZKUL | Single-image bayesian restoration and multi-image super-resolution restoration for b-mode ultrasound images using an accurate system model | |
Watkin et al. | Three-dimensional reconstruction and enhancement of arbitrarily oriented and positioned 2D medical ultrasonic images | |
JPS62150480A (ja) | 立体画像表示装置 | |
Lai et al. | Designs and Implementation of Three Dimensional Nuchal Translucency | |
Lu et al. | A Novel Algorithm for Ultrasonic Stereovision Data Enhancement |
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: 20160713 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. |