CN104159520B - 超声波诊断装置 - Google Patents

超声波诊断装置 Download PDF

Info

Publication number
CN104159520B
CN104159520B CN201380012073.6A CN201380012073A CN104159520B CN 104159520 B CN104159520 B CN 104159520B CN 201380012073 A CN201380012073 A CN 201380012073A CN 104159520 B CN104159520 B CN 104159520B
Authority
CN
China
Prior art keywords
frame
filter
frame data
processing portion
depth
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.)
Expired - Fee Related
Application number
CN201380012073.6A
Other languages
English (en)
Other versions
CN104159520A (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.)
Hitachi Ltd
Original Assignee
Hitachi Aloka Medical Ltd
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
Priority claimed from JP2012045282A external-priority patent/JP5636384B2/ja
Priority claimed from JP2012045280A external-priority patent/JP5367106B2/ja
Application filed by Hitachi Aloka Medical Ltd filed Critical Hitachi Aloka Medical Ltd
Publication of CN104159520A publication Critical patent/CN104159520A/zh
Application granted granted Critical
Publication of CN104159520B publication Critical patent/CN104159520B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • 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/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo 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
    • 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/52077Details 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details

Landscapes

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

Abstract

在超声波图像中,在浅表区域到包含心脏的深区域有很多静止回声。本发明中,HPF处理部(20)通过将已经根据帧内的深度而设定其特性的高通滤波器而应用在所述深度处的帧数据上来对帧数据滤波。例如,可以用数字滤波器实现高通滤波器,并且通过滤波设定部(22)来调整高通滤波器的特性,所述滤波设定部设定所述数字滤波器的滤波系数。也就是说,HPF处理部(20)中的滤波系数是由该滤波设定部(22)控制以使得区域越深,设定越高的偏移级。

Description

超声波诊断装置
技术领域
本发明涉及一种超声波诊断装置,更具体地涉及提高超声波图像的图像质量的技术。
背景技术
通常广泛用于器官、血管和胎儿的诊断的超声波诊断装置已经在例如在心脏状况的诊断中发挥了重要的作用。使用超声波诊断装置,例如使得能够在超声波运动图像上观察跳动的心脏。
在心脏状况的超声波诊断期间,由于胸壁、肋骨、心包等相对于心脏存在于浅表侧(朝向探头),在超声波图像中生成固定回声(静止回声)。这些静止回声导致针对心脏的超声波诊断中的超声波图像的质量下降。因此,通常已经提出了减少这些静止回声的方法。
例如,专利文件1提出的技术是对通过发射和接收超声波而获得的帧数据,沿多个帧的排列方向应用高通滤波器(ahigh-passfilter,HPF),也就是在帧之间应用HPF处理(帧间HPF处理,inter-frameHPFprocessing),由此来减少其关于多个帧的排列方向(时间方向)的波动小于心脏波动的静止回声。
然而,通过简单的帧间HPF处理,由于仅通过完全地移除静止回声难以完全地维持来自心脏的回声,所以部分地生成残留的静止回声,或发生来自心脏的回声的损失,结果在运动图像中产生闪烁。进一步地,相应于在帧间HPF处理中的多个帧之间的加权求和的处理所导致的问题是高速运动的心瓣等在运动图像中双重或三重地出现,并且还导致的问题是心脏(心肌)的界限变得模糊。
例如,专利文件2提出的技术是,为了减少上述的闪烁,在帧间HPF处理后进一步进行帧间LPF(low-passfilter,低通滤波)处理。也就是,通过LPF处理执行帧之间的平滑化(平整化),由此减少闪烁。然而,该技术遭受的问题是,当对整个图像简单地进行帧间LPF处理时,使整个图像模糊。
例如,专利文件3提出的技术是,增加在存在大量静止回声的浅表区域里的HPF处理后的帧数据的权重并增加在心脏存在的深区域里的HPF处理前的帧数据的权重,以使加权求和被应用到HPF之前和之后的帧数据,来抑制上述的闪烁和模糊。通过该技术,减少浅表区域中的静止回声,同时在深区域中抑制与心脏有关的闪烁和模糊。然而,该技术还导致有关电路规模的问题,例如需要应用加权求和的电路结构的问题。
现有技术文件
专利文件
专利文件1:JPH8-107896A
专利文件2:JP2000-139909A
专利文件3:JP2005-288021A
发明内容
技术问题
检视前述的背景技术,本发明的发明人经过多年研究开发了提高包含诸如心脏的待诊断的运动对象的超声波图像的质量的技术,并特别关注于沿着多个帧的排列方向的滤波处理。
本发明产生于这样的研发过程并且目的在于实现关于沿着多个帧的排列方向的滤波处理的改进技术。
解决问题的方案
为了实现上述目标,提供了一种超声波诊断装置,包括:探头,其发射和接收超声波;发射器/接收器部,其控制所述探头从而获得超声波接收信号;滤波处理部,其沿着多个帧的排列方向对基于所述超声波接收信号而获得的帧数据应用滤波处理;以及图像形成部,其基于已经受到所述滤波处理的帧数据而形成超声波图像,其中所述滤波处理部将具有依照帧内的深度设定的偏移级的高通滤波器应用至处于该深度的帧数据,从而执行所述滤波处理。
在上述的超声波诊断装置中,将具有依照帧内的深度而设定的偏移级的高通滤波器应用至处于对应的深度的帧数据。高通滤波器的偏移涉及提高高通滤波器的增益特性,并且偏移级是由在作为高通滤波器的处理目标的频带下限(例如0Hz)的滤波器的增益值定义的。
在数字滤波器中偏移级可相对容易地被调整到期望的级别,例如,通过改变滤波系数等等。因此,根据上述的超声波诊断装置,为了通过沿着多个帧的排列方向的滤波处理来提高超声波图像的质量,能够依照深度来相对容易地调整高通滤波器的特性。
在一个优选的具体实例里,滤波处理部使包括待诊断的运动对象的深区域的帧数据经过具有与浅表区域的帧数据的偏移级相比设定为更高级的偏移级的高通滤波器。在一个优选的具体实例里,滤波处理部应用具有随帧内的深度增加而设定为更高级的偏移级并且具有随帧内的深度增加而更接近全通滤波器特性的特性的高通滤波器。
在一个优选的具体实例里,超声波诊断装置还包括:后处理部,其沿着所述多个帧的的排列方向将低通滤波器或中值滤波器应用至被所述滤波处理部处理过的帧数据。
在一个优选的具体实例里,在应用低通滤波器的情况下,所述后处理部将具有依照帧内的深度而设定的偏移级的低通滤波器应用至处于该深度的帧数据。
在一个优选的具体实例里,在应用中值滤波器的情况下,所述后处理部将中值滤波器唯一地应用至与包括待诊断的运动目标的深区域不同的浅表区域的帧数据。
在一个优选的具体实例里,超声波诊断装置还包括判定部,其确认沿着所述多个帧的排列方向在帧内的各个位置上的被所述滤波处理部处理的帧数据的变化,从而判定在相应位置处是否需要由所述后处理部执行的处理。
在一个优选的具体实例里,所述滤波处理部对帧数据应用具有依照帧数据的帧率和帧数据在帧内的深度而设定的滤波系数的数字高通滤波器。
本发明的有益效果
根据本发明,实现与沿着多个帧的排列方向的滤波处理相关的改进技术。例如,依照本发明的一个优选方案,能够依照该深度相对简单的调整高通滤波器的特性。进一步,根据本发明的另一个优选方案,具有依照帧率而设定滤波系数的数字高通滤波器的应用能够抑制该数字高通滤波器依照帧率的处理结果的波动。
附图说明
图1是示出本发明的实施例中优选的超声波诊断装置的整体结构的示意图。
图2是示出数字滤波器的具体实例的示意图。
图3是示出由3分接头的FIR滤波器形成的HPF的频率特性的示意图。
图4是示出依照深度的滤波系数k0的设定实例的图。
图5是示出由2分接头FIR滤波器形成的HPF的频率特性的示意图。
图6是示出由2分接头FIR滤波器形成的LPF的频率特性的示意图。
图7是示出中值滤波器的具体实例的示意图。
图8是示出判定部的具体实例的示意图。
图9是示出依照帧率的HPF的频率特性的变化的示意图。
图10是示出依照帧率的LPF的频率特性的变化的示意图。
图11是示出依照帧率和深度的滤波系数k0的设定实例(HPF)的图。
图12是示出依照帧率和深度的滤波系数k0的设定实例(LPF)的图。
具体实施方式
图1是示出本发明的实施例中优选的超声波诊断装置的整体结构的示意图(之后称为本超声波诊断装置)。探头10将超声波发射至包括待诊断的组织的诊断区域,并接收从该诊断区域反射的超声波。探头10包括多个用于发射和接收超声波的换能器元件,多个换能器元件通过形成发射波束的发射器/接收器部12进行发射控制。多个换能器元件还接收从该诊断区域反射的超声波。将通过接收超声波所获得的信号输出到发射器/接收器部12,发射器/接收器部12随后形成接收波束。
发射器/接收器部12分别输出与探头10的多个换能器元件对应的发射信号,从而形成超声波的发射波束并且扫描该发射波束。进一步,发射器/接收器部12将整相求和处理(phasealignmentandsummationprocessing)以及其他处理应用至从探头10的相应的多个换能器元件获得的接收信号,由此形成与被扫描的发射波束对应的接收波束,并且输出沿着接收波束获得的回声数据(接收信号)。
发射器/接收器部12在二维平面内扫描超声波束(发射波束和对应的接收波束)并且收集回声数据,从而获得形成对应于该二维平面的帧的帧数据。发射器/接收器部12在该二维平面内重复扫描超声波束从而获得多个帧的帧数据。将如此获得的帧数据存储在例如存储器中,并且将其从存储器读出以例如用于下游侧的HPF处理部20的处理。
本超声波诊断装置适合诊断诸如心脏的运动组织,并且提供有提高包括诸如心脏的待诊断的对象的超声波图像的质量的功能。更具体地,本超声波诊断装置减少与相对于心脏存在于浅表侧(朝向探头10)的胸壁、肋骨、心包等关联的固定回声(静止回声),从而提高超声波图像的质量。
HPF(高通滤波器)处理部20为了减少静止回声而执行对于多个帧所获得的帧数据的滤波处理。HPF处理部20沿着多个帧的排列方向将高通滤波器应用至帧数据,也就是应用帧间HPF处理,从而减少在多个帧的排列方向(时间方向)上的波动小于心脏的波动的静止回声。
在超声波图像上,浅表区域中存在的静止回声多于包含心脏的深区域中存在的静止回声。因此,HPF处理部20通过将具有依照帧内的深度而设定的特性的高通滤波器应用至处于相应深度的帧数据。例如,高通滤波器能够通过数字滤波器来实施,滤波设定部22设定数字滤波器的滤波参数等,从而调整高通滤波器的特性。
图2是示出数字滤波器的具体实例的示意图。帧数据受到沿着多个帧的排列方向的滤波处理。更具体地,关于在帧内某位置(坐标)处的帧数据,在相同位置的帧数据通过多个帧依次受到滤波处理。
图2<A>示出了具有3分接头的FIR滤波器(3分接头FIR滤波器)。多个帧的帧数据被依次输入到该3分接头FIR滤波器。在图2<A>中,能够从三个连续的帧获得的xn-1,xn和xn+1表示在该帧内相同位置(坐标)处的帧数据。
将输入到该3分接头FIR滤波器的一个帧的帧数据存储在帧存储器1中。之后,当输入一个新帧的帧数据时,存储在帧存储器1中的帧数据现在将存储在下一级的帧存储器2中,并将一个帧的新输入的帧数据存储在前一级的帧存储器1中。以此种方式,每次输入一个新帧的帧数据,存储在前一级的帧存储器1中的帧数据被依次地转移到下一级的帧存储器2中。在图2<A>中,xn表示被关注的帧的帧数据;xn-1表示从帧数据xn的帧的前一帧获得的帧数据,并且xn+1表示从帧数据xn的帧的后一帧获得的帧数据。
3分接头FIR滤波器将从三个连续帧获得的三帧数据中的每一个帧数据与滤波系数相乘,并进一步将一起相乘后的三数据项相加。更具体地,在图2<A>中,帧数据xn-1,xn,和xn+1中的每个分别乘以滤波系数k-1,k0,和k1,相乘后作为结果的三数据项被进一步地相加以为被关注的帧输出滤波处理的结果yn=k-1xn-1+k0xn+k1xn+1
图2<B>示出具有2个分接头的FIR滤波器(2分接头FIR滤波器)。在2分接头FIR滤波器的情况下,与3分接头FIR滤波器的情况类似,关于多个帧的帧数据被依次地输入并且将一个帧的帧数据存储在帧存储器中。在图2<B>中,xn表示被关注的帧的帧数据,xn和xn+1表示在该帧中相同位置的帧数据,其能够从两个连续的帧获得。
在2分接头FIR滤波器的情况下,如3分接头FIR滤波器的情况,将帧数据xn和xn+1中的每个与每个滤波系数k0和k1相乘并且在相乘之后作为结果的两数据项被进一步相加以输出关于被关注的帧的滤波处理的结果yn=k0xn+k1xn+1
例如,通过本超声波诊断装置,用图2中示出的3分接头FIR滤波器或2分接头FIR滤波器实现高通滤波器,并且通过调整该滤波系数来调整该高通滤波器的特性。同时图2示出作为具体的数字滤波器的实例的2分接头FIR滤波器和3分接头FIR滤波器,还可以使用具有4个或更多分接头的滤波器。还可以使用IIR滤波器代替FIR滤波器。
关于图2中示出的3分接头FIR滤波器,如果所述滤波系数被设定以满足滤波系数是对称的(k-1=k1),增益特性关于频率单调增加,并且在最高频率处的增益是0[dB]的以下的条件(1-1)和(1-2),则获得展示在图3中示出的特性的高通滤波器(HPF)。
[表达式1]
-0.25≤k1≤0,0.5≤k0≤1(1-1)
|k-1|+|k0|+|k1|=1(1-2)
图3示出由3分接头FIR滤波器(图2)形成的HPF的频率特性。在图3中,横轴显示频率,其中Fs表示帧率,即帧数据的采样频率。纵轴显示滤波器的增益。在图3中,k-1,k0和k1是3分接头FIR滤波器(图2)的滤波系数。
根据在图3中示出的频率特性,HPF的偏移级依照滤波系数而变化。HPF的偏移涉及提高HPF的增益特性。在图3中,将作为待通过HPF处理的频带的下限的0Hz处的滤波器的增益值定义为偏移级。
在图3中,当中心分接头的滤波系数k0是0.5时,偏移级是0且最小,并且在低频侧的抑制效率,即,抑制静止回声的效果为最大。随着滤波系数k0从0.5增加,偏移级也增加并且静止回声的抑制效果逐渐地降低。然后,当滤波系数k0是1.0时,偏移级是1且最大,并且在较低频率侧不能够获得抑制效果,即,不能够获得静止回声的抑制效果,这显示输入信号原样被输出的全通滤波器的滤波特性。
根据本超声波诊断装置,HPF的偏移级依照帧内的帧数据的深度而被设定。特别地,图1中的滤波设定部22控制HPF处理部20中的滤波系数k0,使得深度越深设定的偏移级越高。
图4示出依照深度的滤波系数k0的设定的实例。在图4中,横轴显示深度并且纵轴显示滤波系数k0的值。所述深度涉及帧内的将要被该滤波器处理的帧数据的深度。
在图4示出的设定实例中,在从深度0(接近探头10)到深度d1的浅表部分中,滤波系数k0的值被设定为p(0.5≤p<1.0)。随着深度从d1增加到d2,滤波系数k0的值被增加,并且在深度d2或更深深度的深部分中,滤波系数k0的值被设定为q(p<q≤1.0)。滤波系数k0的值可能在深度d1及其附近和在深度d2及其附近平滑地变化。进一步,滤波系数k0的值可以依需要而在从深度0到深度d1的浅表部分还有在深度d2和更深深度的深部分被改变。还可能适当的调整深度d1和深度d2的位置,以依照过陡的或适中的深度改变滤波系数k0
进一步地,如果期望将重点放在减少浅表区域中的静止回声的效果上,在图4中示出的滤波系数k0的特性能够被转移到整体上更深深度的方向或整体上在更小k0值的方向上。当然,在两个方向的转移能够被一起执行。而且,如果期望将重点放在抑制闪烁和模糊的效果上,在图4中示出的滤波系数k0的特性能够被转移到整体上更浅深度的方向或整体上在更大k0值的方向上。当然,在两个方向上的转移能够被一起执行。
如果2分接头FIR滤波器代替3分接头FIR滤波器使用,则电路结构等能够被进一步简化。例如,关于在图2中示出的2分接头FIR滤波器,如果该滤波系数被设定以满足增益特性关于频率单调增加,并且在最高频率处的增益是0[dB]的以下的条件(2-1)和(2-2),则获得展示在图5中示出的特性的高通滤波器(HPF)。
[表达式2]
0.5≤k0≤1,-0.5≤k1≤0(2-1)
|k0|+|k1|=1(2-2)
图5示出由2分接头FIR滤波器(图2)形成的HPF的频率特性。如在图3中示出的的情况下,在图5中,横轴类似地显示频率,其中Fs表示帧率,即,帧数据的采样频率。纵轴类似地显示滤波器的增益。在图5中,k0和k1是2分接头FIR滤波器(图2)的滤波系数。在图5中示出的频率特性中,如在图3的情况下,所述HPF的偏移级依照滤波系数而改变。
更具体地,在图5中,假设将在作为要通过HPF处理的频带的下限的0Hz的滤波器的增益值定义为偏移级,如图4,当该滤波系数k0是0.5时,偏移级是0且最小,并且在低频率侧的抑制效果,即,抑制静止回声的效果最大。随着该滤波系数k0从0.5增加,偏移级也增加并且静止回声的抑制效果逐渐降低。然后,当滤波系数k0是1.0,偏移级是1且最大,并且在较低频率侧不能够获得抑制效果,即,不能够获得静止回声的抑制效果,这显示全通滤波器的滤波特性。
在2分接头FIR滤波器的情况下,如在3分接头FIR滤波器的情况下,HPF的偏移级依照在帧内的帧数据的深度而被设定。特别地,图1中的滤波设定部22控制滤波系数k0,以使深度越深设定的偏移级越高。对该控制,能够应用用图4描述的滤波系数k0的设定实例。
使用已经参照图2至图5描述的高通滤波器,有可能用相对简单的结构,也就是不需要为例如加权求和处理提供复杂的电路结构(例如,见专利文件3),同时达到静止回声的减少和闪烁及模糊的抑制。
回头参照图1,已经受到HPF处理部20中的帧间HPF处理的多个帧的帧数据首先受到检波处理部30中的检波处理,并且之后被发送到后处理部40。在检波处理部30,执行已知的检波处理。
已经通过本超声波诊断装置,凭借HPF处理部20中的处理而实现了减少静止回声和抑制图像中的闪烁及模糊,有可能通过后处理部40进一步抑制图形中的闪烁。后处理部40沿着多个帧的排列方向将低通滤波器或中值滤波器应用至帧数据。
如果该后处理部40使用低通滤波器(LPF),则能够通过在图2中示出的3分接头FIR滤波器或2分接头FIR滤波器实现所述低通滤波器。当然,可以使用具有4个或更多分接头的滤波器或者使用IIR滤波器代替FIR滤波器。
例如,关于在图2中示出的该2分接头FIR滤波器,如果滤波系数被设定以满足增益特性关于频率单调下降,并且DC分量(频率0Hz分量)的增益是0[dB]的以下条件(3-1)和(3-2),获得展示在图6中示出的特性的低通滤波器(LPF)。
[表达式3]
0.5≤k0≤1(3-1)
k0+k1=1(3-2)
图6示出由2分接头FIR滤波器(图2)形成的LPF的频率特性。在图6中,横轴显示频率,其中Fs表示帧率,即,帧数据的采样频率。纵轴显示滤波器的增益。在图6中,k0和k1是2分接头FIR滤波器(图2)的滤波系数。
在图6示出的频率特性里,偏移级依照滤波系数而变化。LPF的偏移涉及提高LPF的增益特性。在图6中,将在作为要被LPF处理的频带的上限的Fs/2Hz的滤波器的增益值定义为偏移级,也就是,将在采样频率一半的频率的滤波器的增益值定义为偏移级。
在图6中,当该滤波系数k0是0.5时,该偏移级是0且最小,并且在高频率侧的抑制效率,即,平滑效果为最大。随着在多个帧的排列方向上的平滑效果的增加,闪烁的抑制效果也增加。然后,随着滤波系数k0从0.5增加,偏移级增加和平滑效果逐渐下降。然后,当滤波系数k0是1.0,偏移级是1且最大,并且在更高频侧不能够获得抑制效果时,即,不能够获得平滑效果,这显示输入信号被原样输出的全通滤波器的滤波特性。
在本超声波诊断装置,LPF的偏移级依照在帧内的帧数据的深度而设定。具体地,图1中的滤波设定部42控制后处理部40中的滤波系数k0,以使深度越深,设定的偏移级越高。对此控制,能够应用已经使用图4描述的滤波系数k0的设定实例。
如果后处理部40使用中值滤波器,则使用具有在例如图7中示出的结构的中值滤波器。图7示出具体实例的中值滤波器。将多个帧的帧数据依次输入到该中值滤波器。在图7中,xn是被关注的帧数据。进一步,能够从三个连续的帧获得的xn-1,xn和xn+1表示在该帧内的相同位置(坐标)处的帧数据。
将输入到中值滤波器的一个帧的帧数据存储在帧存储器1中。然后,当输入一个新帧的帧数据时,存储在该帧存储器1中的帧数据在现在将存储在下一级的帧存储器2中,并且将一个帧的新输入的帧数据存储在前一级的帧存储器1中。以此种方式,每次输入一个新帧的帧数据,存储在前一级的帧存储器1中的帧数据被依次地转移到下一级的帧存储器2中。在图7中,xn-1表示从帧数据xn的帧的前一帧获得的帧数据,并且xn+1表示从帧数据xn的帧的后一帧获得的帧数据。
中值滤波器选择并输出从三个连续的帧获得的该三个帧数据的中间值。更具体地,在图7中,中值滤波器选择对应于xn-1,xn和xn+1的中间值的数据,并将所选择的数据作为被关注的帧中的处理结果yn而输出。
通过中值滤波器,能够在抑制闪烁的同时抑制图像中的模糊。这里,比将中值滤波器应用至整个图像更期望的是将中值滤波器局部地应用至帧中的浅表部分。更具体地,可期望是通过不将中值滤波器应用至心脏存在的深区域而给出防止在心瓣等的n图像中模糊的产生的考虑。例如,有可能设定深度的阈值并防止将中值滤波器应用至比阈值更深的位置处的帧数据。
进一步,图1中的判定部50可以判定后处理部40中的处理是否必要。在此情况下,关于在HPF处理部20中已经被处理后通过检波处理部30获得的帧数据,判定部50确认在该帧内在每个位置上在多个帧的排列方向上的帧数据的变化,由此判定在后处理部40的处理在该位置上是否必要。
图8示出判定部50的具体实例。多个帧的帧数据被依次输入到判定部50。在图8中,xn-1,xn和xn+1表示能够从三个连续的帧获得的在帧内相同位置(坐标)处的帧数据。
将输入到判定部50的一个帧的帧数据存储在帧存储器1中。然后,当输入一个新帧的帧数据时,存储在帧存储器1的帧数据现在将存储在下一级的帧存储器2中,并且将新输入的一个帧的帧数据存储在前一级的帧存储器1中。以此种方式,每次输入一个新帧的帧数据,存储在前一级的帧存储器1中的帧数据在被依次地转移到下一级的帧存储器2中。在图8中,xn是一个被关注的帧的帧数据,xn-1表示从帧数据xn的帧的前一帧获得的帧数据,并且xn+1表示从帧数据xn的帧的后一帧获得的帧数据。
平均值计算部使用下面的等式计算被关注的帧的之前和之后的帧数据的绝对值的平均值An,等式中M表示被关注的帧之前和之后的帧的数量,在图8的具体实例里M=1。
[表达式4]
A n = 1 2 M &Sigma; m = 1 M ( | x n - m | + | x n + m | ) - - - ( 4 )
判定处理部将被关注的帧的帧数据的绝对值与在平均值计算部计算的该平均值An进行比较,并基于下面的表达式中的条件,为被关注的帧的帧数据判定是否需要进行抑制闪烁的处理。
[表达式5]
||xn|-An|≥TH(5-1)
||xn|-An|/An≥TH(5-2)
在条件(5-1)和(5-2)中,TH是用于判定闪烁的阈值。例如,如果条件(5-1)成立,则判定为被关注的帧的帧数据xn需要闪烁抑制处理。进一步,可以使用被所述平均值An标准化的条件(5-2)代替条件(5-1)。
返回参考图1,后处理部40对在判定部50中被判定为需要进行闪烁抑制处理的被关注的帧的帧数据执行闪烁抑制处理,并且不对被判定为不需要进行闪烁抑制处理的被关注的帧的帧数据执行闪烁抑制处理。当执行闪烁抑制处理时,例如,应用已经参照图6描述的LPF或已经参照图7描述的中值滤波器。进一步,作为对被关注的帧的帧数据的闪烁抑制处理的结果,在图1中的后处理部40可以输出通过表达式4获得的平均值An
对数压缩处理部60对由后处理部40获得的多个帧的帧数据执行已知的对数压缩处理。已经受到对数压缩处理的多个帧的帧数据随后被发送到显示部70,并且将作为例如表现跳动的心脏的运动图像的超声波图像显示在显示部70上。该后处理部40可以被设置在对数压缩处理部60的下游。
如上所述,根据本超声波诊断装置,能够在HPF处理部20中实现静止回声的减少和图像及闪烁与模糊的抑制以及在后处理部40进一步抑制图像的闪烁。当执行该处理时,HPF处理部20使用依照帧中的深度设定其特性的高通滤波器,并且后处理部40使用依照帧内的深度设定其特性的低通滤波器。换句话说,依照帧内的深度来调整该高通滤波器和该低通滤波器的特性。
当高通滤波器和低通滤波器通过数字滤波器来实现时,虽然通过调整滤波系数等的设置而确定数字滤波器的滤波特性,但是数字滤波器的特性还取决于待处理的数据的采样频率。在沿着多个帧的排列方向的滤波处理中,将待处理的帧数据的采样频率依据帧率而变化。因此,当数字滤波器被用在沿着多个帧的排列方向的滤波处理中时,有必要考虑数字滤波器的特性对帧率的依赖性。
本超声波诊断装置应用依照帧率设定的滤波系数的数字滤波器,由此依照帧率来抑制滤波处理的结果的波动。
图9示出依照帧率的HPF的频率特性的变化,并具体地显示了由3分接头FIR滤波器(图2)形成的HPF的频率特性依照帧率而变化。在图9中,横轴显示频率并且纵轴显示滤波器的增益。在图中,Fs代表帧率(帧频率),即,帧数据的采样频率。进一步,在图2中的3分接头FIR滤波器的滤波系数k-1,k0和k1分别被设定为-0.25,0.50和-0.25时,在每个帧率Fs下获得图9所示的频率特性。
另一方面,图10示出依照帧率的LPF的频谱特性的变化,并且具体地显示了由3分接头FIR滤波器(图2)形成的LPF的频率特性依照帧率而变化。在图10中,横轴表示频率并且纵轴表示滤波器的增益。在图中,Fs代表帧率(帧频率),即,帧数据的采样频率。进一步,在图2中的前述3分接头FIR滤波器的滤波系数k-1,k0和k1分别被设置为0.25,0.50和0.25时,在每个帧率Fs下获得图10所示的频率特性。
如图9和图10所示,在数字滤波器中,如果滤波系数固定,则频率特性依照帧率Fs的级别而变化。具体地,随着帧率Fs增加(变高),频率特性趋向于朝向频率轴方向扩展,并且截止频率趋向于变高。
根据本超声波诊断装置,例如依照通过使用者操作而设定的帧率Fs的级别来调整例如图2中示出的数字滤波器的滤波系数。
当分接头的数量固定时,例如由图2示出的3分接头FIR滤波器形成的HPF的频率特性如图3所示。图3显示随着滤波系数k0从0.5增加,HPF的偏移级也增加。进一步,随着偏移级增加,该HPF的截止频率降低。由此,有可能通过增加偏移级来降低该HPF的截止频率。
另一方面,如已经参照图9描述的,随着帧率Fs增加,HPF的频率特性趋向于朝向频率轴方向扩展并且截止频率增加。也就是,截止频率趋向于随着帧率Fs增加而增加。
因此,为了抑制随着帧率Fs增加截止频率也增加的趋向,调整滤波系数以使得,随着帧率Fs增加,偏移级增加并且截止频率下降。该调整允许控制来尽可能地防止HPF的截止频率依照帧率Fs的变化而变化。通过滤波设定部22调整用于图1中的HPF处理部20中的HPF的滤波系数。
而且,当分接头的数量固定时,例如,在图2中示出的由2分接头FIR滤波器形成的LPF的频率特性如图6所示,图6显示随着滤波系数k0从0.5增加,LPF的偏移级也增加。进一步,随着偏移级增加,LPF的截止频率增加。由此,有可能通过增加偏移级来增加LPF的截止频率。
另一方面,如已经参照图10描述的,随着帧率Fs增加,LPF的频率特性趋向于朝向频率轴方向扩展并且截止频率增加。也就是,截止频率趋向于随着帧率Fs的增加而增加。
因此,为了抑制随着帧率Fs增加截止频率也增加的趋向,调整滤波系数以使得,随着帧率Fs增加,偏移级下降并且截止频率下降。该调整允许控制来尽可能地防止LPF的截止频率依照帧率Fs的变化而变化。通过滤波设定部42调整用于图1中的后处理部40中的LPF的滤波系数。
进一步,除了所述滤波系数,也可以调整滤波器的分接头数量。例如,如果在30Hz的帧率下的频率特性有效地在图9中的减少静止回声并抑制图像的闪烁和模糊,则即使当帧率为40Hz或20Hz时,滤波系数和分接头数量也被调整来实现接近于在30Hz的帧率下的那些频率特性的频率特性。具体地,例如,如果帧率从30Hz增加到40Hz,调整滤波系数和分接头数量以使在40Hz的截止频率不高(以使截止频率接近在30Hz的截止频率),同时在40Hz的频率特性不缓和(频率特性接近在30Hz的频率特性)。
类似地,如果在30Hz的帧率下的频率特性有效地抑制图10中的闪烁,则即使当帧率为40Hz或20Hz时,滤波系数和分接头数量也被调整以实现接近于在30Hz的帧率下的那些频率特性的频率特性。
进一步,为了在HPF中均建立依照深度对滤波系数的控制和依照帧率对滤波系数的控制,能够考虑图11示出的设定实例。例如,在依照深度的滤波系数k0的特性被设定为参考帧率(例如,30Hz)的情况中,当帧率被增加(例如,到40Hz)时,图11中示出的HPF的滤波系数k0的特性整体上在更浅的方向上或在滤波系数k0的值整体上变大的方向上转移。当然,在两个方向的转移可以一起执行。另一方面,当上述帧率被降低(例如,到20Hz)时,图11中示出的HPF的滤波系数k0的特性整体上在更深的方向上或在滤波系数k0的值整体上变小的方向上被转移。当然,在两个方向上的转移可以一起执行。
进一步,为了在LPF中均建立依照深度对滤波系数的控制和依照帧率对滤波系数的控制,能够考虑图12示出的设定实例。例如,在依照深度的滤波系数k0的特性被设定为参考帧率(例如,30Hz)的情况中,当帧率被增加(例如,到40Hz)时,图12示出的LPF的滤波系数k0的特性整体上在更深的方向上或在滤波系数k0的值整体上变小的方向上被转移。当然,在两个方向的转移可以一起执行。另一方面,当上述帧率被降低(例如,到20Hz),图12示出的LPF的滤波系数k0的特性整体上在更浅的方向上或在滤波系数k0的值整体上变大的方向上被转移。当然,在两个方向上的转移可以一起执行。
凭借上述的结构,能够在HPF处理部20中实现静止回声的减少和图像的闪烁及模糊的抑制,并且能够在后处理部40进一步抑制图像的闪烁抑制,同时抑制与帧率的变化相关的特性的变化。
上面描述的本发明的优选实施例的目的仅仅是说明,并且不限制本发明的范围。本发明可以包括不离开本发明本质的各种改进。以下也是本发明的优选实施例。
(实施例1)一种超声波诊断装置,包括探头,其发射和接收超声波;发射器/接收器部,其控制探头而获得超声波接收信号;滤波处理部,其沿着多个帧的排列方向对基于超声波接收信号而获得的帧数据应用滤波处理;以及图像形成部,其基于已经受到滤波处理的帧数据来形成超声波图像,其中滤波处理部将具有依照帧数据的帧率设定的滤波系数的高通滤波器应用至帧数据,从而执行滤波处理。
(实施例2)根据实施例1的所述的超声波诊断装置,其中滤波处理部应用具有设定以使帧率越高,偏移级就越高的滤波系数的数字高通滤波器。
(实施例3)根据实施例1或2的所述的超声波诊断装置,其中滤波处理部应用随着帧率增加而具有更多数量的分接头的数字高通滤波器。
(实施例4)根据实施例1到3中的任意实施例所述的超声波诊断装置,其中滤波处理部将具有依照帧数据的帧率和帧数据在帧内的深度设定的滤波系数的数字高通滤波器应用至帧数据。
(实施例5)根据实施例1到4中的任意一个所述的超声波诊断装置,还包括后处理部,其对在滤波处理部中已被处理的帧数据应用具有依照沿着多个帧的排列方向的帧率设定的滤波系数的数字低通滤波器。
(实施例6)根据实施例5所述的超声波诊断装置,其中后处理部应用具有设定以使帧率越低,偏移级越高的滤波系数的数字低通滤波器。
(实施例7)根据实施例6所述的超声波诊断装置,其中后处理部应用随着帧率增加而具有更多分接头的数字低通滤波器。
(实施例8)根据实施例5到7中的任意一个实施例所述的超声波诊断装置,其中后处理部对帧数据应用具有依照帧数据的帧率和帧数据在帧内的深度而设定的滤波系数的数字低通滤波器。
附图标记说明
10,探头,12发射器/接收器部,20,HPF处理部,40,后处理部,50判定部。

Claims (14)

1.一种超声波诊断装置,包括:
探头,其发射和接收超声波;
发射器/接收器部,其控制所述探头从而获得超声波接收信号;
滤波处理部,其沿着多个帧的排列方向对基于所述超声波接收信号而获得的帧数据应用滤波处理;以及
图像形成部,其基于已经受到所述滤波处理的帧数据而形成超声波图像,
其中所述滤波处理部将具有依照帧内的深度设定的偏移级的高通滤波器应用至处于该深度的帧数据,从而执行所述滤波处理,并且
其中所述滤波处理部对帧数据应用具有依照在该帧数据的帧率和帧数据在帧内的深度而设定的滤波系数的数字高通滤波器。
2.根据权利要求1所述的超声波诊断装置,其中
所述滤波处理部应用具有与浅表区域的帧数据的偏移级相比,将包括待诊断的运动对象的深区域的帧数据设定为更高级的偏移级的高通滤波器。
3.根据权利要求1所述的超声波诊断装置,其中
所述滤波处理部应用具有随帧内的深度增加而设定为更高级的偏移级并且具有随帧内的深度增加而更接近全通滤波器特性的特性的高通滤波器。
4.根据权利要求2所述的超声波诊断装置,其中
所述滤波处理部应用具有随帧内的深度增加而被设定为更高级的偏移级并且具有随帧内的深度增加而更接近全通滤波器特性的特性的高通滤波器。
5.根据权利要求1所述的超声波诊断装置,还包括:
后处理部,其沿着所述多个帧的的排列方向将低通滤波器或中值滤波器应用至被所述滤波处理部处理过的帧数据。
6.根据权利要求5所述的超声波诊断装置,其中
在应用低通滤波器的情况下,所述后处理部将具有依照帧内的深度而设定的偏移级的低通滤波器应用至处于该深度的帧数据。
7.根据权利要求5所述的超声波诊断装置,其中
在应用中值滤波器的情况下,所述后处理部将中值滤波器唯一地应用至与包括待诊断的运动目标的深区域不同的浅表区域的帧数据。
8.根据权利要求5所述的超声波诊断装置,还包括:
判定部,其确认沿着所述多个帧的排列方向在帧内的各个位置上的被所述滤波处理部处理的帧数据的变化,从而判定在相应位置处是否需要由所述后处理部执行的处理。
9.根据权利要求1所述的超声波诊断装置,其中
所述滤波处理部应用具有被设定以使偏移级随帧率更高而更高的滤波系数的数字高通滤波器。
10.根据权利要求1所述的超声波诊断装置,其中
所述滤波处理部应用具有随帧率更高而分接头数量就更多的数字高通滤波器。
11.根据权利要求9所述的超声波诊断装置,其中
所述滤波处理部应用具有随帧率更高而分接头数量就更多的数字高通滤波器。
12.根据权利要求5所述的超声波诊断装置,其中
在应用低通滤波器的情况下,所述后处理部将具有依照帧数据的帧率和帧数据在帧内的深度而设定的滤波系数的数字低通滤波器应用至该帧数据。
13.根据权利要求12所述的超声波诊断装置,其中
所述后处理部应用具有被设定以使随帧率更低而偏移级就更高的滤波系数的数字低通滤波器。
14.根据权利要求12所述的超声波诊断装置,其中
所述后处理部应用具有随帧率更高而分接头数量就更多的数字低通滤波器。
CN201380012073.6A 2012-03-01 2013-02-19 超声波诊断装置 Expired - Fee Related CN104159520B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
JP2012-045282 2012-03-01
JP2012-045280 2012-03-01
JP2012045282A JP5636384B2 (ja) 2012-03-01 2012-03-01 超音波診断装置
JP2012045280A JP5367106B2 (ja) 2012-03-01 2012-03-01 超音波診断装置
PCT/JP2013/054015 WO2013129185A1 (ja) 2012-03-01 2013-02-19 超音波診断装置

Publications (2)

Publication Number Publication Date
CN104159520A CN104159520A (zh) 2014-11-19
CN104159520B true CN104159520B (zh) 2016-04-27

Family

ID=49082385

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380012073.6A Expired - Fee Related CN104159520B (zh) 2012-03-01 2013-02-19 超声波诊断装置

Country Status (3)

Country Link
US (1) US9514521B2 (zh)
CN (1) CN104159520B (zh)
WO (1) WO2013129185A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7453040B2 (ja) * 2020-04-01 2024-03-19 富士フイルムヘルスケア株式会社 超音波撮像装置、および、画像処理装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5671744A (en) * 1994-08-15 1997-09-30 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus
CN1695561A (zh) * 2004-05-11 2005-11-16 株式会社东芝 超声波诊断装置和超声波诊断方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5479926A (en) * 1995-03-10 1996-01-02 Acuson Corporation Imaging system display processor
JP4313869B2 (ja) * 1998-11-12 2009-08-12 株式会社東芝 超音波診断装置
JP2000166924A (ja) * 1998-12-10 2000-06-20 Toshiba Iyo System Engineering Kk 超音波画像診断装置および超音波画像の形成方法
JP2000316846A (ja) * 1999-05-11 2000-11-21 Hitachi Medical Corp 超音波診断装置
JP2005288021A (ja) * 2004-04-05 2005-10-20 Toshiba Corp 超音波診断装置及びその診断方法
WO2006088094A1 (ja) 2005-02-17 2006-08-24 Matsushita Electric Industrial Co., Ltd. 超音波ドプラ血流計
JP2006271557A (ja) * 2005-03-28 2006-10-12 Shimadzu Corp 超音波診断装置
JP4928801B2 (ja) 2006-02-23 2012-05-09 株式会社東芝 超音波診断装置
JP5248050B2 (ja) * 2007-06-26 2013-07-31 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波撮像装置
JP5367106B2 (ja) 2012-03-01 2013-12-11 日立アロカメディカル株式会社 超音波診断装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5671744A (en) * 1994-08-15 1997-09-30 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus
CN1695561A (zh) * 2004-05-11 2005-11-16 株式会社东芝 超声波诊断装置和超声波诊断方法

Also Published As

Publication number Publication date
US20150016744A1 (en) 2015-01-15
WO2013129185A1 (ja) 2013-09-06
CN104159520A (zh) 2014-11-19
US9514521B2 (en) 2016-12-06

Similar Documents

Publication Publication Date Title
JP5449852B2 (ja) 超音波診断装置
US6454714B1 (en) Ultrasonic harmonic flash suppression
JP5832737B2 (ja) 超音波診断装置及び超音波画像処理装置
US9782146B2 (en) Ultrasonic diagnostic scanner and method for processing ultrasonic signal
WO2015080006A1 (ja) 超音波診断装置
JP2010068987A (ja) 超音波診断装置
WO2019206709A1 (en) Ultrasound imaging system for high resolution wideband harmonic imaging
JP6300728B2 (ja) 医用画像診断装置
CN104159520B (zh) 超声波诊断装置
JP5367106B2 (ja) 超音波診断装置
JP2011101721A (ja) 超音波診断装置
WO2018168066A1 (ja) 超音波診断装置及びプログラム
JP4276532B2 (ja) 超音波診断装置
JP5636384B2 (ja) 超音波診断装置
JP2019025247A (ja) 超音波撮像装置、および、画像処理装置
JP6045866B2 (ja) 超音波画像処理装置
JP7533326B2 (ja) 医用画像生成装置、医用画像生成方法、およびプログラム
JP4105452B2 (ja) 超音波診断装置
JP7447680B2 (ja) 超音波診断装置、超音波診断装置の制御プログラム、及び、超音波診断装置の制御方法
JPH119595A (ja) 超音波診断装置
US20210055398A1 (en) Ultrasound system with improved noise performance by persistence processing
JP5634817B2 (ja) 超音波診断装置
KR101610877B1 (ko) 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법
CN115462827A (zh) 超声波诊断装置以及图像处理方法
CN115886872A (zh) 血流显影方法、装置、数字扫描变换器和超声成像系统

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: 20170216

Address after: Tokyo, Japan, Japan

Patentee after: Hitachi Ltd.

Address before: Tokyo, Japan

Patentee before: Hitachi Aloka Medical Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160427

Termination date: 20190219