JPH1128214A - 三次元イメージング・システムおよび方法 - Google Patents
三次元イメージング・システムおよび方法Info
- Publication number
- JPH1128214A JPH1128214A JP10123449A JP12344998A JPH1128214A JP H1128214 A JPH1128214 A JP H1128214A JP 10123449 A JP10123449 A JP 10123449A JP 12344998 A JP12344998 A JP 12344998A JP H1128214 A JPH1128214 A JP H1128214A
- Authority
- JP
- Japan
- Prior art keywords
- velocity data
- data
- image
- velocity
- volume
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03B—APPARATUS OR ARRANGEMENTS FOR TAKING PHOTOGRAPHS OR FOR PROJECTING OR VIEWING THEM; APPARATUS OR ARRANGEMENTS EMPLOYING ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ACCESSORIES THEREFOR
- G03B42/00—Obtaining records using waves other than optical waves; Visualisation of such records by using optical means
- G03B42/06—Obtaining records using waves other than optical waves; Visualisation of such records by using optical means using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- 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/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
-
- 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/52034—Data rate converters
-
- 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52068—Stereoscopic displays; Three-dimensional displays; Pseudo 3D displays
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
- A61B8/543—Control of the diagnostic device involving acquisition triggered by a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Abstract
(57)【要約】
【課題】 動いている超音波散乱体を含む被検体ボリュ
ームの三次元イメージング方法およびシステムを提供す
る。 【解決手段】 被検体ボリューム内の多数のサンプル・
ボリュームから反射された超音波エコーを検出し、超音
波散乱体によって反射された超音波エコーから少なくと
も部分的に導き出された速度データを取得し、多数のサ
ンプル・ボリュームの各々に対して速度データを記憶
し、該記憶されている速度データから、所定の範囲内に
ある速度データのみを含むソース速度データ組を検索
し、ソース速度データ組の中の速度データを第１の像平
面へ投影して、第１の投影像を表す投影速度データ組を
形成し、第１の投影像を表示する。
ームの三次元イメージング方法およびシステムを提供す
る。 【解決手段】 被検体ボリューム内の多数のサンプル・
ボリュームから反射された超音波エコーを検出し、超音
波散乱体によって反射された超音波エコーから少なくと
も部分的に導き出された速度データを取得し、多数のサ
ンプル・ボリュームの各々に対して速度データを記憶
し、該記憶されている速度データから、所定の範囲内に
ある速度データのみを含むソース速度データ組を検索
し、ソース速度データ組の中の速度データを第１の像平
面へ投影して、第１の投影像を表す投影速度データ組を
形成し、第１の投影像を表示する。
Description
【０００１】
【発明の分野】本発明は、一般的には、医学診断のため
の人体の解剖学的構造の超音波イメージングに関するも
のである。特に、本発明は、人体内の動いている流体ま
たは組織から反射された超音波エコーのドップラー偏移
を検出することによって該動いている流体または組織の
三次元イメージングを行う方法および装置に関するもの
である。
の人体の解剖学的構造の超音波イメージングに関するも
のである。特に、本発明は、人体内の動いている流体ま
たは組織から反射された超音波エコーのドップラー偏移
を検出することによって該動いている流体または組織の
三次元イメージングを行う方法および装置に関するもの
である。
【０００２】
【発明の背景】通常の超音波スキャナは、画素の輝度が
エコー信号の強度に基づいて定められた、組織の二次元
Ｂモード像を作成する。カラー流れイメージングでは、
血液の流れまたは組織の動きを映像化することが出来
る。ドップラー効果を利用して、心臓および血管内の血
液の流れを測定することはよく知られている。後方散乱
された超音波の周波数偏移を使用することにより、組織
からの後方散乱体または血液の速度を測定することが出
来る。後方散乱された超音波の周波数の変化すなわち周
波数偏移は、血液がトランスジューサの方へ向かって流
れているときは増加し、また血液がトランスジューサか
ら遠ざかる向きに流れているときは減少する。このドッ
プラー偏移は、流れの速度および方向を表すために異な
るカラーを使用して表示することが出来る。カラー流れ
モードは、数百の隣接したサンプル・ボリュームを同時
に表示し、これらは全て各々のサンプル・ボリュームの
速度を表すためにカラー符号化されている。カラー流れ
像はＢモード像に重畳することが出来る。
エコー信号の強度に基づいて定められた、組織の二次元
Ｂモード像を作成する。カラー流れイメージングでは、
血液の流れまたは組織の動きを映像化することが出来
る。ドップラー効果を利用して、心臓および血管内の血
液の流れを測定することはよく知られている。後方散乱
された超音波の周波数偏移を使用することにより、組織
からの後方散乱体または血液の速度を測定することが出
来る。後方散乱された超音波の周波数の変化すなわち周
波数偏移は、血液がトランスジューサの方へ向かって流
れているときは増加し、また血液がトランスジューサか
ら遠ざかる向きに流れているときは減少する。このドッ
プラー偏移は、流れの速度および方向を表すために異な
るカラーを使用して表示することが出来る。カラー流れ
モードは、数百の隣接したサンプル・ボリュームを同時
に表示し、これらは全て各々のサンプル・ボリュームの
速度を表すためにカラー符号化されている。カラー流れ
像はＢモード像に重畳することが出来る。
【０００３】本発明は、４つの主要なサブシステム（図
１参照）、すなわちビーム形成装置２、処理装置４、走
査変換器／表示制御器６および主制御器８で構成されて
いる超音波イメージング・システムに用いられる。シス
テムの制御は主制御器８に集中しており、主制御器はオ
ペレータ・インターフェース（図示していない）を介し
てオペレータ入力を受け入れて、種々のサブシステムを
制御する。主制御器はまたシステム・タイミングおよび
制御信号を発生し、これらの信号がシステム制御母線１
０および走査制御母線（図示していない）を介して分配
される。
１参照）、すなわちビーム形成装置２、処理装置４、走
査変換器／表示制御器６および主制御器８で構成されて
いる超音波イメージング・システムに用いられる。シス
テムの制御は主制御器８に集中しており、主制御器はオ
ペレータ・インターフェース（図示していない）を介し
てオペレータ入力を受け入れて、種々のサブシステムを
制御する。主制御器はまたシステム・タイミングおよび
制御信号を発生し、これらの信号がシステム制御母線１
０および走査制御母線（図示していない）を介して分配
される。
【０００４】主データ路が、トランスジューサからビー
ム形成装置へのディジタル化されたＲＦ入力で始まる。
ビーム形成装置は、２つの加算されたディジタル・ベー
スバンド受信ビームを出力する。ベースバンド・データ
はＢモード処理装置４Ａおよびカラー流れ処理装置４Ｂ
に入力され、そこで取得モードに従って処理されて、走
査変換器／表示処理装置６へ処理済み音響ベクトル（ビ
ーム）データとして出力される。走査変換器／表示処理
装置６がこの処理済み音響データを受け取って、ラスタ
走査フォーマットの像に対するビデオ表示信号をカラー
表示モニタ１２へ出力する。走査変換器／表示処理装置
６は、主制御器８と協働して、表示用の多数の像、表示
の注釈、グラフィック・オーバーレイ、並びにシネ（ｃ
ｉｎｅ；動画）ループおよび記録された時間線データの
リプレイ（再生）をフォーマッティングする。
ム形成装置へのディジタル化されたＲＦ入力で始まる。
ビーム形成装置は、２つの加算されたディジタル・ベー
スバンド受信ビームを出力する。ベースバンド・データ
はＢモード処理装置４Ａおよびカラー流れ処理装置４Ｂ
に入力され、そこで取得モードに従って処理されて、走
査変換器／表示処理装置６へ処理済み音響ベクトル（ビ
ーム）データとして出力される。走査変換器／表示処理
装置６がこの処理済み音響データを受け取って、ラスタ
走査フォーマットの像に対するビデオ表示信号をカラー
表示モニタ１２へ出力する。走査変換器／表示処理装置
６は、主制御器８と協働して、表示用の多数の像、表示
の注釈、グラフィック・オーバーレイ、並びにシネ（ｃ
ｉｎｅ；動画）ループおよび記録された時間線データの
リプレイ（再生）をフォーマッティングする。
【０００５】Ｂモード処理装置４Ａは、ビーム形成装置
からのベースバンド・データを対数圧縮した信号包絡線
へ変換する。Ｂ機能は、信号の包絡線の時間変化振幅
を、各画素に対して８ビットの出力を使用してグレース
ケールで映像化する。ベースバンド信号の包絡線は、ベ
ースバンド・データが表すベクトルの大きさである。血
管や心室などの内部から反射された音波の周波数は血球
の速度に比例してシフトすなわち偏移する。すなわち、
血球がトランスジューサへ向かって動いている場合は正
に偏移し、また血球がトランスジューサから離れる向き
に動いている場合は負に偏移する。カラー流れ（ＣＦ）
処理装置４Ｂは、イメージング平面内における血液の速
度の実時間二次元像を作成するために使用される。血液
の速度は、特定のレンジ・ゲートにおいてファイアリン
グ（ｆｉｒｉｎｇ）相互の間での位相シフト（移相）を
測定することによって計算される。像内の１つのレンジ
・ゲートでドップラー・スペクトルを測定する代わり
に、各々のベクトルに沿った多数のレンジ・ゲートおよ
び多数のベクトル位置から平均血液速度が計算され、こ
の情報から二次元像が作成される。カラー流れ処理装置
の構造および作用が、ここに引用する米国特許第５，５
２４，６２９号に記載されている。
からのベースバンド・データを対数圧縮した信号包絡線
へ変換する。Ｂ機能は、信号の包絡線の時間変化振幅
を、各画素に対して８ビットの出力を使用してグレース
ケールで映像化する。ベースバンド信号の包絡線は、ベ
ースバンド・データが表すベクトルの大きさである。血
管や心室などの内部から反射された音波の周波数は血球
の速度に比例してシフトすなわち偏移する。すなわち、
血球がトランスジューサへ向かって動いている場合は正
に偏移し、また血球がトランスジューサから離れる向き
に動いている場合は負に偏移する。カラー流れ（ＣＦ）
処理装置４Ｂは、イメージング平面内における血液の速
度の実時間二次元像を作成するために使用される。血液
の速度は、特定のレンジ・ゲートにおいてファイアリン
グ（ｆｉｒｉｎｇ）相互の間での位相シフト（移相）を
測定することによって計算される。像内の１つのレンジ
・ゲートでドップラー・スペクトルを測定する代わり
に、各々のベクトルに沿った多数のレンジ・ゲートおよ
び多数のベクトル位置から平均血液速度が計算され、こ
の情報から二次元像が作成される。カラー流れ処理装置
の構造および作用が、ここに引用する米国特許第５，５
２４，６２９号に記載されている。
【０００６】カラー流れ処理装置は、（８ビットの）速
度信号、（４ビットの）分散（乱れ）信号および（８ビ
ットの）パワー信号を発生する。オペレータが、速度お
よび分散またはパワーを走査変換器へ出力するかどうか
選択する。出力信号が、ビデオ処理装置２２内に含まれ
ているクロミナンス制御ルックアップ・テーブルに入力
される。ルックアップ・テーブル内の各々のアドレスは
２４ビットを記憶する。作成される像中の各々の画素に
対して、８ビットが赤の強度を制御し、８ビットが緑の
強度を制御し、８ビットが青の強度を制御する。これら
のビット・パターンは、流れの速度の方向または大きさ
が変わったときに各々の位置の画素の色が変わるように
予め選定されている。例えば、トランスジューサへ向か
う流れは赤で示され、トランスジューサから離れる向き
の流れは青で示される。流れが速くなると、色はより明
るくなる。
度信号、（４ビットの）分散（乱れ）信号および（８ビ
ットの）パワー信号を発生する。オペレータが、速度お
よび分散またはパワーを走査変換器へ出力するかどうか
選択する。出力信号が、ビデオ処理装置２２内に含まれ
ているクロミナンス制御ルックアップ・テーブルに入力
される。ルックアップ・テーブル内の各々のアドレスは
２４ビットを記憶する。作成される像中の各々の画素に
対して、８ビットが赤の強度を制御し、８ビットが緑の
強度を制御し、８ビットが青の強度を制御する。これら
のビット・パターンは、流れの速度の方向または大きさ
が変わったときに各々の位置の画素の色が変わるように
予め選定されている。例えば、トランスジューサへ向か
う流れは赤で示され、トランスジューサから離れる向き
の流れは青で示される。流れが速くなると、色はより明
るくなる。
【０００７】走査変換器／表示制御器６のＢモード音響
線メモリ１４Ａおよびカラー音響線メモリ１４Ｂが、処
理装置４Ａおよび４Ｂからの処理済みディジタル・デー
タをそれぞれ受け取って、Ｂモード・データおよびカラ
ー流れデータを、極座標（Ｒ−θ）セクター・フォーマ
ットまたはデカルト座標線形配列から、適切にスケーリ
ングしたデカルト座標表示画素データへ座標変換する。
この画素データはＸ−Ｙ表示メモリ１８に記憶される。
Ｂモードでは、強度データがＸ−Ｙ表示メモリ１８に記
憶され、各々のアドレスに３つの８ビット画素が記憶さ
れる。また、カラー流れモードでは、データはメモリに
次のように記憶される。すなわち強度データ（８ビッ
ト）、速度またはパワー・データ（８ビット）および分
散（乱れ）データ（４ビット）が記憶される。
線メモリ１４Ａおよびカラー音響線メモリ１４Ｂが、処
理装置４Ａおよび４Ｂからの処理済みディジタル・デー
タをそれぞれ受け取って、Ｂモード・データおよびカラ
ー流れデータを、極座標（Ｒ−θ）セクター・フォーマ
ットまたはデカルト座標線形配列から、適切にスケーリ
ングしたデカルト座標表示画素データへ座標変換する。
この画素データはＸ−Ｙ表示メモリ１８に記憶される。
Ｂモードでは、強度データがＸ−Ｙ表示メモリ１８に記
憶され、各々のアドレスに３つの８ビット画素が記憶さ
れる。また、カラー流れモードでは、データはメモリに
次のように記憶される。すなわち強度データ（８ビッ
ト）、速度またはパワー・データ（８ビット）および分
散（乱れ）データ（４ビット）が記憶される。
【０００８】カラー流れまたはＢモード・データの多数
の相次ぐフレームが、先入れ先出し形式でシネ・メモリ
（ｃｉｎｅ ｍｅｍｏｒｙ）２４に記憶される。シネ・
メモリは、背後で動作していて、使用者に実時間で表示
される像データを捕獲するサーキュラー像バッファに類
似するものである。使用者がシステムを停止（ｆｒｅｅ
ｚｅ）したとき、使用者はシネ・メモリに前に捕獲され
た像データを見ることが出来る。表示された像上にグラ
フィック・オバーレイを作成するためのグラフィック・
データが、時間線／グラフィック処理装置及び表示メモ
リ２０において作成されて記憶される。ビデオ処理装置
２２が、グラフィック・データと像データと時間線デー
タとの間でマルチプレクシングを行って、ビデオ表示モ
ニタ１２上にラスタ走査形式で表示させる最終的なビデ
オ出力を作成する。更に、ビデオ処理装置は様々なグレ
ースケールおよびカラー・マップを提供すると共に、グ
レースケールおよびカラー像を組み合わせる。
の相次ぐフレームが、先入れ先出し形式でシネ・メモリ
（ｃｉｎｅ ｍｅｍｏｒｙ）２４に記憶される。シネ・
メモリは、背後で動作していて、使用者に実時間で表示
される像データを捕獲するサーキュラー像バッファに類
似するものである。使用者がシステムを停止（ｆｒｅｅ
ｚｅ）したとき、使用者はシネ・メモリに前に捕獲され
た像データを見ることが出来る。表示された像上にグラ
フィック・オバーレイを作成するためのグラフィック・
データが、時間線／グラフィック処理装置及び表示メモ
リ２０において作成されて記憶される。ビデオ処理装置
２２が、グラフィック・データと像データと時間線デー
タとの間でマルチプレクシングを行って、ビデオ表示モ
ニタ１２上にラスタ走査形式で表示させる最終的なビデ
オ出力を作成する。更に、ビデオ処理装置は様々なグレ
ースケールおよびカラー・マップを提供すると共に、グ
レースケールおよびカラー像を組み合わせる。
【０００９】通常の超音波イメージング・システムは、
Ｂモードまたはカラー流れモードの像をシネ・メモリ２
４に連続して収集する。シネ・メモリ２４は、単一像検
討および多重像ループ検討のための常駐ディジタル像記
憶並びに様々な制御機能を行う。単一像シネ・リプレイ
の際に表示される関心のある領域は、像の取得の際に使
用されたものである。シネ・メモリはまた、主制御器８
を介してディジタル大容量記憶装置へ像を転送するため
のバッファとして作用する。
Ｂモードまたはカラー流れモードの像をシネ・メモリ２
４に連続して収集する。シネ・メモリ２４は、単一像検
討および多重像ループ検討のための常駐ディジタル像記
憶並びに様々な制御機能を行う。単一像シネ・リプレイ
の際に表示される関心のある領域は、像の取得の際に使
用されたものである。シネ・メモリはまた、主制御器８
を介してディジタル大容量記憶装置へ像を転送するため
のバッファとして作用する。
【００１０】二次元の超音波像は、観察者が走査してい
る解剖学的構造を具体的に思い浮かべることが出来ない
ので、解釈するのがしばしば困難である。しかしなが
ら、超音波プローブが関心のある領域にわたって掃引さ
れて、二次元像を累積して三次元像を形成できれば、熟
練した観察者および未熟な観察者が共に解剖学的構造を
容易に思い浮かべることが出来る。特に、動いている流
体または組織の三次元超音波イメージングは有利であろ
う。
る解剖学的構造を具体的に思い浮かべることが出来ない
ので、解釈するのがしばしば困難である。しかしなが
ら、超音波プローブが関心のある領域にわたって掃引さ
れて、二次元像を累積して三次元像を形成できれば、熟
練した観察者および未熟な観察者が共に解剖学的構造を
容易に思い浮かべることが出来る。特に、動いている流
体または組織の三次元超音波イメージングは有利であろ
う。
【００１１】
【発明の概要】本発明は、速度のドップラー偏移データ
を使用して、動いている流体または組織の三次元像を表
示するための方法および装置である。イメージング平面
に投影される画素強度データのボリュームを制限するこ
とによって、像のセグメント化が改善される。本発明に
よれば、関心のあるボリュームが動いている流体または
組織を含んでいる場合、動いている流体または組織から
反射された超音波中に存在するドップラー偏移を検出
し、このドップラー偏移を使用することにより、投影さ
れるべき画素データの量を制限することが出来る。
を使用して、動いている流体または組織の三次元像を表
示するための方法および装置である。イメージング平面
に投影される画素強度データのボリュームを制限するこ
とによって、像のセグメント化が改善される。本発明に
よれば、関心のあるボリュームが動いている流体または
組織を含んでいる場合、動いている流体または組織から
反射された超音波中に存在するドップラー偏移を検出
し、このドップラー偏移を使用することにより、投影さ
れるべき画素データの量を制限することが出来る。
【００１２】本発明を用いる装置は、患者の心臓サイク
ル波形内の所定の特徴（例えば、ピーク）の発生のよう
な周期的な外部のトリガ事象に応答して、カラー流れモ
ードの像をシネ・メモリ内に収集する超音波スキャナを
有する。本発明によれば、各々の取得された画素は、各
々のサンプル・ボリュームから反射されたエコー信号中
の検出されたドップラー偏移から導き出した８ビットの
速度値を有する。この速度情報は、主制御器（または別
の専用の処理装置）によってイメージング平面に投影さ
れるべき速度データを識別するために処理される。速度
データは、所定の範囲内にある大きさを持つ各々の速度
データ、例えば非ゼロの大きさを持つ速度データを検索
することによって処理される。本発明の広義の概念によ
れば、最小閾値よりも大きい速度値、又は最大閾値より
も小さい速度値、或いは最小閾値と最大閾値との間の速
度値を、シネ・メモリから投影のために検索していよ
い。この範疇内の画素は、投影像を再構成する際に使用
するためのソース画素データ・ボリュームを画成するた
めに用いられる。ソース画素データ・ボリューム内に含
まれている画素は、種々のイメージング平面への投影を
再構成する際に使用される。ソース・データ・ボリュー
ムの外側の画素は、再構成に使用されない。その結果、
セグメント化が改善された投影速度像が得られる。
ル波形内の所定の特徴（例えば、ピーク）の発生のよう
な周期的な外部のトリガ事象に応答して、カラー流れモ
ードの像をシネ・メモリ内に収集する超音波スキャナを
有する。本発明によれば、各々の取得された画素は、各
々のサンプル・ボリュームから反射されたエコー信号中
の検出されたドップラー偏移から導き出した８ビットの
速度値を有する。この速度情報は、主制御器（または別
の専用の処理装置）によってイメージング平面に投影さ
れるべき速度データを識別するために処理される。速度
データは、所定の範囲内にある大きさを持つ各々の速度
データ、例えば非ゼロの大きさを持つ速度データを検索
することによって処理される。本発明の広義の概念によ
れば、最小閾値よりも大きい速度値、又は最大閾値より
も小さい速度値、或いは最小閾値と最大閾値との間の速
度値を、シネ・メモリから投影のために検索していよ
い。この範疇内の画素は、投影像を再構成する際に使用
するためのソース画素データ・ボリュームを画成するた
めに用いられる。ソース画素データ・ボリューム内に含
まれている画素は、種々のイメージング平面への投影を
再構成する際に使用される。ソース・データ・ボリュー
ムの外側の画素は、再構成に使用されない。その結果、
セグメント化が改善された投影速度像が得られる。
【００１３】動いている流体または組織を三次元で映像
化するため、主制御器が、レイ・キャスティング（ｒａ
ｙ−ｃａｓｔｉｎｇ）法を使用して、ソース・データ・
ボリューム内の画素速度データを、複数の回転された像
平面に投影するアルゴリズムを実行する。各々の投影か
ら生じる投影データはシネ・メモリに戻される。投影デ
ータは、表示像の中心領域に表示するために記憶され
る。表示像の周辺部は、オプションとして、システム・
オペレータによるシステムの停止（ｆｒｅｅｚｉｎｇ）
の前にＸ−Ｙメモリから読み出された最後のフレームか
らの未投影の速度データを含んでいてよい。この代わり
に、投影速度像は、背景フレームからの未投影のデータ
なしに、表示することが出来る。
化するため、主制御器が、レイ・キャスティング（ｒａ
ｙ−ｃａｓｔｉｎｇ）法を使用して、ソース・データ・
ボリューム内の画素速度データを、複数の回転された像
平面に投影するアルゴリズムを実行する。各々の投影か
ら生じる投影データはシネ・メモリに戻される。投影デ
ータは、表示像の中心領域に表示するために記憶され
る。表示像の周辺部は、オプションとして、システム・
オペレータによるシステムの停止（ｆｒｅｅｚｉｎｇ）
の前にＸ−Ｙメモリから読み出された最後のフレームか
らの未投影の速度データを含んでいてよい。この代わり
に、投影速度像は、背景フレームからの未投影のデータ
なしに、表示することが出来る。
【００１４】本発明の方法によれば、主制御器がシネ・
メモリから速度データを検索する。速度データは、患者
の心臓サイクル波形内の所定の特徴（例えば、ピーク）
の発生に応答して、シネ・メモリに記憶される。１つの
取得モードでは、毎サイクル当り１フレームがＸ−Ｙメ
モリからシネ・メモリへ転送される。次いで、関心のあ
るボリュームからの速度データが、前に述べた閾値を用
いる方式に従って処理されて、種々のイメージング平面
に投影されるべき速度のソース・データ・ボリュームを
取得する。別の取得モードでは、多数の相次ぐフレーム
がトリガ事象に応答して音響フレーム速度でシネ・メモ
リに記憶される。次いで、速度データを前に述べたよう
に処理して、投影像を再構成し、この投影像をシネ・メ
モリに戻す。シネ・メモリに記憶されたこれらの再構成
されたフレームは、次いでシステム・オペレータによっ
て選択的に表示することが出来る。
メモリから速度データを検索する。速度データは、患者
の心臓サイクル波形内の所定の特徴（例えば、ピーク）
の発生に応答して、シネ・メモリに記憶される。１つの
取得モードでは、毎サイクル当り１フレームがＸ−Ｙメ
モリからシネ・メモリへ転送される。次いで、関心のあ
るボリュームからの速度データが、前に述べた閾値を用
いる方式に従って処理されて、種々のイメージング平面
に投影されるべき速度のソース・データ・ボリュームを
取得する。別の取得モードでは、多数の相次ぐフレーム
がトリガ事象に応答して音響フレーム速度でシネ・メモ
リに記憶される。次いで、速度データを前に述べたよう
に処理して、投影像を再構成し、この投影像をシネ・メ
モリに戻す。シネ・メモリに記憶されたこれらの再構成
されたフレームは、次いでシステム・オペレータによっ
て選択的に表示することが出来る。
【００１５】
【好ましい実施態様の説明】図２を参照して説明する
と、主制御器は中央処理装置（ＣＰＵ）４２およびラン
ダム・アクセス・メモリ４４を有する。ＣＰＵ４２はそ
の中に読出し専用メモリ（ＲＯＭ）を含み、該メモリ
は、取得した速度および乱れデータのボリュームを、異
なる角度で取った多数の三次元投影像に変換するのに使
用されるルーチンを記憶している。ＣＰＵ４２は、シス
テム制御母線１０を介してＸ−Ｙ表示メモリ１８および
シネ・メモリ４２を制御する。具体的に述べると、ＣＰ
Ｕ４２はＸ−Ｙ表示メモリ１８からビデオ処理装置２２
およびシネ・メモリ２４へのデータの流れを制御すると
共に、シネ・メモリ２４からビデオ処理装置２２および
ＣＰＵ４２自身へのデータの流れを制御する。超音波イ
メージング・システムがカラー流れモードで動作してい
るとき、被検体の多数の平行な走査またはスライスの内
の１つを表す各フレームのカラー流れデータが、Ｘ−Ｙ
表示メモリ１８に記憶されて、次のサイクルでビデオ処
理装置２２およびシネ・メモリ２４へ伝送される。走査
された被検体ボリュームを表す１スタック（ｓｔａｃ
ｋ）のフレームが、シネ・メモリ２４内の一区分２４Ａ
に記憶される。初期化の際（図３のステップ２６を参
照）、ＣＰＵ４２はシネ・メモリの区分２４Ａから関心
のある被検体ボリュームに対応するカラー流れデータの
みを検索する。これは、各々の選択されたフレームか
ら、関心のある領域内のカラー流れデータのみを検索す
ることにより達成される。多数の選択されたフレームの
各々からの関心のある領域に対応するカラー流れデータ
は、関心のあるソース・データ・ボリュームを形成す
る。
と、主制御器は中央処理装置（ＣＰＵ）４２およびラン
ダム・アクセス・メモリ４４を有する。ＣＰＵ４２はそ
の中に読出し専用メモリ（ＲＯＭ）を含み、該メモリ
は、取得した速度および乱れデータのボリュームを、異
なる角度で取った多数の三次元投影像に変換するのに使
用されるルーチンを記憶している。ＣＰＵ４２は、シス
テム制御母線１０を介してＸ−Ｙ表示メモリ１８および
シネ・メモリ４２を制御する。具体的に述べると、ＣＰ
Ｕ４２はＸ−Ｙ表示メモリ１８からビデオ処理装置２２
およびシネ・メモリ２４へのデータの流れを制御すると
共に、シネ・メモリ２４からビデオ処理装置２２および
ＣＰＵ４２自身へのデータの流れを制御する。超音波イ
メージング・システムがカラー流れモードで動作してい
るとき、被検体の多数の平行な走査またはスライスの内
の１つを表す各フレームのカラー流れデータが、Ｘ−Ｙ
表示メモリ１８に記憶されて、次のサイクルでビデオ処
理装置２２およびシネ・メモリ２４へ伝送される。走査
された被検体ボリュームを表す１スタック（ｓｔａｃ
ｋ）のフレームが、シネ・メモリ２４内の一区分２４Ａ
に記憶される。初期化の際（図３のステップ２６を参
照）、ＣＰＵ４２はシネ・メモリの区分２４Ａから関心
のある被検体ボリュームに対応するカラー流れデータの
みを検索する。これは、各々の選択されたフレームか
ら、関心のある領域内のカラー流れデータのみを検索す
ることにより達成される。多数の選択されたフレームの
各々からの関心のある領域に対応するカラー流れデータ
は、関心のあるソース・データ・ボリュームを形成す
る。
【００１６】本発明の好ましい実施態様によれば、関心
のあるソース・データ・ボリュームは、所定の範囲内の
速度成分、例えば非ゼロの（すなわち、ゼロでない）速
度値を持つ画素を有する。ソース・データ・ボリューム
内の速度データは、異なる視角で取った投影像を再構成
するために使用される。速度の投影が、米国特許第５，
２２６，１１３号明細書に開示されているレイ・キャス
ティング・アルゴリズムを使用して一連の変換を実行す
るＣＰＵ４２で、再構成される。相次ぐ変換は、所定の
角度範囲内（例えば、＋９０°乃至−９０°の範囲内）
で、所定の角度増分（例えば、１０°の間隔）で作られ
た最大の又は最小の、或いは平均化された速度投影を表
す。しかしながら、角度増分は１０°である必要はな
く、また本発明が特定の角度範囲に制限されるものでも
ない。
のあるソース・データ・ボリュームは、所定の範囲内の
速度成分、例えば非ゼロの（すなわち、ゼロでない）速
度値を持つ画素を有する。ソース・データ・ボリューム
内の速度データは、異なる視角で取った投影像を再構成
するために使用される。速度の投影が、米国特許第５，
２２６，１１３号明細書に開示されているレイ・キャス
ティング・アルゴリズムを使用して一連の変換を実行す
るＣＰＵ４２で、再構成される。相次ぐ変換は、所定の
角度範囲内（例えば、＋９０°乃至−９０°の範囲内）
で、所定の角度増分（例えば、１０°の間隔）で作られ
た最大の又は最小の、或いは平均化された速度投影を表
す。しかしながら、角度増分は１０°である必要はな
く、また本発明が特定の角度範囲に制限されるものでも
ない。
【００１７】本発明で使用されるレイ・キャスティング
法によれば、サンプル５０（図４参照）の立体的に表現
された投影像が、超音波トランスジューサ・アレイを使
用して被検体ボリュームを走査することによって、任意
の視角から、例えば角度パラメータ（θ、φ）で表記さ
れる球面投影角から表示される。ここで、θは視線５８
の射影線（ｅｘｔｅｎｓｉｏｎ）５８’がＸ−Ｙ平面上
に作る角度であり、φは視線５８が射影線５８’に対し
て作る角度である。サンプル・ボリューム５２は、一連
の積み重なった隣接したスライスまたはシートＯＳ1 、
ＯＳ2 、・・・、ＯＳk を作成するように走査され、各
々のスライスは同じ数の被検体ボリューム要素（ボクセ
ル）ＯＶを含む。各々のボクセルはシートの平面（例え
ば、Ｘ−Ｙ平面）内に矩形の輪郭を持つ。この輪郭が正
方形になるように相補的な辺は等しい長さＳであってよ
いが、シートの厚さＴは一般にいずれの辺の長さにも等
しくない。従って、第１の被検体スライスＯＳ1 が第１
の多数の被検体ボクセルＯＶi,j,1 を含み、ここでｉお
よびｊはボクセルのそれぞれのＸ軸位置およびＹ軸位置
である。同様に、第２の被検体スライスＯＳ2 が第２の
多数の被検体ボクセルＯＶi,j,2 を含む。任意の被検体
スライスＯＳk が多数の被検体ボクセルＯＶ i,j,k を含
む。ここで、ｋはそのボクセルのＺ軸位置である。
法によれば、サンプル５０（図４参照）の立体的に表現
された投影像が、超音波トランスジューサ・アレイを使
用して被検体ボリュームを走査することによって、任意
の視角から、例えば角度パラメータ（θ、φ）で表記さ
れる球面投影角から表示される。ここで、θは視線５８
の射影線（ｅｘｔｅｎｓｉｏｎ）５８’がＸ−Ｙ平面上
に作る角度であり、φは視線５８が射影線５８’に対し
て作る角度である。サンプル・ボリューム５２は、一連
の積み重なった隣接したスライスまたはシートＯＳ1 、
ＯＳ2 、・・・、ＯＳk を作成するように走査され、各
々のスライスは同じ数の被検体ボリューム要素（ボクセ
ル）ＯＶを含む。各々のボクセルはシートの平面（例え
ば、Ｘ−Ｙ平面）内に矩形の輪郭を持つ。この輪郭が正
方形になるように相補的な辺は等しい長さＳであってよ
いが、シートの厚さＴは一般にいずれの辺の長さにも等
しくない。従って、第１の被検体スライスＯＳ1 が第１
の多数の被検体ボクセルＯＶi,j,1 を含み、ここでｉお
よびｊはボクセルのそれぞれのＸ軸位置およびＹ軸位置
である。同様に、第２の被検体スライスＯＳ2 が第２の
多数の被検体ボクセルＯＶi,j,2 を含む。任意の被検体
スライスＯＳk が多数の被検体ボクセルＯＶ i,j,k を含
む。ここで、ｋはそのボクセルのＺ軸位置である。
【００１８】各々の被検体ボクセルＯＶi,j,k が分析さ
れて、そのデータ値（強度、速度またはパワー）がデー
タ・ボリューム５４の対応するデータ・ボクセルＤＶ
i,j,kに置かれる。データ・ボリュームＤＶi,j,k は、
各々の被検体スライスＯＳk の厚さおよび各々の被検体
ボクセルの面寸法（Ｘ−Ｙ平面におけるボクセルの大き
さ）が一般に同じでなくても、簡単な立方ｉ，ｊ，ｋ格
子である。すなわち、被検体ボリュームは各々のボクセ
ルに対して異なるＸ、ＹおよびＺ寸法を持っていてよい
ばかりではなく、任意の次元におけるボクセルの総数が
同じである必要もない。例えば、典型的な超音波三次元
走査では、各々のスライスが２５６×２５６行列のボク
セルを持つ、１２８個のスライスを得ることが出来る。
れて、そのデータ値（強度、速度またはパワー）がデー
タ・ボリューム５４の対応するデータ・ボクセルＤＶ
i,j,kに置かれる。データ・ボリュームＤＶi,j,k は、
各々の被検体スライスＯＳk の厚さおよび各々の被検体
ボクセルの面寸法（Ｘ−Ｙ平面におけるボクセルの大き
さ）が一般に同じでなくても、簡単な立方ｉ，ｊ，ｋ格
子である。すなわち、被検体ボリュームは各々のボクセ
ルに対して異なるＸ、ＹおよびＺ寸法を持っていてよい
ばかりではなく、任意の次元におけるボクセルの総数が
同じである必要もない。例えば、典型的な超音波三次元
走査では、各々のスライスが２５６×２５６行列のボク
セルを持つ、１２８個のスライスを得ることが出来る。
【００１９】ＣＰＵ４２によって用いられる公知の技術
に従って、各々のデータ・ボクセルＤＶi,j,k 内の格子
点から像平面５６へのレイ・キャスティングすなわち射
線（ｒａｙ）の投射によって被検体５０の像が投影され
る（図３のステップ３４）。便宜のため、格子点は、例
えばデータ・ボリュームの原点に最も近いデータ・ボク
セルの頂点であってよい。投射された射線６２は被検体
ボリューム５２をみる球面角度パラメータ（θ，φ）か
ら変換された球面角度パラメータ（α，β）を持つ投影
角でデータ・ボリューム５４を出て行く。これらの２つ
の角度は、非立方体の被検体ボリューム５２に対して立
方体のデータ・ボリューム５４を使用したことによる幾
何学的歪みにより、同じではない。しかし、投射された
射線６２はバーＸ−バーＹ平面の射影線６２’を持ち
（ここで、バーＸおよびバーＹは、図示のようにＸおよ
びＹのそれぞれの頭に横棒を引いた記号を表す）、射影
線６２’はデータ・ボリュームのバーＸ軸に対して角度
αを作り、また射線６２はＺ軸と角度βを作る。角度α
およびβは回転プロセス（以下に説明する）によって決
定されて、（球面座標における操作を仮定すると）所望
の視角（θ，φ）でサンプル・ボリューム５２を見るこ
とに対応する。各々の射線６２はデータ・ボリュームの
ボクセル格子点から像平面へ向けて投射される。
に従って、各々のデータ・ボクセルＤＶi,j,k 内の格子
点から像平面５６へのレイ・キャスティングすなわち射
線（ｒａｙ）の投射によって被検体５０の像が投影され
る（図３のステップ３４）。便宜のため、格子点は、例
えばデータ・ボリュームの原点に最も近いデータ・ボク
セルの頂点であってよい。投射された射線６２は被検体
ボリューム５２をみる球面角度パラメータ（θ，φ）か
ら変換された球面角度パラメータ（α，β）を持つ投影
角でデータ・ボリューム５４を出て行く。これらの２つ
の角度は、非立方体の被検体ボリューム５２に対して立
方体のデータ・ボリューム５４を使用したことによる幾
何学的歪みにより、同じではない。しかし、投射された
射線６２はバーＸ−バーＹ平面の射影線６２’を持ち
（ここで、バーＸおよびバーＹは、図示のようにＸおよ
びＹのそれぞれの頭に横棒を引いた記号を表す）、射影
線６２’はデータ・ボリュームのバーＸ軸に対して角度
αを作り、また射線６２はＺ軸と角度βを作る。角度α
およびβは回転プロセス（以下に説明する）によって決
定されて、（球面座標における操作を仮定すると）所望
の視角（θ，φ）でサンプル・ボリューム５２を見るこ
とに対応する。各々の射線６２はデータ・ボリュームの
ボクセル格子点から像平面へ向けて投射される。
【００２０】全ての射線６２は像平面のある部分に突き
当たるが、考慮中の像平面画素６０ａ内に入る射線のみ
が該像平面画素に対するデータに寄与することが出来
る。従って、被検体ボリューム５２の一部分を選び且つ
この選ばれた被検体ボリュームを見る視角（θ，φ）を
選択すると、データ・ボリュームの対応する部分の各ボ
クセル内のデータ値が、像平面５６へある角度（α，
β）で投影される（被検体ボリュームに対して歪んだデ
ータ・ボリュームを見ることに対応する）。従って、第
１のボクセル（例えば、ＤＶi,1,k ）内のデータ値が、
選ばれた角度θおよびφに従って、射線６２ａに沿って
逆投影される。この射線６２ａは画素６０ａ内の衝突位
置６４で像平面５６に突き当たる。これはこの画素に突
き当たる最初の射線であるので、入射データの強度、速
度またはパワー値が所望の画素６０ａに帰する（記憶さ
れる）。データ・ボリューム内の次の第２のボクセル
（例えば、ＤＶi,2,k ）はそのボクセルの格子点から同
じ角度（α，β）で投射される射線６２ｂが関係し、そ
れは像平面５６に衝突位置６４ｂで突き当たる。衝突位
置６４ｂが所望の画素６０ａ内にあると仮定すると、第
２の投影された値が現在記憶されている第１の値と（最
大画素投影のために）比較され、そのうちの大きい方の
値が画素６０ａに対して記憶される。ここで、平均値投
影の場合には、現在の投影されているデータ・ボクセル
の値が、その投影の射線の突き当たる像平面の画素に既
に記憶されている和の値に加算され、次いでその結果の
和が最終的にその画素に突き当たるこのような射線の計
数値で割算されることが理解されよう。選択されたデー
タ・ボリューム内の各々のボクセルが逐次的にエントリ
されて像平面５６に投影されるとき、あるデータ・ボリ
ューム・ボクセル（例えば、ＤＶi,3,k ）はその関連す
る射線６２ｐに沿って投影されるが、所望の画素６０ａ
内に突き当たらず、従ってそのデータ値（例えば、強
度）が画素６０ａに対して現在記憶されているデータ値
と比較されない。特定の三次元の視角（θ，φ）におけ
るデータのその投影に対して、画素画素６０ａに対する
最大データ値がそのとき確立される。しかし、射線６２
ｐが事実上、別の像平面画素（例えば、画素画素６０
ｂ）内にある衝突位置６４ｐを持ち、その画素に記憶さ
れているデータ値と比較されて、比較後の大きい方の値
がその画素に対する記憶装置に戻される。全てのデータ
値は、新しい投影が取られるときにゼロにリセットされ
る。従って、像平面画素の各々は像投影手順の開始時に
リセットされ、（選択された被検体ボリューム５２の部
分によって設定されるような、全空間または選択された
部分内の）データ・ボリューム・ボクセルの全ては個別
に且つ逐次的に走査される。各々データ・ボクセルＤＶ
内のデータ値は、その１つの画素６０内で像平面５６に
突き当たる関連の射線６２により投影される。各々の画
素内の最大値は射線投射されたデータ・ボリューム・ボ
クセルの現在値との間で比較されて、その内の大きい方
の値が決定される。この大きい方の値は、次いで、最大
値像の一部分として記憶される。実際には、最大画素投
影の場合、新しく投射されたデータ・ボクセル値が、新
しく投射された射線が突き当たる像平面画素に対して既
に記憶されているデータ値よりも大きいときだけ、記憶
された最大値が変更される。
当たるが、考慮中の像平面画素６０ａ内に入る射線のみ
が該像平面画素に対するデータに寄与することが出来
る。従って、被検体ボリューム５２の一部分を選び且つ
この選ばれた被検体ボリュームを見る視角（θ，φ）を
選択すると、データ・ボリュームの対応する部分の各ボ
クセル内のデータ値が、像平面５６へある角度（α，
β）で投影される（被検体ボリュームに対して歪んだデ
ータ・ボリュームを見ることに対応する）。従って、第
１のボクセル（例えば、ＤＶi,1,k ）内のデータ値が、
選ばれた角度θおよびφに従って、射線６２ａに沿って
逆投影される。この射線６２ａは画素６０ａ内の衝突位
置６４で像平面５６に突き当たる。これはこの画素に突
き当たる最初の射線であるので、入射データの強度、速
度またはパワー値が所望の画素６０ａに帰する（記憶さ
れる）。データ・ボリューム内の次の第２のボクセル
（例えば、ＤＶi,2,k ）はそのボクセルの格子点から同
じ角度（α，β）で投射される射線６２ｂが関係し、そ
れは像平面５６に衝突位置６４ｂで突き当たる。衝突位
置６４ｂが所望の画素６０ａ内にあると仮定すると、第
２の投影された値が現在記憶されている第１の値と（最
大画素投影のために）比較され、そのうちの大きい方の
値が画素６０ａに対して記憶される。ここで、平均値投
影の場合には、現在の投影されているデータ・ボクセル
の値が、その投影の射線の突き当たる像平面の画素に既
に記憶されている和の値に加算され、次いでその結果の
和が最終的にその画素に突き当たるこのような射線の計
数値で割算されることが理解されよう。選択されたデー
タ・ボリューム内の各々のボクセルが逐次的にエントリ
されて像平面５６に投影されるとき、あるデータ・ボリ
ューム・ボクセル（例えば、ＤＶi,3,k ）はその関連す
る射線６２ｐに沿って投影されるが、所望の画素６０ａ
内に突き当たらず、従ってそのデータ値（例えば、強
度）が画素６０ａに対して現在記憶されているデータ値
と比較されない。特定の三次元の視角（θ，φ）におけ
るデータのその投影に対して、画素画素６０ａに対する
最大データ値がそのとき確立される。しかし、射線６２
ｐが事実上、別の像平面画素（例えば、画素画素６０
ｂ）内にある衝突位置６４ｐを持ち、その画素に記憶さ
れているデータ値と比較されて、比較後の大きい方の値
がその画素に対する記憶装置に戻される。全てのデータ
値は、新しい投影が取られるときにゼロにリセットされ
る。従って、像平面画素の各々は像投影手順の開始時に
リセットされ、（選択された被検体ボリューム５２の部
分によって設定されるような、全空間または選択された
部分内の）データ・ボリューム・ボクセルの全ては個別
に且つ逐次的に走査される。各々データ・ボクセルＤＶ
内のデータ値は、その１つの画素６０内で像平面５６に
突き当たる関連の射線６２により投影される。各々の画
素内の最大値は射線投射されたデータ・ボリューム・ボ
クセルの現在値との間で比較されて、その内の大きい方
の値が決定される。この大きい方の値は、次いで、最大
値像の一部分として記憶される。実際には、最大画素投
影の場合、新しく投射されたデータ・ボクセル値が、新
しく投射された射線が突き当たる像平面画素に対して既
に記憶されているデータ値よりも大きいときだけ、記憶
された最大値が変更される。
【００２１】上記の技術の別の面によれば、データ投影
がスケーリングされ（図３のステップ３６）、被検体ボ
リュームと像平面との間の非等方性が、逆投影の完了後
の一組の計算だけによって除かれる。ここで図５を参照
して説明すると、被検体ボリュームが実際のボリューム
（容積）であるのに対して、データ・ボリュームが抽象
的な概念であるので、第１の平面において、任意の視方
向６６が被検体ボリューム５２およびデータ・ボリュー
ム５４の両方に対して位置決めされる角度ψよりも異な
る角度γで立方体のデータ・ボリューム格子５４を表す
ことに起因するデータ投影の歪みの量を決定することが
必要である。各々のボクセルの見かけの寸法は、有効な
立て方向の角度（仰角）ψおよびγが変わるにつれて変
わろうとする。アスペクト比Ａ（被検体ボリューム５２
内の実際のスライスの厚さＴと同じ被検体ボリューム５
２内の実際の画素の大きさＳとの比として定義される）
が１でない（すなわち、被検体ボクセルがデータ・ボリ
ューム５４におけるように立方体のボクセルではないと
き、１より大きいか又は小さい）場合、立て方向の角度
ψおよびγが異なり、データ・ボリューム内の有効な立
て方向角度ψが、被検体ボリューム内の実際の立て方向
角度γとは異なる。データは、次式で得られる被検体の
立て方向角度に従って回転される。
がスケーリングされ（図３のステップ３６）、被検体ボ
リュームと像平面との間の非等方性が、逆投影の完了後
の一組の計算だけによって除かれる。ここで図５を参照
して説明すると、被検体ボリュームが実際のボリューム
（容積）であるのに対して、データ・ボリュームが抽象
的な概念であるので、第１の平面において、任意の視方
向６６が被検体ボリューム５２およびデータ・ボリュー
ム５４の両方に対して位置決めされる角度ψよりも異な
る角度γで立方体のデータ・ボリューム格子５４を表す
ことに起因するデータ投影の歪みの量を決定することが
必要である。各々のボクセルの見かけの寸法は、有効な
立て方向の角度（仰角）ψおよびγが変わるにつれて変
わろうとする。アスペクト比Ａ（被検体ボリューム５２
内の実際のスライスの厚さＴと同じ被検体ボリューム５
２内の実際の画素の大きさＳとの比として定義される）
が１でない（すなわち、被検体ボクセルがデータ・ボリ
ューム５４におけるように立方体のボクセルではないと
き、１より大きいか又は小さい）場合、立て方向の角度
ψおよびγが異なり、データ・ボリューム内の有効な立
て方向角度ψが、被検体ボリューム内の実際の立て方向
角度γとは異なる。データは、次式で得られる被検体の
立て方向角度に従って回転される。
【００２２】Ψ＝ｔａｎ-1［（１／Ａ）ｔａｎ（γ）］ その後、投影されたデータは、（回転が水平軸の周りに
なされた場合）被検体ボリューム内で正しい高さを持つ
ように、全ての投影されたデータの高さに立て方向スケ
ーリング係数を乗算することによって、スケーリングす
ることが出来る。古い投影像の高さＨを有効なスケーリ
ング係数ＥS により補正することが出来る。ここで、 ＥS ＝［（Ａｃｏｓγ）2＋ｓｉｎ2γ］1/2 であり、新しい高さＨ’はＨ’＝Ｈ・ＥS である。上記
と同じことが、回転が垂直軸の周りになされるときの幅
について当てはまる。
なされた場合）被検体ボリューム内で正しい高さを持つ
ように、全ての投影されたデータの高さに立て方向スケ
ーリング係数を乗算することによって、スケーリングす
ることが出来る。古い投影像の高さＨを有効なスケーリ
ング係数ＥS により補正することが出来る。ここで、 ＥS ＝［（Ａｃｏｓγ）2＋ｓｉｎ2γ］1/2 であり、新しい高さＨ’はＨ’＝Ｈ・ＥS である。上記
と同じことが、回転が垂直軸の周りになされるときの幅
について当てはまる。
【００２３】上記の関係を利用して、データ・ボリュー
ムの角度（α，β）を回転すると角度（θ，φ）にな
り、歪みが１つの軸に沿っているだけであるので、角度
θは角度αに等しい。３×３回転マトリクス［Ｍ］の要
素を決定することができ、２つの関係する回転角度が与
えられていると、これらの関係を使用してデータ・ボリ
ュームから像平面への変換が決定される。
ムの角度（α，β）を回転すると角度（θ，φ）にな
り、歪みが１つの軸に沿っているだけであるので、角度
θは角度αに等しい。３×３回転マトリクス［Ｍ］の要
素を決定することができ、２つの関係する回転角度が与
えられていると、これらの関係を使用してデータ・ボリ
ュームから像平面への変換が決定される。
【００２４】Ｘ’＝Ｍ１Ｘ＋Ｍ２Ｙ＋Ｍ３Ｚ＋ＸＯ Ｙ’＝Ｍ４Ｘ＋Ｍ５Ｙ＋Ｍ６Ｚ＋ＹＯ ここで、Ｍ１−Ｍ６は回転マトリクスの最初の２行（す
なわち、Ｍ１＝−ｓｉｎθ、Ｍ２＝ｃｏｓθｓｉｎψ、
Ｍ３＝０、Ｍ４＝−ｃｏｓθｓｉｎψ２、Ｍ５＝−ｓｉ
ｎθｓｉｎψ、Ｍ６＝ｃｏｓψ）であり、Ｘ’および
Ｙ’は投影された点の像平面上の位置であり、ＸＯおよ
びＹＯは像平面ＸおよびＹのオフセット（それぞれＸお
よびＹ最低値点を基準としている）であって、そこから
像平面の選択された部分が始まるオフセットである。デ
ータが像平面５６上に投影された後、等方性でない被検
体ボクセルの効果を補正するために像がスケーリングさ
れる。回転マトリクスの係数Ｍ１−Ｍ６を投影（所与の
θおよびφ）の始めに予め計算して（図３のステップ３
２）、全ての回転の計算のために使用することが出来る
ことが理解されよう。
なわち、Ｍ１＝−ｓｉｎθ、Ｍ２＝ｃｏｓθｓｉｎψ、
Ｍ３＝０、Ｍ４＝−ｃｏｓθｓｉｎψ２、Ｍ５＝−ｓｉ
ｎθｓｉｎψ、Ｍ６＝ｃｏｓψ）であり、Ｘ’および
Ｙ’は投影された点の像平面上の位置であり、ＸＯおよ
びＹＯは像平面ＸおよびＹのオフセット（それぞれＸお
よびＹ最低値点を基準としている）であって、そこから
像平面の選択された部分が始まるオフセットである。デ
ータが像平面５６上に投影された後、等方性でない被検
体ボクセルの効果を補正するために像がスケーリングさ
れる。回転マトリクスの係数Ｍ１−Ｍ６を投影（所与の
θおよびφ）の始めに予め計算して（図３のステップ３
２）、全ての回転の計算のために使用することが出来る
ことが理解されよう。
【００２５】図６は、主制御器８（または別の専用の処
理装置）内に設けられている上述のレイ・キャスティン
グ法を実行する手段を示す。このような手段は、シネ・
メモリ２４からデータ入力７０ａに受け取ったスライス
・データを記憶するための三次元データ・メモリ手段７
０を有する。各々の被検体ボクセルに関するデータが、
ＣＰＵ７４からボクセル・アドレス入力７０ｂに受け取
ったボクセル・アドレス入力情報に応答して、そのボク
セルのアドレスに記憶される。三次元データ・メモリ手
段７０が満たされたとき（被検体ボリューム５２からデ
ータ・ボリューム５４への全ての要求されたデータの転
送に対応する）、関心のある被検体ボリューム部分が選
択されて、そのＸ、ＹおよびＺ方向における開始コーナ
ーおよび範囲を設定するデータがＣＰＵ７４からアドレ
ス作成手段７２の入力７２ａへ送られる。アドレス作成
手段７２は、アドレス出力７２ｂに、選択された被検体
ボリューム内の各々のボクセルのＸ，Ｙ，Ｚアドレスを
逐次的に供給する。出力７２ｂは三次元データ・メモリ
手段７０の出力データ・アドレス入力７０ｃに接続され
ていて、その１つのボクセルに対する記憶された強度デ
ータがアドレスされて三次元データ・メモリ手段の出力
７０ｄから出力されるようにする。ボクセルのＸ，Ｙ，
Ｚアドレスはまた逐次的に回転パラメータ計算手段７６
の第１の入力７６ａにも供給される。回転パラメータ計
算手段７６は、ＣＰＵ７４を介して角度（α，β）情報
を、計算されたマトリクス要素Ｍ１−Ｍ６値として受け
取って、出力７６ｃに、選択された視角（θ，φ）で見
たときの被検体のＸ，Ｙ，Ｚ画素に対応する像平面画素
のアドレスＸ’，Ｙ’を供給する。視角（θ，φ）情報
はシステムに入力されて、ＣＰＵ７４によって処理され
る。その結果は視方向マトリクス手段７８の入力７８ｂ
および７８ｃに入力されて、その出力７８ａから回転パ
ラメータ計算手段７６へマトリクス要素Ｍ１−Ｍ６が供
給される。像平面画素のアドレスＸ’，Ｙ’は、像平面
メモリ手段８０として作用するフレーム・バッファのア
ドレス入力８０ａに現れる。同時に、データ・ボリュー
ムから像平面に投影された強度データが、三次元データ
・メモリ手段の出力７０ｄから像平面メモリ手段の新デ
ータ入力８０ｂに現れる。このデータはまた、データ比
較手段８２の新データ入力８２ａにも現れる。入力８０
ａのそのアドレスに対して像平面メモリ手段８０に前に
記憶されていた強度データが旧データ出力８０ｃに現
れ、従ってデータ比較手段の旧データ入力８２ｂに現れ
る。入力８２ａおよび８２ｂのそれぞれの新データおよ
び旧データがデータ比較手段８２内で比較され、入力８
２ａの新データが入力８２ｂの旧データよりも大きい場
合は、出力８２ｃが選定された論理状態（例えば、高論
理レベル）に作動される。出力８２ｃは像平面メモリ手
段８０の置換制御データ入力８０ｄに接続されている。
これにより、置換制御データ入力８０ｄが選定された論
理レベルにある場合、入力８０ａによって制御されたア
ドレスに記憶されるデータが、入力８０ｂの新データを
受け入れるように変更される。従って、記憶されている
データは（ＣＰＵ７４からの）データ／制御ポート８０
ｅの信号等によって最初にリセットされ、そして新デー
タが前に記憶された旧データよりも大きいという比較結
果に応答して、各々の像平面画素位置Ｘ’，Ｙ’に対し
て最大値のデータが記憶される。選択されたアドレスの
全てがアドレス作成手段７２によって逐次的に走査され
た後、像平面メモリ手段８０に記憶されているデータが
ＣＰＵ７４においてスケーリングされ、このスケーリン
グされたデータは表示、永久記憶または同様な目的のた
めに像平面メモリ手段８０から取り出すことが出来る。
理装置）内に設けられている上述のレイ・キャスティン
グ法を実行する手段を示す。このような手段は、シネ・
メモリ２４からデータ入力７０ａに受け取ったスライス
・データを記憶するための三次元データ・メモリ手段７
０を有する。各々の被検体ボクセルに関するデータが、
ＣＰＵ７４からボクセル・アドレス入力７０ｂに受け取
ったボクセル・アドレス入力情報に応答して、そのボク
セルのアドレスに記憶される。三次元データ・メモリ手
段７０が満たされたとき（被検体ボリューム５２からデ
ータ・ボリューム５４への全ての要求されたデータの転
送に対応する）、関心のある被検体ボリューム部分が選
択されて、そのＸ、ＹおよびＺ方向における開始コーナ
ーおよび範囲を設定するデータがＣＰＵ７４からアドレ
ス作成手段７２の入力７２ａへ送られる。アドレス作成
手段７２は、アドレス出力７２ｂに、選択された被検体
ボリューム内の各々のボクセルのＸ，Ｙ，Ｚアドレスを
逐次的に供給する。出力７２ｂは三次元データ・メモリ
手段７０の出力データ・アドレス入力７０ｃに接続され
ていて、その１つのボクセルに対する記憶された強度デ
ータがアドレスされて三次元データ・メモリ手段の出力
７０ｄから出力されるようにする。ボクセルのＸ，Ｙ，
Ｚアドレスはまた逐次的に回転パラメータ計算手段７６
の第１の入力７６ａにも供給される。回転パラメータ計
算手段７６は、ＣＰＵ７４を介して角度（α，β）情報
を、計算されたマトリクス要素Ｍ１−Ｍ６値として受け
取って、出力７６ｃに、選択された視角（θ，φ）で見
たときの被検体のＸ，Ｙ，Ｚ画素に対応する像平面画素
のアドレスＸ’，Ｙ’を供給する。視角（θ，φ）情報
はシステムに入力されて、ＣＰＵ７４によって処理され
る。その結果は視方向マトリクス手段７８の入力７８ｂ
および７８ｃに入力されて、その出力７８ａから回転パ
ラメータ計算手段７６へマトリクス要素Ｍ１−Ｍ６が供
給される。像平面画素のアドレスＸ’，Ｙ’は、像平面
メモリ手段８０として作用するフレーム・バッファのア
ドレス入力８０ａに現れる。同時に、データ・ボリュー
ムから像平面に投影された強度データが、三次元データ
・メモリ手段の出力７０ｄから像平面メモリ手段の新デ
ータ入力８０ｂに現れる。このデータはまた、データ比
較手段８２の新データ入力８２ａにも現れる。入力８０
ａのそのアドレスに対して像平面メモリ手段８０に前に
記憶されていた強度データが旧データ出力８０ｃに現
れ、従ってデータ比較手段の旧データ入力８２ｂに現れ
る。入力８２ａおよび８２ｂのそれぞれの新データおよ
び旧データがデータ比較手段８２内で比較され、入力８
２ａの新データが入力８２ｂの旧データよりも大きい場
合は、出力８２ｃが選定された論理状態（例えば、高論
理レベル）に作動される。出力８２ｃは像平面メモリ手
段８０の置換制御データ入力８０ｄに接続されている。
これにより、置換制御データ入力８０ｄが選定された論
理レベルにある場合、入力８０ａによって制御されたア
ドレスに記憶されるデータが、入力８０ｂの新データを
受け入れるように変更される。従って、記憶されている
データは（ＣＰＵ７４からの）データ／制御ポート８０
ｅの信号等によって最初にリセットされ、そして新デー
タが前に記憶された旧データよりも大きいという比較結
果に応答して、各々の像平面画素位置Ｘ’，Ｙ’に対し
て最大値のデータが記憶される。選択されたアドレスの
全てがアドレス作成手段７２によって逐次的に走査され
た後、像平面メモリ手段８０に記憶されているデータが
ＣＰＵ７４においてスケーリングされ、このスケーリン
グされたデータは表示、永久記憶または同様な目的のた
めに像平面メモリ手段８０から取り出すことが出来る。
【００２６】本発明によれば、図３に示されている方法
は、シネ・メモリから検索された関心のあるデータ・ボ
リュームに対して、カラー流れ速度データに適用され
る。投影像内の各々の画素は、所与の像平面上への投影
によって導き出されたそれぞれの変換された速度データ
を含む。更に、シネ・メモリがオペレータによって停止
（ｆｒｅｅｚｅ）されたときに、ＣＰＵ４２はシネ・メ
モリ２４の区分２４Ｂ内の多数の相次ぐアドレスにＸ−
Ｙメモリ１８からの最後のフレームを記憶させる。第１
の投影視角にたいする投影像データがシネ・メモリの区
分２４Ｂ内の第１のアドレスに書き込まれ、これにより
関心のある領域内の投影像データが背景フレーム上に重
畳される。このプロセスは、全ての投影像がシネ・メモ
リの区分２４Ｂ内に記憶されるまで、各々の角度増分に
対して繰り返される。各々の投影像フレームは、変換さ
れたデータを含む関心のある領域で構成され、またオプ
ションとして、関心のある領域を囲んでいて、関心のあ
る領域の変換されたデータによってオーバーライト（上
書き）されない背景フレーム・データより成る背景周辺
部を含む。背景の像は、各々の表示された投影がどの場
所から見たものであるかを一層明確にする。そこで、オ
ペレータは投影像の内の任意の１つを表示のために選択
することが出来る。更に、一連の投影像を表示モニタ上
にリプレイ（再生）して、被検体ボリュームをあたかも
観察者の前で回転しているかのように表示することが出
来る。
は、シネ・メモリから検索された関心のあるデータ・ボ
リュームに対して、カラー流れ速度データに適用され
る。投影像内の各々の画素は、所与の像平面上への投影
によって導き出されたそれぞれの変換された速度データ
を含む。更に、シネ・メモリがオペレータによって停止
（ｆｒｅｅｚｅ）されたときに、ＣＰＵ４２はシネ・メ
モリ２４の区分２４Ｂ内の多数の相次ぐアドレスにＸ−
Ｙメモリ１８からの最後のフレームを記憶させる。第１
の投影視角にたいする投影像データがシネ・メモリの区
分２４Ｂ内の第１のアドレスに書き込まれ、これにより
関心のある領域内の投影像データが背景フレーム上に重
畳される。このプロセスは、全ての投影像がシネ・メモ
リの区分２４Ｂ内に記憶されるまで、各々の角度増分に
対して繰り返される。各々の投影像フレームは、変換さ
れたデータを含む関心のある領域で構成され、またオプ
ションとして、関心のある領域を囲んでいて、関心のあ
る領域の変換されたデータによってオーバーライト（上
書き）されない背景フレーム・データより成る背景周辺
部を含む。背景の像は、各々の表示された投影がどの場
所から見たものであるかを一層明確にする。そこで、オ
ペレータは投影像の内の任意の１つを表示のために選択
することが出来る。更に、一連の投影像を表示モニタ上
にリプレイ（再生）して、被検体ボリュームをあたかも
観察者の前で回転しているかのように表示することが出
来る。
【００２７】本発明の好ましい実施態様では、超音波イ
メージング・システムは複数の異なる投影モードを有す
る。例えば、投影は最大または最小値の画素を含んでい
てよい。別のモードでは、面の表現（ｒｅｎｄｅｒｉｎ
ｇ）を行うようにレイ・キャスティング法を用いてもよ
い。速度のソース・データ・ボリュームを形成すると
き、そこから速度データを取ろうとするフレームまたは
走査を識別するために２つのタイプのゲート動作を使用
することが出来る。システム・オペレータが患者の心臓
サイクル内のある特定の点における血流に関心がある場
合、システムは、患者に接続された心臓モニタからの出
力を受け取るように接続される。各サイクル毎に、モニ
タは心臓サイクル波形内の所定の特徴の発生に応答して
信号を発生する。モニタからの各々の出力に応答して、
主制御器が、トリガ事象の発生時またはトリガ事象から
所定の遅延期間後にＸ−Ｙメモリ内に存在しているフレ
ームをシネ・メモリに記憶させる。従って、毎サイクル
当り１フレームがシネ・メモリに記憶される。この代わ
りに、心臓サイクル波形内の所定の特徴の発生に応答し
て、音響フレーム速度で多数の相次ぐフレームがシネ・
メモリに記憶される。
メージング・システムは複数の異なる投影モードを有す
る。例えば、投影は最大または最小値の画素を含んでい
てよい。別のモードでは、面の表現（ｒｅｎｄｅｒｉｎ
ｇ）を行うようにレイ・キャスティング法を用いてもよ
い。速度のソース・データ・ボリュームを形成すると
き、そこから速度データを取ろうとするフレームまたは
走査を識別するために２つのタイプのゲート動作を使用
することが出来る。システム・オペレータが患者の心臓
サイクル内のある特定の点における血流に関心がある場
合、システムは、患者に接続された心臓モニタからの出
力を受け取るように接続される。各サイクル毎に、モニ
タは心臓サイクル波形内の所定の特徴の発生に応答して
信号を発生する。モニタからの各々の出力に応答して、
主制御器が、トリガ事象の発生時またはトリガ事象から
所定の遅延期間後にＸ−Ｙメモリ内に存在しているフレ
ームをシネ・メモリに記憶させる。従って、毎サイクル
当り１フレームがシネ・メモリに記憶される。この代わ
りに、心臓サイクル波形内の所定の特徴の発生に応答し
て、音響フレーム速度で多数の相次ぐフレームがシネ・
メモリに記憶される。
【００２８】どの取得モードが用いられるかに拘わら
ず、ソース・データ・ボリュームは、各フレーム内の関
心のある領域に対応する画素データをシネ・メモリから
検索し、次いで画素データを処理して、所定の閾値範囲
内にある速度成分（例えば、非ゼロの速度成分）を持つ
画素データのみを取得することによって、形成される。
この速度情報は、次いで、表示のための投影速度像を再
構成するために、種々のイメージング平面に投影され
る。
ず、ソース・データ・ボリュームは、各フレーム内の関
心のある領域に対応する画素データをシネ・メモリから
検索し、次いで画素データを処理して、所定の閾値範囲
内にある速度成分（例えば、非ゼロの速度成分）を持つ
画素データのみを取得することによって、形成される。
この速度情報は、次いで、表示のための投影速度像を再
構成するために、種々のイメージング平面に投影され
る。
【００２９】上記の好ましい実施態様は例示の目的で開
示された。超音波イメージングまたはコンピュータ・グ
ラフィックスの分野における当業者には種々の変更およ
び変形を容易になし得よう。このような全ての変更およ
び変形は特許請求の範囲に包含されるものである。
示された。超音波イメージングまたはコンピュータ・グ
ラフィックスの分野における当業者には種々の変更およ
び変形を容易になし得よう。このような全ての変更およ
び変形は特許請求の範囲に包含されるものである。
【図１】実時間超音波イメージング・システム内の主要
な機能のサブシステムを示すブロック図である。
な機能のサブシステムを示すブロック図である。
【図２】本発明に従った画素強度データの相次ぐ立体的
投影よりなるフレームを再構成する手段を示すブロック
図である。
投影よりなるフレームを再構成する手段を示すブロック
図である。
【図３】本発明に従った画素強度データの相次ぐ立体的
投影よりなるフレームを再構成するためのアルゴリズム
のステップを示す流れ図である。
投影よりなるフレームを再構成するためのアルゴリズム
のステップを示す流れ図である。
【図４】従来技術に従ってレイ・キャスティングの逆投
影を立体的に表現する際に関係するサンプリングされた
関心のある被検体ボリューム、関連のデータ・ボリュー
ムおよび像投影平面を示す概略線図である。
影を立体的に表現する際に関係するサンプリングされた
関心のある被検体ボリューム、関連のデータ・ボリュー
ムおよび像投影平面を示す概略線図である。
【図５】被検体データおよびデータ・ボリュームの同様
なビュー（ｖｉｅｗ）に対応すると共に、三次元超音波
イメージングにおける必要なスケーリング定数を定める
のに有用である一対の幾何学的二次元形状を示す概略線
図である。
なビュー（ｖｉｅｗ）に対応すると共に、三次元超音波
イメージングにおける必要なスケーリング定数を定める
のに有用である一対の幾何学的二次元形状を示す概略線
図である。
【図６】三次元超音波イメージングにおける最大強度投
影を行う手段の概略ブロック図である。
影を行う手段の概略ブロック図である。
２ ビーム形成装置 ４ 処理装置 ４Ａ Ｂモード処理装置 ４Ｂ カラー流れ処理装置 ６ 走査変換器／表示制御器 ８ 主制御器 １２ 表示モニタ １４Ａ Ｂモード音響線メモリ １４Ｂ カラー流れ音響線メモリ １８ Ｘ−Ｙ表示メモリ ２０ 時間線／グラフィック処理装置及び表示メモリ ２２ ビデオ処理装置 ２４ シネ・メモリ ４２ 中央処理装置 ４４ ランダム・アクセス・メモリ ５０ サンプル（被検体） ５２ サンプル（被検体）ボリューム ５４ データ・ボリューム ５６ 像平面 ５８ 視線 ５８’ 射影線 ６０ 画素 ６２ 射線 ６４ 衝突位置 ７０ 三次元データ・メモリ手段 ７２ アドレス作成手段 ７４ ＣＰＵ ７６ 回転パラメータ計算手段 ７８ 視方向マトリクス手段 ８０ 像平面メモリ手段 ８２ データ比較手段
───────────────────────────────────────────────────── フロントページの続き (72)発明者 トッド・マイケル・ティルマン アメリカ合衆国、ウィスコンシン州、ウェ スト・ミルウォーキー、サウス・54ティー エイチ・ストリート、1514番 (72)発明者 パトリシア・アン・シューバート アメリカ合衆国、ウィスコンシン州、ミル ウォーキー、サウス・ショア・ドライブ、 2704エイ
Claims (9)
- 【請求項１】 動いている超音波散乱体を含む被検体ボ
リュームの三次元イメージング・システムにおいて、 超音波ビームを送信して、被検体ボリューム内の多数の
サンプル・ボリュームから反射された超音波エコーを検
出する超音波トランスジューサ・アレイ、 超音波散乱体によって反射された超音波エコーから少な
くとも部分的に導き出された速度データを取得する手段
であって、各々の速度データが前記多数のサンプル・ボ
リュームの内のそれぞれの１つのサンプル・ボリューム
に対応している手段、 前記多数のサンプル・ボリュームの各々に対して速度デ
ータを記憶するメモリ手段、 前記メモリ手段に記憶されている速度データから、所定
の範囲内にある速度データのみを含むソース速度データ
組を検索する手段、 前記ソース速度データ組の中の速度データを第１の像平
面へ投影して、第１の投影像を表す投影速度データ組を
形成する手段、 表示モニタ、および前記第１の投影像を前記表示モニタ
に表示させる手段、を含んでいることを特徴とする三次
元イメージング・システム。 - 【請求項２】 更に、前記ソース速度データ組の中の速
度データを、前記第１の像平面に対して回転されている
第２の像平面に投影して、第２の投影像を表す投影速度
データ組を形成する手段、および前記第２の投影像を前
記表示モニタに表示させる手段を含んでいる請求項１記
載のシステム。 - 【請求項３】 患者の心臓サイクル波形内の所定の特徴
に応答して、少なくとも１つのフレームの速度データが
前記メモリ手段に記憶される請求項１記載のシステム。 - 【請求項４】 動いている超音波散乱体を含む被検体ボ
リュームの三次元イメージング方法において、 被検体ボリューム内に超音波ビームを送信するステッ
プ、 被検体ボリューム内の多数のサンプル・ボリュームから
反射された超音波エコーを検出するステップ、 超音波散乱体によって反射された超音波エコーから少な
くとも部分的に導き出された速度データを取得するステ
ップであって、各々の速度データが前記多数のサンプル
・ボリュームの内のそれぞれ１つのサンプル・ボリュー
ムに対応しているステップ、 前記多数のサンプル・ボリュームの各々に対して速度デ
ータを記憶するステップ、 前記の記憶されている速度データから、所定の範囲内に
ある速度データのみを含むソース速度データ組を検索す
るステップ、 前記ソース速度データ組の中の速度データを第１の像平
面へ投影して、第１の投影像を表す投影速度データ組を
形成するステップ、および前記第１の投影像を表示する
ステップ、を含んでいることを特徴とする三次元イメー
ジング・システム。 - 【請求項５】 更に、前記ソース速度データ組の中の速
度データを、前記第１の像平面に対して回転されている
第２の像平面に投影して、第２の投影像を表す投影速度
データ組を形成するステップ、および前記第２の投影像
を表示するステップを含んでいる請求項４記載のシステ
ム。 - 【請求項６】 患者の心臓サイクル波形内の所定の特徴
に応答して、少なくとも１つのフレームの速度データが
記憶される請求項４記載のシステム。 - 【請求項７】 動いている超音波散乱体を含む被検体ボ
リュームの三次元イメージング方法において、 被検体ボリューム内に超音波ビームを送信するステッ
プ、 被検体ボリューム内の多数のサンプル・ボリュームから
反射された超音波エコーを検出するステップ、 超音波散乱体によって反射された超音波エコーから少な
くとも部分的に導き出された速度データを取得するステ
ップであって、各々の速度データが前記多数のサンプル
・ボリュームの内のそれぞれ１つのサンプル・ボリュー
ムに対応しているステップ、 前記多数のサンプル・ボリュームの各々に対して速度デ
ータを記憶するステップ、 前記の記憶されている速度データから、所定の範囲内に
ある速度データを含み且つ該所定の範囲の外にある速度
データを排除したソース速度データ組を検索するステッ
プ、 前記ソース速度データ組の中の速度データを第１の像平
面へ投影して、第１の投影像を表す投影速度データ組を
形成するステップ、および前記第１の投影像を表示する
ステップ、を含んでいることを特徴とする三次元イメー
ジング・システム。 - 【請求項８】 更に、前記ソース速度データ組の中の速
度データを、前記第１の像平面に対して回転されている
第２の像平面に投影して、第２の投影像を表す投影速度
データ組を形成するステップ、および前記第２の投影像
を表示するステップを含んでいる請求項７記載のシステ
ム。 - 【請求項９】 患者の心臓サイクル波形内の所定の特徴
に応答して、少なくとも１つのフレームの速度データが
記憶される請求項７記載のシステム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/852,613 US6126603A (en) | 1997-05-07 | 1997-05-07 | Method and apparatus for segmenting color flow mode data using velocity information in three-dimensional ultrasound imaging |
| US08/852613 | 1997-05-07 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH1128214A true JPH1128214A (ja) | 1999-02-02 |
Family
ID=25313787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10123449A Pending JPH1128214A (ja) | 1997-05-07 | 1998-05-06 | 三次元イメージング・システムおよび方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6126603A (ja) |
| JP (1) | JPH1128214A (ja) |
| DE (1) | DE19819793A1 (ja) |
| IL (1) | IL124170A (ja) |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6110120A (en) * | 1997-04-11 | 2000-08-29 | Acuson Corporation | Gated ultrasound imaging apparatus and method |
| US6171244B1 (en) * | 1997-12-31 | 2001-01-09 | Acuson Corporation | Ultrasonic system and method for storing data |
| US6245027B1 (en) * | 1998-04-10 | 2001-06-12 | Noam Alperin | Method of measuring intracranial pressure |
| JP3537358B2 (ja) * | 1999-08-30 | 2004-06-14 | 松下電器産業株式会社 | 超音波診断装置 |
| US6520915B1 (en) * | 2000-01-28 | 2003-02-18 | U-Systems, Inc. | Ultrasound imaging system with intrinsic doppler capability |
| EP1123687A3 (en) * | 2000-02-10 | 2004-02-04 | Aloka Co., Ltd. | Ultrasonic diagnostic apparatus |
| US20020045196A1 (en) * | 2000-05-12 | 2002-04-18 | Walt Mahoney | Methods of isolating trophoblast cells from maternal blood |
| US6599248B1 (en) * | 2001-03-20 | 2003-07-29 | Aloka | Method and apparatus for ultrasound diagnostic imaging |
| US7092558B2 (en) * | 2002-08-14 | 2006-08-15 | General Electric Company | Automated optimization of medical 3D visualizations |
| US7003493B2 (en) * | 2003-01-22 | 2006-02-21 | First Data Corporation | Direct payment with token |
| EP1803402A4 (en) * | 2004-10-20 | 2008-10-29 | Toshiba Kk | ULTRASONIC DIAGNOSTIC EQUIPMENT AND CONTROL METHOD THEREFOR |
| DE102005051783A1 (de) * | 2005-10-28 | 2007-05-03 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren und Vorrichtung zur bildgebenden Ultraschallprüfung an einem dreidimensionalen Werkstück |
| EP2005890A2 (en) * | 2006-03-31 | 2008-12-24 | National University Corporation Kyoto Institute of Technology | Image processing device, ultrasonic imaging device using the same, and image processing method |
| JP4881112B2 (ja) * | 2006-09-19 | 2012-02-22 | 株式会社東芝 | 超音波診断装置及び画像データ生成方法 |
| JP5160825B2 (ja) * | 2007-07-17 | 2013-03-13 | 日立アロカメディカル株式会社 | 超音波診断装置及び画像処理プログラム |
| JP5226978B2 (ja) * | 2007-07-17 | 2013-07-03 | 日立アロカメディカル株式会社 | 超音波診断装置及び画像処理プログラム |
| WO2014139032A1 (en) * | 2013-03-15 | 2014-09-18 | Colibri Technologies Inc. | Data display and processing algorithms for 3d imaging systems |
| EP2991556B1 (en) * | 2013-05-03 | 2018-08-15 | Sunnybrook Health Sciences Centre | Systems and methods for super-resolution ultrasound imaging |
| AU2019355381A1 (en) | 2018-10-02 | 2021-05-20 | Sony Corporation | Image processing device and method |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5226113A (en) * | 1989-10-30 | 1993-07-06 | General Electric Company | Method and apparatus for volumetric projection rendering using reverse ray casting |
| JP3187148B2 (ja) * | 1991-08-26 | 2001-07-11 | 株式会社東芝 | 超音波診断装置 |
| JP3248001B2 (ja) * | 1992-03-19 | 2002-01-21 | 株式会社日立メディコ | 三次元カラードプラ画像表示方法及びその装置 |
| US5622174A (en) * | 1992-10-02 | 1997-04-22 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and image displaying system |
| US5365929A (en) * | 1993-10-04 | 1994-11-22 | Advanced Technology Laboratories, Inc. | Multiple sample volume spectral Doppler |
| US5474073A (en) * | 1994-11-22 | 1995-12-12 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic scanning for three dimensional display |
| US5485842A (en) * | 1994-11-30 | 1996-01-23 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic scan conversion for three dimensional display processing |
| US5655535A (en) * | 1996-03-29 | 1997-08-12 | Siemens Medical Systems, Inc. | 3-Dimensional compound ultrasound field of view |
| US6059727A (en) * | 1995-06-15 | 2000-05-09 | The Regents Of The University Of Michigan | Method and apparatus for composition and display of three-dimensional image from two-dimensional ultrasound scan data |
| US5582173A (en) * | 1995-09-18 | 1996-12-10 | Siemens Medical Systems, Inc. | System and method for 3-D medical imaging using 2-D scan data |
-
1997
- 1997-05-07 US US08/852,613 patent/US6126603A/en not_active Expired - Lifetime
-
1998
- 1998-04-22 IL IL12417098A patent/IL124170A/en not_active IP Right Cessation
- 1998-05-04 DE DE19819793A patent/DE19819793A1/de not_active Ceased
- 1998-05-06 JP JP10123449A patent/JPH1128214A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| IL124170A (en) | 2002-11-10 |
| US6126603A (en) | 2000-10-03 |
| DE19819793A1 (de) | 1999-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4204095B2 (ja) | 被検体ボリュームの三次元イメージング・システムおよび方法 | |
| US5904653A (en) | Method and apparatus for three-dimensional ultrasound imaging combining intensity data with color flow velocity or power data | |
| US5779641A (en) | Method and apparatus for three-dimensional ultrasound imaging by projecting filtered pixel data | |
| JP4155618B2 (ja) | 超音波散乱媒体の三次元イメージング・システムおよび方法 | |
| US5934288A (en) | Method and apparatus for displaying 3D ultrasound data using three modes of operation | |
| JP4424799B2 (ja) | ３次元イメージング方法及びシステム | |
| US5895358A (en) | Method and apparatus for mapping color flow velocity data into display intensities | |
| US6413219B1 (en) | Three-dimensional ultrasound data display using multiple cut planes | |
| EP1046929B1 (en) | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering | |
| US5329929A (en) | Ultrasonic diagnostic apparatus | |
| JPH1128214A (ja) | 三次元イメージング・システムおよび方法 | |
| EP0563566B1 (en) | On-line acoustic densitometry tool for use with an ultrasonic imaging system | |
| US5865750A (en) | Method and apparatus for enhancing segmentation in three-dimensional ultrasound imaging | |
| US20040081340A1 (en) | Image processing apparatus and ultrasound diagnosis apparatus | |
| JP4767377B2 (ja) | 物体ボリュームの３次元イメージングのためのシステム及び方法 | |
| JPH11313822A (ja) | ３次元イメ―ジング・システム並びに走査平面の動きを追跡する方法及び装置 | |
| KR20000035616A (ko) | 평균 또는 중간 화소 투영을 이용한 속도 및 전력데이터의 3차원 초음파 화상 진단 | |
| EP1008111A1 (en) | 3d imaging from 2d scans | |
| Hottier et al. | 3D echography: status and perspective | |
| US6458082B1 (en) | System and method for the display of ultrasound data | |
| US6306092B1 (en) | Method and apparatus for calibrating rotational offsets in ultrasound transducer scans | |
| Watkin et al. | Three-dimensional reconstruction and enhancement of arbitrarily oriented and positioned 2D medical ultrasonic images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050428 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070829 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070904 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071109 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080226 |