JP2883584B2 - 超音波画像処理装置及び方法 - Google Patents
超音波画像処理装置及び方法Info
- Publication number
- JP2883584B2 JP2883584B2 JP8185781A JP18578196A JP2883584B2 JP 2883584 B2 JP2883584 B2 JP 2883584B2 JP 8185781 A JP8185781 A JP 8185781A JP 18578196 A JP18578196 A JP 18578196A JP 2883584 B2 JP2883584 B2 JP 2883584B2
- Authority
- JP
- Japan
- Prior art keywords
- voxel
- image processing
- value
- opacity
- ultrasonic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
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/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/8918—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being linear
-
- 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/8934—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration
- G01S15/8945—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration using transducers mounted for linear mechanical movement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Image Analysis (AREA)
Description
法に関し、特に、立体画像と投影画像の両者の性質を合
わせもったような新しい三次元画像(特に、三次元超音
波画像)を形成する装置及び方法に関し、またそのよう
な三次元画像の処理を迅速に行うことができる装置及び
方法に関する。
モード表示では、生体内の組織の断面が白黒の濃淡画像
として表示される。レントゲン像が透視像であるのに対
し、この超音波断層画像は単純な断面の像であり、空間
的な把握が困難であるという面がある。すなわち、レン
トゲン像では、皮膚を通して患者の体内を透かして見た
ような画像が得られ、例えば、観察対象となった臓器の
全体構造やその前後の組織との関係などを把握・認識し
ながら、当該臓器の観察を行えるという利点がある。一
方、単なる超音波断層画像では、超音波ビームが走査さ
れた走査面のみの情報しか画像化されず、その臓器の全
体構造やその前後の他の組織との関係などは画像として
現れない。それゆえ、他の部位を観察したければ超音波
探触子を移動・傾斜させて、走査面自体を変更させる必
要があるが、いずれにしても三次元的に(空間的に)組
織を把握するのは困難である。
移動させて三次元領域内でエコーデータを取り込み、画
像処理により三次元画像を形成することが提案されてい
る。しかし、この従来方法では、三次元領域内のエコー
データをいったん三次元エコーメモリに格納した上で所
定の順序でデータを読み出して処理を行うことが必要と
され、画像処理のためにはデータ処理に時間がかかり、
リアルタイムで画像を形成することが全く困難であっ
た。また、表面抽出により臓器の表面を立体的に濃淡表
示するものであるために、臓器表面と臓器内部の情報を
併せて画像化するようなことはできなかった。従来、エ
コーデータを積算して積算画像を形成することも提案さ
れているが、その積算画像は奥行き感のないものであ
る。
おいては、三次元領域内の各エコーデータがいったん三
次元エコーデータメモリに書き込まれた後に、読み出さ
れて画像処理されており、そのようなデータの書き込み
・読み出しのために、データの取り込みと同期して画像
処理を行うことができなかった。また、従来において
は、組織表面と内部を同時に画像化することや所望の透
明度をもって組織内部を立体的に観察することはでき
ず、このため新しい画像処理法が要望されていた。
超音波診断装置に限られず、超音波などの測定波を三次
元領域に送受波し、その三次元領域全体又はその領域内
の物体を画像化する他の装置でも指摘されている。
ものであり、その目的は、奥行き感(立体感)が得ら
れ、同時にレントゲン像のような性質(透過性)をもっ
た画像を形成できる画像処理装置及び方法を提供するこ
とにある。
迅速に形成できる画像処理装置及び方法を提供すること
にある。
奥行き感と透明度のバランスを自在に調整できる画像処
理装置及び方法を提供することにある。
知られているボリューム・レンダリング(Volume Rende
ring)法を改良・発展させ、三次元計測に適合する新し
い画像処理法を提供することにある。
三次元領域内において処理範囲を適切に設定でき、目的
とする物体を的確に画像化できる画像処理装置及び方法
を提供することにある。
臓器をあたかも光に照らされているかのように表現でき
る超音波画像処理装置及び方法を提供することにある。
処理で必要とされていた三次元データメモリを排除で
き、迅速な処理を行える画像処理装置及び方法を提供す
ることにある。
ームに沿って開始ボクセルから終了ボクセルまでボクセ
ル処理を順次行うことにより、その超音波ビームに対応
する画素の画素値を求める超音波画像処理装置であっ
て、エコー値ei に基づきボクセルi の不透明度αi を
演算する不透明度演算手段と、エコー値eiに基づきボ
クセルi の透明度βi を演算する透明度演算手段と、エ
コー値ei に不透明度αi を乗算し、ボクセルi の発光
量を演算する発光量演算手段と、1つ前のボクセルi-1
の出力光量にボクセルi の透明度βi を乗算し、ボクセ
ルi の透過光量を演算する透過光量演算手段と、前記発
光量と前記透過光量とを加算し、ボクセルi の出力光量
を求める光量加算手段と、を含み、終了ボクセルの出力
光量を画素値に対応させて超音波画像を形成することを
特徴とする。
ビームが三次元的に走査され、その三次元領域内の各位
置でエコー値ei (すなわちエコーデータ)が取り込ま
れる。ここで、その三次元領域は、多数の「ボクセル」
の集合体として仮定される。このモデルでは、各ボクセ
ルは標本点(サンプル・ボリューム)に相当し、各ボク
セルは「ボクセル値」としてエコー値ei を有する。ま
た、各ボクセルには、後述のように「不透明度」及び
「透明度」が定義される。
素に対応する仮想的な透視線が複数設定される。そし
て、各透視線ごとにボクセル処理がなされ、すなわち、
所定の開始ボクセルから所定の終了ボクセルまで、各ボ
クセルのエコー値ei が逐次的にボクセル処理され、最
終的な処理結果(終了ボクセルの出力光量)が画素の画
素値(例えば、輝度値、色相など)に変換される。
ビーム方向に設定され、換言すれば、超音波ビームに沿
ってボクセル処理が実行される。これにより、順次取り
込まれるエコーデータを時系列順で逐次的にリアルタイ
ム処理でき、また、従来装置において必要であった三次
元データメモリを不要にすることができる。すなわち、
本発明によれば、取り込まれたエコーデータはその取り
込み順序で処理され、三次元データメモリにいったんす
べてのエコーデータを格納させなくても、データ処理を
十分に行える。なお、超音波ビームと透視線とを一致さ
せない場合には、三次元領域の外部に視点を設定し、そ
の視点から複数の透視線を放射状に設定することにな
る。
クセル(通常は、原点に最も近い最初のボクセル)から
所定の終了条件が満たされる終了ボクセルまで順次行わ
れる。なお、透視線が超音波ビームとして定義される場
合には、開始ボクセルは、通常、最も送受波器(超音波
探触子など)に近いボクセルとなる。ただし、ノイズの
影響を避けるために、送受波器近傍のボクセルについて
は処理対象としないことも可能であり、この場合には所
定の深さのボクセルを開始ボクセルとし、そこからボク
セル処理を行う。
まず、ボクセルi のエコー値ei に基づき不透明度αi
と透明度βi が定義される。そして、エコー値ei に不
透明度αi が乗算されてボクセルiの発光量が演算さ
れ、また、1つ前のボクセルi-1 の出力光量にボクセル
i の透明度βi が乗算されてボクセルi の透過光量が演
算される。
いての周囲への超音波の拡散・散乱の度合いに関わるも
ので、発光量は、ボクセルi の音源(光源)としての強
さを表すものと思われる。一方、透明度βi は、超音波
の透過率に関わるもので、透過光量は、ボクセルi を伝
達媒体として見た場合にその伝達率に相当するものと思
われる。このような発光量と透過光量とが加算されてボ
クセルi の出力光量が演算される。ここで、出力光量は
ボクセルi の画素値への寄与度を表すものである。この
出力光量は、次のボクセルのボクセル処理(透過光量の
演算)に引き渡される。
すると、その最終ボクセルの出力光量が画素値に変換さ
れる。そして、各透視線の画素値が求まれば、それらの
画素値の集合として1枚の超音波画像が形成される。
と立体画像としての性格とを併せて有することが実験に
より確認されている。すなわち、生体内の組織をレント
ゲン写真のように透かして表現でき、その一方、超音波
三次元画像のような奥行き感をもって表現できる。よっ
て、例えば胎児の表面と内部の同時観察などを行うこと
ができ、疾病診断に当たって組織の三次元的な把握を容
易に行うことができる。
化させることによって、所望の画質の超音波画像を構成
でき、例えば透明感を強調したり、または立体感を強調
したりすることができる。あるいは、組織表面を強調し
たり、または組織内部を強調することができる。このよ
うな設定は、実験結果に基づいて自動的に行うことがで
き、あるいは画像表示を確認しながらオペレータがつま
みをリアルタイムで操作することにより実現することも
できる。
明度αi は、0以上1以下の範囲内で設定され、前記透
明度βi は(1−αi )で設定される。
明度αi は、エコー値ei をパラメータとする関数によ
り定義される。ここで、前記関数はエコー値ei の増大
に対して不透明度αi が増大する非線形関数であること
が望ましい。
に含まれる係数を可変させる係数調整手段を含む。かか
る構成によれば、例えば、表示させる深さを調節するこ
とができ、また表面画像から透過画像まで自在に調整で
きる。
ボクセルを判定する終了判定手段を含む。ここで、前記
終了判定手段は、望ましくは、設定された深さまでボク
セル処理が到達した場合に前記終了ボクセルを判定す
る。また、前記判定終了手段は、望ましくは、開始ボク
セルから各ボクセルごとにその不透明度αi を加算し、
その加算値が所定の終了判定値となった場合に終了ボク
セルを判定する。ここで、前記終了判定値は望ましくは
1である。
予め設定された深さのボクセル(あるいは最後のボクセ
ル)に到達した場合や、逐次加算した不透明度αi が終
了判定値(通常は1)に達した場合に終了を判定するも
のである。よって、逐次加算される各不透明度αi の値
が透視線全体にわたって大きければ、比較的早い段階で
処理が終了することになり、例えば組織の表面まで透視
して画像表現が終了することになる。組織表面は比較的
エコー値が大きいので、終了判定条件を適宜設定すれば
組織表面でボクセル処理を終了させて、組織の表面画像
を形成できる。すなわち、しきい値法などを利用するこ
と無く表面抽出を行える。
ボクセルを設定する手段を有する。ここで、開始ボクセ
ルは例えば透視線(超音波ビーム)上の最初のボクセル
であり、あるいは、人為的に又は自動的に設定された深
さのボクセルである。
明度演算手段は、各エコー値ei に対する不透明度αi
を対応付けたテーブルで構成され、前記透明度演算手段
は、各エコー値ei に対する透明度βi を対応付けたテ
ーブルで構成される。このようなテーブルを利用するこ
とにより処理時間を短縮できる。
波の送受波により得られた受信信号を対数増幅する対数
増幅器を有する。実験によれば、このような対数増幅後
にボクセル処理を行うことにより、より良好な超音波画
像を形成できることが確認されている。
本発明は、三次元領域に対して測定波を送受波し、三次
元領域を構成する各ボクセルごとにボクセル値ei を取
り込む送受波手段と、前記三次元領域を通過する仮想的
な透視線を設定し、その透視線に沿って開始ボクセルか
ら終了ボクセルまでボクセル処理を順次行うことによ
り、その透視線に対応する画素の画素値を求める画像処
理手段と、を含む装置であって、前記画像処理手段は、
ボクセル値ei に基づきボクセルi の不透明度αiを演
算する不透明度演算手段と、ボクセル値ei に基づきボ
クセルi の透明度βi を演算する透明度演算手段と、ボ
クセル値ei と不透明度αi とに基づいて、ボクセルi
の発光量を演算する発光量演算手段と、1つ前のボクセ
ルi-1 の出力光量とボクセルi の透明度βi とに基づい
て、ボクセルi の透過光量を演算する透過光量演算手段
と、前記発光量と前記透過光量とに基づいて、ボクセル
i の出力光量を求める出力光量演算手段と、を含み、終
了ボクセルの出力光量を画素値に対応させて前記三次元
領域の画像を形成することを特徴とする。
明は、生体の超音波診断を行う場合の他に、生体以外の
物体の超音波診断を行う場合(例えば金属探傷、ソナ
ー)や魚群探知機にも適用できる。
に、本発明は、三次元領域を通過する仮想的な透視線上
で、開始ボクセルから終了ボクセルまでボクセル処理を
順次行って、透視線に対応する画素の画素値を求める画
像処理方法において、ボクセル値ei によりボクセルi
の不透明度αi (0≦αi ≦1)を演算する工程と、ボ
クセルi の透明度(1−αi )を演算する工程と、ボク
セル値ei に不透明度αiを乗算し、ボクセルi の発光
量を演算する工程と、1つ前のボクセルi-1 の出力光量
にボクセルiの透明度(1−αi )を乗算し、ボクセル
i の透過光量を演算する工程と、前記発光量と前記透過
光量とを加算し、ボクセルi の出力光量を求める工程
と、を含み、終了ボクセルの出力光量を画素値とするこ
とを特徴とする。ここで、望ましくは、前記不透明度α
i は、k1及びk2を係数として、αi =k1・ei k2により
演算される。
に、本発明は、三次元領域を通過する仮想的な透視線上
で、開始ボクセルから終了ボクセルまでボクセル処理を
順次行うことにより、透視線に対応する画素の画素値を
求める超音波画像処理装置において、ボクセル値ei に
基づきボクセルi の不透明度αi を演算する不透明度演
算手段と、前記ボクセル値ei 、前記不透明度αi 、及
び、1つ前のボクセルi-1 の出力光量に相当する入力光
量CINi に基づいて、ボクセルi の出力光量COU Tiを演
算する出力光量演算手段と、を含み、終了ボクセルの出
力光量を画素値に対応させて画像を形成することを特徴
とする。
音波ビームに相当する。また、本発明の好適な態様で
は、前記出力光量演算手段は、CINi +αi ・(ei −
CINi)の演算により、出力光量COUTiを演算する。
に、本発明は、超音波ビームを三次元領域で走査するこ
とにより得られたエコー値に基づき前記三次元領域を画
像化した超音波画像を形成する画像処理手段を含む超音
波画像処理装置において、前記画像処理手段は、前記超
音波ビームに沿った時系列順で入力される各エコー値に
対して逐次処理を行って、当該超音波ビーム上の複数の
エコー値が反映された画素値を演算し、各超音波ビーム
ごとに前記画素値が演算され、それらの画素値の集合と
して前記三次元領域を画像化した超音波画像が形成され
ることを特徴とする。
込み順序で逐次的に処理され、これにより超音波ビーム
ごとに画素値が演算される。すなわち、この構成によれ
ば、各エコー値を三次元データメモリにいったん格納す
ることなく、その時系列順で処理できるので、超音波の
送受波と同期したデータ処理を実現できる。従来から、
三次元領域を画像化した超音波画像(積算画像、投影画
像を含む)の実時間処理が強く要望されていたが、本発
明によれば、以上のような簡単な構成によってその要望
に応えることができる。
次処理では、前のエコー値の演算結果が次のエコー値の
演算で利用され、終了点のエコー値に対する演算結果と
して画素値が決定される。
ームに沿って各エコー値を参照することにより、各超音
波ビームごとに前記終了点を判定する終了判定手段を含
む。また、本発明の好適な態様では、前記終了判定手段
が有する終了判定条件を設定するための終了条件設定手
段を有する。
ームの走査により取り込まれた各エコー値が三次元デー
タメモリを介することなく前記画像処理手段に順次入力
され、前記画像処理手段は各エコー値の入力タイミング
に同期して前記逐次処理を行う。
波画像は、前記三次元領域を平面上に投影した画像であ
る。ここで、その超音波画像には、上記の非透明度など
を利用したボクセル処理による立体投影画像、エコー値
を積算することにより形成される積算画像、エコー値を
投影面に投影することにより形成される投影画像、など
が含まれる。
リューム・レンダリング(Volume Rendering)法を基礎と
し、リアルタイムの画像処理(特に、超音波画像処理)
にその手法を発展させたものである。その際には、特有
の条件が加味されている。そこで、まず本実施形態に係
る画像処理法の原理について、図1〜図3を用いて説明
する。
音波ビームがX方向に走査されると、走査面10が形成
される。この走査面10をZ方向に移動させると、周知
のように三次元エコーデータ取込み空間12が形成され
る。この三次元エコーデータ取込み空間12に対して、
各超音波ビームに沿って本実施形態に係るボクセル処理
を行い、投影面16上に三次元エコーデータ取込み空間
12を投影したものが、図1(B)の超音波画像100
である。超音波画像100では、そのX方向の1ライン
100aが1つの走査面10に相当する。換言すれば、
超音波ビーム(透視線)1本が超音波画像100内の1
画素に相当する。
列順でそのエコーデータに対して以下に詳述するボクセ
ル処理が行われるので、各エコーデータを三次元エコー
データメモリにいったん蓄積して画像形成に必要な順序
でエコーデータを読み出す必要はなく、データ取り込み
と同期したデータ処理が可能となる。
概念が示されている。1つのボクセルは、受信信号をA
/D変換して得られた1つのエコーデータに相当し、換
言すれば、そのA/D変換レートの1周期に相当するボ
リューム(標本点)に相当するものである。すなわち、
超音波ビームは、多数のボクセルの集合体として仮定さ
れる。図2には各ボクセルがi−1からLLASTまで示さ
れている。最初のボクセルから順次処理を行って得られ
た値が超音波画像を構成する1画素の輝度値P(x,
y)に対応する。
透明度β[本実施形態ではβ=(1−α)]を定義する
ことにする。不透明度αは、図3に示すようにボクセル
の周囲への自発的な発光に相当するものである。透明度
(1−α)は1つ前のボクセルからの光に対する当該ボ
クセル中の透過度合いに相当するものである。不透明度
αは0≦α≦1の範囲に設定され、本実施形態におい
て、その不透明度はエコーデータ(エコー値)の関数と
して定義される。具体的には、例えば、
り、またk1は定数(係数)である。k2としては望ましく
は1よりも大きい数値が代入され、例えばk2=2又は3
である。すなわち、エコーデータの値eに対してαは非
線形に変化する。なお、定数k1は可変できるように構成
するのが望ましい。
は、入力光量CINi と出力光量COU Tiとが定義され、そ
の入力光量CINi は1つ前のボクセルi−1の出力光量
COU Ti-1に等しい。すなわち、
ボクセルにおいてはCIN 1 =0である。
度(1−α)に基づいて、発光量と透過光量が定義され
る。すなわち、ボクセルiの発光量は、不透明度とエコ
ーデータの積として定義され、αi ・ei である。ボク
セルiの透過光量は透明度と入力光量の積として定義さ
れ、(1−αi )・CINi である。
その発光量と透過光量は以下のように加算され、当該ボ
クセルの出力光量COUTiが決定される。
わち、1つ前のボクセルでの計算結果が次のボクセルの
計算に利用される。
セルへ、そして、その次のボクセルへと順次行っていく
間において、各ボクセルの不透明度αi を加算し、その
加算値Σαi が1に到達した時点で、処理を終了させ
る。ただし、処理が最後又は設定された深さのボクセル
LLASTとなった場合にも処理を終了させる。すなわち、
処理の終了条件は、
が1に到達した時点で処理を停止させることを意味し、
もちろん、条件に応じて上記第4式の条件、特にαi の
最大加算値(終了判定値)を変更させてもよい。
ル(最終ボクセル)の出力光量COU T が、対応する画素
の輝度P(x,y)として利用される。そして、このよ
うな超音波ビームごとの画素値演算がすべての超音波ビ
ームについて行われると、超音波画像を構成するすべて
の画素の画素値を得られる。すなわち、1枚の超音波画
像が形成される。
(x,y)には、開始ボクセルから終了ボクセルまでの
すべてのエコーデータの値が反映されている。しかし、
それは従来のように単なる単純積算でなく、各ボクセル
での超音波の散乱と吸収の両方を反映したものとなって
いる。よって、あたかも光源から光が出て、各ボクセル
での散乱及び吸収を経て透過した光によって形成される
像のような奥行き感(立体感)と透明感の両者の性質を
もった超音波画像を構成できる。
が(1−αi )で定義され、すなわち不透明度αi によ
って透明度を表すことができるので、演算式中から透明
度の概念を見掛け上消去することができる。よって、以
下のように第3式を式変形することにより、同じ原理に
基づいて、出力光量COUTiを演算できる。
2項をΔi で置き換えると、第5−2式が得られる。す
なわち、ボクセルiの出力光量COUTiは、入力光量C
INi に修正光量Δi を加算したものとして定義できる。
この5−2式においても、上記の式変形の過程を見れば
明らかなように、透明度(1−αi )の概念は内包され
ており、原理上異ならない。
実施形態につき図面を用いて説明する。
置の全体構成が示されており、図5はそのブロック図で
ある。
音波振動子を有し、そのリニアアレイ型の超音波振動子
を電子走査(リニア走査、セクタ走査)することによっ
て、X−Y面内に走査面10が形成される。この走査面
10を例えば機械的にZ方向へ走査することによって、
図1に示したような三次元エコーデータ取込み空間12
が形成される。ここで、超音波探触子22のZ方向への
機械的な走査は駆動部24によって行われており、ロー
タリーエンコーダ26はその超音波探触子22のZ方向
の座標を検出する。この実施形態では駆動部24によっ
て自動的にZ方向へ超音波探触子22が走査されている
が、もちろん手動により超音波探触子22を走査させて
もよい。この場合においてもZ方向の座標の取込みが必
要となる。なお、Z座標は後述するデジタルスキャンコ
ンバータ(DSC)28に送られる。
子22に対して送信信号を供給すると共に、超音波探触
子22から出力された受信信号を処理するものである。
その送受信部30から出力された受信信号は、まずアン
プ32によって増幅された後、LOGアンプ34におい
て対数増幅され、そしてA/D変換器36にてデジタル
信号に変換された後に、画像処理部38に送られる。画
像処理部38は、受信信号の入力に同期して上述した図
4に示すような処理を行って、最終的に各超音波ビーム
ごとに、対応する画素の輝度値を決定し、DSC28に
出力する。すなわち、このDSC28においては、それ
に内蔵された表示メモリに超音波画像1枚分の画像デー
タが格納される。この場合、各画素の画素データは超音
波ビームのXアドレスと上述したZアドレスとによって
特定されており、それらが表示メモリ上のX、Yアドレ
スに対応づけられ書き込まれる。DSC28から読み出
された画像データはD/A変換器40にてアナログ信号
に戻された後、表示器42に送られ、この表示器42に
て本実施形態に係る三次元超音波画像が表示される。
ィ)調整部44にて決定された係数k1が入力されてお
り、このk1によって各ボクセルの不透明度を調整するこ
とができる。
部38の前段に対数増幅を行うLOGアンプ34が設け
られているが、このようなLOGアンプ34を設けるこ
とによって、より超音波画像を良好なものにできること
が実験により確認されている。
具体的な構成例が示されている。
タei は、ラッチ回路50においていったんラッチされ
た後、そのエコーデータei が発光量演算部52,不透
明度演算部54及び透明度演算部56に入力される。こ
こで、発光量演算部52は、乗算器58で構成され、そ
の乗算器58は、ei ・αi の計算を行って発光量を求
める。不透明度演算部54は、この実施形態においてα
ROMで構成され、すなわちエコーデータei と不透明
度αとの関係を示すテーブルで構成されている。これと
同様に、透明度演算部56は、(1−α)ROMで構成
され、すなわちエコーデータei と透明度(1−α)と
の関係を示すテーブルで構成されている。
部54に入力されると、その出力から不透明度αi が出
力され、一方、そのエコーデータei が透明度演算部5
6に入力されると、その出力から透明度(1−αi )が
出力される。
部56,ラッチ回路50,ラッチ回路60及びラッチ回
路62には、ANDゲート64を介してサンプリングク
ロックが供給されている。このサンプリングクロック
は、図5に示したA/D変換器36に入力されるA/D
変換のためのサンプリングクロックである。すなわち、
図6に示す各回路はA/D変換器36の動作に同期して
動作する。
から出力された透明度(1−αi )に対してラッチ回路
62にてラッチされた一つ前のボクセルの出力光量C
OUTi-1を乗算するものであり、具体的には乗算器68で
構成される。すなわち、この透過光量演算部66から当
該ボクセルの入力光量に透明度を乗算した透過光量が出
力される。
いて、発光量と透過光量とを加算し、当該ボクセルの出
力光量COUTiを出力する。具体的には、この光量加算部
70は加算器72により構成されている。光量加算部7
0から出力された出力光量COUTiは、ゲート機能を有す
るラッチ回路74と上述したラッチ回路62に送られて
いる。すなわち、ラッチ回路62を介して次のボクセル
の演算のために当該ボクセルの出力光量が帰還されてい
る。終了判定部77は、超音波ビームに沿って行われる
上述のボクセル処理の終了判定を行うものであり、具体
的には加算器76,ラッチ回路60及びコンパレータ7
8で構成される。加算器76には、当該ボクセルの不透
明度αi が順次入力されており、一方、その加算器76
の出力はラッチ回路60を介して帰還されて加算器76
の他方の入力に入力されており、これによって加算器7
6の出力として不透明度を順次積算した積算値が出力さ
れることになる。コンパレータ78は、その積算値と所
定の係数Kとを比較し、両者が一致した時点で終了判定
パルスを出力する。具体的には、Kとして1が設定され
ており、終了判定部77は不透明度を積算しその積算値
が1に到達した時点で終了判定パルスを出力する。その
終了判定パルスは、ANDゲート64に入力されサンプ
リングクロックの通過を停止させ、またラッチ回路74
に送られて光量加算部70から出力された出力光量C
OUTiのラッチすなわち通過を許容させる。その出力光量
COUTiは当該超音波ビームに対応する画素の輝度値P
(x,y)となる。ここで、出力光量COUTiを色相に対
応させることもできる。
ていないが、設定された深さのボクセルまで処理が到達
した場合には、上述同様にエコーデータの処理を停止さ
せる。または、超音波ビームの最後のボクセルまで処理
が進行した場合にはその処理を停止させる。超音波ビー
ム上の開始ボクセルは、図6には示されていないが、開
始ボクセルの最初に存在するボクセルとして設定され、
あるいは自動的に又は人為的に設定された深さのボクセ
ルとして設定される。この場合、開始ボクセルの条件を
設定する開始ボクセル設定器を設けるのが望ましい。
によれば、各超音波ビームごとにその超音波ビームに沿
って開始ボクセルからボクセル処理が開始され、そのボ
クセル処理が各ボクセルについて行われ、その際、光量
加算部70の出力が次のボクセル処理で利用される。そ
して、最終ボクセルの出力光量が当該超音波ビームに対
応する画素の輝度値Pとして利用される。1つの超音波
ビームに対する処理が終了した後、次の超音波ビームに
ついて処理がなされ、最終的に超音波画像1枚分につい
てのエコーデータの処理が完了する。その後、DSCか
らその超音波画像が読み出されて表示器42に表示され
る。
述した第1式に示される関数に従って決定され、その第
1式におけるk1は不透明度調整部44によって調整可能
である。すなわち、この不透明度調整部44によって、
表現できる深さを調整でき、また透明感や立体感などを
調整できる。
は、レントゲン像のような透明感を有し、かつ臓器に立
体感を持たせた表現を行うことができ、診断上有用な超
音波画像を構成できる。また、従来のボリューム・リン
ダリング処理ではリアルタイムで画像処理を行うことが
きわめて困難であったが、図5及び図6に示した実施形
態によれば、超音波ビームに沿って各ボクセルの処理が
行われ、しかもテーブルが有効に利用されているため、
リアルタイムでその処理を行うことができ、実用性に優
れる。
ームに沿って各ボクセルの処理を行ったが、場合によっ
ては三次元エコーデータ取込み空間に対してその外部に
視点を設定しその視点から放射状に複数の透視線を設定
し、各透視線に沿ってボクセルの処理を行うこともでき
る。また、上記実施態様では、不透明度演算部54と透
明度演算部56をそれぞれ別個のROMで構成したが、
2進法回路では1つのROMを共用しつつ、ROMの出
力を反転して1の補数を求めることによって透明度を求
める構成を採用することもできる。さらに、上記実施形
態においては、係数Kとして1が設定されていたが、こ
の値を可変させる手段を設けることによって、リアルタ
イムで超音波画像を観察しながら超音波画像の画質を調
整してもよい。
処理部38の他の実施形態が示されている。ROM82
はエコーデータei から不透明度αi を出力している。
一方、減算器84は(ei −CINi )を演算している。
そして、乗算器86では上記5−3式に相当するΔi の
演算が行われ、上記5−2式に基づいて加算器88では
CINi +Δi の演算が行われている。遅延器90では次
のボクセルについての演算のために、加算器88の出力
COUTiが1ボクセル分遅延されている。その遅延器90
の出力は次のボクセルについての入力光量CINi とな
り、減算器84及び加算器88に入力されている。
構成することも可能であるが、ソフトウエアで構成する
こともできる。
生体内部を透明感をもって及び立体感をもって表現する
ことができ、診断上有用な超音波画像を提供できる。ま
た、本発明によれば、そのような画像をリアルタイムで
構成できるという利点がある。また、本発明によれば、
観察可能な範囲を自在に設定することが可能となる。
を示す図である。
示す図である。
る。
を示すブロック図である。
す図である。
図である。
処理部、44 不透明度調整部、52 発光量演算部、
54 不透明度演算部、56 透明度演算部、66 透
過光量演算部、70 光量加算部、77 終了判定部。
Claims (26)
- 【請求項1】 超音波ビームに沿って開始ボクセルから
終了ボクセルまでボクセル処理を順次行うことにより、
その超音波ビームに対応する画素の画素値を求める超音
波画像処理装置であって、 エコー値ei に基づきボクセルi の不透明度αi を演算
する不透明度演算手段と、 エコー値ei に基づきボクセルi の透明度βi を演算す
る透明度演算手段と、 エコー値ei に不透明度αi を乗算し、ボクセルi の発
光量を演算する発光量演算手段と、 1つ前のボクセルi-1 の出力光量にボクセルi の透明度
βi を乗算し、ボクセルi の透過光量を演算する透過光
量演算手段と、 前記発光量と前記透過光量とを加算し、ボクセルi の出
力光量を求める光量加算手段と、 を含み、 終了ボクセルの出力光量を画素値に対応させて超音波画
像を形成することを特徴とする超音波画像処理装置。 - 【請求項2】 請求項1記載の装置において、 前記不透明度αi は、0以上1以下の範囲内で設定さ
れ、 前記透明度βi は(1−αi )で設定されることを特徴
とする超音波画像処理装置。 - 【請求項3】 請求項2記載の装置において、 前記不透明度αi は、エコー値ei をパラメータとする
関数により定義されることを特徴とする超音波画像処理
装置。 - 【請求項4】 請求項3記載の装置において、 前記関数はエコー値ei の増大に対して不透明度αi が
増大する非線形関数であることを特徴とする超音波画像
処理装置。 - 【請求項5】 請求項3記載の装置において、 前記関数に含まれる係数を可変する係数調整手段を含む
ことを特徴とする超音波画像処理装置。 - 【請求項6】 請求項1記載の装置において、 前記終了ボクセルを判定する終了判定手段を含むことを
特徴とする超音波画像処理装置。 - 【請求項7】 請求項6記載の装置において、 前記終了判定手段は、設定された深さまでボクセル処理
が到達した場合に前記終了ボクセルを判定することを特
徴とする超音波画像処理装置。 - 【請求項8】 請求項6記載の装置において、 前記判定終了手段は、開始ボクセルから各ボクセルごと
にその不透明度αi を加算し、その加算値が所定の終了
判定値となった場合に前記終了ボクセルを判定すること
を特徴とする超音波画像処理装置。 - 【請求項9】 請求項8記載の装置において、 前記不透明度αi は、0以上1以下の範囲内で設定さ
れ、 前記透明度βi は(1−αi )で設定され、 前記終了判定値は1であることを特徴とする超音波画像
処理装置。 - 【請求項10】 請求項1記載の装置において、 前記開始ボクセルを設定する手段を有することを特徴と
する超音波画像処理装置。 - 【請求項11】 請求項1記載の装置において、 前記不透明度演算手段は、各エコー値ei に対する不透
明度αi を対応付けたテーブルで構成されたことを特徴
とする超音波画像処理装置。 - 【請求項12】 請求項1記載の装置において、 前記透明度演算手段は、各エコー値ei に対する透明度
βi を対応付けたテーブルで構成されたことを特徴とす
る超音波画像処理装置。 - 【請求項13】 請求項1記載の装置において、 超音波の送受波により得られた受信信号を対数増幅する
対数増幅器を有することを特徴とする超音波画像処理装
置。 - 【請求項14】 三次元領域に対して測定波を送受波
し、三次元領域を構成する各ボクセルごとにボクセル値
ei を取り込む送受波手段と、 前記三次元領域を通過する仮想的な透視線を設定し、そ
の透視線に沿って開始ボクセルから終了ボクセルまでボ
クセル処理を順次行うことにより、その透視線に対応す
る画素の画素値を求める画像処理手段と、 を含む装置であって、 前記画像処理手段は、 ボクセル値ei に基づきボクセルi の不透明度αi を演
算する不透明度演算手段と、 ボクセル値ei に基づきボクセルi の透明度βi を演算
する透明度演算手段と、 ボクセル値ei と不透明度αi とに基づき、ボクセルi
の発光量を演算する発光量演算手段と、 1つ前のボクセルi-1 の出力光量とボクセルi の透明度
βi とに基づき、ボクセルi の透過光量を演算する透過
光量演算手段と、 前記発光量と前記透過光量とに基づき、ボクセルi の出
力光量を求める出力光量演算手段と、 を含み、 終了ボクセルの出力光量を画素値に対応させて画像を形
成することを特徴とする装置。 - 【請求項15】 請求項14記載の装置において、 前記透視線を測定波ビームに一致させたことを特徴とす
る装置。 - 【請求項16】 三次元領域を通過する仮想的な透視線
上で、開始ボクセルから終了ボクセルまでボクセル処理
を順次行うことにより、透視線に対応する画素の画素値
を求める画像処理方法において、 ボクセル値ei によりボクセルi の不透明度αi (0≦
αi ≦1)を演算する工程と、 ボクセルi の透明度(1−αi )を演算する工程と、 ボクセル値ei に不透明度αi を乗算し、ボクセルi の
発光量を演算する工程と、 1つ前のボクセルi-1 の出力光量にボクセルi の透明度
(1−αi )を乗算し、ボクセルi の透過光量を演算す
る工程と、 前記発光量と前記透過光量とを加算し、ボクセルi の出
力光量を求める工程と、を含み、 終了ボクセルの出力光量を輝度値とすることを特徴とす
る画像処理方法。 - 【請求項17】 請求項16記載の方法において、 前記不透明度αi は、k1及びk2を係数として、αi =k1
・ei k2により演算されることを特徴とする画像処理方
法。 - 【請求項18】 三次元領域を通過する仮想的な透視線
上で、開始ボクセルから終了ボクセルまでボクセル処理
を順次行うことにより、透視線に対応する画素の画素値
を求める超音波画像処理装置であって、 ボクセル値ei に基づきボクセルi の不透明度αi を演
算する不透明度演算手段と、 前記ボクセル値ei 、前記不透明度αi 、及び、1つ前
のボクセルi-1 の出力光量に相当する入力光量CINi に
基づいて、ボクセルi の出力光量COUTiを演算する出力
光量演算手段と、 を含み、 終了ボクセルの出力光量を画素値に対応させて超音波画
像を形成することを特徴とする超音波画像処理装置。 - 【請求項19】 請求項18記載の装置において、 前記透視線は超音波ビームに相当することを特徴とする
超音波画像処理装置。 - 【請求項20】 請求項18記載の装置において、 前記出力光量演算手段は、CINi +αi ・(ei −C
INi )の演算により、出力光量COUTiを演算することを
特徴とする超音波画像処理装置。 - 【請求項21】 超音波ビームを三次元領域で走査する
ことにより得られたエコー値に基づき前記三次元領域を
画像化した超音波画像を形成する画像処理手段を含む超
音波画像処理装置において、 前記画像処理手段は、前記超音波ビームに沿った時系列
順で入力される各エコー値に対して逐次処理を行って、
当該超音波ビーム上の複数のエコー値が反映された画素
値を演算し、 各超音波ビームごとに前記画素値が演算され、それらの
画素値の集合として前記三次元領域を画像化した超音波
画像が形成されることを特徴とする超音波画像処理装
置。 - 【請求項22】 請求項21記載の装置において、 前記逐次処理では、前のエコー値の演算結果が次のエコ
ー値の演算で利用され、 終了点のエコー値に対する演算結果として画素値が決定
されることを特徴とする超音波画像処理装置。 - 【請求項23】 請求項22記載の装置において、 超音波ビームに沿って各エコー値を参照することによ
り、各超音波ビームごとに前記終了点を判定する終了判
定手段を含むことを特徴とする超音波画像処理装置。 - 【請求項24】 請求項23記載の装置において、 前記終了判定手段が有する終了判定条件を設定するため
の終了条件設定手段を有することを特徴とする超音波画
像処理装置。 - 【請求項25】 請求項21記載の装置において、 超音波ビームの走査により取り込まれた各エコー値が三
次元データメモリを介することなく前記画像処理手段に
順次入力され、 前記画像処理手段は各エコー値の入力タイミングに同期
して前記逐次処理を行うことを特徴とする超音波画像処
理装置。 - 【請求項26】 請求項21記載の装置において、 前記超音波画像は、前記三次元領域を平面上に投影した
画像であることを特徴とする超音波画像処理装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8185781A JP2883584B2 (ja) | 1995-07-17 | 1996-07-16 | 超音波画像処理装置及び方法 |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP7-180107 | 1995-07-17 | ||
JP18010795 | 1995-07-17 | ||
JP12534196 | 1996-05-21 | ||
JP8-125341 | 1996-05-21 | ||
JP8185781A JP2883584B2 (ja) | 1995-07-17 | 1996-07-16 | 超音波画像処理装置及び方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH1033538A JPH1033538A (ja) | 1998-02-10 |
JP2883584B2 true JP2883584B2 (ja) | 1999-04-19 |
Family
ID=27315101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP8185781A Expired - Fee Related JP2883584B2 (ja) | 1995-07-17 | 1996-07-16 | 超音波画像処理装置及び方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2883584B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006271594A (ja) * | 2005-03-29 | 2006-10-12 | Aloka Co Ltd | 超音波診断装置 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6102861A (en) * | 1999-04-23 | 2000-08-15 | General Electric Company | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering |
JP3944059B2 (ja) * | 2002-11-14 | 2007-07-11 | アロカ株式会社 | 超音波診断装置 |
JP3872424B2 (ja) | 2002-12-20 | 2007-01-24 | アロカ株式会社 | 超音波診断装置 |
JP2004267506A (ja) * | 2003-03-10 | 2004-09-30 | Aloka Co Ltd | 超音波診断装置 |
JP4703144B2 (ja) * | 2004-08-23 | 2011-06-15 | 株式会社東芝 | 超音波診断装置、及び画像処理方法 |
US7810987B2 (en) * | 2005-07-27 | 2010-10-12 | Cargill, Incorporated | Automated solution maker apparatus |
JP4827670B2 (ja) * | 2006-09-20 | 2011-11-30 | 株式会社日立製作所 | 超音波検査装置 |
JP4909137B2 (ja) * | 2007-03-15 | 2012-04-04 | 日立アロカメディカル株式会社 | ボリュームデータ処理装置及び方法 |
JP2009011468A (ja) | 2007-07-03 | 2009-01-22 | Aloka Co Ltd | 超音波診断装置 |
KR101117035B1 (ko) * | 2009-03-24 | 2012-03-15 | 삼성메디슨 주식회사 | 볼륨 데이터에 표면 렌더링을 수행하는 초음파 시스템 및 방법 |
JP5632680B2 (ja) | 2010-08-25 | 2014-11-26 | 日立アロカメディカル株式会社 | 超音波画像処理装置 |
US9655592B2 (en) * | 2014-11-21 | 2017-05-23 | General Electric Corporation | Method and apparatus for rendering an ultrasound image |
-
1996
- 1996-07-16 JP JP8185781A patent/JP2883584B2/ja not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006271594A (ja) * | 2005-03-29 | 2006-10-12 | Aloka Co Ltd | 超音波診断装置 |
JP4653536B2 (ja) * | 2005-03-29 | 2011-03-16 | アロカ株式会社 | 超音波診断装置 |
Also Published As
Publication number | Publication date |
---|---|
JPH1033538A (ja) | 1998-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0754957B1 (en) | Ultrasound image processing apparatus and method for producing an ultrasound image | |
US5766129A (en) | Ultrasound diagnostic apparatus and method of forming an ultrasound image by the apparatus | |
US5911691A (en) | Ultrasound image processing apparatus and method of forming and displaying ultrasound images by the apparatus | |
CN103648400B (zh) | 超声波诊断装置以及方法 | |
US6544178B1 (en) | Methods and systems for volume rendering using ultrasound data | |
JP2883584B2 (ja) | 超音波画像処理装置及び方法 | |
CN1196446C (zh) | 超声图象的相干滤波的方法和装置 | |
US4176658A (en) | Ultrasonic echogram display | |
EP0952462A2 (en) | Method and apparatus for improving visualization of biopsy needle in ultrasound imaging | |
JPH10277030A (ja) | 超音波診断装置 | |
JPH0227631B2 (ja) | ||
JP3015728B2 (ja) | 超音波診断装置 | |
JP3431156B2 (ja) | 超音波走査のための方法及び装置 | |
JP5996268B2 (ja) | 超音波診断装置、画像処理装置、及びプログラム | |
JP3015727B2 (ja) | 超音波診断装置 | |
CN104546005B (zh) | 一种超声非线性成像方法及系统 | |
JP3050374B2 (ja) | 超音波診断装置 | |
JP3522488B2 (ja) | 超音波診断装置 | |
JP3394648B2 (ja) | 超音波診断装置 | |
JPH0751134B2 (ja) | 超音波診断画像処理装置 | |
JP7286025B2 (ja) | 胎盤を評価するためのシステム及び方法 | |
JP3505281B2 (ja) | 超音波3次元画像形成装置 | |
JP3776523B2 (ja) | 超音波診断装置 | |
JP3500014B2 (ja) | 超音波診断装置 | |
EP0809119B1 (en) | Ultrasound image processing apparatus and method of forming and displaying ultra sound images by the apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100205 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100205 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120205 Year of fee payment: 13 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120205 Year of fee payment: 13 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140205 Year of fee payment: 15 |
|
LAPS | Cancellation because of no payment of annual fees |