JP6596013B2 - 肋骨間空間を用いたコヒーレントな複合による音響撮像のためのシステム及び方法 - Google Patents

肋骨間空間を用いたコヒーレントな複合による音響撮像のためのシステム及び方法 Download PDF

Info

Publication number
JP6596013B2
JP6596013B2 JP2016559411A JP2016559411A JP6596013B2 JP 6596013 B2 JP6596013 B2 JP 6596013B2 JP 2016559411 A JP2016559411 A JP 2016559411A JP 2016559411 A JP2016559411 A JP 2016559411A JP 6596013 B2 JP6596013 B2 JP 6596013B2
Authority
JP
Japan
Prior art keywords
transducer array
acoustic transducer
transmit
acoustic
imaging region
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2016559411A
Other languages
English (en)
Other versions
JP2017509416A (ja
JP2017509416A5 (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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2017509416A publication Critical patent/JP2017509416A/ja
Publication of JP2017509416A5 publication Critical patent/JP2017509416A5/ja
Application granted granted Critical
Publication of JP6596013B2 publication Critical patent/JP6596013B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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/52047Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • 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/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/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
    • G01S15/8927Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array using simultaneously or sequentially two or more subarrays or subapertures
    • 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/8977Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • 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
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection

Landscapes

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

Description

本発明は、音響撮像装置、システム、及び方法に関する。
音響(例えば超音波)撮像(acoustic imaging)、又はその他のアクティブ音響映像法においては、まずトランスデューサのアレイが、通常は集束ビームを発生させるために適した遅延で、パルスを媒体内に送信し、次いでエコーを受信する。ビーム形成は媒体の画像を再構築するために用いられる。画像の品質は、点像分布関数(PSF:point−spread function)すなわち単一ポイントに対する撮像システムの応答に関連している。PSFはアレイのアパーチャに関連している。遠距離近似では、PSFの横プロファイル(lateral profile)は、アパーチャ振幅分布のフーリエ変換の2乗である。アポディゼーション(apodization)を用いない場合、これはsinc関数となる。
しかしながら、アパーチャの一部が1つ又はいくつかの障害物によって遮られることがある。医療用超音波では、この障害物は肋骨である場合がある。例えば胸部による心臓の超音波撮像(経胸壁超音波撮像)は、アパーチャが大きすぎる場合にその一部を遮る肋骨が存在すると、複雑になる。遮られたアパーチャ部分がアパーチャ縁部に沿っている場合、結果として解像度が損なわれる。その部分がアパーチャ内にある場合、結果として画像品質を大きく低下させるグレーティングローブ(grating lobe)が生じる。これと同じ現象はアパーチャのギャップからも起こり、これは、機能していないトランスデューサアレイ要素によって、又はギャップで分離された2つの超音波アレイによりコヒーレントに撮像を行おうとした場合に生じる。一方、2本の肋骨間の空間(以降、これを「肋骨間ギャップ」と呼ぶ)に収まる小さいアパーチャを用いる場合、視野及び超音波画像の解像度が限定される。
従って、いくつかの肋骨間空間にまたがる大きいコヒーレントなアパーチャを用いることができる超音波撮像のための方法及びシステムを提供することが望ましい。
本発明の1つの態様において、方法は、音響トランスデューサアレイを用いて、音響トランスデューサアレイと撮像領域の少なくとも一部との間に1つ以上の障害物がある撮像領域の画像データを生成することと、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償することと、補償した画像データから撮像領域の画像を生成することと、を有する。
いくつかの実施形態では、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して1つ以上の障害物による撮像領域の画像データの欠落を補償すること、及び補償した画像データから撮像領域の画像を生成することが、音響トランスデューサアレイのための逆フィルタを求めることであって、逆フィルタを、1つ以上の障害物が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトルで乗算すると、1つ以上の障害物がないと存在するはずの撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトルが生成される、逆フィルタを求めることと、音響トランスデューサアレイの各送信/受信要素対によって生成した信号を、送信/受信対の角度周波数に対応する逆フィルタの値で重み付けすることと、全ての送信受信対の重み付けした信号を合計することによって音響画像を生成することと、を有する。
いくつかの実施形態では、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償すること、及び補償した画像データから撮像領域の画像を生成することが、少なくとも第1及び第2のアポディゼーション関数を用いて少なくとも第1及び第2の送信及び受信動作を実行することであって、音響トランスデューサアレイを用いて、第1の音響波を撮像領域に送信し、撮像領域から第1の音響エコーを受信し、そこから第1の画像データを生成することによって第1の送信及び受信動作を実行することであって、第1のアポディゼーション関数を音響トランスデューサアレイに適用して第1の送信及び受信動作のための第1の送信アパーチャ及び第1の受信アパーチャを生成する、第1の送信及び受信動作を実行することと、音響トランスデューサアレイを用いて、第2の音響波を撮像領域に送信し、撮像領域から第2の音響エコーを受信し、そこから第2の画像データを生成することによって第2の送信及び受信動作を実行することであって、第2のアポディゼーション関数を音響トランスデューサアレイに適用して第2の送信及び受信動作のための第2の送信アパーチャ及び第2の受信アパーチャを生成する、第2の送信及び受信動作を実行することと、を含む、第1及び第2の送信及び受信動作を実行することと、第1の画像データを第2の画像データと組み合わせることによって撮像領域の画像を生成することと、を有する。少なくとも第1及び第2のアポディゼーション関数は、撮像領域に対する音響トランスデューサアレイの逆フィルタから求められ、逆フィルタを、1つ以上の障害物が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトルで乗算すると、1つ以上の障害物がないと存在するはずの撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトルが生成される。
いくつかの実施形態では、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償すること、及び補償した画像データから撮像領域の画像を生成することが、RFデータ(検出前)に基づく又は画像に基づくデコンボリューションアルゴリズムを実行することを有する。
本発明の別の態様では、撮像領域を撮像するための装置は、音響トランスデューサアレイと撮像領域の少なくとも一部との間に1つ以上の障害物がある撮像領域の画像データを生成する音響トランスデューサアレイと、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償し、補償した画像データから撮像領域の画像を生成する1つ以上のプロセッサと、を備える。
いくつかの実施形態では、1つ以上のプロセッサが、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して1つ以上の障害物による撮像領域の画像データの欠落を補償し、補償した画像データから撮像領域の画像を生成することが、音響トランスデューサアレイのための逆フィルタを求めることであって、逆フィルタを、1つ以上の障害物が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトルで乗算すると、1つ以上の障害物がないと存在するはずの撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトルが生成される、逆フィルタを求めることと、音響トランスデューサアレイの各送信/受信要素対によって生成した信号を、送信/受信対の角度周波数に対応する逆フィルタの値で重み付けすることと、及び、全ての送信受信対の重み付けした信号を合計することによって音響画像を生成することと、によって行われる。
いくつかの実施形態では、1つ以上のプロセッサが、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償すること、及び補償した画像データから撮像領域の画像を生成することが、少なくとも第1及び第2のアポディゼーション関数を用いて少なくとも第1及び第2の送信及び受信動作を実行することであって、音響トランスデューサアレイを用いて、第1の音響波を撮像領域に送信し、撮像領域から第1の音響エコーを受信し、そこから第1の画像データを生成することによって第1の送信及び受信動作を実行することであって、第1のアポディゼーション関数を音響トランスデューサアレイに適用して第1の送信及び受信動作のための第1の送信アパーチャ及び第1の受信アパーチャを生成する、第1の送信及び受信動作を実行することと、音響トランスデューサアレイを用いて、第2の音響波を撮像領域に送信し、撮像領域から第2の音響エコーを受信し、そこから第2の画像データを生成することによって第2の送信及び受信動作を実行することであって、第2のアポディゼーション関数を音響トランスデューサアレイに適用して第2の送信及び受信動作のための第2の送信アパーチャ及び第2の受信アパーチャを生成する、第2の送信及び受信動作を実行することと、を含む、第1及び第2の送信及び受信動作を実行することと、第1の画像データを第2の画像データと組み合わせることによって撮像領域の画像を生成することと、によって行われる。少なくとも第1及び第2のアポディゼーション関数は、撮像領域に対する音響トランスデューサアレイの逆フィルタから求められ、逆フィルタを、1つ以上の障害物が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトルで乗算すると、1つ以上の障害物がないと存在するはずの撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトルが生成される。
いくつかの実施形態では、1つ以上のプロセッサが、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償すること、及び補償した画像データから撮像領域の画像を生成することが、RFデータ(検出前)に基づく又は画像に基づくデコンボリューションアルゴリズムを実行することによって行われる。
音響トランスデューサと撮像領域の少なくとも一部との間の肋骨の存在により生成される音響画像アパーチャにおける肋骨間ギャップを示す。 図1に示すような状況において生成される音響画像においてメインローブの両側に顕著なグレーティングローブが観察されることを示す。 音響撮像装置の一実施形態を示す。 撮像デバイスと撮像領域の少なくとも一部との間に1つ以上の障害物が存在する撮像領域を撮像する方法400の一実施形態を示すフローチャートである (1)1つ以上の障害物がないと存在するはずの、撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトル、(2)1つ以上の障害物がある撮像領域に対する音響トランスデューサアレイの角度スペクトルの例を示す。 撮像デバイスと撮像領域の少なくとも一部との間に1つ以上の障害物がある撮像領域を撮像する方法の別の実施形態を示すフローチャートである。 肋骨モデルを示し、どのように光線追跡アルゴリズムを用いてどの画像ポイントがトランスデューサアレイ内のどのトランスデューサ要素を見るかを検出するかを説明する。 2つの送信及び受信動作で用いられる2組の送信/受信アパーチャのためのアポディゼーション関数を示す。 撮像デバイスと撮像領域の少なくとも一部との間に1つ以上の障害物がある撮像領域を撮像する方法の更に別の実施形態を示すフローチャートである。 補償アルゴリズムを肋骨間空間に用いた場合の、図2に示したものと同じ状況での点像分布関数を示す。 補償なしの場合の、大きいアパーチャにおいて肋骨間ギャップが存在するシミュレーション画像を示す。 補償なしの場合の、大きいアパーチャにおいて肋骨間ギャップが存在するシミュレーション画像を示す。 2つの妨害物(blockage)(例えば2本の肋骨)をシミュレーションした場合の、大きいアパーチャを用いて生成された経胸壁心エコー図を示す。 2つの妨害物(例えば2本の肋骨)をシミュレーションし、肋骨間ギャップ補償アルゴリズムを用いた場合の、大きいアパーチャを用いて生成された経胸壁心エコー図を示す。 図12Aと同様であるが妨害物が存在しない場合の、大きいアパーチャを用いて生成された基準経胸壁心エコー図を示す。 図12A〜図12Cを生成するために用いたものよりも小さいアパーチャを用いて生成された経胸壁心エコー図を示す。 肋骨が存在する生成されたコヒーレント画像を示す。 「肋骨間ギャップ補償」アルゴリズムを用いるコヒーレントな複合(coherent compounding)を利用した場合の、図13に示すものと同様の肋骨が存在する生成された画像を示す。 肋骨が存在しないと生成されるはずの、図13及び図14Aに示すものと同様の画像を示す。 図8に示したモデルにおける遠距離中心のポイントの角度スペクトルを示す。 図8に示したモデルにおける画像の縁部上の、音響トランスデューサアレイに近いポイントの角度スペクトルを示す。 図8の例について、画像の様々なポイントについてアルゴリズムにより再構築され得る最大往復アパーチャサイズを示す。 肋骨が肋骨間空間と同じくらい大きい場合の補償アルゴリズムの結果を示す。
これより、本発明の好適な実施形態が示された添付図面を参照して本発明について更に充分に説明する。しかしながら本発明は異なる形態で具現化することも可能であり、本明細書に記載する実施形態に限定されるものとして解釈するべきではない。これらの実施形態は本発明の教示の例として与えられる。
図1は、音響トランスデューサと撮像領域の少なくとも一部との間に配置された障害物(例えば肋骨)の存在によって生じ得る、音響画像アパーチャにおける肋骨間ギャップを示す。
具体的には、図1は音響プローブ10(例えば超音波プローブ)を示し、これは、例えば心臓のような関心領域(ROI:region of interest)5を含む撮像領域7を撮像するための音響トランスデューサアレイ10’を含む。撮像領域7は、例えば肋骨のような複数の障害物15−1、15−2、15−3も含み、これらは肋骨間空間11−1、11−2により相互に分離され離間している。一般に、音響プローブ10及び音響トランスデューサアレイ10’は、音響撮像装置(例えば超音波撮像システム)の構成要素である。音響トランスデューサアレイ10’は複数の音響要素(例えば超音波要素)を含む。音響トランスデューサアレイ10’の音響要素は、例えば音響撮像装置のマイクロビームフォーマから電気信号を受信し、これに応答して対応する音響信号を撮像領域7に送信することによって、送信動作を実行するように構成されている。また、音響トランスデューサアレイ10’の音響要素は、撮像領域7から音響信号(例えば音響エコー)を受信し、これに応答して対応する電気信号を音響撮像装置の信号処理回路に供給することによって、受信動作を実行するように構成することができる。音響プローブ10、音響トランスデューサアレイ10’、音響プローブ10と音響トランスデューサアレイ10’とを含む音響撮像装置の例示的な実施形態のこれ以上の詳細については、図3の説明を参照して以下に与える。
動作において、アパーチャ110−1及び110−2が音響トランスデューサアレイ10’によって形成され、撮像領域7の、特にROI5の対応するセクタスキャン120及び122を形成する。図1からわかるように、この例では、音響トランスデューサアレイ10’とROI5との間に障害物15−1、15−2、15−3が配置されている。従って、スキャン120は、障害物15−1によって遮られるか又は影になる第1の部分120−1、肋骨間空間11−1を通りアパーチャ110−1の視野を含む第2の部分120−2、障害物15−2によって遮られるか又は影になる第3の部分120−3を含む。同様に、スキャン122は、障害物15−2によって遮られるか又は影になる第1の部分122−1、肋骨間空間11−2を通りアパーチャ110−2の視野を含む第2の部分122−2、障害物15−3によって遮られるか又は影になる第3の部分122−3を含む。
音響撮像装置はビーム形成を用いて撮像領域7の画像を再構築する。画像の品質は、点像分布関数(PSF)すなわち単一ポイントに対する音響撮像装置の応答に関連している。PSFは、音響トランスデューサアレイ10’のアパーチャ(複数のアパーチャ)110−1、110−2等に関連している。遠距離近似では、PSFの横プロファイルはアパーチャ振幅分布のフーリエ変換の2乗である。アパーチャ(複数のアパーチャ)にアポディゼーションを適用しない場合、PSFはsinc関数となる。しかしながら、図1に示すように、アパーチャ(複数のアパーチャ)の一部は1つ又はいくつかの障害物によって遮られている。遮られたアパーチャ部分がアパーチャ縁部に沿っている場合、結果として解像度が損なわれる。その部分がアパーチャ内にある場合、結果として画像品質を大きく低下させるグレーティングローブが生じる。これと同じ現象はアパーチャのギャップからも起こり、これは、機能していない要素によって、又はギャップで分離された2つのプローブアレイによりコヒーレントに撮像を行おうとした場合に生じる。
図2は、図1に示し上述したような、1つ以上のアパーチャによる撮像領域の1つ以上の部分の撮像を障害物が遮っている状況において、音響プローブ10により生成された音響画像200のメインローブ202の両側に顕著なグレーティングローブ204が観察されることを示す。具体的には、図2は、200要素アパーチャを3つの42要素サブアパーチャに分割し、2つの37要素ギャップで分離した場合のPSFシミュレーションを示す(この場合、不活性の(inactive)要素によりシミュレーションされている)。シミュレーションされたファントムのポイントは、4cm、8cm、12cm、及び16cmの深さにあり、画像の上から下に深くなる。
本発明者らは、ギャップ又は妨害物が大きいアパーチャをいくつかのサブアパーチャに分離している場合、ギャップがサブアパーチャよりも小さい限り、補償アルゴリズムを用いることによって完全な非遮蔽アパーチャの解像度を回復できることを発見し認識している。以下に記載する様々な実施形態では、補償アルゴリズムは往復信号の冗長性を利用する。具体的には、以下に記載する様々な実施形態は、往復超音波信号の冗長性の概念に基づいている[1]。NTxとNRxをそれぞれ、音響トランスデューサアレイにおける送信要素と受信要素の指数とする。遠距離において(かつ位相収差が存在しない場合)、これらの和(NTx+NRx)が一定であるような全てのTx/Rx要素の対の信号は同一である。従って、これらの対のいくつかが肋骨に遮られたために欠落した場合、遮られていない同一の対によって置換することができる。
図3は、往復超音波信号の冗長性を利用してトランスデューサアレイアパーチャ(複数のアパーチャ)のギャップ及び妨害物を補償するための1つ以上の補償方法又はアルゴリズムを実施することができる音響撮像装置300の一実施形態を示す。音響撮像装置300では、音響トランスデューサアレイ10’(例えば静電容量型超音波トランスデューサ(CMUT:capacitive micromachined ultrasonic transducer)アレイ)は、超音波を送信しエコー情報を受信するための音響プローブ10に設けられている。あるいは、トランスデューサアレイ10’は、PZT又はPVDFのような材料で形成された圧電トランスデューサ要素を含み得る。トランスデューサアレイ10’は、2D平面で及び/又は3D撮像のため3次元でのスキャンが可能な、トランスデューサ要素の1次元又は2次元アレイである。トランスデューサアレイ10’は、音響プローブ10において、CMUTアレイセル又は圧電要素によって信号の送信及び受信を制御するマイクロビームフォーマ12に結合されている。マイクロビームフォーマは、トランスデューサ要素のグループ又は「パッチ」により受信された信号の少なくとも部分的なビーム形成を行うことができる。これについては、米国特許第5,997,479号(Savord等)、第6,013,032号(Savord)、第6,623,432号(Powers等)に記載されている。マイクロビームフォーマ12は、プローブケーブルによって送受信(T/R)スイッチ16に結合されている。T/Rスイッチ16は、送信と受信との間で切り換えを行い、マイクロビームフォーマ12が使用されずにトランスデューサアレイ10’が直接メインシステムビームフォーマ20により動作される場合にメインビームフォーマ20を高エネルギ送信信号から保護する。マイクロビームフォーマ12の制御下でのトランスデューサアレイ10からの超音波ビーム送信は、T/Rスイッチによってマイクロビームフォーマに結合されたトランスデューサコントローラ18と、ユーザインタフェース又は制御パネル38のユーザ動作による入力を受信するメインシステムビームフォーマ20とによって誘導される。トランスデューサコントローラ18が制御する機能の1つは、ビームをステアリングし集束させる方向である。ビームは、トランスデューサアレイ10’から(これに直交して)まっすぐに、又はより広い視野のために様々な角度にステアリングされる。トランスデューサコントローラ18は、CMUTアレイのためDCバイアス制御45を制御するように結合することができる。DCバイアス制御45は、CMUTセルに印加されるDCバイアス電圧(複数の電圧)を制御する。
受信されてマイクロビームフォーマ12により部分的にビーム形成された信号はメインビームフォーマ20に結合され、ここで、個々のトランスデューサ要素パッチからの部分的にビーム形成された信号は組み合わされて、完全にビーム形成された信号となる。例えば、メインビームフォーマ20は128のチャネルを有し、その各々が、CMUTトランスデューサセル又は圧電要素のパッチから部分的にビーム形成された信号を受信する。いくつかの実施形態では、各パッチは数十又は数百のCMUTトランスデューサセル又は圧電要素を含む。このように、トランスデューサアレイ10’の多数のトランスデューサ要素(例えば数十、数百、又は数千の要素)によって受信された信号は、ビーム形成された単一の信号を効率的に生成することができる。
ビーム形成された信号は、信号プロセッサ22に結合される。信号プロセッサ22は、受信したエコー信号を様々な手法で処理する。その手法は例えば、バンドパスフィルタリング、デシメーション、I成分とQ成分の分離、高調波信号分離である。高調波信号分離は、線形信号と非線形信号を分離して、組織及びマイクロバブルから戻った非線形(基本周波数の高調波)エコー信号の識別を可能とするように作用する。信号プロセッサは、スペックル低減、信号複合、及びノイズ除去のような追加の信号向上を実行してもよい。信号プロセッサのバンドパスフィルタは、エコー信号が深い場所から受信されると通過帯域が高周波帯域から低周波帯域へとスライドするトラッキングフィルタとし、これによって、解剖学上の情報が存在しない深い場所からの高周波数のノイズを阻止することができる。
処理された信号は、Bモードプロセッサ26及びドップラープロセッサ28に結合される。Bモードプロセッサ26は、身体の臓器の組織や血管のような身体内の構造を撮像するため、受信した音響信号の振幅の検出を用いる。身体の構造のBモード画像は、高調波画像モード又は基本画像モード又はこれら双方の組み合わせにおいて形成される。これは、米国特許第6,283,919号(Roundhill等)及び米国特許第6,458,083号(Jago等)に記載されている。ドップラープロセッサ28は、画像フィールドにおける血球の流れのような物質の動きを検出するため、組織の動き及び血流からの時間的に別個の信号を処理する。ドップラープロセッサは典型的に、身体内の選択されたタイプの物質から戻ったエコーを通過及び/又は阻止するように設定されたパラメータのウォールフィルタ(wall filter)を含む。例えばウォールフィルタは、高速の物質からの比較的低い振幅の信号を通過させると共に低速又はゼロ速度の物質からの比較的強い信号を阻止する通過帯域特性を有するように設定することができる。この通過帯域特性は、流れている血液からの信号を通過させると共に、心壁のようなほぼ静止しているか又はゆっくりと動く物体からの信号を阻止する。逆の特性では、動いている心臓組織からの信号を通過させると共に血流信号を阻止して、組織ドップラー撮像と称されるものを行い、組織の動きを検出して描写する。ドップラープロセッサは、画像フィールド内の様々なポイントから時間的に離散したエコー信号のシーケンスを受信及び処理する。特定ポイントからのエコーのシーケンスをアンサンブルと称する。比較的短い期間に高速で連続して受信されたエコーのアンサンブルを用いて、流れている血液のドップラーシフト周波数を推定することができ、ドップラー周波数と速度との対応が血流速度を示す。より長い時間期間に受信したエコーのアンサンブルを用いて、ゆっくり流れる血液又はゆっくり動く組織の速度を推定する。
Bモードプロセッサ及びドップラープロセッサにより生成された構造信号及び動き信号は、スキャン変換器32及び多面リフォーマッタ(multiplanar reformatter)44に結合される。スキャン変換器は、受信時の空間関係にあるエコー信号を所望の画像フォーマットに配列する。例えば、スキャン変換器はエコー信号を、2次元(2D)セクタ形状フォーマット又はピラミッド形状の3次元(3D)画像に配列する。スキャン変換器は、Bモード構造画像に、ドップラー推定速度を用いた画像フィールド内のポイントでの動きに対応した色を重ね合わせて、画像フィールド内の組織の動きと血流を表すカラードップラー画像を生成することができる。多面リフォーマッタは、身体の体積領域(volumetric region)の共通平面内のポイントから受信されたエコーを、その平面の超音波画像に変換する。これは米国特許第6,443,896号(Detmer)に記載されている。体積レンダラー(renderer)42は、3Dデータセットのエコー信号を、所与の基準点から見るような投影3D画像に変換する。これは米国特許第6,530,885号(Entrekin等)に記載されている。2D又は3D画像は、スキャン変換器32、多面リフォーマッタ44、及び体積レンダラー42から画像プロセッサ30に結合されて、更に画像を向上させたり、画像ディスプレイ40に表示するためのバッファリングや一時的な記憶を行ったりする。撮像に用いることに加えて、ドップラープロセッサ28により生成した血流値及びBモードプロセッサ26により生成した組織構造情報は、定量化プロセッサ34に結合される。定量化プロセッサは、血流の容積流量のような様々なフロー条件の測定値や臓器の大きさ及び妊娠期間のような構造測定値を生成する。定量化プロセッサは、ユーザ制御パネル38から、測定を行う画像の解剖学的組織内のポイント等の入力を受信する。定量化プロセッサからの出力データはグラフィックプロセッサ36に結合され、ディスプレイ40上の画像によって測定グラフィック及び測定値の再生が行われる。グラフィックプロセッサ36は、音響画像と共に表示するためのグラフィックオーバーレイも発生させることができる。これらのグラフィックオーバーレイは、患者の氏名、画像の日時、撮像パラメータ等の標準的な識別情報を含むことができる。これらの目的のため、グラフィックプロセッサは、患者の氏名等の入力をユーザインタフェース38から受信する。ユーザインタフェースは送信コントローラ18にも結合されて、トランスデューサアレイ10’からの音響信号の発生を制御し、従ってトランスデューサアレイ及び音響システムが生成する画像の発生を制御する。ユーザインタフェースは多面リフォーマッタ44にも結合され、多面リフォーマットされる(MPR)画像の面の選択及び制御が行われる。これを用いて、MPR画像の画像フィールドにおいて定量化測定を実行する。
音響撮像装置300は閉塞回避プロセッサ50も含む。これは、トランスデューサコントローラ18、ビームフォーマッタ20、信号プロセッサ22、及び画像プロセッサ30と協働して、音響トランスデューサアレイ10’の音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償すると共に、補償した画像データから撮像領域の画像を生成する。閉塞回避プロセッサ50によるこのような動作の例示的な実施形態については、以下で更に詳しく説明する。
図4は、撮像デバイスと撮像領域の少なくとも一部との間に1つ以上の障害物が存在する撮像領域を撮像する方法の一実施形態を示すフローチャートである。
動作410では、音響トランスデューサアレイを用いて、音響トランスデューサアレイと撮像領域の少なくとも一部との間に1つ以上の障害物が存在する撮像領域の画像データを生成する。
動作420では、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用して、1つ以上の障害物による撮像領域の画像データの欠落を補償する。
動作430では、補償した画像データから撮像エリアの画像を生成する。
冗長性は、フラウンホーファー近似が当てはまる遠距離において特に明らかである。その場合、フィールド及びアパーチャにおける振幅分布は、単純なフーリエ変換によって関連付けられる。例えば遠距離では、単一の要素によって発生するフィールドは平面波である。ITx及びIRxを、それぞれ送信要素と受信要素の指数とする。遠距離では、対応するフィールドはexp(j.k.x(ITx))(又はexp(j.k.x(IRx))である。ここで、kは波数であり、x(ITx))はアレイ内の送信要素の物理座標である。往復フィールドは、送信フィールドと受信フィールドの積であり、exp(j.k.(x(ITx)+x(IRx)))として表される。このため、x(ITx)+x(IRx)が同一である全ての要素対は同一の往復フィールドを発生する。これらの要素対は冗長であると言われる。
フィールドの角度スペクトルは、その横空間周波数への分解である。ここで、「角度スペクトル」という言葉は、PSFの横空間スペクトルに関連して用いられる。これは、PSFの横プロファイルのフーリエ変換として表される。遠距離では、先に検討したフラウンホーファー近似に従って、一方向(送信又は受信のみ)の角度スペクトルはアパーチャ振幅分布に比例している。アクティブなアパーチャが矩形である場合、一方向PSFはsinc関数である。往復点像分布関数はsinc関数であり、往復PSFの角度スペクトルは三角関数である。往復角度スペクトルは、送信アパーチャ及び受信アパーチャの振幅分布の畳み込み積として計算される。
図5は、(1)1つ以上の障害物がないと存在するはずの、撮像領域に対する200要素を有する矩形の音響トランスデューサアレイの理想的な角度スペクトル510(破線)と、(2)1つ以上の障害物がある撮像領域に対する同じ音響トランスデューサアレイの角度スペクトル520(実線)の例を示す。
図5からわかるように、障害物(例えば肋骨15)が肋骨間空間よりも大きくない限り、各周波数の振幅は非ゼロである。従って、障害物(例えば肋骨15)が存在するために小さく表される角度周波数に重みを加えることによって、真の角度スペクトル(及び点像分布関数)を回復させることができる。角度スペクトルにおける各周波数の重みは、障害物がない理想的な振幅(その周波数におけるスペクトル510の値)を、障害物が存在することによる実際の振幅(そのポイントにおけるスペクトル520の値)で除算した比である。音響トランスデューサアレイの全ての角度周波数におけるこれらの重み又は係数の集合を、ここで逆フィルタ(inverse filter)と称する。
この方法を実施するための1つの直接的な手法は、合成アパーチャ取得を用いることである。いったん障害物(複数の障害物)(例えば肋骨15)の位置(複数の位置)がわかったら、1つ以上の障害物(例えば肋骨15)が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトル520を計算することができ、角度スペクトル520と理想的な角度スペクトル510の比を算出することで各送信/受信要素対の重み又は係数が得られる。ここで、障害物(複数の障害物)の位置(複数の位置)又は場所(複数の場所)は、音響撮像装置300により生成された撮像エリア7の画像を用いる肋骨検出アルゴリズム及びフィーチャ認識アルゴリズムの使用を含む多種多様な技法によって決定される。2013年12月9日に出願された米国仮特許出願第61/913,576号は、肋骨検出アルゴリズムの実施形態を開示している。米国仮特許出願第61/913,576号の内容は引用により本願にも含まれるものとする。利用可能な他の方法は、X線データ、磁気共鳴映像法(MRI)データのような他の撮像様式から得られた画像データの使用を含む。しかしながら、障害物(複数の障害物)(例えば肋骨15)の位置(複数の位置)の特定にどの方法を用いるとしても、いったん場所(複数の場所)が明らかになったら、上述したように、1つ以上の障害物(例えば肋骨15)が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトル520を計算することができ、角度スペクトル520と理想的な角度スペクトル510の比を算出することで各送信/受信要素対の重み又は係数が得られる。合成アパーチャ加算の間、各送信/受信要素対を対応する重み係数で重み付けする。この方法の1つの欠点は、単一要素の送信を用いるために信号対雑音比(SNR)が低くなる可能性があることである。
図6は、撮像デバイスと撮像領域の少なくとも一部との間に1つ以上の障害物(例えば肋骨15)がある撮像領域を撮像する方法600の一実施形態を示すフローチャートである。方法600の様々な動作は、図6に示し以下に説明するものとは異なる順序で実行し得ることは理解されよう。
動作610では、障害物(複数の障害物)がないと存在するはずの、撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトルを求める。
動作620では、障害物(複数の障害物)が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトルを求める。
動作630では、音響トランスデューサアレイの逆フィルタを求める。
いくつかの実施形態では、逆フィルタを得るには、上述のように求められる1つ以上の障害物がないと存在するはずの撮像領域に対する音響トランスデューサアレイの理想的な角度スペクトル(例えば図5の510)を、1つ以上の障害物が存在する撮像領域に対する音響トランスデューサアレイの角度スペクトル(例えば図5の520)で除算すればよい。
動作640では、音響トランスデューサアレイを用いて撮像領域7の画像データを生成する。
動作650では、撮像動作において音響トランスデューサアレイの各送信/受信要素対によって生成した信号を、送信/受信対の角度周波数に対応する逆フィルタの値で重み付けする。
動作660では、補償した画像データから撮像エリアの画像を生成する。
上記の説明は、障害物の影響を簡略化し、それらが音響トランスデューサアレイのすぐ前に位置する薄い幕のようなフィーチャであるようにモデル化している。実際は、障害物(例えば肋骨15)は典型的に厚く、音響トランスデューサアレイから一定の距離に位置している。結果として、アパーチャ上の障害物の「影」は、撮像ポイントの位置と共に変化する。このため、有効な送信及び受信アパーチャ、並びに角度スペクトルは、撮像領域の各ポイントで変化する。しかしながら、障害物(複数の障害物)の場所(複数の場所)が明らかになったら、光線追跡アルゴリズムを用いて、各ポイントによって「見られる」アパーチャを計算し、各ポイントの角度スペクトルを計算し、上述したアルゴリズムを使用することができる。
従って、上述した方法600のいくつかの実施形態では、音響トランスデューサアレイの逆フィルタを求めること、音響トランスデューサアレイの各送信/受信要素対によって生成した信号を送信/受信対の角度周波数に対応する逆フィルタの値で重み付けすること、全ての送信/受信対の重み付けした信号を合計することで音響画像を生成することは、撮像領域内の対象となる複数ポイントの各々についてアルゴリズムを実行することを含み得る。このアルゴリズムは、ポイントに対する音響トランスデューサアレイの逆フィルタを求めることであって、逆フィルタを、1つ以上の障害物が存在するポイントに対する音響トランスデューサアレイの角度スペクトルで乗算すると、1つ以上の障害物がないと存在するはずのポイントに対する音響トランスデューサアレイの理想的な角度スペクトルが生成される、ことと、ポイントからの音響トランスデューサアレイの各送信/受信要素対によって生成した信号を、送信/受信対の角度周波数に対応する逆フィルタの値で重み付けすることと、ポイントからの音響トランスデューサアレイの各送信/受信要素対によって生成した重み付けした信号を合計することによってポイントにおける音響画像の強度を求めることと、を含み得る。
いくつかの実施形態では、各ポイントに対する音響トランスデューサアレイの逆フィルタは、光線追跡アルゴリズムを用いて、ポイントによって見られる有効アパーチャを計算することにより、求めることができる。
図7は肋骨モデルを示し、障害物15が存在する場合に、図6の方法600において、どのように光線追跡アルゴリズムを用いてどの画像ポイントがトランスデューサアレイ内のどのトランスデューサ要素を見るかを検出するかを説明する。具体的には、図7は、特定の要素702が見ることができる第1のエリア722と、要素702からの視野が障害物15によって遮られる第2のエリア724と、を示す。
柔軟性を保ちながら、方法600で用いられる合成アパーチャ撮像の低いSNRを軽減するための1つの解決策は、仮想トランスデューサ手法を用いることである。このような手法では、1組の集束した送信及び受信対を生成する。各集束ビームは、ビーム焦点に位置する仮想トランスデューサ(しかしながら単一要素よりも強度が高い)として見ることができ、合成アパーチャアルゴリズムは、この仮想トランスデューサアレイを用いて適用可能である。仮想トランスデューサの理想的な深さは障害物(例えば肋骨15)の深さである。
従って、上述した方法600のいくつかの実施形態では、音響トランスデューサアレイの逆フィルタを求めること、音響トランスデューサアレイの各送信/受信要素対によって生成した信号を送信/受信対の角度周波数に対応する逆フィルタの値で重み付けすること、及び全ての送信/受信対の重み付けした信号を合計することで音響画像を生成することは、仮想トランスデューサを用いるアルゴリズムを実行することを含み得る。いくつかの実施形態では、このアルゴリズムは、音響トランスデューサアレイから仮想トランスデューサアレイを生成することと、仮想トランスデューサアレイのための逆フィルタを求めることであって、逆フィルタを、1つ以上の障害物が存在する撮像領域に対する仮想トランスデューサアレイの角度スペクトルで乗算すると、1つ以上の障害物がないと存在するはずの撮像領域に対する仮想トランスデューサアレイの理想的な角度スペクトルが生成される、逆フィルタを求めることと、仮想トランスデューサアレイの各送信/受信要素対によって生成した信号を、送信/受信対の角度周波数に対応する逆フィルタの値で重み付けすることと、全ての送信受信対の重み付けした信号を合計することによって音響画像を生成することと、を含む。
他の実施形態は、音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を、ビーム空間において直接に利用することができる。そのような実施形態は従来のビーム形成を用い、空間内のいくつかのポイントからのエネルギを1つの位置に集束できるので、図6の実施形態に比べてSNRを上昇させ得る。これらの実施形態では、単一要素の送信の代わりに集束送信を用いる。この目的のため、計算される往復重みを、送信アパーチャのための重みと受信アパーチャのための重みとに分割しなければならない。これは、特異値分解、固有値分解、又は他の分解のような分解アルゴリズムを用いて実行可能である。しかしながら、所望の往復点像分布関数は単一の送信−受信動作で達成することができず、2つ以上の送信−受信動作の合計を必要とする(例えば、第2の送信/受信点像分布関数は実際はサイドローブのみを含み、メインローブとサイドローブの双方を含む第1のものから減じることができる)。
例えば、往復重みの行列をWとする。例えばWi,jは、対応する送信/受信対の角度スペクトルの値により与えられる送信iと受信jから成る対に対応した重みである。
U(又はV)を、送信(又は受信)におけるアポディゼーションベクトルとする(Uの各要素は送信の場合にアレイ内の要素に適用される重みである)。所与のUとVについて、WはW=UV*によって与えられる。ここで*は行列転置を表す。
ここで、興味深いのは逆の問題である。すなわち、Wは既知であるが、対応する1組の送信及び受信アパーチャを見つけることは常に可能かということである。Wのランクが1である場合、Wは1送信アパーチャ及び1受信アパーチャの積に分解することができる。しかしながらWのランクが1より大きいと、ギャップ補償アルゴリズムの場合のように、Wは、特異値分解(SVD)によって、いくつかの送信及び受信アパーチャの積の和に分解できるだけである。
W=USV* ここで、Sは特異値の体格行列であり、U(又は)は、列が様々な送信(又は受信)アポディゼーションベクトルである行列である。
実際には、2つのみの特異値が有意である(他のものは極めて小さい)ので、Wは、2つの送信−受信ビームの和として表される極めて良好な近似であり得ることがわかっている。
W=S(1,1).U1.+S(2,2).U2. (1)
Wの対称性のため、送信及び受信アポディゼーションは同一である。すなわちU=Vである。
図8は、図2に関して上述したシミュレーションしたアパーチャのアポディゼーション810(U1)及び820(U2)を示す。双方のアポディゼーション810及び820は、サブアパーチャの縁部の近傍に大きい重みを有する。アポディゼーション810のみが大きいサイドローブを生成する。しかしながら、アポディゼーション810及び820の各々を用いて生成した点像分布関数が上記の式(1)に従って組み合わされた場合、結果として得られる点像分布関数からグレーティングローブは消える。実際、アポディゼーション820が生成する点像分布関数では、メインローブは存在せずグレーティグローブのみが存在し、このグレーティングローブは、アポディゼーション810を用いて生成した第1の点像分布関数から減じられる。
いくつかの実施形態では、第1のアポディゼーション810(U)を送信及び受信動作の双方で用いて第1の画像データセットを生成することができる。同様に、送信及び受信動作で第2のアポディゼーション820(U)を用いて第2の画像データセットを発生させる。次いで、2つの合計データセットを、式1に示すようなそれぞれの特異値で重み付けし、これらを組み合わせることで、最終的な画像データを得る。
図9は、撮像デバイスと撮像領域の少なくとも一部との間に1つ以上の障害物(例えば肋骨15)がある撮像領域を撮像する方法900の一実施形態を示すフローチャートである。方法900の様々な動作は、図9に示し以下に説明するものとは異なる順序で実行し得ることは理解されよう。
動作910では、撮像領域に対する音響トランスデューサアレイの逆フィルタから、第1及び第2のアポディゼーション関数を求める。
図8は、図9の方法900で用いられる2組の送信/受信アパーチャの第1及び第2のアポディゼーション関数810及び820の例を示す。
動作920では、第1のアポディゼーション関数を音響トランスデューサアレイに適用する。
動作930では、第1のアポディゼーションを用いて第1の送信−受信動作を実行して、第1の画像データを生成する。
動作940では、第2のアポディゼーション関数を音響トランスデューサアレイに適用する。
動作950では、第2のアポディゼーションを用いて第2の送信−受信動作を実行して、第2の画像データを生成する。
動作960では、第1の画像データと第2の画像データとを組み合わせることで、撮像領域の画像を生成する。
音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用する方法の他の実施形態は、画像に基づくデコンボリューションを用いる。上述のように、遠距離では、所与の周波数における点像分布関数のフーリエ変換はアパーチャに比例する。広帯域信号では、k空間として知られる点像分布関数の2D空間フーリエ変換を求めることができる。所与の波数kにおけるスペクトルは、上述のように導出される往復角度スペクトルの倍率を変えたものである。倍率は波数k自体である。これはフラウンホーファー近似から導出することができる。上述のことからわかるように、アパーチャの1つの要素は平面波exp(j.k.x(ITx))に関連付けられる。空間周波数はkに比例する。
従って、逆フィルタは、合成アパーチャ往復対でなくk空間において直接実施することができる。
いくつかの実施形態では、画像のk空間は、合計RFデータの2Dフーリエ変換を求めることによって計算される(包絡線検出の前)。逆フィルタの重みを上述のように計算し、kの各値について倍率を変え、スペクトルを乗算する。この時点で、信号の時間的な帯域幅を考慮に入れることができる。最後に、2D逆フーリエ変換を計算して、補償された画像を得る。あるいは、k空間における乗算が画像ドメインにおけるコンボリューションと同等である場合、デコンボリューションを画像ドメインで直接実行することができる。その場合、逆フィルタの2Dフーリエ変換を計算して、画像でコンボリューションされるカーネルを得る(合計RFデータ)。カーネルのサイズは、概ね点像分布関数のサイズである。
画像に基づくデコンボリューション使用の基礎にある想定は、上述したような合成アパーチャアルゴリズム使用の基礎にある想定ほどロバストでないことが観察されている。遠距離においても、k空間スペクトルは、ポイントの横位置の関数として大きく変動する。しかしながらこれは、空間依存逆フィルタを用いることで対処可能である。フィルタの重みは、点像分布関数のシミュレーションによって計算することができる。コンボリューションカーネル手法は、カーネルが位置の関数として変化し得るので、空間依存フィルタリングに特に適している。
図10は、図6に関して説明したような補償アルゴリズムを肋骨間空間に用いた場合の、図2に示したものと同じ状況での点像分布関数を示す。ここでも、4ポイント散乱体が深さ4cm、8cm、12cm、16cmに位置している。12cmと16cmにある深い方の2つのポイントでは、画像は図2に比べて大幅に改善されている。グレーティングローブは抑制され、点像分布関数は、ギャップのない完全なアパーチャで得られるものと同様に見える。最も浅い4cmのポイントでは、遠距離近似は成立せず、再構築は良好でない。8cmのポイントは遷移領域である。
図11A及び図11Bは、大きいアパーチャにおいて肋骨間ギャップの存在をシミュレーションした画像を、補償あり(図11B)及び補償なし(図11A)の場合で示す。具体的には、1つの80要素アパーチャを用い、それらの要素の一部を遮ることで「小型肋骨」を人工的に導入した。2つの13要素妨害物を用い、3つの18要素アパーチャを残した。図11Aは顕著なグレーティングローブを示し、図11Bは、本明細書に記載したギャップ補償アルゴリズムを用いることによってグレーティングローブが低減したこと、場合によっては除去されたことを示す。
図12Aから図12Dは、2つの妨害物(例えば2本の肋骨15による)をシミュレーションした場合の、大きいアパーチャを用いて生成された経胸壁心エコー図のシミュレーションを示す。図11Aから図11Bにおいてと同様、1つの80要素アパーチャが、2つの13要素妨害物で分離された3つの18要素アパーチャに分割されている。
図12Aは、2つの妨害物(例えば2本の肋骨による)をシミュレーションした場合の、大きいアパーチャを用いて生成された経胸壁心エコー図のシミュレーションを示す。
図12Bは、2つの妨害物(例えば2本の肋骨による)をシミュレーションし、肋骨間ギャップ補償アルゴリズムを用いた場合の、大きいアパーチャを用いて生成された経胸壁心エコー図を示す。
図12Cは、図12Aと同様であるが妨害物が存在しない、大きいアパーチャを用いて生成された基準経胸壁心エコー図を示す。
図12Dは、図8Aから図8Cを生成するために用いたものよりも小さいアパーチャを用いて生成された経胸壁心エコー図を示す。
図13は、障害物(例えば15)が存在する生成されたコヒーレント画像を示す。図2及び図10においてと同様、4ポイント散乱体が深さ4cm、8cm、12cm、16cmに位置している。グレーティングローブ1304はメインローブ1302に隣接して示されている。
図14Aは、上述したような「肋骨間ギャップ補償」アルゴリズムを用いるコヒーレントな複合を利用した場合の、図13に示すものと同様の肋骨が存在する生成された画像を示す。図14Bは、障害物が存在しないと生成されるはずの、図13及び図14Aに示すものと同様の画像を示す。図13、図14A、及び図14Bを比較することによって、10cmを超える深さでは、補償アルゴリズムを用いた場合の障害物(例えば肋骨)を介した画像は、障害物が存在しない場合とほぼ同じくらい良好であることがわかる。
図15Aは、図8に示したモデルにおける遠距離中心のポイントの角度スペクトルを示す。図15Bは、図8に示したモデルにおける画像の縁部上の、音響トランスデューサアレイに近いポイントの角度スペクトルを示す。図15Bに示すように、角度スペクトルは、送信/受信往復経路における冗長性を利用した上述のアルゴリズムによって埋めることができないギャップ1505を有する。これは、最小振幅がゼロであるか又は極めて小さいからである。矢印で示すスペクトルの部分1510に、コヒーレント複合アルゴリズムを用いることができる。スペクトルの他の部分はインコヒーレントに追加される。
図16は、図8の例について、画像の様々なポイントについてこのアルゴリズムにより再構築され得る最大往復アパーチャサイズを示す。障害物は左の縁部に沿って位置している。
図17は、障害物(例えば肋骨)が肋骨間空間と同じくらい大きい場合の補償アルゴリズムの結果を示す。ここでは、ギャップを埋めることができないので、各ポイントに1つだけ肋骨間空間が用いられる。インコヒーレントな複合を用いてこの結果を改善することができる。
本明細書において好適な実施形態を詳細に開示したが、本発明の概念及び範囲内に収まる多くの変形が可能である。そのような変形は、本明細書、図面、及び特許請求の範囲を精査した後に、当業者には明らかとなろう。従って本発明は、添付の特許請求の範囲内にある場合を除いて制約を受けないものとする。

Claims (15)

  1. 撮像領域を撮像する方法であって、
    音響トランスデューサアレイを用いて、前記音響トランスデューサアレイと前記撮像領域の少なくとも一部との間に1つ以上の障害物がある前記撮像領域の画像データを生成することと、
    (i)前記音響トランスデューサアレイのための逆フィルタを求めることであって、前記1つ以上の障害物がないと存在するはずの前記撮像領域に対する前記音響トランスデューサアレイの理想的な角度スペクトルを、前記1つ以上の障害物が存在する前記撮像領域に対する前記音響トランスデューサアレイの角度スペクトルで除算することで、前記逆フィルタを求めることと、
    (ii)前記音響トランスデューサアレイの各送信/受信要素対によって生成した信号を、前記送信/受信対の角度周波数に対応する前記逆フィルタの値で重み付けすることと、
    によって、前記1つ以上の障害物による前記撮像領域の画像データの欠落を補償するために、前記音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用することと、
    全ての前記送信/受信対の前記重み付けした信号を合計することによって、補償した前記画像データから前記撮像領域の画像を生成することと、
    を有する、方法。
  2. 音響トランスデューサアレイを用いて前記撮像領域の画像データを生成することが、合成アパーチャ取得を実行することを有する、請求項1に記載の方法。
  3. 前記音響トランスデューサアレイの前記逆フィルタを求めること、前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記信号を前記送信/受信対の前記角度周波数に対応する前記逆フィルタの前記値で重み付けすること、及び全ての前記送信/受信対の前記重み付けした信号を合計することで音響画像を生成することが、
    前記撮像領域内の対象となる複数ポイントの各々について、
    前記ポイントに対する前記音響トランスデューサアレイの逆フィルタを求めることであって、前記逆フィルタを、前記1つ以上の障害物が存在する前記ポイントに対する前記音響トランスデューサアレイの角度スペクトルで乗算すると、前記1つ以上の障害物がないと存在するはずの前記ポイントに対する前記音響トランスデューサアレイの理想的な角度スペクトルが生成される、前記逆フィルタを求めることと、
    前記ポイントからの前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記信号を、前記送信/受信対の角度周波数に対応する前記逆フィルタの値で重み付けすることと、
    前記ポイントからの前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記重み付けした信号を合計することによって前記ポイントにおける前記音響画像の強度を求めることと、
    を有する、請求項1に記載の方法。
  4. 前記ポイントに対する前記音響トランスデューサアレイの前記逆フィルタを求めることが、光線追跡アルゴリズムを用いて前記ポイントによって見られる有効アパーチャを計算することを含む、請求項3に記載の方法。
  5. 前記1つ以上の障害物による前記撮像領域の画像データの欠落を補償するために、前記音響トランスデューサアレイの前記音響トランスデューサ間の送信/受信対経路における冗長性を利用すること、及び前記補償した画像データから前記撮像領域の画像を生成することが、
    少なくとも第1及び第2のアポディゼーション関数を用いて少なくとも第1及び第2の送信及び受信動作を実行することであって、
    前記音響トランスデューサアレイを用いて、第1の音響波を前記撮像領域に送信し、前記撮像領域から第1の音響エコーを受信し、そこから第1の画像データを生成することによって前記第1の送信及び受信動作を実行することであって、前記第1のアポディゼーション関数を前記音響トランスデューサアレイに適用して前記第1の送信及び受信動作のための第1の送信アパーチャ及び第1の受信アパーチャを生成する、前記第1の送信及び受信動作を実行することと、
    前記音響トランスデューサアレイを用いて、第2の音響波を前記撮像領域に送信し、前記撮像領域から第2の音響エコーを受信し、そこから第2の画像データを生成することによって前記第2の送信及び受信動作を実行することであって、前記第2のアポディゼーション関数を前記音響トランスデューサアレイに適用して前記第2の送信及び受信動作のための第2の送信アパーチャ及び第2の受信アパーチャを生成する、前記第2の送信及び受信動作を実行することと、
    を含む、第1及び第2の送信及び受信動作を実行することと、
    前記第1の画像データを前記第2の画像データと組み合わせることによって前記撮像領域の画像を生成することと、
    を有し、前記少なくとも第1及び第2のアポディゼーション関数が前記撮像領域に対する前記音響トランスデューサアレイの逆フィルタから求められる、請求項1に記載の方法。
  6. 前記少なくとも第1及び第2のアポディゼーション関数が、分解アルゴリズムを用いることによって、前記撮像領域に対する前記音響トランスデューサアレイの前記逆フィルタから求められる、請求項に記載の方法。
  7. 前記音響トランスデューサアレイの前記音響トランスデューサ間の送信/受信対経路における冗長性を利用して、前記1つ以上の障害物による前記撮像領域の画像データの欠落を補償すること、及び前記補償した画像データから前記撮像領域の画像を生成することが、RFデータに基づく又は画像に基づくデコンボリューションアルゴリズムを実行することを有する、請求項1に記載の方法。
  8. 撮像領域を撮像するための装置であって、前記装置が、
    音響トランスデューサアレイと前記撮像領域の少なくとも一部との間に1つ以上の障害物がある前記撮像領域の画像データを生成する音響トランスデューサアレイと、
    (i)前記音響トランスデューサアレイのための逆フィルタを求めることであって、前記1つ以上の障害物がないと存在するはずの前記撮像領域に対する前記音響トランスデューサアレイの理想的な角度スペクトルを、前記1つ以上の障害物が存在する前記撮像領域に対する前記音響トランスデューサアレイの角度スペクトルで除算することで、前記逆フィルタを求めることと、
    (ii)前記音響トランスデューサアレイの各送信/受信要素対によって生成した信号を、前記送信/受信対の角度周波数に対応する前記逆フィルタの値で重み付けすることと、
    によって、前記1つ以上の障害物による前記撮像領域の画像データの欠落を補償するために、前記音響トランスデューサアレイの音響トランスデューサ間の送信/受信対経路における冗長性を利用する1つ以上のプロセッサと、
    を備え、前記1つ以上のプロセッサが更に、
    全ての前記送信/受信対の前記重み付けした信号を合計することによって、前記補償した画像データから前記撮像領域の画像を生成する、装置。
  9. 前記1つ以上のプロセッサが更に、
    前記音響トランスデューサアレイに対する前記1つ以上の障害物の1つ以上の位置を求め、
    前記求めた1つ以上の位置を用いて、前記1つ以上の障害物が存在する前記撮像領域に対する前記音響トランスデューサアレイの前記角度スペクトルを求める、請求項8に記載の装置。
  10. 前記装置が合成アパーチャ取得を実行することによって前記画像データを生成する、請求項9に記載の装置。
  11. 前記1つ以上のプロセッサが、前記音響トランスデューサアレイの前記逆フィルタを求め、前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記信号を前記送信/受信対の前記角度周波数に対応する前記逆フィルタの前記値で重み付けし、全ての前記送信/受信対の前記重み付けした信号を合計することで音響画像を生成することが、
    前記撮像領域内の対象となる複数ポイントの各々について、
    前記ポイントに対する前記音響トランスデューサアレイの逆フィルタを求めることであって、前記逆フィルタを、前記1つ以上の障害物が存在する前記ポイントに対する前記音響トランスデューサアレイの角度スペクトルで乗算すると、前記1つ以上の障害物がないと存在するはずの前記ポイントに対する前記音響トランスデューサアレイの理想的な角度スペクトルが生成される、前記逆フィルタを求めることと、
    前記ポイントからの前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記信号を、前記送信/受信対の角度周波数に対応する前記逆フィルタの値で重み付けすることと、
    前記ポイントからの前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記重み付けした信号を合計することによって前記ポイントにおける前記音響画像の強度を求めることと、
    によって行われ、
    前記1つ以上のプロセッサが、光線追跡アルゴリズムを用いて前記ポイントによって見られる有効アパーチャを計算することによって、前記ポイントに対する前記音響トランスデューサアレイの前記逆フィルタを求める、請求項8に記載の装置。
  12. 前記1つ以上のプロセッサが、前記音響トランスデューサアレイの前記逆フィルタを求めること、前記音響トランスデューサアレイの各送信/受信要素対によって生成した前記信号を前記送信/受信対の前記角度周波数に対応する前記逆フィルタの前記値で重み付けすること、及び全ての前記送信/受信対の前記重み付けした信号を合計することで音響画像を生成することが、
    前記音響トランスデューサアレイから仮想トランスデューサアレイを生成することと、
    前記仮想トランスデューサアレイのための逆フィルタを求めることであって、前記逆フィルタを、前記1つ以上の障害物が存在する前記撮像領域に対する前記仮想トランスデューサアレイの角度スペクトルで乗算すると、前記1つ以上の障害物がないと存在するはずの前記撮像領域に対する前記仮想トランスデューサアレイの理想的な角度スペクトルが生成される、前記逆フィルタを求めることと、
    コーから前記仮想トランスデューサアレイの各送信/受信要素対によって生成した信号を、前記送信/受信対の角度周波数に対応する前記逆フィルタの値で重み付けすることと、
    前記エコーからの全ての前記送信受信対の前記重み付けした信号を合計することによって音響画像を生成することと、
    によって行われる、請求項8に記載の装置。
  13. 前記1つ以上のプロセッサが、前記音響トランスデューサアレイの前記音響トランスデューサ間の送信/受信対経路における冗長性を利用して、前記1つ以上の障害物による前記撮像領域の画像データの欠落を補償すること、及び前記補償した画像データから前記撮像領域の画像を生成することが、
    少なくとも第1及び第2のアポディゼーション関数を用いて少なくとも第1及び第2の送信及び受信動作を実行することであって、
    前記音響トランスデューサアレイを用いて、第1の音響波を前記撮像領域に送信し、前記撮像領域から第1の音響エコーを受信し、そこから第1の画像データを生成することによって前記第1の送信及び受信動作を実行することであって、前記第1のアポディゼーション関数を前記音響トランスデューサアレイに適用して前記第1の送信及び受信動作のための第1の送信アパーチャ及び第1の受信アパーチャを生成する、前記第1の送信及び受信動作を実行することと、
    前記音響トランスデューサアレイを用いて、第2の音響波を前記撮像領域に送信し、前記撮像領域から第2の音響エコーを受信し、そこから第2の画像データを生成することによって前記第2の送信及び受信動作を実行することであって、前記第2のアポディゼーション関数を前記音響トランスデューサアレイに適用して前記第2の送信及び受信動作のための第2の送信アパーチャ及び第2の受信アパーチャを生成する、前記第2の送信及び受信動作を実行することと、
    を含む、第1及び第2の送信及び受信動作を実行することと、
    前記第1の画像データを前記第2の画像データと組み合わせることによって前記撮像領域の画像を生成することと、
    によって行われ、前記少なくとも第1及び第2のアポディゼーション関数が前記撮像領域に対する前記音響トランスデューサアレイの前記逆フィルタから求められる、請求項8に記載の装置。
  14. 前記1つ以上のプロセッサが、分解アルゴリズムを用いることによって、前記撮像領域に対する前記音響トランスデューサアレイの前記逆フィルタから前記少なくとも第1及び第2のアポディゼーション関数を求める、請求項13に記載の装置。
  15. 前記1つ以上のプロセッサが、前記音響トランスデューサアレイの前記音響トランスデューサ間の送信/受信対経路における冗長性を利用して、前記1つ以上の障害物による前記撮像領域の画像データの欠落を補償すること、及び前記補償した画像データから前記撮像領域の画像を生成することが、RFデータに基づく又は画像に基づくデコンボリューションアルゴリズムを実行することによって行われる、請求項8に記載の装置。
JP2016559411A 2014-03-31 2015-03-18 肋骨間空間を用いたコヒーレントな複合による音響撮像のためのシステム及び方法 Expired - Fee Related JP6596013B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201461972896P 2014-03-31 2014-03-31
US61/972,896 2014-03-31
PCT/IB2015/051972 WO2015150954A1 (en) 2014-03-31 2015-03-18 System and method for acoustic imaging with coherent compounding using intercostal spaces

Publications (3)

Publication Number Publication Date
JP2017509416A JP2017509416A (ja) 2017-04-06
JP2017509416A5 JP2017509416A5 (ja) 2018-04-26
JP6596013B2 true JP6596013B2 (ja) 2019-10-23

Family

ID=52814157

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016559411A Expired - Fee Related JP6596013B2 (ja) 2014-03-31 2015-03-18 肋骨間空間を用いたコヒーレントな複合による音響撮像のためのシステム及び方法

Country Status (6)

Country Link
US (1) US10746859B2 (ja)
EP (1) EP3126872B1 (ja)
JP (1) JP6596013B2 (ja)
CN (1) CN106132308B (ja)
RU (1) RU2676560C2 (ja)
WO (1) WO2015150954A1 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11125866B2 (en) * 2015-06-04 2021-09-21 Chikayoshi Sumi Measurement and imaging instruments and beamforming method
CN109313604B (zh) * 2016-06-27 2023-05-23 英特尔公司 用于压缩虚拟存储器的动态配置的计算系统、装置和方法
CN106680825B (zh) * 2016-12-05 2019-06-04 中国科学院声学研究所 一种声学阵列成像系统与方法
EP3382423A1 (en) * 2017-03-27 2018-10-03 Koninklijke Philips N.V. Methods and systems for filtering ultrasound image clutter
US11147538B2 (en) * 2018-01-15 2021-10-19 International Business Machines Corporation Efficient computation of spatially varying ultrasound analytical point spread functions
JP7401459B2 (ja) * 2018-01-23 2023-12-19 コーニンクレッカ フィリップス エヌ ヴェ 針挿入誘導を提供する超音波イメージングシステム
CN112119327A (zh) * 2018-05-15 2020-12-22 皇家飞利浦有限公司 具有声速像差映射的合成透射聚焦超声系统
WO2019232451A1 (en) * 2018-05-31 2019-12-05 Matt Mcgrath Design & Co, Llc Method of medical imaging using multiple arrays
CN110068833B (zh) * 2019-05-05 2021-10-29 中国科学院电子学研究所 一种合成孔径激光雷达成像方法、仪器及系统
JP2023540566A (ja) * 2020-09-08 2023-09-25 ディープサイト テクノロジー インコーポレイテッド 混合アレイを使用する合成開口撮像システムおよび方法
CN113081051A (zh) * 2021-03-27 2021-07-09 四川维思模医疗科技有限公司 一种tee检查装置、检查系统及成像方法
CN113281729B (zh) * 2021-05-31 2022-02-01 中国科学院声学研究所 一种基于多帧空间谱联合处理的目标自动检测方法及系统
WO2022260112A1 (ja) * 2021-06-11 2022-12-15 株式会社 Integral Geometry Science 映像化装置及び映像化方法
US20240295506A1 (en) * 2021-06-17 2024-09-05 Kenjiro Kimura Imaging device and imaging method
CN113391311B (zh) * 2021-06-21 2022-08-09 电子科技大学 一种分布式雷达广义孔径合成方法

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5278757A (en) * 1991-11-15 1994-01-11 The Trustees Of The University Of Pennsylvania Synthetic aperture ultrasonic imaging system using a minimum or reduced redundancy phased array
JP3059042B2 (ja) * 1994-02-22 2000-07-04 フクダ電子株式会社 超音波診断装置
US6283919B1 (en) 1996-11-26 2001-09-04 Atl Ultrasound Ultrasonic diagnostic imaging with blended tissue harmonic signals
US6458083B1 (en) 1996-11-26 2002-10-01 Koninklijke Philips Electronics N.V. Ultrasonic harmonic imaging with adaptive image formation
US6013032A (en) 1998-03-13 2000-01-11 Hewlett-Packard Company Beamforming methods and apparatus for three-dimensional ultrasound imaging using two-dimensional transducer array
US5997479A (en) 1998-05-28 1999-12-07 Hewlett-Packard Company Phased array acoustic systems with intra-group processors
US6511426B1 (en) * 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6530885B1 (en) 2000-03-17 2003-03-11 Atl Ultrasound, Inc. Spatially compounded three dimensional ultrasonic images
US6443896B1 (en) 2000-08-17 2002-09-03 Koninklijke Philips Electronics N.V. Method for creating multiplanar ultrasonic images of a three dimensional object
US6468216B1 (en) 2000-08-24 2002-10-22 Kininklijke Philips Electronics N.V. Ultrasonic diagnostic imaging of the coronary arteries
US7611465B2 (en) * 2003-07-15 2009-11-03 Board Of Regents, The University Of Texas System Rapid and accurate detection of bone quality using ultrasound critical angle reflectometry
US7330594B2 (en) * 2003-07-31 2008-02-12 Kabushiki Kaisha Toshiba Image enhancement or correction software, method, apparatus and system for substantially minimizing blur in the scanned image
US20050124883A1 (en) * 2003-11-20 2005-06-09 Hunt Thomas J. Adaptive parallel artifact mitigation
US9310475B2 (en) * 2003-11-21 2016-04-12 General Electric Company Method and apparatus for transmitting multiple beams
US8234923B2 (en) * 2004-09-20 2012-08-07 Innervision Medical Technologies Inc. Systems and methods for ultrasound imaging
US8133180B2 (en) * 2004-10-06 2012-03-13 Guided Therapy Systems, L.L.C. Method and system for treating cellulite
EP2088932B1 (en) * 2006-10-25 2020-04-08 Maui Imaging, Inc. Method and apparatus to produce ultrasonic images using multiple apertures
JP5069022B2 (ja) * 2007-03-06 2012-11-07 ゼネラル・エレクトリック・カンパニイ 超音波撮像で使用するための正確な時間遅延推定の方法及びシステム
EP3563768A3 (en) * 2010-10-13 2020-02-12 Maui Imaging, Inc. Concave ultrasound transducers and 3d arrays
US20120302927A1 (en) * 2011-05-23 2012-11-29 University Of Washington Methods for characterizing nonlinear fields of a high-intensity focused ultrasound source and associated systems and devices
EP2833791B1 (en) * 2012-03-26 2022-12-21 Maui Imaging, Inc. Methods for improving ultrasound image quality by applying weighting factors
JP2014023934A (ja) * 2012-07-27 2014-02-06 Samsung Electronics Co Ltd 画像処理モジュール及び画像生成方法
CN103584886B (zh) * 2013-11-20 2015-07-15 无锡祥生医学影像有限责任公司 一种基于相位相干信息的自适应变迹方法

Also Published As

Publication number Publication date
EP3126872B1 (en) 2021-05-12
EP3126872A1 (en) 2017-02-08
JP2017509416A (ja) 2017-04-06
WO2015150954A1 (en) 2015-10-08
RU2016142405A3 (ja) 2018-11-01
US20170184713A1 (en) 2017-06-29
CN106132308A (zh) 2016-11-16
RU2676560C2 (ru) 2019-01-09
RU2016142405A (ru) 2018-05-07
CN106132308B (zh) 2020-03-03
US10746859B2 (en) 2020-08-18

Similar Documents

Publication Publication Date Title
JP6596013B2 (ja) 肋骨間空間を用いたコヒーレントな複合による音響撮像のためのシステム及び方法
CN110023782B (zh) 用于对超声图像杂波滤波的方法和系统
Vignon et al. Capon beamforming in medical ultrasound imaging with focused beams
US20190083058A1 (en) Method and apparatus to produce ultrasonic images using multiple apertures
Szasz et al. Beamforming through regularized inverse problems in ultrasound medical imaging
Bendjador et al. The SVD beamformer: Physical principles and application to ultrafast adaptive ultrasound
US8818064B2 (en) Time-domain estimator for image reconstruction
US10064602B2 (en) Coherence ultrasound imaging with broad transmit beams
CN112513675B (zh) 超声方法和装置
Zhang et al. Ultrafast ultrasound imaging using combined transmissions with cross-coherence-based reconstruction
US20190090850A1 (en) High-resolution and/or high-contrast 3-d and/or 4-d ultrasound imaging with a 1-d transducer array
WO2018099867A1 (en) Methods and systems for filtering ultrasound image clutter
JP6932200B2 (ja) 超音波画像クラッターをフィルタ処理するための方法及びシステム
CN110891492B (zh) 用于处理超声图像的方法和系统
JP6998477B2 (ja) カラードップラー超音波イメージングを行うための方法及びシステム
Morgan et al. Synthetic aperture focusing for multi-covariate imaging of sub-resolution targets
Jakovljevic et al. Short-lag spatial coherence imaging on matrix arrays, Part II: Phantom and in vivo experiments
CN110869799B (zh) 用于处理超声图像的方法和系统
Peralta et al. Extension of coherent multi-transducer ultrasound imaging with diverging waves
Ahmed et al. Large-array deep abdominal imaging in fundamental and harmonic mode
Guo et al. Wavenumber Beamforming with Sub-Nyquist Sampling for Focus-Beam Ultrasound Imaging
JP7345678B2 (ja) 3dベクトルフロー場を取得するための方法及びシステム
Gong et al. Simulation studies of filtered spatial compounding (FSC) and filtered frequency compounding (FFC) in synthetic transmit aperture (STA) imaging
Vrålstad et al. Universal REFoCUS beamforming with spatial weighting
Dozier-Muhammad Improved image uniformity using optimized diverging-wave acquisition sequence for high frame rate pulse-echo ultrasound

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180315

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180315

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20181130

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20181218

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20190314

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190610

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190927

R150 Certificate of patent or registration of utility model

Ref document number: 6596013

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees