Nothing Special   »   [go: up one dir, main page]

JP4065327B2 - 投影像表示方法および装置 - Google Patents

投影像表示方法および装置 Download PDF

Info

Publication number
JP4065327B2
JP4065327B2 JP51736898A JP51736898A JP4065327B2 JP 4065327 B2 JP4065327 B2 JP 4065327B2 JP 51736898 A JP51736898 A JP 51736898A JP 51736898 A JP51736898 A JP 51736898A JP 4065327 B2 JP4065327 B2 JP 4065327B2
Authority
JP
Japan
Prior art keywords
image
projection
display
value
setting
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
JP51736898A
Other languages
English (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
Application granted granted Critical
Publication of JP4065327B2 publication Critical patent/JP4065327B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/08Volume rendering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/001Texturing; Colouring; Generation of texture or colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

技術分野
本発明は、断層像が得られる医用画像装置、例えばX線CT装置、MRI装置あるいは超音波断層装置等からの複数の断層像により投影像を得る投影像構成表示方法及び装置に関するものである。
背景技術
従来のこの種の方法としては、下記文献(1),(2)に示すような最大値投影表示法及び積分値投影表示法等がある。
(1)P.J.Keller,B.P.Drayer,E.K.Fram,K.D.Williams,C.L.Dumoulin,S.P.Souza“MR Angiography with Two-dimensional Acquisition and Three-dimensional Display”;Radiology,November 1989 527−532頁
(2)永井純、安達秀雄編著「三次元CT血管造影法」1994年5月30日(株)メディカル・サイエンス・インターナショナル発行 38頁
しかし、上記最大値投影表示法では濃度の最大値の投影像しか表示されず、臓器の分離表示に劣るという問題点があった。また上記積分値投影表示法では、図21に示すように、投影線Lの通過する距離が小さい臓器d2の陰影は、投影線Lの通過する距離が大きい臓器d1の陰影に比較して薄くなり、小さい臓器d2の陰影が大きな臓器d1の陰影に埋もれてしまう等の問題が生じ、臓器d1とd2を識別して表示することが困難であるという問題点があった。なお図21において、eは視点、Lは視線(投影線)、Pは投影面、11〜14は断層像である。上記視点e、視線L、投影面Pの関係式は特開平8−16813号公報に記述されている。
本発明の目的は、臓器等、断層像中の観察対象部位の分離表示に優れ、各部位が鮮明に分離された投影像が得られる投影像構成表示方法及び装置を提供することにある。
発明の開示
本発明は、MR等の三次元計測データから再構成されたボリューム画像や複数の断層像を用い、これらの断層像に対して任意の位置に視点を設定し、この視点位置から前記各断層像を投影面に投影して投影像を構成する投影像構成表示方法において、特定の投影線が、前記複数の断層像を通過する際の各断層像上の通過点での各画素値に対応付けられた画素値計数メモリの、その通過点が投影される前記投影面の画素位置と同一画素位置に、その通過毎に一定値、例えば1を加算して行く処理を、前記投影面の全画素位置を網羅する全投影線について行う画素値別振分け計数処理ステップと、この画素値別振分け計数処理ステップによる記録計数値について画素値対応で所望の重み付けを行う重み付け処理ステップと、この重み付け処理ステップの処理結果データをモニタ画面に投影像として表示する画像表示ステップとを備えている。
また本発明は、同上投影像構成表示方法において、特定の投影線が、前記複数の断層像を通過する際の各断層像上の通過点での各画素値に対応付けられた画素値標識ビットメモリの、その通過点が投影される投影面の画素位置と同一画素位置に、その通過があったことを標識して行く処理を、前記投影面の全画素位置を網羅する全投影線について行う標識記録処理ステップと、各画素値標識ビットメモリの同一画素位置における標識有りの数の合計を全画素位置について求め、その計数値について画素値対応で所望の重み付けを行う重み付け処理ステップと、この重み付け処理ステップの処理結果データをモニタ画面に投影像として表示する画像表示ステップとを備えている。
メモリは各画素値に対応付けられており、それらメモリの同一画素位置における値の和をとり、それに所望の重み付けを行い、その処理結果データをモニタ画面に投影像として表示するので、臓器等、断層像中の観察対象部位の分離表示に優れ、各部位が鮮明に分離された投影像が得られる。
【図面の簡単な説明】
図1は、本発明の一実施例における画像処理の原理を示す説明図である。
図2A、図2Bおよび図3は、図1に示す実施例の実行手順を示すフローチャートである。
図4は、本発明の他の実施例における画像処理の原理を示す説明図である。
図5は、図4に示す実施例の実行手順を示すフローチャートである。
図6は、本発明の他の実施形態を示す説明図である。
図7Aおよび図7Bは、本発明に係る方法を実行することにより得られた投影像のモニタ表示の一例を示す図である。
図8Aから図8Dは、それぞれ本発明に係る重み係数曲線の例を示す図である。
図9Aは、本発明方法で得られた投影像のモニタ表示の他の例を示す図である。
図9B、図9Cおよび図9Dは、それぞれ図9Aに示す3つの表示を実現するための重み係数曲線を示した図である。
図10は、本発明を応用して動画表示を行うための動作原理を説明するための図である。
図11は、本発明方法で得られた投影像のカラー表示の原理を示す説明図である。
図12は、本発明方法で得られた投影像のカラー表示の他の例を示す説明図である。
図13は、本発明方法で得られた投影像と従来の投影像とを並べてモニタ表示した例を示す図である。
図14Aは、図13に示す像にカラーを入れたモニタ表示の例を示す図である。
図14Bは、図14Aに示す図を得るための動作を示すフローチャートである。
図15は、本発明方法で得られた2種類の投影像のモニタ表示の例を示す図である。
図16Aは、図15に示す投影像の1部にカラーを加えたモニタ表示の例を示す図である。
図16Bは、図16Aの投影像を実現するための手順を示すフローチャートである。
図17は、本発明方法で得られた投影像のモニタ表示の他の例を示す図である。
図18は、図17で示した投影像のモニタ表示の手順の一例を示すフローチャートである。
図19は、本発明方法で得られた透視図の観察部を明示したモニタ表示の例を示す図である。
図20は、本発明に係る投影像構成表示装置の一実施例を示すブロック図である。
図21は、従来方法の説明図である。
発明の実施の形態
以下、図面を参照して本発明の実施形態を説明する。
図1は、本発明による投影像構成表示方法の一実施形態の説明図である。この図1において、eは視点、Lは視線投影線、Pは投影面である。
1はX線CT装置による複数のCT像、MRI装置による複数のMRI像、超音波断層装置による複数の断層像、または所要の厚みを持った立体領域における立体的像を同時に再構成して得たボリューム画像である。以下では、複数のCT像を例にとって説明する。これらCT像11−1nは、視点eと投影面Pとの間において積層状に積み重ねられて仮想立体像空間を形成している。
ここで、視点位置、投影線方向、投影面位置及び複数のCT像位置は適宜設定される。また、CT像は原画像だけでなく、補間により求めたCT画像も含む。
11A、11B及び11Cは、各々CT像11〜1n中の同一画素値領域、すなわち同一CT値領域を示すもので、これら同一CT値領域は、単一のCT像のみに二次元的に存在する場合、あるいは隣接する複数のCT像にわたって三次元的に存在する場合がある。図示例における同一CT値領域11A,11B及び11Cは、CT値5、10及び15の3領域を例示している。それら各領域11A,11B及び11C中の同一投影線L上における目盛り状の短い横線は、各CT像の位置を表わすと共に、その投影線Lの各CT像の通過点を表わす。同一CT値領域中の同一投影線L上における目盛り状の短い横線が複数あることは、同一CT値領域が、隣接する複数のCT像にわたって三次元的に存在していることを意味する。
MAはCT値別に設けられたCT値計数メモリで、MA1は第1CT値を持つ画素の計数用、MA2は第2CT値を持つ画素の計数用、…MAnは第nCT値を持つ画素の計数用である。いま、第1CT値をCT値1とし、第nCT値をCT値1000とすると、例えば第1CT値計数メモリMA1はCT値1計数メモリ、第2CT値計数メモリMA2はCT値2計数メモリ、…第nCT値計数メモリMAnはCT値1000計数メモリとなる。CT値計数メモリの枚数や各CT値計数メモリに記録されるべきCT値は、表示しようとする対象のCT値範囲や表示しようとする対象の数などに応じて適宜設定する。
一般に、骨、皮、臓器等の人体の各組織は、それぞれほぼ均一のCT値で、あるいは一定のCT値範囲内でCT像を形成する。従って、同一のまたは一定の範囲のCT値毎に別々にCT値計数メモリに記録すれば、人体の各組織毎に分離された投影像を形成することができる。
これらのCT値係数メモリMA1〜MAnは、図示するように物理的に分かれていなければならないものではなく、1つのメモリ中の記憶領域を適宜割り当てて設定してもよい。
MBは後述のようにCT値計数メモリMA1〜MAnに基づいて各画素毎に算出された、表示されるべき画像データを記憶する積和メモリ(表示メモリ)である。このメモリMBも、CT値係数メモリMAから物理的に分かれていなければならないものではない。
図2Aおよび図2Bは図1に示す実施例における本発明方法の手順の一例を示すフローチャートであって、図2Aは重み付け前の各CT像1における各画素位置を、そのCT値別に振分け計数する処理(振分け計数処理)を、図2BはCT値係数メモリMAから読み出したデータに対する重み付け処理を示す。
図1,図2Aから分かるように、この実施例においては、予め各CT値毎にそのCT値計数用のメモリMA1〜MAnを用意しておく。そして特定の投影線Lが、CT像11〜1nを通過する際の各CT像1上の通過点でのCT値に該当するCT値計数メモリMAの、その通過点が投影される投影面Pの画素位置と同一画素位置に、その通過毎に一定値、ここでは1を加算して行く処理を行う。投影線は、投影面Pの全画素位置を網羅するように順次に移動され、各移動毎に、即ち各画素毎に上記加算処理が実行される。この実行手順を図2Aを参照して更に詳細に説明する。
すなわち図2Aにおいて、ステップ41では全てのCT値計数メモリMA1〜MAnをクリアし、ステップ42では最初の投影線L1を設定する。通常は、後述モニタ画面の左上隅部に相当する投影面Pの位置に向かう投影線L1を設定する。
ステップ43では上記投影線L1が通過するCT像11〜1nの各通過点での画素値、即ちCT値を読み出す。読み出されたCT値が所要のCT値範囲の上限と下限を表わすしきい値の条件に適合するもののみ、それらCT値に対応する計数メモリの、その通過点(画素位置)が投影される投影面Pの画素位置と同じ画素位置に、その通過点に1を加える。すなわち、投影線L1が通過する画素の値(画素値)別の通過回数をメモリMA上に計数して行くもので、全てのCT像11〜1n通過後、ステップ44へ移る。なお、各CT値計数メモリ毎に所定の範囲内の読み出しCT値を計数するようにすることができる。また、ステップ43の動作は特開平8−16813号公報にて知られている。この文献は本願の開示にとり込まれる。
ステップ44では全ての投影線L1〜Lmを設定し終えたか否かを判定する。YESならば処理を終了し、NOならばステップ45へ跳ぶ。
ステップ45では次の投影線L2を設定し、ステップ43へ戻る。以後ステップ44においてYESと判定されるまでステップ43〜45を繰り返す。なお、ステップ44においてYESと判定されるのは、最終の投影線Lm、通常は、モニタ画面の右下隅部に相当する投影面Pの位置に向かう投影線が設定され、その投影線Lmについてのステップ43の処理が終了した時点である。
以上で各CT像1における各画素位置についてのCT値別の振分け計数処理を終え、全てのCT値計数メモリMA1〜MAnの全ての画素位置についての計数値記録が得られる。
なお、上述説明は投影線Lの単位ステップ長、即ち投影線Lが投影面Pに向かって進み、画素値を読み込んで行く場合の画素値読込み間隔がCT像間を進む距離と同一の場合について述べたが、それらは必ずしも全て一致するものではない。一致しない場合があるときは、一致するCT像1の画素値についてのみ上記計数処理を行うことも考えられるが、一致しない場合、すなわちステップ位置に予め積み上げ設定されたCT像がない場合には、その部分の画素値を隣接するCT像の画素値から補間法(本出願人によるJP−A−8−335278等を参照)により求め、その求めた画素値について上述計数処理を行うようにしてもよい。また、上記単位ステップ長は1画素長でもよく、また1画素長より小さくてもよい。
図2Aに示す手順によって得られたCT値別にCT値計数メモリMA1〜MAnに記録された画像データは、それぞれCT値で特定される人体の個々の組織や臓器を表わすものである。従って、観察されるべき所望の臓器を強調しあるいはそれのみを分離して表示することができる。
図2Aに示す計数処理を終えた後、特定の臓器等の所要の観察対象部位を分離表示または強調表示するために、特定のCT値または特定のCT値範囲を記憶したCT値計数メモリMAのデータのみが表示され、または強調して表示されるように、CT値計数メモリ単位で重み付けがされる。図2Bは、説明を簡単にするために最初の投影線L1のみに対応する各CT値計数メモリMA1〜MAnについて積和処理をする手順の一例を示す図である。
図2Bにおいて、まずステップ46で、主メモリ中の一時記憶用の画素データメモリsumをゼロにクリアする。次いで、ステップ47で最初のCT値計数メモリMA1を指定する。
ステップ48ではそのメモリMA1での最初の投影線L1についての、すなわちモニタ画面の左上隅部に相当する投影面Pの対応位置でのメモリ値(記録計数値)が0か否かを判定する。0ならばステップ50へ跳ぶ。0でなければステップ49へ移る。このステップ48は後述ビットマップを用いた発明方法の場合に特に有用である。図1に示す実施例ではステップ48はなくてもよいがあっても処理に支障を来すものではない。
ステップ49では上記メモリ値を重み係数と掛け算してその積を画素データメモリsumに加え、ステップ50へ移る。なお、重み付けについては更に後述する。
ステップ50では全てのCT値計数メモリMA1〜MAnを指定したか否かを判定する。YESならばステップ52に移り、NOならばステップ51へ跳ぶ。
ステップ51では次のCT値計数メモリMA、ここではCT値2計数メモリMA2を指定し、ステップ48へ跳ぶ。以後ステップ50においてYESと判定されるまでステップ51,48〜50を繰り返す。
ステップ52では画素データメモリsumの画素データを表示メモリ(積和メモリ)MBの対応する記憶領域に格納する。
以上は投影線1本当たりの重み付け処理について述べたもので、以後全ての投影線L1〜Lmについて上述重み付け処理を繰り返す。
なお上述実施例では、全ステップを図2A,図2Bに示すような2段階処理で実行することとしたが、図3に示すように連続処理で実行するようにしてもよい。
この図3において、図2A,図2Bと同一又は相当のステップについては同一符号を付してその説明を省略する。図2A,図2Bと比較すれば分かるように、図3では投影線Lの変更設定を最終ステップ44で行うようにしたものである。すなわち、ここでは特定の投影線Lに関し、各CT像1における各画素位置についてのCT値別の振分け計数処理を行い、その記録計数値について重み付け処理をした後に、次の投影線Lを設定し、同様に計数処理、重み付け処理を続けて行い、以下これを最終の投影線Lmに関してまで繰り返すものである。
このように、各処理,各ステップの順序は適宜変更可能であり、最終的に全ての投影線L1〜Lm(投影面Pの全画素)に関して、計数処理,重み付け処理され、表示メモリ(積和メモリ)MBにその結果が格納されればよい。
表示を避けたいCT値の像に対しては重み付け計数を零にすれば、そのCT値の像を消去することができ、強調して観察したいCT値の像に対しては高い重み付け係数を与え、弱めたいCT値の像に対しては低い重み付け係数を与えるようにすればよい。
この様にして表示メモリMBに得られた画像データを表示装置で表示すれば、所要のCT値の像を強調して表示することができ、不要なCT値の像は消去されることができる。また、CT像11〜1nにおいて、積層方向に厚みを持った同一CT値の像はCT値計数メモリでの記録値が大きくなるので、表示装置に表示される像は積層方向に厚みの大きい部分ほど濃く表示され、投影像としての特徴をよく表現することができる。
図4は投影像構成表示方法の第2の実施例の説明図である。この図4において、図1と同一符号は同一又は相当部分を示す。
MC1〜MCnはCT値別に設けられたCT値標識ビットメモリで、MC1は第1CT値標識用ビットメモリ、MC2は第2CT値標識用ビットメモリ、…MCnは第nCT値標識用ビットメモリである。いま、第1CT値をCT値1とし、第nCT値をCT値1000とすると、第1CT値標識ビットメモリMC1はCT値1標識ビットメモリ、第2CT値標識ビットメモリMC2はCT値2標識ビットメモリ、そして第nCT値標識ビットメモリMCnはCT値1000標識ビットメモリとなる。
それらCT値標識ビットメモリMC1〜MCnの番号に対応付けられるCT値は、用いられるCT値標識ビットメモリの枚数や構成表示しようとする対象のCT値範囲等に応じて適宜設定される。これらのメモリMC1〜MCnは、図示するように物理的に分かれていなければならないものではなく、1つのメモリ中の記憶領域を適宜割り当てて設定してもよい。このことは表示メモリMBについても同様である。
図5は、図4に示す本発明方法の手順の一例を示すフローチャートであって、各投影線が通過する各CT像11〜1nにおける各画素位置のCT値を読み出し、各投影線毎にCT値別に投影線の通過の有無を標識としてCT値標識ビットメモリMC1〜MCnに記録する処理をするものである。
即ち、図4と図5から分かるように、予め各CT値毎にそのCT値標識ビットメモリMC1〜MCnを用意しておく。そして特定の投影線Lが、CT像11〜1nを通過する際の各CT像11〜1n上の通過点でのCT値に該当するCT値標識ビットメモリMCの、その通過点が投影される投影面Pの画素位置と同一画素位置に、その通過があったことを標識して記録して行く処理を、表示メモリ(積和メモリ)MBの大きさに対応する投影面Pの全画素位置を網羅する全投影線について行うものである。
図5において、ステップ41′では全てのCT値標識ビットメモリMC1〜MCnをクリアし、ステップ42では最初の投影線を設定する。最初の投影線は、通常は、モニタ画面の左上隅部に相当する投影面Pの位置に向かう投影線L1である。
ステップ43′では上記投影線L1が通過するCT像11〜1nの各通過点でのCT値について、それが観察をする必要があるCT値範囲を規定するしきい値(上限、下限等)の条件を満たす場合に、その画素値に対応するCT値の標識ビットメモリMCの、その通過点(画素位置)が投影される投影面Pの画素位置と同じ画素位置に、その通過があったことの標識“1”を立てる。
従って、投影線L1が通過する画素の値(CT値)別にその通過があったことを標識して行くもので、全てのCT像11〜1n通過後、ステップ44へ移る。なお、画素値とCT値の対応は適宜定められるが、各CT値標識ビットメモリはそれぞれ所定範囲のCT値を検出した時に標識“1”を立てるようにしてもよい。またCT像11〜1nの通過にあたり、同じ画素値を複数回通過しても、最初の画素値の通過で立てられた標識“1”に何ら変化は生じない。
ステップ44では全ての投影線L1〜Lmを設定し終えたか否かを判定する。YESならば処理を終了し、NOならばステップ45へ跳ぶ。
ステップ45では次の投影線L2を設定し、ステップ43′へ戻る。以後ステップ44においてYESと判定されるまでステップ43′,44,45を繰り返す。なお、ステップ44においてYESと判定されるのは、最終の投影線Lm、通常は、モニタ画面の右下隅部に相当する投影面Pの位置に向かう投影線が設定され、その投影線Lmについてのステップ43′の処理が終了した時点である。
以上で各CT像1における各画素位置についてのCT値別の通過の有無の標識処理を終え、全てのCT値標識ビットメモリMC1〜MCnの全ての画素位置(投影面Pの全画素位置に対応する画素位置全て)についての標識“1”の有無の記録が得られる。
上述標識記録処理を終えると、その標識記録値に対して特定のCT値またはCT値の範囲、換言すれば臓器等、CT像1中の観察対象部位を分離表示、更には強調表示をするために、図2Bと同様の重み付けをする。この場合、図2B中の「メモリ」は「ビットメモリ」と読み替えられる。
以上に述べた第2の実施例においては、各投影線上に所定のCT値を持った画素の有無を記録するのみであるから、図2Bに示す処理を全ての投影線について実施して画素メモリに画像処理データを記録して、これをモニタに表示すれば、選択表示されたCT値の像は、視点eから見た被観察対象の平面的形状が表示されることになる。従って、観察対象部位がCT像の積層方向に厚みを持っていても、その厚みの大きさは投影像に全く影響を与えない。
なお、上述の第2の実施例においても、投影線Lの単位ステップ長が一致しない場合のCT像(画素値)の補間や単位ステップ長の長さ等については、図1に示す第1の実施例の例と同様である。
図1、図4の例における各CT像11〜1nの各画素位置の、投影面P上の画素位置への変換(座標変換)は、例えば本出願人による特開平7−210704号公報および特開平7−296184号公報等に記載の方法が適用される。
また、図1、図4の例では、複数のCT像11〜1nの投影面Pへの投影法に中心投影法を用いたが、図6に示すように平行投影法を用いて投影してもよい。この図6は、図4の実施例において視点eを無限遠点へ持っていったのと同様な平行投影法を適用した例を示しているが、図1の実施例に平行投影法を適用してもよい。図6において図4と同一又は相当部分に同一符号を付すことでその説明を省略する。
以上の処理の結果、表示メモリ(積和メモリ)MBには投影像が構成され、これをD/A変換してモニタに与えることにより、投影像71がモニタ表示される。
図7A,図7Bは、モニタ画面24aに投影像その他の関連情報等の表示例を示す図である。
図7Aにおいて、投影像71aは視線方向(透視角度)0°の表示例であり、投影像71bは同像を視線方向90°から見た場合の表示例である。図7Aには、同一の被観察対象を異なる視線方向から見た2つの投影像71a,71bが並べて表示されている。
図1から図3に示す第1の実施例に基づきCT値毎に計数処理がされておれば、どのCT値の部位または臓器を分離し、または強調して表示するかは重み係数をCT値毎に別々に設定すればよい。この重み係数の設定方法は、図7Aに示すように、横軸にCT値、縦軸に重み係数を取った重み係数曲線72を表示し、この重み係数曲線72の形状を修正することにより、各CT値に対して任意に重み付けの設定をすることができる。具体的には、予め設定された3種類の重み付け係数曲線がメモリに登録されており、表示された3つの重み係数曲線選択ボタン73a,73b,73cの内の1つをカーソルでポイントして選択すると、それに対応する重み係数曲線72が図7Aの曲線72のように表示される。
重み係数曲線72に修正を加える場合は、その曲線の各辺に表示されている小さな四角のマーカーをカーソルで指示して、その指示辺を上下または左右に移動させることにより曲線の一部分を変化させることができる。
なお、図7Aでは特性曲線を矩形波状に示したが、連続的な波形の曲線を用いることも可能で、その曲線上に多数のマーカーを表示して、そのマーカーをカーソルで上下移動させて曲線に修正を加えるようにすることができる。
以上の様なカーソルによる曲線の表示および表示曲線の移動修正はコンピュータの分野で一般によく知られた技術をそのまま応用したものである。また、横軸と縦軸のそれぞれの物理量に対応して特性曲線を表示する技術もコンピュータ分野で一般に行われている技術を応用する。
図7Bは、上記投影像71aと、その元となった複数のCT像11〜1nのうち適宜選択された1つのCT像1aとを並べて表示した例である。ここでは、重み係数選択ボタン73a,73bまたは73cで選択される各重み係数曲線72は山の頂部が1つとなるように設定されている。投影像71aはマウスカーソル26aで指示するCT像1a中の任意箇所の濃度値を山の頂部として重み付けされた時に得られた投影像を表示した例を示す。なお、山の頂部は例えばX線CT像ならばX線吸収が大きいことに対応する。
なお、図7Aや図7Bに示すような複数画像の表示は、それら画像に対応する画像データをそれぞれ画像毎に表示メモリに番地を指定して一時的に記憶して、それをそのままモニタに表示するようにした一般的な画像処理技術を用いて行われる。
図8Aから図8Dは重み係数曲線72の変更例を示す。上述したように、これら例示の曲線72を標準形状として、モニタ画面24a上でマウスカーソル26aを操作することによりその重み係数曲線72の形状を修正可能にしてもよい。また、修正された重み係数曲線72を磁気ディスク等に登録し、再利用可能にしてもよい。
図9Aは同一の被観察部位に対して相互に異なった重み付け曲線を与えた場合の3通りの投影像を並べてモニタ表示した例を示す図である。表示像71cは図9Bに示す重み付け曲線に従って、皮膚と造影剤入りの血管の部分のCT値のみに重み付け係数が与えられた時に得られる投影像である。表示像71dは図9Cに示すように造影剤入り血管の部分のみのCT値に重み付けがされた時の投影像である。また、表示像71eは図9Dに示すように皮膚部のCT値のみに重み付けが与えられた時の投影像である。
図7A、図8Aに示すように、重み係数曲線72のモニタ表示時に、モニタ画面横方向にとられたCT値に対応する位置に肝臓、皮膚、骨等、CT像中の観察対象部位の名前を付加表示すれば、所望の部位を分離、強調表示させる際に、視覚に従って重み係数曲線72の形状を容易に修正することができる。この様な部位の名前の付加表示は従来周知の画像処理技術を用いて行うことができる。
各CT値計数メモリMA,CT値標識ビットメモリMCと対で、もう1組の距離用メモリ(Zバッファ)を追加し、各投影線L上で視点eから最も近い画素位置までの距離を各々記録する。そうすることによって、投影面P上の同一位置に2つ以上の画素が重なって存在する場合に、距離の近い画素を表示して遠い画素を消去する隠面消去(Zバッファ法)の適用に有効である。また、表示される画素の投影線L上の視点eからの距離を各々記録しておけば、その距離が短い画素程、明るく表示する陰影付け(デブス法)も適用可能になる。この様な表示はJP−A−8−16813に記載された技術を用いて実現することができる。この文献は本願の開示に取り入れられる。
図1に示す実施例においては、CT値計数メモリMA1〜MAnのおのおのと各CT値とを1対1に対応付け、すなわち、第1CT値計数メモリMA1はCT値1計数メモリ、第2CT値計数メモリMA2はCT値2計数メモリ…としたが、これを1対複数としてもよい。例えば、第1CT値計数メモリMA1はCT値1〜10計数メモリ、第2CT値計数メモリMA2はCT値11〜20計数メモリ…としてもよく、これによればCT値計数メモリMAの数が少なくてすみ、メモリの節約になる。
また図1に示す実施例においては、各CT値計数メモリMAはそのCT値の領域の厚みデータと見ることもできることから、その厚みデータに基づいて濃度勾配を計算して陰影付けした画像を得ることもできる。この濃度勾配は複数のCT値計数メモリ間の濃度勾配であってもよい。この様にして表示される画像は勾配の大きな部分が濃く表われる。
図1および図4に示す実施例においては、各CT像1の各画素値(CT値)に対してしきい値処理(図2A,図3,図5中のステップ43,43′参照)を行ったが、このしきい値処理は省略してもよい。その場合には、画素値(CT値)の全範囲(全域)が構成表示対象となり、通常のX線投影像に相当する投影像が得られる。
また図1および図4に示す実施例において、重み付けせずにまたは均一の重み付けで、所望の1又は複数のメモリMAまたはMCを選択し、それのみで図2A,図3または図5に示す処理を行った投影像を構成表示するようにしてもよい。
更に図1および図4に示す実施例においては、視点eが1つの場合を説明したが、視点eを所定距離ずらした位置に一対設定し、各々について図2A,図3,図5に示す処理を行えば視差のある2つの投影像が構成表示でき、ステレオ視することも可能になる。例えば、図20に示すように、2つの投影像を30Hzの繰返し周波数で交互に同一表示部に表示し、観察者は液晶シャッタメガネ55を掛けて、そのメガネをCRTモニタ24の垂直同期信号に同期して左右の目を交互に見えるようにする。
所望の1若しくは複数の、又は全てのCT像1の所望箇所に事前に特定の画素値を書き込んでおき、それを図2A,図2B,図3または図5に示すフローに従って処理して投影像を構成、表示するようにしてもよい。例えば、CT像1を画像メモリに読み出してモニタに表示し、マウスの操作により特定の臓器の周辺部に高濃度で表示される画素値を書き込んで、これを元のCT像メモリ部へ記憶させておけば、その部分は実際の画素値を無視して高濃度で表示されることになり、輪郭強調の効果が得られる。また、そのような事前処理されたCT像1による投影像を複数得て、後述のような動画表示をするようにしてもよい。
図4および図6に示す実施例において、投影像を表示する際に、所望のCT値標識ビットメモリMC上に標識“1”が記録された画素位置に対応する位置だけを例えばカラーで、若しくは高輝度に、又は点滅させて、表示するようにしてもよい。これによれば、投影像中の所望のCT値部分が他のCT値部分とは区別されて明示的に表示される。なお、ここで、カラー表示方法については後述する。また高輝度に表示する方法は所望のCT値標識ビットメモリMCからの読み出し画像データに対する重み付けを大きくすればよい。更に、点滅表示をする方法は、所望のCT値標識ビットメモリの出力データを例えば1秒間隔で表示と非表示を繰返すようにする。
図10に示すように、重み係数をNo.1〜No.nと相違させ、各重み係数毎に視点位置および/または視線方向、ここでは視点位置を一定方向に順番に変えて構成した複数の投影像群(画像群FG1〜画像群FGn)を記憶媒体121に格納しておく。一方、CRTモニタ画面24aには上記画像群FG1〜画像群FGnの表示用ボタンSW1〜SWnを表示しておく。そして、画像群FGk(k=1〜n)表示用ボタンの1つSWk(k=1〜n)をマウス(図示せず)でクリックすると、重み係数No.k(k=1〜n)の画像群FGkをなす複数の投影像が順番に画像メモリに読み出されCRTモニタ画面24aに動画像122として表示される。これによれば、重み係数の大きさに応じた部位、例えば骨部、軟部組織あるいは血管等が強調された投影像122を動画で観察できる。
具体例として、いま、原画像である複数のCT像11〜1nが人体の腰部の積み上げ画像であると仮定して、視点位置はその体軸回りを所定の方向に所定角度間隔で変化させ、視線方向は全て複数CT画像の軸中心(原点)方向に共通として骨部及び軟部組織の重みを大きくした重み係数でそれぞれ画像群FG1及びFG2を作成したとする。この場合には、画像群FG1表示用ボタンをマウスでクリックすると、腰部の骨部分の、体軸を中心軸とした回転像がCRTモニタ画面24aに表示され、また、画像群FG2表示用ボタンがマウスでクリックされれば、腰部の肉や皮膚等の軟部組織部分の、体軸を中心軸とした回転像がCRTモニタ画面24aに表示されるもので、腰部の骨部、軟部組織部分の比較診断に有効である。なお、複数の回転像等の動画、上述例では腰部の骨部分及び軟部組織部分の回転像がCRTモニタ画面24a上に並べて表示されるようにしてもよい。また、上記回転像の中心軸を体軸から移動可能(偏心可能)として視線が原点以外を向くようにしてもよく、これによれば、診断したい部分の近傍に上記中心軸を置くことでその部分の回転に伴う振れを少なくすることができ、観察がしやすくなる。
図7Aに示す2つの投影像71aと71b、図7Bに示す投影像71aと断層像26a、あるいは図9に示す3つの投影像71c〜71e等のように、視線(投影線)方向や視点位置を異にして得られた複数の投影像、投影像と断層像あるいは重み係数曲線72の形状を異にして得られた複数の投影像等、複数の画像を表示する場合、図10に示す動画表示の実施例を利用して少なくとも1つの表示画像を動画表示するようにすることができる。
図11に示すように、例えば、メモリ131には軟部組織についての重み係数のみで重み付けした投影像を、メモリ132には骨部についての重み係数のみで重み付けした投影像を、メモリ133には血管部についての重み係数のみで重み付けした投影像を、各々格納しておく。
次に、メモリ131〜133に格納された投影像を次のようなカラーマップ134に従って加算、合成する。図11に示すように、カラーマップ134は0〜255の256階調で、そのうち0〜t1が白黒(モノクローム)、(t1+1)〜t2が赤色、(t2+1)〜255が青色に設定されている。メモリ131の値は各画素の濃淡に従い上記0〜t1の範囲に規格化された値で、メモリ132の値は同じく上記(t1+1)〜t2の範囲に規格化された値で、メモリ3の値は同じく上記(t2+1)〜255の範囲に規格化された値で、メモリ135に加算される。これによれば、軟部組織が白黒の濃淡、骨部が赤色の濃淡、血管部が青色の濃淡で各々表示された合成投影像136が得られる。なお、各色共、256階調の0に近いほど薄い色となり、255に近いほど濃い色となる。
ここで、投影像表示(白黒表示を含む)する場合、その表示色はカラーマップに依存する。したがって、表示メモリ(積和メモリ)MBの格納値(積和結果)は一定の範囲、例えば0〜255の範囲に分布させる必要がある。これには、上記格納値に所定の定数を掛ければよい。この定数は、処理時間の節約のため、表示メモリ(積和メモリ)MB上の数ラインだけサンプリングしてその平均値を計算し、これを平均値と推定し、
平均値×定数=範囲内の中央値
定数=範囲内の中央値÷平均値
となるようにすることで決めてもよい。
すなわち、階調範囲が0から255とすると、その範囲の中央値は128であり、サンプリングで求めた平均値を200とすれば、上記式中の定数は0.64となる。従って、各画素のCT値に0.64を掛算をして表示階調を求めることができる。計算された各画素の表示階調が0から250の範囲を超える時はその画素は最小値0または最大値250に設定される。
図11の合成投影像136に示すような各部位毎の色分け表示は、上記のように加算し合成せずとも可能である。すなわち、図12に示すように、所望のCT値計数メモリMA上の画像を所望の色に対応させて表示することによっても可能である。例えば、第5CT値計数メモリMA5の値は各画素の濃淡に従いt1〜t2の範囲(黄の濃淡)に、第2CT値計数メモリMA2の値は同じくt3〜t4の範囲(青の濃淡)に、第15CT値計数メモリMA15の値は同じくt5〜255の範囲(赤の濃淡)に各々規格化された後に、合成用メモリ(図示せず)上でOR(論理和)をとり、CRTモニタ画面24aで表示する。
これによれば、所望のCT値計数メモリMA上の画像毎に、換言すれば各部位毎に色分け表示された合成投影像(図示せず)が得られる。なお、ここで用いるカラーマップ141も0〜255の256階調であり、各色共、256階調の0に近いほど薄い色となり、255に近いほど濃い色となる。
図13に示すように、図1または図4に示す実施例において得られた投影像(重み付け投影画像)151と、この投影像151と同様の対象を同方向から得た、可能ならば同じ視点位置、視線方向から得た、従来の最大値投影表示法等による濃度の最大値をとった画像(MIP画像)152とをCRTモニタ画面24a上に並べて表示されるようにしてもよい。
これによれば、従来のMIP画像152による投影線上の濃度最大値の画像情報と、それ以外の所望の部位に重み付けした本発明方法による投影像151による画像情報と、を相補的に用いて画像を観察、診断することができる。
図13において、投影像151上で所望のCT値をもつ画素位置に対応するMIP画像上の画素位置をカラー表示してもよい。MIP画像152上のカラー表示する位置は、例えばMIP画像152と図4に示すCT値標識ビットメモリMCのうち、所望のCT値標識用のビットメモリMCとを対照し、図14Bのフローチャートに示すように、MIP画像152上の画素を順番に走査しながら各画素位置毎に所望のCT値標識用のビットメモリMCを参照し、そのメモリMC上の対応画素位置に標識“1”が立っていれば、MIP画像152のその画素位置をカラー化する。これをMIP画像152の全画素位置について行う。図14Aはその一例を示す図で、MIP画像152中の斜線部161がカラー化された所望のCT値をもった領域である。なお図14Aにおいて、図13と同一符号は同一又は相当部分を示す。
このようにすれば、投影線上の濃度最大値の画像情報のみからなるMIP画像152中の所望のCT値部分が他の部分とは区別されて明示的に表示される。明示的表示が可能ならば、カラー表示に代えて高輝度に、又は点滅させて、表示するようにしてもよい。
図7Bでは、本発明方法による重み付け投影画像71aと、その元となった複数のCT像11〜1nのうち適宜選択されたCT像とをCRTモニタ画面24a上に並べて表示した例を示したが、図15に示すように、投影画像71aと三次元的に表示した画像(以下、三次元画像という。)171とを並べて表示してもよい。この三次元画像は複数のCT像11〜1nを用いて陰影付け三次元画像171として再構成される。この様な三次元画像の構成はJP−A−8−16813に記載されている。この場合、三次元画像171の構成時の視点位置、視線方向を、上記投影像71aの構成時と同様に設定すれば、両像71a,171の比較対照に極めて有効である。
図15において、投影像71a上で所望のCT値をもつ画素位置に対応する三次元画像171上の画素位置をカラー表示してもよい。三次元画像171上のカラー表示する位置の設定は図16Bに示すフローチャートに従って行われる。三次元画像171と図4に示すCT値標識ビットメモリMCのうち、所望のCT値標識用のビットメモリMCとを対照し、三次元画像171上の所定の画素位置に対応する所望のCT値標識用のビットメモリMC上の画素位置に標識“1”が立っていれば三次元画像171のその画素位置をカラー化する。これを三次元画像171の全画素位置について行う。図16Aはこれを実行するためのフローチャートの一実施例を示す図で、三次元画像171中の斜線部181がカラー化された所望のCT値をもった領域である。なお図16Aにおいて、図15と同一符号は同一又は相当部分を示す。
これによれば、三次元画像171中の所望のCT値部分が他の部分とは区別されて明示的に表示される。明示的表示が可能ならば、カラー表示に代えて高輝度に、又は点滅させて、表示するようにしてもよい。
図17に示すようにCRTモニタ画面24aに表示する重み付け投影画像が血管像191(一般には臓器像)であって、その血管内にカテーテル等の血管内移動医用具を挿入する場合には、その血管内移動医用具像192を上記CRTモニタ画面24a上の血管像191に合成して表示するようにしてもよい。ここでは、上記合成表示に加えて、その合成像の回転をも行っている。図17では縦軸を中心に(a),(b)---(c),(d)と順に回転している例を示しているが、横軸を中心として回転するようにしてもよい。また、縦軸中心に回転する画像と横軸中心に回転する画像を並べて表示するようにしてもよく、更に両者を切り替えて表示するようにしてもよい。なお、これらの表示方法は図10を参照して説明した実施例を適用することによって実現することができる。
このような画像表示の手順の一例を図18に示す。
ステップ201では血管内移動医用具像、ここではカテーテル像192の座標値変更のフラグ領域をクリアする。
ステップ202ではカテーテル座標値変更のフラグがオンか否かを判定し、オンならばステップ203に進み、そうでなければステップ204に飛ぶ。
ステップ203ではカテーテル座標値をカテーテル座標検出手段(図示せず)から受け取り、カテーテル座標値変更のフラグをオフにする。
ステップ204では受け取ったカテーテル座標値を上記血管像191の座標系での対応座標値に変換する。
ステップ205では座標値変換されたカテーテル像192を上記血管像191に合成し、表示する。
ステップ206では上記血管像191の視点位置を所定の回転軸を中心として所定の方向に所定角度進ませた位置に変更する。
以下、ステップ202〜206が繰り返され、図17(a)〜(d)に示すように回転する合成像(血管像191とカテーテル像192)がCRTモニタ画面24aに表示される。
なお、ステップ207、208はカテーテル座標検出手段側のフローで、ステップ207では、カテーテルの挿入座標が変化したか否かを判定する。変化した場合にはステップ208に進み、そうでなければステップ207を繰り返す。
ステップ208では、上記カテーテル像192の座標値変更フラグをオンにし、その座標値をステップ203で取り込む。
このような血管像191とカテーテル像192の合成表示によれば、カテーテル等の血管内移動医用具の挿入に役立つ。
図19は本発明の一応用例を示す図であって、画像210には血管内腔面211が表示されている。しかし、この画面210は全体像212の中のどの部分を診ているのか一般には明確ではない。そこで画像212中に、画像210を得るために、本発明による実施に係る視点213、視野および投影面215を破線で示すような三角形または四角錘(図示せず)で表示する。この視点、視野および投影面の表示により、画像210が画像212中のどの部分をどんな方向に見た時のどの面上における投影像であるかが明確に認識することができる。
なお、視点、視野および投影面を画面上に表示する方法、およびそのような方法により指定された部分の投影像を再構成して表示する方法については前記のJP−A−8−16813に開示されている。
以上の説明においては、X線CT像について述べたが、本発明においては、これに限られるものではなく、MRI断層像、三次元的なMRIボリューム像または超音波断層像にも適用することができる。なお、MRI像や超音波像ではCT値に対応する画素値として濃度値が用いられる。
図20は、発明による投影像構成表示装置の概略構成の一実施例を示すブロック図である。図20において、20は磁気ディスク、21は主メモリ、22はCPU(中央処理装置)、23は表示メモリ、24はCRTモニタ、25はキーボード、26はマウス、27はマウスコントローラ、28は共通バスである。
この場合、磁気ディスク20は、構成対象であるCT像等の複数の断層像またはボリューム画像、上述計数処理又は標識記録処理及び重み付け処理並びに画素位置座標変換等の実行プログラム、オペレーションプログラム、重み係数曲線等を格納するものである。またこの磁気ディスク20は、各種プログラムの実行により得られる各種データ、例えば構成後の投影像等も格納する。主メモリ21はプログラム実行時に種々のプログラム、データを一時的に格納するものである。なお、磁気ディスクに代えて光ディスクを用いてもよい。CPU22はプログラム実行時に必要な演算、処理及び装置各部の制御を行うものである。
表示メモリ23はCT像、投影像、重み係数曲線等のデータ、マウスカーソルの表示用データ等を格納する。CRTモニタ24は表示メモリ23に格納されたデータを画像として表示するものである。すなわちCRTモニタ24には、各種プログラムの実行に伴って、CT像や投影像、重み係数曲線あるいはマウスカーソル等が表示される。
キーボード25は各種の操作指令やデータ等を入力するものである。マウス26は例えば視点、投影面、断層像等の位置入力、視線(投影線)方向、重み係数曲線の選択や修正等を行うのに用いられるものであり、マウスコントローラ27はそのマウス26を制御するものである。共通バス28は以上の装置各部を接続するものである。なお、液晶シャッタメガネ55は、前述のように、CRTモニタ24の垂直同期信号に同期してメガネの左右液晶を交互に遮光して、交互に表示される視点の異なる2つの同一対象像を立体像として見るようにするものである。
以上の構成において、CPU22はキーボード25やマウス26の操作により図2A,図3および図5に示す処理フローを演算し実行し、投影像を構成し表示する。
産業上の利用分野
以上に述べたように、本発明に係る複数の断層像から投影像を構成し表示する方法および装置は、X線、MRIおよび超音波画像処理装置および方法への適用に有用である。

Claims (30)

  1. 複数の断層像を積層したまたは少なくとも1つのボリューム画像からなる被観察体の仮想立体像を形成する立体像形成工程と、
    前記形成された仮想立体像に対して所要の位置に視点を設定し、その設定された視点から前記仮想立体像を投影する投影面を設定する設定工程と、
    前記視点から前記投影面に向かう投影線によって横切られる前記仮想立体像の画素を該投影線上の所要間隔毎に画素値データとして読み出す読出工程と、
    該読み出された画素値を画素値毎に複数の画素値メモリに記録する記録工程と、
    前記画素値メモリに記録された画素値データから投影像を構成する投影像構成工程と、
    該構成された投影像を表示装置に表示する表示工程と、
    を含むことを特徴とする投影像表示方法。
  2. 請求項1に記載の投影像表示方法において、前記複数の画素値メモリから読み出された画素データ毎に所要の重み係数を設定するための重み係数情報を前記表示工程に表示させる表示制御工程を更に含み、前記設定工程は前記表示制御工程によって表示制御された重み係数情報を参照して重み係数を設定することを含む、方法。
  3. 請求項2に記載の投影像表示方法において、前記投影像構成工程は異なる視線方向からの複数の投影像を構成し、前記表示工程は異なる視線方向からの投影像をそれぞれ表示し、前記設定工程は前記重み係数をそれぞれ表示された投影像を参照して設定する、ことを含む方法。
  4. 請求項2に記載の投影像表示方法において、前記表示工程は前記投影像構成工程によって構成された投影像と断層像をそれぞれ表示し、前記設定工程は前記重み係数をそれぞれ表示された投影像と前記断層像を参照して設定する、ことを含む方法。
  5. 請求項2から4の何れか一項に記載の投影像表示方法において、前記重み係数の設定に際し、前記表示制御工程は画素値の意味する部位情報を付加表示する、ことを含む方法。
  6. 請求項1に記載の投影像表示方法において、前記断層像またはボリューム画像を構成する画素値(CT値)毎に投影像を前記投影像構成工程によって構成し、これら構成された画素値毎の投影像を前記表示工程によって表示し、前記設定工程は前記表示された画素値毎の投影像を参照して重み係数を設定する、ことを含む方法。
  7. 請求項6に記載の投影像表示方法において、前記重み係数の設定に際し、複数種類の重み係数から選択する重み係数選択工程、を更に含む方法。
  8. 請求項2に記載の投影像表示方法において、前記重み係数の設定に際し、複数種類の重み係数を有し、それらの重み係数の番号と前記表示装置に表示されるソフトスイッチとを対応づけて重み係数を記憶媒体に記憶する重み係数記憶工程を更に含み、前記表示制御工程は前記記憶媒体から読み出された重み係数を前記ソフトスイッチの選択により順次変化させて前記投影像を動画的に表示させることを含む、方法。
  9. 請求項1に記載の投影像表示方法において、前記投影像構成工程は前記投影像に対しその濃度値により白黒または色を対応づけて投影像を構成することを含む、方法。
  10. 請求項2に記載の投影像表示方法において、前記複数の断層像または少なくとも1つのボリューム画像から最大値投影法によりMIP画像を構成するMIP画像構成工程を更に含み、前記表示工程は前記MIP画像と重み付けされた投影像とを同一画面に表示することを含む、方法。
  11. 請求項10に記載の投影像表示方法において、前記MIP画像構成工程はそのMIP画像を構成するための断層像の画素値により特定される特定部をカラー化することを含む、方法。
  12. 請求項2に記載の投影像表示方法において、前記複数の断層像または少なくとも1つのボリューム画像から陰影付け3D画像を構成する陰影付け3D画像構成工程を更に含み、前記表示工程は前記陰影付け3D画像と重み付けされた投影像とを同一画面に表示することを含む、方法。
  13. 請求項12に記載の投影像表示方法において、前記陰影付け3D画像構成工程はその陰影付け3D画像を構成するための断層像の画素値により特定される特定部をカラー化することを含む、方法。
  14. 請求項1に記載の投影像表示方法において、前記投影像構成工程は前記投影像と血管内移動医用具像を合成し、前記表示工程は前記投影像と血管内移動医用具像とを合成表示する、ことを含む方法。
  15. 請求項1に記載の投影像表示方法において、前記複数の断層像または少なくとも1つのボリューム画像から中心投影法によって血管内腔面像を構成する工程を更に含み、前記表示工程は前記血管内腔面像と前記投影像とを同一画面に表示することを含む、方法。
  16. 複数の断層像を積層したまたは少なくとも1つのボリューム画像からなる被観察体の仮想立体像を形成する立体像形成手段と、
    前記形成された仮想立体像に対して所要の位置に視点を設定し、その設定された視点から前記仮想立体像を投影する投影面を設定する設定手段と、
    前記視点から前記投影面に向かう投影線によって横切られる前記仮想立体像の画素を該投影線上の所要間隔毎に画素値データとして読み出す読出手段と、
    該読み出された画素値を画素値毎に複数の画素値メモリに記録する記録手段と、
    前記画素値メモリに記録された画素値データから投影像を構成する投影像構成手段と、
    該構成された投影像を表示装置に表示する表示手段と、
    を含むことを特徴とする投影像表示装置。
  17. 請求項16に記載の投影像表示装置において、前記複数の画素値メモリから読み出された画素データ毎に所要の重み係数を設定するための重み係数情報を前記表示手段に表示させる表示制御手段を更に含み、前記設定手段は前記表示制御手段によって表示制御された重み係数情報を参照して重み係数を設定する手段を含む、装置。
  18. 請求項17に記載の投影像表示装置において、前記投影像構成手段は異なる視線方向からの複数の投影像を構成し、前記表示手段は異なる視線方向からの投影像をそれぞれ表示し、前記設定手段は前記重み係数をそれぞれ表示された投影像を参照して設定する、ようにした装置。
  19. 請求項17に記載の投影像表示装置において、前記表示手段は前記投影像構成手段によって構成された投影像と断層像をそれぞれ表示し、前記設定手段は前記重み係数をそれぞれ表示された投影像と前記断層像を参照して設定する、ようにした装置。
  20. 請求項17から19の何れか一項に記載の投影像表示装置において、前記重み係数の設定に際し、前記表示制御手段は画素値の意味する部位情報を付加表示する、ようにした装置。
  21. 請求項16に記載の投影像表示装置において、前記断層像またはボリューム画像を構成する画素値(CT値)毎に投影像を前記投影像構成手段によって構成し、これら構成された画素値毎の投影像を前記表示手段によって表示し、前記設定手段は前記表示された画素値毎の投影像を参照して重み係数を設定する、ようにした装置。
  22. 請求項21に記載の投影像表示装置において、前記重み係数の設定に際し、複数種類の重み係数から選択する重み係数選択手段、を更に備えた装置。
  23. 請求項17に記載の投影像表示装置において、前記重み係数の設定に際し、複数種類の重み係数を有し、それらの重み係数の番号と前記表示装置に表示されるソフトスイッチとを対応づけて重み係数を記憶媒体に記憶する重み係数記憶手段を更に備え、前記表示制御手段は前記記憶媒体から読み出された重み係数を前記ソフトスイッチの選択により順次変化させて前記投影像を動画的に表示させる手段を含む、装置。
  24. 請求項16に記載の投影像表示装置において、前記投影像構成手段は前記投影像に対しその濃度値により白黒または色を対応づけて投影像を構成する手段を含む、装置。
  25. 請求項17に記載の投影像表示装置において、前記複数の断層像または少なくとも1つのボリューム画像から最大値投影法によりMIP画像を構成するMIP画像構成手段を更に含み、前記表示手段は前記MIP画像と重み付けされた投影像とを同一画面に表示する手段を含む、装置。
  26. 請求項25に記載の投影像表示装置において、前記MIP画像構成手段はそのMIP画像を構成するための断層像の画素値により特定される特定部をカラー化する手段を含む、装置。
  27. 請求項17に記載の投影像表示装置において、前記複数の断層像または少なくとも1つのボリューム画像から陰影付け3D画像を構成する陰影付け3D画像構成手段を更に含み、前記表示手段は前記陰影付け3D画像と重み付けされた投影像とを同一画面に表示する手段を含む、装置。
  28. 請求項27に記載の投影像表示装置において、前記陰影付け3D画像構成手段はその陰影付け3D画像を構成するための断層像の画素値により特定される特定部をカラー化する手段を含む、装置。
  29. 請求項16に記載の投影像表示装置において、前記投影像構成手段は前記投影像と血管内移動医用具像を合成し、前記表示手段は前記投影像と血管内移動医用具像とを合成表示する手段を更に含む装置。
  30. 請求項16に記載の投影像表示装置において、前記複数の断層像または少なくとも1つのボリューム画像から中心投影法によって血管内腔面像を構成する手段を更に含み、前記表示手段は前記血管内腔面像と前記投影像とを同一画面に表示する手段を含む、装置。
JP51736898A 1996-10-08 1997-09-09 投影像表示方法および装置 Expired - Fee Related JP4065327B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP26729296 1996-10-08
JP15668997 1997-06-13
PCT/JP1997/003174 WO1998015226A1 (fr) 1996-10-08 1997-09-09 Procede et appareil de formation et d'affichage d'une image a partir d'une pluralite d'images partielles

Publications (1)

Publication Number Publication Date
JP4065327B2 true JP4065327B2 (ja) 2008-03-26

Family

ID=26484374

Family Applications (1)

Application Number Title Priority Date Filing Date
JP51736898A Expired - Fee Related JP4065327B2 (ja) 1996-10-08 1997-09-09 投影像表示方法および装置

Country Status (4)

Country Link
US (1) US6990231B2 (ja)
EP (1) EP0968683B1 (ja)
JP (1) JP4065327B2 (ja)
WO (1) WO1998015226A1 (ja)

Families Citing this family (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4473358B2 (ja) * 1999-01-21 2010-06-02 株式会社東芝 診断装置
JP4313910B2 (ja) * 1999-10-21 2009-08-12 株式会社日立メディコ 画像表示装置
JP4597331B2 (ja) * 2000-04-05 2010-12-15 株式会社日立メディコ 医用画像表示装置
JP3712348B2 (ja) 2000-07-18 2005-11-02 独立行政法人科学技術振興機構 心房粗動および心房細動の心臓磁界診断装置およびその作動方法
JP3712350B2 (ja) 2000-07-18 2005-11-02 独立行政法人科学技術振興機構 心室遅延電位の心臓磁界診断装置およびその作動方法
JP3712349B2 (ja) 2000-07-18 2005-11-02 独立行政法人科学技術振興機構 生存心筋診断装置およびその作動方法
JP4795527B2 (ja) * 2000-11-24 2011-10-19 株式会社東芝 X線画像診断システム
US20050049494A1 (en) * 2003-08-29 2005-03-03 Arthur Gritzky Method and apparatus for presenting multiple enhanced images
DE102004011155A1 (de) * 2004-03-08 2005-08-18 Siemens Ag Verfahren zur Visualisierung von mit einem bildgebenden, endoluminalen Instrument aufgezeichneten 2D-Bilddaten eines Hohlkanals
US7035371B2 (en) 2004-03-22 2006-04-25 Siemens Aktiengesellschaft Method and device for medical imaging
JP2006180987A (ja) * 2004-12-27 2006-07-13 Hitachi Medical Corp 医用画像表示装置
US9510732B2 (en) 2005-10-25 2016-12-06 Intuitive Surgical Operations, Inc. Methods and apparatus for efficient purging
US8137333B2 (en) 2005-10-25 2012-03-20 Voyage Medical, Inc. Delivery of biological compounds to ischemic and/or infarcted tissue
US10064540B2 (en) 2005-02-02 2018-09-04 Intuitive Surgical Operations, Inc. Visualization apparatus for transseptal access
US20080015569A1 (en) 2005-02-02 2008-01-17 Voyage Medical, Inc. Methods and apparatus for treatment of atrial fibrillation
US11478152B2 (en) 2005-02-02 2022-10-25 Intuitive Surgical Operations, Inc. Electrophysiology mapping and visualization system
US8934962B2 (en) 2005-02-02 2015-01-13 Intuitive Surgical Operations, Inc. Electrophysiology mapping and visualization system
US7480402B2 (en) * 2005-04-20 2009-01-20 Visionsense Ltd. System and method for producing an augmented image of an organ of a patient
KR20080038280A (ko) * 2005-06-02 2008-05-06 미티테크 엘엘씨 객체 또는 생물 형태 상에 특성을 디스플레이하기 위한방법 및 장치
US7764818B2 (en) * 2005-06-20 2010-07-27 Siemens Medical Solutions Usa, Inc. Surface parameter adaptive ultrasound image processing
US9055906B2 (en) * 2006-06-14 2015-06-16 Intuitive Surgical Operations, Inc. In-vivo visualization systems
US20080097476A1 (en) 2006-09-01 2008-04-24 Voyage Medical, Inc. Precision control systems for tissue visualization and manipulation assemblies
US10004388B2 (en) 2006-09-01 2018-06-26 Intuitive Surgical Operations, Inc. Coronary sinus cannulation
JP4922734B2 (ja) * 2006-11-13 2012-04-25 株式会社Aze 医用画像生成装置、方法およびプログラム
US8244015B2 (en) * 2006-11-22 2012-08-14 General Electric Company Methods and apparatus for detecting aneurysm in vasculatures
US8160395B2 (en) * 2006-11-22 2012-04-17 General Electric Company Method and apparatus for synchronizing corresponding landmarks among a plurality of images
US9226648B2 (en) 2006-12-21 2016-01-05 Intuitive Surgical Operations, Inc. Off-axis visualization systems
US8657805B2 (en) 2007-05-08 2014-02-25 Intuitive Surgical Operations, Inc. Complex shape steerable tissue visualization and manipulation catheter
US7991450B2 (en) 2007-07-02 2011-08-02 General Electric Company Methods and systems for volume fusion in diagnostic imaging
JP5130002B2 (ja) * 2007-07-30 2013-01-30 富士フイルム株式会社 階調画像表示システム
DE102007046938A1 (de) * 2007-09-28 2009-05-20 Siemens Ag Verfahren zur kombinierten Bilddarstellung eines im Herzbereich eines Patienten eingebrachten Katheters mit elektrophysiologischen Daten des Herzens
WO2009073009A1 (en) * 2007-12-07 2009-06-11 Craig Schwartz Method and apparatus for projecting viewable data onto an imaged object
US8131044B2 (en) * 2008-04-08 2012-03-06 General Electric Company Method and apparatus for determining the effectiveness of an image transformation process
US8938104B2 (en) * 2008-08-29 2015-01-20 Varian Medical Systems International Ag Systems and methods for adaptive filtering
CN102222352B (zh) * 2010-04-16 2014-07-23 株式会社日立医疗器械 图像处理方法和图像处理装置
JP5050094B2 (ja) * 2010-12-21 2012-10-17 株式会社東芝 映像処理装置及び映像処理方法
WO2012118109A1 (ja) * 2011-03-03 2012-09-07 株式会社 日立メディコ 医用画像処理装置及び医用画像処理方法
US8526692B2 (en) * 2011-06-30 2013-09-03 Wisconsin Alumni Research Foundation Reduction of transitivity errors in radiotherapy image registration
JP5835975B2 (ja) * 2011-07-19 2015-12-24 株式会社東芝 画像処理システム、装置、方法及び医用画像診断装置
US9208559B1 (en) * 2014-07-25 2015-12-08 Siemens Aktiengesellschaft Method and apparatus for gastric artery chemical embolization

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57115238A (en) * 1981-01-08 1982-07-17 Shimadzu Corp Image treating apparatus of computer tomogram photograph apparatus
US4577152A (en) * 1983-04-11 1986-03-18 Albert Macovski Selective material projection imaging system using nuclear magnetic resonance
JPH031845A (ja) * 1989-05-31 1991-01-08 Shimadzu Corp X線画像表示装置
US5128864A (en) * 1989-08-09 1992-07-07 W. L. Systems, Inc. Method for computing tomographic scans
JPH04259454A (ja) * 1991-02-15 1992-09-16 Hitachi Medical Corp X線ct装置
GB9122843D0 (en) * 1991-10-28 1991-12-11 Imperial College Method and apparatus for image processing
US5430783A (en) * 1992-08-07 1995-07-04 General Electric Company Reconstruction method for helical scanning computed tomography apparatus with multi-row detector array employing overlapping beams
US5748696A (en) 1993-11-26 1998-05-05 Kabushiki Kaisha Toshiba Radiation computed tomography apparatus
JP3730274B2 (ja) * 1993-11-26 2005-12-21 東芝医用システムエンジニアリング株式会社 手術支援システム
US5694530A (en) 1994-01-18 1997-12-02 Hitachi Medical Corporation Method of constructing three-dimensional image according to central projection method and apparatus for same
JP3411665B2 (ja) 1994-04-27 2003-06-03 株式会社日立メディコ 三次元画像構成装置
JPH07210704A (ja) 1994-01-18 1995-08-11 Hitachi Medical Corp 三次元画像の構成方法
JP3632862B2 (ja) * 1994-06-24 2005-03-23 株式会社日立メディコ 三次元画像表示方法及び装置
US5900878A (en) 1994-01-18 1999-05-04 Hitachi Medical Corporation Method of constructing pseudo-three-dimensional image for obtaining central projection image through determining view point position by using parallel projection image and apparatus for displaying projection image
JPH0822537A (ja) * 1994-07-05 1996-01-23 Hitachi Medical Corp Ct画像からの透視画像作成方法及び装置
CN1164904A (zh) * 1994-09-06 1997-11-12 纽约州州立大学研究基金会 实时立体显象设备与方法
JP3556292B2 (ja) * 1994-10-28 2004-08-18 株式会社東芝 画像処理装置
JP3411745B2 (ja) 1995-04-05 2003-06-03 株式会社日立メディコ 断層像補間法及び装置
US5920660A (en) 1995-04-05 1999-07-06 Hitachi Medical Corporation Tomogram interpolation method for executing interpolation calculation by using pixel values on projection line
US5720291A (en) 1996-03-22 1998-02-24 Advanced Technology Laboratories, Inc. Three dimensional medical ultrasonic diagnostic image of tissue texture and vasculature
US5673300A (en) * 1996-06-11 1997-09-30 Wisconsin Alumni Research Foundation Method of registering a radiation treatment plan to a patient

Also Published As

Publication number Publication date
US20040165766A1 (en) 2004-08-26
US6990231B2 (en) 2006-01-24
EP0968683A1 (en) 2000-01-05
WO1998015226A1 (fr) 1998-04-16
EP0968683B1 (en) 2011-05-25
EP0968683A4 (en) 2008-04-23

Similar Documents

Publication Publication Date Title
JP4065327B2 (ja) 投影像表示方法および装置
JP4421016B2 (ja) 医用画像処理装置
US7127091B2 (en) Method and apparatus for visualizing a limited part of a 3D medical image-point-related data set, through basing a rendered image on an intermediate region between first and second clipping planes, and including spectroscopic viewing of such region
JP5495357B2 (ja) 画像表示方法及び医用画像診断システム
US20110254845A1 (en) Image processing method and image processing apparatus
US9491443B2 (en) Image processing method and image processing apparatus
JP4122463B2 (ja) 医療用可視画像の生成方法
CN109427059A (zh) 解剖结构的平面可视化
JP5289966B2 (ja) インターベンショナル手順中、シルエットレンダリング及び画像を表示するための画像処理システム及び方法
Williams et al. Volumetric curved planar reformation for virtual endoscopy
CN106030657A (zh) 医学4d成像中的运动自适应可视化
JP4786307B2 (ja) 画像処理装置
JPH0838433A (ja) 医用画像診断装置
US11941765B2 (en) Representation apparatus for displaying a graphical representation of an augmented reality
JP2002544604A (ja) 医用画像用半透明オーバーレイ
JPH1176228A (ja) 三次元画像構成装置
JP2000105838A (ja) 画像表示方法及び画像処理装置
EP2082250A1 (en) User interface for displaying mri images
JP4010034B2 (ja) 画像作成装置
JP2006506142A (ja) ボリュームデータセット内に撮像された対象の表示方法
WO1991014397A1 (en) Three-dimensional graphics simulation and actual imaging data composite display
JPH0526231B2 (ja)
Kase et al. Effective three-dimensional representation of internal structures in medical imaging visualization
JP2774512B2 (ja) 画像表示装置
JPH03232077A (ja) 画像処理装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040909

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040909

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080104

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

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20110111

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20120111

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20130111

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20140111

Year of fee payment: 6

LAPS Cancellation because of no payment of annual fees