[go: up one dir, main page]

JP6000196B2 - 超音波診断装置、音速決定方法およびプログラム - Google Patents

超音波診断装置、音速決定方法およびプログラム Download PDF

Info

Publication number
JP6000196B2
JP6000196B2 JP2013145442A JP2013145442A JP6000196B2 JP 6000196 B2 JP6000196 B2 JP 6000196B2 JP 2013145442 A JP2013145442 A JP 2013145442A JP 2013145442 A JP2013145442 A JP 2013145442A JP 6000196 B2 JP6000196 B2 JP 6000196B2
Authority
JP
Japan
Prior art keywords
sound speed
element data
ultrasonic
transmission
accuracy
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
JP2013145442A
Other languages
English (en)
Other versions
JP2014079568A (ja
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.)
Fujifilm Corp
Original Assignee
Fujifilm 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 Fujifilm Corp filed Critical Fujifilm Corp
Priority to JP2013145442A priority Critical patent/JP6000196B2/ja
Priority to PCT/JP2013/075485 priority patent/WO2014050742A1/ja
Priority to CN201380049587.9A priority patent/CN104703543A/zh
Publication of JP2014079568A publication Critical patent/JP2014079568A/ja
Priority to US14/669,646 priority patent/US10786228B2/en
Application granted granted Critical
Publication of JP6000196B2 publication Critical patent/JP6000196B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/58Testing, adjusting or calibrating the diagnostic device
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Clinical applications
    • A61B8/0825Clinical applications for diagnosis of the breast, e.g. mammography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4245Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4416Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4488Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • 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/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5261Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/18Methods or devices for transmitting, conducting or directing sound
    • G10K11/26Sound-focusing or directing, e.g. scanning
    • G10K11/34Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
    • G10K11/341Circuits therefor
    • G10K11/346Circuits therefor using phase variation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/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

Landscapes

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

Description

本発明は、超音波診断において、検査対象物内の音速決定機能を有する超音波診断装置、検査対象物内の音速の決定方法、および、検査対象物の音速を決定するプログラムに関する。
従来から、医療分野において、超音波画像を利用した超音波画像診断装置等の超音波診断装置が実用化されている。
一般に、この種の超音波診断装置は、複数の素子(超音波トランスデューサ)を内蔵した超音波プローブ(超音波探触子 以下、プローブとも言う)と、このプローブに接続された装置本体とを有している。超音波診断装置では、プローブの複数の素子から所定の焦点(送信焦点)を形成するように被検体(検査対象物)に向けて超音波ビームを送信し、被検体からの超音波エコーをプローブで受信して、受信した超音波エコーの受信信号を装置本体で電気的に処理することにより超音波画像が生成される。
ここで、超音波診断装置では、一回の超音波ビームの送信による被検体からの超音波エコーを、複数の素子で受信する。従って、同一の反射体で反射された超音波エコーでも、各素子の位置に応じて、超音波エコーの受信時間は遅延する。
そのため、超音波診断装置では、各素子で受信した超音波エコーの受信信号をA/D(アナログ/デジタル)変換してデジタルの受信信号(以下、素子データとする)とした後、この素子データを受信フォーカス処理、即ち、受信信号の遅延時間に応じた遅延補正を行い、位相を合わせ整相加算して受信データ(音線信号)を生成し、この受信データによって、超音波画像を生成している。
従来の超音波診断装置では、被検体内における超音波の音速は一定であると仮定し、超音波の音速を、予め設定した或る値に固定して、受信フォーカス処理を行っていた。
しかし、生体内の脂肪層、筋肉層等の組織の違いにより音速は異なるので、被検体内における超音波の音速は一様ではない。また、太った被検者と、やせた被検者とでは、脂肪層や筋肉層の厚さが異なる。すなわち、超音波の音速には、個人差がある。
そのため、超音波の音速を固定していた従来の超音波診断装置では、被検体内の実際の音速と、設定した音速とが異なると、被検体内で反射されて素子に至る超音波エコーの到達時間が、設定した遅延時間とズレてしまう。
その結果、適正な位相整合を行うことができず、受信フォーカスが不適正になってしまい、得られる超音波画像の画質が劣化するという問題が有る。また、得られる超音波画像が、実際の被検体に対して歪んでしまうという問題も有る。
このような問題に対し、超音波診断装置では、被検体内における音速を決定(算出)し、この音速を用いて、受信フォーカス処理を行うことが行われている。
例えば、特許文献1には、被検体における診断領域を、撮影しようとする超音波画像上で分割してなる着目領域を設定して、超音波の送受信を行い、適宜、設定した複数の音速(設定音速)を用いて、得られた素子データに対して受信フォーカス処理を行って、各着目領域において、予め設定した複数の音速(設定音速)毎にフォーカス指標を算出し、算出したフォーカス指標を用いて、各着目領域における音速(環境音速)を決定する超音波診断装置が記載されている。
フォーカス指標としては、コントラストや輝度等が例示され、例えば、フォーカス指標としての輝度が最も高かった設定音速を、着目領域における音速として決定している。
特開2011−92686号公報
特許文献1に記載される超音波診断装置においては、決定した音速に応じて遅延時間や遅延パターンを補正あるいは選択して、受信フォーカス処理を行って、超音波画像を生成する。これにより、被検体内の場所や被検体毎の音速の違いによらず、安定して、適正な超音波画像を出力することを可能にしている。
ここで、この音速の決定方法は、焦点など、超音波エコー(反射波)の波面形状が綺麗な場所では、正確な音速を求めることができる。しかしながら、焦点から離れた位置などでは、正確な音速を求めることができない場合も、多い。
そのため、被検体内(超音波画像)の位置によっては、決定した音速を用いても受信フォーカス処理を適正に行うことができず、高画質な超音波画像が得られない場合がある。
一方、被検体内のあらゆる場所に対応して、正確な音速を決定するためには、深さや方位方向の位置が異なる多数の焦点に対応して、超音波の送信(送受信)を行えばよい。
しかしながら、多数の焦点に対応して超音波の送受信を行うと、音速の更新等に伴って、超音波画像のフレームレートが低下してしまうという問題も有る。
本発明の目的は、このような従来技術の問題点を解決することにあり、超音波診断において、深度方向の多数の焦点に対応して超音波の送受信を行わなくても、検査対象物内のあらゆる場所において、正確に音速を決定することができる超音波診断装置、音速決定方法およびプログラムを提供することにある。
この目的を達成するために、本発明の超音波診断装置は、超音波ビームを用いて検査対象物を検査する超音波診断装置であって、
超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子と、
探触子に、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行わせる送信手段と、
個々の超音波ビームの送信に対応して、複数の素子が出力したアナログ素子信号を受け、所定の処理を施す受信手段と、
受信手段が処理したアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとするAD変換手段と、
複数の第1の素子データから、第1の素子データのいずれかに対応する第2の素子データを生成するデータ処理手段と、
第2の素子データを用いて、検査対象物内における音速を決定する音速決定手段とを有することを特徴とする超音波診断装置を提供する。
このような本発明の超音波診断装置において、送信手段は、中心となる素子の変更、および、超音波ビームの送信方向の変更の少なくとも一方を行って、探触子に複数回の超音波ビームの送信を行わせるのが好ましい。
また、データ処理手段は、中心となる素子が互いに異なる超音波ビームの送信で得られた複数の第1の素子データ、および、送信方向が互いに異なる超音波ビームの送信で得られた複数の第1の素子データの、少なくとも一方を用いて、第2の素子データを生成するのが好ましい。
また、データ処理手段は、送信領域が重なり合う超音波ビームの送信で得られた、複数の第1の素子データから、第2の素子データを生成するのが好ましい。
また、データ処理手段は、複数の第1の素子データを、素子が超音波エコーを受信した受信時間および素子の位置に応じて重ね合わせて、第2の素子データを生成するのが好ましい。
また、決定した音速の精度の判定を行う判定手段を有するのが好ましい。
また、音速の精度が、所定精度未満である場合には、この音速が所定精度未満の位置に対して、送信焦点を変更して、再度、超音波ビームの送信を行い、この位置の音速を再決定するのが好ましい。
また、再決定した音速の精度を判定し、この音速の再決定および判定を、音速の精度が、所定精度以上となるまで繰り返し行うのが好ましい。
また、第1の素子データを記憶する記憶手段、および、検査対象物内における位置が、送信焦点の位置の近傍か否かを判定する位置判定手段を有し、音速の決定に対応する第1の素子データを記憶手段に記憶しておき、音速の精度が所定精度未満である場合には、位置判定手段によって、音速を決定した位置が、送信焦点の近傍か否かを判定し、送信焦点の近傍である場合には、超音波ビームの送信を行わずに、第1の素子データを用いて音速の再決定を行い、送信焦点の近傍ではない場合には、音速が所定精度未満の位置に対応して、送信焦点を変更して、再度、超音波ビームの送信を行い、この位置の音速を再決定するのが好ましい。
また、再決定した音速の精度を判定し、再決定した音速の精度が、所定精度未満である場合には、音速が所定精度未満の位置に対して、送信焦点を変更して、再度、超音波ビームの送信を行い、この位置の音速を再決定し、再度、音速の精度を判定することを、音速の精度が、所定精度以上となるまで繰り返し行うのが好ましい。
また、検査対象物内の複数の位置において音速を決定し、各位置で音速の精度の判定を行うのが好ましい。
また、音速の精度の判定を、音速のバラツキ、音速の標準偏差、音速の最大値と最小値との差、および、音速の平均値の少なくとも1以上を用いて行うのが好ましい。
また、少なくとも1つの超音波画像に対応する全ての第1の素子データを記憶する素子データ記憶手段を有するのが好ましい。
また、第2の素子データを用いて超音波画像を形成するのが好ましい。
さらに、音速決定手段が決定した音速を用いて、超音波画像を形成するための整合加算を行い、超音波画像を形成するのが好ましい。
また、本発明の音速算出方法は、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子を用いる超音波診断において、検査対象物の音速を決定するに際し、
探触子によって、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の素子にアナログ素子信号を出力し、
このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとし、
複数の第1の素子データから、第1の素子データのいずれかに対応する第2の素子データを生成し、
第2の素子データを用いて、検査対象物内における音速を決定することを特徴とする音速決定方法を提供する。
このような本発明の音速算出方法において、中心となる素子の変更、および、超音波ビームの送信方向の変更の少なくとも一方を行って、探触子に複数回の超音波ビームの送信を行わせるのが好ましい。
また、中心となる素子が互いに異なる超音波ビームの送信で得られた複数の第1の素子データ、および、送信方向が互いに異なる超音波ビームの送信で得られた複数の第1の素子データの、少なくとも一方を用いて、第2の素子データを生成するのが好ましい。
また、送信領域が重なり合う超音波ビームの送信で得られた、複数の第1の素子データから、第2の素子データを生成するのが好ましい。
また、複数の第1の素子データを、素子が超音波エコーを受信した受信時間および素子の位置に応じて重ね合わせて、第2の素子データを生成するのが好ましい。
また、決定した音速の精度の判定を行うのが好ましい。
また、音速の精度が、所定精度未満である場合には、音速が所定精度未満の位置に対して、送信焦点を変更して、再度、超音波ビームの送信を行い、この位置の音速の再決定を行うのが好ましい。
また、再決定した音速の精度を判定し、音速の再決定および判定を、音速の精度が、所定精度以上となるまで繰り返し行うのが好ましい。
また、音速の決定に対応する第1の素子データを記憶しておき、音速の精度が所定精度未満である場合には、音速を決定した位置が、送信焦点の近傍か否かを判定し、送信焦点の近傍である場合には、超音波ビームの送信を行わずに、第1の素子データを用いて音速の再決定を行い、送信焦点の近傍ではない場合には、音速が所定精度未満の位置に対応して、送信焦点を変更して、再度、超音波ビームの送信を行い、この位置の音速を再決定するのが好ましい。
また、再決定した音速の精度を判定し、再決定した音速の精度が、所定精度未満である場合には、音速が所定精度未満の位置に対して、送信焦点を変更して、再度、超音波ビームの送信を行い、この位置の音速を再決定し、再度、音速の精度を判定することを、音速の精度が、規定値以上となるまで繰り返し行うのが好ましい。
また、検査対象物内の複数の位置において音速を決定し、各位置で音速の精度の判定を行うのが好ましい。
さらに、音速の精度の判定を、音速のバラツキ、音速の標準偏差、音速の最大値と最小値との差、および、音速の平均値の少なくとも1以上を用いて行うのが好ましい。
さらに、本発明のプログラムは、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子に、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、中心となる素子を変更して、複数回、行わせて、かつ、個々の超音波ビームの送信に対応して、複数の素子にアナログ素子信号を出力するステップ、
このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとするステップ、
複数の第1の素子データから、第1の素子データのいずれかに対応する第2の素子データを生成するステップ、
および第2の素子データを用いて、検査対象物内における音速を決定するステップを、コンピュータに実行させるプログラムを提供する。
このような本発明のプログラムにおいて、複数の素子にアナログ素子信号を出力させるステップでは、中心となる素子の変更、および、超音波ビームの送信方向の変更の少なくとも一方を行って、探触子に複数回の超音波ビームの送信を行わせるのが好ましい。
また、第2の素子データを生成するステップでは、中心となる素子が互いに異なる超音波ビームの送信で得られた複数の第1の素子データ、および、送信方向が互いに異なる超音波ビームの送信で得られた複数の第1の素子データの、少なくとも一方を用いて、第2の素子データを生成するのが好ましい。
また、第2の素子データを生成するステップでは、送信領域が重なり合う超音波ビームの送信で得られた、複数の第1の素子データから、第2の素子データを生成するのが好ましい。
さらに、第2の素子データを生成するステップでは、複数の第1の素子データを、素子が超音波エコーを受信した受信時間および素子の位置に応じて重ね合わせて、第2の素子データを生成するのが好ましい。
このような本発明によれば、超音波診断において、深度方向の多数の焦点に対応して超音波の送受信を行わなくても、検査対象物内のあらゆる場所において、正確に音速を決定することができる。
そのため、本発明によれば、フレームレートを低下することなく、正確な音速で適正な受信フォーカス処理を行った、高画質な超音波画像を生成することができる。
本発明の超音波診断装置の構成の一例を概念的に示すブロック図である。 図1に示す超音波診断装置における受信フォーカス処理の一例を説明するための概念図である。 図1に示す超音波診断装置の素子データ処理部の構成の一例を概念的に示すブロック図である。 (a)および(c)は、それぞれ理想的な超音波ビームによる超音波の送受信を説明するための概念図であり、(b)および(d)は、それぞれの超音波の送受信で得られる素子データを示す概念図である。 (a)および(c)は、それぞれ実際の超音波ビームによる超音波の送受信を説明するための概念図であり、(b)および(d)は、それぞれの超音波の送受信で得られる素子データを示す概念図である。 (a)および(b)は、同じ反射点に対して、互いに異なる中心素子による超音波の送受信を行った場合の音波の経路を説明するための概念図であり、(c)は、複数の素子で得られる素子データを、(d)は、(c)に示す素子データの遅延時間を、それぞれ説明するための概念図である。 (a)、(b)および(c)は真の信号の、(d)、(e)および(f)はゴーストの、それぞれ素子データ、それらの遅延時間および素子データの重ね合わせ状態を説明するための概念図であり、(g)は、複数の素子に対応する素子データの重ね合わせ状態を、(h)は、(g)における素子データの重ね合わせの結果を、それぞれ説明するための概念図である。 図1に示す超音波診断装置の音速決定部の構成の一例を概念的に示すブロック図である。 図1に示す超音波診断装置の音速決定処理の一例を説明するためのフローチャートである。 図9のフローチャートにおける音速決定方法を説明するためのフローチャートである。 図1に示す超音波診断装置の音速決定処理の別の例を説明するためのフローチャートである。
以下、本発明の超音波診断装置、音速決定方法およびプログラムについて、添付の図面に示される好適実施例を基に、詳細に説明する。
図1に、本発明の音速決定方法を実施する、本発明の超音波診断装置の一例を、ブロック図によって概念的に示す。
図1に示すように、超音波診断装置10は、超音波プローブ12と、超音波プローブ12に接続される送信部14および受信部16と、A/D変換部18と、素子データ記憶部20と、素子データ処理部22と、音速決定部23と、画像生成部24と、表示制御部26と、表示部28と、制御部30と、操作部32と、格納部34とを有する。
図示例においては、送信部14、受信部16、A/D変換部18、素子データ記憶部20、素子データ処理部22、音速決定部23、画像生成部24、表示制御部26、表示部28、制御部30、操作部32、および、格納部34が、超音波診断装置10の装置本体を構成する。
超音波プローブ(超音波探触子)12は、通常の超音波診断装置に用いられる、公知の超音波プローブである。
超音波プローブ12(以下、プローブ12とする)は、超音波トランスデューサを一次元的または二次元的に配列してなる振動子アレイ36を有する。
超音波トランスデューサは、検査対象物(以下、被検体という)の超音波画像の撮像の際に、それぞれ送信部14から供給される駆動信号に従って超音波ビームを被検体に送信すると共に、被検体で反射された超音波エコーを受信して、受信した超音波の強さに応じた受信信号を出力する。
各超音波トランスデューサは、例えば、PZT(チタン酸ジルコン酸鉛)に代表される圧電セラミックや、PVDF(ポリフッ化ビニリデン)に代表される高分子圧電素子、PMN−PT(マグネシウムニオブ酸・チタン酸鉛固溶体)に代表される圧電単結晶等からなる圧電体の両端に電極を形成した、振動子で構成される。
このような振動子の電極に、パルス状または連続波状の電圧を印加すると、印加された電圧に応じて圧電体が伸縮し、それぞれの振動子からパルス状または連続波の超音波を発生する。また、各振動子から発生した超音波は、各振動子の駆動の遅延に応じて、設定された焦点に収束して合成されて(すなわち送信フォーカスされて)、超音波ビームが形成される。
また、振動子は、被検体内で反射された超音波エコーが入射することで伸縮し、この伸縮の大きさに応じた電気信号を発生する。この電気信号が、受信信号として、受信部16に出力される。
送信部14は、例えば、複数のパルサを有し、プローブ12の各超音波トランスデューサ(振動子)に、駆動信号を供給する(駆動電圧を印加する)。
送信部14は、制御部30が選択した送信遅延パターンに基づいて、所定数(複数)の超音波トランスデューサが送信する超音波が、設定した焦点に収束する目的とする超音波ビームを形成するように、駆動信号の遅延量(駆動電圧の印加タイミング)を調節する送信フォーカスを行って、駆動信号を超音波トランスデューサに供給する。なお、この送信遅延パターンは、後述する環境音速に応じて補正されたものでもよい。
これにより、プローブ12(振動子アレイ36)から被検体に、目的とする超音波ビームが送信される。
受信部16は、制御部30からの制御信号に応じて、1回の超音波ビームの送信に対応して、所定数(複数)の超音波トランスデューサが出力した受信信号を受け取り、増幅等の所定の処理を施して、A/D変換部18に供給する。
なお、本発明の超音波診断装置10において、超音波の送受信の方法は、基本的に、公知の超音波診断装置と同様である。
従って、1回の超音波の送受信(1本の超音波ビームの送信、および、この送信に対応する超音波エコーの受信)において、超音波を発生する超音波トランスデューサの数(送信開口の数)、および、超音波を受信(受信部16が受信信号を受け取る)する超音波トランスデューサの数(受信開口の数)は、共に、複数であれば、限定は無い。また、1回の送受信において、送信と受信とで、開口数は、同じでも異なってもよい。
また、少なくとも方位方向(アジマス方向(超音波トランスデューサの配列方向))に隣接する超音波ビームで、送信領域が重複していれば、1つの超音波画像を形成するための超音波の送受信の回数(音線数)や、送受信の中心となる超音波トランスデューサ(中心素子)の間隔(すなわち走査線/音線の密度)にも、限定は無い。従って、超音波で走査する領域に対応する全ての超音波トランスデューサを中心素子として超音波の送受信を行ってもよく、2個置きや4個置きなどの所定間隔の超音波トランスデューサを中心素子として超音波の送受信を行ってもよい。
A/D変換部18は、受信部16から供給されたアナログの受信信号を、アナログ/デジタル変換して、デジタルの受信信号である素子データ(第1の素子データ)とする。
A/D変換部18は、A/D変換した素子データを素子データ記憶部20に供給する。
素子データ記憶部20は、A/D変換部18から供給された素子データを、順次、記憶する。また、素子データ記憶部20は、制御部30から入力されるフレームレートに関する情報(例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータ)を、各素子データに関連付けて格納する。
好ましくは、素子データ記憶部20は、少なくとも1つの超音波画像(1フレームの超音波画像)に対応する全ての素子データを記憶し、かつ、少なくとも超音波画像の表示を終了するまでは、表示中および表示前の超音波画像の素子データを消去しない。
素子データ処理部22は、本発明の特徴とする部位で、素子データを重ね合わせて、各素子データに対応する処理済素子データ(第2の素子データ)を生成する部位である。
具体的には、素子データ処理部22は、制御部30による制御に基づいて、素子データ記憶部20に記憶された素子データのうち、中心となる超音波トランスデューサ(中心となる素子(中心素子))が異なり、かつ、超音波ビームの送信領域が重なり合う、所定数(複数)の超音波ビームの送信で得られた素子データを、各超音波トランスデューサが超音波エコーを受信した時間、および、超音波トランスデューサの位置に応じて重ね合わせて、素子データ(後述する注目素子の素子データ)に対応する処理済素子データを生成する。
素子データ処理部22は、生成した処理済素子データを、音速決定部23および画像生成部24に送る。
音速決定部23は、素子データ処理部22が生成した処理済素子データを用いて、被検体内における超音波の音速(環境音速)を決定する部位である。
素子データ処理部22および処理済素子データ、ならびに、音速決定部23および環境音速に関しては、後に詳述する。
画像生成部24は、制御部30による制御に基づいて、素子データ処理部22から供給された処理済素子データから受信データ(音線信号)を生成し、この受信データから超音波画像を生成するものである。
画像生成部24は、整相加算部38、検波処理部40、DSC42、画像処理部44、および、画像メモリ46を有する。
整相加算部38は、素子データ処理部22が生成した処理済素子データを整合加算して受信フォーカス処理を行い、受信データを生成する。
前述のように、プローブ12の振動子アレイ36は、複数の素子(超音波トランスデューサ)を一次元的あるいは二次元的に配列してなるものである。従って、被検体内の1つの反射点との間の距離は、各超音波トランスデューサで異なる。そのため、同じ反射点で反射された超音波エコーであっても、各超音波トランスデューサに超音波エコーが到達する時間が異なる。整相加算部38は、制御部30が選択した受信遅延パターンに応じて、各超音波トランスデューサ毎の超音波エコーの到達時刻の差(遅延時間)に相当する分、各処理済素子データを遅延して、遅延時間を与えた処理済素子データを整合加算することにより、デジタル的に受信フォーカス処理を行い、受信データを生成する。
整相加算部38は、生成した受信データを、検波処理部40に供給する。
ここで、整相加算部38は、音速決定部23によって、被検体内における超音波の音速(環境音速)が決定されて、供給されている場合には、この環境音速を用いて、遅延時間や受信遅延パターン等の補正を行って、受信フォーカス処理を行う。
なお、環境音速が決定されていない場合には、整相加算部38は、受信遅延パターンを用いた公知の方法で、受信フォーカス処理を行う。
図2に、環境音速を用いた受信フォーカス処理の一例を示す。
ここで、図2は、プローブ12が有する複数の超音波トランスデューサが、同図中左右方向に一列に配列されている、リニアプローブの場合である。しかしながら、コンベックスプローブの場合もプローブ形状が違うだけで、考え方は同じでよい。
方位方向における各々の超音波トランスデューサの幅をLとすると、方位方向の中心の超音波トランスデューサから端部に向かってn番目の超音波トランスデューサまでの距離はnLとなる。
同図に示すように、超音波の反射点が、中心の超音波トランスデューサから配列方向に対して垂直な距離(深さ)dの位置にあるとすると、n番目の超音波トランスデューサと反射点との間の距離(長さ)dnは、式(1)により算出される。
n=((nL)2+d21/2 … (1)
従って、環境音速Vaを用いて、超音波エコーが反射点からn番目の超音波トランスデューサに到達(受信)する時間tは、式(2)により算出される。
n=dn/Va=((nL)2+d21/2/Va … (2)
前述のように、超音波トランスデューサと反射点との間の距離は、各超音波トランスデューサ毎に異なる。そのため、この例の場合、同図上部のグラフに示すように、超音波エコーの到達時間tnは、配列方向の端部側の超音波トランスデューサほど、長くなる。
具体的には、n番目の超音波トランスデューサで受信される超音波は、超音波が反射点から中心の超音波トランスデューサで受信されるまでの時間をt1とすると、中心の超音波トランスデューサで受信される超音波に対して、時間Δt=tn−t1だけ遅れる。本例では、この遅延時間Δtが、すなわち、受信遅延パターンである。
整相加算部38は、各々の超音波トランスデューサに対応する受信データについて、上記時間Δtで表される遅延時間を用いて整相加算を行い、受信フォーカス処理を行う。
なお、本発明において、環境音速に応じた受信フォーカス処理は、この方法に限定はされず、公知の方法が、各種、利用可能である。
例えば、制御部30が、環境音速に応じた受信遅延パターンを選択して、これに応じた制御信号を整相加算部38に供給してもよい。あるいは、制御部30が環境音速に応じて受信遅延パターンを補正して、補正した受信遅延パターンに応じた制御信号を整相加算部38に供給してもよい。あるいは、整相加算部38が、環境音速に応じて、制御部30から供給された制御信号を補正して、受信フォーカス処理を行ってもよい。
検波処理部40は、整相加算部38が生成した受信データに対し、超音波の反射位置の深度に応じて距離による減衰の補正を施した後、包絡線検波処理を施すことにより、被検体内における断層の画像情報(輝度画像情報)であるBモード画像データを生成する。
DSC(digital scan converter)42は、検波処理部40で生成されたBモード画像データを、通常のテレビジョン信号の走査方式に対応する画像データに変換(ラスター変換)する。
画像処理部44は、DSC42から入力されるBモード画像データに、階調処理等の各種の必要な画像処理を施して、表示に供するためのBモード画像データとする。画像処理部44は、画像処理済のBモード画像データを、表示のために表示制御部26に出力し、および/または、画像メモリ46に格納する。
画像メモリ46は、画像処理部44が処理したBモード画像データを格納する、公知の記憶手段(記憶媒体)である。画像メモリ46に格納されたBモード画像データは、必要に応じて、表示部28で表示するために表示制御部26に読み出される。
表示制御部26は、画像処理部44によって所定の画像処理が施されたBモード画像データを用いて、表示部28に超音波画像を表示させる。
表示部28は、例えば、LCD等のディスプレイ装置を含んでおり、表示制御部26の制御の下で、超音波画像を表示する。
制御部30は、操作者により操作部32から入力された指令に基づいて超音波診断装置10の各部の制御を行う部位である。
また、制御部30は、操作部32を用いて操作者によって入力された各種の情報を、必要な部位に供給する。例えば、操作部32に、素子データ処理部22および画像生成部24の整相加算部38で用いられる遅延時間算出に必要な情報、ならびに、素子データ処理部22における素子データ処理に必要な情報の入力が行われた場合には、これらの情報を、必要に応じて、送信部14、受信部16、素子データ記憶部20、素子データ処理部22、画像生成部24および表示制御部26等の各部に供給する。
操作部32は、操作者が入力操作を行うためのものであり、キーボード、マウス、トラックボール、タッチパネル等から形成することができる。
また、操作部32は、操作者が、必要に応じて各種の情報を入力するための、入力機能を備えている。例えば、操作部32は、プローブ12(超音波トランスデューサ)の情報、プローブ12(振動子アレイ36)における送信開口および受信開口、重ね合わせる素子データ数や方法などの処理済素子データの生成に関する情報、超音波ビームの焦点位置等を入力するための、入力機能を備えている。
これらは、例えば、撮影部位(診察部位)の選択、画質の選択、撮影する超音波画像の深度の選択等によって、入力される。
格納部34は、制御部30が超音波診断装置10の各部の制御を実行するための動作プログラム、送信遅延パターンおよび受信遅延パターン、処理済素子データの生成に関する情報、さらには、操作部32から入力されたプローブ12の情報、送信開口および受信開口、焦点位置の情報など、制御部30が超音波診断装置の動作や制御を行うための必要な情報等を格納するものである。
格納部34には、ハードディスク、フレキシブルディスク、MO、MT、RAM、CD−ROM、DVD−ROM等の公知の記録媒体を用いることができる。
なお、超音波診断装置10において、素子データ処理部22、音速決定部23、整相加算部38、検波処理部40、DSC42、画像処理部44、および表示制御部26等は、CPUと、CPUに各種の処理を行わせるための動作プログラムから構成される。しかしながら、本発明においては、これらの部位をデジタル回路で構成してもよい。
前述のように、素子データ処理部22は、素子データ記憶部20に記憶された素子データ(未処理素子データ)のうち、中心となる超音波トランスデューサ(中心素子)が異なり、かつ、超音波ビームの送信領域が重なる、所定数(複数)の超音波ビームの送信で得られた素子データを、各超音波トランスデューサが受信した時間および超音波トランスデューサの位置に応じて重ね合わせて、処理済素子データを生成する部位である。
なお、以下の説明では、超音波トランスデューサのことを、単に『素子』とも言う。
図3に、素子データ処理部22の構成を、ブロック図で概念的に示す。
図3に示すように、素子データ処理部22は、遅延時間算出部48と、重ね合わせ処理部49とを有する。
遅延時間算出部48は、操作部32から入力された、もしくは、操作部32から入力されて格納部34に格納されているプローブ12(超音波トランスデューサ(素子))、超音波ビームの焦点位置、プローブ12の送信開口および受信開口などに関する情報を事前に取得しておく。
また、遅延時間算出部48は、超音波ビームを送信(生成)するために超音波を発振する送信開口の素子と、被検体からの超音波エコーを受信する受信開口の素子との幾何学的な位置に基づいて、受信開口の素子で受信される超音波エコーすなわち素子データの遅延時間を算出する。
重ね合わせ処理部49は、操作部32から入力された、もしくは、操作部32から入力されて格納部34に格納されている、重ね合わせる素子データの数および重ね合わせ処理方法等の素子データ処理に関する情報に基づいて、素子データ記憶部20に記憶されている素子データから、重ね合わせを行う素子データ(中心素子が異なり、かつ、送信領域が重なる超音波ビームで得られた素子データ(2以上の対象領域毎に生成された2以上の素子データ))を読み出す。
さらに、重ね合わせ処理部49は、遅延時間算出部48で算出された、それぞれの素子データに対応する遅延時間に基づいて、2以上の素子データを、受信時間上で、即ち時間を合わせて、かつ、受信された探触子の素子の絶対的な位置を合わせて、重ね合わせて、処理済素子データを生成する。
以下に、素子データ処理部22で行う素子データの処理について、詳細に説明する。
初めに、超音波プローブ12において、送信開口すなわち超音波ビームを送信するために超音波を発信する素子(以下、単に送信素子という)から、被検体に超音波ビームを送信し、被検体との間の相互作用によって発生された超音波エコーを、受信開口すなわち超音波エコーの受信を行う素子(以下、単に受信素子という)で受信して、素子データを得る場合において、送信素子からの超音波ビームと受信素子で得られる素子データとの関係について説明する。
一例として、図4(a)に示すように、3つの素子52c〜52eを送信素子として超音波ビームを送信し、7つの素子52a〜52gを受信素子として超音波エコーを受信する。次いで、図4(c)に示すように、1素子分、素子を方位方向に移動(以下、シフトとも言う)して、3つの素子52d〜52fを送信素子として超音波ビームを送信し、7つの素子52b〜52hを受信素子として超音波エコーを受信して、それぞれ、素子データを取得する。
すなわち、図4(a)に示す例では、中心素子(中心となる素子)は、素子52dであり、図4(b)に示す例では、中心素子は素子52eである。
この際において、反射点54を含む検査対象領域に送信する超音波ビーム56が、焦点58で収束して、素子間隔以下に絞れている理想的な場合を考える。
図4(a)のように、反射点54の真上(反射点と焦点とを結ぶ直線上)にある素子52dを中心素子として、送信素子である素子52c〜52eから超音波ビーム56を送信し、受信素子である素子52a〜52gで超音波エコーを受信して素子データを取得すると、超音波ビーム56の焦点58は、中心素子である素子52dと反射点54とを結ぶ一直線上にある。この場合、超音波ビーム56は、反射点54まで送信されるので、反射点54から反射される超音波エコーが生成される。
反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子である素子52a〜52gに受信され、素子52a〜52gによって、図4(b)に示すような素子データ62が得られる。なお、図4(b)において、縦軸は時間で、横軸は図4(a)に一致する方位方向の位置(素子の位置)である(図4(d)も同じ)。
これに対し、図4(c)に示すように、中心素子を1素子分、シフトさせた場合には、反射点54の真上にある素子52dの隣の素子52eが、中心素子となる。
素子52eを中心素子として、送信素子である素子52d〜52fから超音波ビーム56を送信し、受信素子である素子52b〜52hで超音波エコーを受信する。この際に、同様に超音波ビーム56が理想的であれば、超音波ビーム56の送信方向、即ち、中心素子52eと焦点58とを結ぶ直線上に反射点54が存在しない。従って、この超音波ビーム56は、反射点54に送信されない。
そのため、反射点54から反射される超音波エコーは生成されず、受信素子である素子52b〜52hは、超音波エコーを受信しないので、図4(d)に示すように、反射点54からの反射信号は得られないことになる(素子データの信号強度が『0』になる)。
しかしながら、実際の超音波ビームは、図5(a)および(c)に示す超音波ビーム64のように、焦点58で収束した後に拡散するので、素子間隔より幅が広い。
ここで、図4(a)と同様に、図5(a)のように、反射点54の真上にある素子52dを中心素子として、素子52c〜52eを送信素子として超音波ビーム64を送信した場合には、超音波ビーム64が幅広であっても、その焦点58は、素子54dと反射点54とを結ぶ一直線上にある。従って、超音波ビーム64は、反射点54で反射され、超音波エコーが生成される。
その結果、図4(a)の場合と同様に、反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子である素子52a〜52gに受信され、同様に、図5(b)に示すような真の素子データ66が得られる。
次いで、図4と同様、図5(c)に示すように、中心素子を1素子分、シフトして、隣の素子52eを中心素子として、素子52d〜52fを送信素子として超音波ビーム64を送信し、素子52b〜52hを受信素子として超音波エコーを受信する。この場合でも、超音波ビーム64は幅広であるため、その超音波の送信方向、即ち、中心素子である素子52eと焦点58とを結ぶ直線上に反射点54が存在していなくても、超音波ビーム64は、反射点54に送信される(到達する)。
そのため、反射点54から、超音波ビームの送信方向には、本来、存在しない超音波エコー、所謂ゴーストの反射エコーが発生する。この反射点54からのゴーストの反射エコーは、図5(c)に示すように、所定角度に拡がる受信経路60を通って受信素子である素子52b〜52hに受信される。その結果、素子52b〜52hによって、図5(d)に示すようなゴーストの素子データ68が得られることになる。
このようなゴーストの素子データ68は、素子データから生成される超音波画像の精度を低下させる原因となる。
素子データ処理部22は、素子データに対応する遅延時間を遅延時間算出部48で算出し、重ね合わせ処理部49が、2以上の素子データを、この遅延時間および素子の絶対的な位置に応じて重ね合わせることで、真の素子データを強調して、ゴーストの素子データを減衰させた、高精度な素子データである処理済素子データを生成するものである。
前述のように、遅延時間算出部48は、受信素子(受信開口)の各素子で受信される素子データの遅延時間を算出する。
すなわち、図5(c)に示す超音波ビーム64の伝播距離は、超音波ビーム64が中心素子である素子52eから焦点58を経由して反射点54に至る送信経路と、ゴーストの反射エコーが反射点54から受信素子である素子52b〜52hの個々に至る受信経路との和となる。
この図5(c)に示す超音波ビーム64の伝播距離は、図5(a)に示す超音波ビーム64の伝播距離、すなわち、超音波ビーム64が中心素子52dから焦点58を経由して反射点54に至る送信経路と、真の超音波エコーが反射点54から受信素子である素子52a〜52gに至る受信経路との和より長くなる。
そのため、図5(d)に示すようなゴーストの素子データ68は、図5(b)に示すような真の素子データ66に対して、遅延することになる。
素子データ処理部22の遅延時間算出部48において、真の素子データに対するゴーストの素子データの時間差、即ち遅延時間は、音速、送信素子、超音波ビームの焦点、被検体の反射点、および受信素子の幾何学的配置から算出される。
したがって、遅延時間の計算には、プローブ12の形状(素子間隔、リニア、コンベックスなど)、音速、焦点の位置、送信開口、受信開口などの情報が必要である。遅延時間算出部48では、操作部32によって入力された、若しくは格納部34に格納されたこれらの情報を取得して、遅延時間の計算を行う。なお、音速は、固定値(例えば、1540m/sec)を用いてもよく、あるいは、後述する音速決定部23が決定した音速(環境音速)を用いてもよく、あるいは、操作者が入力できるようにしてもよい。
遅延時間は、例えば、送信素子、超音波ビームの焦点、被検体の反射点、および受信素子の幾何学的配置から算出される、送信素子から焦点を経て反射点に至る超音波ビームの送信経路および反射点から受信素子に至る真の反射超音波エコーまたはゴーストの反射信号の受信経路の合計長さ(伝播距離)と、音速によって算出される伝播時間の差から算出することができる。
本発明では、例えば、図6(a)および図6(b)に示すようにして、真の超音波エコーとゴーストの反射エコーの場合の超音波ビームの送信経路および受信経路の長さを求めることができる。なお、図6において、x方向は方位方向で、y方向は深度方向である。
また、図6(a)は、図5(a)と同様の超音波の送受信を行い、図6(b)は、図5(c)と同様の超音波の送受信を行うものである。
真の超音波エコーの場合、図6(a)(図5(a))に示すように、中心素子である素子52dと、焦点58と、反射点54とは、一直線上に位置している(方位方向の位置が一致している)。すなわち、中心素子52dの真下に、焦点58および反射点54が位置している。
従って、中心素子である素子52dの位置をx−yの2次元座標上の座標(x0、0)とすると、焦点58および反射点54のx座標も『x0』となる。以下、この送信における焦点58の位置を座標(x0、df)、反射点54の位置を座標(x0、z)とし、さらに、素子の間隔をLeとする。
この際において、中心素子である素子52dから焦点58を経て反射点54に至る超音波ビームの送信経路61の長さ(送信経路距離)Lta、および、反射点54から素子52dに至る真の反射超音波エコーの受信経路60の長さ(受信経路距離)Lraは、Lta=Lra=zによって算出できる。
したがって、真の超音波エコーの場合、超音波エコーの伝播距離Luaは、Lua=Lta+Lra=2zとなる。
次いで、図6(b)に示すように、送信素子および受信素子をx方向(方位方向)に1素子分ずらして(図中右方向にシフトして)、中心素子を素子52eとして送受信を行う。図5(c)で示したように、この場合には、反射点54で反射されるのは、ゴーストの反射エコーとなる。
反射点54は、素子52dの真下(方位方向の同位置)に位置している。従って、図6(b)に示すように、この送受信では、中心素子である素子52eと、反射点54とのx方向の位置は、1素子分すなわちLeだけ、x方向にずれる。
反射点54とx方向の位置が一致する素子52dの座標が(x0、0)であるので、中心素子である素子52eの座標は(x0+Le、0)、この送信における焦点58の座標は(x0+Le、df)となる。なお、前述のように、反射点54の座標は(x0、z)である。
従って、中心素子である素子52eから焦点58を経て、反射点54に至る超音波ビームの送信経路61の長さ(送信経路距離)Ltbは、Ltb=df+√{(z−df)2+Le2}で算出できる。他方、反射点54から、直下(x方向=方位方向の同位置)の素子52dに至るゴーストの反射信号の受信経路60の長さ(受信経路距離)Lrbは、Lrb=zによって算出できる。
したがって、ゴーストの反射エコーの場合の超音波の伝播距離Lubは、Lub=Ltb+Lrb=df+√{(z−df)2+Le2}+zとなる。
こうして、図6(a)に示す幾何学配置で求めた送信経路61の距離Ltaと受信経路60の距離Lraを合計した超音波の伝播距離Luaを、音速で割った値が、真の超音波エコーの伝播時間となる。また、図6(b)に示す幾何学配置で求めた送信経路61の距離Ltbと受信経路60の距離Lrbを合計した超音波の伝播距離Lubを、音速で割った値がゴーストの反射エコーの伝播時間となる。
遅延時間は、反射点54と中心素子とのx座標が一致している時の真の超音波エコーの伝播時間と、反射点54と中心素子とのx座標を1素子間隔ずつずらしたときのゴーストの反射エコーの伝播時間の差から求められる。
なお、図6(a)および図6(b)の幾何学モデルでは、送信経路61が焦点58を経由したモデルになっているが、本発明はこれに限定されず、例えば、焦点58を経由せずに、直接、反射点54に至る経路であっても良い。
また、図6(a)および図6(b)の幾何学モデルはリニアプローブの場合であるが、これに限らず他のプローブにおいても、プローブの形状から同様の幾何学計算を行うことができる。
例えば、コンベックスプローブの場合、プローブの半径と素子間隔の角度から幾何学モデルを設定して、同じように計算することができる。
また、ステア送信の場合には、送信角度などの情報を考慮した幾何学モデルを用い、送信素子と反射点との位置関係から、真の素子データおよびその周辺のゴーストの素子データの遅延時間を算出することができる。
さらに、幾何学モデルよって遅延時間を算出する方法に限らず、あらかじめ装置の計測条件に合わせて高輝度反射点を計測した計測結果から、計測条件毎に遅延時間を求めておき、その遅延時間を装置内に記憶しておくことで、同じ計測条件の遅延時間を読み出すようにしておいてもよい。
図6(c)に、真の素子データ66およびゴーストの素子データ68を示す。
図6(c)において、方位方向の中央が、真の素子データ66、すなわち、中心素子と反射点54とでx方向の位置が一致している送受信によって得られた素子データ(図示例では、素子52dを中心素子とした素子データ)である。また、中央の両側が、ゴーストの素子データ、すなわち、中心素子と反射点54とでx方向の位置が一致していない送受信によって得られた素子データ(図示例では、素子52cや素子52e等を中心素子とした素子データ)である。
また、図6(d)に、上述の幾何学計算から得られた真の素子データ66に対する、ゴーストの素子データ68の遅延時間の一例を示す。真の素子データ66を中心に、ゴーストの信号の素子データ68は、x方向すなわち方位方向に対称的に時間が遅れることが示されている。
なお、こうして、素子データ処理部22の遅延時間算出部48において算出された遅延時間は、整相加算部38における遅延補正に用いることもできる。
後に詳述するが、本発明においては、或る注目素子を中心素子とする超音波ビームの送信(注目素子の送受信)で得られた素子データに、中心素子が異なり、かつ、超音波ビームの少なくとも一部が重複する超音波ビームの送信で得られた素子データを、超音波エコーの受信時間と素子の位置とを合わせて重ね合わせることで、注目素子の処理済素子データ(第2の素子データ)を生成する(注目素子の素子データを再構築する)。
図6において、反射点54は、注目素子の真下(方位方向の同位置/注目素子と焦点とを結ぶ直線上)に位置する或るサンプリングポイントの位置(素子データの出力位置)を示している。本発明では、注目素子の送受信におけるサンプリングポイントへの送受信経路を真の素子データの送受信経路と見なし、中心素子が異なる超音波の送受信(周辺素子からの送受信)における同じサンプリングポイントへの送受信経路をゴーストの送受信経路と見なして、両送信経路の差から、遅延時間を算出して、この遅延時間を用いて素子データの時間を合わせて、重ね合わせを行う。言い換えれば、注目素子の送受信で得られた素子データを真の素子データ、中心素子が異なる送受信で得られた素子データをゴーストの素子データと仮定して、遅延時間を算出し、素子データの重ね合わせを行う。
本発明では、全てのサンプリングポイント(全ての素子データの出力位置)に対応して、同様の考え方で遅延時間を算出して、素子データの重ね合わせを行い、各素子の処理済素子データを生成する。
ここで、実際には、方位方向(x方向)にサンプリングポント(反射点)の位置をズラしても、受信経路の長さ(受信経路距離Lrb)は変わらない。従って、各注目素子に関しては、深さ方向(y方向)の各サンプリングポイント毎に、中心素子が異なる送受信による素子データとの遅延時間の算出を行えばよい。
また、この重ね合わせ処理においては、真の素子データがどの素子データであるかを知っている必要はない。すなわち、後に図7を用いて詳述するが、この重ね合わせ処理では、注目素子の素子データが真の素子データであれば、自動的に強調されて素子データが残り、ゴーストであれば素子データは打ち消される。すなわち、注目素子の素子データが真の素子データである場合には、遅延時間による処理が一致して信号が強調され、注目素子の素子データがゴーストの素子データである場合には、遅延時間による処理が一致せずに、信号が打ち消される。
次に、本発明の素子データ処理部22の重ね合わせ処理部49においては、こうして遅延時間算出部48において算出された遅延時間を用いて、素子データの重ね合わせ処理を行う。
なお、重ね合わせ処理部49における重ね合わせ処理では、重ね合わせる時の重ね合わせ素子データ数と重ね合わせ処理方法の情報が必要になるが、これらは、予め、操作部32によって入力しておいても良いし、格納部34に格納しておいても良い。
図7(a)〜(h)に、重ね合わせ処理部49で行われる、重ね合わせ処理の一例を示す。なお、図7に示す例は、素子データ数が5つ、重ね合わせ素子データ数が3つの場合である。
図7(a)は、5回の超音波の送受信によって得られた5つの素子データを横に並べて表示している。また、図7(a)は、素子データ毎に、超音波ビームを送信して、超音波エコーを受信した様子を表している。各素子データの横軸は、受信素子を表しており、それぞれの素子データにおいて超音波ビームの送受信における中心素子を中心にして表示している。縦軸は、受信時間を表す。この例では、例えば前記素子52b〜52fなど、中心素子を、1素子ずつ、ずらして、5回の超音波の送受信を行っている。
図7では、中央の素子データにおける中心素子の真下にのみ、1つの反射点が存在している状態を示す。すなわち、5つの素子データのうち、真中の素子データでは、超音波の送受信において、反射点からの真の超音波エコーが受信されている。つまり、真中の素子データは、真の素子データである。
真中の素子データ以外の両側2つの素子データについては、超音波の送受信の中心素子の真下には反射点は存在していない。しかしながら、送信した超音波ビームの広がりによって、真中の素子データの送信素子の真下に存在する反射点に超音波ビームが当たることで生じた反射エコーの素子データ、即ちゴーストの素子データが写り込んでいる。
ゴーストの素子データは、真の素子データから離れるほど、反射点までの超音波の伝播時間が長くなるため、真の素子データよりも受信時間が遅くなる。また、反射点からの超音波エコーが初めに受信される受信素子の位置は、反射点の真上の素子(反射点と方位方向の位置が一致する素子)である。
ここで、図7の各素子データの横軸は、超音波ビームの送信時における中心素子を中心にしている。従って、図7に示す例では、素子データ毎に、この中心素子を1素子ずつずらして送信していることから、各素子データにおいて方位方向の素子の絶対位置は、1素子ずつずれている。つまり、真中の素子データでは、反射点からの反射信号が初めに受信される受信素子は中心素子であるが、両隣の素子データにおいては、真中の素子データよりも1素子ずれており、右側の素子データでは左に1素子ずれ、左側の素子データでは右に1素子ずれている。さらに、両端の素子データでは、真中の素子データよりも2素子ずれており、右端の素子データでは左に2素子ずれ、左端の素子データでは右に2素子ずれている。このように、ゴーストの信号は、真の信号に対して、受信時間が遅れるだけでなく、受信素子の方向に対してもずれを生じている。
図7(b)に、図7(a)に示す5つの素子データのうちの真中の素子データに対する受信時間の遅延時間の一例を示す。
重ね合わせ処理部49では、図7(b)に示す遅延時間を用いて、真中の素子データを注目素子の素子データとした場合に、注目素子の素子データを中心に、重ね合わせる素子データ数分、図示例では3素子データ分だけ遅延時間補正を行うと共に、各素子データを注目素子との素子位置の差(中心素子の位置の差)に応じて、図示例では両側に1素子分だけ方位方向にシフトさせて、即ち位相を合わせて3素子データ分の未処理素子データを重ね合わせ、注目素子の素子データの1つの重ね合わせ処理済素子データとして求める。
すなわち、本例においては、注目素子を中心素子とする超音波の送受信によって得られた素子データ(以下、注目素子の素子データとも言う)に、注目素子の隣の素子を中心素子とする超音波の送受信によって得られた素子データ(以下、隣の素子の素子データとも言う)を重ね合わせて、注目素子の素子データの処理済素子データを生成している。
こうして得られた注目素子の素子データの重ね合わせ処理済素子データを図7(c)に示す。
前述のように、図7(a)に示す注目素子の素子データは、中心素子(すなわち注目素子)の真下に反射点が存在する、真の素子データである。また、注目素子に隣接する素子を中心素子とする送受信によって得られた素子データも、反射点に入射して、反射された超音波エコーのデータである。
従って、注目素子の両側の隣の素子の素子データに遅延時間補正および方位方向のシフトを行って位相合わせを行うと、図7(c)に示すように、隣の素子の素子データと、注目素子の素子データとは、位相が合うので高輝度位置で重なり合う。そのため、これらの素子データを、例えば加算すると素子データ値は大きな値(高輝度値)を示し、例えば、平均して平均値を求めても強調された値(高輝度値)を示す。
これに対し、図7(d)は、図7(a)と同じ素子データであるが、真中の素子データの左隣の素子データを、注目素子の素子データとした場合の一例を示す。すなわち、この例は、真下に反射点が存在しない素子を中心素子とする超音波の送受信の、中心素子を注目素子とした場合の一例を示す。従って、この素子を中心素子とする素子データは、ゴーストの素子データである。
図7(e)は、図7(b)と同じものであり、図7(a)に示す5つの素子データの注目素子の素子データに対する受信時間の遅延時間の一例を示す。即ち、図7(a)と図7(d)は同じ素子データであるので、図7(d)に示す5つの素子データの注目素子の素子データに対する受信時間の遅延時間とも同じである。
重ね合わせ処理部49では、図7(e)(即ち、図7(b)と同じ)に示す遅延時間を用いて、注目素子の素子データを中心に、重ね合わせ素子データ数分、図示例では3素子データ分だけ遅延時間補正を行うと共に、各素子データを注目素子との素子位置の差(中心素子の位置の差)に応じて、図示例では両側に1素子分だけ方位方向にシフトさせて、3つの未処理素子データを重ね合わせ、注目素子の素子データの1つの重ね合わせ処理済素子データとして求める。
こうして得られた注目素子の素子データの重ね合わせ処理済素子データを図7(f)に示す。
図7(d)に示す注目素子の素子データは、ゴーストの素子データである。そのため、注目素子の素子データの両側の隣接素子データの未処理素子データに遅延時間補正および方位方向のシフトを行って位相合わせを行っても、図7(f)に示すように、隣接素子データの各素子データと注目素子の素子データとは、互いに位相が合わないので重なり合わない。このため、これらの3つの素子データを、例えば加算しても、位相が合っていないために、位相が反転している信号などは信号が打ち消しあうため、加算値は大きくならず、例えば、平均して平均値を求めると小さな値を示すことになる。
他の素子データに関しても、注目素子の素子データとして同様の遅延時間補正および方位方向のシフトを行った結果、図示例の5素子データそれぞれについての隣接する3つの素子データの重なり状態を図7(g)に示し、これらに対して、重ね合わせ処理として、例えば、加算処理、若しくは平均処理した結果を図7(h)に示す。
図7(h)に示すように、図7(a)に示す直下に反射点が存在している中心素子を注目素子とした素子データの場合には、真の信号の素子データが高輝度値を持つ重ね合わせ処理済素子データとして求められる。これに対して、その両側の各2素子データの全4素子データでは、ゴーストの素子データは互いに位相が合わない素子データを加算し、または平均化する。そのため、素子データ同士が、互いに打ち消し合うことになるため、ゴーストの重ね合わせ処理済素子データは、その値が真の信号の素子データである高輝度値を持つ重ね合わせ処理済素子データに対して小さくなり、真の素子データに対してゴーストの素子データの影響を低減させることができ、または、その影響を無視できる程、小さくすることができる。
すなわち、或る素子を注目素子として、この注目素子を中心素子とする超音波ビームの送信によって得られた素子データ(注目素子の素子データ)に、中心素子が異なり、かつ、超音波ビームの送信領域が重なり合う超音波の送受信によって得られた素子データを、1以上、時間および方位方向の位置合わせを行って重ね合わせて、注目素子の素子データに対応する処理済素子データを生成することにより(言い換えれば、超音波ビームの少なくとも一部が重なり、かつ、中心素子が異なる送受信による素子データを用いた、注目素子の素子データの再構築(補正)を行うことにより)、真の素子データを高輝度化して、かつ、ゴーストの素子データを小さくできる。
そのため、後述するように、処理済素子データを用いて音速の決定を行う本発明によれば、ゴーストの影響を無くして、送信する音線上の多数の点で焦点を結んだ場合に等しい素子データ、いわば仮想的なマルチ焦点での超音波の送信で得られた素子データ(受信データ(超音波画像データ))を用いて、高精度に被検体内の音速を決定できる。
また、処理済素子データに整相加算や検波処理を行って、受信データを生成して、超音波画像を生成することにより、同様に、ゴーストの影響を無くし、すなわち音線上の全ての点で焦点を結んだのに等しい素子データで超音波画像生成できるので、高輝度で、鮮鋭性に優れた、高画質な超音波画像を生成することができる。
なお、以下の説明では、この処理済素子データの生成を、マルチライン処理とも言う。
本発明において、中心素子とは、送信の開口数(超音波の送信を行う素子数)が奇数の場合には、方位方向の中央の素子である。
他方、開口数が偶数の場合には、方位方向の中央の素子のいずれかを中心素子とし、あるいは、方位方向の真中に素子が有ると仮定して、中心素子とする。すなわち、開口数が偶数の場合には、開口の真中のライン上に焦点が有るものとして計算を行ってもよい。
なお、重ね合わせ処理部49における重ね合わせ処理方法としては、単に、加算するだけでなく、平均値や中央値をとってもよいし、係数を掛け合わせた上で加算してもよい。なお、平均値や中央値を取ることは、素子データレベルでの平均化フィルタやメディアンフィルタを掛けることに相当すると考えられるが、平均化フィルタやメディアンフィルタの代わりに、通常の画像処理で行われる逆フィルタなども適用してもよい。
あるいは、重ね合わせる各素子データ同士を比較し、類似している場合には最大値、類似していない場合には平均値、分布の偏りがある場合には中間値をとるなど、これに限らず、重ね合わせる各素子データの特徴量に基づいて重ね合わせ処理を変えてもよい。
また、注目素子の素子データに重ね合わせる素子データの数は、図示例の2つに限定はされず、1つでもよく、あるいは、3つ以上でもよい。すなわち、注目素子の素子データに重ね合わせる素子データの数は、要求される処理速度(フレームレート等)や画質などに応じて、適宜、設定すればよい。
ここで、注目素子の素子データに重ね合わせる素子データ数は、超音波ビームのビーム幅の広がり程度に合わせた方が望ましい。従って、深さによってビーム幅が変わる場合には、重ね合わせる素子データ数も深さによって変更してもよい。
また、ビーム幅は送信開口数に依存することから、送信開口数に応じて重ね合わせる素子データの数を変更してもよい。あるいは、画像の輝度値などの特徴量に基づいて重ね合わせ素子データ数を変更してもよいし、重ね合わせ素子データ数を複数パターンを変えて作成した画像から最適な重ね合わせ素子データ数を選択してもよい。
なお、以上のマルチライン処理では、中心素子が異なり、かつ、超音波ビームの送信方向が平行(角度が同一)である、複数の超音波ビームの送信によって得られた素子データを重ね合わせることにより、注目素子の素子データの処理済素子データを生成したが、本発明は、これに限定はされない。
例えば、中心素子を同一として、送信方向(角度)が異なる複数の超音波ビームの送信によって得られた素子データを重ね合わせることにより、処理済素子データを生成してもよい。この際において、何れの超音波ビームの送信で得られた素子データの処理済素子データを生成するか(すなわち、どの方向の音線の処理済素子データを生成するか)は、診察部位やプローブの種類等に応じてデフォルトで設定されていてもよく、あるいは、操作者が選択するようにしてもよい。
また、中心素子が異なり、平行な超音波ビームの送信で得られた素子データと、中心素子を同一として、送信方向が異なる超音波ビームの送信で得られた素子データとの両方を用いて、処理済素子データを生成してもよい。
前述のように、素子データ処理部22は、生成した処理済素子データを、画像生成部24(整相加算部38)に送る。また、被検体の音速の決定(更新)を行う際には、素子データ処理部22は、生成した処理済素子データを、音速決定部23あるいはさらに画像生成部24に送る。
処理済素子データを供給された画像生成部24では、前述のように、整相加算部38が処理済素子データを整相加算して受信フォーカス処理を行って受信データを生成し、検波処理部40が、受信データに減衰補正および包絡線検波処理を施すことにより、Bモード画像データを生成する。
画像生成部24では、さらに、DSC48が、Bモード画像データを通常のテレビジョン信号の走査方式に対応する画像データにラスター変換し、画像処理部44で階調処理等の所定の処理を施す。
画像処理部44は、生成したBモード画像データを画像メモリ46に格納し、および/または、表示制御部26に送って、被検体のBモード画像を表示部28が表示する。
他方、音速決定部23は、供給された処理済素子データを用いて、被検体内の超音波の音速を決定(音速を算出)する。
図8に、音速決定部23の構成をブロック図で概念的に示す。
図8に示すように、音速決定部23は、着目領域設定部70、送信フォーカス制御部72、設定音速指定部74、フォーカス指標算出部76、環境音速決定部78および精度判定部80を有する。
着目領域設定部70は、制御部30からの指示に応じて、Bモード画像上(超音波画像上)において着目領域を設定するものである。
音速決定部23において、被検体の音速は、この着目領域毎に決定する。
本実施形態において、着目領域設定部70は、Bモード画像の画面全体を格子状に分割して、その個々を着目領域とする。
この分割の数(格子の数)は、予めデフォルトで設定されていても良く、操作者が方位方向および/または深度方向で任意に設定できるようにしてもよい。分割数がデフォルトで設定されている場合には、画像サイズ毎や観察部位毎に、異なる設定を有してもよい。さらに、予め複数の分割数を設定しておき、操作者が選択できるようにしてもよい。
なお、本発明において、着目領域は、Bモード画像を格子状に分割した各領域に限定はされない。
例えば、受信データ(Bモード画像データ)生成する全ての画素(全ての画素に対応する位置(領域))を、着目領域としてもよい。言い換えれば、前述の画面を分割する態様において、受信データを生成する全画素に対応して、画面を格子状に分割してもよい。
あるいは、画面全体ではなく、予め設定された、もしくは、複数の選択肢から選択された、画面の一部分について格子状に分割して、その個々を着目領域としてもよい。また、画面全体ではなく、操作者が設定したROIに対応して、着目領域を設定してもよい。なお、画面の一部分やROI内で着目領域を設定する場合にも、分割は、前記画面全体と同様に行えばよい。また、画面全体での着目領域の設定と、ROI内での着目領域の設定とを、操作者が選択できるようにしてもよい。
また、分割する形状は格子状に限らず、例えば、コンベックスプローブによる超音波画像のような扇形のBモード画像であれば、分割する形状もこれに合わせて扇形にしても良い。この場合も、上記の各態様は、全て利用可能である。
なお、着目領域は、画像が大きく変動した場合(画像特徴量の変動値が閾値を超えた場合など)、観察倍率の変更や観察深度の変更などの観察条件の変更が行われた場合等に、変更あるいは更新してもよく、着目領域の変更あるいは更新を、操作者が指示できるようにしてもよい。
着目領域設定部70は、さらに、設定した着目領域に対して、音速の決定に対応する超音波の送信(送信フォーカス)を行うための焦点(焦点の位置)を設定する。
焦点は、観察部位、音線数、送受信の開口数、プローブ12の種類等に応じて、予めデフォルトで設定されていてもよく、操作者が選択や入力指示してもよく、デフォルトでの設定と操作者による指示都を選択できるようにしてもよい。
なお、前述のように、素子データの重ね合わせを行った処理済素子データを用いて音速の決定を行う本発明は、仮想的なマルチ焦点による送信を行うことができる。そのため、焦点の位置は、1つの音線(1つの中心素子に対応する同方向の超音波の送受信/1つの走査線)に対して、複数を設定してもよいが、基本的に、1つの音線に対して、焦点は1カ所でいい。従って、本発明によれば、フレームレートを低減することなく、音速の決定を行うことができる。
また、焦点の位置は、全ての音線で同じでもよく、あるいは、焦点が異なる音線が混在してもよい。
着目領域設定部70は、設定した着目領域および焦点(焦点の位置)の情報を、送信フォーカス制御部72に送る。
送信フォーカス制御部72は、着目領域設定部70が設定した着目領域および焦点に応じて送信部14が送信フォーカスを実行するように、制御部30に送信フォーカス指示を行うものである。
設定音速指定部74は、制御部30の制御に基づき、環境音速の決定において、受信データに対して受信フォーカスを実行するための設定音速を指定するものである。
フォーカス指標算出部76は、素子データ処理部22が生成した処理済素子データを用いて、設定音速指定部74が指定した複数の設定音速毎に受信データに対して受信フォーカスを行って、受信データのフォーカス指標を算出するものである。
環境音速決定部78は、複数の設定音速毎のフォーカス指標に基づき、着目領域の環境音速を決定するものである。
精度判定部80は、環境音速決定部78が決定した音速の精度を判定するものである。
精度判定部は、決定した音速が所定の精度以上の精度を有する場合には、決定した音速を画像生成部24(整相加算部38)に供給するように、環境音速決定部78に指示を出す。逆に、決定した音速が所定の精度に満たない場合には、着目領域設定部70に、焦点の再設定を行うように指示を出す。
本発明において、この精度判定部80は、好ましい態様として設けられる。
以下、図9に示すフローチャートを参照して、超音波診断装置10における音速の決定方法(本発明の音速決定方法)について、詳細に説明する。
本発明のプログラムは、超音波診断装置10が有するコンピュータに、以下の音速決定方法を実行させるプログラムである。
超音波診断装置10において、環境音速を決定する際には、まず、前述のように、制御部30からの指示に応じて、着目領域設定部70が、着目領域および焦点を設定する。
なお、本発明において、環境音速を決定するタイミング(環境音速の更新タイミング)には、特に限定はなく、公知の超音波診断装置と同様に行えば良い。例えば、測定開始の指示に応じて1回だけ行ってもよく、画像が大きく変動した場合(画像特徴量の変動値が閾値を超えた場合など)に環境音速の決定を行ってもよく、適宜決定した所定フレーム数毎や所定時間の経過毎に環境音速の決定を行ってもよく、操作者の入力指示に応じて環境音速の決定を行ってもよく、これらの音速決定のタイミングの2以上を、適宜、選択できるようにしてもよい。
いずれのタイミングで環境音速を決定する場合でも、マルチライン処理を行う本発明によれば、基本的に、1つの音線に対して1つの焦点での送信を行えば良いので、従来に比して、音速の決定に起因するフレームレートの低下を避けることができる。
この着目領域の設定に応じて、送信フォーカス制御部72は、設定された着目領域および焦点に対して送信部14が送信フォーカスを実行して目的とする超音波ビームを送信するように、制御部30に送信フォーカス指示を行う。
これに応じて、送信部14が、プローブ12(振動子アレイ36の対応する超音波トランスデューサ(素子))を駆動して被検体に超音波ビームを送信し、被検体で反射された超音波エコーが超音波トランスデューサ(素子)によって受信され、アナログの受信信号が受信部16に出力される。
受信部16は、アナログの受信信号に、増幅等の所定の処理を施して、A/D変換部18に供給する。
A/D変換部18は、受信部16から供給されたアナログの受信信号をA/D変換して、デジタルの受信信号である素子データとする。
素子データは、素子データ記憶部20に記憶される。
素子データが素子データ記憶部に記憶されると、素子データ処理部22が、前述のマルチライン処理を行って、処理済素子データを生成する。
すなわち、前述の図7に示すように、素子データ処理部22は、例えば、注目素子と、その両隣の素子とに対して、注目素子の素子データに対する、両隣の素子の素子データの遅延時間を算出し、隣の素子の素子データの遅延時間補正および方位方向のシフトを行い、注目素子の素子データに、両側の隣の素子の素子データを重ね合わせて、注目素子の処理済素子データを生成する。
素子データ処理部22は、生成した処理済素子データを、音速決定部23(フォーカス指標算出部76)に供給する。なお、素子データ処理部22は、音速決定のために生成した処理済素子データを、画像生成部24にも供給し、画像生成部24は、この処理済素子データを用いて、超音波画像(Bモード画像データ)を生成してもよい。また、音速決定のために生成した処理済素子データ(素子データ)は、音速決定のみに使用してもよい。
音速決定部23は、供給された処理済素子データを用いて、被検体内における超音波の音速を決定する。
図10に、音速決定部23における音速決定方法の一例のフローチャートを示す。なお、本発明において、音速決定部23での音速決定方法は、この方法に限定はされず、超音波診断装置で行われている各種の音速決定方法(音速の算出方法)が利用可能である。
処理済素子データが供給されると、音速決定部23では、必要に応じて処理済素子データを所定の部位に記憶させると共に、まず、設定音速指定部74が、設定音速Vの開始音速Vstと終了音速Vendを設定し、さらに、設定音速Vに開始音速Vstをセットする。
開始音速Vstおよび終了音速Vendを含む設定音速は、予めデフォルトで設定されていてもよく、あるいは、開始音速Vstおよび終了音速Vendのみを操作者が任意に入力して、間の刻み幅(所定ステップ音速量ΔV)のみがデフォルトで設定されていてもよく、あるいは、操作者が任意に入力するようにしてもよい。また、設定音速や、その刻み幅がデフォルトで設定される場合には、観察部位や性別等に応じて、複数種類の設定音速が設定され、操作者が、適宜、選択できるようにしてもよい。
本例においては、一例として、開始音速Vstとして1410m/secが、終了音速Vendとして1570m/secが設定され、それに応じて、所定の刻み幅として、40m/secの間隔で設定音速が設定されたとする。
次いで、フォーカス指標算出部76が、各着目領域に対応して、設定音速指定部74が指定した複数の設定音速毎に処理済素子データに対して受信フォーカスして、受信データのフォーカス指標を算出する。
具体的には、フォーカス指標算出部76は、着目領域における受信データ(超音波画像データ/超音波画像)の積分値、2乗積分値、ピーク値、鮮鋭度(シャープネス)、コントラスト、輝度値、半値幅、周波数スペクトル積分、最大値や直流成分で規格化された周波数スペクトル積分値や2乗積分値、自己相関値等をフォーカス指標として算出する。
次いで、音速決定部23は、設定音速指定部74にて、設定音速Vが終了音速Vendに達したかどうか判定し、設定音速Vが終了音速Vend未満ならば(No)、所定ステップ音速量ΔV、すなわち本例では40m/secを設定音速Vに加算して、着目領域のフォーカス指標を算出する。
このルーチンを繰り返し、設定音速Vが終了音速Vendに達したと判定すると(Yes)、環境音速決定部78にて、複数の設定音速毎のフォーカス指標に基づき、最も高いフォーカス指標の設定音速を着目領域の環境音速とするなどして、着目領域の環境音速を決定する。例えば、超音波画像の輝度をフォーカス指標として、着目領域において、最も高輝度な超音波画像が得られた音速を、その着目領域の環境音速とする。
すなわち、本例における環境音速とは、プローブ12(振動子アレイ36(超音波トランスデューサ))から、或る着目領域までの音速が一定であると仮定した際の、超音波プローブ12と着目領域との間の領域の平均的な音速である。
前述のように、音速決定部23は、設定した全ての着目領域において、このような環境音速の決定を行う。
前述のように、マルチライン処理で生成した処理済素子データは、ゴーストの影響を無くして、送信する音線上の多数の点で焦点を結んだ場合に等しい素子データ、いわば仮想的なマルチ焦点での超音波の送信で得られた素子データである。
そのため、処理済素子データを用いて音線の決定を行う本発明によれば、1音線1焦点の超音波の送信でも、1音線多焦点の超音波の送信を行った場合と同等以上の高い精度で、音速を決定することができる。また、1音線1焦点の超音波の送信で、高精度に音速を決定できるので、音速の決定(音速の更新)に伴うフレームレートの低下も防止できる。
ここで、環境音速決定部78は、着目領域の環境音速を決定したら、決定した環境音速を、精度判定部80に送る。
前述のように、マルチライン処理によって得られる処理済素子データは、1焦点での超音波の送信でも、ゴーストが大幅に減衰された、仮想的なマルチ焦点での超音波の送信で得られた素子データである。そのため、多焦点での超音波送信で決定した音速と同等の、高精度な音速の決定が可能である。
しかしながら、やはり、処理済素子データを用いても、不適性な音速が決定される可能性は、ゼロでは無い(例えば、後述する送信の焦点近傍など)。
これに対して、図示例の超音波診断装置10では、好ましい態様として、精度判定部80を設け、決定した環境音速の精度を判定する。
その結果、環境音速の精度が、所定の精度以上(『OK』 以下、『音速OK』とも言う)の場合には、精度判定部80は、その着目領域の環境音速を、画像生成部24に送るように、環境音速決定部78に指示を出す(エンド)。
これに対し、環境音速が所定の精度未満(『NG』 以下、『音速NG』とも言う)の場合には、精度判定部80は、着目領域設定部70に、音速NGの着目領域に関して、先とは異なる焦点を設定するように、指示を出し、『着目領域・焦点の設定』のステップに戻って、焦点を変更して、再度、超音波の送受信を行って、環境音速の再決定を行い、さらに、環境音速の精度の判定を行う。
超音波診断装置10においては、この環境音速の再決定を、対応する着目領域が音速OKとなるまで、繰り返し、行い、全ての着目領域において、音速OKとする。
これにより、マルチライン処理を利用することで、従来の音速決定よりもフレームレートの低下を避けつつ、同等以上の精度で音速を決定することが可能となる。
環境音速の精度判定方法には限定はなく、公知の判定方法が、各種、利用可能である。
一例として、画像を複数の着目領域からなる所定領域に分割(例えば、プローブがリニア形であれば正方格子状、コンベック型であれば扇形の格子状など)し、所定領域における環境音速の標準偏差を算出して、標準偏差が予め設定した閾値に満たない場合に、この所定領域の着目領域は音速NGと判定する方法が例示される。
また、同様の所定領域において、環境音速の最大値と最小値を検出し、その差が所定の閾値を超えた場合に、この所定領域の着目領域は音速NGと判定する方法も利用可能である。また、同様の所定領域において、環境音速の平均値を算出し、平均値が所定範囲外である場合に、この所定領域の着目領域は音速NGと判定する方法も利用可能である。また、同様の所定領域において、環境音速の頻度分布を算出し、分布のバラツキが所定の閾値を超えた場合に、この所定領域の着目領域は音速NGと判定する方法も利用可能である。
あるいは、同様の所定領域において、環境音速の平均値を算出し、この平均値との差が所定の閾値を超えている着目領域は、音速NGと判定する方法も利用可能である。
さらに、所定領域の環境音速の積分値、2乗積分値、ピーク値、コントラスト、周波数スペクトル積分、最大値や直流成分で規格化された周波数スペクトル積分値や2乗積分値、自己相関値等を用いて、音速OKおよび音速NGを判定してもよい。
なお、このような環境音速の精度判定と、各着目領域の環境音速の決定(再決定)とは、並行して行ってもよいのは、もちろんである。
環境音速の再決定は、再決定を行う着目領域に対応して、焦点を変更した超音波の送信を行えば、先の音速決定と同じ方法で行っても、異なる方法で行ってもよい。
例えば、環境音速が所定精度に至らない着目領域に対応して、焦点を変更した超音波の送信を行って、同様にマルチライン処理を行って得られた処理済素子データを用いて環境音速を再決定してもよい。あるいは、環境音速の再決定を行う着目領域に焦点を合わせた超音波の送信を行って、マルチライン処理を行わない通常の素子データを用いて、音速を決定してもよい。また、フォーカス指針の変更など、環境音速決定の演算方法を変更してもよい。あるいは、これらの方法を選択して、環境音速の再決定を実行するようにしてもよい。
前述のように、マルチライン処理を行うと、仮想的なマルチ焦点のような超音波の送信を行うことができるので、1音線に対して1焦点の超音波の送信でも、多焦点の送信と同等以上の高精度な環境音速の決定を行うことができる。
ここで、前述の遅延時間の算出方法からも分かるように、マルチライン処理では、焦点を理想的な点音源として見なして、素子データの処理を行う。そのため、焦点が理想的な点音源とは見なせない場合には、焦点近傍で、十分な精度の環境音速を決定できない場合も有る。
これに応じて、超音波診断装置10では、図11のフローチャートに示すように、環境音速の精度を判定して、着目領域の環境音速が所定の精度未満(音速NG)の場合には、その着目領域が焦点近傍(焦点を含む)か否かを判定して、この判定結果に応じて、その後の環境音速の再決定の処理方法を変更してもよい。
図11のフローチャートに示す環境音速の決定では、音速決定の開始(スタート)から環境音速の精度判定までは、前述の図9に示す環境音速の決定と同様に処理が行われる。
環境音速の精度判定の結果、音速OKの場合には、先と同様に、精度判定部80は、その着目領域の環境音速を画像生成部24に送るように、環境音速決定部78に指示を出す(エンド)。
これに対し、音速NGの場合には、精度判定部80は、その着目領域が超音波ビーム送信の焦点の近傍か否かを判定する。
なお、本発明において、焦点の近傍とは、一例として、焦点の位置から10mm以内の位置である。
その結果、音速NGの着目領域が、焦点近傍ではない場合(NO)には、先と同様に、着目領域設定部70に、その着目領域に関して異なる焦点を設定するように、指示を出し、『着目領域・焦点の設定』のステップに戻って、焦点を変更して、再度、超音波の送受信を行って、環境音速の再決定を行い、精度の判定を行う。
これに対し、音速NGの着目領域が、焦点近傍である場合には、環境音速決定部78に、素子データ記憶部20から、この着目領域に対応する素子データを読み出し、素子データを用いて、環境音速を決定するように、指示を出す。
前述のように、処理済素子データではなく、通常の素子データを用いて環境音速を決定する場合には、焦点などの反射波の波面形状が綺麗な領域では、正確に環境音速を決定することができる。従って、音速NGの場合に、その着目領域が焦点近傍か否かを判定して、焦点近傍の場合には、素子データを用いて環境音速を再決定することにより、再度の超音波の送受信を行うことなく、高精度に環境音速を再決定できる。
なお、素子データを用いた環境音速の決定は、処理済素子データを用いた環境音速の決定と同様、公知の方法で行えばよい。また、素子データを用いた環境音速の決定は、処理済素子データと同様に行ってもよく、異なる方法で環境音速を決定(算出)してもよい。
図11に示す例においても、全ての着目領域で音速OKとなるまで、音速の再決定および精度判定を繰り返し行う。
また、焦点近傍において再決定した音速がNGである場合には、音速NGの着目領域内で焦点の再設定を行って、超音波の送受信を行って、素子データを用いて音速の再決定を行ってもよく、あるいは、焦点の位置を着目領域から離れた位置に再設定して、マルチライン処理を行った処理済素子データを用いて音速の再決定を行ってもよく、あるいは、いずれかを選択できるようにしてもよい。
あるいは、周辺で音速OKとなった環境音速を用い、補間によって音速NGの着目領域の環境音速を算出してもよく、あるいは、周辺で音速OKとなった環境音速を、そのまま、音速NGの着目領域の環境音速として用いてもよい。
以上、本発明の超音波診断装置、音速決定方法およびプログラムに関して詳細に説明したが、本発明は、上述の例に限定はされず、本発明の要旨を逸脱しない範囲において、各種の改良や変更を行ってもよいのは、もちろんである。
例えば、図示例の超音波診断装置10では、画像生成部24は、好ましい態様として、マルチライン処理によって生成した処理済素子データを用いて超音波画像を生成しているが、本発明は、これに限定はされない。
すなわち、本発明の超音波診断装置は、被検体内における超音波の音速決定を、マルチライン処理によって生成した処理済素子データを用いて行えば、超音波画像の生成は、通常の素子データを用いて行ってもよい。
また、1画像分の素子データを記憶する素子データ記憶部20を有さずに、後述するマルチライン処理を行うために、1つの注目素子に対応して、毎回、必要な回数の超音波の送受信を行うようにしてもよい。
また、図示例の超音波診断装置10では、着目領域設定部70が、画面全体あるいは一部を格子状に分割して、その個々を着目領域とし、各領域の最適音速を個別に求める構成としたが、本発明は、これに限定はされず、画面全体を1つの着目領域としてもよい。すなわち、従来のように、画面全体の最適音速を1つ求める構成として、画面全体の最適音速を求める際に、処理済素子データを用いる構成としてもよい。
10 超音波診断装置
12 (超音波)プローブ
14 送信部
16 受信部
18 A/D変換部
20 素子データ記憶部
22 素子データ処理部
23 音速決定部
24 画像生成部
26 表示制御部
28 表示部
30 制御部
32 操作部
34 格納部
36 振動子アレイ
38 整相加算部
40 検波処理部
42 DSC
44 画像処理部
46 画像メモリ
48 遅延時間算出部
49 重ね合わせ処理部
70 着目領域設定部
72 送信フォーカス制御部
74 設定音速指定部
76 フォーカス指標算出部
78 環境音速決定部
80 精度判定部

Claims (29)

  1. 超音波ビームを用いて検査対象物を検査する超音波診断装置であって、
    前記超音波ビームを送信し、かつ、前記検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子と、
    前記探触子に、複数の前記素子を用い、所定の送信焦点を形成するように前記超音波ビームを送信させることを、複数回、行わせる送信手段と、
    個々の前記超音波ビームの送信に対応して、複数の前記素子が出力したアナログ素子信号を受け、所定の処理を施す受信手段と、
    前記受信手段が処理したアナログ素子信号をA/D変換して、素子位置と深度と強度との関係を示すデジタル素子信号である第1の素子データとするAD変換手段と、
    複数の前記第1の素子データから、前記第1の素子データのいずれかに対応する、素子位置と深度と強度との関係を示すデータである第2の素子データを生成するデータ処理手段と、
    前記第2の素子データを用いて、前記検査対象物内における音速を決定する音速決定手段とを有し、
    前記データ処理手段は、事前に取得された前記探触子、前記検査対象物内の音速、前記超音波ビームの焦点位置、前記送信手段による前記探触子の送信開口および前記受信手段による前記探触子の受信開口に関する少なくとも1つの情報に基づいて、前記第1の素子データの遅延時間を算出する遅延時間算出手段と、
    予め設定された、複数の前記第1の素子データのうち重ね合わせる第1の素子データの数および重ね合わせ処理方法に基づいて前記複数の第1の素子データを取得し、それぞれの前記第1の素子データに対応する前記遅延時間および受信された前記素子の位置に基づいて重ね合わせ、少なくとも1つの第2の素子データを生成する重ね合わせ処理手段とを有することを特徴とする超音波診断装置。
  2. 前記送信手段は、中心となる素子の変更、および、超音波ビームの送信方向の変更の少なくとも一方を行って、前記探触子に前記複数回の超音波ビームの送信を行わせる請求項1に記載の超音波診断装置。
  3. 前記データ処理手段は、中心となる素子が互いに異なる前記超音波ビームの送信で得られた複数の前記第1の素子データ、および、送信方向が互いに異なる前記超音波ビームの送信で得られた複数の前記第1の素子データの、少なくとも一方を用いて、前記第2の素子データを生成する請求項2に記載の超音波診断装置。
  4. 前記データ処理手段は、送信領域が重なり合う前記超音波ビームの送信で得られた、複数の前記第1の素子データから、前記第2の素子データを生成する請求項1〜3のいずれか1項に記載の超音波診断装置。
  5. 前記決定した音速の精度の判定を行う判定手段を有する請求項1〜4のいずれか1項に記載の超音波診断装置。
  6. 前記音速の精度が、所定精度未満である場合には、この音速が所定精度未満の位置に対して、前記送信焦点を変更して、再度、前記超音波ビームの送信を行い、この位置の音速を再決定する請求項5に記載の超音波診断装置。
  7. 前記再決定した音速の精度を判定し、この音速の再決定および判定を、音速の精度が、前記所定精度以上となるまで繰り返し行う請求項6に記載の超音波診断装置。
  8. 前記第1の素子データを記憶する記憶手段、および、検査対象物内における位置が、前記送信焦点の位置の近傍か否かを判定する位置判定手段を有し、
    前記音速の決定に対応する前記第1の素子データを前記記憶手段に記憶しておき、
    前記音速の精度が所定精度未満である場合には、前記位置判定手段によって、前記音速を決定した位置が、前記送信焦点の近傍か否かを判定し、
    前記送信焦点の近傍である場合には、前記超音波ビームの送信を行わずに、前記第1の素子データを用いて音速の再決定を行い、
    前記送信焦点の近傍ではない場合には、音速が所定精度未満の位置に対応して、前記送信焦点を変更して、再度、前記超音波ビームの送信を行い、この位置の音速を再決定する、請求項5に記載の超音波診断装置。
  9. 前記再決定した音速の精度を判定し、再決定した音速の精度が、前記所定精度未満である場合には、音速が所定精度未満の位置に対して、前記送信焦点を変更して、再度、前記超音波ビームの送信を行い、この位置の音速を再決定し、再度、音速の精度を判定することを、音速の精度が、前記所定精度以上となるまで繰り返し行う請求項8に記載の超音波診断装置。
  10. 前記検査対象物内の複数の位置において音速を決定し、各位置で音速の精度の判定を行う請求項5〜9のいずれか1項に記載の超音波診断装置。
  11. 前記音速の精度の判定を、音速のバラツキ、音速の標準偏差、音速の最大値と最小値との差、および、音速の平均値の少なくとも1以上を用いて行う請求項10に記載の超音波診断装置。
  12. 少なくとも1つの超音波画像に対応する全ての前記第1の素子データを記憶する素子データ記憶手段を有する請求項1〜11のいずれか1項に記載の超音波診断装置。
  13. 前記第2の素子データを用いて超音波画像を形成する請求項1〜12のいずれか1項に記載の超音波診断装置。
  14. 前記音速決定手段が決定した音速を用いて、超音波画像を形成するための整合加算を行い、超音波画像を形成する請求項1〜13のいずれか1項に記載の超音波診断装置。
  15. 超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子を用いる超音波診断において、前記検査対象物の音速を決定するに際し、
    前記探触子によって、複数の前記素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行い、かつ、個々の前記超音波ビームの送信に対応して、複数の前記素子がアナログ素子信号を出力し、
    このアナログ素子信号をA/D変換して、素子位置と深度と強度との関係を示すデジタル素子信号である第1の素子データとし、
    事前に取得された前記探触子、前記検査対象物内の音速、前記超音波ビームの焦点位置、前記探触子の送信開口および前記探触子の受信開口に関する少なくとも1つの情報に基づいて、前記第1の素子データの遅延時間を算出し、
    予め設定された、複数の前記第1の素子データのうち重ね合わせる第1の素子データの数および重ね合わせ処理方法に基づいて前記複数の第1の素子データを読み出し、それぞれの前記第1の素子データに対応する前記遅延時間および受信された前記素子の位置に基づいて重ね合わせ、前記第1の素子データのいずれかに対応する少なくとも1つの、素子位置と深度と強度との関係を示すデータである第2の素子データを生成し、
    前記第2の素子データを用いて、前記検査対象物内における音速を決定することを特徴とする音速決定方法。
  16. 中心となる素子の変更、および、超音波ビームの送信方向の変更の少なくとも一方を行って、前記探触子に前記複数回の超音波ビームの送信を行わせる請求項15に記載の音速決定方法。
  17. 中心となる素子が互いに異なる前記超音波ビームの送信で得られた複数の前記第1の素子データ、および、送信方向が互いに異なる前記超音波ビームの送信で得られた複数の前記第1の素子データの、少なくとも一方を用いて、前記第2の素子データを生成する請求項16に記載の音速決定方法。
  18. 送信領域が重なり合う前記超音波ビームの送信で得られた、複数の前記第1の素子データから、前記第2の素子データを生成する請求項15〜17のいずれか1項に記載の音速決定方法。
  19. 前記決定した音速の精度の判定を行う請求項14〜18のいずれか1項に記載の音速決定方法。
  20. 前記音速の精度が、所定精度未満である場合には、この音速が所定精度未満の位置に対して、前記送信焦点を変更して、再度、前記超音波ビームの送信を行い、この位置の音速の再決定を行う請求項19に記載の音速決定方法。
  21. 前記再決定した音速の精度を判定し、この音速の再決定および判定を、音速の精度が、前記所定精度以上となるまで繰り返し行う請求項20に記載の音速決定方法。
  22. 前記音速の決定に対応する前記第1の素子データを記憶しておき、
    前記音速の精度が所定精度未満である場合には、前記音速を決定した位置が、前記送信焦点の近傍か否かを判定し、
    前記送信焦点の近傍である場合には、前記超音波ビームの送信を行わずに、前記第1の素子データを用いて音速の再決定を行い、
    前記送信焦点の近傍ではない場合には、音速が所定精度未満の位置に対応して、前記送信焦点を変更して、再度、前記超音波ビームの送信を行い、この位置の音速を再決定する請求項19に記載の音速決定方法。
  23. 前記再決定した音速の精度を判定し、再決定した音速の精度が、所定精度未満である場合には、音速が所定精度未満の位置に対して、前記送信焦点を変更して、再度、前記超音波ビームの送信を行い、この位置の音速を再決定し、再度、音速の精度を判定することを、音速の精度が、前記所定精度以上となるまで繰り返し行う請求項22に記載の音速決定方法。
  24. 前記検査対象物内の複数の位置において音速を決定し、各位置で音速の精度の判定を行う請求項19〜23のいずれか1項に記載の音速決定方法。
  25. 前記音速の精度の判定を、音速のバラツキ、音速の標準偏差、音速の最大値と最小値との差、および、音速の平均値の少なくとも1以上を用いて行う請求項24に記載の音速決定方法。
  26. 超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子に、複数の前記素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行わせて、かつ、個々の前記超音波ビームの送信に対応して、複数の前記素子にアナログ素子信号を出力させるステップ、
    このアナログ素子信号をA/D変換して、素子位置と深度と強度との関係を示すデジタル素子信号である第1の素子データとするステップ、
    事前に取得された前記探触子、前記検査対象物内の音速、前記超音波ビームの焦点位置、前記探触子の送信開口および前記探触子の受信開口に関する少なくとも1つの情報に基づいて、前記第1の素子データの遅延時間を算出するステップ、
    予め設定された、複数の前記第1の素子データのうち重ね合わせる第1の素子データの数および重ね合わせ処理方法に基づいて前記複数の第1の素子データを読み出し、それぞれの前記第1の素子データに対応する前記遅延時間および受信された前記素子の位置に基づいて重ね合わせ、前記第1の素子データのいずれかに対応する少なくとも1つの、素子位置と深度と強度との関係を示すデータである第2の素子データを生成するステップ、
    および、前記第2の素子データを用いて、前記検査対象物内における音速を決定するステップを、コンピュータに実行させるプログラム。
  27. 前記複数の前記素子にアナログ素子信号を出力させるステップでは、中心となる素子の変更、および、超音波ビームの送信方向の変更の少なくとも一方を行って、前記探触子に前記複数回の超音波ビームの送信を行わせる請求項26に記載のプログラム。
  28. 前記第2の素子データを生成するステップでは、中心となる素子が互いに異なる前記超音波ビームの送信で得られた複数の前記第1の素子データ、および、送信方向が互いに異なる前記超音波ビームの送信で得られた複数の前記第1の素子データの、少なくとも一方を用いて、前記第2の素子データを生成する請求項27に記載のプログラム。
  29. 前記第2の素子データを生成するステップでは、送信領域が重なり合う前記超音波ビームの送信で得られた、複数の前記第1の素子データから、前記第2の素子データを生成する請求項26〜28のいずれか1項に記載のプログラム。
JP2013145442A 2012-09-27 2013-07-11 超音波診断装置、音速決定方法およびプログラム Active JP6000196B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2013145442A JP6000196B2 (ja) 2012-09-27 2013-07-11 超音波診断装置、音速決定方法およびプログラム
PCT/JP2013/075485 WO2014050742A1 (ja) 2012-09-27 2013-09-20 超音波診断装置、音速決定方法およびプログラム
CN201380049587.9A CN104703543A (zh) 2012-09-27 2013-09-20 超声波诊断装置、声速决定方法及程序
US14/669,646 US10786228B2 (en) 2012-09-27 2015-03-26 Ultrasound diagnostic apparatus, sound velocity determining method, and recording medium

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012215258 2012-09-27
JP2012215258 2012-09-27
JP2013145442A JP6000196B2 (ja) 2012-09-27 2013-07-11 超音波診断装置、音速決定方法およびプログラム

Publications (2)

Publication Number Publication Date
JP2014079568A JP2014079568A (ja) 2014-05-08
JP6000196B2 true JP6000196B2 (ja) 2016-09-28

Family

ID=50388146

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013145442A Active JP6000196B2 (ja) 2012-09-27 2013-07-11 超音波診断装置、音速決定方法およびプログラム

Country Status (4)

Country Link
US (1) US10786228B2 (ja)
JP (1) JP6000196B2 (ja)
CN (1) CN104703543A (ja)
WO (1) WO2014050742A1 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6006249B2 (ja) * 2014-03-24 2016-10-12 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
CN107847218B (zh) * 2015-07-09 2020-11-06 奥林巴斯株式会社 超声波观测装置、超声波观测装置的工作方法及存储介质
JP2017046811A (ja) * 2015-08-31 2017-03-09 セイコーエプソン株式会社 超音波デバイス、超音波モジュール、及び超音波測定機
EP3395257A4 (en) * 2015-12-24 2019-08-07 Olympus Corporation ULTRASONIC OBSERVATION DEVICE, METHOD OF OPERATING THE ULTRASONIC OBSERVATION DEVICE AND A PROGRAM FOR OPERATING THE ULTRASONIC OBSERVATION DEVICE
CN106770647B (zh) * 2016-11-17 2019-07-09 四川大学 脂肪含量测量系统及方法
CN107997784B (zh) * 2017-11-29 2021-01-12 声泰特(成都)科技有限公司 一种基于声速自适应修正的超声波束合成方法和系统
FR3099580B1 (fr) 2019-08-02 2024-11-01 Centre Nat Rech Scient Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5570691A (en) * 1994-08-05 1996-11-05 Acuson Corporation Method and apparatus for real-time, concurrent adaptive focusing in an ultrasound beamformer imaging system
JP3645347B2 (ja) * 1996-02-23 2005-05-11 株式会社東芝 超音波診断装置及び遅延時間最適化方法
US6547732B2 (en) * 1998-10-01 2003-04-15 Koninklijke Philips Electronics N.V. Adaptive image processing for spatial compounding
JP2005046193A (ja) * 2003-07-29 2005-02-24 Toshiba Corp 超音波診断装置及び超音波画像データ生成方法
CN100435740C (zh) * 2003-08-14 2008-11-26 松下电器产业株式会社 超声波诊断设备
WO2007133879A2 (en) * 2006-05-12 2007-11-22 Koninklijke Philips Electronics, N.V. Ultrasonic synthetic transmit focusing with motion compensation
JP5179963B2 (ja) * 2007-09-04 2013-04-10 富士フイルム株式会社 超音波診断装置及びその作動方法、並びに、画像処理プログラム
JP5460144B2 (ja) * 2008-08-11 2014-04-02 キヤノン株式会社 超音波受信ビーム成形装置
JP5486257B2 (ja) * 2009-09-28 2014-05-07 富士フイルム株式会社 超音波診断装置及び弾性指標算出方法
JP5389722B2 (ja) 2009-09-30 2014-01-15 富士フイルム株式会社 超音波診断装置及びその作動方法
CN102834060B (zh) * 2010-03-31 2014-12-24 富士胶片株式会社 焦点信息确定方法及装置以及环境声速取得方法及装置
CN103124521A (zh) 2010-09-27 2013-05-29 富士胶片株式会社 超声波图像生成装置、超声波图像生成方法、以及程序

Also Published As

Publication number Publication date
US20150196283A1 (en) 2015-07-16
WO2014050742A1 (ja) 2014-04-03
US10786228B2 (en) 2020-09-29
JP2014079568A (ja) 2014-05-08
CN104703543A (zh) 2015-06-10

Similar Documents

Publication Publication Date Title
JP6000197B2 (ja) 超音波診断装置、超音波画像生成方法およびプログラム
JP6006249B2 (ja) 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
JP5905856B2 (ja) 超音波検査装置
JP6000196B2 (ja) 超音波診断装置、音速決定方法およびプログラム
JP5905808B2 (ja) 超音波検査装置、超音波画像データ生成方法およびプログラム
JP5946427B2 (ja) 超音波検査装置、超音波検査方法、プログラム及び記録媒体
JP6074299B2 (ja) 超音波診断装置、超音波診断装置の信号処理方法およびプログラム
JP6165089B2 (ja) 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
US10980515B2 (en) Acoustic wave processing apparatus, signal processing method, and program for acoustic wave processing apparatus
JP5873412B2 (ja) 超音波診断装置、音速決定方法およびプログラム
JP5964774B2 (ja) 超音波診断装置、超音波診断装置の信号処理方法およびプログラム
US10383601B2 (en) Acoustic wave processing apparatus, signal processing method, and program for acoustic wave processing apparatus
JP6047041B2 (ja) 超音波診断装置、超音波診断装置の信号処理方法およびプログラム
WO2014050847A1 (ja) 超音波診断装置、超音波画像データ生成方法およびプログラム
WO2014050756A1 (ja) 超音波検査装置、超音波画像データ生成方法およびプログラム
WO2014050897A1 (ja) 超音波検査装置、超音波画像データ生成方法およびプログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141023

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151027

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151215

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160607

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160727

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160816

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160830

R150 Certificate of patent or registration of utility model

Ref document number: 6000196

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250