JP2014079565A - 超音波検査装置、超音波検査方法、プログラム及び記録媒体 - Google Patents

超音波検査装置、超音波検査方法、プログラム及び記録媒体 Download PDF

Info

Publication number
JP2014079565A
JP2014079565A JP2013125337A JP2013125337A JP2014079565A JP 2014079565 A JP2014079565 A JP 2014079565A JP 2013125337 A JP2013125337 A JP 2013125337A JP 2013125337 A JP2013125337 A JP 2013125337A JP 2014079565 A JP2014079565 A JP 2014079565A
Authority
JP
Japan
Prior art keywords
sound speed
element data
unit
ultrasonic
sound
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2013125337A
Other languages
English (en)
Other versions
JP5946427B2 (ja
Inventor
Hiroaki Yamamoto
拓明 山本
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 JP2013125337A priority Critical patent/JP5946427B2/ja
Priority to CN201380050604.0A priority patent/CN104684487A/zh
Priority to PCT/JP2013/075666 priority patent/WO2014050797A1/ja
Publication of JP2014079565A publication Critical patent/JP2014079565A/ja
Priority to US14/671,236 priority patent/US10687786B2/en
Application granted granted Critical
Publication of JP5946427B2 publication Critical patent/JP5946427B2/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/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/13Tomography
    • A61B8/14Echo-tomography
    • A61B8/145Echo-tomography characterised by scanning multiple planes
    • 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/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • 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
    • 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • 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
    • 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/8997Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
    • 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/52036Details of receivers using analysis of echo signal for target characterisation
    • 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
    • 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • 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

Landscapes

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

Abstract

【課題】マルチライン処理の計算を含めた上で音速補正を行うことにより、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができる超音波検査装置、超音波検査方法、プログラム及び記録媒体を提供する。
【解決手段】複数の素子を備える探触子と、超音波ビームを焦点を形成するように複数の素子から検査対象物に送信させる送信部と、検査対象物からの超音波エコー信号を複数の素子を介して受信する受信部と、検査対象物内の音速値を決定する音速決定部と、その音速値を用いて受信部によって得られる2つ以上の第1の素子データから第1の素子データのいずれかに対応する第2の素子データを作成する素子データ処理部とを有し、音速決定部は、素子データ処理部において第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として求めることにより、上記課題を解決する。
【選択図】図2

Description

本発明は、超音波ビームを送受信することにより生体内の臓器等の検査対象物の撮像を行って、検査対象物の検査や診断のために用いられる超音波画像を生成する超音波検査装置、超音波検査方法、プログラム及び記録媒体に関する。
従来から、医療分野において、超音波画像を利用した超音波画像診断装置等の超音波検査装置が実用化されている。一般に、この種の超音波検査装置は、複数の素子(超音波トランスデューサ)を内蔵した超音波探触子(超音波プローブ)と、この超音波探触子に接続された装置本体とを有しており、超音波探触子の複数の素子から検査対象物(被検体)に向けて超音波ビームを送信し、被検体からの超音波エコーを超音波探触子で受信して、その受信した超音波エコー信号を装置本体で電気的に処理することにより超音波画像が生成される。
超音波検査装置においては、超音波画像を生成するとき、被検体の検査対象領域、例えば、生体内の臓器やその臓器内の病巣等に探触子の複数の素子から焦点を合わせて超音波ビームを送信し、検査対象領域の反射体、例えば臓器や病巣等の表面や界面からの超音波エコーを複数の素子を介して受信しているが、同一の反射体で反射された超音波エコーを複数の素子で受信するので、送信素子から送信された超音波ビームの焦点位置に位置する反射体で反射され、送信素子で受信された超音波エコー信号に対して、同一の反射体で反射され、送信素子と異なるその他の素子で受信された超音波エコー信号は遅延することになる。このため、超音波検査装置では、複数の素子で受信した超音波エコー信号をA/D(アナログ/デジタル)変換して素子データにとした後、素子データを受信フォーカス処理して、即ち遅延補正して位相を合わせ整相加算して音線信号を生成し、こうして得られた音線信号に基づいて超音波画像を生成している。
このような超音波検査技術において、超音波画像の画質を向上させるために、複数の異なる焦点に収束する複数の超音波ビームの送信で得られた信号を足し合わせることで、従来よりも信号の質を改善することが行われている。
例えば、特許文献1には、送信用振動素子群を構成する複数個の振動素子から放射される送信超音波を送信集束点に集束させることにより仮想的な点音源を形成し、そして、この点音源から放射される送信超音波によって連続した複数の観測点から反射する受信超音波を、受信用振動素子群を構成する複数個の振動素子によって受信し、得られたチャンネル分の受信信号に対し観測点が受信集束点となるような受信整相加算を行う超音波診断装置が開示されている。更に、この超音波診断装置では、受信用振動素子群と振動素子の配列方向に順次シフトさせた送信用振動素子群の各々を用いて得られた受信信号に対しても同様の受信整相加算を行い、これら受信整相加算後の受信信号に対し各々の送信集束点から観測点までの伝搬距離の差異に起因した送信遅延を補正する送信整相加算を行う。
特許文献1に開示の超音波診断装置では、複数の振動素子から得られた受信信号に対して受信整相加算と送信整相加算を行なうことにより、被検体の深さ方向に対してほぼ一様な細いビーム幅を有した送信ビーム及び受信ビームを高精度かつ高感度で形成することができる。このため、特許文献1は、超音波診断装置において、空間分解能、コントラスト分解能及びS/Nに優れた画像データの生成と表示が可能となることを開示している。
特に、特許文献1(その図5に示す実施例参照)には、図8に示すように、3つの送信用振動素子群がそれぞれの送信収束点(以下、焦点とする)Ft1〜Ft3で焦点を結び、各焦点を仮想的な点音源(以下、仮想点音源とする)と見なした場合に、各仮想点音源からの波面が同一の観測点(反射点)Pxの位置で反射することによって得られた反射波を受信用振動素子群で受信することが示されている。特許文献1では、これにより、1つの点を複数の異なる焦点で送信した信号を足し合わせることで、従来よりも信号の質を改善させることができるとしている。
特開2009−240700号公報
しかしながら、特許文献1に開示の技術では、従来技術よりも高画質な画像が得られるが、1ラインのデータを作るのに複数の送信ビームを、送信位置を変えて発生させる必要があり、従来技術よりも送信回数が増えるためフレームレートが低下し、リアルタイム性が悪くなるという問題があった。
また、特許文献1に開示の技術では、生体組織内の音速分布にムラがあると、各々の仮想点音源からの波面にずれが生じるため、精度が悪くなってしまうという問題があった。
本発明の目的は、上記従来技術の問題点を解消し、超音波プローブの複数の素子からの超音波ビームを送信する際に、異なる2以上の素子からの送信で取得した素子データを重ねるマルチライン処理を行うことで、送信超音波ビームの広がりによる影響を低減し、SN比を上げ、解像度を上げることができ、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高い解像度で、最適な空間分解能を持つシャープな超音波画像を得ることができると共に、マルチライン処理の計算を含めた上で音速補正を行うことにより、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができる超音波検査装置、超音波検査方法、プログラム及び記録媒体を提供することにある。
本発明者は、本出願人の先の出願に係る特願2012−158009号明細書において、超音波プローブの複数の素子からの超音波ビームを送受信する際にマルチライン処理を行うことで画質の良い超音波画像を得ることができる超音波検査装置を提案している。上記目的を達成するために、更に、鋭意研究を重ねた結果、先の出願に開示の技術には、音速分布のムラを補正する音速補正について記載されているが、反射点(上記特許文献1の図6の観測点Px)からプローブまでの音速しか補正されておらず、各々の仮想的な点音源を形成する送信素子から反射点(Pxの点)までの音速は補正されていないため、音速の補正が不十分であったことを知見し、本発明に至ったものである。
すなわち、本発明の超音波検査装置は、超音波ビームを用いて検査対象物を検査する超音波検査装置であって、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子と、探触子に、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行わせる送信部と、個々の超音波ビームの送信に対応して、複数の素子が出力したアナログ素子信号を受け、所定の処理を施す受信部と、受信部が処理したアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとするAD変換部と、検査対象物内の音速値を決定する音速決定部と、検査対象物内の音速値を用いて2つ以上の第1の素子データから、2つ以上の第1の素子データのいずれかに対応する第2の素子データを生成する素子データ処理部と、を有し、音速決定部は、素子データ処理部において2つ以上の第1の素子データから第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として求めることを特徴とする。
ここで、音速決定部は、素子データ処理部を有し、素子データ処理部における2つ以上の第1の素子データの遅延時間補正計算に用いる音速値を最適化して最適な音速として求めることが好ましい。
また、音速決定部は、検査対象物内を複数の領域に分割したときの各領域における音速値を決定することが好ましい。
また、更に、2つ以上の第1の素子データを含む一部の第1の素子データ、あるいは全ての第1の素子データを保持する素子データ保持部を有することが好ましい。
また、素子データ保持部は、送信部により、検査対象物内における少なくとも2つ以上の重なり合う対象領域毎において、超音波ビームを送信し、受信部により、超音波ビームによって検査対象物内における少なくとも2つ以上の重なり合う対象領域毎に発生した超音波エコー信号を受信して生成された、素子データ毎に受信した各素子における受信時間毎の受信データを含む2以上の第1の素子データを保持し、素子データ処理部は、素子データ毎の超音波ビームを送信したときの素子の幾何学的な配置及び素子データ毎の受信データの各素子の幾何学的な配置に基づいて、2以上の第1の素子データから第2の素子データを生成することが好ましい。
また、送信部は、超音波ビームの送信を、中心となる素子を変更して、複数回、行うことが好ましい。
また、素子データ処理部は、検査対象物内の音速値を用いて2つ以上の第1の素子データの各遅延時間を算出する遅延時間算出部と、2つ以上の第1の素子データを算出された各遅延時間及び受信された探触子の素子の位置に基づいて重ね合わせ、第2の素子データを生成する重ね合わせ処理部とを有することが好ましい。
また、更に、素子データ処理部で生成された第2の素子データに基づいて超音波画像を生成する画像生成部と、を有し、素子データ処理部は、画像生成部で生成されるべき超音波画像上の1ラインの第2の素子データを作成することが好ましい。
また、音速決定部は、素子データ処理部と、素子データ処理部における遅延時間算出に用いられる検査対象物の音速値を変更する音速変更部と、音速変更部で変更された音速値を用いて第2の素子データを作成する素子データ処理部と、素子データ処理部により作成された第2の素子データを、音速変更部で変更された音速値を用いて整相加算して音線信号を生成する整相加算部と、整相加算部で生成された音線信号に基づいて超音波画像を生成する画像生成部と、画像生成部で生成された超音波画像の画質指標を算出する画質指標算出部と、画質指標算出部で算出された超音波画像の画質指標を判定する判定部と、判定部による画質指標の判定結果に応じて最適な音速値を設定する音速設定部と、を有することが好ましい。
また、音速決定部は、音速変更部が、音速値を初期音速値から最終音速値まで所定音速間隔で変更する度毎に、素子データ処理部が変更された当該音速値を用いて第2の素子データを作成し、整相加算部が同じ当該音速値を用いて第2の素子データを整相加算して音線信号を生成し、画像生成部が音線信号に基づいて超音波画像を生成し、画質指標算出部が超音波画像の画質指標を算出することを繰り返して、音速変更部によって変更された全ての複数の音速値にそれぞれ対応する複数の画質指標を算出しておき、判定部が、画質指標算出部で全ての複数の音速値毎に算出された複数の画質指標の中から最適な画質指標を判定し、音速設定部が、判定部で判定された最適な画質指標に対応する音速値を最適な音速値として設定することが好ましい。
さらに、音速決定部よりも前段に、第1の素子データを用いて最適な音速値を算出する第2の音速決定部を有し、音速決定部は、第2の音速決定部で算出した音速値を初期値として用いることが好ましい。
また、上記目的を達成するために、本発明の超音波検査方法は、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子を用いて検査対象物を検査する超音波検査方法であって、探触子によって、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の素子がアナログ素子信号を出力し、このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとし、検査対象物内の音速値を用いて、2つ以上の第1の素子データから、第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として決定することを特徴とする。
ここで、最適な音速を決定するステップは、音速値を初期音速値から最終音速値まで所定音速間隔で変更する毎に、変更された当該音速値を用いて第2の素子データを作成し、同じ当該音速値を用いて第2の素子データを整相加算して音線信号を生成し、生成された音線信号に基づいて超音波画像を生成し、画質指標算出部が超音波画像の画質指標を算出することを繰り返して、変更された全ての複数の音速値毎にそれぞれ対応する複数の画質指標を算出しておき、全ての複数の音速値毎に算出された複数の画質指標の中から最適な画質指標を判定し、判定された最適な画質指標に対応する音速値を最適な音速値として設定することが好ましい。
また、上記目的を達成するために、本発明のプログラムは、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して検査対象物を検査するために、複数の素子が配列された探触子によって、複数の素子を用い、所定の送信焦点を形成するように、超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の素子がアナログ素子信号を出力し、このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データを複数、取得するステップと、検査対象物内の音速値を用いて、得られた2つ以上の第1の素子データから第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として決定するステップとを手順としてコンピュータに実行させることを特徴とする。
また、上記目的を達成するために、本発明の記録媒体は、上記記載のプログラムを記録したコンピュータ読み取り可能な記録媒体である。
本発明によれば、超音波プローブの複数の素子からの超音波ビームを送信する際に、異なる2以上の素子からの送信で取得した素子データを重ねるマルチライン処理を行うことで、送信超音波ビームの広がりによる影響を低減し、SN比を上げ、解像度を上げることができ、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高い解像度で、最適な空間分解能を持つシャープな超音波画像を得ることができると共に、マルチライン処理の計算を含めた上で音速補正を行うことにより、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができる。
本発明に係る超音波検査装置の構成の一例を概念的に示すブロック図である。 図1に示す本発明の超音波検査装置の音速補正部を含む要部を処理フローに沿って示すブロック図であり、音速補正部の一例の詳細を示す。 図1に示す超音波検査装置の素子データ処理部の構成の一例を概念的に示すブロック図である。 (a)及び(c)は、それぞれ理想的な超音波ビームを被検体の反射点の真上の素子及び真上でない素子から送信する場合の説明図であり、(b)及び(d)は、それぞれ得られる素子データを示す説明図である。 (a)及び(c)は、それぞれ実際の超音波ビームを被検体の反射点の真上の素子及び真上でない素子から送信する場合の説明図であり、(b)及び(d)は、それぞれ得られる素子データを示す説明図である。 (a)及び(b)は、それぞれ真の反射超音波エコーとゴーストの反射信号の場合の超音波ビームの送信経路及び受信経路の距離を説明する説明図であり、(c)及び(d)は、それぞれ複数の素子で得られる素子データ及びそれらの遅延時間を示す説明図である。 (a)、(b)及び(c)並びに(d)、(e)及び(f)は、それぞれ真の信号の場合並びにゴーストの場合の複数の素子で得られる素子データ、それらの遅延時間及び素子データの重ね合わせ状態を示す説明図であり、(g)及び(h)は、それぞれ複数の素子に対応する素子データの重ね合わせ状態及びその結果を示す説明図である。 従来技術の超音波診断装置における超音波の送受信を説明する説明図である。 本発明に係る超音波検査装置の実施例2の構成の一例を概念的に示すブロック図である。 図9に示す本発明の超音波検査装置の音速補正部を含む要部を処理フローに沿って示すブロック図であり、音速補正部の一例の詳細を示す。
本発明に係る超音波検査装置、超音波検査方法、プログラム及び記録媒体を添付の図面に示す好適実施形態に基づいて以下に詳細に説明する。
図1は、本発明の超音波検査装置の構成の一実施例を概念的に示すブロック図である。
同図に示すように、超音波検査装置10は、超音波プローブ12と、超音波プローブ12に接続される送信部14及び受信部16と、A/D変換部18と、素子データ保持部20と、素子データ処理部22を備える音速補正部24と、画像生成部26と、表示制御部27と、表示部28と、制御部30と、操作部32と、格納部34とを有する。
超音波プローブ12は、通常の超音波検査装置に用いられる探触子36を有する。
探触子36は、1次元又は2次元アレイ状に配列された複数の素子、即ち超音波トランスデューサを有している。これらの超音波トランスデューサは、検査対象物(以下、被検体という)の超音波画像の撮像の際に、それぞれ送信部14から供給される駆動信号に従って超音波ビームを被検体に送信すると共に、被検体からの超音波エコーを受信して受信信号を出力する。本実施形態では、探触子36の複数の超音波トランスデューサの内の一組を成す所定数の超音波トランスデューサの各々は、1つの超音波ビームの各成分を発生し、一組の所定数の超音波トランスデューサは、被検体に送信する1つの超音波ビームを発生する。
各超音波トランスデューサは、例えば、PZT(チタン酸ジルコン酸鉛)に代表される圧電セラミックや、PVDF(ポリフッ化ビニリデン)に代表される高分子圧電素子、PMN−PT(マグネシウムニオブ酸・チタン酸鉛固溶体)に代表される圧電単結晶等からなる圧電体の両端に電極を形成した素子、即ち振動子によって構成される。即ち、探触子36は、複数の超音波素子として複数の振動子が1次元又は2次元アレイ状に配列された振動子アレイということができる。
このような振動子の電極に、パルス状又は連続波状の電圧を印加すると、圧電体が伸縮し、それぞれの振動子からパルス状又は連続波状の超音波が発生して、それらの超音波の合成により超音波ビームが形成される。また、それぞれの振動子は、伝搬する超音波を受信することにより伸縮して電気信号を発生し、それらの電気信号は、超音波の受信信号として出力される。
送信部14は、例えば、複数のパルサを含んでおり、制御部30からの制御信号に応じて選択された送信遅延パターンに基づいて設定される音速または音速の分布に従い、探触子36の一組の所定数の超音波トランスデューサ(以下、超音波素子という)から送信される超音波ビーム成分が1つの超音波ビームを形成するようにそれぞれの駆動信号の遅延量を調節して組を成す複数の超音波素子に供給する。
受信部16は、制御部30からの制御信号に応じて、探触子36の各超音波素子によって超音波ビームと被検体との間の相互作用によって発生された超音波エコーを被検体から受信して、受信信号、即ち超音波素子毎のアナログ素子信号を増幅して出力し、増幅されたアナログ素子信号をA/D変換部18に供給する。
A/D変換部18は、受信部16に接続され、受信部16から供給されたアナログ素子信号を、デジタル素子データに変換する。A/D変換部18は、A/D変換されたデジタル素子データを素子データ保持部20に供給する。
素子データ保持部20は、A/D変換部18から出力されるデジタル素子データを順次格納する。また、素子データ保持部20は、制御部30から入力されるフレームレートに関する情報(例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータ)を上記のデジタル素子データ(以下、単に素子データという)に関連付けて格納する。
素子データ保持部20は、制御部30からの制御信号に応じて、2次元以上の位置座標上の対象領域において重なり合う少なくとも2つ以上の対象領域を検査する際に、この2つ以上の対象領域毎に受信部16で受信した超音波エコーから2つ以上の対象領域毎に生成された2以上の素子データを記憶保持するものであり、素子データ保持部20に記憶保持される素子データは、素子データ毎に受信した各素子における受信時間毎の受信データを含む、2つ以上の素子データである。即ち、素子データ保持部20は、2つ以上の素子データを含む一部の素子データ、あるいは全ての第1の素子データを保持する。
また、素子データ処理部22は、音速補正部24の内に設けられ、制御部30による制御に基づいて、被検体の検査対象領域内の所定の計算領域の計算座標(所定領域内のサンプリング点や観測点の集まり)について、その計算座標の音速(以下、環境音速ともいう)を用いてマルチライン処理を行って、即ち、素子データ保持部20によって記憶保持された2つ以上の対象領域毎に生成された2以上の素子データ(以下、未処理素子データという)を受信時間上で時間を合わせて重ね合わせて重ね合わせ処理後素子データ(以下、処理済素子データという)を生成するものである。なお、本発明においては、計算座標の音速又は環境音速とは、探触子36の超音波素子から送信された超音波ビームが被検体の検査対象領域内の所定の計算座標までの送信経路及び計算座標で反射した超音波エコーが探触子36の超音波素子に戻るまでの受信経路における平均音速、即ち、探触子36の超音波素子と計算座標との間の平均音速を言う。
音速補正部24は、素子データ保持部20と画像生成部26の整相加算部38との間に設けられ、内部にマルチライン処理を行う素子データ処理部22を備えるもので、制御部30による制御に基づいて、被検体の検査対象領域内の所定の計算領域の計算座標毎に、素子データ処理部22によるマルチライン処理を行った処理済データに基づいて、整相加算部38で用いられる送受信の最適な音速を求め、被検体の検査対象領域内の環境音速、したがって、被検体の検査対象領域の超音波画像に最も適した環境音速を求めると共に、素子データ処理部22によって最適な環境音速を用いるマルチライン処理が行われた処理済データを得るものである。
なお、素子データ処理部22及び音速補正部24の詳細については、後述する。
画像生成部26は、制御部30による制御下で、素子データ処理部22で生成され、音速補正部24から供給された処理済素子データから音線信号(受信データ)を生成し、この音線信号から超音波画像を生成するものである。
画像生成部26は、整相加算部38、検波処理部40、DSC42、画像作成部44、および、画像メモリ46を有する。
整相加算部38は、制御部30において設定された受信方向に応じて、音速補正部24で求められた最適な環境音速またはその分布に従い、素子データ処理部22で生成された処理済素子データにそれぞれの遅延を与えて加算することにより、受信フォーカス処理を行う。この受信フォーカス処理により、超音波エコーの焦点が絞り込まれた受信データ(音線信号)が生成される。なお、音速補正部24において最適な環境音速またはその分布が求まっていない場合には、予め記憶されている複数の受信遅延パターンの中から1つの受信遅延パターンを選択し、選択された受信遅延パターンに基づいて設定される環境音速またはその分布を用いても良い。
整相加算部38は、受信データを検波処理部40に供給する。
検波処理部40は、整相加算部38で生成された受信データに対し、超音波の反射位置の深度に応じて距離による減衰の補正を施した後、包絡線検波処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像データを生成する。
DSC(digital scan converter)42は、検波処理部40で生成されたBモード画像データを通常のテレビジョン信号の走査方式に従う画像データに変換(ラスター変換)する。
画像作成部44は、DSC42から入力されるBモード画像データに階調処理等の各種の必要な画像処理を施して検査や表示に供するためのBモード画像データを作成した後、作成された検査用又は表示用Bモード画像データを表示のために表示制御部27に出力する、或いは画像メモリ46に格納する。
画像メモリ46は、画像作成部44で作成された検査用Bモード画像データを一旦格納する。画像メモリ46に格納された検査用Bモード画像データは、必要に応じて、表示部28で表示するために表示制御部27に読み出される。
表示制御部27は、画像作成部44によって画像処理が施された検査用Bモード画像信号に基づいて、表示部28に超音波画像を表示させる。
表示部28は、例えば、LCD等のディスプレイ装置を含んでおり、表示制御部27の制御の下で、超音波画像を表示する。
制御部30は、操作者により操作部32から入力された指令に基づいて超音波検査装置10の各部の制御を行う。
ここで、制御部30は、種々の情報、特に、音速補正部24で用いられる最適環境音速の算出に必要な情報、素子データ処理部22及び画像生成部26の整相加算部38で最適環境音速と共に用いられる遅延時間算出に必要な情報及び素子データ処理部22における素子データ処理に必要な情報が、操作者によって操作部32を介して入力された際に、かつ/又は格納部34から読み出された際に、操作部32から入力された、かつ/又は格納部34から読み出された上述の種々の情報を、必要に応じて、送信部14、受信部16、素子データ保持部20、素子データ処理部22、音速補正部24、画像生成部26及び表示制御部27等の各部に供給する。
操作部32は、操作者が入力操作を行うためのものであり、キーボード、マウス、トラックボール、タッチパネル等から形成することができる。
また、操作部32は、操作者が必要に応じて各種の情報、特に上述の遅延時間算出に用いられるプローブ12の探触子36の複数の超音波素子、被検体の検査対象領域、その計算領域(計算座標X1〜Xend)、音速補正(最適環境音速決定)のための音速探索範囲(初期音速値Vst、最終音速値Vend、音速値変更幅(音速間隔)ΔV)、超音波ビームの焦点位置、探触子36の送信開口及び受信開口等に関する情報、並びに重ね合わせ素子データ数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報等を入力操作するための入力装置を備えている。
格納部34は、操作部32から入力された各種の情報、特に、上述のプローブ12、計算領域、音速探索範囲、焦点位置、送信開口及び受信開口等に関する情報、重ね合わせ素子データ数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報等、及び送信部14、受信部16、素子データ保持部20、素子データ処理部22、音速補正部24、画像生成部26及び表示制御部27等の制御部30で制御される各部の処理や動作に必要な情報、並びに、各部の処理や動作を実行させるための動作プログラムや処理プログラム等を格納するもので、ハードディスク、フレキシブルディスク、MO、MT、RAM、CD−ROM、DVD−ROM等の記録媒体を用いることができる。
なお、素子データ処理部22、音速補正部24、画像生成部26の整相加算部38、検波処理部40、DSC42並びに画像作成部44、及び表示制御部27は、CPUと、CPUに各種の処理を行わせるための動作プログラムから構成されるが、それらをデジタル回路で構成してもよい。
次に、図1に示す素子データ処理部22を備える音速補正部24について説明する。
図2は、図1に示す本発明の超音波検査装置の音速補正部の一例の詳細含む超音波検査装置の要部を処理フローに沿って示すブロック図である。
音速補正部24は、設定音速Vを検査対象領域内の音速として推定される所定音速範囲全域を所定音速値変更幅(音速間隔)で変化させて、例えば、人体等の生体を対象とする場合、具体的には、1400m/sから1650m/s又は1700m/sまで、1m/s〜50m/sの所定刻み(ステップ)、例えば、10m/s刻みで変化させて、検査対象領域内の画像の輝度値、鮮鋭度等の画質指標(画質評価指標)や音速の収束度合い等に基づいて網羅的に最適音速を求めるものである。
このような音速補正部24による音速補正方法は、例えば、特開2011−92686号公報に開示の音速補正方法を用いることができる。なお、この公報に開示の音速補正方法では、着目領域毎に送信フォーカスを結ぶ必要がある。しかしながら、本発明では、音速補正部24内の素子データ処理部22の重ね合わせ処理部の重ね合わせ処理によって擬似的に如何なる深さにおいても送信フォーカスを結ぶことができている処理済素子データが得られているので、計算領域や着目領域毎に送信フォーカスを設定する必要はない。ただし、本発明においても、この公報に開示の音速補正方法の通りに、計算領域や着目領域毎に送信フォーカスを設定するようにしても良いのは勿論である。
図2に示す実施形態においては、音速補正部24は、計算領域設定部70と、初期計算座標設定部72と、最適音速決定部23と、最終計算座標判定部88と、計算座標更新部90とを有し、最適音速決定部23は、音速探索範囲設定部74と、初期音速値設定部76と、整相加算部38aと、検波処理部40aと、画像作成部44aと、画質指標算出部78と、最終音速値判定部80と、音速値更新部82と、最適音速判定部84と、音速設定部86と、を有する。
音速補正部24では、図2に示すように、検査対象領域内の計算座標X1〜Xendにおける設定音速Vを変化させ、それぞれの設定速度Vにおいて設定速度Vに基づいて、素子データ保持部20から供給された未処理素子データを素子データ処理部22でマルチライン処理した処理済素子データを用いて整相加算部38aにおいて受信フォーカス処理をして整相加算を行い音線信号を生成し、この音線信号から画像作成部44aにおいて超音波画像を形成し、それぞれの設定音速Vでの検査対象領域内の超音波画像の輝度値、鮮鋭度等の画質指標を算出する。
音速補正部24では、各設定音速Vでの超音波画像の輝度値、鮮鋭度等の画質指標を最適音速判定部84において比較して、例えば、得られた輝度値、鮮鋭度等の画質指標の値が最も高い設定音速Vを判定し、判定された設定音速Vを音速設定部86において最適な音速値(最適音速値)として設定する。
なお、音速補正部24においては、音速値の変更及び最適音速の設定は、図2に示すように網羅的に行っても良いし、探索的に行っても良い。
また、音速補正部24における音速値の変更及び最適音速の設定の際に、音速値を設定する領域や計算座標は、特に制限的ではなく、被写体全体や検査対象領域全体を同一の音速値として設定してもよいが、被写体や検査対象領域を微小領域ごとに区切って、微小領域ごとに音速値を設定してもよく、画素単位で1点ずつ音速値を設定しても良い。
まず、図2に示す音速補正部24の素子データ処理部22について、図3に基づいて詳細に説明する。
同図に示すように、素子データ処理部22は、マルチライン処理を行うためのもので、遅延時間算出部48と、重ね合わせ処理部50とを有する。
なお、ここでは、素子データ処理部22でマルチライン処理を行う被検体の検査対象領域内の計算座標は、予め、初期計算座標設定部72によって初期計算座標X1に、又は計算座標更新部90によって更新された計算座標Xi(i=2〜end:以下、計算座標Xi(i=1〜end)で代表する)に設定されているものとし、その音速値は、初期音速値設定部76によって初期音速値Vstに、又は音速値更新部82によって更新された音速値V(以下、音速値Xで代表する)に予め設定されているものとする。また、素子データ処理部22におけるマルチライン処理に必要である、計算座標Xiに関わる未処理素子データは、素子データ保持部20から読み出されて、供給されているものとする。
遅延時間算出部48は、操作部32から入力された、もしくは、操作部32から入力されて格納部34に格納されているプローブ12の探触子36の複数の超音波素子、超音波ビームの焦点位置、探触子36の送信開口及び受信開口等に関する情報を事前に取得しておき、被検体の検査対象領域内の計算座標Xiの音速は、初期音速値設定部76又は音速値更新部82によって予め設定された音速値V(Vst)を用いて、超音波ビームを形成し送信する、送信開口の超音波素子(送信素子)と、被検体からの、超音波ビームによる超音波エコーを受信する、受信開口の超音波素子との幾何学的配置に基づいて、受信開口の各超音波素子で受信される素子データの遅延時間を算出する。
重ね合わせ処理部50は、操作部32から入力された、もしくは、操作部32から入力されて格納部34に格納されている、重ね合わせる素子データの数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報に基づいて、素子データ保持部20によって記憶保持された2以上の対象領域毎に生成された2以上の未処理素子データを読み出し、遅延時間算出部48でそれぞれ算出された遅延時間に基づいて、2以上の未処理素子データを受信時間上で、即ち時間を合わせて、かつ、受信された探触子の素子の絶対的な位置を合わせて、重ね合わせて処理済素子データを生成する。
次に、素子データ処理部22で行う素子データのマルチライン処理について詳細に説明する。
まず、超音波プローブ12の探触子36の送信用超音波素子(以下、単に送信素子という)から超音波ビーム(以下、送信ビームという)を被検体に送信し、被検体との間の相互作用によって発生された超音波エコーを探触子36の受信用超音波素子(以下、単に受信素子という)で受信して素子データを得る場合において、送信素子からの送信ビームと受信素子で得られる素子データとの関係について説明する。
図4(a)及び(c)に示すように、それぞれ7つの超音波素子(以下、単に、素子ともいう)52a〜52g及び52b〜52hを受信素子として超音波エコーを受信して素子データを取得する時、反射点54を含む検査対象領域に送信する送信ビーム56が理想的に素子間隔以下に絞れている理想的な場合には、図4(a)のように、検査対象領域内の反射点54の真上にある、素子52a〜52gの中心の素子52dを送信素子として送信ビーム56を送信し、受信素子52a〜52gで超音波エコーを受信して素子データを取得する場合、送信ビーム56の焦点58は、素子54dと反射点54とを結ぶ一直線上にあり、送信ビーム56は、反射点54まで送信されるので、反射点54から反射される超音波エコーが生成される。反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子52a〜52gに受信され、受信素子52a〜52gによって、図4(b)に示すような素子データ62が得られることになる。
これに対し、図4(c)に示すように、送信素子の中心を、反射点54に対して1素子分、素子の方向(図中右方向)にシフトさせ、反射点54の真上にある素子52dに隣接する素子52eを送信素子として送信ビーム56を送信し、受信素子52b〜52hで超音波エコーを受信する場合には、送信ビーム56の送信方向、即ち、送信素子52eと焦点58とを結ぶ直線上に反射点54が存在しないため、送信ビーム56は、反射点54に送信されない。このため、反射点54から反射される超音波エコーは生成されず、受信素子52b〜52hでは、超音波エコーを受信しないので、図4(d)に示すように、素子データが得られないことになる。
しかしながら、図5(a)及び(c)に示すように、実際の送信ビーム64は、素子間隔より幅が広い。
ここで、図5(a)のように、反射点54の真上にある素子52dを送信素子として送信ビーム64を送信した場合には、図4(a)の場合と同様に、送信ビーム64が幅広であっても、その焦点58は、素子54dと反射点54とを結ぶ一直線上にあり、送信ビーム64は、反射点54で反射され、超音波エコーが生成される。その結果、図4(a)の場合と同様に、反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子52a〜52gに受信され、受信素子52a〜52gによって、図5(b)に示すような真の素子データ66が得られることになる。
一方、図5(c)に示すように、図4(c)の場合と同様に、送信素子の中心を、反射点54に対して1素子分、素子の方向(図中右方向)にシフトさせ、反射点54の真上にある素子52dに隣接する素子52eを送信素子として送信ビーム64を送信し、受信素子52b〜52hで超音波エコーを受信する場合、送信ビーム64は幅広であるため、その送信方向、即ち、送信素子52eと焦点58とを結ぶ直線上に反射点54が存在していなくても、送信ビーム64は、反射点54に送信されることになる。このため、反射点54から本来存在しない超音波エコー、所謂ゴーストの反射信号が発生し、反射点54からのゴーストの反射信号は、所定角度に拡がる受信経路60を通って受信素子52b〜52hに受信され、受信素子52b〜52hによって、図5(d)に示すようなゴーストの素子データ68が得られることになる。
このようなゴーストの素子データ68は、素子データから生成される超音波画像の精度を低下させる原因となる。
ここで、図5(c)に示す送信ビーム64が送信素子52eから焦点58を経由して反射点54に至る送信経路とゴーストの反射信号が反射点54から各受信素子52b〜52hに至る受信経路との和(伝播距離)は、それぞれ図5(a)に示す送信ビーム64が送信素子52dから焦点58を経由して反射点54に至る送信経路と真の反射超音波エコーが反射点54から各受信素子52a〜52gに至る受信経路との和(伝播距離)より長くなるため、図5(d)に示すようなゴーストの素子データ68は、図5(b)に示すような真の素子データ66に対して遅延することになる。
本発明の素子データ処理部22の遅延時間算出部48においては、真の素子データに対するゴーストの素子データの時間差、即ち遅延時間は、送信素子、超音波ビームの焦点、被検体の反射点、及び受信素子の幾何学的配置から算出される。したがって、遅延時間の計算には、超音波プローブ12の形状(素子間隔、リニア、コンベックスなど)、焦点位置、送信開口、受信開口などの情報と、被検体の検査対象領域の音速値とが必要であり、遅延時間算出部48では、操作部32によって入力された、若しくは格納部34に格納されたこれらの情報と、音速値更新部82等で予め設定された音速値Vを取得して遅延時間の計算を行う。遅延時間は、音速値Vを用いて、例えば送信素子、超音波ビームの焦点、被検体の反射点、及び受信素子の幾何学的配置から算出される、送信素子から焦点を経て反射点に至る送信ビームの送信経路及び反射点から受信素子に至る真の反射超音波エコー又はゴーストの反射信号の受信経路の合計長さ(伝播距離)とその音速によって算出される伝播時間の差から算出することができる。
本発明では、例えば、図6(a)及び図6(b)に示すようにして、真の反射超音波エコーとゴーストの反射信号の場合の送信ビームの送信経路及び受信経路の長さを求めることができる。
真の反射超音波エコーの場合、図6(a)に示すように、送信素子52dと受信素子52d(受信素子52a〜52gの中心)とが一致し、その真下に、焦点58及び反射点54が配置されており、反射点54の真上の素子52dの位置をxy2次元座標上の座標(x0、0)とし、素子間隔をLe、焦点58の位置を座標(x0、df)、反射点54の位置を座標(x0、z)とする時、送信素子52dの位置も反射点54の真上の素子52dと同じく座標(x0、0)となり、送信素子52dから焦点58を経て反射点54に至る送信ビームの送信経路61の長さ(送信経路距離)Ltaは、及び、反射点54から受信素子52dに至る真の反射超音波エコーの受信経路60の長さ(受信経路距離)Lraは、Lta=Lra=zによって算出することができる。
したがって、真の反射超音波エコーの場合の超音波の伝播距離Luaは、Lua=Lta+Lra=2zとなる。
ゴーストの反射信号の場合、図6(b)に示すように、図6(a)の場合と比べて、送信素子52eの位置が反射点54に対して1素子分横(x方向:図中右方向)にずれて、焦点58は送信素子52eの真下に配置されるが、反射点54は受信素子52dの真下に配置されており、反射点54の真上の受信素子52dの位置を、図6(a)の場合と同じくxy2次元座標上の座標(x0、0)とし、素子間隔をLe、反射点54の位置を座標(x0、z)とすると、送信素子52eの位置は座標(x0+Le、0)、焦点58の位置は座標(x0+Le、df)となるので、送信素子52eから焦点58を経て反射点54に至る送信ビームの送信経路61の長さ(送信経路距離)Ltbは、Ltb=df+√{(z−df)+Le}によって算出することができ、反射点54から受信素子52dに至るゴーストの反射信号の受信経路60の長さ(受信経路距離)Lrbは、Lrb=zによって算出することができる。
したがって、ゴーストの反射信号の場合の超音波の伝播距離Lubは、Lub=Ltb+Lrb=df+√{(z−df)+Le}+zとなる。
こうして、図6(a)に示す幾何学配置で求めた送信経路61の距離Ltaと受信経路60の距離Lraを合計した超音波の伝播距離Luaを音速で割った値が真の反射超音波エコーの伝播時間となり、図6(b)に示す幾何学配置で求めた送信経路61の距離Ltbと受信経路60の距離Lrbを合計した超音波の伝播距離Lubを音速で割った値がゴーストの反射信号の伝播時間となる。遅延時間の算出は、反射点54と送信素子52(52d)のx座標が一致している時の真の超音波エコーの伝播時間から、反射点54と送信素子52(52e)のx座標を1素子間隔ずつずらしたときのゴーストの反射信号の伝播時間の差から遅延時間を求める。
なお、図6(a)及び図6(b)の幾何学モデルでは、送信経路61が焦点58を経由したモデルになっているが、本発明はこれに限定されず、例えば、焦点58を経由せずに直接反射点54に至る経路であっても良い。
また、図6(a)及び図6(b)の幾何学モデルはリニアプローブの場合であるが、これに限らず他のプローブにおいても、プローブの形状から同様の幾何学計算を行うことができる。例えば、コンベックスプローブの場合、プローブの半径と素子間隔の角度から幾何学モデルを設定して同じように計算することができる。
また、ステア送信の場合には、送信角度などの情報を考慮した幾何学モデル(図示せず)を用い、送信素子と反射点との位置関係から真の超音波エコーの素子データ及びその周辺のゴーストの素子データの遅延時間を算出することができる。
さらに、幾何学モデルよって遅延時間を算出する方法に限らず、あらかじめ装置の計測条件に合わせて高輝度反射点を計測した計測結果から、計測条件毎に遅延時間を求めておき、その遅延時間を装置内に記憶しておくことで、同じ計測条件の遅延時間を読み出すようにしておいてもよい。
図6(c)に、中央に真の信号の素子データ66及び真の信号の周辺のゴーストの素子データ68を示し、図6(d)に、上述の幾何学計算から得られた素子データ66に対するゴーストの素子データ68の遅延時間の一例を示す。真の信号の素子データ66を中心にゴーストの信号の素子データ68は、対称的に時間が遅れることが示されている。
なお、こうして、素子データ処理部22の遅延時間算出部48において算出された遅延時間を整相加算部38における遅延補正に用いることもできる。
次に、本発明の素子データ処理部22の重ね合わせ処理部50においては、こうして遅延時間算出部48において算出された遅延時間を用いて真の超音波エコーの素子データ及びその周辺のゴーストの素子データの重ね合わせ処理を行う。
重ね合わせ処理部50における重ね合わせ処理では、重ね合わせる時の重ね合わせ素子データ数と重ね合わせ処理方法の情報が必要になるが、これらは、予め、操作部32によって入力しておいても良いし、格納部34に格納しておいても良い。
なお、重ね合わせ処理部50において重ね合わせる時の重ね合わせ素子データの範囲は、計算領域設定部70で設定される音速補正計算の対象とする領域内にとどまる範囲にある素子データであることが望ましいが、計算領域設定部70で設定される領域によって制限されるものではない。例えば、重ね合わせる素子データが隣接する領域同士でオーバーラップするような範囲で重ね合わせてもよい。
図7(a)〜(h)に、重ね合わせ処理部50で行われる、素子データ数が5つ、重ね合わせ素子データ数が3つの場合の重ね合わせ処理の一具体例を示す。
図7(a)は、5つの素子データを横に並べて表示しており、素子データ毎に、超音波ビームを送信し、反射信号を受信した様子を表している。各素子データの横軸は、受信素子を表しており、それぞれの素子データにおいて超音波ビームの送信時における中心の素子を中心にして表示している。縦軸は、受信時間を表す。
5つの素子データのうち、真中の素子データでは、素子データの中心の素子(受信素子の中心の素子)、即ち、送信時における中心の素子(送信素子)の真下に反射点が存在しており、反射点からの反射信号が受信されている。つまり、この反射信号は真の信号であり、真中の素子データは、真の信号を表す。
真中の素子データ以外の両側2つの素子データについては、送信時における中心の素子の真下には反射点は存在していないが、送信した超音波ビームの広がりによって、真中の素子データの送信素子の真下に存在する反射点に超音波ビームが当たることで生じた反射信号、即ちゴーストが写り込んでいる。ゴーストは、真の信号から離れるほど反射点までの超音波の伝播時間が長くなるため、真の信号よりも受信時間が遅くなる。また、反射点からの反射信号が初めに受信される受信素子の位置は、反射点の真上の素子であるが、素子データの横軸は超音波ビームの送信時における中心の素子を中心にしているため、素子データ毎にこの中心の素子を1素子ずつずらして送信していることから、各素子データにおいて素子の絶対位置は1素子ずつずれている。つまり、真中の素子データでは、反射点からの反射信号がはじめに受信される受信素子は真中の素子であるが、両隣の素子データにおいては、真中の素子データよりも1素子ずれており、右側の素子データでは左に1素子ずれ、左側の素子データでは右に1素子ずれている。更に、両端の素子データでは、真中の素子データよりも2素子ずれており、右端の素子データでは左に2素子ずれ、左端の素子データでは右に2素子ずれている。このように、ゴーストの信号は、真の信号に対して、受信時間が遅れるだけでなく、受信素子の方向に対してもずれを生じている。
図7(b)に、図7(a)に示す5つの素子データのうちの真中の素子データに対する受信時間の遅延時間の一例を示す。
重ね合わせ処理部50では、図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つの素子データの真中の素子データに対する受信時間の遅延時間とも同じである。
重ね合わせ処理部50では、図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素子データでは、ゴーストの素子データは互いに位相が合わない素子データを加算し、又は平均するので、互いに打ち消し合うことになるため、ゴーストの重ね合わせ処理済素子データは、その値が真の信号の素子データである高輝度値を持つ重ね合わせ処理済素子データに対して小さくなり、真の信号の素子データに対してゴーストの素子データの影響を低減させることができ、又は、その影響を無視できる程、小さくすることができる。
そのため、処理済素子データに整相加算や検波処理を行って、受信データを生成して、超音波画像を生成することにより、同様に、ゴーストの影響を無くし、すなわち音線上の全ての点で焦点を結んだのに等しい素子データで超音波画像生成できるので、高輝度で、鮮鋭性に優れた、高画質な超音波画像を生成することができる。
なお、本発明においては、この処理済素子データの生成を、マルチライン処理と言う。
本発明において、中心素子とは、送信の開口数(超音波の送信を行う素子数)が奇数の場合には、方位方向の中央の素子である。
他方、開口数が偶数の場合には、方位方向の中央の素子のいずれかを中心素子とし、あるいは、方位方向の真中に素子が有ると仮定して、中心素子とする。すなわち、開口数が偶数の場合には、開口の真中のライン上に焦点が有るものとして計算を行ってもよい。
なお、重ね合わせ処理部50における重ね合わせ処理方法としては、単に、加算するだけでなく、平均値や中央値をとってもよいし、係数を掛け合わせた上(重み付して)で加算してもよい。なお、平均値や中央値を取ることは、素子データレベルでの平均化フィルタやメディアンフィルタを掛けることに相当すると考えられるが、平均化フィルタやメディアンフィルタの代わりに、通常の画像処理で行われる逆フィルタなども適用してもよい。或いは、重ね合わせる各素子データ同士を比較し、類似している場合には最大値、類似していない場合には平均値、分布の偏りがある場合には中間値をとるなど、これに限らず、重ね合わせる各素子データの特徴量に基づいて重ね合わせ処理を変えてもよい。
また、注目素子の素子データに重ね合わせる素子データの数は、図示例の2つに限定はされず、1つでもよく、あるいは、3つ以上でもよい。すなわち、注目素子の素子データに重ね合わせる素子データの数は、要求される処理速度(フレームレート等)や画質などに応じて、適宜、設定すればよい。
また、重ね合わせる素子データ数は、超音波ビームのビーム幅の広がり程度に合わせた方が望ましい。従って、深さによってビーム幅が変わる場合には、重ね合わせ素子データ数も深さによって変更してもよい。また、ビーム幅は送信開口数に依存することから、送信開口数に応じて重ね合わせ素子データ数を変更してもよい。或いは、画像の輝度値などの特徴量に基づいて重ね合わせ素子データ数を変更してもよいし、重ね合わせ素子データ数を複数パターン変えて作成した画像から最適な重ね合わせ素子データ数を選択してもよい。
重ね合わせた結果、上述したように、真の信号の素子データでは信号の位相が合うが、ゴーストでは信号の位相が合わないため、加算などの重ね合わせ処理の結果、様々な位相の信号がお互いに打ち消し合い、信号が弱くなる。結果的に、真の信号は、有効な値を持つ、例えば高輝度の素子データとして残り、ゴーストの信号は、減弱した値を持つ、例えば低輝度の素子データとして得ることができる。
本発明に用いられる素子データ処理部は、基本的に以上のように構成される。
再び、図2に戻って、音速値の変更及び最適音速の設定を網羅的に行う音速補正部24の各部についての説明を続ける。
図2に示す音速補正部24において、計算領域設定部70は、操作者による操作部32(図1参照)からの入力等に応じて、検査対象領域内の音速補正計算の対象とする領域、例えば、着目領域ROIを設定し、設定された領域のxy座標上の計算対象とする全ての点を抽出し、計算座標X1〜Xendとして設定する。ここで、計算座標X1〜Xendのそれぞれで表される各領域は、環境音速が同一であると考えられる領域であれば、特に制限的ではなく、被写体全体や検査対象領域全体であってもよいが、被写体や検査対象領域を区切った微小領域であってもよく、画素単位で1点ずつであっても良い。
初期計算座標設定部72は、計算座標Xiを、計算領域設定部70で設定された計算座標X1〜Xendの中の初期計算座標X1に設定する。初期計算座標設定部72で設定された計算座標X1のデータは、最適音速決定部23に送られる。
なお、計算領域設定部70は、操作部32からの入力に応じて計算領域を設定する構成に限定はされず、画像生成部26(図1参照)の画像作成部44で生成された超音波画像(Bモード画像信号)を解析して、例えば着目領域ROIを設定し、計算領域とする構成としてもよい。
例えば、計算領域設定部70は、画像作成部44で生成された超音波画像において、隣接する画素同士の輝度値の差が所定の値以上となる位置(画素)を構造領域(組織や病変)の辺縁部(エッジ)と判定して、エッジに囲まれた領域、あるいは、エッジを多く含む領域を構造領域、即ち、対象部位Pと判定して抽出し、対象部位を含み、かつ、所定の形状、大きさとなる領域を、例えば着目領域ROIに設定し、計算領域に設定してもよい。
なお、超音波画像を解析して計算領域を設定する場合には、構造領域を1種類の領域とみなしても良く、あるいは、エッジの連続性や、輝度値のレベル等に応じて2種以上の領域とみなして、各領域を対象部位と判定して抽出し、対象部位毎に着目領域ROIを設定し、計算領域としても良い。
最適音速決定部23は、初期計算座標設定部72で設定された初期計算座標X1、又は計算座標更新部90で更新された計算座標Xiにおける最適な環境音速値(以下、最適音速値という)を算出して設定する部位である。
ここで、最適音速値とは、設定音速(設定環境音速)Vを種々変更し、それぞれの設定音速Vに基づいて、受信フォーカス処理をして超音波画像を形成した場合に、画像のコントラストおよび/またはシャープネスが最も高くなる設定音速Vの値であり、例えば、特開平8−317926号公報に記載のように、画像のコントラスト、スキャン方向の空間周波数、分散等に基づいて最適音速値の判定を行うことができる。
最適音速決定部23は、素子データ保持部20から読み出された、計算座標Xiに関わる未処理素子データを入力値として用いて最適な環境音速値(以下、最適音速値ともいう)を出力値として網羅的に求めて設定するためのものである。
最適音速決定部23としては、上述したように、所定の音速探索範囲を網羅的に探索して最適音速値を設定するものを挙げることができ、音速探索範囲設定部74と、初期音速値設定部76と、素子データ処理部22と、整相加算部38aと、検波処理部40aと、画像作成部44aと、画質指標算出部78と、最終音速値判定部80と、最適音速判定部84と、音速値更新部82と、音速設定部86と、を有する。
音速探索範囲設定部74は、計算座標Xiにおける環境音速の探索範囲(初期音速値及び最終音速値)及び音速値変更幅(音速間隔)を設定するためのものである。音速探索範囲設定部74は、例えば、音速の探索開始の初期音速値Vst、音速の探索終了の最終音速値Vend、設定音速Vを変更するステップ幅(音速変更幅:音速間隔)ΔVを設定する。一例として、生体等の場合、音速探索範囲設定部74は、初期音速値Vstを1400m/sに、最終音速値Vendを1700m/sに、音速変更幅ΔVを10m/sに設定することができる。
なお、音速探索範囲設定部74が設定する初期音速値Vst、最終音速値Vend及び音速変更幅ΔVは、上述の値に限定されず、被検体やその検査対象領域やその計算座標Xi等に応じて決定されるものであればよい。また、これらの初期音速値Vst、最終音速値Vend及び音速変更幅ΔVは、音速探索範囲設定部74に、被検体やその検査対象領域やその計算座標Xi等に応じて、操作者によって操作部32を介して入力されるものであっても良いし、予め、格納部34に格納されているものの中から読み出されるものであっても良い。
初期音速値設定部76は、設定音速Vを音速探索範囲設定部74が設定した初期音速値Vstに設定するものである。上述の生体の場合、初期音速値設定部76は、設定音速Vを初期音速値Vstとして、1400m/sに設定する。初期音速値設定部76で設定された初期音速値Vstのデータは、後段の素子データ処理部22に送られる。
整相加算部38aは、整相加算部38と全く同様の構成及び機能を有し、設定音速Vを用いて素子データ処理部22で求められた計算座標Xiに関わる処理済素子データに対して、音速補正部24で設定された最適音速に基づいて遅延を与えて加算する受信フォーカス処理を行って音線信号を生成する。
検波処理部40aは、検波処理部40aと全く同様の構成及び機能を有し、整相加算部38aで生成された音線信号に対して、包絡線検波処理等の処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像データを生成する。
画像作成部44aは、画像作成部44と全く同様の構成及び機能を有し、検波処理部40aで生成されたBモード画像データに画像処理を施して超音波画像を作成する。
画質指標算出部78は、画像作成部44aで作成された超音波画像から設定音速Vにおける画質指標、即ち画像の輝度値や鮮鋭度等の画質指標、例えば、各設定音速Vでの計算座標Xi(例えば、着目領域ROI)の画像の鮮鋭度(シャープネス)等を算出する。算出された画質指標は、設定音速Vに対応付けられて画質指標算出部78内に保持される、例えば、画質指標算出部78内のメモリ等に保持される。
最終音速値判定部80は、設定音速Vを最終音速値Vendと比較して、両者が等しいか否かの判定を行い、設定音速Vが最終音速値Vendに等しくない(V≠Vend)場合には、音速値更新部82に移行させ、設定音速Vが最終音速値Vendに等しい(V=Vend)場合には、最適音速判定部84に移行させる。
音速値更新部82は、現在の設定音速Vに音速変更幅ΔVを加算して新しい設定音速V(V=V+ΔV)を設定する。音速値更新部82で更新設定された設定音速V(V=V+ΔV)のデータは、前段の素子データ処理部22に送られる。
最適音速判定部84は、画質指標算出部78に保持されている各音速値V(Vst〜Vend)毎に求められた画質指標を読み出し、各設定音速Vでの画質指標、例えば、超音波画像の鮮鋭度(シャープネス値)を比較し、最も画質がよい時、例えば、超音波画像の鮮鋭度(シャープネス値)が最も高い時の設定音速Vを判定するものである。
音速設定部86は、最適音速判定部84において生成された超音波画像の画質が最もよいと判定された設定音速Vを最適音速値として設定するものである。
なお、図2に示す最適音速決定部23では、設定音速Vの変更を繰り返して処理しているが、本発明はこれに限定されず、GPUなどのマルチスレッドを用いて、全探索音速値を並列計算し、一度に全ての音速値における画質指標を求めることも可能である。
最適音速決定部23は、音速設定部86において計算座標Xiに対応する最適音速値を設定した後、最終計算座標判定部88に移行する。
最適音速決定部23は、基本的に以上のように構成され、計算座標Xiに係る最適音速値を算出する。
最終計算座標判定部88は、計算座標Xiを最終計算座標Xendと比較して、両者が等しいか否かの判定を行い、計算座標Xiが最終計算座標Xendに等しくない(X≠Xend)場合には、計算座標更新部90に移行させ、計算座標Xiが最終計算座標Xendに等しい(X=Xend)場合には、画像生成部26の整相加算部38に移行させ、計算座標X1〜Xendの各々において算出された最適音速値及び最適音速値を用いて素子データ処理部22で処理された処理済素子データを計算座標X1〜Xendの各々の座標データと共に整相加算部38に供給する。
計算座標更新部90は、現在の計算座標Xiを新しい次の計算座標Xi+1に更新する。計算座標更新部90で更新された新しい次の計算座標Xi+1のデータは、前段の最適音速決定部23の音速探索範囲設定部74に送られる。
音速補正部24は、基本的に以上のように構成され、全ての計算座標Xi(i=1〜end)各々に係る最適音速値及び処理済素子データを算出する。
なお、上述したように、本実施形態の超音波検査装置10の音速補正部24による音速補正における最適音速を探索するループの中にマルチライン処理を行う素子データ処理部22を組み込み、マルチライン処理時の遅延時間補正を計算する際の音速値に、音速補正部24で探索する音速値を用いることで、反射点(サンプリング点、観測点)から超音波プローブ12の探触子36の受信素子までの音速だけでなく、各々の仮想的な点音源を形成する探触子36の送信素子から反射点までの音速も含めて最適な環境音速値を求めることができる。
本発明の超音波検査装置は、基本的に以上のように構成されるが、以下に、図1及び図2を参照して、本発明の超音波検査装置の動作、作用及び本発明の超音波検査方法について説明する。
操作者が、図1に示す超音波プローブ12を被検体の表面(検査対象領域)に当接させて、測定を開始すると、送信部14から供給される駆動信号に従って探触子36の複数の素子から超音波ビームが送信され、被検体からの超音波エコーを、探触子36の複数の素子が受信し、受信信号としてアナログ素子信号を出力する。
受信部16は、アナログ素子信号を増幅してA/D変換部18に供給し、A/D変換部18は、アナログ素子信号をデジタル素子データに変換して素子データ保持部20に供給して、未処理素子データとして記憶保持させる。
次に、図2に示す音速補正部24において、操作者による操作部32(図1参照)からの入力等に応じて、計算領域設定部70によって計算座標X1〜Xendが設定された後、初期計算座標設定部72によって計算座標Xiが、初期計算座標X1に設定され、設定された計算座標X1のデータが、最適音速決定部23に送られる。
最適音速決定部23では、まず、音速探索範囲設定部74が、被検体やその検査対象領域やその計算座標Xi等に応じて、操作者によって操作部32を介して入力された、又は予め、格納部34に格納されている音速探索範囲の情報の中から読み出された初期音速値Vst、最終音速値Vend及び音速変更幅ΔVを設定し、例えば、上述したように、生体では、初期音速値Vstを1400m/sに、音速値Vendを1700m/sに、音速変更幅ΔVを10m/sに設定する。
続いて、初期音速値設定部76が、設定音速Vを初期音速値Vst(例えば、1400m/s)に設定して、素子データ処理部22に移行する。
次に、素子データ処理部22は、遅延時間算出部48(図3参照)において、真の信号の未処理素子データに対する周辺のゴーストの信号の未処理素子データの遅延時間(例えば、図7(b)、図7(e)、どちらも同じものである)を、送信素子、焦点、反射点、及び受信素子の幾何学的配置、及び予め初期音速値設定部76で設定された設定音速V等から算出する(例えば、図6の幾何学モデルを用いて算出する)。
次に、素子データ処理部22は、素子データ保持部20から計算座標X1に関わる未処理素子データを素子データ保持部20から読み出し、処理を行う素子データを注目素子データとし、重ね合わせ処理部50(図3参照)において、遅延時間算出部48で算出された遅延時間を用いて、注目素子データとその周辺の未処理素子データとを位相合わせて重ね合わせてマルチライン処理を行い、処理済素子データを求める。これにより、真の信号を含む未処理素子データであれば、強調された処理済素子データが求まると共に、ゴーストの未処理素子データであれば減弱した処理済素子データが求まる。
素子データ処理部22は、こうして求めた処理済素子データを整相加算部38aに供給する。
次に、整相加算部38aは、初期音速値値Vstである設定音速Vに基づいて、素子データ処理部22から供給された計算座標X1に対応する処理済素子データを用いて受信フォーカス処理をして整相加算を行い音線信号を生成する。次いで、検波処理部40aは、生成された音線信号に対して検波処理を行ってBモード画像信号を生成する。そして、画像作成部44aは、生成されたBモード画像信号から超音波画像を作成する。
その後、画質指標算出部78は、作成された超音波画像から設定音速Vにおける超音波画像の輝度値や鮮鋭度等の画質指標、例えば、各設定音速VでのXi(X1)の画像の鮮鋭度等を算出する。算出された画質指標は、設定音速値Vに対応付けられて画質指標算出部78内に保持される。
次に、最終音速値判定部80は、設定音速Vを最終音速値Vendと比較して、両者が等しいか否かの判定を行う。最終音速値判定部80の判定において、設定音速Vが最終音速値Vendに等しい(V=Vend)場合には、最適音速判定部84に移行し、設定音速Vが最終音速値Vendに等しくない(V≠Vend)場合には、音速値更新部82に移行する。
音速値更新部82に移行した場合には、音速値更新部82が、現在の設定音速Vに音速変更幅ΔVを加算して新しい設定音速V(V=V+ΔV)を設定する。
この後、素子データ処理部22に戻り、新しい設定音速に基づいて素子データ処理部22によるマルチライン処理、整相加算部38aによる整相加算、検波処理部40aによる検波処理、画像作成部44aによる画像作成、画質指標算出部78による画質指標の算出及び保持、並びに最終音速値判定部80による判定が行われる。
最終音速値判定部80による判定結果がV≠Vend場合には、上述したように、音速値更新部82に移行し、音速値更新部82による新しい設定音速Vの設定、新しい設定音速Vに基づく素子データ処理部22によるマルチライン処理、整相加算部38aによる整相加算、検波処理部40aによる検波処理、画像作成部44aによる画像作成、画質指標算出部78による画質指標の算出及び保持、並びに最終音速値判定部80による判定が、V=Vendになるまで繰り返される。即ち、これらの処理を、音速範囲の終了の最終音速値Vendになるまで設定音速の値をΔVずつ変更しながら、例えば、設定音速Vを、1400m/sから1700m/sまで、10m/s刻みで変化させながら繰り返し行う。
一方、判定部80による判定結果がV=Vendである場合には、音速の探索範囲内の全ての設定音速Vに関して画質指標が求められたことになるので、上述したように、最適音速判定部84に移行し、最適音速判定部84は、画質指標算出部78に保持されている各音速値毎に求められた画質指標を読み出し、最適音速判定部84は、各設定音速Vでの画質指標、例えば、超音波画像の鮮鋭度(シャープネス値)を比較し、最も画質がよい時、例えば、超音波画像の鮮鋭度(シャープネス値)が最も高い時の設定音速Vを判定し、音速設定部86は、最適音速判定部84にて判定された設定音速Vを最適音速値として採用し、設定する。
こうして、最適音速決定部23は、計算座標Xi(X1)に対応する最適音速値を算出する。
次に、最適音速決定部23による計算座標Xi(X1)に対応する最適音速値の算出が終了すると、最終計算座標判定部88に移行する。
最終計算座標判定部88では、計算座標Xiを最終計算座標Xendと比較して、両者が等しいか否かの判定を行う。その判定の結果、計算座標Xiが最終計算座標Xendに等しくない(X≠Xend)場合には、計算座標更新部90に移行する。
この後、計算座標更新部90では、現在の計算座標Xiを新しい次の計算座標Xi+1に更新し、前段の最適音速決定部23に移行する。
この後、最適音速決定部23では、上述した設定音速Vの網羅的な探索を行い、計算座標Xi+1に対応する最適音速値を算出し、再び、最終計算座標判定部88に移行して、最終計算座標Xendの判定が行われる。
最終計算座標判定部88による判定結果がXi≠Xend場合には、上述したように、計算座標更新部90に移行し、計算座標更新部90による新しい計算座標Xiの設定、最適音速決定部23における設定音速Vの網羅的な探索による計算座標Xiに対応する最適音速値の算出、最終計算座標判定部88の判定が、Xi=Xendになるまで繰り返される。
一方、最終計算座標判定部88における判定の結果、計算座標Xiが最終計算座標Xendに等しい(X=Xend)場合には、音速補正部24による最適音速値の設定が終了し、全ての計算座標Xi(i=1〜end)各々に対応する最適音速値及び処理済素子データが算出され、画像生成部26の整相加算部38に移行する。こうして、計算座標X1〜Xendの各々において算出された最適音速値及び最適音速値を用いて素子データ処理部22で処理された処理済素子データは、計算座標X1〜Xendの各々の座標データと共に整相加算部38に供給される。
次に、画像生成部26の整相加算部38では、素子データに受信フォーカス処理を施して受信データ(音線信号)を生成し、検波処理部40に供給する。検波処理部40では、音線信号を処理してBモード画像信号を生成する。Bモード画像信号を、DSC42がラスター変換し、画像作成部44が画像処理を施し、超音波画像が生成される。生成された超音波画像は、画像メモリ46に格納されると共に、表示制御部27により超音波画像が表示部28に表示される。
このように、本発明の超音波検査装置は、素子データ処理部でマルチライン処理され、ゴーストの信号が弱減され、真の信号が強調された処理済素子データを用いて作成した超音波画像を評価する画質指標を用いて最適な環境音速を求め、求められた最適な環境音速を用いてマルチライン処理された処理済素子データを用いて超音波画像を作成するので、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができ、また、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高解像度で、SN比が高く、最適な空間分解能を持つシャープな超音波画像を得ることができる。
次に、図9及び図10に基づいて、本発明の実施例2の超音波検査装置について説明する。
図9は、本発明の実施例2に係る超音波検査装置の構成を概念的に示すブロック図であり、図10は、図9に示す超音波検査装置の要部の一例を処理フローに沿って概念的に示すブロック図である。
なお、図9及び図10に示す超音波検査装置100は、図1及び図2に示す超音波検査装置10と、第2の音速補正部102を備えている点で異なる以外は、全く同様の構成を有するものであるので、同一の構成要素には、同一の参照符号を付し、その詳細な説明は省略する。
図9に示すように、超音波検査装置100は、超音波プローブ12と、超音波プローブ12に接続される送信部14及び受信部16と、A/D変換部18と、素子データ保持部20と、第2の音速補正部102と、素子データ処理部22を備える音速補正部24と、画像生成部26と、表示制御部27と、表示部28と、制御部30と、操作部32と、格納部34とを有する。
超音波検査装置100では、音速補正部24による音速補正を行う前に、第2の音速補正部102にて、未処理素子データ(第1の素子データ)を用いて音速補正を行って暫定的な音速値を求め、この音速値を、音速補正部24で音速補正を行う際の初期値として利用して最適音速値を求めるものである。
第2の音速補正部102は、素子データ保持部20と音速補正部24との間に接続されており、素子データ保持部20に保存された未処理素子データ(第1の素子データ)を用いて、音速補正を行い算出した音速値を音速補正部24に供給する。
第2の音速補正部102が行う音速補正には、特に限定はないが、例えば、特開2011−92686号に記載の方法が利用できる。この方法は、複数の設定音速毎に未処理素子データに対して受信フォーカス処理を行って設定音速毎にフォーカス指標を算出し、このフォーカス指標に基づいて最適な音速値を決定するものである。すなわち、未処理素子データを最も適切に受信フォーカス処理できる設定音速を最適な音速値とする。
なお、第2の音速補正部102での音速補正は、音速補正部24における音速補正の際の音速探索の初期値として大まかに求めればよいので、必ずしも計算領域や着目領域毎に送信フォーカスを設定する必要はない。また、計算領域や着目領域毎に行うものに限定はされず、画面全体を1つの着目領域として音速値を求めてもよい。
図10に示すように、第2の音速補正部102は、求めた音速値を音速補正部24の最適音速決定部23の音速探索範囲設定部74に供給する。
最適音速決定部23の音速探索範囲設定部74は、第2の音速補正部102から供給された音速値を初期音速値Vstとして音速値を設定する。
音速値更新部82は、初期音速値Vst及びこの前後の音速値(Vst±ΔV)を順次設定して、画質指標算出部78で画質指標を算出し、3つの音速値での結果に基づいて、次の設定音速Vを設定する。例えば、設定音速VがVst+ΔVの場合が最も良い結果であれば、音速値更新部82は、次の設定音速VとしてVst+2×ΔVを設定し、設定音速VがVstの場合が最も良い結果であれば、Vstを最適音速値として採用する。すなわち、最適音速決定部23は、第2の音速補正部102から供給された音速値を初期音速値として探索的に最適音速を求める。
第2の音速補正部102が未処理素子データを用いて音速補正を行い、音速補正部24での音速補正の初期値として用いる構成とすることにより、最適音速値の測定・演算にかかる時間を短縮することができる。
以上、本発明について詳細に説明したが、本発明は上記実施例に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
例えば、図1及び図2に示す各構成要素をハードウエアとして構成いても良いし、コンピュータ等で実行されるソフトウエアとして構成しても良い。
更に、本発明の各実施形態について説明した超音波検査装置の各機能に対応して、例えば、上述した音速補正部に含まれる素子データ処理部を始めとする各構成要素、及び画像生成部として機能するようにコンピュータを動作させる超音波検査ためのプログラムも、また、上記超音波検査方法の各ステップ、例えば、上述したマルチライン処理の各ステップを含む音速補正の各ステップ、及び画像生成の各ステップを手順としてコンピュータに実行させる超音波検査のためのプログラムも、本発明の実施の形態の1つである。また、そのようなプログラムを記録したコンピュータ読み取り可能な記憶媒体も、本発明の実施の形態の1つである。
10、100 超音波検査装置
12 超音波プローブ
14 送信部
16 受信部
18 A/D変換部
20 素子データ保持部
22 素子データ処理部
23 最適音速決定部
24 音速補正部
26 画像生成部
27 表示制御部
28 表示部
30 制御部
32 操作部
34 格納部
36 探触子
38、38a 整相加算部
40、40a 検波処理部
42 DSC
44、44a 画像作成部
46 画像メモリ
48 遅延時間算出部
50 重ね合わせ処理部
70 計算領域設定部
72 初期計算座標設定部
74 音速探索範囲設定部
76 初期音速値設定部
78 画質指標算出部
80 最終音速値判定部
82 音速値更新部
84 最適音速判定部
86 音速設定部
88 最終計算座標判定部
90 計算座標更新部
102 第2の音速補正部

Claims (15)

  1. 超音波ビームを用いて検査対象物を検査する超音波検査装置であって、
    前記超音波ビームを送信し、かつ、前記検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子と、
    前記探触子に、複数の前記素子を用い、所定の送信焦点を形成するように前記超音波ビームを送信させることを、複数回、行わせる送信部と、
    個々の前記超音波ビームの送信に対応して、複数の前記素子が出力したアナログ素子信号を受け、所定の処理を施す受信部と、
    前記受信部が処理したアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとするAD変換部と、
    前記検査対象物内の音速値を決定する音速決定部と、
    前記検査対象物内の音速値を用いて2つ以上の前記第1の素子データから、前記2つ以上の第1の素子データのいずれかに対応する第2の素子データを生成する素子データ処理部と、を有し、
    前記音速決定部は、前記素子データ処理部において前記2つ以上の第1の素子データから前記第2の素子データを作成する時に用いる前記音速値を最適化して最適な音速値として求めることを特徴とする超音波検査装置。
  2. 前記音速決定部は、前記素子データ処理部を有し、前記素子データ処理部における前記2つ以上の第1の素子データの遅延時間補正計算に用いる前記音速値を最適化して最適な音速として求める請求項1に記載の超音波検査装置。
  3. 前記音速決定部は、前記検査対象物内を複数の領域に分割したときの各領域における音速値を決定する請求項1又は2に記載の超音波検査装置。
  4. 更に、前記2つ以上の第1の素子データを含む一部の第1の素子データ、あるいは全ての第1の素子データを保持する素子データ保持部を有する請求項1〜3のいずれか1項に記載の超音波検査装置。
  5. 前記素子データ保持部は、前記送信部により、前記検査対象物内における少なくとも2つ以上の重なり合う対象領域毎において、前記超音波ビームを送信し、前記受信部により、前記超音波ビームによって前記検査対象物内における少なくとも2つ以上の重なり合う対象領域毎に発生した前記超音波エコー信号を受信して生成された、素子データ毎に受信した各素子における受信時間毎の受信データを含む前記2以上の第1の素子データを保持し、
    前記素子データ処理部は、前記素子データ毎の前記超音波ビームを送信したときの素子の幾何学的な配置及び前記素子データ毎の前記受信データの各素子の幾何学的な配置に基づいて、前記2以上の第1の素子データから前記第2の素子データを生成する請求項4に記載の超音波検査装置。
  6. 前記送信部は、前記超音波ビームの送信を、中心となる素子を変更して、複数回、行う請求項5に記載の超音波診断装置。
  7. 前記素子データ処理部は、前記検査対象物内の音速値を用いて前記2つ以上の第1の素子データの各遅延時間を算出する遅延時間算出部と、前記2つ以上の第1の素子データを算出された各遅延時間及び受信された前記探触子の素子の位置に基づいて重ね合わせ、前記第2の素子データを生成する重ね合わせ処理部とを有する請求項1〜6のいずれか1項に記載の超音波検査装置。
  8. 更に、前記素子データ処理部で生成された前記第2の素子データに基づいて超音波画像を生成する画像生成部と、を有し、
    前記素子データ処理部は、前記画像生成部で生成されるべき超音波画像上の1ラインの前記第2の素子データを作成する請求項1〜7のいずれか1項に記載の超音波検査装置。
  9. 前記音速決定部は、
    前記素子データ処理部と、
    前記素子データ処理部における前記遅延時間算出に用いられる前記検査対象物の音速値を変更する音速変更部と、
    前記音速変更部で変更された音速値を用いて前記第2の素子データを作成する前記素子データ処理部と、
    前記素子データ処理部により作成された前記第2の素子データを、前記音速変更部で変更された音速値を用いて整相加算して音線信号を生成する整相加算部と、
    前記整相加算部で生成された前記音線信号に基づいて超音波画像を生成する画像生成部と、
    前記画像生成部で生成された前記超音波画像の画質指標を算出する画質指標算出部と、
    前記画質指標算出部で算出された前記超音波画像の画質指標を判定する判定部と、
    前記判定部による前記画質指標の判定結果に応じて前記最適な音速値を設定する音速設定部と、を有する請求項1〜8のいずれか1項に記載の超音波検査装置。
  10. 前記音速決定部は、
    前記音速変更部が、前記音速値を初期音速値から最終音速値まで所定音速間隔で変更する度毎に、前記素子データ処理部が変更された当該音速値を用いて前記第2の素子データを作成し、前記整相加算部が同じ当該音速値を用いて前記第2の素子データを整相加算して前記音線信号を生成し、前記画像生成部が前記音線信号に基づいて前記超音波画像を生成し、前記画質指標算出部が前記超音波画像の前記画質指標を算出することを繰り返して、前記音速変更部によって変更された全ての複数の前記音速値にそれぞれ対応する複数の前記画質指標を算出しておき、
    前記判定部が、前記画質指標算出部で全ての複数の前記音速値毎に算出された複数の前記画質指標の中から最適な画質指標を判定し、前記音速設定部が、前記判定部で判定された前記最適な画質指標に対応する音速値を前記最適な音速値として設定する請求項9に記載の超音波検査装置。
  11. さらに、前記音速決定部よりも前段に、前記第1の素子データを用いて最適な音速値を算出する第2の音速決定部を有し、
    前記音速決定部は、前記第2の音速決定部で算出した音速値を初期値として用いる請求項1〜10のいずれか1項に記載の超音波検査装置。
  12. 超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子を用いて検査対象物を検査する超音波検査方法であって、
    前記探触子によって、複数の前記素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行い、かつ、個々の前記超音波ビームの送信に対応して、複数の前記素子がアナログ素子信号を出力し、
    このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとし、
    前記検査対象物内の音速値を用いて、2つ以上の前記第1の素子データから、前記第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる前記音速値を最適化して最適な音速値として決定することを特徴とする超音波検査方法。
  13. 前記最適な音速を決定するステップは、
    前記音速値を初期音速値から最終音速値まで所定音速間隔で変更する毎に、変更された当該音速値を用いて前記第2の素子データを作成し、同じ当該音速値を用いて前記第2の素子データを整相加算して前記音線信号を生成し、生成された前記音線信号に基づいて前記超音波画像を生成し、前記画質指標算出部が前記超音波画像の前記画質指標を算出することを繰り返して、変更された全ての複数の前記音速値毎にそれぞれ対応する複数の前記画質指標を算出しておき、
    全ての複数の前記音速値毎に算出された複数の前記画質指標の中から最適な画質指標を判定し、判定された前記最適な画質指標に対応する音速値を前記最適な音速値として設定する請求項12に記載の超音波検査方法。
  14. 超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して検査対象物を検査するために、複数の素子が配列された探触子によって、複数の前記素子を用い、所定の送信焦点を形成するように、超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の前記素子がアナログ素子信号を出力し、このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データを複数、取得するステップと、
    前記検査対象物内の音速値を用いて、得られた2つ以上の前記第1の素子データから前記第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる前記音速値を最適化して最適な音速値として決定するステップとを手順としてコンピュータに実行させることを特徴とするプログラム。
  15. 請求項14に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
JP2013125337A 2012-09-28 2013-06-14 超音波検査装置、超音波検査方法、プログラム及び記録媒体 Active JP5946427B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2013125337A JP5946427B2 (ja) 2012-09-28 2013-06-14 超音波検査装置、超音波検査方法、プログラム及び記録媒体
CN201380050604.0A CN104684487A (zh) 2012-09-28 2013-09-24 超声波检查装置、超声波检查方法、程序及记录介质
PCT/JP2013/075666 WO2014050797A1 (ja) 2012-09-28 2013-09-24 超音波検査装置、超音波検査方法、プログラム及び記録媒体
US14/671,236 US10687786B2 (en) 2012-09-28 2015-03-27 Ultrasound inspection apparatus, ultrasound inspection method and recording medium

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012215721 2012-09-28
JP2012215721 2012-09-28
JP2013125337A JP5946427B2 (ja) 2012-09-28 2013-06-14 超音波検査装置、超音波検査方法、プログラム及び記録媒体

Publications (2)

Publication Number Publication Date
JP2014079565A true JP2014079565A (ja) 2014-05-08
JP5946427B2 JP5946427B2 (ja) 2016-07-06

Family

ID=50388196

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013125337A Active JP5946427B2 (ja) 2012-09-28 2013-06-14 超音波検査装置、超音波検査方法、プログラム及び記録媒体

Country Status (4)

Country Link
US (1) US10687786B2 (ja)
JP (1) JP5946427B2 (ja)
CN (1) CN104684487A (ja)
WO (1) WO2014050797A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102015004318A1 (de) 2014-04-08 2015-10-08 Fanuc Corporation Servo-Steuergerät zur Verringerung eines Synchronisierungsfehlers bei einer synchronen Bearbeitung

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9798856B2 (en) * 2012-03-21 2017-10-24 Koninklijke Philips N.V. Clinical workstation integrating medical imaging and biopsy data and methods using same
JP6006249B2 (ja) * 2014-03-24 2016-10-12 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
KR20150118732A (ko) * 2014-04-15 2015-10-23 삼성전자주식회사 초음파 장치 및 그 제어 방법
CN104188684B (zh) * 2014-09-15 2016-08-31 声泰特(成都)科技有限公司 一种自适应医疗超声成像声速优化和信号修正方法及系统
WO2019044075A1 (ja) * 2017-08-31 2019-03-07 富士フイルム株式会社 画像生成装置および作動方法
CN110680398A (zh) * 2019-10-11 2020-01-14 深圳开立生物医疗科技股份有限公司 一种超声回波成像方法、系统及超声设备
CN115867202A (zh) * 2020-07-09 2023-03-28 深圳迈瑞生物医疗电子股份有限公司 超声声速确定方法和超声成像系统
WO2022117393A1 (en) * 2020-12-01 2022-06-09 Koninklijke Philips N.V. Retrospective transmit focusing using a transmit speed systems, devices, and methods
US20230061869A1 (en) * 2021-08-26 2023-03-02 GE Precision Healthcare LLC System and methods for beamforming sound speed selection
US11867807B2 (en) * 2021-09-01 2024-01-09 GE Precision Healthcare LLC System and methods for beamforming sound speed selection

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004522515A (ja) * 2001-01-26 2004-07-29 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 空間的複合のための適応画像処理
JP2009090102A (ja) * 2007-09-18 2009-04-30 Fujifilm Corp 超音波診断方法及び装置
JP2010082190A (ja) * 2008-09-30 2010-04-15 Fujifilm Corp 超音波信号処理装置及び方法
US20120238874A1 (en) * 2011-03-18 2012-09-20 Fujifilm Corporation Ultrasound diagnostic apparatus and ultrasound image producing method

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4893284A (en) * 1988-05-27 1990-01-09 General Electric Company Calibration of phased array ultrasound probe
JPH0595946A (ja) * 1991-10-09 1993-04-20 Yokogawa Medical Syst Ltd 超音波の音速測定法及び音速測定手段を備えた超音波診断装置
EP1974672B9 (en) * 2007-03-28 2014-04-16 Kabushiki Kaisha Toshiba Ultrasonic imaging apparatus and ultrasonic velocity optimization method
JP2009090104A (ja) 2007-09-18 2009-04-30 Fujifilm Corp 超音波診断方法及び装置
JP2009089940A (ja) * 2007-10-10 2009-04-30 Toshiba Corp 超音波診断装置
JP5355924B2 (ja) 2008-03-31 2013-11-27 株式会社東芝 超音波診断装置
JP2010099452A (ja) * 2008-09-25 2010-05-06 Fujifilm Corp 超音波診断装置及び超音波診断方法
JP5566773B2 (ja) * 2009-06-30 2014-08-06 株式会社東芝 超音波診断装置および音速設定方法
JP5486257B2 (ja) * 2009-09-28 2014-05-07 富士フイルム株式会社 超音波診断装置及び弾性指標算出方法
JP5292440B2 (ja) * 2011-06-03 2013-09-18 富士フイルム株式会社 超音波診断装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004522515A (ja) * 2001-01-26 2004-07-29 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 空間的複合のための適応画像処理
JP2009090102A (ja) * 2007-09-18 2009-04-30 Fujifilm Corp 超音波診断方法及び装置
JP2010082190A (ja) * 2008-09-30 2010-04-15 Fujifilm Corp 超音波信号処理装置及び方法
US20120238874A1 (en) * 2011-03-18 2012-09-20 Fujifilm Corporation Ultrasound diagnostic apparatus and ultrasound image producing method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102015004318A1 (de) 2014-04-08 2015-10-08 Fanuc Corporation Servo-Steuergerät zur Verringerung eines Synchronisierungsfehlers bei einer synchronen Bearbeitung

Also Published As

Publication number Publication date
US10687786B2 (en) 2020-06-23
CN104684487A (zh) 2015-06-03
US20150196274A1 (en) 2015-07-16
JP5946427B2 (ja) 2016-07-06
WO2014050797A1 (ja) 2014-04-03

Similar Documents

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

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141016

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151013

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151202

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160531

R150 Certificate of patent or registration of utility model

Ref document number: 5946427

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