CN101396287B - 超声波诊断装置以及超声波诊断图像生成方法 - Google Patents

超声波诊断装置以及超声波诊断图像生成方法 Download PDF

Info

Publication number
CN101396287B
CN101396287B CN200810161771.3A CN200810161771A CN101396287B CN 101396287 B CN101396287 B CN 101396287B CN 200810161771 A CN200810161771 A CN 200810161771A CN 101396287 B CN101396287 B CN 101396287B
Authority
CN
China
Prior art keywords
mentioned
image
filtering operation
coefficient
frames
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
CN200810161771.3A
Other languages
English (en)
Other versions
CN101396287A (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 CN101396287A publication Critical patent/CN101396287A/zh
Application granted granted Critical
Publication of CN101396287B publication Critical patent/CN101396287B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54—Control of the diagnostic device
    • 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995—Combining images from different aspect angles, e.g. spatial compounding
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023—Details of receivers
    • G01S7/52025—Details of receivers for pulse systems
    • G01S7/52026—Extracting wanted echo signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053—Display arrangements
    • G01S7/52057—Cathode ray tube displays
    • G01S7/5206—Two-dimensional coordinated display of distance and direction; B-scan display
    • G01S7/52065—Compound scan display, e.g. panoramic imaging
    • 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]

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供一种超声波诊断装置以及程序。在超声波诊断装置中，具备进行如下处理的图像合成部(14)。即，图像合成部(14)对多个帧的图像，分别通过规定的变换处理进行多重分辨率分析，针对通过该多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理，对该滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像。

Description

超声波诊断装置以及超声波诊断图像生成方法
相关申请的交叉引用 
本申请是基于2007年9月28日提交的在先的日本专利申请No.2007-256338，并且要求该专利申请的优先权，该专利申请的全部内容以引用的方式并入本文中。 
技术领域
本发明涉及超声波诊断装置以及程序，特别涉及对执行超声波扫描而取得的多个帧的图像进行合成处理的超声波诊断装置以及程序。 
背景技术
超声波诊断装置是向被检体内放射从内置于超声波探测器中的压电振子发生的超声波，利用上述压电振子接收由于被检体组织的声阻抗的差异而产生的反射信号并显示在监视器上的装置。根据本超声波诊断装置，可以仅通过将超声波探测器接触到体表的简单操作而容易地得到实时的图像数据。因此，超声波诊断装置被广泛用于脏器的功能诊断和形态诊断。 
由于超声波脉冲反射法和超声波多普勒法这2个大的技术开发，利用来自被检体的组织或血球的反射波得到生物体信息的超声波诊断方法得到了急速的进步。使用这些技术得到的B模式图像数据和彩色多普勒图像数据在当今的超声波诊断中成为不可欠缺的一部分。 
另外，在通过超声波诊断装置得到的超声波图像数据中，发生由于来自被检体内的多个反射体的反射超声波的随机干涉而产生的粒状性的噪声(所谓散斑噪声)。公知该散斑噪声与超声波诊断装置的系统噪声一起妨碍超声波诊断。 
作为降低这样的散斑噪声的方法，例如可以举出复合扫描(compound scan)法。该复合扫描法进一步被分类成空间复合扫描和频率复合扫描。 
上述空间复合扫描是对被检体的同一部位从多个不同的方向进行超声波发送接收，并对所得到的多个图像数据进行加法合成而生成显示用的图像数据的技术。另一方面，上述频率复合扫描是对被检体的同一部位利用不同的超声波频率收集多张图像数据，并对这些收集到的图像数据进行加法合成而生成显示用图像数据的技术。 
作为如此对图像数据进行合成处理的技术，例如在特开平8-280688号公报中公开有如下那样的技术。即、在特开平8-280688号公报中，公开出一种发生所复合的超声波影像的方法，包括：将各个连续的影像帧分割成多个子影像区域，对该子影像区域的局部移动向量进行估计评价的阶段；根据所估计评价的局部移动向量估计评价总括性影像移动的阶段；以及根据所估计评价的总括性影像移动显示所合成的XFOV影像的阶段。 
根据该特开平8-280688号公报公开的方法，可以发生被表现成所扩大的视场(XFOV)影像的大的复合超声波影像。 
另外，作为降低上述散斑噪声的方法，例如在特开2005-296331号公报中公开出如下那样的技术。即，在特开2005-296331号公报中，公开出一种超声波诊断装置，具备：图像数据生成单元，对被检体的多个方向进行超声波发送接收而生成图像数据；小波变换单元，对上述图像数据进行小波变换而计算出各分辨率等级下的低频小波变换系数和多个高频小波变换系数；边缘信息检测单元，根据上述高频小波变换系数对上述图像数据的边缘信息进行检测；系数处理单元，根据所检测出的上述边缘信息对上述高频小波变换系数进行加权运算；小波逆变换单元，通过使用了加权运算后的上述高频小波变换系数的小波逆变换而生成被边缘强调的图像数据；以及显示单元，显示上述被边缘强调的图像数据。 
根据该特开2005-296331号公报公开的超声波诊断装置，可以在超声波图像数据中实现降低有效的散斑噪声和强调边缘。
另外，在通过超声波诊断装置来显示组织像时，作为对执行超声波扫描而得到的图像数据合成并显示时的合成方法，主要公知以下3种合成方法。 
(合成方法1)空间复合扫描 
在空间复合扫描中，如图8所示，固定保持超声波探测器101，针对每个帧从不同的方向进行扫描。如果像这样从不同的方向进行超声波扫描，则即使是来自相同位置的回波，干涉条件也不同，所以散斑的相关变少。因此，通过针对所重复的区域(第1扫描图像103与第2扫描图像105的重复区域107)进行平均值运算或最大值运算，可以减少散斑。 
另外，与超声波波束垂直的组织边界通过镜面反射，发生更清楚的边界回波。因此，通过从多个方向照射超声波波束，可以清楚地显示来自各种角度的组织边界的回波。 
(合成方法2)全景图(panoramic view) 
所谓全景图是特开平8-280688号公报中示出的那样的成像法，是如下的技术：如图9所示使超声波探测器101沿着扫描断面的方向移动，根据所得到的图像数据对超声波探测器101的动作进行检测，使图像数据并行移动且旋转移动来进行定位而执行合成处理，从而显示更宽的视场。当前，一般进行利用最新的图像来改写的合成处理。 
(合成方法)组合对焦(combination focus) 
所谓组合对焦是指，如图10所示，从同一方向而并非从不同方向，改变发送聚焦而进行多次超声波波束的发送接收，对通过这些发送接收取得的图像进行合成的技术。即、在组合对焦中，例如关于近距离区域的图像111A利用近距离聚焦的超声波波束111来取得，关于远距离区域的图像113A利用远距离聚焦的超声波波束113来取得，并合成这些图像。 
但是，上述图像数据的合成方法分别具有如下的问题。 
首先，关于空间复合扫描，在合成图像数据时，如图8所示在与该合成相关的图像数据(第1扫描图像103与第2扫描图像105)之 间，产生所重复的区域107和不重复的区域109。此处，关于不重复的区域109，直接显示与该合成相关的图像数据。另一方面，关于重复的区域107，例如利用平均值运算或最大值运算等图像合成方法来进行图像合成处理。 
此处，在作为图像合成方法采用最大值运算的情况下，所重复的区域107与不重复的区域109之间的连接点(边界)处的显示成为自然的显示。但是，如果在一个图像数据中存在灰度高的伪(artifact)像，则该伪像被直接显示。另一方面，在作为图像合成方法采用平均值运算的情况下，所重复的区域107与不重复的区域109的连接点(边界)处的显示看起来不自然。进而，作为图像合成方法，不论在采样最大值运算的情况下，还是在采用平均值运算的情况下，该合成处理后的图像数据都成为某处模糊的显示，具有成为对观察者造成上述那样的印象的图像数据这样的问题。 
另外，在使用了全景图情况下，并产生了所重复的区域的情况下，如上所述对该区域的图像数据进行改写。因此，不产生起因于合成处理的图像的模糊。但是，尽管存在所重复的区域的图像数据，但不将其利用时，无法说充分地享受到通过对图像数据进行合成处理而得到的优点。换言之，通过利用所重复的区域的图像数据，可以得到与通过上述空间复合扫描得到的效果同样的效果。但是，还同时造成与空间复合扫描所具有的问题同样的问题，所以无论如何必需解决空间复合扫描所具有的问题。 
进而，例如如使用了全景图的情况那样，在从不同的位置扫描而取得的图像数据中，起因于分别取得的位置，而在该图像数据中存在清晰的区域和不清晰的区域。 
例如如图11所示，如果针对在关心区域131的右侧存在不清晰区域的第1图像133、和在关心区域131的左侧存在不清晰区域的第2图像135，执行平均值运算或最大值运算来进行合成处理，则仅可以得到不清晰的区域扩大的图像数据。 
因此，在这样的情况下，优选通过不使用与该合成处理相关的图 像数据中存在的不清晰的区域，而仅使用清晰的区域来进行合成处理，得到仅具有清晰的区域的图像数据137。但是，当前，自动地进行这样的合成处理的技术仍未被公开过。 
而且，在使用组合聚焦的情况下，由于将图像数据分割而合成，所以存在在该合成处理后的图像数据上，观察到连接点(边界)这样的问题。尽管还存在针对连接点(边界)依次改变权重来相加这样的方法，但更优选为与改变权重来相加的方法相比对分辨率提高的贡献更大那样的合成处理方法。 
另外，特开平8-280688号公报以及特开2005-296331号公报公开的技术当然不能够解决参照图8至图11说明的上述问题。 
发明内容
本发明是鉴于上述情况而完成的，其目的在于提供一种超声波图像诊断装置以及程序，在通过合成处理制成的图像数据中，实现与该合成处理相关的图像数据彼此重复的区域与不重复的区域之间的边界不能观察成不连续、高亮度的伪像不被强调、以及即使在与该合成处理相关的图像数据中包含具有不清晰的区域的图像数据的情况下也不产生起因于此的不清晰的区域。 
为了达成上述目的，本发明的第1方式的超声波诊断装置以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：分析部，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；滤波运算部，对通过由上述分析部执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及图像生成部，对由上述滤波运算部执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像。 
为了达成上述目的，本发明的第2方式的超声波诊断图像生成方法以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫 描，根据接收回波生成多个帧的图像，其特征在于，具备：分析步骤，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；滤波运算步骤，对通过由上述分析步骤执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及图像生成步骤，对由上述滤波运算步骤执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像。 
根据本发明，可以提供一种超声波图像诊断装置以及程序，在通过合成处理制成的图像数据中，实现与该合成处理相关的图像数据彼此重复的区域与不重复的区域之间的边界不能观察成不连续、高亮度的伪像不被强调、以及即使在与该合成处理相关的图像数据中包含具有不清晰的区域的图像数据的情况下也不产生起因于此的不清晰的区域。 
本发明的优点将在随后的描述中阐述，并且部分地将从所述描述显而易见，或者可通过本发明的实践而得知。可通过尤其是后面指出的手段和组合而实现和获得本发明的优点。 
被并入说明书中并构成说明书一部分的附图示出本发明的实施例，并且与上面给出的一般描述以及下面给出的对实施例的详细描述一起用于解释本发明的原理。 
附图说明
图1是示出本发明的一个实施方式的超声波诊断装置的结构例子的图。 
图2是示出图像合成单元的输入输出的图。 
图3是示出由图像合成单元进行的图像数据的合成处理的流程的图。 
图4是示出本发明的一个实施方式中的图像数据的合成处理的原理的图。 
图5是示出三维扫描的概念的图。
图6是示出插入散斑降低和强调边缘的滤波处理时的由图像合成单元进行的图像数据的合成处理的流程的图。 
图7是示出2级组合聚焦法的概念的图。 
图8是示出空间复合扫描的概念的图。 
图9是示出全景图的概念的图。 
图10是示出组合聚焦的概念的图。 
图11是示出以往技术具有的问题点的图。 
具体实施方式
以下，参照附图对本发明的一个实施方式进行说明。 
图1是示出该实施方式的超声波诊断装置的结构例子的图。如该图所示，该实施方式的超声波诊断装置1具备超声波探测器2、发送电路4、接收电路6、B模式处理单元8、彩色多普勒处理单元10、坐标变换单元12、图像合成单元14、监视器16、控制单元18。 
上述超声波探测器2是承担在与被检体之间照射、反射的超声波信号的发送接收波的设备(探触子)，由作为电气/机械可逆性变换元件的压电陶瓷等压电元件形成。该超声波探测器2例如由在前端部具备阵列状地排列的多个压电元件的相控阵列类型的部件构成。由此，超声波探测器2将所供给的脉冲驱动电压变换成超声波脉冲信号而沿着被检体的扫描区域内的期望方向发送，并且将从被检体反射来的超声波信号变换成与其对应的电压的回波信号。 
上述发送电路4根据基于上述控制单元18的控制信号，以针对每个发送通道附加了规定的发送延迟时间的定时，向上述超声波探测器2的各压电元件发送驱动信号。由此从上述超声波探测器2的各压电元件朝向被检体内发送超声波信号。 
上述接收电路6将包括根据超声波信号的发送波在被检体内的声阻抗的失配面反射，而由组织内的散乱体散乱的分量等的超声波回波信号，经由上述超声波探测器2的各压电元件接收而作为与该信号对应的电压量的回波信号，对该回波信号实施接收延迟、加法处理而生 成接收信号并分别输出到上述B模式处理单元8以及彩色多普勒处理单元10。 
上述B模式处理单元8对来自上述接收电路6的输出信号进行包络线检波，将其检波出的信号作为承担被检体内的组织的形态信息的信号，输出到上述坐标变换单元12。 
上述彩色多普勒处理单元10根据来自上述接收电路6的输出信号对速度信息进行频率分析，将该分析结果作为承担被检体内的血流或组织的移动速度信息的信号，输出到上述坐标变换单元12。 
上述坐标变换单元12将来自上述B模式处理单元8以及上述彩色多普勒处理单元10的输出信号向显示坐标进行坐标变换后输出到上述图像合成单元14。 
上述图像合成单元14接收来自上述坐标变换单元12的信号，在控制单元18的控制下，进行与B模式像以及CFM(彩色血流成像，Color Flow Mapping)像相关的各种图像处理以及合成处理(在后面详述)，进而根据这些图像进行各种定量分析和测量后进行将表示其结果的信息附加到图像上等图像处理，将通过这些图像处理得到的图像信号变换成TV用的扫描信号而输出到监视器16。 
上述控制单元18具有作为具备与未图示的内部总线连接的CPU(处理器)、存储器(RAM/ROM)、硬盘装置、可移动体介质(CD-ROM、软盘(注册商标)、存储卡等)的驱动装置、及其外围装置的计算机的功能。上述控制单元18按照在检查时预先编程的步骤对该超声波诊断装置1整体的动作进行控制。 
通过以上说明的结构，在上述监视器16上，显示超声波图像(包括与B模式像以及CFM像相关的各种图像、与其测量、分析结果相关的信息等)。 
另外，在空间复合扫描中，通常针对每个帧改变上述超声波探测器2的倾斜扫描的方向(角度)。但是，针对每个帧改变扫描仅设对B模式像，而CFM像中针对每个帧设为相同的扫描。即，此处仅对B模式进行空间复合扫描以及合成处理。以下，仅对B模式像的处理进 行详细说明。 
以下，对作为本第1实施方式的超声波诊断装置1的主要特征之一的上述图像合成单元14所执行的合成处理进行详细说明。 
首先，关于由上述图像合成单元14执行的合成处理的输入输出关系例如成为图2所示那样的输入输出关系。即，通过N帧(或大于等于N帧)的运算来合成图像，输出与输入相同帧数的图像数据。更具体而言，如图2所示合成图像数据A1与图像数据B1而生成图像数据F1，合成图像数据B1与图像数据A2而生成图像数据F2，合成图像数据A2与图像数据B2而生成图像数据F3。 
另外，图2所示的例子是进行2个方向的合成处理的例子，但当然也可以由本一个实施方式的超声波诊断装置进行N个方向的合成处理。 
接下来，参照图3所示的流程图，对由上述图像合成单元14执行的合成处理的流程进行说明。此处，以对(k+1)帧的图像数据进行合成处理的情况为例子进行说明。另外，参照图4在后面对各步骤中的处理进行详细说明。 
首先，上述图像合成单元14如果被输入最新的图像数据(第n帧的图像)(步骤S1)，则进行L等级的离散二维小波变换(步骤S2)。 
接下来，上述图像合成单元14去除包括与合成处理相关的图像彼此重叠的区域中存在的边缘分量的边缘区域(步骤S3)。通常，如果在图像的边缘部中进行小波变换，则在高频侧发生边缘分量。而且，在合成处理后的图像中，上述图像的边缘部当然早已并非边缘部。因此，无需上述边缘部中的边缘分量。因此，在步骤S3中去除边缘分量。 
之后，上述图像合成单元14进行后述的规定的运算处理(针对(k+1)张的帧的各小波变换系数的每一个进行运算的处理)(步骤S4)，对该运算处理结果进行二维离散逆小波变换(步骤S5)，得到合成处理后的图像数据、即A(0)。 
以下，参照图4对本一个实施方式中的图像数据的合成处理的原理进行详细说明。此处，以对图像a以及图像b这2个图像进行合成 处理的情况为例子进行说明。具体而言，上述图像合成单元14如下所述进行图像a以及图像b的多重分辨率分析(此处为等级2的二维离散小波变换)(步骤S2)。 
即、上述图像合成单元14首先通过等级1的二维离散小波变换，针对图像a以及图像b，分解成A1(Approximation)、H1(Horizontaldetail)、V1(Vertical detail)、D1(Diagonal detail)的小波系数。接下来，仅对上述A1，进一步进行二维离散小波变换。由此，如图4所示，得到A2、H2、V2、D2。针对图像a以及图像b这双方分别进行这些处理。 
接下来，进行去除在图4中用椭圆形状51包围示出的区域、即边缘区域(将该边缘区域中的高频侧的系数设为0)的处理(步骤S3)。 
另外，该边缘区域是在合成处理后的图像中不成为边缘部的边缘区域。换言之，该边缘区域对于合成处理前的某1张图像来说是该图像的边缘部，但在合成处理后的图像中是不成为边缘部的区域。另外，该区域在该合成处理前成为既知的区域。同样地，在进行了小波变换的情况下，该图像的边缘部造成影响的范围的区域也是在该合成处理前既知的区域。 
根据以上说明的步骤S1至步骤S3的处理，当然，从相同图像可以得到相同结果。因此，通过依次保存当前的结果，无需每次再次运算，而可以得到在图3中用虚线包围的框内所示的(n-k)帧至(n-1)帧的离散二维小波变换的值。 
而且，上述步骤S4中的运算处理以及上述步骤S5中的二维离散逆小波变换例如是进行如下那样的运算的处理。 
此处，用func(a1、a2、...、aN、M、Level)表示运算函数，用idwt2表示离散二维逆小波变换。其中，a1、a2、...、aN表示输入数据，M＝0表示Approximation，M＝1表示Horizontal detail，M＝2表示Vertical detail，M＝3表示Diagonal detail，Level表示小波展开数。
A(L)＝func(A(L，n)，A(L，n-1)，···，A(L，n-k)，O，L) 
H(L)＝func(H(L，n)，H(L，n-1)，···，H(L，n-k)，1，L) 
V(L)＝func(V(L，n)，V(L，n-1)，···，V(L，n-k)，2，L) 
D(L)＝func(D(L，n)，D(L，n-1)，···，D(L，n-k)，3，L) 
H(L-1)＝func(H(L-1，n)，H(L-1，n-1)，···，H(L-1，n-k)，1，L-1) 
V(L-1)＝func(V(L-1，n)，V(L-1，n-1)，···，V(L-1，n-k)，2，L-1) 
D(L-1)＝func(D(L-1，n)，D(L-1，n-1)，···，D(L-1，n-k)，3，L-1) 
········ 
H(1)＝func(H(L，n)，H(L，n-1)，···，H(L，n-k)，1，1) 
V(1)＝func(V(L，n)，V(L，n-1)，···，V(L，n-k)，2，1) 
D(1)＝func(D(L，n)，D(L，n-1)，···，D(L，n-k)，3，1) 
A(L-1)＝idwt2(A(L)，H(L)，V(L)，D(L)) 
A(L-2)＝idwt2(A(L-1)，H(L-1)，V(L-1)，D(L-1)) 
········ 
A(0)＝idwt2(A(1)，H(1)，V(1)，D(1)) 
此处，上述A(0)是通过合成处理最终制成的图像数据。 
另外，作为上述的运算函数的一个例子，在全部Level中，可以举出下式。 
func(a1，a2，…，aN，M，Level)＝Mean(a1，a2，…，aN)M＝0 
func(a1，a2，…，aN，M，Level)＝AbsMax(a1，a2，…，aN)M＝1，2，3 
其中，Mean(a，b)是计算a与b的平均值的运算，是在a＝0的情况下输出b的值，在b＝0的情况下输出a的值的运算。另外，AbsMax(a，b)是输出a以及b中的绝对值大的一方的运算。 
另外，更具体而言，在上述步骤S4中，例如进行下式那样的运算。
A2＝Mean(A2a，A2b) 
H2＝AbsMax(H2a，H2b) 
V2＝AbsMax(V2a，V2b) 
D2＝AbsMax(D2a，D2b) 
H1＝AbsMax(H1a，H1b) 
V1＝AbsMax(V1a，V1b) 
D1＝AbsMax(D1a，D1b) 
然后，在进行了上述的运算之后，上述图像合成单元14针对上述的系数A2、H2、V2、D2进行二维离散逆小波变换而计算出A1，针对A1、H1、V1、D1进行二维离散逆小波变换(步骤S5)，得到合成处理后的图像数据、即A(0)。 
以下，对通过进行上述步骤S4中的运算处理而得到的效果进行说明。 
首先，通过针对上述Approximation的系数进行平均值运算，可以得到平均亮度不变化这样的效果、以及用与合成处理相关的各图像间的平均值来表示图像的大致构造这样的效果。 
另外，通过对上述H、V、D的detail系数进行最大值运算，选择出与图像处理相关的各图像中的具有最大分辨率的图像，针对各像素的每一个自动地选择出具有最大的分辨率的图像。因此，作为整体，重合的边界也观察不到，而被观察成自然的重合，在各像素处显示出所重合的图像中的具有最大分辨率的图像。因此，可以得到分辨率提高地进行观察这样的效果。 
即、可以达成图像数据彼此重复的区域与不重复的区域之间的边界不被观察成不连续、高亮度的伪像不被强调、以及即使在图像数据中包含具有不清晰的区域的图像数据的情况下也不产生起因于此的不清晰的区域这样的目的。 
通过这样的运算，达成在复合扫描中期待的通过镜面反射实现的图像的清晰化。此处，为了进一步得到散斑去除这样的效果，利用进行上述那样的小波变换的做法，对小波收缩(Wavelet Shrinkage)这 样的图像处理技术并用即可。等级低的detail系数多数情况下为噪声，所以通过进行小波收缩，可以降低或去除散斑。 
另外，作为进行这样的处理时的上述func函数，例如可以举出下式的函数。 
func(a1，a2，…，aN，M，Level)＝AbsMax(a1，a2，…，aN)*c(M，Level)M＝1，2，3 
此处，c(M，Level)表示小波系数的权重，例如在Level为1那样小的情况下取近似0的值，在Level大的情况下取近似1的值。另外，c(M，Level)的值也可以根据通过M指定的Horizontal、Vertical、Diagonal，分别改变值。进而，为了强调组织的边缘等边缘，c(M，Level)的值也可以设为比1大的值。 
另外，以上说明的由上述图像合成单元14执行的处理例如也可以适用于图5所示那样的相互重复的通过多次三维扫描而取得的图像数据。 
图5是示出一边进行三维扫描一边移动上述超声波探测器2而从不同的窗口观察同一组织的样子、即重复临时将上述超声波探测器2从被检体离开后放置到该被检体的其他位置而再次取得图像这样的步骤的样子的图。 
根据该方法，对被检体中的同一区域，从不同的位置实施多次三维扫描，三维地进行通过那些扫描取得的图像的定位，从而可以得到完全同一断面或同一体素的数据。在该情况下，关于同一区域存在多个体素数据，对那些多个体素数据或所切出的多个二维图像数据进行合成处理。而且，即使在采用这样的方法的情况下，也当然可以应用由上述图像合成单元14执行的上述合成处理。 
另外，由上述图像合成单元14执行的上述一连串的合成处理通过程序化或在程序化后将该程序读入存储介质，从而易于作为与超声波诊断装置1独立的软件产品单体来销售和发布，并且还可以将本一个实施方式涉及的技术利用在其他硬件上。 
如上述说明，根据本一个实施方式，在通过合成处理制成的图像数据中，可以提供实现了如下目的的超声波图像诊断装置以及程序： 与该合成处理相关的图像数据彼此重复的区域与不重复的区域之间的边界不会被观察成不连续、高亮度的伪像不被强调、以及即使在与该合成处理相关的图像数据中包含具有不清晰的区域的图像数据的情况下也不产生起因于此的不清晰的区域。 
另外，在通过上述的由上述图像合成单元14执行的合成处理而得到的图像数据中，该合成处理前的图像数据彼此重叠的区域与不重叠的区域之间的边界自然且分辨率提高。分辨率提高的理由为，该合成处理后的图像数据具有合成处理前的图像数据中的高频分量中的最大分量，所以在分辨率最高的状态下被合成处理。 
以上，根据一个实施方式说明了本发明，但本发明不限于上述的实施方式，当然可以在本发明的要旨的范围内实现各种变形以及应用。 
(第1变形例) 
作为合成处理中的运算方法，除了上述的运算方法以外，例如也可以在所有情况下采用使用Mean运算的方法、或在所有情况下使用AbsMax运算的方法。 
在所有的情况下使用Mean运算的方法是与不进行小波变换而仅取Mean的处理基本上等同的处理。但是，在本一个实施方式中，如上所述通过小波系数进行所重复的一侧的边缘区域去除处理，所以所重复的区域中的边界不会被看起来不连续。 
在所有情况下使用AbsMax运算的方法与不进行小波变换而仅取Max的处理不同，针对各像素的每一个输出具有最大的分辨率的图像。但是，有时平均亮度上升，或残留高亮度的伪像。 
(第2变形例) 
在上述那样的合成处理的运算方法中，仅使用各像素的每一个的系数来进行了合成处理，但当然也可以使用该像素的附近的像素的系数来进行合成处理。 
例如，作为2张图像的位置(x，y)处的通过AbsMax进行的计算方法，计算出第1张图像中的包括位置(x，y)的周围5×5点的平均值a1，并且计算出第2张图像中的包括位置(x，y)的周围5×5点 的平均值a2，将上述a1以及上述a2中的绝对值大的一方的值设为位置(x，y)处的输出值。通过该处理，可以降低噪声被强调而观察到的情况。 
(第3变形例) 
在上述的合成处理的运算方法中，在多重分辨率分析中使用了离散小波变换。在一般的离散小波变换中，对该图像数据实施LPF(LowPass Filter，低通滤波)和HPF(High Pass Filter，高通滤波)，将各输出欠采样成一半。因此，信息量在该变换前和变换后没有变化。 
但是，如果在通过小波收缩处理了系数之后，进行逆小波变换，则有时在该图像上会观察到块状的伪像。为了防止该现象，也可以代替离散小波变换，而例如进行平稳(stationary)小波变换。 
在平稳小波变换中不进行欠采样。因此，在对二维图像进行处理的情况下，信息量通过1次变换成为4倍。但是，即使在进行了小波收缩之后进行逆变换的情况下，也可以得到在该图像上不能观察到上述块状的伪像这样的效果。 
(第4变形例) 
在合成处理中的上述的运算方法中，在多重分辨率分析中使用了离散小波变换。但是，在多重分辨率变换中，也可以代替小波变换，而例如使用拉普拉斯(Laplacian)金字塔、RoL(Ratio of Low pass，比率低通)金字塔、或梯度(Gradient)金字塔这样的金字塔变换(Pyramid transform)。 
小波变换与金字塔变换之间的相异点在于，小波变换为正交变换，与其相对金字塔变换未必是正交变换。但是，对于小波变换以及金字塔变换这双方，在都可以通过那些变换进行多重分辨率分析的点上是共通的。 
(第5变形例) 
如图6所示，也可以将用于散斑降低和强调边缘的非线性扩散滤波处理(步骤S11)插入到上述步骤S4中的运算处理之前或上述步骤S4中的运算处理之后的任意一方或两方。另外，在图6中，为便于说 明，示出成在上述步骤S4的前后双方插入了非线性扩散滤波处理(步骤S11)。 
即，如例如特开2006-116307中公开的那样，通过在合成处理之前插入将多重分辨率分析与非线性扩散滤波组合的处理，可以得到散斑少且边缘被强调的易于观察的合成图像。在进行这样的处理的情况下，由于已经进行了多重分辨率分析，所以可以节省多重分辨率分析所需的时间。 
(第6变形例) 
在以上说明的例子中，假设使用了复合扫描的情况，但不限于复合扫描，例如即使在一边进行全景图或三维扫描一边移动超声波探测器2而从不同的窗口观察同一组织的情况下也可以同样地应用上述一个实施方式。 
另外，在全景图的情况下，有时需要从图像估计出移动量，而使该图像并行移动和/或旋转移动。在上述特开平8-280688号公报中示出其详细内容。即使在从其他位置进行三维扫描的情况下，通过三维地并行移动和/或旋转移动，可以将断面定位到期望的断面。在这样操作而使断面一致之后，可以进行在上述一个实施方式中说明的图像的合成处理。 
(第7变形例) 
也可以在通过组合聚焦法执行的合成处理中应用上述一个实施方式。在该情况下，使与合成处理相关的图像的边界交叠，进行在上述一个实施方式中说明的由图像合成单元14执行的合成处理。 
图7是示出2级的组合聚焦法的概念的图。使通过近距离聚焦的超声波波束111取得的图像51A与通过远距离聚焦的超声波波束53取得的图像53A之间的边界交叠而制成交叠区域55A，针对该交叠区域55A，使用在上述一个实施方式中说明的合成处理方法。 
由此，可以期待提高合成处理后的图像中的分辨率且降低散斑。在该情况下，也可以交换由上述图像合成单元14执行的处理与由上述坐标变换单元12执行的处理之间的顺序。
(第8变形例) 
也可以在通过频率复合扫描执行的合成处理中应用上述一个实施方式。另外，基波图像与2次高次谐波图像之间的合成处理也是频率复合的一种。根据频率复合扫描，通常由于各扫描范围一致，所以可以与组合聚焦法同样地应用上述一个实施方式。 
本领域的设计人员很容易想到附加的优点和修改。因此，从更广义的方面来说，本发明并不限于本文中所示出和描述的具体细节和代表性实施例。从而，在不脱离由所附权利要求及其等同形式所限定的总体发明构思的精神和范围的情况下可以进行各种修改。

Claims (9)

1.一种超声波诊断装置，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析部，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算部，对通过由上述分析部执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及
图像生成部，对由上述滤波运算部执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像，
由上述分析部进行的上述规定的变换处理是指，基于小波变换或金字塔变换的变换处理，
针对上述各系数的每一个进行的所述滤波运算处理是使用各像素的系数和该各像素的附近像素的系数来进行的处理。
2.一种超声波诊断装置，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析部，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算部，对通过由上述分析部执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及
图像生成部，对由上述滤波运算部执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像，
由上述滤波运算处理部执行的上述滤波运算处理是指，针对通过上述多重分辨率分析取得的上述系数中的低频侧的系数取平均值的运算，针对高频侧的系数输出绝对值大的值的运算处理。
3.一种超声波诊断装置，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析部，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算部，对通过由上述分析部执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及
图像生成部，对由上述滤波运算部执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像，
由上述滤波运算处理部执行的上述滤波运算处理是指，针对所有上述系数输出绝对值大的值的运算处理。
4.一种超声波诊断装置，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析部，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算部，对通过由上述分析部执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及
图像生成部，对由上述滤波运算部执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像，
上述滤波运算处理部在决定上述多个帧的图像中的规定位置的系数的输出值时，还利用该规定位置的附近的系数来进行上述滤波运算处理。
5.一种超声波诊断装置，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析部，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算部，对通过由上述分析部执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；
图像生成部，对由上述滤波运算部执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像；以及
边缘分量去除部，在由上述分析部实施了上述多重分辨率分析之后，且在由上述滤波运算处理部实施上述滤波运算处理之前，在上述多个帧的图像中的相互重复的区域的边缘部中，去除起因于由上述分析部执行的上述规定的变换处理而发生的边缘分量。
6.一种超声波诊断图像生成方法，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析步骤，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算步骤，对通过由上述分析步骤执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及
图像生成步骤，对由上述滤波运算步骤执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像，
上述滤波运算步骤中执行的上述滤波运算处理是指，针对通过上述多重分辨率分析取得的上述系数中的低频侧的系数取平均值的运算，针对高频侧的系数输出绝对值大的值的运算处理。
7.一种超声波诊断图像生成方法，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析步骤，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算步骤，对通过由上述分析步骤执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；以及
图像生成步骤，对由上述滤波运算步骤执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像，
上述滤波运算处理步骤中执行的上述滤波运算处理是指，针对所有上述系数输出绝对值大的值的运算处理。
8.一种超声波诊断图像生成方法，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析步骤，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算步骤，对通过由上述分析步骤执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；
图像生成步骤，对由上述滤波运算步骤执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像；以及
上述滤波运算步骤是在决定上述多个帧的图像中的规定位置的系数的输出值时，还利用该规定位置的附近的系数来进行上述滤波运算处理的步骤。
9.一种超声波诊断图像生成方法，以使被检体中的同一区域至少一部分重复的方式进行多次超声波扫描，根据接收回波生成多个帧的图像，其特征在于，具备：
分析步骤，对上述多个帧的图像，分别通过规定的变换处理进行多重分辨率分析；
滤波运算步骤，对通过由上述分析步骤执行的多重分辨率分析取得的各分辨率的各系数的每一个，使用上述多个帧的图像间的对应的系数彼此分别进行滤波运算处理；
图像生成步骤，对由上述滤波运算步骤执行的滤波运算处理结果，实施上述规定的变换处理的逆变换处理而生成1个帧的图像；以及
边缘分量去除步骤，在由上述分析步骤实施了上述多重分辨率分析之后，且由上述滤波运算处理步骤实施上述滤波运算处理之前，在上述多个帧的图像中的相互重复的区域的边缘部中，去除起因于上述分析步骤中执行的上述规定的变换处理而发生的边缘分量。
CN200810161771.3A 2007-09-28 2008-09-26 超声波诊断装置以及超声波诊断图像生成方法 Active CN101396287B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2007-256338 2007-09-28
JP2007256338 2007-09-28
JP2007256338A JP5002397B2 (ja) 2007-09-28 2007-09-28 超音波診断装置及びプログラム

Publications (2)

Publication Number Publication Date
CN101396287A CN101396287A (zh) 2009-04-01
CN101396287B true CN101396287B (zh) 2012-02-08

Family

ID=40509163

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810161771.3A Active CN101396287B (zh) 2007-09-28 2008-09-26 超声波诊断装置以及超声波诊断图像生成方法

Country Status (3)

Country Link
US (1) US8647275B2 (zh)
JP (1) JP5002397B2 (zh)
CN (1) CN101396287B (zh)

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5269517B2 (ja) * 2008-08-14 2013-08-21 株式会社東芝 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP5449852B2 (ja) 2009-05-08 2014-03-19 株式会社東芝 超音波診断装置
US9255914B2 (en) 2009-08-13 2016-02-09 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus and program
US9173629B2 (en) 2009-11-18 2015-11-03 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus
KR101120820B1 (ko) 2009-11-19 2012-03-22 삼성메디슨 주식회사 초음파 공간 합성 영상을 제공하는 초음파 시스템 및 방법
JP2012071115A (ja) * 2010-08-30 2012-04-12 Toshiba Corp 超音波診断装置
JP5832737B2 (ja) * 2010-11-01 2015-12-16 株式会社東芝 超音波診断装置及び超音波画像処理装置
JP5283725B2 (ja) * 2011-03-28 2013-09-04 富士フイルム株式会社 超音波診断装置
US8900148B2 (en) 2011-03-09 2014-12-02 Fujifilm Corporation Ultrasound diagnostic apparatus
KR101797038B1 (ko) 2011-04-22 2017-11-13 삼성전자주식회사 진단영상을 생성하는 방법, 이를 수행하는 장치, 진단시스템 및 의료영상시스템
JP5984260B2 (ja) * 2011-09-20 2016-09-06 東芝メディカルシステムズ株式会社 画像処理装置及び医用画像診断装置
MX2014005408A (es) * 2011-11-02 2015-02-12 Seno Medical Instr Inc Sistema de imagenología de modalidad doble para mapeo funcional y anatómico corregistrado.
US20130289381A1 (en) 2011-11-02 2013-10-31 Seno Medical Instruments, Inc. Dual modality imaging system for coregistered functional and anatomical mapping
CN104137144B (zh) * 2012-02-29 2017-03-08 独立行政法人科学技术振兴机构 图像处理用数字滤波器及字符串倾斜错觉生成装置
WO2013154079A1 (ja) * 2012-04-11 2013-10-17 株式会社東芝 超音波診断装置
US9453908B2 (en) * 2012-06-12 2016-09-27 University Of Rochester Superresolution imaging of scatterers in pulse-echo imaging
JP6415852B2 (ja) * 2013-07-12 2018-10-31 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用画像処理装置及び医用画像処理方法
JP6114663B2 (ja) * 2013-08-27 2017-04-12 富士フイルム株式会社 超音波診断装置および超音波画像生成方法
JP6274495B2 (ja) * 2013-10-25 2018-02-07 キヤノンメディカルシステムズ株式会社 画像処理装置および超音波診断装置
JP6386739B2 (ja) * 2014-02-13 2018-09-05 キヤノンメディカルシステムズ株式会社 超音波診断装置
CN104546008B (zh) * 2015-02-02 2016-09-21 声泰特（成都）科技有限公司 一种基波/谐波融合与空间复合相结合的成像方法
JP6596907B2 (ja) * 2015-05-01 2019-10-30 コニカミノルタ株式会社 超音波画像診断装置
CN106725607A (zh) * 2017-02-17 2017-05-31 深圳华声医疗技术股份有限公司 超声多普勒参数优化方法与超声多普勒装置
JP6767902B2 (ja) * 2017-03-17 2020-10-14 株式会社日立製作所 超音波診断装置及びプログラム
WO2019057981A1 (en) * 2017-09-25 2019-03-28 Koninklijke Philips N.V. METHODS AND SYSTEMS FOR PARASITIC ACOUSTIC ECHO FILTERING AND RANDOM NOISE
CN114052786B (zh) * 2021-10-22 2023-02-07 武汉联影医疗科技有限公司 一种超声波发射方法和系统
CN114287963A (zh) * 2021-12-30 2022-04-08 重庆海扶医疗科技股份有限公司 图像处理方法、装置、电子设备及计算机可读介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1739623A1 (en) * 2005-07-01 2007-01-03 Medison Co., Ltd. Hierarchical motion estimation method and ultrasound imaging system using the same

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06154218A (ja) * 1992-11-24 1994-06-03 Matsushita Electric Ind Co Ltd 超音波診断装置
JP3974946B2 (ja) * 1994-04-08 2007-09-12 オリンパス株式会社 画像分類装置
US5619998A (en) * 1994-09-23 1997-04-15 General Electric Company Enhanced method for reducing ultrasound speckle noise using wavelet transform
US5575286A (en) 1995-03-31 1996-11-19 Siemens Medical Systems, Inc. Method and apparatus for generating large compound ultrasound image
KR19990005283A (ko) * 1997-06-30 1999-01-25 배순훈 웨이브렛 변환 영역에서의 계층적 움직임 추정 방법
US6511426B1 (en) * 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6042545A (en) * 1998-11-25 2000-03-28 Acuson Corporation Medical diagnostic ultrasound system and method for transform ultrasound processing
ITFI20020034A1 (it) * 2002-02-27 2003-08-27 Actis Active Sensors S R L Metodo e dispositivo per l'analisi spettrale di un segnale ecografico
US6676605B2 (en) * 2002-06-07 2004-01-13 Diagnostic Ultrasound Bladder wall thickness measurement system and methods
JP4362327B2 (ja) * 2002-09-24 2009-11-11 パナソニック株式会社 画像符号化方法および画像符号化装置
EP1593087A4 (en) * 2003-01-30 2006-10-04 Chase Medical Lp METHOD AND SYSTEM FOR IMAGE PROCESSING AND CONTOUR EVALUATION
US20050277835A1 (en) * 2003-05-30 2005-12-15 Angelsen Bjorn A Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties
US20050053305A1 (en) * 2003-09-10 2005-03-10 Yadong Li Systems and methods for implementing a speckle reduction filter
JP4632685B2 (ja) * 2004-04-12 2011-02-16 株式会社東芝 超音波診断装置及び画像データ処理装置
KR100646715B1 (ko) * 2004-10-18 2006-11-23 한국과학기술원 후처리를 통한 2차원 초음파 영상의 화질 개선 방법
JP4590293B2 (ja) * 2005-04-11 2010-12-01 富士フイルム株式会社 超音波観測装置
KR100815977B1 (ko) * 2005-04-20 2008-03-24 주식회사 메디슨 연속되는 영상으로부터 움직임을 추정하는 방법 및 시스템
EP1880364A1 (en) * 2005-05-12 2008-01-23 Bracco Imaging S.P.A. Method for coding pixels or voxels of a digital image and a method for processing digital images
GB0514715D0 (en) * 2005-07-18 2005-08-24 Isis Innovation Combination of images
US20070083114A1 (en) * 2005-08-26 2007-04-12 The University Of Connecticut Systems and methods for image resolution enhancement

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1739623A1 (en) * 2005-07-01 2007-01-03 Medison Co., Ltd. Hierarchical motion estimation method and ultrasound imaging system using the same

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特开2005-296331A 2005.10.27

Also Published As

Publication number Publication date
CN101396287A (zh) 2009-04-01
JP5002397B2 (ja) 2012-08-15
US8647275B2 (en) 2014-02-11
JP2009082469A (ja) 2009-04-23
US20090088638A1 (en) 2009-04-02

Similar Documents

Publication Publication Date Title
CN101396287B (zh) 超声波诊断装置以及超声波诊断图像生成方法
KR101205107B1 (ko) 스페클 감소 필터의 구현 방법, 스페클 감소 필터링 장치 및 초음파 촬상 시스템
US8021301B2 (en) Ultrasonic image processing apparatus, ultrasonic image processing method and ultrasonic image processing program
US7404797B2 (en) Ultrasonic imaging system and ultrasonic signal processing method
CN105455843B (zh) 在超声成像中的阴影抑制
US6872181B2 (en) Compound image display system and method
US9585636B2 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method
JP5586375B2 (ja) 超音波診断装置、及びプログラム
WO2012049124A2 (en) Methods and systems for producing compounded ultrasound images
JPH05228145A (ja) 超音波診断装置
US20120053464A1 (en) Ultrasonic image processing apparatus
WO2007092054A2 (en) Method and apparatus to visualize the coronary arteries using ultrasound
JPH1147132A (ja) 三次元画像形成方法及び超音波画像処理装置
JP7366221B2 (ja) 解剖学的、機能的及び血行動態イメージングのためのトリプルモード超音波イメージング
JP7014044B2 (ja) 超音波診断装置
Jirik et al. High-resolution ultrasonic imaging using two-dimensional homomorphic filtering
US20230355213A1 (en) Ultrasound image processing
JP2763126B2 (ja) カラー超音波診断装置
JP2006212054A (ja) 超音波観測装置、画像処理装置、及び、画像処理プログラム
Biswas et al. An FPGA-based architecture of DSC–SRI units specially for motion blind ultrasound systems
WO2018168066A1 (ja) 超音波診断装置及びプログラム
JP2002186615A (ja) 超音波診断装置
JP5085144B2 (ja) 超音波診断装置
Biswas et al. DSC-SRI ALGORITHMS FOR MOTION BLIND DETECTION OF OBJECTS IN AN ULTRASOUND SYSTEM
JP2021074207A (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.