JPH1128211A - 超音波散乱媒体の三次元イメージング・システムおよび方法 - Google Patents
超音波散乱媒体の三次元イメージング・システムおよび方法Info
- Publication number
- JPH1128211A JPH1128211A JP10123446A JP12344698A JPH1128211A JP H1128211 A JPH1128211 A JP H1128211A JP 10123446 A JP10123446 A JP 10123446A JP 12344698 A JP12344698 A JP 12344698A JP H1128211 A JPH1128211 A JP H1128211A
- Authority
- JP
- Japan
- Prior art keywords
- pixel data
- data
- volume
- pixel
- image
- 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
- 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
-
- 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/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/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
- 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/5206—Two-dimensional coordinated display of distance and direction; B-scan display
-
- 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/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
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 である。上記
と同じことが、回転が垂直軸の周りになされるときの幅
について当てはまる。
【００２７】上記の関係を利用して、データ・ボリュー
ムの角度（α，β）を回転すると角度（θ，φ）にな
り、歪みが１つの軸に沿っているだけであるので、角度
θは角度αに等しい。３×３回転マトリクス［Ｍ］の要
素を決定することができ、２つの関係する回転角度が与
えられていると、これらの関係を使用してデータ・ボリ
ュームから像平面への変換が決定される。
ムの角度（α，β）を回転すると角度（θ，φ）にな
り、歪みが１つの軸に沿っているだけであるので、角度
θは角度αに等しい。３×３回転マトリクス［Ｍ］の要
素を決定することができ、２つの関係する回転角度が与
えられていると、これらの関係を使用してデータ・ボリ
ュームから像平面への変換が決定される。
【００２８】Ｘ’＝Ｍ１Ｘ＋Ｍ２Ｙ＋Ｍ３Ｚ＋ＸＯ Ｙ’＝Ｍ４Ｘ＋Ｍ５Ｙ＋Ｍ６Ｚ＋ＹＯ ここで、Ｍ１−Ｍ６は回転マトリクスの最初の２行（す
なわち、Ｍ１＝−ｓｉｎθ、Ｍ２＝ｃｏｓθｓｉｎψ、
Ｍ３＝０、Ｍ４＝−ｃｏｓθｓｉｎψ２、Ｍ５＝−ｓｉ
ｎθｓｉｎψ、Ｍ６＝ｃｏｓψ）であり、Ｘ’および
Ｙ’は投影された点の像平面上の位置であり、ＸＯおよ
びＹＯは像平面ＸおよびＹのオフセット（それぞれＸお
よびＹ最低値点を基準としている）であって、そこから
像平面の選択された部分が始まるオフセットである。デ
ータが像平面５６上に投影された後、等方性でない被検
体ボクセルの効果を補正するために像がスケーリングさ
れる。回転マトリクスの係数Ｍ１−Ｍ６を投影（所与の
θおよびφ）の始めに予め計算して（図３のステップ３
２）、全ての回転の計算のために使用することが出来る
ことが理解されよう。図６は、主制御器８（または別の
専用の処理装置）内に設けられている上述のレイ・キャ
スティング法を実行する手段を示す。このような手段
は、シネ・メモリ２４からデータ入力７０ａに受け取っ
たスライス・データを記憶するための三次元データ・メ
モリ手段７０を有する。各々の被検体ボクセルに関する
データが、ＣＰＵ７４からボクセル・アドレス入力７０
ｂに受け取ったボクセル・アドレス入力情報に応答し
て、そのボクセルのアドレスに記憶される。三次元デー
タ・メモリ手段７０が満たされたとき（被検体ボリュー
ム５２からデータ・ボリューム５４への全ての要求され
たデータの転送に対応する）、関心のある被検体ボリュ
ーム部分が選択されて、そのＸ、ＹおよびＺ方向におけ
る開始コーナーおよび範囲を設定するデータがＣＰＵ７
４からアドレス作成手段７２の入力７２ａへ送られる。
アドレス作成手段７２は、アドレス出力７２ｂに、選択
された被検体ボリューム内の各々のボクセルのＸ，Ｙ，
Ｚアドレスを逐次的に供給する。出力７２ｂは三次元デ
ータ・メモリ手段７０の出力データ・アドレス入力７０
ｃに接続されていて、その１つのボクセルに対する記憶
された強度データがアドレスされて三次元データ・メモ
リ手段の出力７０ｄから出力されるようにする。ボクセ
ルのＸ，Ｙ，Ｚアドレスはまた逐次的に回転パラメータ
計算手段７６の第１の入力７６ａにも供給される。回転
パラメータ計算手段７６は、ＣＰＵ７４を介して角度
（α，β）情報を、計算されたマトリクス要素Ｍ１−Ｍ
６値として受け取って、出力７６ｃに、選択された視角
（θ，φ）で見たときの被検体のＸ，Ｙ，Ｚ画素に対応
する像平面画素のアドレスＸ’，Ｙ’を供給する。視角
（θ，φ）情報はシステムに入力されて、ＣＰＵ７４に
よって処理される。その結果は視方向マトリクス手段７
８の入力７８ｂおよび７８ｃに入力されて、その出力７
８ａから回転パラメータ計算手段７６へマトリクス要素
Ｍ１−Ｍ６が供給される。像平面画素のアドレスＸ’，
Ｙ’は、像平面メモリ手段８０として作用するフレーム
・バッファのアドレス入力８０ａに現れる。同時に、デ
ータ・ボリュームから像平面に投影された強度データ
が、三次元データ・メモリ手段の出力７０ｄから像平面
メモリ手段の新データ入力８０ｂに現れる。このデータ
はまた、データ比較手段８２の新データ入力８２ａにも
現れる。入力８０ａのそのアドレスに対して像平面メモ
リ手段８０に前に記憶されていた強度データが旧データ
出力８０ｃに現れ、従ってデータ比較手段の旧データ入
力８２ｂに現れる。入力８２ａおよび８２ｂのそれぞれ
の新データおよび旧データがデータ比較手段８２内で比
較され、入力８２ａの新データが入力８２ｂの旧データ
よりも大きい場合は、出力８２ｃが選定された論理状態
（例えば、高論理レベル）に作動される。出力８２ｃは
像平面メモリ手段８０の置換制御データ入力８０ｄに接
続されている。これにより、置換制御データ入力８０ｄ
が選定された論理レベルにある場合、入力８０ａによっ
て制御されたアドレスに記憶されるデータが、入力８０
ｂの新データを受け入れるように変更される。従って、
記憶されているデータは（ＣＰＵ７４からの）データ／
制御ポート８０ｅの信号等によって最初にリセットさ
れ、そして新データが前に記憶された旧データよりも大
きいという比較結果に応答して、各々の像平面画素位置
Ｘ’，Ｙ’に対して最大値のデータが記憶される。選択
されたアドレスの全てがアドレス作成手段７２によって
逐次的に走査された後、像平面メモリ手段８０に記憶さ
れているデータがＣＰＵ７４においてスケーリングさ
れ、このスケーリングされたデータは表示、永久記憶ま
たは同様な目的のために像平面メモリ手段８０から取り
出すことが出来る。
なわち、Ｍ１＝−ｓｉｎθ、Ｍ２＝ｃｏｓθｓｉｎψ、
Ｍ３＝０、Ｍ４＝−ｃｏｓθｓｉｎψ２、Ｍ５＝−ｓｉ
ｎθｓｉｎψ、Ｍ６＝ｃｏｓψ）であり、Ｘ’および
Ｙ’は投影された点の像平面上の位置であり、ＸＯおよ
びＹＯは像平面ＸおよびＹのオフセット（それぞれＸお
よびＹ最低値点を基準としている）であって、そこから
像平面の選択された部分が始まるオフセットである。デ
ータが像平面５６上に投影された後、等方性でない被検
体ボクセルの効果を補正するために像がスケーリングさ
れる。回転マトリクスの係数Ｍ１−Ｍ６を投影（所与の
θおよびφ）の始めに予め計算して（図３のステップ３
２）、全ての回転の計算のために使用することが出来る
ことが理解されよう。図６は、主制御器８（または別の
専用の処理装置）内に設けられている上述のレイ・キャ
スティング法を実行する手段を示す。このような手段
は、シネ・メモリ２４からデータ入力７０ａに受け取っ
たスライス・データを記憶するための三次元データ・メ
モリ手段７０を有する。各々の被検体ボクセルに関する
データが、ＣＰＵ７４からボクセル・アドレス入力７０
ｂに受け取ったボクセル・アドレス入力情報に応答し
て、そのボクセルのアドレスに記憶される。三次元デー
タ・メモリ手段７０が満たされたとき（被検体ボリュー
ム５２からデータ・ボリューム５４への全ての要求され
たデータの転送に対応する）、関心のある被検体ボリュ
ーム部分が選択されて、そのＸ、ＹおよびＺ方向におけ
る開始コーナーおよび範囲を設定するデータがＣＰＵ７
４からアドレス作成手段７２の入力７２ａへ送られる。
アドレス作成手段７２は、アドレス出力７２ｂに、選択
された被検体ボリューム内の各々のボクセルのＸ，Ｙ，
Ｚアドレスを逐次的に供給する。出力７２ｂは三次元デ
ータ・メモリ手段７０の出力データ・アドレス入力７０
ｃに接続されていて、その１つのボクセルに対する記憶
された強度データがアドレスされて三次元データ・メモ
リ手段の出力７０ｄから出力されるようにする。ボクセ
ルのＸ，Ｙ，Ｚアドレスはまた逐次的に回転パラメータ
計算手段７６の第１の入力７６ａにも供給される。回転
パラメータ計算手段７６は、ＣＰＵ７４を介して角度
（α，β）情報を、計算されたマトリクス要素Ｍ１−Ｍ
６値として受け取って、出力７６ｃに、選択された視角
（θ，φ）で見たときの被検体のＸ，Ｙ，Ｚ画素に対応
する像平面画素のアドレスＸ’，Ｙ’を供給する。視角
（θ，φ）情報はシステムに入力されて、ＣＰＵ７４に
よって処理される。その結果は視方向マトリクス手段７
８の入力７８ｂおよび７８ｃに入力されて、その出力７
８ａから回転パラメータ計算手段７６へマトリクス要素
Ｍ１−Ｍ６が供給される。像平面画素のアドレスＸ’，
Ｙ’は、像平面メモリ手段８０として作用するフレーム
・バッファのアドレス入力８０ａに現れる。同時に、デ
ータ・ボリュームから像平面に投影された強度データ
が、三次元データ・メモリ手段の出力７０ｄから像平面
メモリ手段の新データ入力８０ｂに現れる。このデータ
はまた、データ比較手段８２の新データ入力８２ａにも
現れる。入力８０ａのそのアドレスに対して像平面メモ
リ手段８０に前に記憶されていた強度データが旧データ
出力８０ｃに現れ、従ってデータ比較手段の旧データ入
力８２ｂに現れる。入力８２ａおよび８２ｂのそれぞれ
の新データおよび旧データがデータ比較手段８２内で比
較され、入力８２ａの新データが入力８２ｂの旧データ
よりも大きい場合は、出力８２ｃが選定された論理状態
（例えば、高論理レベル）に作動される。出力８２ｃは
像平面メモリ手段８０の置換制御データ入力８０ｄに接
続されている。これにより、置換制御データ入力８０ｄ
が選定された論理レベルにある場合、入力８０ａによっ
て制御されたアドレスに記憶されるデータが、入力８０
ｂの新データを受け入れるように変更される。従って、
記憶されているデータは（ＣＰＵ７４からの）データ／
制御ポート８０ｅの信号等によって最初にリセットさ
れ、そして新データが前に記憶された旧データよりも大
きいという比較結果に応答して、各々の像平面画素位置
Ｘ’，Ｙ’に対して最大値のデータが記憶される。選択
されたアドレスの全てがアドレス作成手段７２によって
逐次的に走査された後、像平面メモリ手段８０に記憶さ
れているデータがＣＰＵ７４においてスケーリングさ
れ、このスケーリングされたデータは表示、永久記憶ま
たは同様な目的のために像平面メモリ手段８０から取り
出すことが出来る。
【００２９】本発明の別の面によれば、表示の前に、ス
ケーリングされた像平面データが所望の輝度およびコン
トラスト範囲を達成するように写像（マッピング）され
る（図３のステップ３８）。三次元再構成のベースとな
るソース・フレームに対して関心のある領域を読み込み
ながら、所与の輝度を持つ画素の数についてのヒストグ
ラムがオプションとして主制御器８において作成され
る。この代わりに、ヒストグラムは投影像を使用して形
成することが出来る。同時に、最大画素強度が決定され
る。各々のビン（ｂｉｎ）内の画素が、全画素数の所与
のパーセントに達するまで計数される。このビン数が画
素閾値になる。次いで、意図した結果に応じて定めた画
素閾値より大きい又は小さい所望の輝度およびコントラ
スト範囲に各々の画素値が写像されるように、マップが
作成される。
ケーリングされた像平面データが所望の輝度およびコン
トラスト範囲を達成するように写像（マッピング）され
る（図３のステップ３８）。三次元再構成のベースとな
るソース・フレームに対して関心のある領域を読み込み
ながら、所与の輝度を持つ画素の数についてのヒストグ
ラムがオプションとして主制御器８において作成され
る。この代わりに、ヒストグラムは投影像を使用して形
成することが出来る。同時に、最大画素強度が決定され
る。各々のビン（ｂｉｎ）内の画素が、全画素数の所与
のパーセントに達するまで計数される。このビン数が画
素閾値になる。次いで、意図した結果に応じて定めた画
素閾値より大きい又は小さい所望の輝度およびコントラ
スト範囲に各々の画素値が写像されるように、マップが
作成される。
【００３０】本発明によれば、シネ・メモリ２４内に記
憶されている各々のスライス又はフレームからの関心の
ある領域内の画素がＣＰＵ４２（図１参照）によってコ
ンボリューション・フィルタリングされ、次いでメモリ
４４（図２参照）に記憶される。シネ・メモリに記憶さ
れている関心のある領域の画素は、シネ・メモリから読
み出して、相次ぐ画素として又は画素の配列としてコン
ボリューション・フィルタに供給することが出来る。Ｃ
ＰＵ４２内に設けられているコンボリューション・フィ
ルタによって実行される演算は、図９を参照して後で説
明する。
憶されている各々のスライス又はフレームからの関心の
ある領域内の画素がＣＰＵ４２（図１参照）によってコ
ンボリューション・フィルタリングされ、次いでメモリ
４４（図２参照）に記憶される。シネ・メモリに記憶さ
れている関心のある領域の画素は、シネ・メモリから読
み出して、相次ぐ画素として又は画素の配列としてコン
ボリューション・フィルタに供給することが出来る。Ｃ
ＰＵ４２内に設けられているコンボリューション・フィ
ルタによって実行される演算は、図９を参照して後で説
明する。
【００３１】本発明による装置は、それぞれの多数の重
み付け係数を受け取るための多数のタップを持つコンボ
リューション・フィルタを有する。この一組の重み付け
係数はランダム・アクセス・メモリ４４に記憶しておい
て、ＣＰＵ４２によって該メモリから検索することが出
来る。また、多数組の重み付け係数をメモり４４内にル
ックアップ・テーブルとして記憶しておいて、システム
・オペレータによる選択に応答して１組を検索するよう
にしてもよい。好ましい実施態様によれば、１組内の重
み付け係数の数は９個であり、各々の重み付け係数は多
ビット（例えば、８ビット）ディジタル数として表すこ
とが出来る。９個の重み付け係数（Ｗ1乃至Ｗ9 ）は３
×３カーネルとして配列され、これはカーネルの中心点
が各々のレンジ・ベクトルを（最初および最後のレンジ
・ベクトルを除き且つ各ベクトルの最初および最後のデ
ータ点を除いて）逐次的にステップ・ダウンされるよう
に関心のある全画素領域にわたって効果的に動かされ
る。
み付け係数を受け取るための多数のタップを持つコンボ
リューション・フィルタを有する。この一組の重み付け
係数はランダム・アクセス・メモリ４４に記憶しておい
て、ＣＰＵ４２によって該メモリから検索することが出
来る。また、多数組の重み付け係数をメモり４４内にル
ックアップ・テーブルとして記憶しておいて、システム
・オペレータによる選択に応答して１組を検索するよう
にしてもよい。好ましい実施態様によれば、１組内の重
み付け係数の数は９個であり、各々の重み付け係数は多
ビット（例えば、８ビット）ディジタル数として表すこ
とが出来る。９個の重み付け係数（Ｗ1乃至Ｗ9 ）は３
×３カーネルとして配列され、これはカーネルの中心点
が各々のレンジ・ベクトルを（最初および最後のレンジ
・ベクトルを除き且つ各ベクトルの最初および最後のデ
ータ点を除いて）逐次的にステップ・ダウンされるよう
に関心のある全画素領域にわたって効果的に動かされ
る。
【００３２】例えば、図７を参照して説明すると、９点
カーネル（図８参照）が最初の３つのベクトル、すなわ
ちＤＶ4,4,1 −ＤＶ4,1,1 、ＤＶ3,4,1 −ＤＶ3,1,1 お
よびＤＶ2,4,1 −ＤＶ2,1,1 に対してレンジを下向きに
動かされ、これにより画素データのこれらの最初の３つ
のベクトルは処理されて第１の新しいフィルタリングさ
れたベクトルを形成し、このベクトルは別の位置に記憶
される。第１のフィルタリング操作において、図７に示
されている配列の左上のコーナー内の３×３の画素配列
に重み付け係数が適用される。重み付け係数Ｗ1 が画素
ＤＶ4,4,1 を重み付けするために使用され、重み付け係
数Ｗ2 が画素ＤＶ3,4,1 を重み付けするために使用さ
れ、重み付け係数Ｗ3 が画素ＤＶ2,4,1 を重み付けする
ために使用され、重み付け係数Ｗ4 が画素ＤＶ4,3,1 を
重み付けするために使用されるとうように適用される。
コンボリューション・フィルタは、画素ＤＶ3,3,1 と置
換される処理された（すなわち、フィルタリングされ
た）画素を作成する。このフィルタリングされた画素Ｄ
Ｖ3,3,1 FILは、次のように計算することが出来る。
カーネル（図８参照）が最初の３つのベクトル、すなわ
ちＤＶ4,4,1 −ＤＶ4,1,1 、ＤＶ3,4,1 −ＤＶ3,1,1 お
よびＤＶ2,4,1 −ＤＶ2,1,1 に対してレンジを下向きに
動かされ、これにより画素データのこれらの最初の３つ
のベクトルは処理されて第１の新しいフィルタリングさ
れたベクトルを形成し、このベクトルは別の位置に記憶
される。第１のフィルタリング操作において、図７に示
されている配列の左上のコーナー内の３×３の画素配列
に重み付け係数が適用される。重み付け係数Ｗ1 が画素
ＤＶ4,4,1 を重み付けするために使用され、重み付け係
数Ｗ2 が画素ＤＶ3,4,1 を重み付けするために使用さ
れ、重み付け係数Ｗ3 が画素ＤＶ2,4,1 を重み付けする
ために使用され、重み付け係数Ｗ4 が画素ＤＶ4,3,1 を
重み付けするために使用されるとうように適用される。
コンボリューション・フィルタは、画素ＤＶ3,3,1 と置
換される処理された（すなわち、フィルタリングされ
た）画素を作成する。このフィルタリングされた画素Ｄ
Ｖ3,3,1 FILは、次のように計算することが出来る。
【００３３】 ＤＶ3,3,1 FIL＝ＤＶ4,4,1・Ｗ1 ＋ＤＶ3,4,1・Ｗ2 ＋ＤＶ2,4,1・Ｗ3 ＋ＤＶ4,3,1・Ｗ4 ＋ＤＶ3,3,1・Ｗ5 ＋ＤＶ2,3,1・Ｗ6 ＋ＤＶ4,2,1・Ｗ7 ＋ＤＶ3,2,1・Ｗ8 ＋ＤＶ2,2,1・Ｗ9 従って、コンボリューション・フィルタによって生じた
フィルタリングされた画素は、処理される画素とその周
囲の画素とを含む３×３の画素の和である。これらの画
素の各々は、重み付け係数Ｗ1 −Ｗ9 のそれぞれの１つ
により重み付けされている。次の段で、図８の９点カー
ネルが、図７の配列の左下のコーナーの３×３の画素配
列に適用されて、ＤＶ3,2,1 FILが求められる。この簡単
な例では、画素ＤＶ3,3,1 FILおよびＤＶ3,2,1 FILは第１
の新しいフィルタリングされたベクトルを形成する。縦
の一列の画素が５個より多い場合、このプロセスは列の
末端から二番目の画素がフィルタリングされるまで繰り
返される。その後、画素データの第２（ＤＶ3,4,1 −Ｄ
Ｖ3,1,1 ）、第３（ＤＶ2,4,1 −ＤＶ2,1,1 ）および第
４（ＤＶ1,4,1 −ＤＶ1,1,1 ）の隣接するベクトルが処
理されて、第２の新しいフィルタリングされたベクトル
が形成され、これは第１の新しいフィルタリングされた
ベクトルの隣に記憶される。このプロセスは、新しいフ
ィルタリングされた画素の像を形成する新しい一組のフ
ィルタリングされたベクトルが取得されて、新しい位置
に記憶されるまで、繰り返される。
フィルタリングされた画素は、処理される画素とその周
囲の画素とを含む３×３の画素の和である。これらの画
素の各々は、重み付け係数Ｗ1 −Ｗ9 のそれぞれの１つ
により重み付けされている。次の段で、図８の９点カー
ネルが、図７の配列の左下のコーナーの３×３の画素配
列に適用されて、ＤＶ3,2,1 FILが求められる。この簡単
な例では、画素ＤＶ3,3,1 FILおよびＤＶ3,2,1 FILは第１
の新しいフィルタリングされたベクトルを形成する。縦
の一列の画素が５個より多い場合、このプロセスは列の
末端から二番目の画素がフィルタリングされるまで繰り
返される。その後、画素データの第２（ＤＶ3,4,1 −Ｄ
Ｖ3,1,1 ）、第３（ＤＶ2,4,1 −ＤＶ2,1,1 ）および第
４（ＤＶ1,4,1 −ＤＶ1,1,1 ）の隣接するベクトルが処
理されて、第２の新しいフィルタリングされたベクトル
が形成され、これは第１の新しいフィルタリングされた
ベクトルの隣に記憶される。このプロセスは、新しいフ
ィルタリングされた画素の像を形成する新しい一組のフ
ィルタリングされたベクトルが取得されて、新しい位置
に記憶されるまで、繰り返される。
【００３４】ソース・データが強度データおよび速度ま
たはパワー・ドップラー・データを有している場合は、
上記のコンボリューション・フィルタ操作は各々の成分
に対して別々に実施することが出来る。本発明の好まし
い実施態様によれば、重み付け係数Ｗ5 ＝４であり、他
の８個の重み付け係数は１に等しい。このようなコンボ
リューション・フィルタは図９に示されるように構成す
ることが出来るが、ディジタル・フィルタに関する当業
者には他の構成を使用し得ることが理解されよう。
たはパワー・ドップラー・データを有している場合は、
上記のコンボリューション・フィルタ操作は各々の成分
に対して別々に実施することが出来る。本発明の好まし
い実施態様によれば、重み付け係数Ｗ5 ＝４であり、他
の８個の重み付け係数は１に等しい。このようなコンボ
リューション・フィルタは図９に示されるように構成す
ることが出来るが、ディジタル・フィルタに関する当業
者には他の構成を使用し得ることが理解されよう。
【００３５】図９を参照すると、コンボリューション・
フィルタが多数の画素遅延素子８６ａ−８６ｉ、多数の
乗算器８８ａ−８８ｉおよび多数の加算器９０ａ−９０
ｉを有するものとして表されている。画素遅延素子はシ
フトレジスタを構成するように接続されていて、各々の
遅延素子はメモリ４４から画素を読み出す速度に等しい
周期の間、画素値を記憶（格納）するようになってい
る。画素遅延素子８６ａ−８６ｃは１つのシフトレジス
タを構成するように接続され、画素遅延素子８６ｄ−８
６ｆは別の１つのシフトレジスタを構成するように接続
され、画素遅延素子８６ｇ−８６ｉは更に別の１つのシ
フトレジスタを構成するように接続されている。これら
のシフトレジスタは、乗算器８８ａ−８８ｉに対して３
×３の画素配列を提供するように使用される。
フィルタが多数の画素遅延素子８６ａ−８６ｉ、多数の
乗算器８８ａ−８８ｉおよび多数の加算器９０ａ−９０
ｉを有するものとして表されている。画素遅延素子はシ
フトレジスタを構成するように接続されていて、各々の
遅延素子はメモリ４４から画素を読み出す速度に等しい
周期の間、画素値を記憶（格納）するようになってい
る。画素遅延素子８６ａ−８６ｃは１つのシフトレジス
タを構成するように接続され、画素遅延素子８６ｄ−８
６ｆは別の１つのシフトレジスタを構成するように接続
され、画素遅延素子８６ｇ−８６ｉは更に別の１つのシ
フトレジスタを構成するように接続されている。これら
のシフトレジスタは、乗算器８８ａ−８８ｉに対して３
×３の画素配列を提供するように使用される。
【００３６】画素遅延素子８６ａ−８６ｃから成るシフ
トレジスタは、メモリ４４から直接に相次ぐ画素を受け
取るように接続されている。画素遅延素子８６ｄ−８６
ｆから成るシフトレジスタは、メモリ４４から読み出さ
れた相次ぐ画素をライン・メモリ９２を介して受け取る
ように接続されている。ライン・メモリ９２は、メモリ
４４から供給された画素に、１つの水平線の期間に等し
い遅延を与える。また、画素遅延素子８６ｇ−８６ｉか
ら成るシフトレジスタは、ライン・メモリ９２に縦続接
続されたライン・メモリ９４に接続されている。ライン
・メモリ９４は、ライン・メモリ９２によって既に遅延
されている画素に、１つの水平線の期間に等しい遅延を
与える。従って、ライン・メモリ９２および９４は、１
つの水平線内の画素が同じ縦の列内にあるが２つ先行す
る線内にある画素と時間的に一致するようにする。例え
ば、図７を参照して説明すると、画素ＤＶ3,2,1 がメモ
り４４から供給されるとき、画素ＤＶ3,3,1 がライン・
メモリ９２の出力に供給され、且つ画素ＤＶ3,4,1 がラ
イン・メモリ９４の出力に供給される。
トレジスタは、メモリ４４から直接に相次ぐ画素を受け
取るように接続されている。画素遅延素子８６ｄ−８６
ｆから成るシフトレジスタは、メモリ４４から読み出さ
れた相次ぐ画素をライン・メモリ９２を介して受け取る
ように接続されている。ライン・メモリ９２は、メモリ
４４から供給された画素に、１つの水平線の期間に等し
い遅延を与える。また、画素遅延素子８６ｇ−８６ｉか
ら成るシフトレジスタは、ライン・メモリ９２に縦続接
続されたライン・メモリ９４に接続されている。ライン
・メモリ９４は、ライン・メモリ９２によって既に遅延
されている画素に、１つの水平線の期間に等しい遅延を
与える。従って、ライン・メモリ９２および９４は、１
つの水平線内の画素が同じ縦の列内にあるが２つ先行す
る線内にある画素と時間的に一致するようにする。例え
ば、図７を参照して説明すると、画素ＤＶ3,2,1 がメモ
り４４から供給されるとき、画素ＤＶ3,3,1 がライン・
メモリ９２の出力に供給され、且つ画素ＤＶ3,4,1 がラ
イン・メモリ９４の出力に供給される。
【００３７】乗算器８８ａが重み付け係数Ｗ1 と遅延素
子８６ａの出力の画素データ値とを受け取って、これら
の入力の積を形成する。同様に、乗算器８８ｂが重み付
け係数Ｗ2 と遅延素子８６ｂの出力の画素データ値とを
受け取って、これらの入力の積を形成する。このように
して、処理している３×３の配列を構成する９個の画素
データ値と９個の重み付け係数とを使用することによ
り、９個の積が作られる。これらの９個の積は加算器９
０ａ−９０ｈによって加算されて、新しいフィルタリン
グされた画素データ値を形成する。このコンボリューシ
ョン・フィルタは、画素データの新しいフィルタリング
されたフレームが取得されるまで連続的に動作する。
子８６ａの出力の画素データ値とを受け取って、これら
の入力の積を形成する。同様に、乗算器８８ｂが重み付
け係数Ｗ2 と遅延素子８６ｂの出力の画素データ値とを
受け取って、これらの入力の積を形成する。このように
して、処理している３×３の配列を構成する９個の画素
データ値と９個の重み付け係数とを使用することによ
り、９個の積が作られる。これらの９個の積は加算器９
０ａ−９０ｈによって加算されて、新しいフィルタリン
グされた画素データ値を形成する。このコンボリューシ
ョン・フィルタは、画素データの新しいフィルタリング
されたフレームが取得されるまで連続的に動作する。
【００３８】本発明はコンボリューション・フィルタリ
ングを使用することに制限されない。代わりに、メジア
ン（ｍｅｄｉａｎ）フィルタまたはエッジ保存フィルタ
を使用して、画素データをフィルタリングしてもよい。
上記の投影法は、シネ・メモリから検索された関心のあ
るデータ・ボリュームに対して、Ｂモード強度データ
に、或いはカラー流れ速度またはパワー・データに適用
することが出来る。投影像内の各々の画素は、所与の像
平面上への投影によって導き出された、変換された強度
データおよび変換された速度またはパワー・データを含
む。更に、シネ・メモリがオペレータによって停止され
たときに、ＣＰＵ４２はオプションとしてシネ・メモリ
２４の区分２４Ｂ内の多数の相次ぐアドレスにＸ−Ｙメ
モリ１８からの最後のフレームを記憶させる。第１の投
影視角にたいする投影像データがシネ・メモリの区分２
４Ｂ内の第１のアドレスに書き込まれ、これにより関心
のある領域内の投影像データが背景のフレーム上に重畳
される。このプロセスは、全ての投影像がシネ・メモリ
の区分２４Ｂ内に記憶されるまで、各々の角度増分に対
して繰り返される。各々の投影像フレームは、変換され
たデータを含む関心のある領域で構成され、オプション
として関心のある領域を囲んでいて、関心のある領域の
変換されたデータによってオーバーライト（上書き）さ
れない背景フレーム・データより成る背景周辺部を含
む。背景の像は、各々の表示された投影がどの場所から
見たものであるかを一層明確にする。そこで、オペレー
タは投影像の内の任意の１つを表示のために選択するこ
とが出来る。更に、一連の投影像を表示モニタ上にリプ
レイ（再生）して、被検体ボリュームをあたかも観察者
の前で回転しているかのように表示することが出来る。
ングを使用することに制限されない。代わりに、メジア
ン（ｍｅｄｉａｎ）フィルタまたはエッジ保存フィルタ
を使用して、画素データをフィルタリングしてもよい。
上記の投影法は、シネ・メモリから検索された関心のあ
るデータ・ボリュームに対して、Ｂモード強度データ
に、或いはカラー流れ速度またはパワー・データに適用
することが出来る。投影像内の各々の画素は、所与の像
平面上への投影によって導き出された、変換された強度
データおよび変換された速度またはパワー・データを含
む。更に、シネ・メモリがオペレータによって停止され
たときに、ＣＰＵ４２はオプションとしてシネ・メモリ
２４の区分２４Ｂ内の多数の相次ぐアドレスにＸ−Ｙメ
モリ１８からの最後のフレームを記憶させる。第１の投
影視角にたいする投影像データがシネ・メモリの区分２
４Ｂ内の第１のアドレスに書き込まれ、これにより関心
のある領域内の投影像データが背景のフレーム上に重畳
される。このプロセスは、全ての投影像がシネ・メモリ
の区分２４Ｂ内に記憶されるまで、各々の角度増分に対
して繰り返される。各々の投影像フレームは、変換され
たデータを含む関心のある領域で構成され、オプション
として関心のある領域を囲んでいて、関心のある領域の
変換されたデータによってオーバーライト（上書き）さ
れない背景フレーム・データより成る背景周辺部を含
む。背景の像は、各々の表示された投影がどの場所から
見たものであるかを一層明確にする。そこで、オペレー
タは投影像の内の任意の１つを表示のために選択するこ
とが出来る。更に、一連の投影像を表示モニタ上にリプ
レイ（再生）して、被検体ボリュームをあたかも観察者
の前で回転しているかのように表示することが出来る。
【００３９】本発明の好ましい実施態様では、超音波イ
メージング・システムは複数の異なる投影モードを有す
る。例えば、投影は最大または最小値の画素を含んでい
てよい。或いは、画素データを逆転して、その最大値を
像平面に投影するような、血管のイメージングに有用な
モードを選択してもよい。更に別のモードでは、面の表
現を行うようにレイ・キャスティング法を用いてもよ
い。
メージング・システムは複数の異なる投影モードを有す
る。例えば、投影は最大または最小値の画素を含んでい
てよい。或いは、画素データを逆転して、その最大値を
像平面に投影するような、血管のイメージングに有用な
モードを選択してもよい。更に別のモードでは、面の表
現を行うようにレイ・キャスティング法を用いてもよ
い。
【００４０】上記の好ましい実施態様は例示の目的で開
示された。超音波イメージングまたはコンピュータ・グ
ラフィックスの分野における当業者には種々の変更およ
び変形を容易になし得よう。このような全ての変更およ
び変形は特許請求の範囲に包含されるものである。
示された。超音波イメージングまたはコンピュータ・グ
ラフィックスの分野における当業者には種々の変更およ
び変形を容易になし得よう。このような全ての変更およ
び変形は特許請求の範囲に包含されるものである。
【図１】実時間超音波イメージング・システム内の主要
な機能のサブシステムを示すブロック図である。
な機能のサブシステムを示すブロック図である。
【図２】本発明の好ましい実施態様に従って強度および
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成する手段を示すブロック図である。
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成する手段を示すブロック図である。
【図３】本発明の好ましい実施態様に従って強度および
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成するためのアルゴリズムのステップ
を示す流れ図である。
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成するためのアルゴリズムのステップ
を示す流れ図である。
【図４】従来技術に従ってレイ・キャスティングの逆投
影を立体的に表現する際に関係するサンプリングされた
関心のある被検体ボリューム、関連のデータ・ボリュー
ムおよび像投影平面を示す概略線図である。
影を立体的に表現する際に関係するサンプリングされた
関心のある被検体ボリューム、関連のデータ・ボリュー
ムおよび像投影平面を示す概略線図である。
【図５】被検体データおよびデータ・ボリュームの同様
なビュー（ｖｉｅｗ）に対応すると共に、三次元超音波
イメージングにおける必要なスケーリング定数を定める
のに有用である一対の幾何学的二次元形状を示す概略線
図である。
なビュー（ｖｉｅｗ）に対応すると共に、三次元超音波
イメージングにおける必要なスケーリング定数を定める
のに有用である一対の幾何学的二次元形状を示す概略線
図である。
【図６】三次元超音波イメージングにおける最大強度投
影を行う手段の概略ブロック図である。
影を行う手段の概略ブロック図である。
【図７】画素データの二次元配列の配列図である。
【図８】本発明に従ってコンボリューション・フィルタ
で使用される３×３の重み係数の組の配列図である。
で使用される３×３の重み係数の組の配列図である。
【図９】本発明の好ましい実施態様によるコンボリュー
ション・フィルタの概略ブロック図である。
ション・フィルタの概略ブロック図である。
２ ビーム形成装置 ４ 処理装置 ４Ａ Ｂモード処理装置 ４Ｂ カラー流れ処理装置 ６ 走査変換器／表示制御器 ８ 主制御器 １２ 表示モニタ １４Ａ Ｂモード音響線メモリ １４Ｂ カラー流れ音響線メモリ １８ Ｘ−Ｙ表示メモリ ２０ 時間線／グラフィック処理装置及び表示メモリ ２２ ビデオ処理装置 ２４ シネ・メモリ ４２ 中央処理装置 ４４ ランダム・アクセス・メモリ ５０ サンプル（被検体） ５２ サンプル（被検体）ボリューム ５４ データ・ボリューム ５６ 像平面 ５８ 視線 ５８’ 射影線 ６０ 画素 ６２ 射線 ６４ 衝突位置 ７０ 三次元データ・メモリ手段 ７２ アドレス作成手段 ７４ ＣＰＵ ７６ 回転パラメータ計算手段 ７８ 視方向マトリクス手段 ８０ 像平面メモリ手段 ８２ データ比較手段 ８６ 画素遅延素子 ８８ 乗算器 ９０ 加算器 ９２ ライン・メモリ ９４ ライン・メモリ
Claims (20)
- 【請求項１】 被検体ボリューム中の超音波散乱媒体の
三次元イメージング・システムにおいて、 超音波ビームを送信して、被検体ボリューム内の多数の
サンプル・ボリュームから反射された超音波エコーを検
出する超音波トランスジューサ・アレイ、 散乱媒体によって反射された超音波エコーから少なくと
も部分的に導き出された画素データを取得する取得手段
であって、各々の画素データが前記多数のサンプル・ボ
リュームの内のそれぞれ１つのサンプル・ボリュームに
対応している取得手段、 前記多数のサンプル・ボリュームの各々に対して前記取
得された画素データを記憶するメモリ手段、 前記メモリ手段から、前記被検体ボリューム内の関心の
あるボリュームに対応する一組の画素データを検索する
手段、 スペックルを低減するために、該一組の画素データをフ
ィルタリングする手段、 前記フィルタリングされた画素データの組を第１の像平
面へ投影して、第１の投影像を表す投影データの組を形
成する手段、 表示モニタ、および前記第１の投影像を前記表示モニタ
に表示させる手段、を含んでいることを特徴とする三次
元イメージング・システム。 - 【請求項２】 前記フィルタリング手段が、ｎ≧２とし
て、ｎ×ｎカーネルを持つコンボリューション・フィル
タである請求項１記載のシステム。 - 【請求項３】 ｎ＝３である請求項２記載のシステム。
- 【請求項４】 前記カーネルが次の重み付け係数：１１
１ １４１ １１１を持つ請求項３記載のシステム。 - 【請求項５】 前記画素データが強度データを有する請
求項１記載のシステム。 - 【請求項６】 前記画素データが流れ速度データを有す
る請求項１記載のシステム。 - 【請求項７】 前記画素データが流れパワー・データを
有する請求項１記載のシステム。 - 【請求項８】 更に、前記前記フィルタリングされた画
素データの組を、前記第１の像平面に対して回転されて
いる第２の像平面に投影して、第２の投影像を表す投影
データの組を形成する手段、および前記第２の投影像を
前記表示モニタに表示させる手段を含んでいる請求項１
記載のシステム。 - 【請求項９】 前記フィルタリング手段がメジアン・フ
ィルタである請求項１記載のシステム。 - 【請求項１０】 前記フィルタリング手段がエッジ保存
フィルタである請求項１記載のシステム。 - 【請求項１１】 被検体ボリューム中の超音波散乱媒体
の三次元イメージング・システムにおいて、 超音波ビームを送信して、被検体ボリュームと交差する
走査平面内にある多数のサンプル・ボリュームから反射
された超音波エコーを検出する超音波トランスジューサ
・アレイ、 散乱媒体によって反射された超音波エコーから少なくと
も部分的に導き出された画素データを取得する取得手段
であって、各々の画素データが前記走査平面内にある前
記多数のサンプル・ボリュームの内のそれぞれ１つのサ
ンプル・ボリュームに対応している取得手段、 多数の相次ぐ走査平面の各々に対して１フレームの画素
データを実時間で記憶するフレーム・メモリ手段、 前記フレーム・メモリ手段から相次いで出力される画素
データの多数のフレームを記憶するシネ・メモリ手段、 前記シネ・メモリ手段から、前記被検体ボリューム内の
関心のあるボリュームに対応する一組の画素データを検
索する手段、 スペックルを低減するために、該一組の画素データをフ
ィルタリングする手段、 前記画素データの組の多数の像平面へのそれぞれの投影
を表す、対応する多数のサブセットのフィルタリングさ
れ投影された画素データを作成する手段、 前記多数のサブセットのフィルタリングされ投影された
画素データを前記シネ・メモリに記憶させる手段、 前記多数の像平面の内の１つを選択する手段、 表示モニタ、および前記選択された像平面に対応する前
記サブセットのフィルタリングされ投影された画素デー
タを表示させる手段、を含んでいることを特徴とする三
次元イメージング・システム。 - 【請求項１２】 前記フィルタリング手段が、ｎ≧２と
して、ｎ×ｎカーネルを持つコンボリューション・フィ
ルタである請求項１１記載のシステム。 - 【請求項１３】 ｎ＝３である請求項１２記載のシステ
ム。 - 【請求項１４】 前記カーネルが重み付け係数：１１１
１４１ １１１を持つ請求項１３記載のシステム。 - 【請求項１５】 前記フィルタリング手段がメジアン・
フィルタである請求項１１記載のシステム。 - 【請求項１６】 前記フィルタリング手段がエッジ保存
フィルタである請求項１１記載のシステム。 - 【請求項１７】 被検体ボリューム中の超音波散乱媒体
の三次元イメージング方法において、 被検体ボリュームと交差する走査平面内に超音波ビーム
を送信するステップ、 前記走査平面内にある多数のサンプル・ボリュームから
反射された超音波エコーを検出するステップ、 前記被検体ボリュームにわたって前記走査平面を走査す
るステップ、 散乱媒体によって反射された超音波エコーから少なくと
も部分的に導き出された画素データを取得するステップ
であって、各々の画素データが前記多数のサンプル・ボ
リュームの内のそれぞれ１つのサンプル・ボリュームに
対応しているステップ、 前記多数のサンプル・ボリュームの各々に対して前記取
得された画素データを記憶するステップ、 前記記憶された画素データから、前記被検体ボリューム
内の関心のあるボリュームに対応する一組の画素データ
を検索するステップ、 スペックルを低減するために、該一組の画素データをフ
ィルタリングするステップ、 前記フィルタリングされた画素データの組を、前記走査
平面ではない第１の像平面へ投影して、第１の投影像を
表す投影データの組を形成するステップ、および前記第
１の投影像を表示するステップ、を含んでいることを特
徴とする三次元イメージング方法。 - 【請求項１８】 前記フィルタリングするステップが、
ｎ≧２として、ｎ×ｎカーネルを持つコンボリューショ
ン・フィルタにより前記画素データの組を処理する請求
項１７記載の方法。 - 【請求項１９】 ｎ＝３である請求項１８記載の方
法。。 - 【請求項２０】 前記カーネルが重み付け係数：１１１
１４１ １１１を持つ請求項１８記載のシステム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/852,773 US5779641A (en) | 1997-05-07 | 1997-05-07 | Method and apparatus for three-dimensional ultrasound imaging by projecting filtered pixel data |
| US08/852773 | 1997-05-07 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH1128211A true JPH1128211A (ja) | 1999-02-02 |
Family
ID=25314177
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10123446A Pending JPH1128211A (ja) | 1997-05-07 | 1998-05-06 | 超音波散乱媒体の三次元イメージング・システムおよび方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5779641A (ja) |
| JP (1) | JPH1128211A (ja) |
| DE (1) | DE19819831A1 (ja) |
| IL (1) | IL124171A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006305337A (ja) * | 2005-03-31 | 2006-11-09 | Toshiba Corp | 超音波診断装置及び超音波画像処理方法 |
| US8834371B2 (en) | 2005-03-31 | 2014-09-16 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus and ultrasound image processing program |
Families Citing this family (51)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5899999A (en) * | 1996-10-16 | 1999-05-04 | Microsoft Corporation | Iterative convolution filter particularly suited for use in an image classification and retrieval system |
| US6102864A (en) * | 1997-05-07 | 2000-08-15 | General Electric Company | Three-dimensional ultrasound imaging of velocity and power data using average or median pixel projections |
| US5967987A (en) * | 1997-12-18 | 1999-10-19 | Acuson Corporation | Ultrasonic system and method for measurement of fluid flow |
| US6262749B1 (en) | 1997-12-31 | 2001-07-17 | Acuson Corporation | Ultrasonic system and method for data transfer, storage and/or processing |
| US6171244B1 (en) | 1997-12-31 | 2001-01-09 | Acuson Corporation | Ultrasonic system and method for storing data |
| US6135960A (en) * | 1998-08-31 | 2000-10-24 | Holmberg; Linda Jean | High-resolution, three-dimensional whole body ultrasound imaging system |
| US5951479A (en) * | 1998-09-29 | 1999-09-14 | General Electric Company | Method and apparatus for synthetic transmit aperture imaging |
| US6283917B1 (en) * | 1998-10-01 | 2001-09-04 | Atl Ultrasound | Ultrasonic diagnostic imaging system with blurring corrected spatial compounding |
| US6210328B1 (en) * | 1998-10-01 | 2001-04-03 | Atl Ultrasound | Ultrasonic diagnostic imaging system with variable spatial compounding |
| US6155978A (en) * | 1998-12-09 | 2000-12-05 | General Electric Company | Three-dimensional imaging by projecting morphologically filtered pixel data |
| EP1034742A1 (de) * | 1999-03-09 | 2000-09-13 | Kreztechnik Aktiengesellschaft | Verfahren zur Untersuchung von Objekten mit Ultraschall |
| US6535623B1 (en) | 1999-04-15 | 2003-03-18 | Allen Robert Tannenbaum | Curvature based system for the segmentation and analysis of cardiac magnetic resonance images |
| US6484191B1 (en) | 1999-07-02 | 2002-11-19 | Aloka Co., Ltd. | Apparatus and method for the real-time calculation of local variance in images |
| US6530885B1 (en) * | 2000-03-17 | 2003-03-11 | Atl Ultrasound, Inc. | Spatially compounded three dimensional ultrasonic images |
| US7615008B2 (en) * | 2000-11-24 | 2009-11-10 | U-Systems, Inc. | Processing and displaying breast ultrasound information |
| US7597663B2 (en) | 2000-11-24 | 2009-10-06 | U-Systems, Inc. | Adjunctive ultrasound processing and display for breast cancer screening |
| US7103205B2 (en) * | 2000-11-24 | 2006-09-05 | U-Systems, Inc. | Breast cancer screening with ultrasound image overlays |
| US7940966B2 (en) | 2000-11-24 | 2011-05-10 | U-Systems, Inc. | Full-field breast image data processing and archiving |
| US7556602B2 (en) | 2000-11-24 | 2009-07-07 | U-Systems, Inc. | Breast cancer screening with adjunctive ultrasound mammography |
| KR100388407B1 (ko) * | 2001-04-27 | 2003-06-25 | 주식회사 메디슨 | 표시 장치의 화소에 대응하는 복셀에서 수신 집속하는 3차원 초음파 영상 시스템 |
| JP3898047B2 (ja) * | 2001-07-09 | 2007-03-28 | セイコーインスツル株式会社 | 血液レオロジー測定装置 |
| DE10141724B4 (de) * | 2001-08-25 | 2006-11-16 | Brüel & Kjaer Vibro GmbH | Verfahren zur Darstellung eines Histogramms auf einer Rasteranzeige |
| US7386153B2 (en) * | 2001-11-23 | 2008-06-10 | Infinitt Co., Ltd. | Medical image segmentation apparatus and method thereof |
| KR100439756B1 (ko) * | 2002-01-09 | 2004-07-12 | 주식회사 인피니트테크놀로지 | 3차원 가상내시경 화면 표시장치 및 그 방법 |
| JP2005536269A (ja) * | 2002-08-21 | 2005-12-02 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 適合可能な空間像合成による超音波撮像装置 |
| JP3944059B2 (ja) * | 2002-11-14 | 2007-07-11 | アロカ株式会社 | 超音波診断装置 |
| US6896657B2 (en) * | 2003-05-23 | 2005-05-24 | Scimed Life Systems, Inc. | Method and system for registering ultrasound image in three-dimensional coordinate system |
| US20050033123A1 (en) * | 2003-07-25 | 2005-02-10 | Siemens Medical Solutions Usa, Inc. | Region of interest methods and systems for ultrasound imaging |
| US20050110793A1 (en) * | 2003-11-21 | 2005-05-26 | Steen Erik N. | Methods and systems for graphics processing in a medical imaging system |
| JP4083719B2 (ja) * | 2004-08-18 | 2008-04-30 | ザイオソフト株式会社 | 画像再構成における逆投影方法を実行する画像処理装置 |
| EP1876959A1 (en) | 2005-04-25 | 2008-01-16 | Koninklijke Philips Electronics N.V. | Targeted additive gain tool for processing ultrasound images |
| US7991242B2 (en) | 2005-05-11 | 2011-08-02 | Optosecurity Inc. | Apparatus, method and system for screening receptacles and persons, having image distortion correction functionality |
| US20090174554A1 (en) | 2005-05-11 | 2009-07-09 | Eric Bergeron | Method and system for screening luggage items, cargo containers or persons |
| US20070070063A1 (en) * | 2005-09-29 | 2007-03-29 | Siemens Medical Solutions Usa, Inc. | Non-photorealistic volume rendering of ultrasonic data |
| US7899232B2 (en) | 2006-05-11 | 2011-03-01 | Optosecurity Inc. | Method and apparatus for providing threat image projection (TIP) in a luggage screening system, and luggage screening system implementing same |
| JP4864554B2 (ja) * | 2006-06-12 | 2012-02-01 | 株式会社東芝 | 超音波診断装置、医用画像処理装置、及び医用画像処理プログラム |
| US8494210B2 (en) | 2007-03-30 | 2013-07-23 | Optosecurity Inc. | User interface for use in security screening providing image enhancement capabilities and apparatus for implementing same |
| US8548568B2 (en) * | 2006-09-08 | 2013-10-01 | General Electric Company | Methods and apparatus for motion compensation |
| JP4881112B2 (ja) * | 2006-09-19 | 2012-02-22 | 株式会社東芝 | 超音波診断装置及び画像データ生成方法 |
| US20080081998A1 (en) * | 2006-10-03 | 2008-04-03 | General Electric Company | System and method for three-dimensional and four-dimensional contrast imaging |
| US8332016B2 (en) | 2009-08-04 | 2012-12-11 | Carl Zeiss Meditec, Inc. | Non-linear projections of 3-D medical imaging data |
| JP6025849B2 (ja) | 2011-09-07 | 2016-11-16 | ラピスカン システムズ、インコーポレイテッド | マニフェストデータをイメージング／検知処理に統合するｘ線検査システム |
| DE102014102367A1 (de) | 2014-02-24 | 2015-08-27 | Universität Rostock | Vorrichtung und Verfahren zur Hohlraumdetektion |
| US10456116B2 (en) * | 2014-09-30 | 2019-10-29 | Siemens Medical Solutions Usa, Inc. | Shadow suppression in ultrasound imaging |
| WO2016118947A1 (en) | 2015-01-23 | 2016-07-28 | The University Of North Carolina At Chapel Hill | Apparatuses, systems, and methods for preclinical ultrasound imaging of subjects |
| WO2016192114A1 (zh) * | 2015-06-05 | 2016-12-08 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体成像方法及超声流体成像系统 |
| CN107847214B (zh) * | 2015-08-04 | 2021-01-01 | 深圳迈瑞生物医疗电子股份有限公司 | 三维超声流体成像方法及系统 |
| WO2017042304A1 (en) * | 2015-09-10 | 2017-03-16 | Koninklijke Philips N.V. | An ultrasound system with wide depth and detailed viewing |
| EP3772702A3 (en) | 2016-02-22 | 2021-05-19 | Rapiscan Systems, Inc. | Methods for processing radiographic images |
| CN108078590B (zh) * | 2018-01-03 | 2021-02-09 | 声泰特（成都）科技有限公司 | 基于超声频谱多普勒的血流动力学可视化方法与系统 |
| US10803612B2 (en) * | 2018-09-25 | 2020-10-13 | General Electric Company | Method and system for structure recognition in three-dimensional ultrasound data based on volume renderings |
Family Cites Families (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5159931A (en) * | 1988-11-25 | 1992-11-03 | Riccardo Pini | Apparatus for obtaining a three-dimensional reconstruction of anatomic structures through the acquisition of echographic images |
| US5226113A (en) * | 1989-10-30 | 1993-07-06 | General Electric Company | Method and apparatus for volumetric projection rendering using reverse ray casting |
| US5282471A (en) * | 1991-07-31 | 1994-02-01 | Kabushiki Kaisha Toshiba | Ultrasonic imaging system capable of displaying 3-dimensional angiogram in real time mode |
| US5241371A (en) * | 1991-10-30 | 1993-08-31 | Pioneer Electronic Corporation | Image display control apparatus and still picture reproducing system |
| US5327894A (en) * | 1993-03-01 | 1994-07-12 | General Electric Company | Wall filter using circular convolution for a color flow imaging system |
| US5601084A (en) * | 1993-06-23 | 1997-02-11 | University Of Washington | Determining cardiac wall thickness and motion by imaging and three-dimensional modeling |
| US5365929A (en) * | 1993-10-04 | 1994-11-22 | Advanced Technology Laboratories, Inc. | Multiple sample volume spectral Doppler |
| US5409007A (en) * | 1993-11-26 | 1995-04-25 | General Electric Company | Filter to reduce speckle artifact in ultrasound imaging |
| 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 |
| WO1997000482A1 (en) * | 1995-06-15 | 1997-01-03 | The Regents Of The University Of Michigan | Method and apparatus for composition and display of three-dimensional image from two-dimensional ultrasound |
| US5566674A (en) * | 1995-06-30 | 1996-10-22 | Siemens Medical Systems, Inc. | Method and apparatus for reducing ultrasound image shadowing and speckle |
| US5582176A (en) * | 1995-08-15 | 1996-12-10 | Medasonics | Methods and apparatus for automatically determining edge frequency in doppler ultrasound signals |
| 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,773 patent/US5779641A/en not_active Expired - Lifetime
-
1998
- 1998-04-22 IL IL12417198A patent/IL124171A/xx not_active IP Right Cessation
- 1998-05-04 DE DE19819831A patent/DE19819831A1/de not_active Ceased
- 1998-05-06 JP JP10123446A patent/JPH1128211A/ja active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006305337A (ja) * | 2005-03-31 | 2006-11-09 | Toshiba Corp | 超音波診断装置及び超音波画像処理方法 |
| US8834371B2 (en) | 2005-03-31 | 2014-09-16 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus and ultrasound image processing program |
Also Published As
| Publication number | Publication date |
|---|---|
| IL124171A (en) | 2000-10-31 |
| US5779641A (en) | 1998-07-14 |
| DE19819831A1 (de) | 1998-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH1128211A (ja) | 超音波散乱媒体の三次元イメージング・システムおよび方法 | |
| JP4155618B2 (ja) | 超音波散乱媒体の三次元イメージング・システムおよび方法 | |
| US5904653A (en) | Method and apparatus for three-dimensional ultrasound imaging combining intensity data with color flow velocity or power data | |
| JP4204095B2 (ja) | 被検体ボリュームの三次元イメージング・システムおよび方法 | |
| JP4424799B2 (ja) | ３次元イメージング方法及びシステム | |
| US8021300B2 (en) | Three-dimensional fly-through systems and methods using ultrasound data | |
| US6102861A (en) | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering | |
| US5895358A (en) | Method and apparatus for mapping color flow velocity data into display intensities | |
| US20040081340A1 (en) | Image processing apparatus and ultrasound diagnosis apparatus | |
| US5329929A (en) | Ultrasonic diagnostic apparatus | |
| US6193660B1 (en) | Medical diagnostic ultrasound system and method for region of interest determination | |
| US5546807A (en) | High speed volumetric ultrasound imaging system | |
| US6511426B1 (en) | Medical diagnostic ultrasound system and method for versatile processing | |
| US5865750A (en) | Method and apparatus for enhancing segmentation in three-dimensional ultrasound imaging | |
| US6102858A (en) | Method and apparatus for three-dimensional ultrasound imaging using contrast agents and harmonic echoes | |
| US6126603A (en) | Method and apparatus for segmenting color flow mode data using velocity information in three-dimensional ultrasound imaging | |
| US20120053464A1 (en) | Ultrasonic image processing apparatus | |
| JP2000152938A (ja) | 物体ボリュ―ムをイメ―ジングするシステム及び方法 | |
| JP2003061956A (ja) | 超音波診断装置、医用診断装置及び画像処理方法 | |
| JP4653324B2 (ja) | 画像表示装置、画像表示プログラム、画像処理装置、及び医用画像診断装置 | |
| US6458082B1 (en) | System and method for the display of ultrasound data | |
| JP3936450B2 (ja) | 投影画像生成装置及び医用画像装置 | |
| US5865752A (en) | Method and apparatus for ultrasound imaging using normalized difference between successive frames | |
| JPH09204513A (ja) | 投影画像処理方法及び投影画像処理装置 | |
| 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: 20070815 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070821 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071011 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071113 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080117 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080304 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080328 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080729 |