JP4767377B2 - 物体ボリュームの3次元イメージングのためのシステム及び方法 - Google Patents
物体ボリュームの3次元イメージングのためのシステム及び方法 Download PDFInfo
- Publication number
- JP4767377B2 JP4767377B2 JP11569099A JP11569099A JP4767377B2 JP 4767377 B2 JP4767377 B2 JP 4767377B2 JP 11569099 A JP11569099 A JP 11569099A JP 11569099 A JP11569099 A JP 11569099A JP 4767377 B2 JP4767377 B2 JP 4767377B2
- Authority
- JP
- Japan
- Prior art keywords
- low
- data
- harmonic mode
- projection
- harmonic
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims description 20
- 238000000034 method Methods 0.000 title claims description 19
- 238000002604 ultrasonography Methods 0.000 claims description 14
- 230000006870 function Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000002592 echocardiography Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 4
- 239000002872 contrast media Substances 0.000 description 22
- 239000000523 sample Substances 0.000 description 21
- 210000001519 tissue Anatomy 0.000 description 12
- 238000010304 firing Methods 0.000 description 9
- 238000012285 ultrasound imaging Methods 0.000 description 8
- 210000003484 anatomy Anatomy 0.000 description 6
- 230000017531 blood circulation Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 230000011218 segmentation Effects 0.000 description 6
- 210000004369 blood Anatomy 0.000 description 5
- 239000008280 blood Substances 0.000 description 5
- 238000005266 casting Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 210000004027 cell Anatomy 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000007792 addition Methods 0.000 description 3
- 210000000601 blood cell Anatomy 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 210000000577 adipose tissue Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008081 blood perfusion Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8993—Three dimensional imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/916—Ultrasound 3-D imaging
Description
【発明の分野】
本発明は、一般的には、医療診断の目的のための人体の解剖学的構造の超音波イメージングに関する。具体的には、本発明は、人体内の運動する流体又は組織に超音波を送信して、これらの流体又は組織から反射された超音波エコーを検出することによって、人体内の運動する流体又は組織を3次元イメージングするための方法及び装置に関する。
【0002】
【発明の背景】
従来の超音波スキャナは、Bモード及びカラー・フロー・モード等の様々なイメージング・モードで動作することが可能である。Bモードでは、表示ピクセルの輝度がエコー反射の強度に基づいているような2次元画像が形成される。カラー・フロー・イメージングでは、血液の流れ又は組織の運動をイメージングすることができる。イメージングされている物体に動きが存在していれば、反射信号に、動きの速度に比例する(送信信号に対する)ドプラ・シフトが生ずる。
【0003】
ドプラ効果に基づいて血流を検出する超音波スキャナは周知である。これらのシステムは、物体内に超音波を送信し且つ物体から後方散乱した超音波エコーを受信するように超音波トランスデューサ・アレイを駆動することにより動作する。血流の特性の測定においては、反射された超音波を周波数基準と比較し、血球等の流動する散乱体によって反射波に対して付与された周波数シフトを決定する。この周波数シフト、即ち位相シフトは、血流の速度として解釈される。血液の速度は、特定のレンジ・ゲートにおけるファイアリング(firing)毎の位相シフトを測定することにより算出される。
【0004】
後方散乱された周波数の変化又はシフトは、血液がトランスデューサに向かって流れているときには増大し、血液がトランスデューサから遠ざかって流れているときには減少する。異なる色を用いて速度を表示し、大きさ及び方向を表現することもできる。また、ドプラ・スペクトルの全パワーを表示することもできる。カラー・フロー画像は、血液等の運動する物質の速度又はパワーのカラー画像を白黒の解剖学的構造のBモード画像に重畳することにより形成される。典型的には、カラー・フロー・モードは、そのすべてがBモード画像上に置かれた状態で各々のサンプル・ボリュームの速度又はパワーを表現するようにカラー符号化されている隣接する何百ものサンプル・ボリュームを同時に表示する。
【0005】
従来の超音波イメージング・システム(図1に示す)では、超音波トランスデューサ・アレイ2が作動されて、同一の送信特性で同一の送信焦点位置にフォーカスされている一連のマルチ・サイクル(典型的には4サイクル〜8サイクル)・トーン・バーストを送信する。これらのトーン・バーストは、一定のパルス繰り返し周波数(PRF)でファイアリングされる。PRFは典型的には、キロヘルツの範囲にある。同一の送信焦点位置にフォーカスされている一連の送信ファイアリングは「パケット(packet)」と呼ばれる。各々の送信ビームは、走査されている物体を通じて伝播し、血球等の物体内の超音波散乱体によって反射される。反射信号は、トランスデューサ・アレイの素子によって検出された後に、ビームフォーマ4によって受信ビームとして形成される。ビームフォーマは、遅延されたチャネル・データを加算し、RFデータか、又は同相及び直角位相(I/Q)データかを出力する。図1には後者の選択肢が示されている。従来のシステムでは、ビームフォーマの出力の周波数は、フィルタ6によってフィルタリングされた後に、動作モードに応じて、Bモード・プロセッサ8か、又は1対のウォール・フィルタ12とパラメータ推定器14とを含んでいるカラー・フロー・プロセッサかに出力される。
【0006】
Bモードでは、Bモード・プロセッサ8に組み込まれている包絡線検出器が、量(I2 +Q2 )1/2 を算出することによりビーム加算後の受信信号の包絡線を形成する。信号の包絡線は、エッジ強調及び対数圧縮のような更なる何らかのBモード処理を施されて表示データを形成し、この表示データがスキャン・コンバータ16に出力される。
【0007】
カラー・フロー・モードでは、I成分及びQ成分は、コーナ・ターナ・メモリ(corner turner memory)10に出力される。その目的は、インタリーブされている可能性のあるファイアリングからのデータをバッファリング(buffer)して、これらのデータを所与のレンジ・セルにおける複数のファイアリングにまたがる点から成るベクトルとして出力することにある。データは、「高速」で、即ち、各回のファイアリング毎に(ベクトルに沿って)ダウン・レンジの順で受信される。コーナ・ターナ・メモリの出力は、「低速」になるように再配列されており、即ち各々のレンジ・セル毎にファイアリングの順に再配列されている。結果として得られる「低速」のI信号及びQ信号のサンプルはカラー・フロー・プロセッサに出力される。カラー・フロー・プロセッサは、静止している組織又は極めて低速で動いている組織に対応するあらゆるクラッタを除去するそれぞれのウォール・フィルタ12を組み込んである。
【0008】
インソニファイ(insonifying )ビームと流れの軸との間の角度をθとすると、速度ベクトルの大きさは下記の標準的なドプラ方程式によって決定することができる。
v=cfd /(2f0 cos θ) (1)
ここで、cは血中での音速であり、f0 は送信周波数であり、fd は後方散乱した超音波における動きにより誘起されたドプラ周波数シフトである。
【0009】
ウォール・フィルタでフィルタリングされた出力はパラメータ推定器14へ供給され、パラメータ推定器14は、レンジ・セルの情報を中間自己相関パラメータN、D及びR(0)へ変換する。N及びDは、自己相関方程式の分子と分母であり、以下のように示される。
【0010】
【数1】
【0011】
ここで、Ii 及びQi はファイアリングiについての入力データであり、Mはパケット内のファイアリングの数である。R(0)は、あるパケットにおけるファイアリングの数の全体にわたる有限の和として近似することができ、以下のようになる。
【0012】
【数2】
【0013】
R(0)は、反射された超音波エコーのパワーを指示する。
パラメータ推定器14のプロセッサは、N及びDを各々のレンジ・セル毎に振幅及び位相へ変換する。用いられる式は、以下の通りである。
【0014】
【数3】
【0015】
パラメータ推定器14は、これらの振幅値及び位相値を処理して、パワー、速度及び乱れ(turbulence)の推定値とする。位相は、後に示すように、速度に比例する平均ドプラ周波数を算出するのに用いられ、R(0)及び|R(T)|(振幅)は、乱れ(ばらつき)を推定するのに用いられる。
平均ドプラ周波数は、N及びDの位相、並びにパルス繰り返し時間Tから得られる。
【0016】
【数4】
【0017】
平均速度は、ドプラ・シフト方程式を用いて算出される。
【0018】
【数5】
【0019】
なお、パラメータ推定器14は、平均ドプラ周波数を中間出力として算出せずに、ルックアップ・テーブルを使用することによりプロセッサの位相出力から直接
【0020】
【外1】
【0021】
一般的には、スキャン・コンバータ16に入力される表示データは、R−θフォーマットを有しており(セクタ走査の場合)、このデータは、スキャン・コンバータによって、ビデオ表示のためのX−Yフォーマットへ変換される。走査変換後のフレームはビデオ・プロセッサ18へ渡され、ビデオ・プロセッサ18は、ビデオ・データをグレイ・スケール及び/又はカラー・マッピングとして写像(マッピング)する。次いで、グレイ・スケール及び/又はカラーの画像フレームは、ビデオ・モニタ20へ送られて表示される。典型的には、速度若しくはパワーのいずれかが単独で表示されるか、又は速度がパワー若しくは乱れのいずれかと共に表示される。システム制御は、ホスト・コンピュータ(図1には示されていない)に集中されており、ホスト・コンピュータは、オペレータ・インタフェイス(例えば、キーボード)を介してオペレータの入力を受け取って、様々なサブシステムを制御する。
【0022】
超音波プローブが関心のある区域の全体を掃引していれば、2次元画像を蓄積して3次元データ・ボリュームを形成することができる。このボリューム内のデータは、ボリューム・レンダリング及び表面投影を含めたいくつかの方法で操作することができる。具体的には、最大、最小、平均又はメジアンのピクセル値をイメージング平面に投影することにより、Bモード(強度)データ、速度データ又はパワー・データの3次元画像を形成することができる。
【0023】
Bモード・データの3次元再構成は、コントラスト及びダイナミック・レンジの不足のためデータをセグメント化するのが困難であることから劣化する。速度データ及びパワー・データの再構成は、付随するクラッタ信号がセグメント化に誤差を生ずることから劣化する。従って、セグメント化が改善されている3次元イメージング手法が必要とされている。
【0024】
【発明の概要】
本発明は、セグメント化が改善された超音波データの3次元画像を表示する方法及び装置である。これは、造影剤を用いて又は用いないで解剖学的構造の調波イメージング(harmonic imaging)を行うことによって達成される。調波イメージングには2つの種類がある。(1)流体中に注入された造影剤から反射された調波のイメージング、及び(2)一般に「組織調波(tissue harmonics)」と呼ばれる自然に生起する調波のイメージングである。
【0025】
組織調波は、一般的には、造影剤からの調波よりも弱いが、サイド・ローブがより小さいこと且つ一般には基本周波数の倍数の高調波を検出するように選択が行われることから、よりよい分解能を提供する。多量の脂肪組織を有する患者は、一般的には、脂肪内での音速のばらつきのため、イメージングが困難である。しかしながら、音波の伝播時に増大した量の調波を発生するのが脂肪組織の特徴であるので、大柄の患者に関して有利になる。
【0026】
造影剤は、送信された超音波パルスの周波数を持つ反射信号を生じさせるばかりでなく、送信周波数の高調波(harmonic)又は低調波(subharmonic )の反射信号も生じさせるという独特の特性を有している。造影剤は、造影剤が到達した区域をイメージングする目的で体内に導入されるものである。従って、造影剤は、流動及び潅流を研究するのに有用である。
【0027】
本発明の好ましい実施例によれば、超音波トランスデューサ・アレイは、造影剤を含む所望のサンプル・ボリュームにフォーカスされている所定の送信中心周波数の超音波パルスによって形成されるビームを送信するように制御される。受信モードでは、受信器は、送信中心周波数の倍数又は約数の周波数で反射されたエコーをビーム加算後の受信信号として形成する。この処理は、物体ボリューム内の関心のある領域に交差する多数の走査平面の各々における各々のサンプル・ボリューム毎に繰り返される。
【0028】
受信信号内の不要な周波数をフィルタリング除去した後に、即ち所望の高(低)調波信号の実質的にすべてを通過させた後に、超音波スキャナは、高(低)調波信号を処理してBモード・データ、速度データ又はパワー・データを取得し、次いで、シネ・メモリにこのデータを収集して、ピクセル表示データのボリュームを形成する。3次元画像を形成するためには、様々な種類の投影(例えば、ピクセル投影、表面投影又は合成投影)を用いることができる。Bモード・データの場合には、最大、最小、平均又はメジアン・ピクセル投影を行うことができる。速度データ又はパワー・データの場合には、好ましくは平均又はメジアン・ピクセル投影が行われる。最終的な値は、各々の射線(ray )毎に取得される。一旦、投影が完了すると(即ち、すべての射線が検査されると)、投影が表示される。
【0029】
本発明によれば、高(低)調波信号は、実質的に造影剤(即ち、マイクルバブル)が存在している区域のみから反射される。これにより、Bモード・データのセグメント化が、高調波又は低調波信号を発生する区域の選択に単純化される。速度信号及びパワー信号は、大幅に増大するので、背景クラッタから関心のあるデータを分離するのがより容易になる。高(低)調波信号に加えて、音響刺激型音波射出(ASAE、acoustically stimulated acoustic emissions)をやはりセグメント化の目的に用いることができる。ASAEは、バブルの破裂に関連しており、小血管の流れのイメージングを可能にする造影剤現象である。
【0030】
【好ましい実施態様の詳しい説明】
本発明による調波イメージングは、所定の基本周波数でパルスを送信し且つ基本周波数の倍数又は約数の周波数で受信するようにシステムを設定することにより実現される。これは、基本周波数(例えば、2MHz)で音波を送信するように送信ビームフォーマの送信メモリをプログラムすること、並びに検出されるべき調波周波数(例えば、4MHz)を中心とする周波数帯域を通過させるようにビームフォーマ4(図1を参照)の各々の受信チャネル内のフィルタ及びフィルタ6を設定することを伴う。フィルタ6は、好ましくは、所望の高調波又は低調波周波数帯域を通過させるように設計された1組のフィルタ係数でプログラムされるFIRフィルタである。
【0031】
後に詳述する理由により、場合によっては、所望の高(低)調波周波数帯域と共に基本周波数帯域の一部を通過させるようにすると有利なことがある。例えば、場合によっては、第2高調波と、これと連続した基本周波数帯域の部分帯域(sub-band)とを包含する帯域の周波数を通過させると望ましいことがある。このことは、通過帯域を基本領域に延長するフィルタ係数でフィルタをプログラムすることにより達成することができる。
【0032】
引き続き図1を参照して述べると、帯域フィルタリングされた信号は、フィルタ6から、Bモード・プロセッサ8及びカラー・フロー・プロセッサ(12、14)の一方又は両方に出力される。次いで、処理後の信号は、走査変換されて、1フレーム分のピクセル表示データを取得する。3次元画像は、以下に詳述するように、シネ・メモリに記憶されている多数の画像フレームを再構成することにより形成される。
【0033】
図2を参照して述べると、システムの制御は、マスタ・コントローラ又はホスト・コンピュータ22に集中されており、ホスト・コンピュータ22は、オペレータ・インタフェイス(図示されていない)を介してオペレータの入力を受け取って、様々なサブシステムを制御する。ホスト・コンピュータ22はまた、システムのタイミング信号及び制御信号を発生する。スキャン・コンバータ16は、音響線メモリ24及びX−Yメモリ26を含んでいる。音響線メモリ24に極座標(R−θ)セクタ・フォーマットで記憶されているデータ(強度、速度又はパワー)は、適当にスケーリング(拡縮)されたデカルト座標ピクセル表示データへ変換され、このデータがX−Yメモリ26に記憶される。ピクセル表示データの多数の相次ぐフレームは、先入れ先出し方式でシネ・メモリ28に記憶されている。シネ・メモリは、ユーザに対して実時間で表示されるピクセル表示データを絶えず取得しながらバックグラウンドで動作している循環画像バッファのようなものである。ユーザがシステムをフリーズさせると、ユーザは、シネ・メモリ28内に過去に取得されているピクセル表示データを見ることが可能になる。
【0034】
従来の超音波イメージング・システムは、連続的にシネ・メモリ28に画像フレームを収集する。シネ・メモリ28は、単一の画像閲覧及び多数の画像ループ閲覧、並びに様々な制御機能のための常駐ディジタル画像記憶手段を提供する。単一画像のシネ再生時に表示される関心のある領域は、画像の取得時に用いられた関心のある領域である。シネ・メモリはまた、ホスト・コンピュータ22を介してディジタル大容量記憶装置(図示されていない)へ画像を転送するためのバッファとしても作用する。
【0035】
中央処理装置(CPU)30は読み出し専用メモリを有しており、該読み出し専用メモリは、取得されたイメージング・データ・ボリュームを、様々な角度でとられる多数の3次元投影画像へ変換するのに用いられるルーチンを記憶している。CPU30は、システム制御バス34を介して、X−Yメモリ26及びシネ・メモリ28を制御する。具体的には、CPU30は、X−Yメモリ26からビデオ・プロセッサ18及びシネ・メモリ28へのデータの流れ、並びにシネ・メモリからビデオ・プロセッサ22及びCPU30自体へのデータの流れを制御する。データの各々のフレームは、被検体を通る多数のスキャン又はスライスの1つを表していて、X−Yメモリ26に記憶されて、次のサイクルでビデオ・プロセッサ18及びシネ・メモリ28へ送信される。フレームのスタックは、走査されている物体ボリュームを表していて、シネ・メモリ28のセクション28Aに記憶される。初期化(図3の工程36)に応答して、CPU30は、シネ・メモリ・セクション28Aから、関心のある物体ボリュームに対応しているピクセル表示データのみを検索する。これは、各々の選択されたフレームから関心のある領域にあるデータのみを検索することにより達成される。多数の選択されたフレームの各々からの関心のある領域に対応しているデータは、関心のあるソース・データ・ボリュームを形成し、これはローカル・メモリに記憶される(図3の工程38)。
【0036】
従来の2次元超音波画像は、基本信号成分と調波信号成分との組合せから形成されており、後者の成分は、組織、又は造影剤を含んでいる血流のような非線形媒体において発生される。本発明によれば、3次元超音波画像のセグメント化が、基本信号成分のすべて又はかなりの部分を抑制すると共に実質的にすべての所望の調波信号成分を通過させることにより改善される。
【0037】
本発明の1つの好ましい実施例によれば、関心のある物体ボリュームへの造影剤の注入が、受信信号の調波帯域フィルタリングと組み合わされる。造影剤は、医療超音波法については、従来イメージングが困難であるとされている脈管の解剖学的構造の診断を助けるために開発されてきた。例えば、造影剤の使用は、Ultrasonics 誌、第29巻、第324頁〜第380頁(1991年)のde Jong 等による「超音波造影剤の原理と最近の発展(Principles and Recent Developments in Ultrasound Contrast Agents)」で議論されている。これらの造影剤は典型的には、直径が1μm〜10μmの範囲にあるマイクロバブルであり、血流内に注入される。マイクロバブルの後方散乱信号は血球の後方散乱信号よりも遥かに大きいので、マイクロバブルは、血流のイメージングを可能にするマーカとして用いられる。これらの造影剤からエコーを更に分離する1つの方法は、造影剤エコーの高(低)調波成分を用いるものであり、この成分は、造影剤が存在していない周囲組織の調波成分よりも遥かに大きい。[例えば、Proc. 1992 IEEE Ultrason. Symp. 誌、第1175頁〜第1177頁のNewhouse等の「第2高調波ドプラ超音波血液潅流測定(Second Harmonic Doppler Ultrasound Blood Perfusion Measurement)」、及びProc. 1994 IEEE Ultrason. Symp. 誌、第1547頁〜第1550頁のBurns 等の「マイクロバブル造影剤を用いた調波パワー・モード・ドプラ:小血管のフロー・イメージングの改良された方法(Harmonic Power Mode Doppler Using Microbubble Contrast Agents: An Improved Method for Small Vessel Flow Imaging )を参照されたい。]高(低)調波信号のフロー・イメージングは大部分が、周波数f0 の狭帯域信号を送信し、周波数2f0 (第2高調波)又は周波数f0 /2(低調波)を中心とする帯域幅で受信してから従来の処理を行うことにより行われている。
【0038】
造影剤の注入に依らない好ましい実施例によれば、高(低)調波を発生する人体組織からの受信信号を帯域フィルタリングして、後続の3次元再構成のために所望の(例えば、第2高)調波を通過させる。
両方の実施例とも、高(低)調波周波数を中心とする帯域フィルタは、所望の高(低)調波成分を実質的に分離する。例えば、フィルタ6(図1を参照)は好ましくは、MタップのFIRフィルタである。このMタップFIRフィルタが受信信号の所望の高(低)調波周波数の大部分を通過させるように、1組の高(低)調波フィルタ係数c0 ,c1 ,…,cM-1 が選択される。送信される中心周波数がf0 にあるならば、組織及び造影剤の非線形性から、kf0 の高調波が発生される。ここで、kは、2よりも大きい又は2に等しい整数である。また、造影剤のバブルの破裂により、周波数f0 /kの低調波が発生され得る。好ましい実施例では、第2高調波信号がイメージングされる。送信にあたっては、送信信号は、第2高調波がトランスデューサの帯域幅の範囲内にも適合可能となるように十分に狭い帯域幅を有していなければならない。
【0039】
3次元超音波イメージングを行うためには、ピクセル投影はCPU30において再構成され、CPU30は米国特許第5,226,113号に開示されているレイ・キャスティング(ray-casting )・アルゴリズムを用いて一連の変換を実行する。相次ぐ変換は、例えば+90°〜−90°の角度範囲内で、例えば10°の間隔の角度増分で行われる最大、最小、平均又はメジアン速度投影を表す。但し、角度増分は10°でなくてもよいし、本発明は、いかなる特定の角度範囲にも限定されていない。
【0040】
本発明に用いられるレイ・キャスティング手法によれば、超音波トランスデューサで物体ボリューム52を走査することにより、サンプル50(図4を参照)の立体的にレンダリングされた投影画像が、任意の視角、例えば角度パラメータ(θ,φ)によって表される球面投影角から表示される。ここで、θは、視線58の補助線(extension) 58′がX−Y平面上で作る角度であり、φは、補助線58′に対する視線58の角度である。サンプル・ボリューム52は、同一個数の物体ボリューム要素(ボクセル)OVをそれぞれ含んでいる積み重なった連続的なスライス又はシートOS1 ,OS2 ,…,OSk の系列を形成するような方式で走査される。各々のボクセルは、シート平面(例えば、X−Y平面)内で矩形の輪郭を有しており、縦横の各辺(complementary sides )は等しい長さSであってよく、その場合はこの輪郭は正方形となり、シート厚Tは一般的にはいずれの辺の長さにも等しくない。このように、第1の物体スライスOS1 は、第1の多数の物体ボクセルOVi,j,1 を含んでおり、ここで、i及びjはそれぞれ、ボクセルのX軸及びY軸での位置である。同様に、第2の物体スライスOS2 は、物体ボクセルOVi,j,2 を含んでいる。任意の物体スライスOSk は、ボクセルOVi,j,k を含んでおり、ここで、kはこのボクセルのZ軸での位置である。
【0041】
各々の物体ボクセルOVi,j,k は分析され、そのデータ値(強度、速度又はパワー)が、データ・ボリューム54の対応するデータ・ボクセルDVi,j,k に配置される。各々の物体スライスOSk の厚み及び各々の物体ボクセル面のサイズ(ボクセルのZ−Y平面におけるサイズ)が一般に同じではなくても、データ・ボリューム54は単純なi,j,k立方格子である。即ち、物体ボリュームは、各々のボクセルについて異なるX、Y及びZの寸法を有していてもよいばかりでなく、任意の次元におけるボクセルの総数も同じである必要はない。例えば、典型的な超音波3次元走査は、各々のスライスがボクセルの256×256のマトリクスを含むようにして、128のスライスを含み得る。
【0042】
CPU30によって用いられる公知の手法によれば、物体50の画像は、データ・ボクセルDVi,j,k の格子点から画像平面56に向かってのレイ・キャスティングにより投影される(図3の工程42)。便宜を図るために、格子点は、例えば、データ・ボリュームの原点に最も近いデータ・ボクセルの頂点としてよい。投射される射線62は、物体ボリューム52を眺める球面角度パラメータ(θ,φ)から変換された球面角度パラメータ(α,β)による投影角度でデータ・ボリューム54を出発する。これら2つの角度は、非立方体である物体ボリューム52に対して立方体であるデータ・ボリューム54を用いることによる幾何学的な歪みのため、同一ではない。しかしながら、投射される射線62は、データ
【0043】
【外2】
【0044】
(後述する)によって、所望の視角(θ,φ)で物体ボリューム52を眺めることに対応するように決定される(球面座標での演算を仮定している)。射線62の各々は、データ・ボリュームのボクセルの格子点から画像平面に向かって投射される。
すべての射線62が画像平面の何らかの部分に入射するが、考察されている画像平面ピクセル60aの範囲内に到達した射線のみが、この画像平面ピクセルのデータに寄与することを許される。このように、物体ボリューム52の眺めている部分と、この選択された物体ボリュームを眺める視角(θ,φ)とを選択すると、データ・ボリュームの対応する部分の各々のボクセルのデータ値が、画像平面56に向かって何らかの角度(α,β)で投射される(物体ボリュームに関して歪んだデータ・ボリュームを眺めることに対応する。)。第1のボクセル(例えば、ボクセルDVi,1,k )のデータ値はこのようにして、選択されたθ及びφの値に従って、射線62aに沿って逆投影される。この射線62aは、ピクセル60a内の位置64aにおいて画像平面56に入射し、このピクセルに入射するのはこれが最初の射線であるので、入射データの値が所望のピクセル60aに帰属される(記憶される)。データ・ボリュームの次のボクセル(例えば、ボクセルDVi,2,k )は、ボクセルの格子点から同一の角度(α,β)構成で投影されるその関連した射線62bを有しており、画像平面67上でのその位置64bが記録される。入射位置64bが、所望のピクセル60aの範囲内にあるとすると、この第2の投影後の値は(最大ピクセル投影の場合には)、現在記憶されている第1の値と比較されて、より大きい方の値がピクセル60aについて記憶装置に配置される。選択されたデータ・ボリュームの各々のボクセルが画像平面56に向かって順に入って投影されるにつれて、データ・ボリュームのボクセル(例えば、ボクセルDVi,3,k )は遂には、その関連する射線62pに沿って投影されると所望のピクセル60aの範囲内に入射しなくなるので、そのデータ値は、ピクセル60aについて現在記憶されているデータ値とは比較されない。この時点で、ピクセル60aについての最大データ値が、特定の3次元視角(θ,φ)におけるデータのこの投影について設定される。但し、射線62pは実際には、入射点64pを有しており、入射点64pは他の画像平面ピクセル(例えば、ピクセル60b)の範囲内に位置しているので、この位置に記憶されているデータ値と比較され、比較の後に、より大きい方の値がこのピクセルについて記憶装置に返される。新たな投影を取得するときには、すべてのデータ値はゼロにリセットされる。このように、画像平面ピクセルの各々が画像投影手順の開始時にリセットされ、データ・ボリュームのボクセルのすべてが(選択された物体ボリューム52の部分によって設定されているボリュームの全体で又は選択された部分で)個別に且つ順に走査される。各々のデータ・ボクセルDVのデータ値が、関連する射線62にとり投影されて、画像平面56内の1つのピクセル60に入射し、ここで、各々のピクセルの最大値を射線投射されたデータ・ボリュームのボクセルの現在の値との間で比較してそのうちの大きい方を決定し、次いで、この大きい方の値が最大値画像の一部として記憶される。実際には、最大ピクセル投影の場合には、新たに投射される射線が入射する画像平面のピクセルにおいて既に記憶されているデータ値よりも、新たに投射されるデータ・ボクセル値の方が大きい場合にのみ、記憶されている最大値が変更される。
【0045】
上述の手法のもう1つの面によれば、データ投影がスケーリングされて(図3の工程44)、物体ボリュームと画像平面との間のあらゆる異方性が、逆投影の完了後にただ単一の組の計算を行うことによって除去される。ここで図5を参照すると、物体ボリューム52は実在のボリュームであるのに対し、データ・ボリューム54は抽象概念であるので、第1の平面内での、異なる角度γにおける立方体のデータ・ボリューム格子54の存在と、任意の視線方向66が物体ボリューム52及びデータ・ボリューム54の両方に関して配置される角度Ψとによるデータ投影の歪みの量を決定することが必要である。各々のボクセルの見かけ上の寸法は、実効の仰角Ψ及びγが変化するにつれて変化してゆく。縦横(アスペクト)比A(物体ボリューム52における実際のスライス厚Tと同じ物体ボリューム52における実際のピクセル・サイズSとの比として定義される)が1でない(即ち、物体ボクセルが、データ・ボリューム54の場合に見られるように立方体ボクセルではないとき、1よりも大きいか又は小さくなる)場合、仰角Ψ及びγは異なることになり、データ・ボリュームにおける実効的な仰角Ψは、物体ボリュームにおける実際の仰角γと異なることになる。データの回転は、以下の式によって得られる物体の仰角に従う。
【0046】
Ψ=tan-1[(1/A)tan(γ)]
この後に、仰角スケーリング・ファクタをすべての投影後のデータの高さに乗算することにより、投影後のデータを物体ボリュームにおいて正しい高さ(回転が水平軸の周りで行われる場合)を有するようにスケーリングすることができる。古い投影後の画像の高さHは、実効スケーリング・ファクタEs によって補正することができ、ここで、
Es =[(Acosγ)2 +sin2 γ]1/2
であり、新たな高さはH′=H・Es である。回転が垂直軸の周りで行われる場合には、幅についても同じことが当てはまる。
【0047】
上述の関係を利用すると、歪みは1つの軸のみに沿ったものであるので、データ・ボリュームの角度(α,β)の回転が角度(θ,φ)となり、すると、角度θが角度αと等しくなる。3×3の回転行列[M]の各要素を決定することができ、2つの関連する回転角度が与えられると、これらの関係を用いて、データ・ボリュームから画像平面への変換が決定される。
【0048】
X′=M1X+M2Y+M3Z=XO
Y′=M4X+M5Y+M6Z+YO
ここで、M1〜M6は回転行列の最初の2行であり(即ち、M1=−sinθ、M2=cosθsinΨ、M3=0、M4=−cosθsinΨ2、M5=−sinθsinΨ、及びM6=cosΨ)、X′及びY′は投影された点の画像平面上での位置であり、XO及びYOは画像平面の選択された部分が開始する位置での画像平面のX及びYのオフセット(それぞれX及びYの最低値の点を基準とする)である。データが画像平面56上へ投影された後に、異方性の物体ボクセルの影響を補正するように画像をスケーリングする。ファクタM1〜M6は、投影の開始時(所与のθ及びφ)に予め算出しておいて(図3の工程40)、すべての回転計算に用いることができることが理解されよう。
【0049】
図6は、ホスト・コンピュータ22(又は独立した専用のプロセッサ)に組み込まれる上述のレイ・キャスティング手法を実行するシステムを示している。このようなシステムは、シネ・メモリ24からのデータ入力70aにおいて受け取られるスライス・データを記憶する3次元データ・メモリ70を含んでいる。各々の物体ボクセルに関連するデータは、CPU74からのボクセル・アドレス入力70bにおいて受け取られるボクセル・アドレス入力情報に応答して、このボクセルのアドレスに記憶されている。一旦、データ・メモリ手段が充填されたら(物体ボリューム52からデータ・ボリューム54へのすべての所要データの転送に対応する)、関心のある物体ボリュームの部分が選択され、その開始頂点と、X、Y及びZ方向での範囲とを設定するデータが、CPU74からアドレス発生器72の入力72aへ送られる。アドレス発生器72は、アドレス出力72bにおいて、選択された物体ボリューム内の各々のボクセルのX,Y,Zアドレスを順に供給する。出力72bは、データ・メモリ70の出力データ・アドレス入力70cに接続されているので、この時点でアドレス指定されたこの1つのボクセルについて記憶されているデータが、データ・メモリ出力70dから出力される。ボクセルのX,Y,Zのアドレスの系列はまた、回転パラメータ・コンピュータ76の第1の入力76aにも供給され、回転パラメータ・コンピュータ76は、CPU74を介して角度(α,β)情報を算出された行列要素M1〜M6の値として受け取って、出力76cにおいて、選択された視角(θ,φ)で眺めたときの物体のこのX,Y,Zのピクセルに対応する画像平面ピクセルのアドレスX′,Y′を供給する。視角(θ,φ)情報は、システムに入って、CPU74によって処理される。この結果は、視角行列78の入力78b及び78cに入って、その出力78aに行列要素M1〜M6を供給し、そこから回転コンピュータ76へ供給する。画像平面のピクセル・アドレスX′,Y′は、画像平面メモリ80として動作するフレーム・バッファのアドレス入力80aに現れる。同時に、データ・ボリュームから投影平面への投影後のピクセル表示データが、3次元データ・メモリの出力70dから、画像平面メモリの新規データ入力80bに現れる。このデータはまた、データ比較器82の新規データ入力82aにも現れる。入力80aのこのアドレスについて画像平面メモリ80に過去に保存されているピクセル表示データは、旧データ出力80cに現れ、そこから比較器の旧データ入力82bに現れる。入力82b及び82aにおける旧データ及び新規データはそれぞれ、比較器82で比較され、入力82aの新規データが入力82bの旧データよりも大きさが大きければ、比較器82の出力82cが、ある選択された論理条件(例えば、高論理レベル)にイネーブルされる。出力82cは、画像平面メモリの置換制御データ入力80dに接続されており、置換データ制御入力80dが選択された論理レベルにある場合には、入力80aによって制御されるアドレスに記憶されているデータが変更されて、入力80bの新たなデータを受け入れるようにする。このように、記憶されているデータは、(CPU74からの)データ/制御ポート80eを介した信号等によって最初にリセットされ、新たなデータが過去に記憶されている古いデータの値を上回ったことを指示する比較に応答して、各々の画像平面ピクセル位置X′,Y′毎に最大値のデータが記憶される。すべての選択されたアドレスがアドレス発生器72によって順に走査された後に、画像平面メモリ80に記憶されているデータをCPU74においてスケーリングして、スケーリング後の画像平面データを、表示、永久保存又は類似の目的のためにメモリ80から引き出すことができる。
【0050】
速度データ及びパワー・データについては、好ましくは、平均又はメジアン・ピクセル投影が用いられる。図7は、平均又はメジアン強度投影を行うシステムを示している。図7に示されているシステムは、加算器84がデータ比較器82に代えて用いられ且つトリガ86が追加されている点で、図6に示すシステムと異なっている。また、図7の理解のためには、「新規データ」とラベルされた線が、図6に示す3次元データ・メモリ70の出力に接続されており、X′,Y′とラベルされた線が図6に示す回転パラメータ・コンピュータ76の出力に接続されているものと理解されたい。
【0051】
平均値投影の場合には、カレントの投影後データ・ボクセルの値が、この投影射線の入射した画像平面のピクセルに既に記憶されている和に加算され、その和は最終的にこのピクセルについてのこのような入射線のカウント数で除算される。選択されたデータ・ボリュームにおける各々のボクセルが画像平面56に向かって順に入って投影されるにつれて、データ・ボリュームのボクセル(例えば、ボクセルDVj,3,k )は遂にはその関連する射線62pに沿って投影されて、所望のピクセル60aの範囲内に入射しなくなり、これにより、そのデータ値は、ピクセル60aについて現在記憶されている和の値には加算されない。次いで、ピクセル60aについての平均データ値が、特定の3次元視角(θ,φ)におけるデータのこの投影について算出される。但し、射線62pは実際には、他の画像平面ピクセル(例えば、ピクセル60b)の範囲内にある入射点64pを有しているので、そこに記憶されている和のデータ値に加算される。新たな投影を取得するときには、すべてのデータ値はゼロにリセットされる。このように、画像平面ピクセルの各々は、画像投影手順の開始時にリセットされ、データ・ボリュームのボクセルのすべて(選択されている物体ボリューム52の部分によって設定されている空間の全体で又は選択された部分で)個別に且つ順に走査される。各々のデータ・ボクセルDVのデータ値が、関連する射線62を介して投影されて、画像平面56の1つのピクセル60において画像平面56に入射する。次いで、各々のピクセルについて平均値若しくはメジアン値又はその両方が算出されて表示される。
【0052】
各々の画像平面ピクセルの平均強度は、以下の方法によって求めることができる。各々の物体ボリューム・ボクセル・アドレスX,Y,Zについて、加算器84のデータ入力84aに新たなデータが印加される。同時に、アドレスX′,Y′が、画像平面メモリ80のアドレス入力80aに現れる。新たなアドレス・データの出現は、トリガ86によって感知されて、入力80fにおいて置換データ信号が供給される。この信号は、画像平面メモリ80のこのアドレスに過去に記憶されているデータが検索されて、画像平面メモリ80のデータ出力80cを介して加算器84のデータ入力84bに印加されるように十分に遅延して発生する。加算器は、新たなデータと記憶されているデータとの和を形成する。この和は、画像平面メモリ80の和データ入力80dに接続されている加算器84の出力84cにおいて利用可能になる。関心のあるボリュームの物体ボリュームのボクセルのすべてが走査されて、画像平面の関連するピクセルに対するそれらのすべてのボクセルの寄与が為されてしまうまで、和データは、画像平面のこのピクセルについて画像平面メモリ80のアドレスX′,Y′に記憶されている。次いで、画像平面ピクセルの和の値は、各々の画像平面ピクセルについての加算の回数(これは、各々のピクセル毎に入力80dの作動を加算することにより得られる数であり、やはり画像平面メモリ80に記憶されている)に従って、画像平面メモリ80のデータ・ポート80eを介してCPU74によって演算され、記憶及び後続の表示のために各々の画像平面ピクセル毎の平均強度が導出される。カラー・フロー・モードでは、平均ピクセル値は、反射したドプラ信号に含まれる平均速度又はフロー・パワーのいずれかを表す。
【0053】
代替的な実施例によれば、画像フレームにおける各々のピクセル毎の投影値のそれぞれの一覧を記憶し、次いで、関心のあるボリュームにおける物体ボリュームのボクセルのすべてが走査された後に、各々のピクセル毎にそれぞれのメジアン・ピクセル値を識別することにより、メジアン・ピクセル投影が実行される。次いで、これらのメジアン値を表示する。メジアン・ピクセル投影については、出力70dからの新たなデータは、画像平面メモリ手段80の新規データ入力80bに入力されることが当業者には理解されよう。これらのデータ値は、CPU74によってメジアン表に配列され、CPU74は後に、表に含まれている値からメジアン値を識別する。
【0054】
更なる変形によれば、各々のピクセル投影について、平均値及びメジアン値の両方が算出され記憶される。次いで、選択アルゴリズムを各々のピクセル毎の平均値及びメジアン値に適用して、表示されるべき投影画像にどちらの値を用いるかを決定する。用途によっては、2つの値のうちの大きい方を選択するような選択アルゴリズムを用いてもよい。また、他の用途では、2つの値のうちの小さい方を選択するような選択アルゴリズムを用いてもよい。
【0055】
本発明によれば、図3に示す方法は、シネ・メモリから検索される関心のあるデータ・ボリュームについての強度データ、速度データ又はフロー・パワー・データに適用される。投影画像の各々のピクセルは、所与の画像平面への投影によって導出されるそれぞれの変換後のデータを含んでいる。加えて、シネ・メモリがオペレータによってフリーズされた時点で、CPU30は、シネ・メモリ28のセクション28B内の多数の連続したアドレスに、X−Yメモリ26からの最新のフレームを記憶している。第1の投影視角についての投影後のピクセル表示データは、シネ・メモリ・セクション28Bの第1のアドレスに書き込まれるので、関心のある領域の投影後のピクセル表示データは背景フレームに重畳される。この処理は、すべての投影後の画像がシネ・メモリ・セクション28Bに記憶されるまで各々の角度増分毎に繰り返される。尚、各々の投影後の画像フレームは、変換後のデータを含んでいる関心のある領域と、選択により、関心のある領域を包囲する背景の周囲辺とから成っており、背景の周囲辺は、関心のある領域の変換後のデータで上書きされていない背景フレーム・データから成っている。背景画像によって、各々の表示投影がどこから眺めたものであるかがより明瞭になる。次いで、オペレータは、表示のために任意の投影後の画像を選択することができる。加えて、投影後の画像の系列を表示モニタ上で再生して、物体ボリュームが観察者の目前で回転しているかのように物体ボリュームを表現することができる。調波分離(帯域)フィルタを介して基本信号の一部が通過し得るようにするにより、調波反射をしない解剖学的標識点を医療診断中にイメージングして観察することができる。
【0056】
画像は、連続的にシネ・メモリ内に収集される。フリー・ハンド式走査法によって、又はプローブを移動させるシステムを用いることにより、プローブで解剖学的構造のある区域の全体を掃引すると、3次元ボリュームが得られる。各々の画像の次の画像に関する位置は、多くの方法で決定することができる。プローブは、プローブを直線的に移動させ、又はある角度にわたってプローブを揺動させ、或いはプローブ面に垂直にある角度にわたってプローブを回転させる等の何らかの方法によって走査させることができる。プローブが既知の距離にわたって、又は既知の角度範囲にわたって一定の速度で平行移動されていれば、各々の画像の間の距離は容易に決定することができる。プローブが平行移動した距離は多数の方法で決定することができ、例えば、適当な測定装置を用いて距離を単純に測定することにより、又は解剖学的構造の上若しくは内部のマーカを利用することにより、或いはプローブに位置センサを取り付けることにより決定することができる。また、プローブを任意の経路を通じて移動させて、プローブに装着された位置センサからのデータを用いて各々の画像の位置を決定することもできる。次いで、これらの画像をデカルト座標系に配置することができる。加えて、プローブは、既知の経路を通してプローブを移動させるような装置に固定してもよい。
【0057】
ここで用いる際には、「高(低)調波モード」という用語は、基本送信周波数の倍数である周波数を有している高調波モードと、基本送信周波数の約数である周波数を有している低調波モードとを包含する。
本発明のいくつかの好ましい特徴についてのみ例示し説明したが、当業者には多くの改変及び変形が想到されよう。従って、特許請求の範囲は、発明の要旨に含まれるこれらのようなすべての改変又は変形を網羅しているものとする。
【図面の簡単な説明】
【図1】実時間超音波イメージング・システム内の主要な機能的サブシステムを示すブロック図である。
【図2】ピクセル表示データの相次ぐ三次元投影を有するフレームを再構成するシステムを示すブロック図である。
【図3】ピクセル表示データの相次ぐ三次元投影を有するフレームを再構成するアルゴリズムの各工程を示す流れ図である。
【図4】逆レイ・キャスティング投影を三次元レンダリングする際における関連するサンプリングされた関心のある物体ボリューム、関連するデータ・ボリューム及び画像投影平面を示す概略線図である。
【図5】物体ボリューム及びデータ・ボリュームの類似した眺めに対応しており、3次元超音波イメージングにおける必要なスケーリング定数を定義するのに有用な1対の2次元幾何構成の概略図である。
【図6】3次元超音波イメージングにおいて最大強度投影を行うシステムのブロック図である。
【図7】3次元超音波イメージングにおいて平均又はメジアン強度投影を行うシステムのブロック図である。
Claims (6)
- 超音波散乱体を含んでいる物体ボリュームの3次元イメージングのためのシステムであって、
基本周波数を有している超音波ビームを送信して、多数のサンプル・ボリュームにおいて前記物体ボリュームから反射された超音波エコーを検出する超音波トランスデューサ・アレイ(2)と、
表示モニタ(20)と、
前記サンプル・ボリュームから反射された前記超音波ビームの超音波エコーから導出される高(低)調波モード・データを取得する高(低)調波モード取得手段(4,6)であって、各々の高(低)調波モード・データが、前記多数のサンプル・ボリュームのそれぞれのサンプル・ボリュームに対応しており、前記高(低)調波モードが前記基本周波数の倍数又は約数である高(低)調波モード取得手段(4,6)と、
前記多数のサンプル・ボリュームの各々について、高(低)調波モード・データを記憶するメモリ(28)と、
前記メモリ手段に記憶されている前記高(低)調波モード・データからソース高(低)調波モード・データ・セット(34)を検索する手段と、
第1の画像フレームの各々のピクセルについて投影値を決定するコントローラであって、これらの投影値が、所定の関数に従って、前記ソース高(低)調波モード・データ・セット内の前記高(低)調波モード・データを第1の画像平面に投影することにより算出されて、これにより第1の投影後のデータ・セットを形成するコントローラ(22)と、
前記表示モニタ上に、前記第1の投影後のデータ・セットの少なくとも一部の関数である第1の投影画像を表示させるプロセッサと、
を具備し、
前記高(低)調波モード・データはBモード強度データであり、
前記高(低)調波モード取得手段は、所望の高(低)調波周波数帯域の実質的にすべてを通過させる帯域フィルタ(6)を有し、
この帯域フィルタ(6)が、前記基本周波数を中心とする基本周波数帯域のうち前記所望の高(低)調波周波数帯域と連続している部分帯域を通過させることを特徴とする3次元イメージングシステム。 - 前記所定の関数は最大化を行うものである請求項1に記載のシステム。
- 前記所定の関数は最小化を行うものである請求項1に記載のシステム。
- 前記コントローラは更に、第2の画像フレームの各々のピクセルについて、前記所定の関数に従って、前記第1の画像平面に対して回転されている第2の画像平面に前記ソース高(低)調波モード・データ・セット内の前記高(低)調波モード・データを投影することにより算出される投影値を決定し、これにより第2の投影後のデータ・セットを形成するように構成されており、前記プロセッサは更に、前記表示モニタ上に、前記第2の投影後のデータ・セットの少なくとも一部の関数である第2の投影画像を表示するように構成されている請求項1に記載のシステム。
- 超音波散乱体を含んでいる物体ボリュームの3次元イメージングのための方法であって、
超音波トランスデューサを送信モードに駆動して超音波ビームを送信する工程と、
前記超音波トランスデューサを受信モードにおいて、前記超音波ビームが多数のサンプル・ボリュームにおいて反射された超音波エコーを検出する工程と、
前記超音波エコーから導出される高(低)調波モード・データを取得する工程であって、各々の高(低)調波モード・データが前記多数のサンプル・ボリュームのそれぞれのサンプル・ボリュームに対応しており、前記高(低)調波モードが前記基本周波数の倍数又は約数である、高(低)調波モード・データ取得工程と、
前記多数のサンプル・ボリュームの各々について、高(低)調波モード・データを記憶する工程と、
該記憶された高(低)調波モード・データからソース高(低)調波モード・データ・セットを検索する工程と、
第1の画像フレームの各々のピクセルについて、所定の関数に従って、第1の画像平面に前記ソース高(低)調波モード・データ・セット内の前記高(低)調波モード・データを投影することにより算出された投影値を決定して、これにより第1の投影後のデータ・セットを形成する工程と、
該第1の投影後のデータ・セットの少なくとも一部の関数である第1の投影画像を表示する工程と、
を具備し、
前記高(低)調波モード・データはBモード強度データであり、
前記高(低)調波モード・データ取得工程は、所望の高(低)調波周波数帯域の実質的にすべてを通過させるように帯域フィルタリングする工程を含み、
この帯域フィルタリングする工程は、前記基本周波数を中心とする基本周波数帯域のうちの、前記所望の高(低)調波周波数帯域と連続している部分帯域を通過させる工程を更に含むことを特徴とする超音波散乱体の3次元イメージング方法。 - 第2の画像フレームの各々のピクセルについて、前記所定の関数に従って、前記第1の画像平面に対して回転されている第2の画像平面に前記ソース高(低)調波モード・データ・セット内の前記高(低)調波モード・データを投影することにより投影値を算出し、これにより第2の投影後のデータ・セットを形成する工程と、
該第2の投影後のデータ・セットの少なくとも一部の関数である第2の投影画像を表示する工程と、
を更に含んでいる請求項5に記載の方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/065,213 US6102858A (en) | 1998-04-23 | 1998-04-23 | Method and apparatus for three-dimensional ultrasound imaging using contrast agents and harmonic echoes |
US09/065213 | 1998-04-23 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000152939A JP2000152939A (ja) | 2000-06-06 |
JP4767377B2 true JP4767377B2 (ja) | 2011-09-07 |
Family
ID=22061104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP11569099A Expired - Fee Related JP4767377B2 (ja) | 1998-04-23 | 1999-04-23 | 物体ボリュームの3次元イメージングのためのシステム及び方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US6102858A (ja) |
EP (1) | EP0952463B1 (ja) |
JP (1) | JP4767377B2 (ja) |
DE (1) | DE69936418T2 (ja) |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6110120A (en) | 1997-04-11 | 2000-08-29 | Acuson Corporation | Gated ultrasound imaging apparatus and method |
US6231512B1 (en) * | 1999-05-28 | 2001-05-15 | General Electric Company | Method and apparatus for parametric harmonic imaging |
US6290647B1 (en) * | 1999-07-02 | 2001-09-18 | Acuson Corporation | Contrast agent imaging with subharmonic and harmonic signals in diagnostic medical ultrasound |
US6533727B1 (en) * | 2000-02-11 | 2003-03-18 | Koninklijke Phillips Electronics N.V. | Ultrasonic system and method employing non-integer harmonic echo signals for imaging |
ITSV20000027A1 (it) * | 2000-06-22 | 2001-12-22 | Esaote Spa | Metodo e macchina per l'acquisizione di immagini ecografiche in particolare di tipo tridimensionale nonche' sonda di acquisizione |
EP1358849B1 (en) * | 2001-02-01 | 2011-04-20 | Hitachi Medical Corporation | Ultrasonic contrast-enhanced imager |
US6626836B2 (en) | 2001-04-04 | 2003-09-30 | Siemens Medical Solutions Usa, Inc. | Adaptive signal processing scheme for contrast agent imaging |
US6638230B2 (en) * | 2001-07-31 | 2003-10-28 | Koninklijke Philips Electronics N.V. | Apparatus and method of frequency compounding to perform contrast imaging |
US6450962B1 (en) * | 2001-09-18 | 2002-09-17 | Kretztechnik Ag | Ultrasonic diagnostic methods and apparatus for generating images from multiple 2D slices |
GB2384310A (en) * | 2002-01-19 | 2003-07-23 | Stichting Tech Wetenschapp | Ultrasonic tissue imaging using high order harmonics |
US6705993B2 (en) | 2002-05-10 | 2004-03-16 | Regents Of The University Of Minnesota | Ultrasound imaging system and method using non-linear post-beamforming filter |
US8221321B2 (en) | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods for quantification and classification of fluids in human cavities in ultrasound images |
GB2391625A (en) * | 2002-08-09 | 2004-02-11 | Diagnostic Ultrasound Europ B | Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams |
US7819806B2 (en) | 2002-06-07 | 2010-10-26 | Verathon Inc. | System and method to identify and measure organ wall boundaries |
US7520857B2 (en) * | 2002-06-07 | 2009-04-21 | Verathon Inc. | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume |
US8221322B2 (en) | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods to improve clarity in ultrasound images |
US6589177B1 (en) | 2002-11-15 | 2003-07-08 | Koninklijke Philips Electronics N.V. | Method and apparatus for obtaining B-flow and B-mode data from multiline beams in an ultrasound imaging system |
JP3964364B2 (ja) * | 2003-07-22 | 2007-08-22 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
JP2007508889A (ja) * | 2003-10-23 | 2007-04-12 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 超音波画像化方法及び装置 |
US20050110793A1 (en) * | 2003-11-21 | 2005-05-26 | Steen Erik N. | Methods and systems for graphics processing in a medical imaging system |
US20060058674A1 (en) * | 2004-08-31 | 2006-03-16 | General Electric Company | Optimizing ultrasound acquisition based on ultrasound-located landmarks |
JP4268976B2 (ja) * | 2006-06-15 | 2009-05-27 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | イメージング装置 |
JP4881112B2 (ja) * | 2006-09-19 | 2012-02-22 | 株式会社東芝 | 超音波診断装置及び画像データ生成方法 |
EP2097009A4 (en) * | 2006-12-29 | 2010-01-06 | Verathon Inc | SYSTEM AND METHOD FOR ULTRASONIC HARMONIC IMAGING |
US8167803B2 (en) | 2007-05-16 | 2012-05-01 | Verathon Inc. | System and method for bladder detection using harmonic imaging |
EP2323559A4 (en) | 2008-08-07 | 2016-09-21 | Verathon Inc | DEVICE, SYSTEM AND METHOD FOR MEASURING THE DIAMETER OF GROUPS OF AQUARIUS ANEUTS |
DE102009015386A1 (de) * | 2009-03-27 | 2010-09-30 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Erzeugung eines funktionellen Datensatzes eines perfundierten Bereichs des menschlichen oder tierischen Körpers |
US8159680B2 (en) * | 2010-02-16 | 2012-04-17 | Massachusetts Institute Of Technology | Single-transducer, three-dimensional laser imaging system and method |
KR101526486B1 (ko) * | 2013-07-24 | 2015-06-10 | 디게이트 주식회사 | 광학 파라미터를 이용한 3차원 객체 렌더링 장치 |
US11166700B2 (en) * | 2015-03-25 | 2021-11-09 | Koninklijke Philips N.V. | Optimal ultrasound-based organ segmentation |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4276885A (en) * | 1979-05-04 | 1981-07-07 | Rasor Associates, Inc | Ultrasonic image enhancement |
FR2653919A1 (fr) * | 1989-10-30 | 1991-05-03 | Gen Electric | Procede et appareil pour la restitution de projections volumetriques utilisant la projection inverse de rayons. |
US5226113A (en) * | 1989-10-30 | 1993-07-06 | General Electric Company | Method and apparatus for volumetric projection rendering using reverse ray casting |
JPH0675576B2 (ja) * | 1991-03-13 | 1994-09-28 | アロカ株式会社 | 超音波画像表示装置 |
JP3046424B2 (ja) * | 1991-11-11 | 2000-05-29 | 株式会社東芝 | 超音波診断装置 |
JP3373268B2 (ja) * | 1993-12-10 | 2003-02-04 | ジーイー横河メディカルシステム株式会社 | 超音波診断装置 |
JP3480608B2 (ja) * | 1994-11-07 | 2003-12-22 | 株式会社日立メディコ | 投影像の作成方法およびそのための装置 |
US5608690A (en) * | 1995-03-02 | 1997-03-04 | Acuson Corporation | Transmit beamformer with frequency dependent focus |
JP3865800B2 (ja) * | 1995-05-15 | 2007-01-10 | 株式会社東芝 | 超音波診断装置 |
JP3492035B2 (ja) * | 1995-07-28 | 2004-02-03 | ジーイー横河メディカルシステム株式会社 | 超音波診断装置 |
US5706819A (en) * | 1995-10-10 | 1998-01-13 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with harmonic contrast agents |
US5833613A (en) * | 1996-09-27 | 1998-11-10 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with contrast agents |
EP1374777A1 (en) * | 1995-10-10 | 2004-01-02 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with contrast agents |
US5720291A (en) * | 1996-03-22 | 1998-02-24 | Advanced Technology Laboratories, Inc. | Three dimensional medical ultrasonic diagnostic image of tissue texture and vasculature |
JPH1071146A (ja) * | 1996-08-30 | 1998-03-17 | Ge Yokogawa Medical Syst Ltd | 動作諸元表示方法および装置並びに超音波映像装置 |
JP4382884B2 (ja) * | 1996-11-08 | 2009-12-16 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 高調波による超音波画像処理方法および装置 |
US5928151A (en) * | 1997-08-22 | 1999-07-27 | Acuson Corporation | Ultrasonic system and method for harmonic imaging in three dimensions |
-
1998
- 1998-04-23 US US09/065,213 patent/US6102858A/en not_active Expired - Lifetime
-
1999
- 1999-04-23 DE DE69936418T patent/DE69936418T2/de not_active Expired - Lifetime
- 1999-04-23 JP JP11569099A patent/JP4767377B2/ja not_active Expired - Fee Related
- 1999-04-23 EP EP99303174A patent/EP0952463B1/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
DE69936418D1 (de) | 2007-08-16 |
EP0952463B1 (en) | 2007-07-04 |
EP0952463A3 (en) | 2003-02-05 |
EP0952463A2 (en) | 1999-10-27 |
DE69936418T2 (de) | 2008-03-20 |
JP2000152939A (ja) | 2000-06-06 |
US6102858A (en) | 2000-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4767377B2 (ja) | 物体ボリュームの3次元イメージングのためのシステム及び方法 | |
EP1004896B1 (en) | Three-dimensional ultrasound imaging of velocity and power data | |
EP1046929B1 (en) | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering | |
US5779641A (en) | Method and apparatus for three-dimensional ultrasound imaging by projecting filtered pixel data | |
JP4204095B2 (ja) | 被検体ボリュームの三次元イメージング・システムおよび方法 | |
US5429137A (en) | Acoustic scan conversion method and apparatus for velocity flow | |
EP1008864B1 (en) | Three-dimensional imaging by projecting morphologically filtered pixel data | |
US5954653A (en) | Method and apparatus for automatically enhancing contrast in projected ultrasound image | |
US5904653A (en) | Method and apparatus for three-dimensional ultrasound imaging combining intensity data with color flow velocity or power data | |
US6511426B1 (en) | Medical diagnostic ultrasound system and method for versatile processing | |
US8002704B2 (en) | Method and system for determining contact along a surface of an ultrasound probe | |
US5895358A (en) | Method and apparatus for mapping color flow velocity data into display intensities | |
US20030114756A1 (en) | Method and system for ultrasound blood flow imaging and volume flow calculations | |
US6126603A (en) | Method and apparatus for segmenting color flow mode data using velocity information in three-dimensional ultrasound imaging | |
KR20190044758A (ko) | 3차원 초음파 영상 복원 방법 및 그 초음파 영상 장치 | |
US4252025A (en) | Ultrasonic diagnostic apparatus | |
Hottier et al. | 3D echography: status and perspective | |
US20220280138A1 (en) | Ultrasonic imaging of acoustic attenuation coefficients with elevation compounding | |
EP4208735A1 (en) | Time-aligned plane wave compounding of ultrasound data | |
US6135962A (en) | Method and apparatus for adaptive filtering by counting acoustic sample zeroes in ultrasound imaging | |
JP2001061839A (ja) | 超音波診断装置及び超音波画像処理装置 | |
Bylund | Spatio-temporal filtering of ultrasound image sequences | |
Sato et al. | Analysis of posterior echoes using reconstructed vertical ultrasound images | |
Le Croissette | Signal Processing In Ultrasound | |
Bijnens et al. | Reconstruction of echocardiographic sector images from digital radiofrequency data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060418 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090226 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090317 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20090617 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20090617 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20090617 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090716 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20090707 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090916 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091117 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100217 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100222 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100615 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100915 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100921 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101215 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110125 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110414 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20110426 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110531 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110615 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140624 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |