JPH1147134A - 超音波散乱媒体の三次元イメージング・システムおよび方法 - Google Patents
超音波散乱媒体の三次元イメージング・システムおよび方法Info
- Publication number
- JPH1147134A JPH1147134A JP10123450A JP12345098A JPH1147134A JP H1147134 A JPH1147134 A JP H1147134A JP 10123450 A JP10123450 A JP 10123450A JP 12345098 A JP12345098 A JP 12345098A JP H1147134 A JPH1147134 A JP H1147134A
- Authority
- JP
- Japan
- Prior art keywords
- intensity
- data
- pixel
- range
- ranges
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
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
-
- 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/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/52073—Production of cursor lines, markers or indicia by electronic means
-
- 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/52046—Techniques for image enhancement involving transmitter or receiver
-
- 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
-
- 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/52074—Composite displays, e.g. split-screen displays; Combination of multiple images or of images and alphanumeric tabular information
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
- Image Generation (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 である。上記
と同じことが、回転が垂直軸の周りになされるときの幅
について当てはまる。
【００２５】上記の関係を利用して、データ・ボリュー
ムの角度（α，β）を回転すると角度（θ，φ）にな
り、歪みが１つの軸に沿っているだけであるので、角度
θは角度αに等しい。３×３回転マトリクス［Ｍ］の要
素を決定することができ、２つの関係する回転角度が与
えられていると、これらの関係を使用してデータ・ボリ
ュームから像平面への変換が決定される。
ムの角度（α，β）を回転すると角度（θ，φ）にな
り、歪みが１つの軸に沿っているだけであるので、角度
θは角度αに等しい。３×３回転マトリクス［Ｍ］の要
素を決定することができ、２つの関係する回転角度が与
えられていると、これらの関係を使用してデータ・ボリ
ュームから像平面への変換が決定される。
【００２６】Ｘ’＝Ｍ１Ｘ＋Ｍ２Ｙ＋Ｍ３Ｚ＋ＸＯ Ｙ’＝Ｍ４Ｘ＋Ｍ５Ｙ＋Ｍ６Ｚ＋ＹＯ ここで、Ｍ１−Ｍ６は回転マトリクスの最初の２行（す
なわち、Ｍ１＝−ｓｉｎθ、Ｍ２＝ｃｏｓθｓｉｎψ、
Ｍ３＝０、Ｍ４＝−ｃｏｓθｓｉｎψ２、Ｍ５＝−ｓｉ
ｎθｓｉｎψ、Ｍ６＝ｃｏｓψ）であり、Ｘ’および
Ｙ’は投影された点の像平面上の位置であり、ＸＯおよ
びＹＯは像平面ＸおよびＹのオフセット（それぞれＸお
よびＹ最低値点を基準としている）であって、そこから
像平面の選択された部分が始まるオフセットである。デ
ータが像平面５６上に投影された後、等方性でない被検
体ボクセルの効果を補正するために像がスケーリングさ
れる。回転マトリクスの係数Ｍ１−Ｍ６を投影（所与の
θおよびφ）の始めに予め計算して（図３のステップ３
２）、全ての回転の計算のために使用することが出来る
ことが理解されよう。図６は、主制御器８（または別の
専用の処理装置）内に設けられている上述のレイ・キャ
スティング法を実行する手段を示す。このような手段
は、シネ・メモリ２４からデータ入力７０ａに受け取っ
たスライス・データを記憶するための三次元データ・メ
モリ手段７０を有する。各々の被検体ボクセルに関する
データが、ＣＰＵ７４からボクセル・アドレス入力７０
ｂに受け取ったボクセル・アドレス入力情報に応答し
て、そのボクセルのアドレスに記憶される。三次元デー
タ・メモリ手段７０が満たされたとき（被検体ボリュー
ム５２からデータ・ボリューム５４への全ての要求され
たデータの転送に対応する）、関心のある被検体ボリュ
ーム部分が選択されて、そのＸ、ＹおよびＺ方向におけ
る開始コーナーおよび範囲を設定するデータがＣＰＵ７
４からアドレス作成手段７２の入力７２ａへ送られる。
アドレス作成手段７２は、アドレス出力７２ｂに、選択
された被検体ボリューム内の各々のボクセルのＸ，Ｙ，
Ｚアドレスを逐次的に供給する。出力７２ｂは三次元デ
ータ・メモリ手段７０の出力データ・アドレス入力７０
ｃに接続されていて、その１つのボクセルに対する記憶
された強度データがアドレスされて三次元データ・メモ
リ手段の出力７０ｄから出力されるようにする。ボクセ
ルのＸ，Ｙ，Ｚアドレスはまた逐次的に回転パラメータ
計算手段７６の第１の入力７６ａにも供給される。回転
パラメータ計算手段７６は、ＣＰＵ７４を介して角度
（α，β）情報を、計算されたマトリクス要素Ｍ１−Ｍ
６値として受け取って、出力７６ｃに、選択された視角
（θ，φ）で見たときの被検体のＸ，Ｙ，Ｚ画素に対応
する像平面画素のアドレスＸ’，Ｙ’を供給する。視角
（θ，φ）情報はシステムに入力されて、ＣＰＵ７４に
よって処理される。その結果は視方向マトリクス手段７
８の入力７８ｂおよび７８ｃに入力されて、その出力７
８ａから回転パラメータ計算手段７６へマトリクス要素
Ｍ１−Ｍ６が供給される。像平面画素のアドレスＸ’，
Ｙ’は、像平面メモリ手段８０として作用するフレーム
・バッファのアドレス入力８０ａに現れる。同時に、デ
ータ・ボリュームから像平面に投影された強度データ
が、三次元データ・メモリ手段の出力７０ｄから像平面
メモリ手段の新データ入力８０ｂに現れる。このデータ
はまた、データ比較手段８２の新データ入力８２ａにも
現れる。入力８０ａのそのアドレスに対して像平面メモ
リ手段８０に前に記憶されていた強度データが旧データ
出力８０ｃに現れ、従ってデータ比較手段の旧データ入
力８２ｂに現れる。入力８２ａおよび８２ｂのそれぞれ
の新データおよび旧データがデータ比較手段８２内で比
較され、入力８２ａの新データが入力８２ｂの旧データ
よりも大きい場合は、出力８２ｃが選定された論理状態
（例えば、高論理レベル）に作動される。出力８２ｃは
像平面メモリ手段８０の置換制御データ入力８０ｄに接
続されている。これにより、置換制御データ入力８０ｄ
が選定された論理レベルにある場合、入力８０ａによっ
て制御されたアドレスに記憶されるデータが、入力８０
ｂの新データを受け入れるように変更される。従って、
記憶されているデータは（ＣＰＵ７４からの）データ／
制御ポート８０ｅの信号等によって最初にリセットさ
れ、そして新データが前に記憶された旧データよりも大
きいという比較結果に応答して、各々の像平面画素位置
Ｘ’，Ｙ’に対して最大値のデータが記憶される。選択
されたアドレスの全てがアドレス作成手段７２によって
逐次的に走査された後、像平面メモリ手段８０に記憶さ
れているデータがＣＰＵ７４においてスケーリングさ
れ、このスケーリングされたデータは表示、永久記憶ま
たは同様な目的のために像平面メモリ手段８０から取り
出すことが出来る。
なわち、Ｍ１＝−ｓｉｎθ、Ｍ２＝ｃｏｓθｓｉｎψ、
Ｍ３＝０、Ｍ４＝−ｃｏｓθｓｉｎψ２、Ｍ５＝−ｓｉ
ｎθｓｉｎψ、Ｍ６＝ｃｏｓψ）であり、Ｘ’および
Ｙ’は投影された点の像平面上の位置であり、ＸＯおよ
びＹＯは像平面ＸおよびＹのオフセット（それぞれＸお
よびＹ最低値点を基準としている）であって、そこから
像平面の選択された部分が始まるオフセットである。デ
ータが像平面５６上に投影された後、等方性でない被検
体ボクセルの効果を補正するために像がスケーリングさ
れる。回転マトリクスの係数Ｍ１−Ｍ６を投影（所与の
θおよびφ）の始めに予め計算して（図３のステップ３
２）、全ての回転の計算のために使用することが出来る
ことが理解されよう。図６は、主制御器８（または別の
専用の処理装置）内に設けられている上述のレイ・キャ
スティング法を実行する手段を示す。このような手段
は、シネ・メモリ２４からデータ入力７０ａに受け取っ
たスライス・データを記憶するための三次元データ・メ
モリ手段７０を有する。各々の被検体ボクセルに関する
データが、ＣＰＵ７４からボクセル・アドレス入力７０
ｂに受け取ったボクセル・アドレス入力情報に応答し
て、そのボクセルのアドレスに記憶される。三次元デー
タ・メモリ手段７０が満たされたとき（被検体ボリュー
ム５２からデータ・ボリューム５４への全ての要求され
たデータの転送に対応する）、関心のある被検体ボリュ
ーム部分が選択されて、そのＸ、ＹおよびＺ方向におけ
る開始コーナーおよび範囲を設定するデータがＣＰＵ７
４からアドレス作成手段７２の入力７２ａへ送られる。
アドレス作成手段７２は、アドレス出力７２ｂに、選択
された被検体ボリューム内の各々のボクセルのＸ，Ｙ，
Ｚアドレスを逐次的に供給する。出力７２ｂは三次元デ
ータ・メモリ手段７０の出力データ・アドレス入力７０
ｃに接続されていて、その１つのボクセルに対する記憶
された強度データがアドレスされて三次元データ・メモ
リ手段の出力７０ｄから出力されるようにする。ボクセ
ルのＸ，Ｙ，Ｚアドレスはまた逐次的に回転パラメータ
計算手段７６の第１の入力７６ａにも供給される。回転
パラメータ計算手段７６は、ＣＰＵ７４を介して角度
（α，β）情報を、計算されたマトリクス要素Ｍ１−Ｍ
６値として受け取って、出力７６ｃに、選択された視角
（θ，φ）で見たときの被検体のＸ，Ｙ，Ｚ画素に対応
する像平面画素のアドレスＸ’，Ｙ’を供給する。視角
（θ，φ）情報はシステムに入力されて、ＣＰＵ７４に
よって処理される。その結果は視方向マトリクス手段７
８の入力７８ｂおよび７８ｃに入力されて、その出力７
８ａから回転パラメータ計算手段７６へマトリクス要素
Ｍ１−Ｍ６が供給される。像平面画素のアドレスＸ’，
Ｙ’は、像平面メモリ手段８０として作用するフレーム
・バッファのアドレス入力８０ａに現れる。同時に、デ
ータ・ボリュームから像平面に投影された強度データ
が、三次元データ・メモリ手段の出力７０ｄから像平面
メモリ手段の新データ入力８０ｂに現れる。このデータ
はまた、データ比較手段８２の新データ入力８２ａにも
現れる。入力８０ａのそのアドレスに対して像平面メモ
リ手段８０に前に記憶されていた強度データが旧データ
出力８０ｃに現れ、従ってデータ比較手段の旧データ入
力８２ｂに現れる。入力８２ａおよび８２ｂのそれぞれ
の新データおよび旧データがデータ比較手段８２内で比
較され、入力８２ａの新データが入力８２ｂの旧データ
よりも大きい場合は、出力８２ｃが選定された論理状態
（例えば、高論理レベル）に作動される。出力８２ｃは
像平面メモリ手段８０の置換制御データ入力８０ｄに接
続されている。これにより、置換制御データ入力８０ｄ
が選定された論理レベルにある場合、入力８０ａによっ
て制御されたアドレスに記憶されるデータが、入力８０
ｂの新データを受け入れるように変更される。従って、
記憶されているデータは（ＣＰＵ７４からの）データ／
制御ポート８０ｅの信号等によって最初にリセットさ
れ、そして新データが前に記憶された旧データよりも大
きいという比較結果に応答して、各々の像平面画素位置
Ｘ’，Ｙ’に対して最大値のデータが記憶される。選択
されたアドレスの全てがアドレス作成手段７２によって
逐次的に走査された後、像平面メモリ手段８０に記憶さ
れているデータがＣＰＵ７４においてスケーリングさ
れ、このスケーリングされたデータは表示、永久記憶ま
たは同様な目的のために像平面メモリ手段８０から取り
出すことが出来る。
【００２７】本発明の別の面によれば、表示の前に、ス
ケーリングされた像平面データが所望の輝度およびコン
トラスト範囲を達成するように写像される（図３のステ
ップ３８）。本発明によれば、投影像のコントラスト
は、主制御器によって、未調節の画素強度データから調
節された画素強度データへの１対１の写像（ｍａｐｐｉ
ｎｇ）を作成することにより調節される。この写像はソ
ース像または投影像のいずれかの画素強度データに基づ
いて作成されるが、好ましい実施態様では、第１の（す
なわち、０゜の）投影像データを使用して作成される。
ケーリングされた像平面データが所望の輝度およびコン
トラスト範囲を達成するように写像される（図３のステ
ップ３８）。本発明によれば、投影像のコントラスト
は、主制御器によって、未調節の画素強度データから調
節された画素強度データへの１対１の写像（ｍａｐｐｉ
ｎｇ）を作成することにより調節される。この写像はソ
ース像または投影像のいずれかの画素強度データに基づ
いて作成されるが、好ましい実施態様では、第１の（す
なわち、０゜の）投影像データを使用して作成される。
【００２８】コントラスト調節写像を作成するため、主
制御器は、多数の規定された強度範囲またはビン（ｂｉ
ｎ）の内の各々の中の強度を持つ画素の数のヒストグラ
ムを編集する。図７は、比較的大きい強度の画素を持つ
投影像に対するこのようなヒストグラムを示す。図９
は、比較的小さい強度の画素を持つ投影像に対するこの
ようなヒストグラムを示す。同時に、主制御器は、ソー
スまたは投影像内の最大画素強度を決定する。最大の画
素数を含むビンから開始して、順次画素数の少ない残り
のビンを加えることにより、各々の強度範囲またはビン
内の画素の数が、像内の全画素数の所定のパーセントに
達するまで加算される。計数されたビンが主として比較
的大きい強度の画素を含んでいるとき（図７参照）、計
数値内に含まれている最も小さい強度のビンの下限が画
素強度閾値となる。逆に、計数されたビンが主として比
較的小さい強度の画素を含んでいるとき（図９参照）、
計数値内に含まれている最も大きい強度のビンの上限が
画素強度閾値となる。いずれの場合も、画素強度閾値よ
り大きいまたは小さい強度値の範囲が、例えば、８ビッ
ト強度値の場合は０〜２５５の拡大した表示強度値に線
形に相関するように、１対１の写像が作成される。図７
のヒストグラムから導き出した写像が図８に示されてお
り、また図９のヒストグラムから導き出した写像が図１
０に示されている。
制御器は、多数の規定された強度範囲またはビン（ｂｉ
ｎ）の内の各々の中の強度を持つ画素の数のヒストグラ
ムを編集する。図７は、比較的大きい強度の画素を持つ
投影像に対するこのようなヒストグラムを示す。図９
は、比較的小さい強度の画素を持つ投影像に対するこの
ようなヒストグラムを示す。同時に、主制御器は、ソー
スまたは投影像内の最大画素強度を決定する。最大の画
素数を含むビンから開始して、順次画素数の少ない残り
のビンを加えることにより、各々の強度範囲またはビン
内の画素の数が、像内の全画素数の所定のパーセントに
達するまで加算される。計数されたビンが主として比較
的大きい強度の画素を含んでいるとき（図７参照）、計
数値内に含まれている最も小さい強度のビンの下限が画
素強度閾値となる。逆に、計数されたビンが主として比
較的小さい強度の画素を含んでいるとき（図９参照）、
計数値内に含まれている最も大きい強度のビンの上限が
画素強度閾値となる。いずれの場合も、画素強度閾値よ
り大きいまたは小さい強度値の範囲が、例えば、８ビッ
ト強度値の場合は０〜２５５の拡大した表示強度値に線
形に相関するように、１対１の写像が作成される。図７
のヒストグラムから導き出した写像が図８に示されてお
り、また図９のヒストグラムから導き出した写像が図１
０に示されている。
【００２９】本発明の方法によれば、各々の投影像のコ
ントラストは、主制御器が各々の投影像をシネ・メモリ
に書き込む前に、画素強度ヒストグラムから作成された
写像を使用して自動的に調節される。例えば、図７のグ
ラフの示した大きい強度の投影像の場合、約１６５〜２
５０の範囲内の強度値を持つ投影画素データが、図８に
傾斜した直線で示されるように、０〜２５５の範囲内に
ある表示強度値に写像または変換される。同様に、図９
のグラフの示した小さい強度の投影像の場合、約０〜８
０の範囲内の強度値を持つ投影画素データが、図１０に
傾斜した直線で示されるように、０〜２５５の範囲内に
ある表示強度値に写像または変換される。写像のための
入力範囲は、特定の像から導き出したそれぞれのヒスト
グラムにおける計数されたビンの範囲に応じて変わるこ
とが理解されよう。写像入力範囲の外側の強度値を持つ
画素データは廃棄される。従って最も関心の大きい画素
強度データのコントラストを増大し且つ最も関心の少な
い画素強度データを廃棄することによって、各々の投影
像の画素強度データは、画素閾値より大きい又は小さい
所望の輝度およびコントラスト範囲へ写像される。コン
トラスト調節された画素強度データは次いでその後の表
示のためにシネ・メモリへ戻される。
ントラストは、主制御器が各々の投影像をシネ・メモリ
に書き込む前に、画素強度ヒストグラムから作成された
写像を使用して自動的に調節される。例えば、図７のグ
ラフの示した大きい強度の投影像の場合、約１６５〜２
５０の範囲内の強度値を持つ投影画素データが、図８に
傾斜した直線で示されるように、０〜２５５の範囲内に
ある表示強度値に写像または変換される。同様に、図９
のグラフの示した小さい強度の投影像の場合、約０〜８
０の範囲内の強度値を持つ投影画素データが、図１０に
傾斜した直線で示されるように、０〜２５５の範囲内に
ある表示強度値に写像または変換される。写像のための
入力範囲は、特定の像から導き出したそれぞれのヒスト
グラムにおける計数されたビンの範囲に応じて変わるこ
とが理解されよう。写像入力範囲の外側の強度値を持つ
画素データは廃棄される。従って最も関心の大きい画素
強度データのコントラストを増大し且つ最も関心の少な
い画素強度データを廃棄することによって、各々の投影
像の画素強度データは、画素閾値より大きい又は小さい
所望の輝度およびコントラスト範囲へ写像される。コン
トラスト調節された画素強度データは次いでその後の表
示のためにシネ・メモリへ戻される。
【００３０】上記の投影法は、シネ・メモリから検索さ
れた関心のあるデータ・ボリュームに対して、Ｂモード
強度データとカラー流れ速度またはパワー・データとに
別々に適用される。投影像内の各々の画素は、所与の像
平面上への投影によって導き出された、変換された強度
データおよび変換された速度またはパワー・データを含
む。更に、シネ・メモリがオペレータによって停止（ｆ
ｒｅｅｚｅ）されたときに、ＣＰＵ４２はオプションと
してシネ・メモリ２４の区分２４Ｂ内の多数の相次ぐア
ドレスにＸ−Ｙメモリ１８からの最後のフレームを記憶
させる。第１の投影視角に対する投影像データがシネ・
メモリの区分２４Ｂ内の第１のアドレスに書き込まれ、
これにより関心のある領域内の投影像データが背景のフ
レーム上に重畳される。このプロセスは、全ての投影像
がシネ・メモリの区分２４Ｂ内に記憶されるまで、各々
の角度増分に対して繰り返される。各々の投影像フレー
ムは、変換されたデータを含む関心のある領域で構成さ
れ、オプションとして関心のある領域を囲んでいて、関
心のある領域の変換されたデータによってオーバーライ
ト（上書き）されない背景フレーム・データより成る背
景周辺部を含む。背景の像は、各々の表示された投影が
どの場所から見たものであるかを一層明確にする。そこ
で、オペレータは投影像の内の任意の１つを表示のため
に選択することが出来る。更に、一連の投影像を表示モ
ニタ上にリプレイ（再生）して、被検体ボリュームをあ
たかも観察者の前で回転しているかのように表示するこ
とが出来る。
れた関心のあるデータ・ボリュームに対して、Ｂモード
強度データとカラー流れ速度またはパワー・データとに
別々に適用される。投影像内の各々の画素は、所与の像
平面上への投影によって導き出された、変換された強度
データおよび変換された速度またはパワー・データを含
む。更に、シネ・メモリがオペレータによって停止（ｆ
ｒｅｅｚｅ）されたときに、ＣＰＵ４２はオプションと
してシネ・メモリ２４の区分２４Ｂ内の多数の相次ぐア
ドレスにＸ−Ｙメモリ１８からの最後のフレームを記憶
させる。第１の投影視角に対する投影像データがシネ・
メモリの区分２４Ｂ内の第１のアドレスに書き込まれ、
これにより関心のある領域内の投影像データが背景のフ
レーム上に重畳される。このプロセスは、全ての投影像
がシネ・メモリの区分２４Ｂ内に記憶されるまで、各々
の角度増分に対して繰り返される。各々の投影像フレー
ムは、変換されたデータを含む関心のある領域で構成さ
れ、オプションとして関心のある領域を囲んでいて、関
心のある領域の変換されたデータによってオーバーライ
ト（上書き）されない背景フレーム・データより成る背
景周辺部を含む。背景の像は、各々の表示された投影が
どの場所から見たものであるかを一層明確にする。そこ
で、オペレータは投影像の内の任意の１つを表示のため
に選択することが出来る。更に、一連の投影像を表示モ
ニタ上にリプレイ（再生）して、被検体ボリュームをあ
たかも観察者の前で回転しているかのように表示するこ
とが出来る。
【００３１】本発明の好ましい実施態様では、超音波イ
メージング・システムは複数の異なる投影モードを有す
る。例えば、投影は最大または最小値の画素を含んでい
てよい。或いは、画素データを逆転して、その最大値を
像平面に投影するような、血管のイメージングに有用な
モードを選択してもよい。更に別のモードでは、面のレ
ンダリング（ｒｎｄｅｒｉｎｇ）を行うようにレイ・キ
ャスティング法を用いてもよい。
メージング・システムは複数の異なる投影モードを有す
る。例えば、投影は最大または最小値の画素を含んでい
てよい。或いは、画素データを逆転して、その最大値を
像平面に投影するような、血管のイメージングに有用な
モードを選択してもよい。更に別のモードでは、面のレ
ンダリング（ｒｎｄｅｒｉｎｇ）を行うようにレイ・キ
ャスティング法を用いてもよい。
【００３２】上記の好ましい実施態様は例示の目的で開
示された。超音波イメージングまたはコンピュータ・グ
ラフィックスの分野における当業者には種々の変更およ
び変形を容易になし得よう。このような全ての変更およ
び変形は特許請求の範囲に包含されるものである。
示された。超音波イメージングまたはコンピュータ・グ
ラフィックスの分野における当業者には種々の変更およ
び変形を容易になし得よう。このような全ての変更およ
び変形は特許請求の範囲に包含されるものである。
【図１】実時間超音波イメージング・システム内の主要
な機能のサブシステムを示すブロック図である。
な機能のサブシステムを示すブロック図である。
【図２】本発明の好ましい実施態様に従って強度および
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成する手段を示すブロック図である。
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成する手段を示すブロック図である。
【図３】本発明の好ましい実施態様に従って強度および
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成するためのアルゴリズムのステップ
を示す流れ図である。
速度またはパワー画素データの相次ぐ立体的投影よりな
るフレームを再構成するためのアルゴリズムのステップ
を示す流れ図である。
【図４】従来技術に従ってレイ・キャスティングの逆投
影を立体的にレンダリングする際に関係するサンプリン
グされた関心のある被検体ボリューム、関連のデータ・
ボリュームおよび像投影平面を示す概略線図である。
影を立体的にレンダリングする際に関係するサンプリン
グされた関心のある被検体ボリューム、関連のデータ・
ボリュームおよび像投影平面を示す概略線図である。
【図５】被検体データおよびデータ・ボリュームの同様
なビュー（ｖｉｅｗ）に対応すると共に、三次元超音波
イメージングにおける必要なスケーリング定数を定める
のに有用である一対の幾何学的二次元形状を示す概略線
図である。
なビュー（ｖｉｅｗ）に対応すると共に、三次元超音波
イメージングにおける必要なスケーリング定数を定める
のに有用である一対の幾何学的二次元形状を示す概略線
図である。
【図６】三次元超音波イメージングにおける最大強度投
影を行う手段の概略ブロック図である。
影を行う手段の概略ブロック図である。
【図７】大きい強度の画素を持つ像に対して所定の範囲
内の強度値を持つ画素の集団を示すヒストグラムであ
る。
内の強度値を持つ画素の集団を示すヒストグラムであ
る。
【図８】本発明の好ましい実施態様に従った図７に示し
た大きい強度の画素データのコントラスト調節写像を示
すグラフである。
た大きい強度の画素データのコントラスト調節写像を示
すグラフである。
【図９】小さい強度の画素を持つ像に対して所定の範囲
内の強度値を持つ画素の集団を示すヒストグラムであ
る。
内の強度値を持つ画素の集団を示すヒストグラムであ
る。
【図１０】本発明の好ましい実施態様に従った図９に示
した小さい強度の画素データのコントラスト調節写像を
示すグラフである。
した小さい強度の画素データのコントラスト調節写像を
示すグラフである。
２ ビーム形成装置 ４ 処理装置 ４Ａ Ｂモード処理装置 ４Ｂ カラー流れ処理装置 ６ 走査変換器／表示制御器 ８ 主制御器 １２ 表示モニタ １４Ａ Ｂモード音響線メモリ １４Ｂ カラー流れ音響線メモリ １８ Ｘ−Ｙ表示メモリ ２０ 時間線／グラフィック処理装置及び表示メモリ ２２ ビデオ処理装置 ２４ シネ・メモリ ４２ 中央処理装置 ４４ ランダム・アクセス・メモリ ５０ サンプル（被検体） ５２ サンプル（被検体）ボリューム ５４ データ・ボリューム ５６ 像平面 ５８ 視線 ５８’ 射影線 ６０ 画素 ６２ 射線 ６４ 衝突位置 ７０ 三次元データ・メモリ手段 ７２ アドレス作成手段 ７４ ＣＰＵ ７６ 回転パラメータ計算手段 ７８ 視方向マトリクス手段 ８０ 像平面メモリ手段 ８２ データ比較手段 ８６ 画素遅延素子 ８８ 乗算器 ９０ 加算器 ９２ ライン・メモリ ９４ ライン・メモリ
Claims (20)
- 【請求項１】 被検体ボリューム中の超音波散乱媒体の
三次元イメージング・システムにおいて、 超音波ビームを送信して、被検体ボリューム内の多数の
サンプル・ボリュームから反射された超音波エコーを検
出する超音波トランスジューサ・アレイ、 散乱媒体によって反射された超音波エコーから導き出さ
れた画素強度データを取得する取得手段であって、各々
の画素強度データが前記多数のサンプル・ボリュームの
内のそれぞれ１つのサンプル・ボリュームに対応してい
る取得手段、 前記多数のサンプル・ボリュームの各々に対して前記取
得された画素強度データを記憶するメモリ手段、 前記メモリ手段から、前記被検体ボリューム内の関心の
あるボリュームに対応する一組の画素強度データを検索
する手段、 前記関心のあるボリュームに対応する前記画素強度デー
タの組を像平面へ投影して、投影像を表す投影画素強度
データの組を形成する手段、 前記投影画素強度データの組を変換して、コントラスト
調節された投影像を表す変換手段、 表示モニタ、および前記コントラスト調節された投影像
を前記表示モニタに表示させる手段、を含んでいること
を特徴とする三次元イメージング・システム。 - 【請求項２】 前記変換手段が、第１の強度範囲内にあ
る前記投影画素強度データの組の内の強度値をゼロにす
る手段、および第２の強度範囲内にある前記投影画素強
度データの組の内の強度値を、線形写像に従って、第３
の強度範囲内にあるコントラスト調節された強度値に変
換する手段を含み、 前記第１および第２の強度範囲は画素強度閾値で連続し
ており、前記第３の強度範囲は前記第２の強度範囲より
も広く且つ前記第２の強度範囲にオーバーラップしてお
り、また前記コントラスト調節された強度値は前記コン
トラスト調節された投影像を表す写像された投影画素強
度データの組を形成する請求項１記載のシステム。 - 【請求項３】 前記第２の強度範囲が、前記第１の強度
範囲内の強度値よりも大きい強度値を含んでいる請求項
２記載のシステム。 - 【請求項４】 前記第２の強度範囲が、前記第１の強度
範囲内の強度値よりも小さい強度値を含んでいる請求項
２記載のシステム。 - 【請求項５】 更に、多数の強度サブレンジの各々の強
度サブレンジ内にある強度範囲を持つ前記投影画素強度
データの組の中の画素の数を計数する手段を含み、前記
第１の強度範囲が前記多数の強度サブレンジの内の第１
組の強度サブレンジを含み、前記第２の強度範囲が前記
多数の強度サブレンジの内の第２組の強度サブレンジを
含んでいるが前記第１組の強度サブレンジを含んでいな
い請求項２記載のシステム。 - 【請求項６】 更に、前記の計数値が所定の数に等しい
ことに応答して、前記画素強度閾値を決定する手段、お
よび前記画素強度閾値の関数として前記線形写像を作成
する手段を含んでいる請求項５記載のシステム。 - 【請求項７】 更に、多数の強度サブレンジの各々の強
度サブレンジ内にある強度範囲を持つ前記画素強度デー
タの組の中の画素の数を計数する手段を含み、前記第１
の強度範囲が前記多数の強度サブレンジの内の第１組の
強度サブレンジを含み、前記第２の強度範囲が前記多数
の強度サブレンジの内の第２組の強度サブレンジを含ん
でいるが前記第１組の強度サブレンジを含んでいない請
求項２記載のシステム。 - 【請求項８】 更に、前記の計数値が所定の数に等しい
ことに応答して、前記画素強度閾値を決定する手段、お
よび前記画素強度閾値の関数として前記線形写像を作成
する手段を含んでいる請求項７記載のシステム。 - 【請求項９】 被検体ボリューム中の超音波散乱媒体の
三次元イメージング方法において、 被検体ボリュームと交差する走査平面内に超音波ビーム
を送信するステップ、前記走査平面内にある多数のサン
プル・ボリュームから反射された超音波エコーを検出す
るステップ、 前記被検体ボリュームにわたって前記走査平面を走査す
るステップ、 散乱媒体によって反射された超音波エコーから導き出さ
れた画素強度データを取得するステップであって、各々
の画素強度データが前記多数のサンプル・ボリュームの
内のそれぞれ１つのサンプル・ボリュームに対応してい
るステップ、 前記多数のサンプル・ボリュームの各々に対して前記取
得された画素強度データを記憶するステップ、 前記記憶された画素強度データから、前記被検体ボリュ
ーム内の関心のあるボリュームに対応する一組の画素強
度データを検索するステップ、 前記関心のあるボリュームに対応する前記画素強度デー
タの組を像平面へ投影して、投影像を表す投影画素強度
データの組を形成するステップ、 前記投影画素強度データの組を変換して、コントラスト
調節された投影像を表す変換ステップ、および前記コン
トラスト調節された投影像を表示するステップ、を含ん
でいることを特徴とする三次元イメージング方法。 - 【請求項１０】 前記変換ステップが、第１の強度範囲
内にある前記投影画素強度データの組の内の強度値をゼ
ロにするステップ、および第２の強度範囲内にある前記
投影画素強度データの組の内の強度値を、線形写像に従
って、第３の強度範囲内にあるコントラスト調節された
強度値に変換するステップを含み、 前記第１および第２の強度範囲は画素強度閾値で連続し
ており、前記第３の強度範囲は前記第２の強度範囲より
も広く且つ前記第２の強度範囲にオーバーラップしてお
り、また前記コントラスト調節された強度値は前記コン
トラスト調節された投影像を表す写像された投影画素強
度データの組を形成する請求項９記載の方法。 - 【請求項１１】 前記第２の強度範囲が、前記第１の強
度範囲内の強度値よりも大きい強度値を含んでいる請求
項１０記載の方法。 - 【請求項１２】 前記第２の強度範囲が、前記第１の強
度範囲内の強度値よりも小さい強度値を含んでいる請求
項１０記載の方法。 - 【請求項１３】 更に、多数の強度サブレンジの各々の
強度サブレンジ内にある強度範囲を持つ前記投影画素強
度データの組の中の画素の数を計数するステップを含
み、前記第１の強度範囲が前記多数の強度サブレンジの
内の第１組の強度サブレンジを含み、前記第２の強度範
囲が前記多数の強度サブレンジの内の第２組の強度サブ
レンジを含んでいるが前記第１組の強度サブレンジを含
んでいない請求項１０記載の方法。 - 【請求項１４】 更に、前記の計数値が所定の数に等し
いことに応答して、前記画素強度閾値を決定するステッ
プ、および前記画素強度閾値の関数として前記線形写像
を作成するステップを含んでいる請求項１３記載の方
法。 - 【請求項１５】 更に、多数の強度サブレンジの各々の
強度サブレンジ内にある強度範囲を持つ前記関心のある
ボリュームに対応する前記画素強度データの組のサブセ
ット内の画素の数を計数するステップを含み、前記第１
の強度範囲が前記多数の強度サブレンジの内の第１組の
強度サブレンジを含み、前記第２の強度範囲が前記多数
の強度サブレンジの内の第２組の強度サブレンジを含ん
でいるが前記第１組の強度サブレンジを含んでいない請
求項１０記載の方法。 - 【請求項１６】 更に、前記の計数値が所定の数に等し
いことに応答して、前記画素強度閾値を決定するステッ
プ、および前記画素強度閾値の関数として前記線形写像
を作成するステップを含んでいる請求項１５記載の方
法。 - 【請求項１７】 被検体ボリューム中の超音波散乱媒体
の三次元イメージング方法において、 被検体ボリュームと交差する走査平面内に超音波ビーム
を送信するステップ、前記走査平面内にある多数のサン
プル・ボリュームから反射された超音波エコーを検出す
るステップ、 前記被検体ボリュームにわたって多数の走査位置で前記
走査平面を走査するステップ、 前記走査平面の多数の走査位置によって定められたそれ
ぞれの平面内の散乱媒体によって反射された超音波エコ
ーから導き出された多数のサブセットの画素強度データ
を取得するステップであって、画素強度データの各々の
サブセットがそれぞれ１つの強度データ・フレームを形
成するステップ、 前記の多数の強度データ・フレームを記憶するステッ
プ、 選択された相次ぐ記憶された強度データ・フレームか
ら、関心のある領域に対応するそれぞれのサブセットの
画素強度データを検索するステップであって、これらの
画素強度データのサブセットが被検体ボリューム内の関
心のあるボリュームに対応する画素強度データの組を形
成するステップ、 前記画素強度データの組を像平面へ投影して、投影像を
表す投影画素強度データの組を形成するステップ、 前記投影画素強度データの組を変換して、コントラスト
調節された投影像を表す変換ステップ、および前記コン
トラスト調節された投影像を前記多数の強度データ・フ
レームの選択された１つの中心領域上に重畳して成る像
を表示するステップ、 を含んでいることを特徴とする三次元イメージング方
法。 - 【請求項１８】 前記変換ステップが、第１の強度範囲
内にある前記投影画素強度データの組の内の強度値をゼ
ロにするステップ、および第２の強度範囲内にある前記
投影画素強度データの組の内の強度値を、線形写像に従
って、第３の強度範囲内にあるコントラスト調節された
強度値に変換するステップを含み、 前記第１および第２の強度範囲は画素強度閾値で連続し
ており、前記第３の強度範囲は前記第２の強度範囲より
も広く且つ前記第２の強度範囲にオーバーラップしてお
り、また前記コントラスト調節された強度値は前記コン
トラスト調節された投影像を表す写像された投影画素強
度データの組を形成する請求項１７記載の方法。 - 【請求項１９】 更に、多数の強度サブレンジの各々の
強度サブレンジ内にある強度範囲を持つ前記投影画素強
度データの組の中の画素の数を計数するステップを含
み、前記第１の強度範囲が前記多数の強度サブレンジの
内の第１組の強度サブレンジを含み、前記第２の強度範
囲は前記多数の強度サブレンジの内の第２組の強度サブ
レンジを含んでいるが前記第１組の強度サブレンジを含
んでいない請求項１８記載の方法。 - 【請求項２０】 更に、前記の計数値が所定の数に等し
いことに応答して、前記画素強度閾値を決定するステッ
プ、および前記画素強度閾値の関数として前記線形写像
を作成するステップを含んでいる請求項１９記載の方
法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/852,772 US5954653A (en) | 1997-05-07 | 1997-05-07 | Method and apparatus for automatically enhancing contrast in projected ultrasound image |
| US08/852772 | 1997-05-07 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1147134A true JPH1147134A (ja) | 1999-02-23 |
| JP4155618B2 JP4155618B2 (ja) | 2008-09-24 |
Family
ID=25314174
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP12345098A Expired - Fee Related JP4155618B2 (ja) | 1997-05-07 | 1998-05-06 | 超音波散乱媒体の三次元イメージング・システムおよび方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5954653A (ja) |
| JP (1) | JP4155618B2 (ja) |
| DE (1) | DE19819800A1 (ja) |
| IL (1) | IL124164A (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000333956A (ja) * | 1999-05-31 | 2000-12-05 | Toshiba Corp | 超音波診断装置及び超音波画像の表示切換方法 |
| JP2001043396A (ja) * | 1999-07-15 | 2001-02-16 | Mitsubishi Electric Inf Technol Center America Inc | ボリュームデータセットのヒストグラムを生成する方法及びその装置 |
| JP2007054504A (ja) * | 2005-08-26 | 2007-03-08 | Hitachi Medical Corp | 超音波診断装置 |
| JP2009153925A (ja) * | 2007-12-27 | 2009-07-16 | Toshiba Corp | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
Families Citing this family (42)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11164833A (ja) * | 1997-09-30 | 1999-06-22 | Toshiba Corp | 医用画像診断装置 |
| JP3218216B2 (ja) * | 1998-03-20 | 2001-10-15 | アロカ株式会社 | 三次元画像処理装置 |
| DE69908713T2 (de) * | 1998-06-22 | 2004-05-13 | Texas Instruments Inc., Dallas | Selektive Helligkeitserweiterung |
| US6734904B1 (en) * | 1998-07-23 | 2004-05-11 | Iteris, Inc. | Imaging system and method with dynamic brightness control |
| US6048313A (en) * | 1998-09-17 | 2000-04-11 | Stonger; Kelly A. | Method and apparatus for fractal-based enhancement of ultrasound imaging |
| US6059729A (en) * | 1998-10-19 | 2000-05-09 | Stonger; Kelly A. | Method and apparatus for edge enhancement in ultrasound imaging |
| US6017309A (en) * | 1998-12-31 | 2000-01-25 | Washburn; Michael J. | Ultrasound color flow display optimization by adjusting color maps |
| US6102861A (en) * | 1999-04-23 | 2000-08-15 | General Electric Company | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering |
| US6413219B1 (en) | 1999-03-31 | 2002-07-02 | General Electric Company | Three-dimensional ultrasound data display using multiple cut planes |
| US6434260B1 (en) * | 1999-07-12 | 2002-08-13 | Biomedicom, Creative Biomedical Computing Ltd. | Facial imaging in utero |
| US6542626B1 (en) * | 1999-11-05 | 2003-04-01 | General Electric Company | Method and apparatus for adapting imaging system operation based on pixel intensity histogram |
| US6579238B1 (en) | 2000-04-24 | 2003-06-17 | Acuson Corporation | Medical ultrasonic imaging system with adaptive multi-dimensional back-end mapping |
| US6398733B1 (en) | 2000-04-24 | 2002-06-04 | Acuson Corporation | Medical ultrasonic imaging system with adaptive multi-dimensional back-end mapping |
| FI20010563A (fi) * | 2001-03-20 | 2002-09-21 | Runtech Systems Oy | Laitteisto ja menetelmä telapintojen ja/tai kudosten käsittelemiseksi |
| 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 |
| US7092558B2 (en) * | 2002-08-14 | 2006-08-15 | General Electric Company | Automated optimization of medical 3D visualizations |
| US20050053305A1 (en) * | 2003-09-10 | 2005-03-10 | Yadong Li | Systems and methods for implementing a speckle reduction filter |
| US7275547B2 (en) * | 2003-10-08 | 2007-10-02 | Boston Scientific Scimed, Inc. | Method and system for determining the location of a medical probe using a reference transducer array |
| US6979295B2 (en) | 2003-11-19 | 2005-12-27 | Ge Medical Systems Global Technology Company, Llc | Automatic color gain adjustments |
| US20050110793A1 (en) * | 2003-11-21 | 2005-05-26 | Steen Erik N. | Methods and systems for graphics processing in a medical imaging system |
| JP4564286B2 (ja) * | 2004-06-14 | 2010-10-20 | 株式会社東芝 | ３次元超音波画像化装置 |
| US7627386B2 (en) * | 2004-10-07 | 2009-12-01 | Zonaire Medical Systems, Inc. | Ultrasound imaging system parameter optimization via fuzzy logic |
| US8834371B2 (en) * | 2005-03-31 | 2014-09-16 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus and ultrasound image processing program |
| US8002705B1 (en) | 2005-07-22 | 2011-08-23 | Zonaire Medical Systems, Inc. | Continuous transmit focusing method and apparatus for ultrasound imaging system |
| US8784318B1 (en) | 2005-07-22 | 2014-07-22 | Zonare Medical Systems, Inc. | Aberration correction using channel data in ultrasound imaging system |
| EP1750143A3 (en) | 2005-07-27 | 2009-09-09 | Medison Co., Ltd. | Ultrasound diagnostic system and method of automatically controlling brightness and contrast of a three-dimensional ultrasound image |
| JP5226181B2 (ja) * | 2005-11-24 | 2013-07-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像診断装置 |
| US7796830B2 (en) * | 2006-08-15 | 2010-09-14 | Nokia Corporation | Adaptive contrast optimization of digital color images |
| JP4881112B2 (ja) * | 2006-09-19 | 2012-02-22 | 株式会社東芝 | 超音波診断装置及び画像データ生成方法 |
| CN101523237B (zh) * | 2006-10-13 | 2015-01-14 | 皇家飞利浦电子股份有限公司 | 采用灰度级反转的3d超声彩色流成像 |
| US20080146922A1 (en) * | 2006-10-24 | 2008-06-19 | Zonare Medical Systems, Inc. | Control of user interfaces and displays for portable ultrasound unit and docking station |
| US20080242988A1 (en) * | 2007-03-28 | 2008-10-02 | Tetsuya Yoshida | Ultrasound diagnostic apparatus, ultrasound image processing apparatus, and ultrasound image processing method |
| US9060669B1 (en) | 2007-12-20 | 2015-06-23 | Zonare Medical Systems, Inc. | System and method for providing variable ultrasound array processing in a post-storage mode |
| US20100036727A1 (en) * | 2008-08-07 | 2010-02-11 | Sarelson Seth H | Method of Tracking the Impact of Paid Search on Offline Sales |
| DE102009058887A1 (de) * | 2009-12-18 | 2011-06-22 | Daimler AG, 70327 | Anzeigeverfahren und Anzeigevorrichtung für ein Fahrzeug und betreffendes Computerprogrammprodukt |
| KR102185726B1 (ko) * | 2014-01-28 | 2020-12-02 | 삼성메디슨 주식회사 | 관심 영역에 대응하는 초음파 영상 표시 방법 및 이를 위한 초음파 장치 |
| US9842424B2 (en) * | 2014-02-10 | 2017-12-12 | Pixar | Volume rendering using adaptive buckets |
| WO2019038097A1 (en) * | 2017-08-23 | 2019-02-28 | Koninklijke Philips N.V. | METHOD FOR DETECTING FETAL MOTION FROM DOPPLER ULTRASONIC SIGNAL |
| EP3466344A1 (en) * | 2017-10-03 | 2019-04-10 | Koninklijke Philips N.V. | A fetal movement detection method from doppler ultrasound signal |
| US11771399B2 (en) | 2018-02-07 | 2023-10-03 | Atherosys, Inc. | Apparatus and method to guide ultrasound acquisition of the peripheral arteries in the transverse plane |
| EP3820375A1 (en) * | 2018-07-11 | 2021-05-19 | Koninklijke Philips N.V. | Ultrasound imaging system with pixel extrapolation image enhancement |
Family Cites Families (13)
| 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 |
| GB9025431D0 (en) * | 1990-11-22 | 1991-01-09 | Advanced Tech Lab | Three dimensional ultrasonic imaging |
| US5469849A (en) * | 1993-06-14 | 1995-11-28 | Kabushiki Kaisha Toshiba | Ultrasound diagnosis apparatus |
| US5456255A (en) * | 1993-07-12 | 1995-10-10 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus |
| 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 |
| JP3406106B2 (ja) * | 1995-02-06 | 2003-05-12 | ジーイー横河メディカルシステム株式会社 | 超音波画像表示方法および超音波診断装置 |
| 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 |
| EP0754957B1 (en) * | 1995-07-17 | 2001-03-21 | Aloka Co. Ltd. | Ultrasound image processing apparatus and method for producing an ultrasound image |
| 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 |
| JP3407169B2 (ja) * | 1995-10-12 | 2003-05-19 | 富士写真光機株式会社 | 超音波画像立体表示装置及び超音波画像立体表示方法 |
-
1997
- 1997-05-07 US US08/852,772 patent/US5954653A/en not_active Expired - Lifetime
-
1998
- 1998-04-22 IL IL12416498A patent/IL124164A/xx not_active IP Right Cessation
- 1998-05-04 DE DE19819800A patent/DE19819800A1/de not_active Ceased
- 1998-05-06 JP JP12345098A patent/JP4155618B2/ja not_active Expired - Fee Related
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000333956A (ja) * | 1999-05-31 | 2000-12-05 | Toshiba Corp | 超音波診断装置及び超音波画像の表示切換方法 |
| JP2001043396A (ja) * | 1999-07-15 | 2001-02-16 | Mitsubishi Electric Inf Technol Center America Inc | ボリュームデータセットのヒストグラムを生成する方法及びその装置 |
| JP2007054504A (ja) * | 2005-08-26 | 2007-03-08 | Hitachi Medical Corp | 超音波診断装置 |
| JP2009153925A (ja) * | 2007-12-27 | 2009-07-16 | Toshiba Corp | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| IL124164A (en) | 2000-06-29 |
| JP4155618B2 (ja) | 2008-09-24 |
| US5954653A (en) | 1999-09-21 |
| DE19819800A1 (de) | 1998-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4155618B2 (ja) | 超音波散乱媒体の三次元イメージング・システムおよび方法 | |
| US5779641A (en) | Method and apparatus for three-dimensional ultrasound imaging by projecting filtered pixel data | |
| US5904653A (en) | Method and apparatus for three-dimensional ultrasound imaging combining intensity data with color flow velocity or power data | |
| JP4204095B2 (ja) | 被検体ボリュームの三次元イメージング・システムおよび方法 | |
| JP4424799B2 (ja) | ３次元イメージング方法及びシステム | |
| 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 | |
| US6413219B1 (en) | Three-dimensional ultrasound data display using multiple cut planes | |
| US5546807A (en) | High speed volumetric ultrasound imaging system | |
| US5329929A (en) | Ultrasonic diagnostic apparatus | |
| US5934288A (en) | Method and apparatus for displaying 3D ultrasound data using three modes of operation | |
| JP4628645B2 (ja) | 複数の２ｄスライスから画像を作成するための超音波診断方法及び装置 | |
| US20040081340A1 (en) | Image processing apparatus and ultrasound diagnosis apparatus | |
| US20050267366A1 (en) | Ultrasonic diagnostic apparatus and image processing method | |
| US6126603A (en) | Method and apparatus for segmenting color flow mode data using velocity information in three-dimensional ultrasound imaging | |
| JP2000152939A (ja) | 物体ボリュ―ムの３次元イメ―ジングのためのシステム及び方法 | |
| US20120053464A1 (en) | Ultrasonic image processing apparatus | |
| JP2003061956A (ja) | 超音波診断装置、医用診断装置及び画像処理方法 | |
| US7302092B1 (en) | Three-dimensional imaging system | |
| US6458082B1 (en) | System and method for the display of ultrasound data | |
| JP3936450B2 (ja) | 投影画像生成装置及び医用画像装置 | |
| US6306092B1 (en) | Method and apparatus for calibrating rotational offsets in ultrasound transducer scans | |
| US6482159B1 (en) | Method for the examination of objects with ultrasound | |
| EP0912905B1 (en) | Three-dimensional imaging system | |
| JP2001515373A (ja) | 高速三次元超音波画像形成システム |
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: 20071106 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080226 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080321 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20080507 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080610 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080708 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110718 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110718 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120718 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120718 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130718 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |