[go: up one dir, main page]

JP4397131B2 - 3D image display device - Google Patents

3D image display device Download PDF

Info

Publication number
JP4397131B2
JP4397131B2 JP2001092944A JP2001092944A JP4397131B2 JP 4397131 B2 JP4397131 B2 JP 4397131B2 JP 2001092944 A JP2001092944 A JP 2001092944A JP 2001092944 A JP2001092944 A JP 2001092944A JP 4397131 B2 JP4397131 B2 JP 4397131B2
Authority
JP
Japan
Prior art keywords
light emission
emission rate
image
volume
projection plane
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
Application number
JP2001092944A
Other languages
Japanese (ja)
Other versions
JP2001351120A5 (en
JP2001351120A (en
Inventor
良洋 後藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2001092944A priority Critical patent/JP4397131B2/en
Priority to US10/471,548 priority patent/US7019745B2/en
Priority to PCT/JP2001/008716 priority patent/WO2002078545A1/en
Publication of JP2001351120A publication Critical patent/JP2001351120A/en
Publication of JP2001351120A5 publication Critical patent/JP2001351120A5/ja
Application granted granted Critical
Publication of JP4397131B2 publication Critical patent/JP4397131B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Image Generation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Processing Or Creating Images (AREA)

Description

【0001】
【発明が属する技術分野】
本発明は三次元画像表示装置に係り、特に医用画像装置、例えば、X線CT装置、MRI装置等から得られる複数の断層像を使用し、ボリュームレンダリング法により陰影付けして三次元画像を構成し表示する三次元画像表示装置に関する。
【0002】
【従来の技術】
従来技術において、X線CT装置やMRI装置を使用すれば、被検体の断層像が得られることはよく知られている。また、この断層像を複数枚積み上げてなるボリューム画像を、投影面上に陰影付けして3次元的な画像、即ち、三次元画像(3D画像)を構成する方法も、ボリュームレンダリング法、サーフエイス法、デプス法などが知られている。
【0003】
尚、ボリュームレンダリング法を扱った文献として『Marc Levoy; “Display of Surfaces from Volume Data”, IEEE Computer Graphics & Applications, pp29-37, 1988』がある。
【0004】
【発明が解決しようとする課題】
しかしながら、従来のボリュームレンダリング法によって構成された3D画像は、特定部位を他の部位から識別しやすいように強調表示することが事実上できないという問題がある。例えば、健康診断時に病巣部分を発見しやすくするために特定の病巣部分を強調表示したり、手術で血管にステントグラフトを挿入したあとの様子を患者さんに説明するためにステントグラフトを強調表示することができなかった。
【0005】
また、上記特定の病巣部分やステントグラフトに色付けして、他の部位との差を明確にする表示も行われていなかった。
【0006】
本発明はこのような事情に鑑みてなされたもので、ボリュームレンダリング法により特定部位を強調表示した三次元画像を構成し表示することができる三次元画像表示装置を提供することを目的とする。
【0007】
【課題を解決するための手段】
前記目的を達成するために、本願請求項1に係る発明は、複数枚の断層像が積み上げられてなるボリューム画像を所定の投影面上に陰影付けして投影し、三次元画像を表示する三次元画像表示装置であって、前記陰影付けのアルゴリズムとして前記ボリューム画像の外に設定した仮想光源から光線を照射し、該光線の前記ボリューム画像の各画素での反射率及び濃度勾配に基づいて前記投影面上の画素値を算出するボリュームレンダリング法を用いた三次元画像表示装置において、前記ボリューム画像中の表示対象物のうちから強調表示する特定の対象物を示す画素に対して該画素に入射する光線を増幅させる発光率を割り当てる発光率割当手段と、前記発光率割当手段によって割り当てられた発光率をボリュームレンダリング法における陰影付け計算のパラメータとして追加し、投影面上の画素値を算出する演算手段と、を備え、前記発光率が割り当てられた特定の対象物を強調表示するようにしたことを特徴としている。
【0008】
また、前記発光率割当手段は、本願請求項2に示すように被検体の各部位の病巣ごとに、その病巣を示す画素の濃度値に対応して発光率が割り当てられた発光率パターンを有し、強調表示させたい部位の病巣を選択指定することにより、該部位の病巣に対応する発光率パターンを選択することを特徴としている。
【0009】
即ち、特定の病巣部分を強調表示させたい場合には、その病巣に対応する発光率パターンを選択するだけでよく、病巣がある場合には病巣部分が発光する(明るくなる)ように表示される。
【0010】
また、本願請求項3に係る発明は、複数枚の断層像が積み上げられてなるボリューム画像を所定の投影面上に陰影複数枚の断層像が積み上げられてなるボリューム画像を所定の投影面上に陰影付けして投影し、三次元画像を表示する三次元画像表示装置であって、前記陰影付けのアルゴリズムとして前記ボリューム画像の外に設定した仮想光源から光線を照射し、該光線の前記ボリューム画像の各画素での反射率及び濃度勾配に基づいて前記投影面上の画素値を算出するボリュームレンダリング法を用いた三次元画像表示装置において、前記ボリューム画像中の表示対象物のうちから強調表示する特定の対象物を示す画素に対して色発光率を割り当てる色発光率割当手段と、前記色発光率割当手段によって割り当てられた色発光率をボリュームレンダリング法における陰影付け計算のパラメータとして追加し、投影面上の画素値を算出する演算手段と、を備え、前記色発光率が割り当てられた特定の対象物を強調表示するようにしたことを特徴としている。
【0011】
即ち、特定の病巣部分を強調表示させたい場合には、その病巣に対応する色発光率パターンを選択するだけでよく、病巣がある場合には病巣部分が色付けされて発光するように表示される。
【0012】
【発明の実施の形態】
以下添付図面に従って本発明に係る三次元画像表示装置の好ましい実施の形態について説明する。
【0013】
まず、本発明に適用されるボリュームレンダリング法による三次元画像(3D画像)の構成方法について説明する。
【0014】
一般に、ボリュームレンダリング法では、図1に示すように複数枚の断層像が積み上げられてなるボリューム画像の各画素を、立方体の形状をもったボクセルと呼ぶ。従って、ボリューム画像は、このボクセルを三次元的に積み上げた画像となる。ボクセルのデータは、X線CT装置であれば、そのボクセルに対応する人体の箇所のX線吸収値に基づいて算出されるCT値、MRI装置であればプロトン密度等の測定値である。
【0015】
また、ボリュームレンダリング法による陰影付けのパラメータとしては、CT値に対して設定した光の不透明度や、CT値のつくる等高面の勾配(以下、濃度勾配という)等が使用される。
【0016】
ボリュームレンダリング法は、上記ボリューム画像のボクセルのデータ(以下、CT値という)に対して適宜設定した閾値条件に合ったCT値を抽出することにより、骨や皮膚等の所望の被検体を区別して抽出するとともに、そのCT値を投影面上の陰影付けのパラメータとして使用することもできる。
【0017】
即ち、図1に示すように仮想光源から出射される仮想光線は、ボリューム画像の各ボクセルを通過するが、その時各ボクセルのCT値で定まる不透明度にしたがって減衰する。従って、各ボクセルにおける反射光の大きさを求めることができ、ボリュームレンダリング法による陰影付けアルゴリズムでは、反射光が大きい程、そのボクセルを明るくする。また、ボクセルのCT値が大きく変化するところに仮想的な表面が存在すると考え、その仮想的な表面が正面を向いている程、即ち、濃度勾配が大きい程、そのボクセルを明るくする。
【0018】
そして、仮想光線が0に減衰するまで、上記ボクセルの明るさの情報を視点から奥行き方向の各ボクセルにわたって累算した情報を、投影面上の各画素の明るさを示す値(以下、画素値という)とし、投影面上の全画素の画素値を求めるようにしている。
【0019】
本発明では、特定部位を強調表示できるように投影面上の画素値を、次式により求める。
【0020】
【数1】

Figure 0004397131
但し、V:CT値、濃度値
ξ:反射率
η:発光率(負の吸収率)
Δ:濃度勾配
δ(V):約0.1〜6.0の実数で、Vに依存する。
入射光線の強度=1と仮定する。
【0021】
上記式1に示すように、投影面上の画素値は、仮想光線の方向に並ぶ各ボクセル(j=0〜k)の明るさの情報を累算した情報である。ここで、kは、仮想光線がボクセルを通過するごとに反射、減衰して0(即ち、式1の{ }内が0)になるまでの累算されるボクセルの数である。
【0022】
さて、本発明では、式1の[ ]内の2項目に、負の値をもつ発光率η(Vi)が追加されている。即ち、一般に仮想光線は各ボクセルを通過するごとに減衰し、増幅されることはないが、本発明では、強調表示させたい部位のボクセルに対し、入射光線の強度に対する比で表してなる発光率を割り当て、結果的に強調表示させたい部位のボクセルでは入射光線の強度を増強するように作用させる。
【0023】
図2はボクセルのCT値に応じて反射率ξが割り当てられた反射率パターンA、及び発光率ηが割り当てられた発光率パターンBの一例を示している。尚、発光率ηは負の値をもち、前述した式1の{ }の値(即ち、入射光線の強度)を増加させるように作用する。また、図2上で、Cは吸収率のパターンを示しており、上記のようにして定義した発光率は、仮想光線の負の吸収率で近似することができる。
【0024】
図3は上記式1に基づいて作成した3D画像等を表示するCRTディスプレイの表示画面の一例を示している。
【0025】
同図に示すように、反射率割当部1では、○印をマウスで移動させることにより、CT値に対して所望の反射率を割り当てることができ、同様に発光率割当部2では、□印をマウスで移動させることにより、CT値に対して所望の発光率を割り当てることができる。このようにして割り当てた反射率及び発光率に基づいて式1より投影面上の画素値を計算し、3D画像3を構成し表示する。
【0026】
また、発光率割当部2では、肝臓の病巣部位、軟骨の病巣部位、骨の病巣部位等をマウスで選択できるようになっており、これらの病巣部位から所望の病巣部位を選択すると、その病巣部位を強調表示するための発光率パターンが自動的に設定できるようにもなっている。
【0027】
尚、立方体内の3D画像3の視点位置を変えるには、前記立方体の角の白丸4をマウスでドラッグし、前記立方体を立方体の中心点を基準に適宜回転させる。これにより、3D画像3の視線方向が決定し、新たに3D画像が構成表示される。上記のようにして構成された3D画像は必要に応じて磁気ディクスに格納される。また、図3上で、5は断層像(CT画像)を示している。
【0028】
図4(a)は脳に対応するCT値に対して発光率を設定して作成された頭部の3D画像10を示し、図4(b)は発光率を設定せずに(全てのCT値の発光率が0)作成された頭部の3D画像12を示している。図4(a)に示す脳(特に眼球部位11)は明るく表示され、図4(b)に示す脳(特に眼球部位13)は暗く表示されている。
【0029】
また、図5(a)はステントグラフト15に対応するCT値に対して発光率を設定して作成された血管の3D画像14を示し、図4(b)は発光率を設定せずに作成された血管の3D画像16とを示している。図5(a)に示すステントグラフト15は明るく表示され、図5(b)に示すステントグラフト17は暗く表示されている。
【0030】
図6は本発明に係る3D画像表示の処理手順を示すフローチャートである。
〔ステップ50〕
反射率、発光率、及び初期値の設定を行う。
〔ステップ52〕
視線方向(仮想光線の方向)の初期値を設定する。
〔ステップ54〕
式1に基づいて投影面上の画素値を計算し、3D画像を構成し表示する。CRTディスプレイに3D画像を表示するには、表示諧調に合わせるためのテーブルを参照して投影面メモリ上の値を変換したあとで、表示メモリに書き込むことにより表示される。
〔ステップ56、58、60〕
視線方向、反射率、及び発光率のいずれかが変更されたか否かを判別する。いずれかが変更されると、ステップ54に飛び、再度3D画像を構成し表示する。
〔ステップ62〕
現在表示されている3D画像を格納するか否かを調べる。格納指示があれば、ステップ64に飛び、なければステップ66に移行する。
〔ステップ64〕
現在表示されている3D画像を磁気ディクスに格納する。
〔ステップ66〕
3D画像の表示等の終了指示があったか否かを判定し、終了指示がなければステップ56に跳び、終了指示があると、処理動作を終了させる。
【0031】
また、図7に示すように、発光させる領域を特定の関心領域としてもよい。その関心領域は、図7(a)に示すように、投影面と発光対象部位との距離で決まるものである。また、表示画像からの関心領域の設定の仕方は、図7(b)に示すように、操作者がマウスを使い発光させたい領域を囲み、奥行き方向の設定は、操作者が所望の発光対象部位と投影面との距離を適宜設定すればよい。
【0032】
また、これまでの説明は白黒表示について主に説明したが、カラー画像として適用できることもいうまでもない。この場合は、白黒画像での発光率をカラー画像の色発光率に代えることで実現できる。
【0033】
具体的には、図8に示すように、赤(R)、緑(G)、青(B)の色成分につき、別々に色発光できるようにしてある。横軸はCT値であるが、例えば脳を赤く発光させたい場合、成分G、Bについては脳を示すCT値範囲の光の反射率(不透明度:オパシティともいう)をマウスの操作により低く抑制し、成分Rの発光率を大きくする。これにより、脳を赤く発光できることとなる。なお、このカラー発光は一例であり、RGBによって決められる発色系によらないものも含まれることはいうまでもないことである。
【0034】
図9は本発明に係る三次元画像表示装置のハードウェア構成を示すブロック図である。この三次元画像表示装置は、例えばX線CT装置等で被検体の対象部位について収集した複数の断層像により、式1に基づいて3D画像等を構成して表示するもので、各構成要素の動作を制御する中央処理装置(CPU)20と、装置の制御プログラムが格納された主メモリ22と、複数の断層像データ及びプログラム等が格納された磁気ディスク24と、表示用の画像データを一時記憶する表示メモリ26と、この表示メモリ26からの画像データに基づいて画像を表示する表示装置としてのCRTモニタ28と、位置入力装置としてのマウス30と、マウス30の状態を検出してCRTモニタ18上のマウスポインタの位置やマウス30の状態等の信号をCPU20に出力するマウスコントローラ32と、各種の操作指令等を入力するキーボード34と、上記各構成要素を接続する共通バス26とから構成される。
【0035】
尚、この実施の形態では、式1に基づいて投影面上の画素値を計算し、3D画像を構成するようにしたが、これに限らず、下記の式2、式3等によって投影面上の画素値を計算してもよい。
【0036】
【数2】
Figure 0004397131
但し、Cは定数
【0037】
【数3】
Figure 0004397131
を加えて反射率を増大させており、強調度合いを更に強めるようにしている。また、上記式3は、入射強度を乗算して減衰させており、式1と比べて入射強度の減衰が遅いため、比較的内部まで透視可能な3D画像を構成することができる。
【0038】
また、本発明には、発光率および発光色をそれぞれ割当てることで説明したが、それらを組み合わせ、その組み合わせの結果として特定部位やステントを強調表示することも含まれる。
【0039】
【発明の効果】
以上説明したように本発明によれば、ボリュームレンダリング法により特定部位を強調表示した三次元画像を構成し表示することができ、健康診断時に病巣部分を発見しやすい等の効果がある。
【図面の簡単な説明】
【図1】ボリュームレンダリング法による陰影付けを説明するために用いた図。
【図2】本発明に係る発光率等を説明するために用いた図。
【図3】本発明に基づいて作成した3D画像等を表示するCRTディスプレイの表示画面の一例を示す図。
【図4】本発明によって表示される3D画像と従来の3D画像との一例を示す図。
【図5】本発明によって表示される3D画像と従来の3D画像との他の例を示す図。
【図6】本発明に係る3D画像表示の処理手順を示すフローチャート。
【図7】本発明の発光率、発光させる色等を設定する例を説明する図。
【図8】本発明の発光させる色を設定する例、およびその表示例を示す図。
【図9】本発明に係る三次元画像表示装置のハードウェア構成を示すブロック図。
【符号の説明】
1…反射率割当部、2…発光率割当部、3、10、12、14、、16…3D画像、15、17…ステントグラフト、20…CPU、22…主メモリ、24…磁気ディスク、26…表示メモリ、28…CRTモニタ、30…マウス、32…コントローラ、34…キーボード、36…共通バス[0001]
[Technical field to which the invention belongs]
The present invention relates to a three-dimensional image display apparatus, and in particular, uses a plurality of tomographic images obtained from a medical image apparatus, for example, an X-ray CT apparatus, an MRI apparatus, etc., and forms a three-dimensional image by shading by a volume rendering method. The present invention relates to a three-dimensional image display device for displaying and displaying.
[0002]
[Prior art]
In the prior art, it is well known that a tomographic image of a subject can be obtained if an X-ray CT apparatus or an MRI apparatus is used. In addition, a volume rendering method and a surf face method are also used as a method of constructing a three-dimensional image, that is, a three-dimensional image (3D image) by shading a volume image formed by stacking a plurality of tomographic images on a projection plane. The depth method is known.
[0003]
As a document dealing with the volume rendering method, there is “Marc Levoy;“ Display of Surfaces from Volume Data ”, IEEE Computer Graphics & Applications, pp29-37, 1988”.
[0004]
[Problems to be solved by the invention]
However, a 3D image constructed by the conventional volume rendering method has a problem that it is practically impossible to highlight a specific part so that it can be easily distinguished from other parts. For example, a specific lesion may be highlighted to make it easier to find the lesion during a health checkup, or the stent graft may be highlighted to explain to the patient what the stent graft has been inserted into the blood vessel during surgery. could not.
[0005]
In addition, the above-mentioned specific lesion part or stent graft is colored to display a difference from other parts.
[0006]
The present invention has been made in view of such circumstances, and an object of the present invention is to provide a three-dimensional image display device capable of constructing and displaying a three-dimensional image in which a specific part is highlighted by a volume rendering method.
[0007]
[Means for Solving the Problems]
In order to achieve the above-mentioned object, the invention according to claim 1 of the present application is a tertiary image displaying a three-dimensional image by projecting a volume image formed by stacking a plurality of tomographic images by shading on a predetermined projection plane. In the original image display device, a light source is irradiated from a virtual light source set outside the volume image as the shading algorithm, and the light beam is based on a reflectance and a density gradient at each pixel of the volume image. In a three-dimensional image display device using a volume rendering method for calculating a pixel value on a projection plane, the pixel is incident on a pixel indicating a specific object to be highlighted from among the display objects in the volume image. A light emission rate assigning means for assigning a light emission rate for amplifying the light rays to be amplified, and the light emission rate assigned by the light emission rate assigning means in the volume rendering method. With added as parameter calculation, and calculation means for calculating a pixel value on the projection plane, comprising a, it is characterized in that so as to highlight the particular object in which the light-emitting rate is allocated.
[0008]
Further, as shown in claim 2 of the present application, the light emission rate assigning means has a light emission rate pattern in which a light emission rate is assigned to each lesion of each part of the subject corresponding to a density value of a pixel indicating the lesion. Then, by selecting and specifying a lesion at a site to be highlighted, a light emission rate pattern corresponding to the lesion at the site is selected.
[0009]
That is, when a specific lesion portion is to be highlighted, it is only necessary to select a light emission rate pattern corresponding to the lesion, and when there is a lesion, the lesion portion is displayed to emit light (become bright). .
[0010]
In the invention according to claim 3 of the present application, a volume image formed by stacking a plurality of tomographic images on a predetermined projection plane and a volume image formed by stacking a plurality of shaded tomographic images on a predetermined projection plane. A three-dimensional image display device that displays a three-dimensional image by shading and projecting, and irradiating light from a virtual light source set outside the volume image as the shading algorithm, and the volume image of the light In the three-dimensional image display device using the volume rendering method that calculates the pixel value on the projection plane based on the reflectance and the density gradient at each pixel, highlighting is performed from among the display objects in the volume image. A color light emission rate assigning unit for assigning a color light emission rate to a pixel indicating a specific object, and a color light emission rate assigned by the color light emission rate assigning unit. An arithmetic means for calculating a pixel value on a projection plane added as a parameter for shading calculation in the rendering method, and highlighting a specific object to which the color light emission rate is assigned. It is said.
[0011]
That is, when a specific lesion portion is to be highlighted, it is only necessary to select a color light emission rate pattern corresponding to the lesion, and when there is a lesion, the lesion portion is displayed to be colored and emit light. .
[0012]
DETAILED DESCRIPTION OF THE INVENTION
Preferred embodiments of a three-dimensional image display device according to the present invention will be described below with reference to the accompanying drawings.
[0013]
First, a method for constructing a three-dimensional image (3D image) by the volume rendering method applied to the present invention will be described.
[0014]
In general, in the volume rendering method, each pixel of a volume image formed by stacking a plurality of tomographic images as shown in FIG. 1 is called a voxel having a cubic shape. Therefore, the volume image is an image in which the voxels are stacked three-dimensionally. The voxel data is a CT value calculated based on an X-ray absorption value of a human body corresponding to the voxel in the case of an X-ray CT apparatus, and a measured value such as proton density in the case of an MRI apparatus.
[0015]
Further, as the parameters for shading by the volume rendering method, the opacity of light set for the CT value, the gradient of the contour plane created by the CT value (hereinafter referred to as density gradient), and the like are used.
[0016]
The volume rendering method distinguishes a desired subject such as a bone or skin by extracting a CT value that meets a threshold condition appropriately set for voxel data (hereinafter referred to as a CT value) of the volume image. In addition to extraction, the CT value can also be used as a parameter for shading on the projection plane.
[0017]
That is, as shown in FIG. 1, the virtual light beam emitted from the virtual light source passes through each voxel of the volume image, but at that time, attenuates according to the opacity determined by the CT value of each voxel. Therefore, the magnitude of the reflected light in each voxel can be obtained. In the shading algorithm based on the volume rendering method, the larger the reflected light, the brighter the voxel. Further, it is considered that a virtual surface exists where the CT value of the voxel greatly changes, and the voxel is brightened as the virtual surface is facing the front, that is, the density gradient is larger.
[0018]
A value indicating the brightness of each pixel on the projection plane (hereinafter referred to as pixel value) is obtained by accumulating the information on the brightness of the voxel over each voxel in the depth direction from the viewpoint until the virtual ray is attenuated to 0. The pixel values of all the pixels on the projection plane are obtained.
[0019]
In the present invention, the pixel value on the projection plane is obtained by the following equation so that the specific part can be highlighted.
[0020]
[Expression 1]
Figure 0004397131
However, V: CT value, density value ξ: reflectance η: light emission rate (negative absorption rate)
Δ: concentration gradient δ (V): a real number of about 0.1 to 6.0, depending on V.
Assume that the intensity of the incident light = 1.
[0021]
As shown in the above formula 1, the pixel value on the projection plane is information obtained by accumulating the brightness information of the voxels (j = 0 to k) arranged in the direction of the virtual ray. Here, k is the number of voxels that are accumulated until a virtual ray passes through the voxel and is reflected and attenuated to 0 (that is, 0 in {} in Equation 1).
[0022]
Now, in the present invention, the light emission rate η (Vi) having a negative value is added to the two items in [] in Equation 1. That is, in general, a virtual ray attenuates and does not amplify as it passes through each voxel. However, in the present invention, a luminous rate expressed by a ratio to the intensity of incident light with respect to a voxel of a portion to be highlighted. , And as a result, the voxel at the part to be highlighted is operated to increase the intensity of the incident light.
[0023]
FIG. 2 shows an example of the reflectance pattern A to which the reflectance ξ is assigned according to the CT value of the voxel, and the light emission pattern B to which the light emission rate η is assigned. Note that the light emission rate η has a negative value, and acts to increase the value of {} in Equation 1 (that is, the intensity of incident light). In FIG. 2, C indicates an absorptance pattern, and the light emission rate defined as described above can be approximated by a negative absorptance of a virtual ray.
[0024]
FIG. 3 shows an example of a display screen of a CRT display that displays a 3D image or the like created based on Equation 1 above.
[0025]
As shown in the figure, in the reflectance assignment unit 1, a desired reflectance can be assigned to the CT value by moving the circle mark with a mouse. Similarly, in the light emission rate assignment unit 2, the square mark Is moved with a mouse, a desired light emission rate can be assigned to the CT value. A pixel value on the projection plane is calculated from Equation 1 based on the reflectance and the light emission rate assigned in this way, and a 3D image 3 is constructed and displayed.
[0026]
Further, the light emission rate allocating unit 2 can select a lesion site of the liver, a lesion site of the cartilage, a lesion site of the bone, etc. with a mouse. When a desired lesion site is selected from these lesion sites, the lesion site is selected. It is also possible to automatically set a light emission rate pattern for highlighting a region.
[0027]
In order to change the viewpoint position of the 3D image 3 in the cube, the white circle 4 at the corner of the cube is dragged with the mouse, and the cube is appropriately rotated based on the center point of the cube. Thereby, the line-of-sight direction of the 3D image 3 is determined, and a new 3D image is configured and displayed. The 3D image configured as described above is stored in the magnetic disk as necessary. In FIG. 3, reference numeral 5 denotes a tomographic image (CT image).
[0028]
FIG. 4 (a) shows a 3D image 10 of the head created by setting the light emission rate with respect to the CT value corresponding to the brain, and FIG. 4 (b) shows that all the CT values are set without setting the light emission rate. The light emission rate of the value is 0) and the created 3D image 12 of the head is shown. The brain (particularly the eyeball part 11) shown in FIG. 4A is displayed brightly, and the brain (particularly the eyeball part 13) shown in FIG. 4B is displayed darkly.
[0029]
FIG. 5A shows a 3D image 14 of a blood vessel created by setting the light emission rate with respect to the CT value corresponding to the stent graft 15, and FIG. 4B is created without setting the light emission rate. A 3D image 16 of a blood vessel is shown. The stent graft 15 shown in FIG. 5A is displayed brightly, and the stent graft 17 shown in FIG. 5B is displayed darkly.
[0030]
FIG. 6 is a flowchart showing a 3D image display processing procedure according to the present invention.
[Step 50]
A reflectance, a light emission rate, and an initial value are set.
[Step 52]
Sets the initial value of the line-of-sight direction (the direction of virtual rays).
[Step 54]
A pixel value on the projection plane is calculated based on Equation 1, and a 3D image is constructed and displayed. In order to display a 3D image on a CRT display, a value on the projection plane memory is converted with reference to a table for adjusting display gradation, and then written into the display memory.
[Steps 56, 58, 60]
It is determined whether any of the line-of-sight direction, the reflectance, and the light emission rate has been changed. If any of them is changed, the process jumps to step 54, and a 3D image is constructed and displayed again.
[Step 62]
It is checked whether or not to store the currently displayed 3D image. If there is a storage instruction, the process jumps to step 64;
[Step 64]
The currently displayed 3D image is stored in the magnetic disk.
[Step 66]
It is determined whether or not there is an end instruction for displaying a 3D image or the like. If there is no end instruction, the process jumps to step 56, and if there is an end instruction, the processing operation is ended.
[0031]
Further, as shown in FIG. 7, a region to emit light may be a specific region of interest. The region of interest is determined by the distance between the projection plane and the light emission target part, as shown in FIG. Further, as shown in FIG. 7B, the method of setting the region of interest from the display image encloses the region where the operator wants to emit light using the mouse, and the depth direction is set by the operator as desired light emission target. What is necessary is just to set the distance of a site | part and a projection surface suitably.
[0032]
In addition, although the description so far has mainly described the monochrome display, it is needless to say that it can be applied as a color image. In this case, it can be realized by replacing the light emission rate of a black and white image with the color light emission rate of a color image.
[0033]
Specifically, as shown in FIG. 8, the color components of red (R), green (G), and blue (B) can be emitted separately. The horizontal axis is the CT value. For example, when you want the brain to emit red light, the light reflectance (opacity: also called opacity) of the CT value range indicating the brain is suppressed to a low level by operating the mouse for components G and B. And increase the emission rate of component R. As a result, the brain can emit red light. Note that this color light emission is an example, and it is needless to say that light emission not depending on the color development system determined by RGB is included.
[0034]
FIG. 9 is a block diagram showing a hardware configuration of the 3D image display apparatus according to the present invention. This three-dimensional image display apparatus is configured to display a 3D image or the like based on Equation 1 by using a plurality of tomographic images collected on a target region of a subject with an X-ray CT apparatus, for example. A central processing unit (CPU) 20 that controls the operation, a main memory 22 that stores a control program for the device, a magnetic disk 24 that stores a plurality of tomographic image data and programs, and image data for display are temporarily stored. A display memory 26 to be stored, a CRT monitor 28 as a display device for displaying an image based on image data from the display memory 26, a mouse 30 as a position input device, and a CRT monitor by detecting the state of the mouse 30 A mouse controller 32 for outputting signals such as the position of the mouse pointer on the mouse 18 and the state of the mouse 30 to the CPU 20; A board 34, and a common bus 26 for connecting the above components.
[0035]
In this embodiment, the pixel value on the projection plane is calculated based on Expression 1, and a 3D image is formed. However, the present invention is not limited to this, and the following Expression 2, Expression 3, etc. May be calculated.
[0036]
[Expression 2]
Figure 0004397131
Where C is a constant
[Equation 3]
Figure 0004397131
To increase the reflectivity and further enhance the degree of emphasis. Further, the above Expression 3 is attenuated by multiplying the incident intensity, and the attenuation of the incident intensity is slower than that of Expression 1, so that a 3D image that can be seen through to the inside can be formed.
[0038]
Further, although the present invention has been described by assigning the light emission rate and the light emission color, it includes combining them and highlighting a specific part or stent as a result of the combination.
[0039]
【The invention's effect】
As described above, according to the present invention, it is possible to construct and display a three-dimensional image in which a specific part is highlighted by the volume rendering method, and it is easy to find a lesion part at the time of a health check.
[Brief description of the drawings]
FIG. 1 is a diagram used for explaining shading by a volume rendering method.
FIG. 2 is a diagram used for explaining the light emission rate and the like according to the present invention.
FIG. 3 is a diagram showing an example of a display screen of a CRT display that displays a 3D image or the like created based on the present invention.
FIG. 4 is a diagram showing an example of a 3D image displayed by the present invention and a conventional 3D image.
FIG. 5 is a diagram showing another example of a 3D image displayed according to the present invention and a conventional 3D image.
FIG. 6 is a flowchart showing a 3D image display processing procedure according to the present invention.
FIG. 7 is a diagram illustrating an example of setting the light emission rate, the color to be emitted, and the like according to the present invention.
FIGS. 8A and 8B are diagrams illustrating an example of setting a light emission color according to the present invention and a display example thereof; FIGS.
FIG. 9 is a block diagram showing a hardware configuration of a 3D image display apparatus according to the present invention.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 ... Reflectance assigning part, 2 ... Luminescence assigning part 3, 10, 12, 14, 16, ... 3D image, 15, 17 ... Stent graft, 20 ... CPU, 22 ... Main memory, 24 ... Magnetic disk, 26 ... Display memory, 28 ... CRT monitor, 30 ... Mouse, 32 ... Controller, 34 ... Keyboard, 36 ... Common bus

Claims (3)

複数枚の断層像が積み上げられてなるボリューム画像を所定の投影面上に陰影付けして投影し、三次元画像を表示する三次元画像表示装置であって、前記陰影付けのアルゴリズムとして前記ボリューム画像の外に設定した仮想光源から光線を照射し、該光線の前記ボリューム画像の各画素での反射率及び濃度勾配に基づいて前記投影面上の画素値を算出するボリュームレンダリング法を用いた三次元画像表示装置において、前記ボリューム画像中の表示対象物のうちから強調表示する特定の対象物を示す画素に対して該画素に入射する光線を増幅させる発光率を割り当てる発光率割当手段と、前記発光率割当手段によって割り当てられた発光率をボリュームレンダリング法における陰影付け計算のパラメータとして追加し、投影面上の画素値を算出する演算手段と、を備え、前記発光率が割り当てられた特定の対象物を強調表示するようにしたことを特徴とする三次元画像表示装置。  A three-dimensional image display device for displaying a three-dimensional image by shading and projecting a volume image formed by stacking a plurality of tomographic images on a predetermined projection plane, wherein the volume image is used as the shading algorithm 3D using a volume rendering method that irradiates a light beam from a virtual light source set outside and calculates a pixel value on the projection plane based on a reflectance and a density gradient at each pixel of the volume image of the light beam In the image display device, a light emission rate assigning unit that assigns a light emission rate for amplifying a light beam incident on a pixel indicating a specific object to be highlighted among the display objects in the volume image, and the light emission The light emission rate assigned by the rate assigning means is added as a parameter for shading calculation in the volume rendering method, and the pixel value on the projection plane Comprising calculating means for calculating, for the emission ratios three-dimensional image display apparatus being characterized in that so as to highlight certain objects assigned. 前記発光率割当手段は、被検体の各部位の病巣ごとに、その病巣を示す画素の濃度値に対応して発光率が割り当てられた発光率パターンを有し、強調表示させたい部位の病巣を選択指定することにより、該部位の病巣に対応する発光率パターンを選択することを特徴とする請求項1の三次元画像表示装置。  The light emission rate assigning means has a light emission rate pattern in which a light emission rate is assigned to each lesion of each part of the subject corresponding to a density value of a pixel indicating the lesion, and a lesion of a part to be highlighted is displayed. 2. The three-dimensional image display device according to claim 1, wherein the light emission rate pattern corresponding to the lesion of the part is selected by selecting and specifying. 複数枚の断層像が積み上げられてなるボリューム画像を所定の投影面上に陰影付けして投影し、三次元画像を表示する三次元画像表示装置であって、前記陰影付けのアルゴリズムとして前記ボリューム画像の外に設定した仮想光源から光線を照射し、該光線の前記ボリューム画像の各画素での反射率及び濃度勾配に基づいて前記投影面上の画素値を算出するボリュームレンダリング法を用いた三次元画像表示装置において、前記ボリューム画像中の表示対象物のうちから強調表示する特定の対象物を示す画素に対して色発光率を割り当てる色発光率割当手段と、前記色発光率割当手段によって割り当てられた色発光率をボリュームレンダリング法における陰影付け計算のパラメータとして追加し投影面上の画素値を算出する演算手段とを備え、前記色発光率が割り当てられた特定の対象物を強調表示するようにしたことを特徴とする三次元画像表示装置。  A three-dimensional image display device for displaying a three-dimensional image by shading and projecting a volume image formed by stacking a plurality of tomographic images on a predetermined projection plane, wherein the volume image is used as the shading algorithm 3D using a volume rendering method that irradiates a light beam from a virtual light source set outside and calculates a pixel value on the projection plane based on a reflectance and a density gradient at each pixel of the volume image of the light beam In the image display device, a color light emission rate assigning unit for assigning a color light emission rate to a pixel indicating a specific target object to be highlighted among the display objects in the volume image, and the color light emission rate assigning unit. And calculating means for calculating the pixel value on the projection plane by adding the color emission rate as a parameter for shading calculation in the volume rendering method. Three-dimensional image display device is characterized in that so as to highlight the particular object in which the color light-emitting rate is allocated.
JP2001092944A 2000-04-03 2001-03-28 3D image display device Expired - Fee Related JP4397131B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2001092944A JP4397131B2 (en) 2000-04-03 2001-03-28 3D image display device
US10/471,548 US7019745B2 (en) 2001-03-28 2001-10-03 Three-dimensional image display device
PCT/JP2001/008716 WO2002078545A1 (en) 2001-03-28 2001-10-03 Three-dimensional image display device

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2000-101094 2000-04-03
JP2000101094 2000-04-03
JP2001092944A JP4397131B2 (en) 2000-04-03 2001-03-28 3D image display device

Publications (3)

Publication Number Publication Date
JP2001351120A JP2001351120A (en) 2001-12-21
JP2001351120A5 JP2001351120A5 (en) 2007-01-25
JP4397131B2 true JP4397131B2 (en) 2010-01-13

Family

ID=26589353

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001092944A Expired - Fee Related JP4397131B2 (en) 2000-04-03 2001-03-28 3D image display device

Country Status (1)

Country Link
JP (1) JP4397131B2 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6692441B1 (en) * 2002-11-12 2004-02-17 Koninklijke Philips Electronics N.V. System for identifying a volume of interest in a volume rendered ultrasound image
EP1985236A4 (en) 2006-02-17 2010-11-17 Hitachi Medical Corp Image display device and program
JP4716903B2 (en) * 2006-03-09 2011-07-06 株式会社日立メディコ Image display device and image processing program
JP5049654B2 (en) * 2006-05-29 2012-10-17 株式会社東芝 Medical image processing apparatus and medical image processing method
US8442282B2 (en) 2006-05-29 2013-05-14 Kabushiki Kaisha Toshiba Computer-aided imaging diagnostic processing apparatus and computer-aided imaging diagnostic processing method
JP6202757B2 (en) * 2012-09-26 2017-09-27 株式会社日立製作所 Ultrasonic diagnostic apparatus and ultrasonic two-dimensional tomographic image generation method
CN109934902B (en) * 2019-03-13 2020-08-11 南京大学 A Gradient Domain Rendering Image Reconstruction Method Using Scene Feature Constraints
EP4200812A1 (en) * 2020-08-20 2023-06-28 Koninklijke Philips N.V. Rendering of two dimensional data sets

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02239383A (en) * 1989-03-14 1990-09-21 Toshiba Corp Three dimensional picture display device
JP2938502B2 (en) * 1990-03-01 1999-08-23 株式会社東芝 Image display device
JP2651787B2 (en) * 1993-09-03 1997-09-10 日本アイ・ビー・エム株式会社 Interactive volume rendering method and apparatus
JP3570576B2 (en) * 1995-06-19 2004-09-29 株式会社日立製作所 3D image synthesis and display device compatible with multi-modality
JPH09259298A (en) * 1996-03-25 1997-10-03 Hitachi Medical Corp Pseudo three-dimensional image constitution method
JPH1021362A (en) * 1996-06-28 1998-01-23 Hitachi Medical Corp Method and device for measuring three-dimensional concerned area
JP2891350B2 (en) * 1997-03-04 1999-05-17 技術研究組合医療福祉機器研究所 3D image processing method
JPH1125287A (en) * 1997-07-08 1999-01-29 Toshiba Mach Co Ltd Method and device for setting voxel opacity
JP2000090283A (en) * 1998-09-09 2000-03-31 Toshiba Iyo System Engineering Kk Volume rendering image display method, image processor and storage medium storing program for the same method
US6356265B1 (en) * 1998-11-12 2002-03-12 Terarecon, Inc. Method and apparatus for modulating lighting with gradient magnitudes of volume data in a rendering pipeline

Also Published As

Publication number Publication date
JP2001351120A (en) 2001-12-21

Similar Documents

Publication Publication Date Title
JP6735016B2 (en) Medical image processor
US9330485B2 (en) Volume rendering of medical images
US10420534B2 (en) Ultrasonic diagnostic device, image processing device, and image processing method
US9741161B2 (en) Method, system, and medical imaging device for shading with multiple light sources
US7019745B2 (en) Three-dimensional image display device
JP2009034521A (en) System and method for volume rendering data in medical diagnostic imaging, and computer readable storage medium
JP6248044B2 (en) Visualization of 3D medical perfusion images
CN109887064A (en) Method and system for X-ray image generation
JP6563177B2 (en) Medical image processing apparatus and medical image processing program
US20170301129A1 (en) Medical image processing apparatus, medical image processing method, and medical image processing system
JP4397131B2 (en) 3D image display device
JP4122314B2 (en) Projection image processing method, projection image processing program, and projection image processing apparatus
JP2004215846A (en) Method and device for volume rendering image processing and program therefor
US10249074B2 (en) Medical image processing device, medical image processing method and computer readable medium for displaying color volume rendered images
JP2010131315A (en) Medical image processor and medical image processing program
US20060152502A1 (en) Method for depicting an object displayed in a volume data set
US20150320507A1 (en) Path creation using medical imaging for planning device insertion
US11259782B2 (en) Medical imaging data processing apparatus and method
JP4716903B2 (en) Image display device and image processing program
JP2018121857A (en) Medical image processing device, medical image processing method, and medical image processing program
JP7613995B2 (en) Ultrasound diagnostic device, image processing device, and program
JP6035057B2 (en) 3D image display apparatus and 3D image display method
WO2011062108A1 (en) Image processing apparatus and image processing method
JP2024000957A (en) Medical image processing device and medical image processing method
WO2024100186A1 (en) System for virtual reality-based preoperative surgical planning

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061204

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061204

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: 20091013

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: 20091020

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121030

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: 20121030

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131030

Year of fee payment: 4

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees