JPH0816813A - 画像形成方法及び装置、表示方法及び表示装置 - Google Patents
画像形成方法及び装置、表示方法及び表示装置Info
- Publication number
- JPH0816813A JPH0816813A JP6143496A JP14349694A JPH0816813A JP H0816813 A JPH0816813 A JP H0816813A JP 6143496 A JP6143496 A JP 6143496A JP 14349694 A JP14349694 A JP 14349694A JP H0816813 A JPH0816813 A JP H0816813A
- Authority
- JP
- Japan
- Prior art keywords
- viewpoint
- image
- projection
- dimensional
- line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Image Generation (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
原画像に対して、内視鏡的な視点の移り変わりにより、
腸や気管支、血管等の内部の様子を擬似三次元画像とし
て表示したい。 【構成】 点としての視点と、この視点からの視線方向
とを、三次元原画像の奥行き方向に更新しながら設定す
る。更新毎に、視線方向に直交するような投影面を設定
する。更新毎に、視点と投影面との間に介在する単位三
次元画像を、視点から投影面に中心投影して中心投影像
を得る。この中心投影像に陰影処理及び隠面処理を行い
擬似三次元画像を得、その結果を表示画面に表示する。
かくして、視点と視線方向とを三次元原画像の奥行き方
向に更新してゆき、その更新の都度、中心投影がなされ
る。この結果、これらを次々に表示することで、内視鏡
的な視点の移り変わりによる擬似三次元画像が観察でき
る。
Description
的映像化を行う疑似三次元画像の形成方法及び装置、表
示方法及び装置に関する。
像、超音波断層像等の種々の画像より成る。疑似三次元
画像を得るには、X線CT画像では断層面の画像を積み
上げるやり方をとる。MRI画像では三次元的な画像計
測を行うことができるため、それらを配列処理すること
で三次元原画像を得る。
影面との間に存在する三次元原画像を、その投影面に、
視点からみたように投影する投影法が存在する。投影法
には、平行投影法と中心投影法とが知られている。平行
投影法は、視点が面又は線で構成されており、この視点
から平行に投影面に投影を行うものであり、三次元画像
化した臓器などの対象物をその外側からみた画像として
構成する映像化には適している。中心投影法は、数学上
では知られているが、疑似三次元画像化した医用画像の
投影法としては、未だ採用されていない。
映像化にあっては、外形的な形状も大切であるが内視鏡
的な観察を行ったような映像化も必要である。平行投影
法ではその実現は困難であった。そこで、本件出願人
は、先に、中心投影法を採用した三次元画像構成法につ
いて出願を行った(特願平6−3492号)。この先願
は、三次元画像に対して視点を面や線ではなく点とし、
且つこの点としての視点を中心として位置づけ、視点と
投影面との間に存在する三次元原画像を、視点を中心位
置として中心投影して投影面に投影しようとしたもので
ある。投影面が二次元表示面とすれば、二次元表示面に
は、視点から中心投影法で得た三次元画像が二次元化さ
れて表示される。
の通りである。三次元原画像が腸内部を含む画像である
とすると、視点をこの腸内部に持ってくる。そして視点
からその腸内部に向けて視線を移せば、視点から奥側に
ある腸内部は拡がりを持って観察できるはずである。こ
うした様相を実現するために、視点を中心位置として、
視点よりも奥側を、中心投影法により、投影面である表
示面に投影すれば、この表示面の投影画像は丁度内視鏡
で腸内部を観察したかの如き映像となる。これが先願の
主たる特徴である。
結果をそのまま表示したのでは、画素の遠近が考慮され
ず、且つ三次元原画像上での視点からみての手前側と奥
側との区別もつかないものとなり、画像のリアルティが
なくなる。そこで、先願では、画素の遠近を考慮すべく
陰影づけアルゴリズム法を開示した。例えば、デプス法
の例では、視点と各画素との距離Rを考慮し、投影後の
画素濃度(階調度)Iを下記の如くとる。Cは定数、I
maxは最大階調度(設定した定数)である。
るくなる。
バッファ法を採用する。Zバッファ法とは、隠れた面を
消去する方法であり、過去には平行投影法で採用されて
いた。平行投影法において、各画素から投影面に垂線を
下し、その交点と画素までの距離Z(これは、垂線方向
がZ軸と定義していたためであり、このZ軸成分が距離
Zとなる。距離Zは奥行き距離とも呼ばれている)を求
める。そして、投影面上の同一位置に2つ以上の画素が
重なって存在する場合には、各画素のZをみて、最小の
Zの値をZバッファメモリに記憶し、その他の画素は陰
れ画素として消去する。かくして、最小のZ、即ち一番
手前側の画素のみが選択され、これを表示することで、
手前側の画素のみが観察できる。この平行投影法のZバ
ッファ法を、中心投影法としても採用できる。即ち、Z
バッファの中味のデータをZとするのではなく、点とし
ての視点と投影面の各画素位置との距離Rとする。これ
により、投影面の画素位置に、2つ以上の画素が重なっ
て投影されている場合には、距離Rの最小の画素を選
び、残りの画素は消去する。この結果を表示すれば、隠
面処理された二次元画像が表示面上に得られる。勿論、
デプス法やZバッファは一例にすぎない。
る。簡単のため、2つの平行な面PL1、PL2を考え
る。手前の面PL1が三次元原画像を構成する1つの面
であり、いわゆる被投影面である。奥側の面PL2が投
影面である。手前の面PL1の更に手前に点としての視
点Eを与える。視点Eから投影面PL2にむけて中心投
影を行うと、視点Eから投影面PL2にむけて放射状に
のびる数多くの中心投影線P(P1、P2、……)が得ら
れる。中心投影線Pと被投影面との交点Q(Q1、Q2、
……)が被投影画素位置となる。そして交点Qをつきぬ
けて投影面PL2と交わる。その交点N(N1、N2、…
…)が、被投影画素位置Qの投影画素位置となる。被投
影画素位置Qと投影画素位Nとは、視点Eと被投影面P
L1と投影面PL2とが幾何学的に定まれば座標変換を示
す関数で規定できる。
一般的には、図3の如く面PL1とPL2とは平行とは限
らず、視点Eの位置も種々である。しかし、いずれの態
様であれ、EとPL1とPL2とは幾何学的に規定できる
ものであり、座標変換関数も、一義的に数式で規定でき
る故、何ら問題はない。そして、投影後の画素位置に関
してはコンピュータを利用して座標変換の処理を行い、
画素濃度については、デプス法やZバッファ法を利用す
ることでリアルティを持たせることで、二次元表示画面
上には、疑似三次元画像を内視鏡的にみた映像がリアル
ティを持って表示できることになった。
もあれば、多階調の画像の例もある。2値化画像とは、
関心領域のみを“1”とし、他は“0”とする如き例で
ある。そして“1”となっている部位のみを(数1)使
って陰影処理を行い、Zバッファ処理を行う。多階調の
画像とは、2値化する前の画像である。
示のための第2の先願である特許出願を行っている(特
願平5−259581号)。この先願は、腸などの折れ
曲がった画像を動画として表示しようとするものであ
り、視点を前進させてその腸に沿って内視鏡であたかも
観察するかの如き表示を行わせる。前進ではなく、後退
して表示させる例も開示する。更に、腸の折れ曲がりに
沿って追跡してゆく例も開示する。これらの動画を得る
ために、三次元原画像を平行投影法で投影して投影画像
を得、動画表示に際してはそれらの投影画像を動画とし
て選択するやり方も開示する。
は、視点位置をマウスやトラックボールを使って人間の
指示によって動かす例を開示する。これによって腸など
の内部を内視鏡を動かしながら見ているような投影画像
が得られる。しかし、第1の先願は、視点位置の移動に
あたって、忠実に腸などの内部を内視鏡で追従させるや
り方については開示していない。
号は、内視鏡で観察する如き投影画像の動画表示につい
て開示し、その投影画像の平行投影は述べているが、中
心投影法による記載はない。
あって、その目的は、点としての視点位置を腸などの内
部に沿って忠実に移動可能にする疑似三次元画像形成法
及び装置を提供するものである。
て視点位置の腸などの内部に沿って忠実に移動可能にす
る疑似三次元画像形成法及び装置を提供するものであ
る。
に対して、視点とこの視点からの視線方向と投影面とを
対にして、その三次元原画像の奥行き方向に向かうよう
に更新し、各対となる視点と投影面毎にその間に介在す
る三次元原画像(以下、単位三次元画像)をその視点位
置及び視線方向から投影面に中心投影して、中心投影画
像を得る画像形成方法を開示する。
この視点からの視線方向と投影面とを対にして、その三
次元原画像の奥行き方向に向かうように更新し、各対と
なる視点と投影面毎にその間に介在する単位三次元画像
をその視点位置及び視線方向から投影面に中心投影して
中心投影画像を得、この中心投影画像をその時の視点と
の距離を基準として陰影処理し且つ隠面処理して擬似三
次元画像を得、各更新毎に得た擬似三次元画像を表示す
る画像表示方法を開示する。
この視点からの視線方向と投影面とを対にして、その三
次元原画像の奥行き方向に向かうように更新し、各対と
なる視点と投影面毎にその間に介在する単位三次元画像
を、その視点位置と視線方向とから、該投影面に中心投
影して中心投影画像を得、視点位置から投影面に下した
垂線位置を中心位置として表示範囲を切り出し、この切
り出した中心投影画像をその時の視点との距離を基準と
する陰影処理し且つ隠面処理して擬似三次元画像を得、
各更新毎に得た擬似三次元画像を表示する画像表示方法
を開示する。
この視点からの視線方向と投影面とを対にして、その三
次元原画像の奥行き方向に向かい且つ投影面がその視点
からの視線方向に対して垂直になるように更新し、各対
となる視点と投影面毎にその間に介在する単位三次元画
像をその視点位置及び視線方向から投影面に中心投影し
て、中心投影画像を得る画像形成方法を開示する。
この視点からの視線方向と投影面とを対にして、その三
次元原画像の奥行き方向に向かい且つ投影面がその視点
からの視線方向に対して垂直になるように更新し、各対
となる視点と投影面毎にその間に介在する単位三次元画
像をその視点位置及び視線方向から投影面に中心投影し
て、中心投影画像を得、この中心投影画像をその時の視
点との距離を基準とする陰影処理し且つ隠面処理して擬
似三次元画像を得、各更新毎に得た擬似三次元画像を表
示する画像表示方法を開示する。
この視点からの視線方向と投影面とを対にして、その三
次元原画像の奥行き方向に向かい且つ投影面がその視点
からの視線方向に対して垂直になるように更新し、各対
となる視点と投影面毎にその間に介在する単位三次元画
像を、その視点位置及び視線方向から投影面に中心投影
して、中心投影画像を得、上記垂線位置を中心位置とし
て投影面から表示範囲を切り出し、この切り出した中心
投影画像をその時の視点との距離を基準とする陰影処理
し且つ隠面処理して擬似三次元画像を得、各更新毎に得
た擬似三次元画像を表示する画像表示方法を開示する。
投影面とを対にして、その三次元原画像の奥行き方向に
向かい且つ投影面がその視点からの視線方向に対して垂
直となるように更新し、各対となる視点と投影面毎にそ
の間に介在する単位三次元画像をその視点位置及び視線
方向から投影面に中心投影して中心投影画像を得、この
中心投影画像の各画素をその視点位置からの距離が遠い
程画素濃度を小さくする陰影化処理及び同一画素位置で
重なりあった画素濃度があった時に視点位置からの距離
の最も小さい位置の画素濃度をその画素位置の濃度とす
る隠面化処理を行い、これによって擬似三次元画像を
得、かくして得た更新毎の単位三次元画像の擬似三次元
画像を表示する画像表示方法において、更新毎に、その
更新前の単位三次元画像の隠面化画像の各画素濃度に含
まれる視点からの距離をみてその時の視点位置から最も
遠い画素位置を見つけ出し、この最も遠い画素位置に近
づくように更新後の視点位置及び視線方向を定めるよう
にした画像表示方法を開示する。
この視点からの視線方向と投影面とを対にして、その三
次元原画像の奥行き方向に向かい且つ投影面がその視点
からの視線方向に対して垂直になるように更新し、各対
となる視点と投影面毎にその間に介在する、関心部位を
2値化した単位三次元画像を、その視点位置及び視線方
向から投影面に中心投影して中心投影画像を得、この中
心投影画像の“1”である各画素をその視点位置からの
距離をが遠い程画素濃度を小さくする陰影化処理及び同
一画素位置で重なりあった画素濃度があったときに視点
からの距離の最も小さい位置の画素濃度をその画素位置
の濃度とする隠面化処理を行い、これによって擬似三次
元画像を得、かくして得た更新毎の擬似三次元画像を表
示画面に表示する画像表示方法において、更新毎に、そ
の更新前の擬似三次元画像の各有意の画素の距離をみて
その時の視点位置から最も遠い画素位置を見つけ出し、
この最も遠い画素位置に近づくように更新後の視点位置
及び視線方向を定めるようにした画像表示方法を開示す
る。
向と投影面とを対にして三次元原画像の奥行き方向に向
うように更新する。更に、この視点と投影面との間に介
在する単位三次元画像を視点位置と視線方向から投影面
に中心投影する。これによって三次元原画像に対して、
奥行き方向に内視鏡的な移動が行われ、画像の更新がな
される。更に本発明によれば、視点と視線方向とは隠面
化処理で得た擬似三次元画像の画素に近づくように更新
する。これによって奥行き方向への内視鏡的な視点の移
動がなされる。
し、具体的には、立体的な各位置と各位置毎の画素濃度
とでデータ化された画像を云う。医用画像にあっては、
複数のCT画像(スライス画像とも云う)を積み上げて
得た積み上げ三次元画像、MRIで得た三次元画像、超
音波で得た三次元画像等がある。これらの三次元画像の
データには、計測して得られたものの他に、計測内容か
ら補間して得られた三次元画素データを含む。 中心投影法…図2で述べたやり方をとる。更に、三次元
原画像に対して、点としての視点とこの視点からの視線
方向と投影面とを対にして更新する。 単位三次元画像…更新毎に定まる視点と投影面との間に
介在する三次元原画像を指す。単位三次元画像は三次元
原画像の一部であり、更新毎に単位三次元画像も更新
し、変更される。
の対にしての更新…これは、三次元原画像に対してその
奥行き方向に点としての視点を移しながら更新する、と
の意味である。例えば、奥行き方向に管路(腸や気管)
が形成されている如き三次元原画像に対して、その管路
の奥側を内視鏡的にみたい場合、奥行き方向に画像を追
って観察できるようにすれば都合がよい。その奥行き方
向へ画像を追うために、視点と視線方向と投影面とを対
にして更新してゆくのである。投影面は視線方向に直角
になるように定めることが好ましい。視線方向に直角と
は、視線方向の中心線が投影面の垂線となる如く、投影
面を定めることである。更新ピッチ距離は大きくとると
粗く進み、小さくとると細かく進む。これらは観察目的
に応じて定める。更新ピッチ時間は短い程よいが、更新
毎に所定のデータ処理(視点と視線方向と投影面との更
新処理、陰影化処理、隠面化処理等)を行うだけの時間
は必要である。 擬似三次元画像…単位三次元画像の中心投影画像を陰影
化処理及び隠面化処理して得た二次元配列の画像であ
る。擬似とは、陰影化処理と隠面化処理により、あたか
も三次元の如き奥行き感のあるリアルティで得た画像が
得られるために使った用語である。
(イ)が操作者とのインターフェースを受け持つプログ
ラムの手順、図1(ハ)が単位三次元画像、擬似三次元
画像を得るためのプログラムの手順を示す。図1(ロ)
は、図1(イ)と図1(ハ)のプログラムの処理の中継
を行う共有メモリ9である。先ず、図1(イ)の各ステ
ップを説明する。 ステップ1…擬似三次元画像の初期画像を表示するステ
ップである。この初期画像の表示にあっては、視点と視
線方向と投影面とは初期設定されたものである。この初
期画面には、「自動アイコン、手動アイコン、終了アイ
コン」のアイコンモードと、「複数ポイント選択アイコ
ン」(ポイントとは点の他に領域を含む)とが表示され
ている。アイコンモードにおける「手動アイコン」と
は、これをONすることで手動によって視点、視線方
向、投影面とを対にして更新可能な状態にする表示スイ
ッチである。「自動アイコン」とは、これをONするこ
とで自動によって視点、視線方向、投影面とを対にして
更新可能な状態にする表示スイッチである。「終了アイ
コン」とは、これをONすることで図1(イ)のプログ
ラムが終了することを指す表示スイッチである。「複数
ポイント選択アイコン」とは、奥行き方向に進むべきポ
イントが複数個(例えばP1、P2)存在するときに、そ
の番号を(例えば1、2)表示しておき、どれか1つを
選択させるアイコンである。進むべきポイントが複数個
存在すると、内視鏡的な自動更新が不可となり、どれか
1つを選択させることにしたのある。この選択は操作者
が行う。図4には、その一例を示す。
次元画像101を表示してある。この画像101には、
2つの異なる奥行き方向のポイントP1とP2とが存在す
る例を示してある。例えば、気管支が2つに分岐する部
位の画像や2つの血管が近接して走っている部位の画像
が相当する。こうした画像の表示に併せて、アイコン1
02、103を表示しておく。アイコン102では3つ
のアイコン「自動、手動、終了」の中のいずれかをON
することで、それに該当するフラグが“1”となる。ア
イコン103では、P1を選択しようとすれば「1」を
ONにし、P2を選択しようとすれば「2」をONにす
る。これらのONにより計算機内では、対応するフラグ
が“1”となる。尚、104はマウスである。
による更新を、自動で行うか手動で行うかの指示を操作
者にさせるステップである。1ビットの自動フラグを与
えておき、これが“1”であれば自動、“0”であれば
手動とのやり方をとる。この手動か自動かの指示は、図
4の手動アイコンか自動アイコンかのいずれかをONす
ることで行う。 ステップ3…自動か手動かの判定を、自動フラグが
“1”か“0”かをみることで行う。手動であればステ
ップ4に移り、自動であればステップ5に移る。 ステップ4…手動による視点、視線方向、投影面の更新
を行って擬似三次元画像を得る処理を行う。この手動に
よる処理とはマウスを使っての対話形式による更新を云
う。
イコンがONしているか否かをチェックする。これは終
了フラグが“1”か“0”かでわかる。“1”であれ
ば、処理を終了する。ステップ4内でも終了アイコンが
ONか否かのチェックを行っていることはいうまでもな
い。 ステップ6…共有メモリ9からのステータスの読み込み
を行う。ステータスとは、図1(ハ)のプログラムによ
ってセットされるものであって、奥行き方向が複数個あ
るか否か、複数個ある場合の状態を示すデータである。
具体的には、選択フラグ(複数個あれば“1”、1個で
あれば“0”となるフラグのこと)、選択数、選択肢
1、2、…より成る。選択肢1、2、…とは、図4のP
1、P2、…に相当する。選択数は総数であり、図4では
2となる。 ステップ7…ステータスを画面に表示させて、どの奥行
き方向を選ぶかの選択を操作者に行わせる。 ステップ8…その選択結果を図1(ロ)の共有メモリ9
に書き込む。
する。 ステップ10…共有メモリ9上の自動フラグを読み込
む。 ステップ11…自動フラグが“1”か否かをチェックす
る。“1”であれば図1(ハ)のプログラムの全体が走
り出す。尚、図1(ハ)の開始条件に、自動フラグが
“1”であることを加味してもよい。 ステップ12…共有メモリ9上の終了フラグを読み込
む。 ステップ13、14…終了フラグが“1”であればステ
ップ14で図1(ハ)のプログラムの走りを停止する。
て、その時の視点から最も遠い距離にある画素位置を自
動的に見つける。この疑似三次元画像とは、初期時では
ステップ1で設定した画面上の画像、それ以降において
は更新毎に得られる疑似三次元画像である。ここで、視
点から最も遠い距離とは、視点から最も遠い画面上の位
置との意味ではなく、疑似三次元画像即ち隠面化処理で
得た画像の画素濃度が視点からの距離で示していること
から(正確な表現としては視点からの距離が大きい程画
素濃度を小さくしていること)、この距離が最も大き
い、隠面化処理画像上の画素位置を選び出すとの意であ
る。このステップ15の処理内容は後述の図5に示す通
りである。
ら最も遠い距離が予め規定した一定値以下か否かを自動
的にチェックする。これは更新の結果、三次元原画像の
最も奥まで更新してしまった場合にはそれ以上の奥行き
には画像は存在しないので、このことを距離が一定値以
内か否かで判定するのである。また、着目した管路があ
る画素位置で終了しそれ以上の先は存在しない時にも、
このステップ16の処理でわかる。一定値以下になれ
ば、処理は終了する。
いるか否かを自動的に見つける。複数領域に分かれてい
るとは図4のP1、P2の如き例である。このP1、P2の
如き例を自動的にみつけるには、ステップ15によって
視点からの距離を、疑似三次元画像の全画素について求
めておき、その中の最大遠距離と思われるものが2個以
上存在するときに、複数領域が存在すると判定する。こ
こで、最大距離とはあるしきい値内の値とすることが好
ましい。しきい値を考慮すれば、最大距離とは極大値で
もよい。複数領域に分かれていればステップ19へ移
り、分かれていなければ(即ち1個の管路のみ)、ステ
ップ20へ移る。
は、それに関するステータスを共有メモリ9にセットす
る。ステータスとは、選択フラグを“1”とし、選択数
を領域数分セットし、選択肢をその番号でセットした内
容を指す。このステータスに対しては図1(イ)のステ
ップ7で操作者が選択する。ステップ19Aで選択結果
を読み込む。 ステップ20…図1(イ)ステップ7で選択したステー
タスの選択結果で選ばれた1つの領域又は、単一の領域
の時のその領域上の最大遠距離画素位置に近づくよう
に、視点、視線方向、投影面を自動的に更新する。この
詳細フローは後述の図5に詳述する。
視点、視線方向、投影面に従って、この視点と投影面と
の間に新しく介在することになった単位三次元画像に対
して中心投影処理を行い、中心投影画像を得る。更に、
中心投影画像に対して陰影化処理及び隠面化処理を行
い、擬似三次元画像を得る。この詳細フローは後述の図
6に示す。 ステップ22…ステップ21で算出した擬似三次元画像
を階調ウィンドレベルに従って表示画面に表示する。こ
こで、階調ウィンドレベルレベルとは表示装置の表示階
調と擬似三次元画像の画素階調との関係を規定したもの
である。表示階調数が64階調とした場合、この64階
調で表示できるように擬似三次元画像の画素階調を変換
する。この変換式は擬似三次元階調を横軸で表示し、縦
軸に表示階調を示した場合、一次関数であったり、折れ
線関数であったり、二次関数であったりの種々の関数と
なる。
更新が行われ、その更新の都度、擬似三次元画像が得ら
れ、表示画面に表示される。そして、更新ピッチ時間が
無視できる程小さい時には、更新した擬似三次元画像を
次々に表示すれば、動画的表示となる。この他の動画的
表示を実現するには、上記更新した擬似三次元画像を補
助メモリに格納しておき、先の第2の先願である特願平
5−259581号の発明を利用するやり方がある。
プ20との具体的処理フローを示す図である。簡単化の
ため省略しているが、図1(ハ)のステップ16〜19
Aは、実際上は図5のステップ31と32との間に挿入
するのが好ましい。図5の各ステップを説明する。 ステップ30…擬似三次元画像を形成し、表示画面に表
示するステップである。この画像は、図1(イ)のステ
ップ1の初期画面の画像又は、図1(ハ)のステップ2
1の処理で得た画像のいずれかである。 ステップ31…ステップ30で表示した画像の全画素を
チェックして、その時の視点から奥行き方向への距離R
の最大値又は極大値Rmaxを求める。これは、表示中の
画像は、陰影化処理し隠面化処理した画像であり、各画
素濃度Iには、I=Imax−C・Rの式に従う距離Rが
保存されていることから、この距離RをR=(Imax−
I)/Cで求めることで、最大値又は極大値を求める。
最大値とは奥行き方向に一画素のみが続いている場合に
利用し、奥行き方向に複数画素が続いている場合や複数
領域(図1(ハ)のステップ18)に分かれている場合
には極大値を利用する。
1(イ)のステップ7)にしてこのステップに入る。ス
テップ32では、画面中心から最大値点又は極大値点ま
での距離dX、dYを算出する。距離dX、dYの基準
位置を画面中心としたのは、この擬似三次元画像を算出
したとき、その視点からの視線方向に対して、投影面を
直角になるように設定したことによる。視線方向に対し
て投影面が直角とは、視線が視点からの中心投影線に相
当していることからこの中心投影線の真中の位置にある
中心線に対して、直角になるように投影面を定めたこと
を意味する。更に、中心線は投影面の中心位置に交わる
ようにしてある。この時の投影面を、表示画面そのもの
に対応させておけば、表示画面の中心が視線からの中心
投影の中心線に一致する。従って、画面中心からの距離
dX、dYを算出することで、視線の中心からどれだけ
最大値点又は極大値点がずれているかをみつける。この
ずれ量は新しい視点と視線方向の算出に利用する。尚、
dX、dYのX、Yとは、投影面の座標系である。
示す図である。CT画像は、現視点e1と投影面との間
に介在する単位三次元画像の一部を示すものである。絶
対座標系としてx、y、zを使用し、投影面がX、Y座
標系としている。現視点e1からの放射状に広がる多数
の投影線の中で、その中央に位置する中心線が、投影面
の中心位置に交わり且つその交わり方が直角になるよう
に、投影面を選ぶ。中心線は視点e1からの垂線であ
る。投影面は表示面に一致させたものとする。こうした
状態で、ステップ31、32で最大値点又は極大値点を
求めてdX、dYを算出したとすると、その様子を図7
に示す。最大値点又は極大値点C1′が奥行きの座標位
置であり、これに向けて視点及び視線方向及び投影面を
更新する。この更新のフローを示したのが後述のステッ
プ33〜37である。図7では更新後の新視点をe2で
示している。尚、現視点e1が投影面の中心位置でない
ように見えるが、実際上は投影面の中心位置に現視点e
1が存在する。
在の投影面は図7に示した通りである。これに対して現
視点e1からの最大値点又は極大値点C1′が投影面の中
心位置からずれた。そこで、視点を現視点e1から最大
値点又は極大値点C1′に近づくように、新視点e2を定
める。更に、この視点e2から新視線方向を決定する。
新視線方向は、新視点e2から最大値点又は極大値点
C1′へ中心投影線の中心線が点C1′を通るようにす
る。更に、この中心線と点C1′で直交するように且つ
この新しい投影面の中心位置が点C1′になるように、
新投影面を設定する。
面とxyz座標系との位置関係を示す数値であり、aは
投影面とx軸との交わる点のx軸成分、bは投影面とy
軸との交わる点のy軸成分、cは投影面とz軸との交わ
り点のz軸成分である。更に、図7に関しては、xyz
座標系の原点から投影面に下した垂線を、z−x面に投
影した線がx軸となす角(これをαとする)、及びこの
垂線がx−z面となす角(これをβとする)を求めてお
くことが必要である。これによって、a、b、cと異な
る次元のα、βで投影面とxyz座標系との関係を律す
ることができる。
系での位置を求める。先ず図8において、旧交点C1か
ら新交点C1′への変位分を求めると、以下となる。
画面中央からdX、dYの変位があった場合のxyz座
標系での変位dx、dy、dzを示す。次に、旧交点C
1の座標を(xc1、yc1、zc1)としたとき、近似的新
交点C1′の座標(xc1′、yc1′、zc1′)を下式で
求める。
る。新視点位置e2は、現視点e1と新交点C1′を結ぶ
直線L上であって新交点C1′に近づくような位置に設
定する。新視点位置e2(x1′、y1′、z1′)は下記
で求める。
は直線Lの長さ(長さがhでないので新交点C1′は一
時的使用のための近似点である。後述するように正確な
新直交点は新a、b、cを使って求める値)、Qは交点
C1′の最大値又は極大値Rmaxに対してある割合率mを
乗算して得た値である。
れによって、e1よりも奥側のmで定まる位置に新視点
e2が定まる。
めるステップである。先ず図7の説明で示した角度α、
βの変化量を下記で求める。
る。dαを|β|が90゜に近づく時に小さくする理由
は、図8の(ロ)の視点1から視点2に近づくにつれ
て、dxが変化してもαに対する効果が少なくなるから
である(小さくしたいこと)。 ステップ36…角度α、βの前回の値をα1、 β1とす
ると、今回の変化分による角度α、βは
方向の中心線に直交する投影面を求める。投影面は図7
に示したようにa、b、cで表わされることから、この
a、b、cを求める処理がステップ37である。即ち、
新視点(x1′、y1′、z1′)から視線方向がα、β
の直線と距離hで直交する投影面のx、y、z軸と交わ
る点a、b、cは下記となる。
(数12)により求まる。前述したように以前に求めた
C1′は一時的使用のための新交点である。
フローである。即ち、視点と視線方向と投影面とが定ま
ると、視点と投影面との間に介在する、新たな三次元原
画像の一部の画像が、新しい単位三次元画像となる。こ
の新しい単位三次元画像をその視点から視線方向に向け
て中心投影し、その中心投影画像を陰影化処理し、隠面
化処理し、新しい擬似三次元画像を得る。これを実現す
るのが図6の処理フローである。 ステップ40…投影面のメモリの全画素濃度をゼロクリ
アする(初期化する)。 ステップ41…投影面から表示画面を切り出す。投影面
のサイズと表示画面のサイズとが同一であれば、投影面
は表示画面と一致する。しかし、投影面のサイズは、表
示画面のサイズより大きくしている例が多く、その場合
には、投影面から表示画面相当分の切り出しが必要とな
る。また、両サイズが一致している場合でも、投影面の
一部の領域のみを表示画面に表示したいこともある。こ
うした場合も、投影面からの切り出しが必要となる。そ
こで、ステップ41では、切り出し中心位置(Xc、
Yc)を与えて、この中心から表示画面の切り出しサイ
ズ(例えば表示画面の大きさの切り出しであれば画面サ
イズ)相当分を切り出すことにした。ここで切り出し中
心位置を交点C1′とすれば、奥行き方向が表示画面の
中心となり、好ましい。このステップ41を設けたこと
で、表示メモリのアドレス(X、Y)が指定可能とな
る。最初に指定するアドレス(X、Y)は、この切り出
しサイズ相当の投影面領域の左上である。アドレス更新
は左から右、上から下へのラスタスキャンに従って行う
(ステップ48、50)。
めのスキャン開始ステップである。以下、単位三次元画
像を、CT画像の積み上げ三次元原画像の一部とした例
で説明する。この例を図9に示す。新視点eに対して投
影面21が与えられ、その間に単位三次元画像が介在す
る。単位三次元画像は、CT画像を積み上げて得た三次
元原画像の一部であり、複数枚のCT画像を積み上げて
得たものである。そして、その中の視点に最も近い1枚
のCT画像23(#1)を先ず選ぶ。この際積み上げC
T画像は、絶対座標系x、y、zに対してy軸に平行に
積み上げたものとしている。この視点に最も近いCT断
層像は、CT断層像番号がわかることから、y成分もC
T断層番号から自動的にわかる。これをy0とする。
の画素の中から、表示メモリのアドレス(X、Y)に投
影される、断層像の対応する画素位置(x0、y0、
z0)を算出する。但し、y0はスライス像をy軸上に並
べたので、予めわかっている(計測時には「スライス間
隔」として相対的な値がわかる)。従って実際にはx0
とz0とを算する。この算出は、中心投影変換式(数
9)によって行う。表示メモリの全アドレス(X、Y)
について、対応する画素位置(x0、y0、z0)を算出
する。ここで、図9において、投影面21上のP点が表
示メモリのアドレス(X、Y)に相当し、断層像23上
のS点が、対応する断層像の画素位置に相当する。但
し、図9では投影面上の任意点Pは、絶対座標系xyz
の座標(x、y、z)で示している。また、C1点(x
c1、yc1、zc1)は投影面23への視点eからの垂線交
点である。その他、α、β、a、b、cは図7の定義と
同じである。中心投影変換式は、大きく2つに分かれ
る。第1は、X、Y座標系からxyz系への変換式であ
り、第2はS点の座標の決定式である。X、Y系からx
yz系への変換が(数9)である。
9)の(2)にyとして代入し、x、zを求め、且つそ
れからyを求めると(数9)の(3)が得られる。更
に、投影面21は、
z)を通る直線22は、
えば視点e(x1、y1、z1)から投影面21に下した
垂線と投影面21との交わる点(この点と視点e間の距
離h)として、
示画面(図示せず)上に、縦512画素×横512画素
で表示するとき、X、Yは−256から+256までの
値を取る。それぞれのX、Yに対して上掲(数9)式に
よりx、y、zが決まる。e点のx1、y1、z1は手動
モードの時、キー入力等により任意に与えるもので、
(数13)式により、y0=d0の断層像上で画素S点の
座標x0、z0が決まる。これがx0、z0の決定式であ
る。
像は複数(#1、#2、…)であって、d0も複数個で
あるので、1組のX、Yに対して複数の投影すべき点x
0、z0が決まる。この中から1つの断層像の投影すべき
点x0、z0を選ぶ。選び方は、Zバッファ法(但し、Z
成分の代わりに距離Rを使う)を使ってもよいが、本実
施例では別のやり方を使う。そこで、ステップ44、4
5、46を設けておく。
では画素点(x0、z0)のしきい値処理をする。これ
は、関心部位を抽出するためであり、しきい値はその抽
出用のしきい値である。関心部位とは、例えば臓器の種
類(腸とか気管支とか)であり、臓器抽出用のしきい値
を与えておくことで、関心臓器の自動抽出がなされる。
その画素点(x0、z0)でしきい値を満足しなければ、
ステップ45で次のCT画像(即ち、視点からみて次に
近いCT画像)#2を指定し、ステップ41で指定した
アドレス(X、Y)に対応するy=y0での画素点
(x0、z0)をステップ43で見つけ出し(y0は指定
CT画像から自動算出)、ステップ44で再びしきい値
以上の画素濃度か否かを判定する。しきい値の範囲外の
画素濃度であれば再びステップ45で次に近いCT画像
を指定する。以下、ステップ41で指示した(X、Y)
に投影する、全CT画像にしきい値範囲内の画素濃度が
なければ、ステップ46を経てステップ47へ移る。し
きい値範囲内の画素濃度があるCT画像#iのy=y0
での画素点(x0、z0)が存在すれば、それ以上の奥側
のCT画像#(i+1)、#(i+2)、…はサーチし
ない。そしてステップ47へと移る。
点(x1′、y1′、z1′)とy=y0での画素点
(x0、z0)との距離Rを下記で求める。
ス法による濃度であるが同時に一種のZバッファで処理
された隠面処理後の濃度である。かくして、1点(X、
Y)についての陰影処理及び隠面処理がなされ、これを
メモリアドレス(X、Y)に擬似三次元画像の画素濃度
として格納する。 ステップ48…表示切り出しアドレスXを+1更新す
る。 ステップ49…Xが表示切り出しアドレスの最大値にな
ったかをチェックする。最大値に達していなければ、ス
テップ42へ戻る。最大値に達していれば、ステップ5
0へ移る。 ステップ50…表示切り出しアドレスのXが最大値に達
したことにより、アドレスを1行更新すべく、表示切り
出しアドレスのYを1だけ更新する。アドレスXも初期
値へと戻す。 ステップ51…表示切り出しアドレスのYが最大値に達
したか否かをチェックし、達していなければステップ4
2に戻り、達していれば全体の処理を終了する。
の断層面23A〜23Eより成る三次元原画像(これは
積み上げ三次元原画像である)の、中心投影を説明する
図である。複数の断層面23A〜23Eの中で中心投影
に利用するのは、断層面23A〜23Dであり、これが
単位三次元画像となる。23Eは視点eよりも外側(手
前側)にあるため、中心投影の対象からはずす。今、断
層面23BにはB1、B2、B3の如き映像があり、断層
面23CにはC1、C2の如き映像があるものとする。こ
こで、映像とは、着目臓器であり、その臓器抽出用のし
きい値を与えて、しきい値範囲内の濃度の画素を着目臓
器として抽出し、その抽出画素に“1”を与えた画像で
ある。抽出臓器以外の画素は“0”にする。視点eから
の視線方向を決定し、これに直交する投影面を決定す
る。視線方向は、視点から投影面へ中心投影する方向で
あり、具体的には、視線方向の中心線0が投影面の垂線
となるように、視線方向を定める。この視線方向と投影
面との関係は、視点と視線方向と投影面とが更新されて
も、その更新毎に成り立つようにしている。即ち、視点
から視線方向を定め、この定まった視線方向の中心線が
投影面の垂線となるように投影面を定めることにしたも
のである。
心投影像は、投影面21上では図のように、視点から放
射状に投影して得られる。B1→B1′、B2→B2′、B
3→B3′、C1→C1′、C2→C2′となって、投影像B
1′、B2′、B3′、C1′、C2が得られる。尚、図で
B3′、C2′とが若干の距離を隔てて表記しているが、
これは、投影像をわかりやすくするためであって、実際
上は重なったものとなる。尚、図10で使用したC1、
C1′は図8、9で使用したものとは意味が異なる。更
に、図10(イ)で断層面23A〜23Eは計測CT断
層像の面の他に、補間等の計算処理で求めた断層面を含
む。
3Eと投影面21とが図3のようにある傾きを持った関
係にある場合の図である。図2のような単純な関係にす
るには、断層面23A〜23Eを、視線方向の中心線0
に直交するような断層面23a〜23eとし、投影面2
1に投影させる。断層面23a〜23eは、断層面23
A〜23E(及び他の断層面も含まれる)からの計算処
理で求める。この計算処理には、断層面の画素位置計算
とその画素位置への画素濃度計算とがある。画素濃度計
算は、補間処理を利用する。
での初期画面の時に用いられる。視点と投影面とが奥行
き方向に更新してゆく本実施例にあっては、この初期画
面以降の更新毎の、視点と投影面と単位三次元画像とは
ほとんどが図10(ロ)の如き関係となる。
奥行き方向に視点と投影面とを更新してゆくことの模式
図である。最初に視点1と視線方向1と投影面1とを与
えて、投影面1に視点1からの中心投影像を得る。この
中心投影像は投影面1に挟まれた単位三次元画像1を中
心投影したものである。次に、視点1よりも奥行き方向
に近づく位置に視点2を与え、視線方向2を求め、対向
する投影面2に中心投影像を得る。以下、同様にして視
点と視線方向と投影面とを、画像の奥行き方向に更新し
てゆく。視点iはその任意の更新位置を示している。更
に、視点jは、管路が2つに分岐した例を示す。これは
図4に示した例に該当する。どちらか一方を選択するこ
とになる。図1で説明した通りである。
である。投影面とは、数学的には投影されるべき面であ
るが、画像処理上では、二次元配列のバッファメモリを
該当させている。一方、表示面とはCRT等の表示画面
であるが、実際上は、二次元配列の表示メモリの内容を
そのままスキャンして表示することから、この二次元配
列の表示メモリが該当するとみてよい。投影面である二
次元配列のバッファメモリと二次元配列の表示メモリと
は、そのサイズ(縦×横の大きさ)によって、以下の如
き使い方がある。 (イ)、第1は、バッファメモリサイズと表示メモリサ
イズとが同一である場合。これは、図12(イ)に示す
例であり、実線の投影面と点線の表示面とが互いに全く
重なった関係にある。この場合には、投影面に得た中心
投影像をすべて表示面に表示させて使う。尚、同一サイ
ズでも表示切取りサイズがバッファメモリの一部であっ
て、これを表示メモリの対応した一部又は対応しない一
部に表示させる例もある。 (ロ)、第2は、バッファメモリサイズが表示メモリサ
イズよりも大きい場合。これは、図12(ロ)に示す例
である。この場合には、表示すべき対象を選択すること
になる。表示面1とは左上側を選択した例、表示面2と
は中央よりも右下側を選択した例である。中央付近を選
択する例もありうる。この選択は、自動も手動もありう
る。
影像は、単位三次元画像をそのまま投影したものではな
く、単位三次元画像に陰影処理を施して得た陰影処理画
像である。この陰影処理とは、視点から投影対象となる
画素位置までの距離Rを求め、この距離Rに逆比例する
ように、その画素位置の投影濃度Iを決定する処理であ
る。例えば(数1)で計算する。
バッファには、その画素位置に、2つ以上の濃度が重な
って投影されることがある。図10(イ)のB3′と
C2′とが一部重複している画素位置の如き場合であ
る。こうした場合、距離Rの手前の画素位置の濃度を残
し、その他の遠い画素位置のものは消去する。これは一
種のZバッファ法である。図10(イ)の例では、重複
している画素位置にあってはC2′が手前の距離のた
め、C2′を残すことになる。こうしたZバッファ法で
処理された二次元配列の画像が最終的に投影面に相当す
るバッファメモリに得られる。こうした最終的に投影面
に相当するバッファメモリに得られた二次元配列画像が
表示対象となり、図12(イ)や(ロ)に従って表示さ
れる。
態様がある。 (1)、投影面の位置…投影面の中心位置を図7のxy
z座標系のどの位置におくかを決めるやり方がある。そ
の他に、投影面の右上端位置等の特徴的な位置をxyz
座標系のどの位置におくかを決めるやり方もある。いず
れでも、新視点とこの位置との距離が一定値になるよう
に決める。かくして更新毎に、その時の新視点が得られ
れば、自動的に投影面の位置が定まる。 (2)、投影面の傾き…どの位置を基準にして傾きを定
義するかが前提である。この基準位置は、(1)で述べ
た投影面の位置を使うことが好ましい。投影面の傾き
は、本来任意の値であってよいが、視線方向に直交する
ような値とすることが好ましい。
(又は極大値点)について追記する。最も遠い部分が点
ではなく一定の広がりを持った領域や面として得られて
いる例がある。こうした場合、最も遠い点がその領域や
面に沿って複数個存在することが考えられる。そこで、
統計的な処理を行って、その領域や面の中の中心位置や
重心位置を求め、これを最も遠い点として決めるやり方
もある。例えば、図4で領域P1とP2に大きな極大値
点があり、領域P1付近をさらに細かくみると極大値点
が複数個(P11、P12、P13、……)存在し、領
域P2付近にも極大値点が複数個(P21、P22、P
23、……)存在する場合(いずれも図示せず)であ
る。極大値点の数を減らす一方法は画像の平滑化であ
る。これを次に記す。
の合計15画素の画素値をくわえて25で割り、5x5
画素領域の中心画素の新たな画素値とする(もちろん、
平滑化まえのメモリと平滑化後のメモリは異なる)こと
で得られる。極大値点の数が減ったところで、極大値点
の間の距離を求める。この相互の距離のうち、予め設定
した値(例えば、気管支を見ている場合は、見たい気管
支の直径の1cm程度)より小さいものがあれば、さら
に平滑化をする。前記のすべての距離が、予め設定した
値より大きければ、そこで改めて極大値点が幾つかある
かを調べる。一つなら自動更新が可能であり、複数個な
ら操作者に問い合わせる。また最初に見つかった任意の
一点のみを選んで最も遠い点とするやり方もある。関心
領域について述べる。関心領域には以下の態様がある。 (1)、臓器そのものを関心領域とする例である。 (2)、病巣部を関心領域とする例である。 (3)、表示画面に切り取る対象を、関心領域とする例
である。
では、三次元原画像を表示させながらマウスやキーボー
ドを使って、操作員の指示に従って行う。例えば、図5
のステップ32のdX、dYを手動で入力する。または
もっと簡単にα、βとhをキー入力する。尚、陰影づけ
のアルゴリズムとしては、ボリュームレンダリング法
(例えば Maec Levoy著「Display of Surfaces from Vo
lume Date」、IEEE Computer Graphics & Applications
May 1988、29〜37P)を用いてもよい。これは、第1の
先願である特願平6−3492号の図6で述べたもので
ある。
示す。図13において、CPU1は、全体の管理及び更
新処理(図1、図5、図6)を行うものであり、主メモ
リ2はそれに必要なプログラムのラッチ及び各種の作業
用データを記憶する。磁気ディスク3は、三次元原画像
を記憶しており、これをCPU1に送り更新処理の対象
とする。また更新の処理後の擬似三次元画像を含む各種
の処理経過画像が再利用すべき各種の画像を記憶する。
マウス8は、コントローラ7を介してCRT上にマウス
の表示を行わせ、自動更新や手動更新に利用される。表
示メモリ5は、CRT6に表示するための画像や各種マ
イコンをラッチする。共通バス4はそれらの共通インタ
ーフェース線である。
図8で説明する。図14(イ)でdYを変化させるとβ
が変化し、断層像(CT像)をG3、G4方向に傾けたと
同じ効果が得られる(β変化する)。同じ効果をdXに
も持たせたい。本来、dXは画像をy軸の回りに回転さ
せる(αが変化する)効果を持つが、全く別の意味を持
たせたい場合がある。即ち、dXを変化させたとき
G3、G4とは直角方向のG1、G2方向に傾ける効果を持
たせたいときがある。しかも画像は、あまり回転させた
くない。これに対する、近似的で便宜的な手段を以下に
説明する。例えば、dY=0、dXがゼロでないとき、
G1、G2方向に傾けるには、投影面をα方向に90度回
転させて(βはそのまま)、dYの大きさをdXと等し
くし、強制的にdX=0と置き換えると、dXはゼロな
のでy軸の回りに回せず(αはそのまま)、しかもdY
は大きさがゼロでなくなり方向がG1、G2の方向をむく
のでG1、G2の方向に傾いた(βが変化する)疑似三次
元画像が得られることになる。但し、上記のように投影
面を回転させて、見る方向を90度回転させたので、疑
似三次元画像の構成後にx=y=z=0の投影点のまわ
りに回転して元に戻さなければならない。通常はdX、
dYともにゼロでないので、以下のように近似的な処理
をする。図14(ロ)のように、dYをsqrt(dX
2+dY2)、dXをゼロと置き換えて、同時にα方向に
−ξだけ回転しておいて(ここで、ξ=arctan
(dX/dY)である)、疑似三次元画像を構成する。
これによりG5、G6方向の傾きが含まれることになる。
但し、−ξだけ角度を回転させているので、疑似三次元
画像を構成したあとで、原点x=y=z=0の投影点の
まわりに+ξだけ回転して元に戻す必要がある。図15
にはこの手順を示す。
更新例としたが、これとは逆に奥から手前の視点の更新
をも、ほぼ同じ考え方で可能であることは云うまでもな
い。
図である。人間の気管支の断層像(CT画像)#1〜#
kを用意しておき、これに気管支の手前から内視鏡的に
近づき、分岐した気管支A、Bの一方の気管支B内を更
に内視鏡的に進ときの様子を示している。視点1→2→
3が断層面に直角方向に視点を更新した例であり、視点
4→5が気管支B内を進むときの視点の更新例を示す。
これらの視点1〜5に合わせて投影面1〜5も更新され
る。
5による具体的な表示画像例を示す。尚、図で補間とは
補間処理で得た画像とのことであり、図に示す個所以外
にも図の精度上から適当に作られて表示に使われる。
が、腸などの内部に沿って忠実に移動なように更新する
ことができるようになった。更に視点からの中心投影を
行わせたことで、腸などの内部であたかも移動するよう
な擬似三次元画像を得ることができた。
る。
ーチャートである。
である。
ある。
る。
Claims (16)
- 【請求項1】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かうように更新し、各対となる視点
と投影面毎にその間に介在する三次元原画像(以下、単
位三次元画像)をその視点位置及び視線方向から投影面
に中心投影して、中心投影画像を得る画像形成方法。 - 【請求項2】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かうように更新し、各対となる視点
と投影面毎にその間に介在する単位三次元画像をその視
点位置及び視線方向から投影面に中心投影して中心投影
画像を得、この中心投影画像をその時の視点との距離を
基準として陰影処理し且つ隠面処理して擬似三次元画像
を得、各更新毎に得た擬似三次元画像を表示する画像表
示方法。 - 【請求項3】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かうように更新し、各対となる視点
と投影面毎にその間に介在する単位三次元画像を、その
視点位置と視線方向とから、該投影面に中心投影して中
心投影画像を得、視点位置から投影面に下した垂線位置
を中心位置として表示範囲を切り出し、この切り出した
中心投影画像をその時の視点との距離を基準とする陰影
処理し且つ隠面処理して擬似三次元画像を得、各更新毎
に得た擬似三次元画像を表示する画像表示方法。 - 【請求項4】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かい且つ投影面がその視点からの視
線方向に対して垂直になるように更新し、各対となる視
点と投影面毎にその間に介在する単位三次元画像をその
視点位置及び視線方向から投影面に中心投影して、中心
投影画像を得る画像形成方法。 - 【請求項5】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かい且つ投影面がその視点からの視
線方向に対して垂直になるように更新し、各対となる視
点と投影面毎にその間に介在する単位三次元画像をその
視点位置及び視線方向から投影面に中心投影して、中心
投影画像を得、この中心投影画像をその時の視点との距
離を基準とする陰影処理し且つ隠面処理して擬似三次元
画像を得、各更新毎に得た擬似三次元画像を表示する画
像表示方法。 - 【請求項6】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かい且つ投影面がその視点からの視
線方向に対して垂直になるように更新し、各対となる視
点と投影面毎にその間に介在する単位三次元画像を、そ
の視点位置及び視線方向から投影面に中心投影して、中
心投影画像を得、上記垂線位置を中心位置として投影面
から表示範囲を切り出し、この切り出した中心投影画像
をその時の視点との距離を基準とする陰影処理し且つ隠
面処理して擬似三次元画像を得、各更新毎に得た擬似三
次元画像を表示する画像表示方法。 - 【請求項7】 請求項2、3、5、6のいずれかの画像
表示方法において、上記陰影化処理は視点と画素位置と
の距離が大きい程に小さい画素濃度とするデプス法によ
って行い、隠面化処理は画素に複数の濃度が重なった場
合、視点と画素位置との距離が最も小さい位置の画素濃
度をその画素位置の濃度として残すZバッファ法によっ
て行うものとした画像表示方法。 - 【請求項8】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かうように更新し、各対となる視点
と投影面毎にその間に介在する単位三次元画像をその視
点位置及び視線方向から投影面に中心投影して、中心投
影像を得る画像形成方法において、各更新にあっては、
奥行き方向が2つ以上存在するときにはその一方を選択
して更新を行わせるようにした画像形成方法。 - 【請求項9】 三次元原画像に対して、視点とこの視点
からの視線方向と投影面とを対にして、その三次元原画
像の奥行き方向に向かい且つ投影面がその視点からの視
線方向に対して垂直となるように更新し、各対となる視
点と投影面毎にその間に介在する単位三次元画像をその
視点位置及び視線方向から投影面に中心投影して中心投
影画像を得、この中心投影画像をその時の視点からの距
離が遠い程画素濃度を小さくする陰影処理を行って陰影
化画像を得る画像形成方法において、 更新毎に、その更新前の単位三次元画像の陰影化画像の
各画素濃度に含まれる視点からの距離をみてその時の視
点位置から最も遠い画素位置を見つけ出し、この最も遠
い画素位置に近づくように更新後の視点位置及び視線方
向を定めるようにした画像形成方法。 - 【請求項10】 三次元原画像に対して、視点とこの視
点からの視線方向と投影面とを対にして、その三次元原
画像の奥行き方向に向かい且つ投影面がその視点からの
視線方向に対して垂直となるように更新し、各対となる
視点と投影面毎にその間に介在する単位三次元画像をそ
の視点位置及び視線方向から投影面に中心投影して中心
投影画像を得、この中心投影画像の関心部位に対してそ
の時の視点からの距離が遠い程画素濃度を小さくする陰
影化処理を行って陰影化画像を得る画像形成方法におい
て、 更新毎に、その更新前の単位三次元画像の陰影化画像の
関心部位の各画素濃度に含まれる視点からの距離をみて
その時の視点位置から最も遠い画素位置を見つけ出し、
この最も遠い画素位置に近づくように更新後の視点位置
及び視線方向を定めるようにした画像形成方法。 - 【請求項11】 三次元原画像に対して、視点と投影面
とを対にして、その三次元原画像の奥行き方向に向い且
つ投影面がその視点からの視線方向に対して垂直となる
ように更新し、各対となる視点と投影面毎にその間に介
在する単位三次元画像をその視点位置及び視線方向から
投影面に中心投影して中心投影画像を得、この中心投影
画像の各画素をその視点位置からの距離が遠い程画素濃
度を小さくする陰影化処理及び同一画素位置で重なりあ
った画素濃度があった時に視点位置からの距離の最も小
さい位置の画素濃度をその画素位置の濃度とする隠面化
処理を行い、これによって擬似三次元画像を得、かくし
て得た更新毎の単位三次元画像の擬似三次元画像を表示
する画像表示方法において、 更新毎に、その更新前の単位三次元画像の陰影化画像の
各画素濃度に含まれる視点からの距離をみてその時の視
点位置から最も遠い画素位置を見つけ出し、この最も遠
い画素位置に近づくように更新後の視点位置及び視線方
向を定めるようにした画像表示方法。 - 【請求項12】 三次元原画像に対して、視点とこの視
点からの視線方向と投影面とを対にして、その三次元原
画像の奥行き方向に向かい且つ投影面がその視点からの
視線方向に対して垂直になるように更新し、各対となる
視点と投影面毎にその間に介在する、関心部位を2値化
した単位三次元画像を、その視点位置及び視線方向から
投影面に中心投影して中心投影画像を得、この中心投影
画像の“1”である各画素をその視点位置からの距離が
遠い程画素濃度を小さくする陰影化処理及び同一画素位
置で重なりあった画素濃度があった時に視点からの距離
の最も小さい位置の画素濃度をその画素位置の濃度とす
る隠面化処理を行い、これによって擬似三次元画像を
得、かくして得た更新毎の擬似三次元画像を表示画面に
表示する画像表示方法において、 更新毎に、その更新前の擬似三次元画像の各有意の画素
の距離をみてその時の視点位置から最も遠い画素位置を
見つけ出し、この最も遠い画素位置に近づくように更新
後の視点位置及び視線方向を定めるようにした画像表示
方法。 - 【請求項13】 三次元原画像に対して、視点とこの視
点からの視線方向と投影面とを対にして、その三次元原
画像の奥行き方向に向かうように更新する手段と、各対
となる視点と投影面毎にその間に介在する三次元原画像
(以下、単位三次元画像)をその視点位置及び視線方向
から投影面に中心投影して、中心投影画像を得る手段
と、より成る画像形成装置。 - 【請求項14】 三次元原画像に対して、視点とこの視
点からの視線方向と投影面とを対にして、その三次元原
画像の奥行き方向に向かうように更新する手段と、各対
となる視点と投影面毎にその間に介在する単位三次元画
像を、その視点位置と視線方向とから、該投影面に中心
投影して中心投影画像を得る手段と、視点位置から投影
面に下した垂線位置を中心位置として表示範囲を切り出
し、この切り出した中心投影画像をその時の視点との距
離を基準とする陰影処理し且つ隠面処理して擬似三次元
画像を得る手段と、各更新毎に得た擬似三次元画像を表
示する手段と、より成る画像表示装置。 - 【請求項15】 三次元原画像に対して、視点とこの視
点からの視線方向と投影面とを対にして、その三次元原
画像の奥行き方向に向かい且つ投影面がその視点からの
視線方向に対して垂直になるように更新する手段と、各
対となる視点と投影面毎にその間に介在する単位三次元
画像をその視点位置及び視線方向から投影面に中心投影
して、中心投影画像を得る手段と、より成る画像形成装
置。 - 【請求項16】 三次元原画像に対して、視点と投影面
とを対にして、その三次元原画像の奥行き方向に向かい
且つ投影面がその視点からの視線方向に対して垂直にな
るように更新する手段と、各対となる視点と投影面毎に
その間に介在する単位三次元画像をその視点位置及び視
線方向から投影面に中心投影して中心投影画像を得る手
段と、この中心投影画像の各画素をその視点位置からの
距離が遠い程画素濃度を小さくする陰影化処理及び同一
画素で重なりあった画素濃度があったときに視点位置か
らの距離の最も小さい位置の画素濃度をその画素位置の
濃度とする隠面化処理を行い、これによって擬似三次元
画像を得る手段と、かくして得た更新毎の単位三次元画
像の擬似三次元画像を表示する手段と、より成る画像表
示装置において、 更新毎に、その更新前の単位三次元画像の隠面化画像の
各画素濃度に含まれる視点からの距離をみてその時の視
点位置から最も遠い画素位置を見つけ出し、この最も遠
い画素位置に近づくように更新後の視点位置及び視線方
向を定めるようにした手段を備えて成る画像表示装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP14349694A JP3632862B2 (ja) | 1994-06-24 | 1994-06-24 | 三次元画像表示方法及び装置 |
US08/374,088 US5694530A (en) | 1994-01-18 | 1995-01-18 | Method of constructing three-dimensional image according to central projection method and apparatus for same |
US08/708,959 US5900878A (en) | 1994-01-18 | 1996-09-06 | 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 |
US08/856,675 US5953013A (en) | 1994-01-18 | 1997-04-29 | Method of constructing three-dimensional image according to central projection method and apparatus for same |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP14349694A JP3632862B2 (ja) | 1994-06-24 | 1994-06-24 | 三次元画像表示方法及び装置 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002262698A Division JP3485913B2 (ja) | 2002-09-09 | 2002-09-09 | 画像表示方法及び装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0816813A true JPH0816813A (ja) | 1996-01-19 |
JP3632862B2 JP3632862B2 (ja) | 2005-03-23 |
Family
ID=15340069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP14349694A Expired - Lifetime JP3632862B2 (ja) | 1994-01-18 | 1994-06-24 | 三次元画像表示方法及び装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3632862B2 (ja) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998015226A1 (fr) * | 1996-10-08 | 1998-04-16 | Hitachi Medical Corporation | Procede et appareil de formation et d'affichage d'une image a partir d'une pluralite d'images partielles |
US5883933A (en) * | 1996-03-29 | 1999-03-16 | Hitachi Medical Corporation | Method and apparatus for displaying three-dimensional image |
JP2001052195A (ja) * | 1999-06-21 | 2001-02-23 | Ge Medical Syst Sa | 3次元画像の一部を視覚化する方法 |
JP2001118086A (ja) * | 1999-10-21 | 2001-04-27 | Hitachi Medical Corp | 画像表示方法および装置 |
JP2001351093A (ja) * | 2000-04-05 | 2001-12-21 | Hitachi Medical Corp | 医用画像表示装置 |
US6411298B1 (en) | 1996-06-25 | 2002-06-25 | Hitachi Medical Corporation | Method and apparatus for determining visual point and direction of line of sight in three-dimensional image construction method |
WO2003005902A1 (fr) | 2001-07-12 | 2003-01-23 | Hitachi Medical Corporation | Procede de saisie d'images endoscopiques et dispositif d'irm l'utilisant |
JP2006174440A (ja) * | 2000-02-21 | 2006-06-29 | Sharp Corp | 画像合成装置 |
EP1681011A1 (en) * | 2003-10-31 | 2006-07-19 | Olympus Corporation | Insertion support system |
JP2006246941A (ja) * | 2005-03-08 | 2006-09-21 | Toshiba Corp | 画像処理装置及び管走行トラッキング方法 |
JP2009022411A (ja) * | 2007-07-18 | 2009-02-05 | Hitachi Medical Corp | 医用画像処理装置 |
-
1994
- 1994-06-24 JP JP14349694A patent/JP3632862B2/ja not_active Expired - Lifetime
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5883933A (en) * | 1996-03-29 | 1999-03-16 | Hitachi Medical Corporation | Method and apparatus for displaying three-dimensional image |
US6411298B1 (en) | 1996-06-25 | 2002-06-25 | Hitachi Medical Corporation | Method and apparatus for determining visual point and direction of line of sight in three-dimensional image construction method |
WO1998015226A1 (fr) * | 1996-10-08 | 1998-04-16 | Hitachi Medical Corporation | Procede et appareil de formation et d'affichage d'une image a partir d'une pluralite d'images partielles |
US6990231B2 (en) | 1996-10-08 | 2006-01-24 | Hitachi Medical Corporation | Method and apparatus for forming and displaying projection image from a plurality of sectional images |
JP2001052195A (ja) * | 1999-06-21 | 2001-02-23 | Ge Medical Syst Sa | 3次元画像の一部を視覚化する方法 |
JP2011216122A (ja) * | 1999-06-21 | 2011-10-27 | Ge Medical Syst Sa | 3次元画像の一部を視覚化する方法 |
JP2001118086A (ja) * | 1999-10-21 | 2001-04-27 | Hitachi Medical Corp | 画像表示方法および装置 |
JP4658787B2 (ja) * | 2000-02-21 | 2011-03-23 | シャープ株式会社 | 画像合成装置 |
JP2006174440A (ja) * | 2000-02-21 | 2006-06-29 | Sharp Corp | 画像合成装置 |
JP2001351093A (ja) * | 2000-04-05 | 2001-12-21 | Hitachi Medical Corp | 医用画像表示装置 |
CN100350875C (zh) * | 2001-07-12 | 2007-11-28 | 株式会社日立医药 | 类内窥镜图像获取方法和使用该方法的磁共振成像设备 |
US7653426B2 (en) | 2001-07-12 | 2010-01-26 | Hitachi Medical Corporation | Endoscopic image pickup method and magnetic resonance imaging device using the same |
WO2003005902A1 (fr) | 2001-07-12 | 2003-01-23 | Hitachi Medical Corporation | Procede de saisie d'images endoscopiques et dispositif d'irm l'utilisant |
EP1681011A4 (en) * | 2003-10-31 | 2009-08-19 | Olympus Corp | EINFÜHRSTÜTZSYSTEM |
EP1681011A1 (en) * | 2003-10-31 | 2006-07-19 | Olympus Corporation | Insertion support system |
US8049777B2 (en) | 2003-10-31 | 2011-11-01 | Olympus Corporation | Insertion support system for specifying a location of interest as an arbitrary region and also appropriately setting a navigation leading to the specified region |
EP2245983A3 (en) * | 2003-10-31 | 2013-12-25 | Olympus Corporation | Insertion support system |
JP2006246941A (ja) * | 2005-03-08 | 2006-09-21 | Toshiba Corp | 画像処理装置及び管走行トラッキング方法 |
JP2009022411A (ja) * | 2007-07-18 | 2009-02-05 | Hitachi Medical Corp | 医用画像処理装置 |
Also Published As
Publication number | Publication date |
---|---|
JP3632862B2 (ja) | 2005-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5953013A (en) | Method of constructing three-dimensional image according to central projection method and apparatus for same | |
JP4200546B2 (ja) | 画像表示装置 | |
US5891030A (en) | System for two dimensional and three dimensional imaging of tubular structures in the human body | |
US9053565B2 (en) | Interactive selection of a region of interest in an image | |
JP4584575B2 (ja) | 3d画像に表示される3dサーフェイスとインタラクトする画像処理方法 | |
US6411298B1 (en) | Method and apparatus for determining visual point and direction of line of sight in three-dimensional image construction method | |
US9430879B2 (en) | Methods and apparatuses for creating orientation markers and 3D ultrasound imaging systems using the same | |
JP3632862B2 (ja) | 三次元画像表示方法及び装置 | |
US5900878A (en) | 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 | |
US20240299098A1 (en) | View angle-independent visual representation of a cut procedure | |
WO2005055126A1 (en) | System and method for vascular visualization using planar reformation of vascular central axis surface with biconvex slab | |
JP4515615B2 (ja) | 画像表示装置 | |
JPH1176228A (ja) | 三次元画像構成装置 | |
JP4738236B2 (ja) | 画像表示装置 | |
JPH0981786A (ja) | 3次元画像処理方法 | |
JP3485913B2 (ja) | 画像表示方法及び装置 | |
CN113256818A (zh) | 基于不连续性检测的可度量细粒度遮挡去除可视化方法 | |
JPH08255261A (ja) | 三次元画像表示方法 | |
JP3325635B2 (ja) | 画像処理方法及び装置 | |
JP3851364B2 (ja) | 投影画像表示装置 | |
JPH10201755A (ja) | 擬似三次元画像における三次元量計測方法及びその装置 | |
JP3411745B2 (ja) | 断層像補間法及び装置 | |
JP3597917B2 (ja) | 擬似三次元画像構成法並びに投影画像表示装置 | |
JP2006178772A (ja) | 画像処理方法および装置並びにプログラム | |
JP3749322B2 (ja) | 擬似三次元画像上で距離計測をする距離計測装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20041110 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041216 |
|
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: 20080107 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090107 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100107 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100107 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110107 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110107 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120107 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130107 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140107 Year of fee payment: 9 |
|
EXPY | Cancellation because of completion of term |