JP2001184491A - 三次元画像表示装置 - Google Patents
三次元画像表示装置Info
- Publication number
- JP2001184491A JP2001184491A JP37003499A JP37003499A JP2001184491A JP 2001184491 A JP2001184491 A JP 2001184491A JP 37003499 A JP37003499 A JP 37003499A JP 37003499 A JP37003499 A JP 37003499A JP 2001184491 A JP2001184491 A JP 2001184491A
- Authority
- JP
- Japan
- Prior art keywords
- image
- dimensional image
- value
- dimensional
- display device
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 abstract description 50
- 210000004204 blood vessel Anatomy 0.000 abstract description 11
- 210000000056 organ Anatomy 0.000 abstract description 5
- 230000002950 deficient Effects 0.000 abstract 1
- 230000015654 memory Effects 0.000 description 18
- 238000010586 diagram Methods 0.000 description 8
- 238000009877 rendering Methods 0.000 description 4
- 102100028780 AP-1 complex subunit sigma-2 Human genes 0.000 description 1
- 101100316860 Autographa californica nuclear polyhedrosis virus DA18 gene Proteins 0.000 description 1
- 102100035893 CD151 antigen Human genes 0.000 description 1
- 102100033040 Carbonic anhydrase 12 Human genes 0.000 description 1
- 101100321669 Fagopyrum esculentum FA02 gene Proteins 0.000 description 1
- 102100036738 Guanine nucleotide-binding protein subunit alpha-11 Human genes 0.000 description 1
- 101100055680 Homo sapiens AP1S2 gene Proteins 0.000 description 1
- 101000946874 Homo sapiens CD151 antigen Proteins 0.000 description 1
- 101000867855 Homo sapiens Carbonic anhydrase 12 Proteins 0.000 description 1
- 101100283445 Homo sapiens GNA11 gene Proteins 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010968 computed tomography angiography Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
- Image Generation (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】
【課題】断層像に基づいて所望臓器の三次元画像を再構
成する際に、その臓器のＣＴ値（しきい値の下限値、上
限値）のような三次元画像再構成パラメータを最適値に
することができる三次元画像表示装置を提供する。 【解決手段】三次元画像構成時にＣＴ値の下限値、上限
値をいくつに設定するかが問題で、範囲が広すぎると血
管以外の臓器も含まれることになり、範囲が狭すぎると
欠けた表示になる。そこで、重み付け投影方法により構
成した画像中の対象の特定濃度部分と前記陰影づけた三
次元画像中の前記対象の大きさを関係づけて三次元画像
再構成パラメータを決定する。
成する際に、その臓器のＣＴ値（しきい値の下限値、上
限値）のような三次元画像再構成パラメータを最適値に
することができる三次元画像表示装置を提供する。 【解決手段】三次元画像構成時にＣＴ値の下限値、上限
値をいくつに設定するかが問題で、範囲が広すぎると血
管以外の臓器も含まれることになり、範囲が狭すぎると
欠けた表示になる。そこで、重み付け投影方法により構
成した画像中の対象の特定濃度部分と前記陰影づけた三
次元画像中の前記対象の大きさを関係づけて三次元画像
再構成パラメータを決定する。
Description
【０００１】
【発明の属する技術分野】本発明は三次元画像表示装置
に係り、特に所望の臓器部分の陰影付け三次元画像を再
構成し、表示する三次元画像表示装置に関する。
に係り、特に所望の臓器部分の陰影付け三次元画像を再
構成し、表示する三次元画像表示装置に関する。
【０００２】
【従来の技術】Ｘ線ＣＴによって収集した断層像（ＣＴ
画像）から三次元画像を再構成する方法として、代表的
なものにサーフェイス法とボリュームレンダリング法が
知られているが、例えば、ＣＴ画像からサーフェイス法
でどの臓器部分の三次元画像を再構成するかを決めるに
は、その関心対象のＣＴ値のしきい値を指定する必要が
ある。例えば、血管部のＣＴ値を含むようなしきい値
（下限値＝ｔｌ、上限値＝ｔ２）を指定するとＣＴ画像
から血管部の三次元画像が再構成できる。
画像）から三次元画像を再構成する方法として、代表的
なものにサーフェイス法とボリュームレンダリング法が
知られているが、例えば、ＣＴ画像からサーフェイス法
でどの臓器部分の三次元画像を再構成するかを決めるに
は、その関心対象のＣＴ値のしきい値を指定する必要が
ある。例えば、血管部のＣＴ値を含むようなしきい値
（下限値＝ｔｌ、上限値＝ｔ２）を指定するとＣＴ画像
から血管部の三次元画像が再構成できる。
【０００３】尚、 サーフェイス法に関する詳細につい
ては、“Ｓｈａｄｉｎｇ ３Ｄ−Ｉｍａｇｅｓ ｆｒｏ
ｍ ＣＴ Ｕｓｉｎｇ Ｇｒａｙ−Ｌｅｖｅｌ Ｇｒａ
ｄｉｅｎｔｓ” ｂｙ ＫＡＲＬ ＨＥＩＮＺ ＨＯＮ
Ｅ ａｎｄ ＲＡＰＨ ＢＥＲＮＳＴＥＩＮ；ＩＥＥＥ
ＴＲＡＮＳＡＣＴＩＯＮ ＯＮ ＭＥＤＩＣＡＬ Ｉ
ＭＡＧＩＮＧ，ＶＯＬ．ＭＩ−５，ＮＯ，１ＭＡＲＣＨ
Ｉ９８６ ４５−４７ページに記載されている。
ては、“Ｓｈａｄｉｎｇ ３Ｄ−Ｉｍａｇｅｓ ｆｒｏ
ｍ ＣＴ Ｕｓｉｎｇ Ｇｒａｙ−Ｌｅｖｅｌ Ｇｒａ
ｄｉｅｎｔｓ” ｂｙ ＫＡＲＬ ＨＥＩＮＺ ＨＯＮ
Ｅ ａｎｄ ＲＡＰＨ ＢＥＲＮＳＴＥＩＮ；ＩＥＥＥ
ＴＲＡＮＳＡＣＴＩＯＮ ＯＮ ＭＥＤＩＣＡＬ Ｉ
ＭＡＧＩＮＧ，ＶＯＬ．ＭＩ−５，ＮＯ，１ＭＡＲＣＨ
Ｉ９８６ ４５−４７ページに記載されている。
【０００４】また、ボリュームレンダリング法に関する
詳細については、“三次元ＣＴ診断”隈崎達夫、他編著
１９９５年 （株）南江堂 発行 ８ページに記載さ
れている。
詳細については、“三次元ＣＴ診断”隈崎達夫、他編著
１９９５年 （株）南江堂 発行 ８ページに記載さ
れている。
【０００５】
【発明が解決しようとする課題】しかしながら、三次元
画像として表示したい関心対象のＣＴ値（しきい値の下
限値、上限値）をいくつに設定するかが問題で、範囲が
広すぎるとその対象（例えば、血管）以外の臓器も含ま
れることになり（図３のｂ−２やｃ−２の４１ａ〜４１
ｃ、４２ａ〜４２ｅ）、範囲が狭すぎると欠けた表示に
なる。
画像として表示したい関心対象のＣＴ値（しきい値の下
限値、上限値）をいくつに設定するかが問題で、範囲が
広すぎるとその対象（例えば、血管）以外の臓器も含ま
れることになり（図３のｂ−２やｃ−２の４１ａ〜４１
ｃ、４２ａ〜４２ｅ）、範囲が狭すぎると欠けた表示に
なる。
【０００６】本発明はこのような事情に鑑みてなされた
もので、前記しきい値の下限値、上限値のような関心対
象の陰影付け三次元画像を再構成するための三次元画像
再構成パラメータを最適な値に設定して、関心対象の陰
影付け三次元画像を適切に表示することができる三次元
画像表示装置を提供することを目的とする。
もので、前記しきい値の下限値、上限値のような関心対
象の陰影付け三次元画像を再構成するための三次元画像
再構成パラメータを最適な値に設定して、関心対象の陰
影付け三次元画像を適切に表示することができる三次元
画像表示装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】本発明は前記目的を解決
するために、複数の二次元断層像及び複数の二次元断層
像に分解可能な三次元原データから三次元画像再構成パ
ラメータに基づいて関心対象を陰影付けて三次元画像を
再構成し、表示器に表示する三次元画像表示装置におい
て、上記再構成した画像中の対象の特定濃度部分と前記
再構成した三次元画像中の前記対象の大きさを設定する
手段と、該設定した特定濃度部分と対象の大きさを関係
付けて前記三次元画像再構成パラメータを決定する手段
とを備え、該パラメータを用いて陰影付け三次元画像を
再構成し、前記表示器に表示することを特徴としてい
る。
するために、複数の二次元断層像及び複数の二次元断層
像に分解可能な三次元原データから三次元画像再構成パ
ラメータに基づいて関心対象を陰影付けて三次元画像を
再構成し、表示器に表示する三次元画像表示装置におい
て、上記再構成した画像中の対象の特定濃度部分と前記
再構成した三次元画像中の前記対象の大きさを設定する
手段と、該設定した特定濃度部分と対象の大きさを関係
付けて前記三次元画像再構成パラメータを決定する手段
とを備え、該パラメータを用いて陰影付け三次元画像を
再構成し、前記表示器に表示することを特徴としてい
る。
【０００８】本発明によれば、重み付け投影方法又はＭ
ＩＰ法により再構成した画像中の対象の特定濃度部分に
基づいて三次元画像再構成パラメータを決定するように
したため、三次元画像再構成パラメータを最適な値に設
定することができる。
ＩＰ法により再構成した画像中の対象の特定濃度部分に
基づいて三次元画像再構成パラメータを決定するように
したため、三次元画像再構成パラメータを最適な値に設
定することができる。
【０００９】重み付け投影方法は、例えば、特願平１０
−５１７３６８号に記載されている。図５に示すよう
に、画面サイズと同一サイズの分類メモリを、観察した
いＣＴ値の範囲の枚数だけ用意し、即ち（ＣＴ値−１０
００用メモリ〜ＣＴ値２０００用メモリの枚数）×画面
サイズの三次元的メモリを用意し、投影線上で分類メモ
リのＣＴ値をとる画像の有無をその分類メモリに記録し
ている。そして、上記メモリに記憶された画像の有無情
報と分類されたＣＴ値ごとの重み係数とを積和演算し、
その演算値に応じた濃度値を表示メモリ（係数加算メモ
リ）に記録することにより、重み係数の大きいＣＴ値を
有する特定部位が強調された透視像を表示する。
−５１７３６８号に記載されている。図５に示すよう
に、画面サイズと同一サイズの分類メモリを、観察した
いＣＴ値の範囲の枚数だけ用意し、即ち（ＣＴ値−１０
００用メモリ〜ＣＴ値２０００用メモリの枚数）×画面
サイズの三次元的メモリを用意し、投影線上で分類メモ
リのＣＴ値をとる画像の有無をその分類メモリに記録し
ている。そして、上記メモリに記憶された画像の有無情
報と分類されたＣＴ値ごとの重み係数とを積和演算し、
その演算値に応じた濃度値を表示メモリ（係数加算メモ
リ）に記録することにより、重み係数の大きいＣＴ値を
有する特定部位が強調された透視像を表示する。
【００１０】図３左側のａ−１，ｂ−１，ｃ−１は重み
付け投影方法により構成した血管の画像を示し、右側の
ａ−２，ｂ−２，ｃ−２は例えばサーフェイス法により
構成した陰影付け三次元画像を示す。ａ−２は最適なし
きい値を設定した場合の陰影付け三次元画像であり、血
管の太さはどちらも正しく描出されている。
付け投影方法により構成した血管の画像を示し、右側の
ａ−２，ｂ−２，ｃ−２は例えばサーフェイス法により
構成した陰影付け三次元画像を示す。ａ−２は最適なし
きい値を設定した場合の陰影付け三次元画像であり、血
管の太さはどちらも正しく描出されている。
【００１１】右側のａ−２，ｂ−２，ｃ−２は陰影付け
三次元画像では、一般的には、しきい値設定の範囲が正
しくないと、ｂ−２，ｃ−２に示すように正しくない部
分（４１ａ〜４１ｃ、４２ａ〜４２ｅ）が描出されてし
まう。
三次元画像では、一般的には、しきい値設定の範囲が正
しくないと、ｂ−２，ｃ−２に示すように正しくない部
分（４１ａ〜４１ｃ、４２ａ〜４２ｅ）が描出されてし
まう。
【００１２】本発明では、図３に示す一つの実験事実に
注目する。図３左側のａ−１，ｂ‐１，ｃ−１は重み付
け投影方法により構成した画像でも正しくない部分は現
れる（４１ａ〜４１ｃ、４２ａ〜４２ｅ）が、血管の正
しい輪郭は依然として高濃度（太線部分）のまま保存さ
れる。本発明は、この設定パラメータによる影響の少な
い高濃度部分を利用する。重み付け投影方法により構成
した画像中の対象の特定濃度部分と前記陰影付けた三次
元画像中の前記対象の大きさを関係づけて陰影づけ三次
元画像再構成パラメータを決定する。
注目する。図３左側のａ−１，ｂ‐１，ｃ−１は重み付
け投影方法により構成した画像でも正しくない部分は現
れる（４１ａ〜４１ｃ、４２ａ〜４２ｅ）が、血管の正
しい輪郭は依然として高濃度（太線部分）のまま保存さ
れる。本発明は、この設定パラメータによる影響の少な
い高濃度部分を利用する。重み付け投影方法により構成
した画像中の対象の特定濃度部分と前記陰影付けた三次
元画像中の前記対象の大きさを関係づけて陰影づけ三次
元画像再構成パラメータを決定する。
【００１３】例えば、重み付け投影方法により構成した
画像中の対象の輪郭に、サーフェイス三次元画像の輪郭
が一致するかどうか（目視のより判定してもよく、コン
ピュータ演算により両者の差をとり、差が最小のとき一
致すると判定してもよい）でしきい値を決定する。
画像中の対象の輪郭に、サーフェイス三次元画像の輪郭
が一致するかどうか（目視のより判定してもよく、コン
ピュータ演算により両者の差をとり、差が最小のとき一
致すると判定してもよい）でしきい値を決定する。
【００１４】
【発明の実施の形態】以下添付図面に従って本発明に係
る三次元画像表示装置の好ましい実施の形態について説
明する。
る三次元画像表示装置の好ましい実施の形態について説
明する。
【００１５】尚、本実施の形態ではしきい値のとりかた
で表示する対象の大きさが大きく影響をうける、サーフ
ェイス法に基づく陰影付け三次元画像を再構成する場合
を例に説明する。
で表示する対象の大きさが大きく影響をうける、サーフ
ェイス法に基づく陰影付け三次元画像を再構成する場合
を例に説明する。
【００１６】陰影付け三次元画像は、Ｘ線ＣＴ装置から
取得された複数枚の断層像が積み上げられてなる三次元
原画像に基づいて再構成されるが、血管等の所要臓器を
表示対象とする場合にはその表示対象に適合したパラメ
ータ（本実施の形態では表示対象のＣＴ値のしきい値）
を設定することが必要である。尚、断層像はＸ線ＣＴ装
置から取得する場合に限らず、ＭＲＩ装置、あるいは超
音波装置等から取得する場合であっても本発明を同様に
適用できる。
取得された複数枚の断層像が積み上げられてなる三次元
原画像に基づいて再構成されるが、血管等の所要臓器を
表示対象とする場合にはその表示対象に適合したパラメ
ータ（本実施の形態では表示対象のＣＴ値のしきい値）
を設定することが必要である。尚、断層像はＸ線ＣＴ装
置から取得する場合に限らず、ＭＲＩ装置、あるいは超
音波装置等から取得する場合であっても本発明を同様に
適用できる。
【００１７】そこで、本発明では、重み付け投影方法に
より構成した画像中の対象の特定濃度部分と前記陰影付
け三次元画像中の前記対象の大きさが一致するように三
次元画像再構成パラメータを設定する。一致するかどう
かの判定には、陰影付け三次元画像中の前記対象の輪郭
部分と特定濃度部分の積和をとって、最大値をとった状
態を一致したと判定するなどの方法をとればよい。
より構成した画像中の対象の特定濃度部分と前記陰影付
け三次元画像中の前記対象の大きさが一致するように三
次元画像再構成パラメータを設定する。一致するかどう
かの判定には、陰影付け三次元画像中の前記対象の輪郭
部分と特定濃度部分の積和をとって、最大値をとった状
態を一致したと判定するなどの方法をとればよい。
【００１８】図１のフローチャートを参照しながら本発
明における陰影付け三次元画像の再構成方法について説
明する。 〔ステップＳ１〕重み付け投影方法におけるしきい値の
上限値ｔ２、下限値ｔ１として、表示対象（本実施の形
態では血管）に応じたデフォルト値を磁気ディスクなど
から読み込んで設定する。 〔ステップＳ２〕重み付け投影方法により、重み付け投
影画像を作成する。作成された重み付け投影画像は図２
のｃ−１で示すようにＣＲＴモニタ１８に表示される。 〔ステップＳ３〕作成画像（重み付け投影画像）中の陰
影の高濃度部分（図２のｃ−１の太線部）を抽出する。
濃度一定値以上の領域をしきい値処理などにより抽出す
る。 〔ステップＳ４〕サーフェイス法におけるしきい値とし
てデフォルト値を磁気ディスクなどから読み込んで設定
し、サーフェイス三次元画像を構成する。作成されたサ
ーフェイス三次元画像は、図２のｃ−２で示すようにＣ
ＲＴモニタ１８に表示される。 〔ステップＳ５〕サーフェイス三次元画像のエッジを抽
出する。濃度一定値以上の領域をしきい値処理などによ
り抽出する。 〔ステップＳ６〕上記ステップＳ３で抽出した陰影の高
濃度部分の抽出部分と上記ステップＳ５で抽出したサー
フェイス三次元画像のエッジ抽出部分とを比較する。例
えば、おのおのの画素間で引き算して絶対値をとる。 〔ステップＳ７〕ステップＳ６の比較の結果、陰影の高
濃度部分の抽出部分とサーフェイス三次元画像のエッジ
抽出部分とが一致するかどうかを判定する。１画像すべ
ての画素についての上記絶対値の和が、所定の値より小
さいかどうかを判定する。目視で判定してもよい。小さ
ければ、一致したとしてステップＳ１０に跳ぶ。小さく
なければステップＳ８に行く。 〔ステップＳ８〕しきい値を変更する。自動的に一定値
を加算／減算してもよく、キー入力によってもよい。 〔ステップＳ９〕ステップＳ８のしきい値に基づいてサ
ーフェイス三次元画像を構成する。その後ステップＳ５
に跳ぶ。 〔ステップＳ１０〕サーフェイス三次元画像を表示し、
その画像と再構成パラメータ（しきい値）を磁気ディス
クに格納する。重み係数の下限値、上限値は、図２
（ｂ）のようになだらかに変化するものであってもよ
い。重み係数を図２（ｂ）のようにした場合でも本質的
なかわりはないが、図４に追加説明する。この場合、ａ
とｂ，ｃとｄの比は磁気ディスクなどに予め書き込んで
おいたものを読み込んでもよく、キー入力してもよい。 〔ステップＳ４１〕重み係数を設定する。 〔ステップＳ４２〕重み付け投影画像を構成する。 〔ステップＳ４３〕陰影の高濃度部分の抽出部分とエッ
ジ抽出部分は一致するかどうかを判定する。１画像すべ
ての画素についての絶対値の和が、所定の値より小さい
かどうかを判定する。目視で判定してもよい。一致すれ
ばステップＳ４５に跳ぶ。 〔ステップＳ４４〕重み係数を変更する。自動でもマウ
ス操作による手動でもよい。 〔ステップＳ４５〕重み係数設定パラメータから三次元
画像再構成パラメータを決定する。 〔ステップＳ４６〕上記三次元画像再構成パラメータに
基づいて三次元画像を再構成、表示し、パラメータとと
もに磁気ディスクに格納する。
明における陰影付け三次元画像の再構成方法について説
明する。 〔ステップＳ１〕重み付け投影方法におけるしきい値の
上限値ｔ２、下限値ｔ１として、表示対象（本実施の形
態では血管）に応じたデフォルト値を磁気ディスクなど
から読み込んで設定する。 〔ステップＳ２〕重み付け投影方法により、重み付け投
影画像を作成する。作成された重み付け投影画像は図２
のｃ−１で示すようにＣＲＴモニタ１８に表示される。 〔ステップＳ３〕作成画像（重み付け投影画像）中の陰
影の高濃度部分（図２のｃ−１の太線部）を抽出する。
濃度一定値以上の領域をしきい値処理などにより抽出す
る。 〔ステップＳ４〕サーフェイス法におけるしきい値とし
てデフォルト値を磁気ディスクなどから読み込んで設定
し、サーフェイス三次元画像を構成する。作成されたサ
ーフェイス三次元画像は、図２のｃ−２で示すようにＣ
ＲＴモニタ１８に表示される。 〔ステップＳ５〕サーフェイス三次元画像のエッジを抽
出する。濃度一定値以上の領域をしきい値処理などによ
り抽出する。 〔ステップＳ６〕上記ステップＳ３で抽出した陰影の高
濃度部分の抽出部分と上記ステップＳ５で抽出したサー
フェイス三次元画像のエッジ抽出部分とを比較する。例
えば、おのおのの画素間で引き算して絶対値をとる。 〔ステップＳ７〕ステップＳ６の比較の結果、陰影の高
濃度部分の抽出部分とサーフェイス三次元画像のエッジ
抽出部分とが一致するかどうかを判定する。１画像すべ
ての画素についての上記絶対値の和が、所定の値より小
さいかどうかを判定する。目視で判定してもよい。小さ
ければ、一致したとしてステップＳ１０に跳ぶ。小さく
なければステップＳ８に行く。 〔ステップＳ８〕しきい値を変更する。自動的に一定値
を加算／減算してもよく、キー入力によってもよい。 〔ステップＳ９〕ステップＳ８のしきい値に基づいてサ
ーフェイス三次元画像を構成する。その後ステップＳ５
に跳ぶ。 〔ステップＳ１０〕サーフェイス三次元画像を表示し、
その画像と再構成パラメータ（しきい値）を磁気ディス
クに格納する。重み係数の下限値、上限値は、図２
（ｂ）のようになだらかに変化するものであってもよ
い。重み係数を図２（ｂ）のようにした場合でも本質的
なかわりはないが、図４に追加説明する。この場合、ａ
とｂ，ｃとｄの比は磁気ディスクなどに予め書き込んで
おいたものを読み込んでもよく、キー入力してもよい。 〔ステップＳ４１〕重み係数を設定する。 〔ステップＳ４２〕重み付け投影画像を構成する。 〔ステップＳ４３〕陰影の高濃度部分の抽出部分とエッ
ジ抽出部分は一致するかどうかを判定する。１画像すべ
ての画素についての絶対値の和が、所定の値より小さい
かどうかを判定する。目視で判定してもよい。一致すれ
ばステップＳ４５に跳ぶ。 〔ステップＳ４４〕重み係数を変更する。自動でもマウ
ス操作による手動でもよい。 〔ステップＳ４５〕重み係数設定パラメータから三次元
画像再構成パラメータを決定する。 〔ステップＳ４６〕上記三次元画像再構成パラメータに
基づいて三次元画像を再構成、表示し、パラメータとと
もに磁気ディスクに格納する。
【００１９】以上の説明のように、陰影づけた三次元画
像の再構成時のパラメータ決定は、操作者が画面を見な
がら決定してもよく自動的決定してもよい。“合成”ボ
タンがクリックされたらｃ−１の高濃度抽出部分とｃ−
２を重ね合わせて表示すれば、一致度は明瞭になる。重
み付け投影方法により構成した画像中の対象の特定濃度
部分を参考にして陰影づけた三次元画像の再構成時のパ
ラメータを決定可能とする装置は、全て本発明に含まれ
る。
像の再構成時のパラメータ決定は、操作者が画面を見な
がら決定してもよく自動的決定してもよい。“合成”ボ
タンがクリックされたらｃ−１の高濃度抽出部分とｃ−
２を重ね合わせて表示すれば、一致度は明瞭になる。重
み付け投影方法により構成した画像中の対象の特定濃度
部分を参考にして陰影づけた三次元画像の再構成時のパ
ラメータを決定可能とする装置は、全て本発明に含まれ
る。
【００２０】次に、重み付け投影画像の高濃度部分の抽
出方法の一例について説明する。図６は、その説明図で
ある。同図に示すように重み付け投影画像とＭＩＰ画像
とを構成し、メモリ１とメモリ２に格納する。尚、ＭＩ
Ｐ画像は、最大値投影法によるもので、ＭＩＰの詳細に
ついては、“三次元ＣＴ血管造影法”永井純編著 １９
９４年 （株）メディカルサイエンスインターナショナ
ル発行 ４７ページに記載されている。
出方法の一例について説明する。図６は、その説明図で
ある。同図に示すように重み付け投影画像とＭＩＰ画像
とを構成し、メモリ１とメモリ２に格納する。尚、ＭＩ
Ｐ画像は、最大値投影法によるもので、ＭＩＰの詳細に
ついては、“三次元ＣＴ血管造影法”永井純編著 １９
９４年 （株）メディカルサイエンスインターナショナ
ル発行 ４７ページに記載されている。
【００２１】そして、これらメモリ１、２に格納した画
像の対応する画素同志で、（重み付け投影画像）×Ｃ１
−（ＭＩＰ画像）×Ｃ２の演算を行う。Ｃ１、Ｃ２は定
数である。（重み付け投影画像）×Ｃ１−（ＭＩＰ画
像）×Ｃ２の正値部分は高濃度部分を含んでいる。この
演算結果を表示メモリに判定画像として格納してＣＲＴ
モニタ１８に表示する。これにより、操作者が画面を見
ながらマウス２０等を使用して手動でＣ２を好適な値に
決定する。
像の対応する画素同志で、（重み付け投影画像）×Ｃ１
−（ＭＩＰ画像）×Ｃ２の演算を行う。Ｃ１、Ｃ２は定
数である。（重み付け投影画像）×Ｃ１−（ＭＩＰ画
像）×Ｃ２の正値部分は高濃度部分を含んでいる。この
演算結果を表示メモリに判定画像として格納してＣＲＴ
モニタ１８に表示する。これにより、操作者が画面を見
ながらマウス２０等を使用して手動でＣ２を好適な値に
決定する。
【００２２】画面では、判定画素＞０の時、判定画素を
カラー表示、又は、判定画素＞０の時、対応する重み付
け投影画像の画素をカラー表示、又は、判定画素＞０の
時、対応するＭＩＰ画像の画素をカラー表示、にするこ
とにより操作者に分かりやくす表示してもよい。
カラー表示、又は、判定画素＞０の時、対応する重み付
け投影画像の画素をカラー表示、又は、判定画素＞０の
時、対応するＭＩＰ画像の画素をカラー表示、にするこ
とにより操作者に分かりやくす表示してもよい。
【００２３】上記Ｃ２を手動ではなく、自動的に決定す
る場合について図７のフローチャートを用いて説明す
る。 〔ステップＳ８０〕Ｃ２を予め決めておいた初期値１、
ｄＣを予め決めておいた初期値２とする。 〔ステップＳ８２〕図６の計算を行う。 〔ステップＳ８４〕Ｒ＝カラー画素数／全画素数を求め
る。 〔ステップＳ８６〕Ｒ＜下限値か否かを判定する。ＹＥ
Ｓと判定した場合にはステップＳ８８に移行し、ＮＯと
判定した場合にはステップＳ９４に移行する。 〔ステップＳ８８〕フラグが−１か否かを判定する。Ｙ
ＥＳと判定した場合はステップＳ９０に移行し、ＮＯと
判定した場合にはステップＳ９２に移行する。 〔ステップＳ９０〕収束させるため、変化分ｄＣをｄＣ
／２としてｄＣを小さくする。 〔ステップＳ９２〕フラグを１とし、Ｃ２をＣ２−ｄＣ
とする。そして、上記ステップＳ８２に戻る。 〔ステップＳ９４〕Ｒ＞上限値か否かを判定する。ＹＥ
Ｓと判定した場合にはステップＳ９６に移行し、ＮＯと
判定した場合には終了する。 〔ステップＳ９６〕フラグが１か否かを判定する。ＹＥ
Ｓと判定した場合はステップＳ９８に移行し、ＮＯと判
定した場合にはステップＳ１００に移行する。 〔ステップＳ９８〕ｄＣをｄＣ／２とする。 〔ステップＳ１００〕フラグを−１とし、Ｃ２をＣ２＋
ｄＣとする。そして上記ステップＳ８２に戻る。
る場合について図７のフローチャートを用いて説明す
る。 〔ステップＳ８０〕Ｃ２を予め決めておいた初期値１、
ｄＣを予め決めておいた初期値２とする。 〔ステップＳ８２〕図６の計算を行う。 〔ステップＳ８４〕Ｒ＝カラー画素数／全画素数を求め
る。 〔ステップＳ８６〕Ｒ＜下限値か否かを判定する。ＹＥ
Ｓと判定した場合にはステップＳ８８に移行し、ＮＯと
判定した場合にはステップＳ９４に移行する。 〔ステップＳ８８〕フラグが−１か否かを判定する。Ｙ
ＥＳと判定した場合はステップＳ９０に移行し、ＮＯと
判定した場合にはステップＳ９２に移行する。 〔ステップＳ９０〕収束させるため、変化分ｄＣをｄＣ
／２としてｄＣを小さくする。 〔ステップＳ９２〕フラグを１とし、Ｃ２をＣ２−ｄＣ
とする。そして、上記ステップＳ８２に戻る。 〔ステップＳ９４〕Ｒ＞上限値か否かを判定する。ＹＥ
Ｓと判定した場合にはステップＳ９６に移行し、ＮＯと
判定した場合には終了する。 〔ステップＳ９６〕フラグが１か否かを判定する。ＹＥ
Ｓと判定した場合はステップＳ９８に移行し、ＮＯと判
定した場合にはステップＳ１００に移行する。 〔ステップＳ９８〕ｄＣをｄＣ／２とする。 〔ステップＳ１００〕フラグを−１とし、Ｃ２をＣ２＋
ｄＣとする。そして上記ステップＳ８２に戻る。
【００２４】以上、上記実施の形態では表示対象を血管
として説明した本発明において表示対象は血管に限定さ
れない。
として説明した本発明において表示対象は血管に限定さ
れない。
【００２５】また、陰影づけ三次元画像再構成パラメー
タの決定時は特定濃度部分と前記陰影付け三次元画像中
の前記対象の大きさが一致するかどうかを判定したが、
特定濃度部分と“最終的な”陰影付け三次元画像中の前
記対象の輪郭部分の大きさは一致する必要はない。表示
時に、経験的なｒ値（ｒ＝０．９〜１．１）でｒ倍に拡
大／縮小して表示することもある。
タの決定時は特定濃度部分と前記陰影付け三次元画像中
の前記対象の大きさが一致するかどうかを判定したが、
特定濃度部分と“最終的な”陰影付け三次元画像中の前
記対象の輪郭部分の大きさは一致する必要はない。表示
時に、経験的なｒ値（ｒ＝０．９〜１．１）でｒ倍に拡
大／縮小して表示することもある。
【００２６】また、重み係数の下限値・上限値をそれぞ
れしきい値の下限値・上限値としてもよい。ｙ＝ａｘ＋
ｂなどの一次式で関係づけてもよい。
れしきい値の下限値・上限値としてもよい。ｙ＝ａｘ＋
ｂなどの一次式で関係づけてもよい。
【００２７】また、これまでの説明ではサーフェイス画
像を対象にしたが、本発明はこれに限定されない。ボリ
ュームレンダリング法に適用し、重み係数の下限値また
は上限値をオパシティ最大の点としてもよい。重み係数
の分布形とオパシティの分布形を相似にしてもよい。
像を対象にしたが、本発明はこれに限定されない。ボリ
ュームレンダリング法に適用し、重み係数の下限値また
は上限値をオパシティ最大の点としてもよい。重み係数
の分布形とオパシティの分布形を相似にしてもよい。
【００２８】また、本発明は、ボリュームレンダリング
に適用し、半透明にした画像に対しても有効である。２
物体が半透明で同時表示されており一見比較できないと
思えるような場合でも、さらに重み付け投影方法により
構成した画像を半透明合成すれば、オパシティの分布形
を変えつつ表示すれば操作者により比較可能となる。
に適用し、半透明にした画像に対しても有効である。２
物体が半透明で同時表示されており一見比較できないと
思えるような場合でも、さらに重み付け投影方法により
構成した画像を半透明合成すれば、オパシティの分布形
を変えつつ表示すれば操作者により比較可能となる。
【００２９】また、重みづけ投影方法による画像をもと
にパラメータを決定したが、ＭＩＰ（最大値投影方法）
の輪郭をもとにパラメータを決定してもよい。
にパラメータを決定したが、ＭＩＰ（最大値投影方法）
の輪郭をもとにパラメータを決定してもよい。
【００３０】図８は、本発明に係る三次元画像表示装置
のハードウェア構成を示すブロック図である。この三次
元画像表示装置は、例えばコーンビームＣＴ装置で被検
体の対象部位について収集した複数の断層像により陰影
付け三次元画像を再構成して表示するもので、各構成要
素の動作を制御する中央処理装置（ＣＰＵ）１０と、装
置の制御プログラムが格納された主メモリ１２と、複数
の断層像データ及びプログラム等が格納された磁気ディ
スク１４と、表示用の画像データ（濃度値）を一時記憶
する表示メモリ１６と、この表示メモリ１６からの画像
データに基づいて画像を表示する表示装置としてのＣＲ
Ｔモニタ１８と、位置入力装置としてのマウス２０と、
マウス２０の状態を検出してＣＲＴモニタ１８上のマウ
スポインタの位置やマウス２０の状態等の信号をＣＰＵ
１０に出力するマウスコントローラ２２と、各種の操作
指令等を入力するキーボード２４と、上記各構成要素を
接続する共通バス２６とから構成される。
のハードウェア構成を示すブロック図である。この三次
元画像表示装置は、例えばコーンビームＣＴ装置で被検
体の対象部位について収集した複数の断層像により陰影
付け三次元画像を再構成して表示するもので、各構成要
素の動作を制御する中央処理装置（ＣＰＵ）１０と、装
置の制御プログラムが格納された主メモリ１２と、複数
の断層像データ及びプログラム等が格納された磁気ディ
スク１４と、表示用の画像データ（濃度値）を一時記憶
する表示メモリ１６と、この表示メモリ１６からの画像
データに基づいて画像を表示する表示装置としてのＣＲ
Ｔモニタ１８と、位置入力装置としてのマウス２０と、
マウス２０の状態を検出してＣＲＴモニタ１８上のマウ
スポインタの位置やマウス２０の状態等の信号をＣＰＵ
１０に出力するマウスコントローラ２２と、各種の操作
指令等を入力するキーボード２４と、上記各構成要素を
接続する共通バス２６とから構成される。
【００３１】
【発明の効果】以上説明したように本発明に係る三次元
画像表示装置によれば、血管等の所望の対象の陰影付け
三次元画像を最適な太さで再構成できる。
画像表示装置によれば、血管等の所望の対象の陰影付け
三次元画像を最適な太さで再構成できる。
【００３２】また、上記再構成法として重み付け投影方
法とＭＩＰ法を挙げたが、例えばデプス法等でも再構成
した画像の特定濃度と対象の大きさを設定できる再構成
法であればよい。
法とＭＩＰ法を挙げたが、例えばデプス法等でも再構成
した画像の特定濃度と対象の大きさを設定できる再構成
法であればよい。
【図１】本発明における陰影付け三次元画像の再構成手
順を示したフローチャート。
順を示したフローチャート。
【図２】（ａ）は本発明に係る三次元画像表示装置の表
示画面例、（ｂ）は重み係数の例を示した図。
示画面例、（ｂ）は重み係数の例を示した図。
【図３】実験事実を示した図であり、パラメータによる
画像の変化例を示した図。
画像の変化例を示した図。
【図４】重み係数から三次元画像再構成パラメータを決
定する場合のフローチャート。
定する場合のフローチャート。
【図５】重み付け投影法の説明図。
【図６】高濃度部分の抽出方法の説明図。
【図７】高濃度部分の抽出方法の説明に用いたフローチ
ャート。
ャート。
【図８】本発明に係る三次元画像表示装置のハードウエ
ア構成図。
ア構成図。
１０…ＣＰＵ、１２…主メモリ、１４…磁気ディスク、
１６…表示メモリ、１８…ＣＲＴモニタ、２０…マウ
ス、２４…キーボード
１６…表示メモリ、１８…ＣＲＴモニタ、２０…マウ
ス、２４…キーボード
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考） Ｇ０６Ｔ 15/00 Ｇ０６Ｆ 15/72 ４５０Ｋ ５Ｂ０８０ Ｆターム(参考） 4C093 CA01 FD09 FD12 FF34 FF43 4C096 AB01 DC33 DC36 4C301 EE07 JB27 JC12 KK12 KK16 5B050 AA02 BA03 BA09 DA02 EA30 FA02 5B057 AA09 BA03 CA02 CA08 CA12 CB02 CB08 CB13 DA16 DB03 DB05 DB09 DC22 5B080 AA19 BA03 GA11
Claims (1)
- 【請求項１】複数の二次元断層像及び複数の二次元断層
像に分解可能な三次元原データから三次元画像再構成パ
ラメータに基づいて関心対象を陰影付けて三次元画像を
再構成し、表示器に表示する三次元画像表示装置におい
て、 上記再構成した画像中の対象の特定濃度部分と前記再構
成した三次元画像中の前記対象の大きさを設定する手段
と、該設定した特定濃度部分と対象の大きさを関係付け
て前記三次元画像再構成パラメータを決定する手段とを
備え、該パラメータを用いて陰影付け三次元画像を再構
成し、前記表示器に表示することを特徴とする三次元画
像表示装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP37003499A JP2001184491A (ja) | 1999-12-27 | 1999-12-27 | 三次元画像表示装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP37003499A JP2001184491A (ja) | 1999-12-27 | 1999-12-27 | 三次元画像表示装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2001184491A true JP2001184491A (ja) | 2001-07-06 |
Family
ID=18495914
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP37003499A Pending JP2001184491A (ja) | 1999-12-27 | 1999-12-27 | 三次元画像表示装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2001184491A (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005261531A (ja) * | 2004-03-17 | 2005-09-29 | Hitachi Medical Corp | 医用画像表示方法および装置 |
| JP2006075390A (ja) * | 2004-09-10 | 2006-03-23 | Toshiba Corp | 医用画像表示装置 |
| CN102018524A (zh) * | 2009-09-09 | 2011-04-20 | 上海西门子医疗器械有限公司 | 一种伪影检测方法和装置 |
| JP2013116132A (ja) * | 2010-11-25 | 2013-06-13 | Toshiba Corp | 超音波診断装置、画像生成方法及び画像処理装置 |
| KR101394274B1 (ko) * | 2013-11-27 | 2014-05-13 | (주) 골프존 | 뎁스 정보 분석을 통한 신체 검출 방법 및 사용자 신체 검출을 위한 뎁스 정보 분석 장치 |
| WO2024034996A1 (ko) * | 2022-08-12 | 2024-02-15 | 아주대학교산학협력단 | 수술 시뮬레이션을 제공하는 방법 및 장치 |
-
1999
- 1999-12-27 JP JP37003499A patent/JP2001184491A/ja active Pending
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005261531A (ja) * | 2004-03-17 | 2005-09-29 | Hitachi Medical Corp | 医用画像表示方法および装置 |
| JP4497965B2 (ja) * | 2004-03-17 | 2010-07-07 | 株式会社日立メディコ | 医用画像表示装置 |
| JP2006075390A (ja) * | 2004-09-10 | 2006-03-23 | Toshiba Corp | 医用画像表示装置 |
| CN102018524A (zh) * | 2009-09-09 | 2011-04-20 | 上海西门子医疗器械有限公司 | 一种伪影检测方法和装置 |
| JP2013116132A (ja) * | 2010-11-25 | 2013-06-13 | Toshiba Corp | 超音波診断装置、画像生成方法及び画像処理装置 |
| KR101394274B1 (ko) * | 2013-11-27 | 2014-05-13 | (주) 골프존 | 뎁스 정보 분석을 통한 신체 검출 방법 및 사용자 신체 검출을 위한 뎁스 정보 분석 장치 |
| WO2015080498A1 (en) * | 2013-11-27 | 2015-06-04 | Golfzon Co.,Ltd. | Method for detecting human body through depth information analysis and apparatus for analyzing depth information for user body detection |
| WO2024034996A1 (ko) * | 2022-08-12 | 2024-02-15 | 아주대학교산학협력단 | 수술 시뮬레이션을 제공하는 방법 및 장치 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4835688A (en) | Three-dimensional image processing apparatus | |
| Stytz et al. | Three-dimensional medical imaging: algorithms and computer systems | |
| JP4209017B2 (ja) | 所望の構造の強調像を作成する方法および装置 | |
| JP6080248B2 (ja) | ３次元画像表示装置および方法並びにプログラム | |
| US4903202A (en) | Three-dimensional object removal via connectivity | |
| JP2004529437A (ja) | オブジェクトオーダーボリュームレンダリングのオクルージョンカリング | |
| US8705821B2 (en) | Method and apparatus for multimodal visualization of volume data sets | |
| US4953087A (en) | Three-dimensional images obtained from tomographic data having unequally spaced slices | |
| JPS6237782A (ja) | ３次元面構造を表示するための装置と方法 | |
| Venu et al. | Comparison of Traditional Method with watershed threshold segmentation Technique | |
| US20070133849A1 (en) | Method and apparatus for visualisation of a tubular structure | |
| JP2008514261A (ja) | 画像処理装置及び方法 | |
| JP4087517B2 (ja) | 領域抽出方法 | |
| CN111047519B (zh) | 用于调整图像观感的方法 | |
| JP4313910B2 (ja) | 画像表示装置 | |
| US8422749B2 (en) | Method and system for myocardial region of interest definition | |
| Ropinski et al. | Internal labels as shape cues for medical illustration. | |
| JP2001184491A (ja) | 三次元画像表示装置 | |
| JP2005020338A (ja) | 異常陰影検出方法および装置並びにプログラム | |
| US20030128890A1 (en) | Method of forming different images of an object to be examined | |
| EP2051205B1 (en) | Method, computer program and workstation for removing undesirable objects for a digital medical image | |
| US20100265252A1 (en) | Rendering using multiple intensity redistribution functions | |
| JP4571378B2 (ja) | 画像処理方法および装置並びにプログラム | |
| JP2002236910A (ja) | ３次元画像作成方法 | |
| Sadowsky et al. | Atlas-assisted tomography: registration of a deformable atlas to compensate for limited-angle cone-beam trajectory |