JPH057554A - デイスプレイリスト表面データを使用する手術シミユレーシヨン装置 - Google Patents
デイスプレイリスト表面データを使用する手術シミユレーシヨン装置Info
- Publication number
- JPH057554A JPH057554A JP4036801A JP3680192A JPH057554A JP H057554 A JPH057554 A JP H057554A JP 4036801 A JP4036801 A JP 4036801A JP 3680192 A JP3680192 A JP 3680192A JP H057554 A JPH057554 A JP H057554A
- Authority
- JP
- Japan
- Prior art keywords
- point
- display list
- image plane
- coordinates
- cut
- 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
- 238000004088 simulation Methods 0.000 title claims description 13
- 238000005520 cutting process Methods 0.000 claims abstract description 20
- 238000013507 mapping Methods 0.000 claims abstract description 16
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000001356 surgical procedure Methods 0.000 claims description 5
- 239000003086 colorant Substances 0.000 claims description 4
- 238000005286 illumination Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 54
- 239000013598 vector Substances 0.000 description 53
- 230000008569 process Effects 0.000 description 25
- 238000003384 imaging method Methods 0.000 description 7
- 238000002591 computed tomography Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000002452 interceptive effect Effects 0.000 description 5
- 210000000988 bone and bone Anatomy 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 210000001835 viscera Anatomy 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241001486234 Sciota Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Image Generation (AREA)
Abstract
(57)【要約】 （修正有）
【目的】立方体分割法や類似の画像発生方法によって作
成されるような三次元の点のディスプレイリストで記述
される表面を切断する効率的な方法の提供。 【構成】対応する法線をそなえた三次元の各点の第一の
ディスプレイリスト１８として表される表面構造に対す
るカット軸に沿った装置に於いて、各点の第一のディス
プレイリストに作用して各点の座標を回転させる変換手
段、各点を受けて、回転した座標を画像平面２４内の画
素座標を持つ画素２６に写像することにより、表面構造
の二次元表現を作成し、カット軸に平行な線に沿った各
点が画像平面内の単一の画素に写像されるようにする写
像手段、写像手段と通信して多角形領域５０を記述し、
また各点を識別する手段、および各点のいくつかの点を
第二のディスプレイリスト１８に記憶する手段を含むこ
とを特徴とする。
成されるような三次元の点のディスプレイリストで記述
される表面を切断する効率的な方法の提供。 【構成】対応する法線をそなえた三次元の各点の第一の
ディスプレイリスト１８として表される表面構造に対す
るカット軸に沿った装置に於いて、各点の第一のディス
プレイリストに作用して各点の座標を回転させる変換手
段、各点を受けて、回転した座標を画像平面２４内の画
素座標を持つ画素２６に写像することにより、表面構造
の二次元表現を作成し、カット軸に平行な線に沿った各
点が画像平面内の単一の画素に写像されるようにする写
像手段、写像手段と通信して多角形領域５０を記述し、
また各点を識別する手段、および各点のいくつかの点を
第二のディスプレイリスト１８に記憶する手段を含むこ
とを特徴とする。
Description
【０００１】
【産業上の利用分野】本発明は一般に表面データの画像
のディスプレイと操作を行うための装置と方法に関する
ものである。
のディスプレイと操作を行うための装置と方法に関する
ものである。
【０００２】
【従来の技術】従来のＸ線装置では、骨と軟部組織の異
なる吸収特性に基いて二次元の陰影画像が作成される。
診断能力の改善がコンピュータ断層撮影装置（ＣＴ装
置）によって与えられてきた。このコンピュータ断層撮
影装置はＸ線に基いている。また、コンピュータ断層撮
影装置を使うことにより、検査している身体、対象、ま
たは患者の横断スライスを描いた二次元ビュー（ｖｉｅ
ｗ）を作成することができる。多数の隣接するスライス
に対するデータを収集することにより、三次元の情報を
得ることができる。
なる吸収特性に基いて二次元の陰影画像が作成される。
診断能力の改善がコンピュータ断層撮影装置（ＣＴ装
置）によって与えられてきた。このコンピュータ断層撮
影装置はＸ線に基いている。また、コンピュータ断層撮
影装置を使うことにより、検査している身体、対象、ま
たは患者の横断スライスを描いた二次元ビュー（ｖｉｅ
ｗ）を作成することができる。多数の隣接するスライス
に対するデータを収集することにより、三次元の情報を
得ることができる。
【０００３】同様に、より新しい磁気共鳴（ＭＲ−ｍａ
ｇｎｅｔｉｃ ｒｅｓｏｎａｎｃｅ）イメージング技術
は内臓の物理的な性質を表す三次元のデータアレーを発
生することができる。更に、ＭＲ装置は骨と軟部組織と
の間だけでなく、種々の組織種別相互の間の識別を行
う。ＭＲイメージング装置は吸収データだけではなく生
理データを発生することもできる。
ｇｎｅｔｉｃ ｒｅｓｏｎａｎｃｅ）イメージング技術
は内臓の物理的な性質を表す三次元のデータアレーを発
生することができる。更に、ＭＲ装置は骨と軟部組織と
の間だけでなく、種々の組織種別相互の間の識別を行
う。ＭＲイメージング装置は吸収データだけではなく生
理データを発生することもできる。
【０００４】どちらの装置からの三次元情報も一連のス
ライス画像として見ることができ、観測者の推理能力を
用いて種々の内臓の三次元が心に描かれる。代案とし
て、陰影と輪郭を付けた三次元画像をコンピュータで作
成してもよい。三次元画像のディスプレイは主として計
算機利用設計（ＣＡＤ−ｃｏｍｐｕｔｅｒ ａｉｄｅｄ
ｄｅｓｉｇｎ）および計算機利用製造（ＣＡＭ−ｃｏ
ｍｐｕｔｅｒ ａｉｄｅｄ ｍａｎｕｆａｃｔｕｒｉｎ
ｇ）産業により推進されてきた。製造された部品にたい
する立体モデルをディスプレイするため、また回転して
異なる方向からこれらの部品を見るための装置が開発さ
れてきた。一つの装置では、個別画素の集合として画像
がスクリーン等のディスプレイ装置上に作成される。個
別画素の強さと色が深さ情報を与える。
ライス画像として見ることができ、観測者の推理能力を
用いて種々の内臓の三次元が心に描かれる。代案とし
て、陰影と輪郭を付けた三次元画像をコンピュータで作
成してもよい。三次元画像のディスプレイは主として計
算機利用設計（ＣＡＤ−ｃｏｍｐｕｔｅｒ ａｉｄｅｄ
ｄｅｓｉｇｎ）および計算機利用製造（ＣＡＭ−ｃｏ
ｍｐｕｔｅｒ ａｉｄｅｄ ｍａｎｕｆａｃｔｕｒｉｎ
ｇ）産業により推進されてきた。製造された部品にたい
する立体モデルをディスプレイするため、また回転して
異なる方向からこれらの部品を見るための装置が開発さ
れてきた。一つの装置では、個別画素の集合として画像
がスクリーン等のディスプレイ装置上に作成される。個
別画素の強さと色が深さ情報を与える。
【０００５】このような装置では、画像がそれから作成
される三次元データをディスプレイ装置の二次元画像平
面に写像（マッピング）しなければならない。この写像
を行うための一つの方法は「射線追跡（ray tracing
）」と呼ばれ、この方法では観測点のまわりに規則正
しく分散配置された射線に沿って存在するデータ点を調
べなければならない。観測点に対して各射線に沿った一
番先のデータ点により、ディスプレーする表面が決ま
る。このような各データ点に対応する画素が画像平面に
作成される。観測点から一番先のデータ点までの距離お
よび一番先のデータ点での表面法線の射線からのずれを
推定することにより、画素のシェージング（shading ）
が決定される。垂直な平面はそのデータ点で表面から垂
直に伸びる線である。
される三次元データをディスプレイ装置の二次元画像平
面に写像（マッピング）しなければならない。この写像
を行うための一つの方法は「射線追跡（ray tracing
）」と呼ばれ、この方法では観測点のまわりに規則正
しく分散配置された射線に沿って存在するデータ点を調
べなければならない。観測点に対して各射線に沿った一
番先のデータ点により、ディスプレーする表面が決ま
る。このような各データ点に対応する画素が画像平面に
作成される。観測点から一番先のデータ点までの距離お
よび一番先のデータ点での表面法線の射線からのずれを
推定することにより、画素のシェージング（shading ）
が決定される。垂直な平面はそのデータ点で表面から垂
直に伸びる線である。
【０００６】射線追跡法の一つの利点はこれが三次元表
面の高分解能と実際的なシェージングを提供するという
ことである。一つの欠点は画像の回転によって生じ得る
ような観測点の相対運動では通常、アレーの各データ点
に対して射線の再計算が必要となる。射線追跡計算を反
復するために必要な時間により、ユーザが画像を操作す
る速度が遅くなる。したがって、ユーザとイメージング
装置との間の相互作用が制限される。
面の高分解能と実際的なシェージングを提供するという
ことである。一つの欠点は画像の回転によって生じ得る
ような観測点の相対運動では通常、アレーの各データ点
に対して射線の再計算が必要となる。射線追跡計算を反
復するために必要な時間により、ユーザが画像を操作す
る速度が遅くなる。したがって、ユーザとイメージング
装置との間の相互作用が制限される。
【０００７】三次元データを画像平面に写像し、各画素
のシェージングを決定する改善された方法が、米国特許
第４，７１９，５８５号に説明されている。そこに開示
された「立方体分割」法では、まず三次元データを分析
することにより、関心のある表面にまたがる隣接した頂
点データ点で構成される体積要素（ボクセル）を識別す
る。関心のある表面はデータ点によつて測定された物理
的性質について定義される。選択されたボクセルの各頂
点の空間座標が各頂点に対して計算された「勾配法線」
とともに「ディスプレイリスト（display list）」に集
められる。勾配法線はその頂点に於ける物理的性質の勾
配の方向を示し、与えられたボクセルの外側にあるが与
えられた頂点に隣接している三つのデータ点の値を調べ
ることにより決定される。立方体分割法は頂点相互の間
に内挿を行うことにより結果の画像の見掛けの分解能を
増大する方法も教示している。
のシェージングを決定する改善された方法が、米国特許
第４，７１９，５８５号に説明されている。そこに開示
された「立方体分割」法では、まず三次元データを分析
することにより、関心のある表面にまたがる隣接した頂
点データ点で構成される体積要素（ボクセル）を識別す
る。関心のある表面はデータ点によつて測定された物理
的性質について定義される。選択されたボクセルの各頂
点の空間座標が各頂点に対して計算された「勾配法線」
とともに「ディスプレイリスト（display list）」に集
められる。勾配法線はその頂点に於ける物理的性質の勾
配の方向を示し、与えられたボクセルの外側にあるが与
えられた頂点に隣接している三つのデータ点の値を調べ
ることにより決定される。立方体分割法は頂点相互の間
に内挿を行うことにより結果の画像の見掛けの分解能を
増大する方法も教示している。
【０００８】立方体分割法は通常の射線追跡に比べて二
つの利点がある。第一に、立方体分割法では画像の操作
と発生をより高速で行うことができる。ディスプレイリ
ストは関心のある構造の表面を記述するために必要なデ
ータだけを保持するので、回転やスケーリングのような
動作が早くなる。また、種々の点で構造の表面法線を近
似するために勾配法線を使うことは演算効率が良く、シ
ェージング情報をより早く発生することができる。
つの利点がある。第一に、立方体分割法では画像の操作
と発生をより高速で行うことができる。ディスプレイリ
ストは関心のある構造の表面を記述するために必要なデ
ータだけを保持するので、回転やスケーリングのような
動作が早くなる。また、種々の点で構造の表面法線を近
似するために勾配法線を使うことは演算効率が良く、シ
ェージング情報をより早く発生することができる。
【０００９】第二に、立方体分割法で作成されるディス
プレイリストの中の点の順番ともとのデータアレーの中
のそれらの点の座標との間に関係を必要としない。した
がって、骨と組織の異なる表面を表す数個のディスプレ
イリストの組み合わせは単に一つのディスプレイリスト
を他方のディスプレイリストの端に連結することにより
行うことができ、どちらのディスプレイリストも分類し
たり整理し直す必要はない。これは数個のディスプレイ
リストから複合画像を発生する際にかなり有利なことで
ある。
プレイリストの中の点の順番ともとのデータアレーの中
のそれらの点の座標との間に関係を必要としない。した
がって、骨と組織の異なる表面を表す数個のディスプレ
イリストの組み合わせは単に一つのディスプレイリスト
を他方のディスプレイリストの端に連結することにより
行うことができ、どちらのディスプレイリストも分類し
たり整理し直す必要はない。これは数個のディスプレイ
リストから複合画像を発生する際にかなり有利なことで
ある。
【００１０】しかし、立方体分割法で作成されるような
ディスプレイリストでは、特定の体積の中のデータを識
別することは多少難しい。これはデータがディスプレイ
リスト全体にわたって任意に分布されるからである。し
たがって、例えば予め定義された体積の中のデータ点を
識別して除去することによる画像の「カット」（切断）
にはディスプレイリストの中の各データ点の座標をその
体積を定義する一組の範囲と照合することが含まれる。
三次元の点に対しては、データ点毎に六つの試験を行わ
なければならない。すなわち、上の範囲と下の範囲に対
して一つの試験が行われ、これが三つのデカルト（直
交）座標の各々に対して繰り返される。このプロセスは
煩わしく、時間のかかるものである。
ディスプレイリストでは、特定の体積の中のデータを識
別することは多少難しい。これはデータがディスプレイ
リスト全体にわたって任意に分布されるからである。し
たがって、例えば予め定義された体積の中のデータ点を
識別して除去することによる画像の「カット」（切断）
にはディスプレイリストの中の各データ点の座標をその
体積を定義する一組の範囲と照合することが含まれる。
三次元の点に対しては、データ点毎に六つの試験を行わ
なければならない。すなわち、上の範囲と下の範囲に対
して一つの試験が行われ、これが三つのデカルト（直
交）座標の各々に対して繰り返される。このプロセスは
煩わしく、時間のかかるものである。
【００１１】予め定義された体積を記述する範囲の定義
は一般に直観的でなく、複雑な計算と多大の試行錯誤を
必要とすることがある。特に、予め定義された体積がも
とのデータアレーのデカルト座標系の軸に対して斜めの
軸に沿っている場合にそうである。それにもかかわら
ず、あるディスプレイリストの中のデータ点の空間位置
に基いてそのディスプレイリストの中のデータに対応す
る画像の部分を識別して分離できることはイメージング
される表面の更に詳細な理解を可能とする上で有益であ
る。例えば、三次元画像を切断して描写された表面の部
分を除去できれば、外科医は全く非侵襲的な診断方法を
使って詳細で複雑な手術手順を計画することができる。
三次元表面の画像を「切断する」ことにより、臓器や大
血管を避けるような適当な軌道を決める手術のシミュレ
ーションを行うことができる。
は一般に直観的でなく、複雑な計算と多大の試行錯誤を
必要とすることがある。特に、予め定義された体積がも
とのデータアレーのデカルト座標系の軸に対して斜めの
軸に沿っている場合にそうである。それにもかかわら
ず、あるディスプレイリストの中のデータ点の空間位置
に基いてそのディスプレイリストの中のデータに対応す
る画像の部分を識別して分離できることはイメージング
される表面の更に詳細な理解を可能とする上で有益であ
る。例えば、三次元画像を切断して描写された表面の部
分を除去できれば、外科医は全く非侵襲的な診断方法を
使って詳細で複雑な手術手順を計画することができる。
三次元表面の画像を「切断する」ことにより、臓器や大
血管を避けるような適当な軌道を決める手術のシミュレ
ーションを行うことができる。
【００１２】しかし、このようなシミュレーションされ
た手術が広く受け入れられるためには、画像のほぼ「実
時間」の操作が行えるように、対応するイメージング装
置の操作が容易で充分高速でなければならない。
た手術が広く受け入れられるためには、画像のほぼ「実
時間」の操作が行えるように、対応するイメージング装
置の操作が容易で充分高速でなければならない。
【００１３】
【発明の概要】本発明は立方体分割法や類似の画像発生
方法によって作成されるような三次元の点のディスプレ
イリストで記述される表面を切断する効率的な方法を提
供する。本発明の方法はディスプレイリストに直接作用
するものであり、それからディスプレイリストを編集し
た原始データアレーを参照する必要はない。更に、本発
明の方法は容易に通常の画像処理手順に融け込ませるこ
とができ、また高速の専用回路で具体化することができ
る。本発明はまた切断すべき体積を記述するための簡単
で直観的な方法を提供する。
方法によって作成されるような三次元の点のディスプレ
イリストで記述される表面を切断する効率的な方法を提
供する。本発明の方法はディスプレイリストに直接作用
するものであり、それからディスプレイリストを編集し
た原始データアレーを参照する必要はない。更に、本発
明の方法は容易に通常の画像処理手順に融け込ませるこ
とができ、また高速の専用回路で具体化することができ
る。本発明はまた切断すべき体積を記述するための簡単
で直観的な方法を提供する。
【００１４】とりわけ、写像装置は点のディスプレイリ
ストを受けて、各点に対する画像平面座標を発生する。
これらの画像平面座標を使って、切断すべき体積の画像
平面の中の輪郭を定める多角形を記憶するカットバッフ
ァ（cutbuffer）のアドレッシングを行う。多角形の内
側の面積はカットバッファの中で、記憶された第一の値
により指示される。
ストを受けて、各点に対する画像平面座標を発生する。
これらの画像平面座標を使って、切断すべき体積の画像
平面の中の輪郭を定める多角形を記憶するカットバッフ
ァ（cutbuffer）のアドレッシングを行う。多角形の内
側の面積はカットバッファの中で、記憶された第一の値
により指示される。
【００１５】ディスプレイリストを画像平面に写像する
間、画像平面座標を使ってカットバッファのアドレッシ
ングを行う。第一の値がカットバッファの指示されたア
ドレスに記憶されていれば、写像される点がカット点の
新しいディスプレイリストに記憶される。本発明の一つ
の目的はディスプレイリストの点をカット点の第二のデ
ィスプレイリストに分離する簡単な方法を提供すること
である。本発明が使用する画像平面座標は予め存在する
画像ディスプレイルーチンから得られる。そしてカット
バッファをアドレッシングすることにより各点が「切
断」されるか否かをチェックする付加的なステップはこ
のプロセスの時間を殆ど増加させない。カットバッファ
に記憶された二次元多角形はカット体積の二次元の簡単
な定義を与え、これはトラックボールのようなカーソル
制御装置により対話的に容易に入力することができる。
間、画像平面座標を使ってカットバッファのアドレッシ
ングを行う。第一の値がカットバッファの指示されたア
ドレスに記憶されていれば、写像される点がカット点の
新しいディスプレイリストに記憶される。本発明の一つ
の目的はディスプレイリストの点をカット点の第二のデ
ィスプレイリストに分離する簡単な方法を提供すること
である。本発明が使用する画像平面座標は予め存在する
画像ディスプレイルーチンから得られる。そしてカット
バッファをアドレッシングすることにより各点が「切
断」されるか否かをチェックする付加的なステップはこ
のプロセスの時間を殆ど増加させない。カットバッファ
に記憶された二次元多角形はカット体積の二次元の簡単
な定義を与え、これはトラックボールのようなカーソル
制御装置により対話的に容易に入力することができる。
【００１６】カット体積の第三次元を定義する、カット
軸に平行な点が画像平面内の単一の画素に写像されるよ
うにディスプレイリストの点の座標を回転させることに
より、まずディスプレイリストを変換することができ
る。したがって、カット体積の第三次元は画像平面に垂
直になっているものと考えることができる。本発明のも
う一つの目的は二次元ディスプレイを使って三次元カッ
ト体積を定義する簡単で直観的な方法を提供することで
ある。ディスプレイリストの回転により、カット体積の
第三の軸を二次元ディスプレイの平面に対して垂直にな
るように定義することができる。画像平面に対して平行
な多角形が他の二つの次元を定義する。多角形を画像に
重畳した状態で、回転した画像を見ることによって、切
断すべき体積が容易に理解される。カット軸に沿った切
断の深さは予め定義されたカット平面の片側の点を消去
することによって制御することができる。
軸に平行な点が画像平面内の単一の画素に写像されるよ
うにディスプレイリストの点の座標を回転させることに
より、まずディスプレイリストを変換することができ
る。したがって、カット体積の第三次元は画像平面に垂
直になっているものと考えることができる。本発明のも
う一つの目的は二次元ディスプレイを使って三次元カッ
ト体積を定義する簡単で直観的な方法を提供することで
ある。ディスプレイリストの回転により、カット体積の
第三の軸を二次元ディスプレイの平面に対して垂直にな
るように定義することができる。画像平面に対して平行
な多角形が他の二つの次元を定義する。多角形を画像に
重畳した状態で、回転した画像を見ることによって、切
断すべき体積が容易に理解される。カット軸に沿った切
断の深さは予め定義されたカット平面の片側の点を消去
することによって制御することができる。
【００１７】本発明の上記および他の目的および利点は
以下の説明から明らかとなる。説明には付図を参照する
が、付図は本明細書の一部を構成し、本発明の一実施例
を図示したものである。しかし、このような実施例は必
ずしも本発明の全範囲を表すものではないので、本発明
の範囲の解釈に当たっては特許請求の範囲を参照すべき
である。
以下の説明から明らかとなる。説明には付図を参照する
が、付図は本明細書の一部を構成し、本発明の一実施例
を図示したものである。しかし、このような実施例は必
ずしも本発明の全範囲を表すものではないので、本発明
の範囲の解釈に当たっては特許請求の範囲を参照すべき
である。
【００１８】
【詳しい説明】データ値の三次元アレーの取得
図１に示すように、ＣＴ装置またはＭＲ装置からのデー
タ値１０の三次元アレーがボクセル１２の各頂点に一つ
ずつデータ値を持つほぼ長方形の横断面の一組のボクセ
ル１２を定義する。明快にするため、ボクセル１２のデ
ータ値１０が黒丸で示され、一つにだけ参照番号を付し
てある。各データ値１０はカーテシアン座標系の三つの
寸法を与えられた、定義された空間ロケーションをそな
え、またその空間位置でイメージング対象の測定された
性質を表す値をそなえている。
タ値１０の三次元アレーがボクセル１２の各頂点に一つ
ずつデータ値を持つほぼ長方形の横断面の一組のボクセ
ル１２を定義する。明快にするため、ボクセル１２のデ
ータ値１０が黒丸で示され、一つにだけ参照番号を付し
てある。各データ値１０はカーテシアン座標系の三つの
寸法を与えられた、定義された空間ロケーションをそな
え、またその空間位置でイメージング対象の測定された
性質を表す値をそなえている。
【００１９】また図２に示すように、データ値１０の位
置での測定された性質の勾配を示すボクセル１２の各デ
ータ値１０に対してスケーリングされた勾配法線１４
（矢印で示してある）が定義される。データ値１０に対
応するこのスケーリングされた勾配法線１４はデータ値
１０に隣接しているがボクセル１２の外側にあり、デカ
ルト座標系のそれぞれの軸に沿ってデータ値１０からず
れたデータ値１１（ｘ）、１１（ｙ）、および１１
（ｚ）（図１で白丸で示してある）から決定される。デ
ータ値１０の値と１１（ｘ）の値との差はデータ値１０
での勾配法線のｘ成分である。同様に、データ値１０の
値と１１（ｙ）の値との差はデータ値１０での勾配法線
のｙ成分である。最後に、データ値１０の値と１１
（ｚ）の値との差はデータ値１０での勾配法線のｚ成分
である。これらの三つの差から形成される合成のベクト
ルが単位の大きさにスケーリングされることにより、以
後単に勾配法線と呼ぶスケーリングされた勾配法線１４
が形成される。
置での測定された性質の勾配を示すボクセル１２の各デ
ータ値１０に対してスケーリングされた勾配法線１４
（矢印で示してある）が定義される。データ値１０に対
応するこのスケーリングされた勾配法線１４はデータ値
１０に隣接しているがボクセル１２の外側にあり、デカ
ルト座標系のそれぞれの軸に沿ってデータ値１０からず
れたデータ値１１（ｘ）、１１（ｙ）、および１１
（ｚ）（図１で白丸で示してある）から決定される。デ
ータ値１０の値と１１（ｘ）の値との差はデータ値１０
での勾配法線のｘ成分である。同様に、データ値１０の
値と１１（ｙ）の値との差はデータ値１０での勾配法線
のｙ成分である。最後に、データ値１０の値と１１
（ｚ）の値との差はデータ値１０での勾配法線のｚ成分
である。これらの三つの差から形成される合成のベクト
ルが単位の大きさにスケーリングされることにより、以
後単に勾配法線と呼ぶスケーリングされた勾配法線１４
が形成される。
【００２０】ディスプレイリストの作成
図２に示すように、表面１６はイメージング対象の測定
された性質の閾値によって定義される。ボクセル１２の
類別は、１）各ボクセル１２のすべてのデータ値１０の
値がこの閾値より大きいか、またはすべてのデータ値１
０がこの閾値より小さいか、もしくは２）いくつかのデ
ータ値１０だけがこの閾値より大きくていくつかのデー
タ値１０だけがこの閾値より小さいか、に従って行われ
る。単一のボクセル１２の規模で、イメージング対象の
測定された性質の変化が空間で実質的に単調であると仮
定すると、この類別はボクセル１２が完全に表面１６の
片側にあるか、またはボクセル１２が表面１６にまたが
っているか判定することと同等である。
された性質の閾値によって定義される。ボクセル１２の
類別は、１）各ボクセル１２のすべてのデータ値１０の
値がこの閾値より大きいか、またはすべてのデータ値１
０がこの閾値より小さいか、もしくは２）いくつかのデ
ータ値１０だけがこの閾値より大きくていくつかのデー
タ値１０だけがこの閾値より小さいか、に従って行われ
る。単一のボクセル１２の規模で、イメージング対象の
測定された性質の変化が空間で実質的に単調であると仮
定すると、この類別はボクセル１２が完全に表面１６の
片側にあるか、またはボクセル１２が表面１６にまたが
っているか判定することと同等である。
【００２１】この閾値定義表面１６はユーザが選択す
る。本発明の医療への応用では、種々の身体構造を選択
的に見るように閾値を選択することができる。例えば、
骨の構造を見たい場合には、特定の閾値が選択される。
しかし、より軟部の組織の表面構造を見たい場合には、
ユーザが異なる閾値を選択する。表面１６にまたがるボ
クセル１２のデータ値はそれらのデータ値に対応する勾
配法線１４とともにディスプレイリスト１８（図３に示
す）に集められる。
る。本発明の医療への応用では、種々の身体構造を選択
的に見るように閾値を選択することができる。例えば、
骨の構造を見たい場合には、特定の閾値が選択される。
しかし、より軟部の組織の表面構造を見たい場合には、
ユーザが異なる閾値を選択する。表面１６にまたがるボ
クセル１２のデータ値はそれらのデータ値に対応する勾
配法線１４とともにディスプレイリスト１８（図３に示
す）に集められる。
【００２２】三次元表面を表す上記の方法は上記米国特
許に述べられている前記「立方体分割法」に教示されて
いる。図４に示されているように、ディスプレイリスト
１８がデータ値１０に対応する点ベクトル２０のリスト
としてコンピュータメモリに記憶される。各点ベクトル
２０は順次リストアドレスｉをそなえている。ディスプ
レイリストの中の点ベクトル２０の各アドレスｉには、
表面１６にまたがるボクセル１２のデータ値１０の空間
座標（ｘi ，ｙi ，ｚi ）とともにスケーリングされた
勾配法線のカーテシアン成分（ｎxi，ｎyi，ｎzi）によ
って表されるようなそのデータ値１０に対応する勾配法
線が含まれている。
許に述べられている前記「立方体分割法」に教示されて
いる。図４に示されているように、ディスプレイリスト
１８がデータ値１０に対応する点ベクトル２０のリスト
としてコンピュータメモリに記憶される。各点ベクトル
２０は順次リストアドレスｉをそなえている。ディスプ
レイリストの中の点ベクトル２０の各アドレスｉには、
表面１６にまたがるボクセル１２のデータ値１０の空間
座標（ｘi ，ｙi ，ｚi ）とともにスケーリングされた
勾配法線のカーテシアン成分（ｎxi，ｎyi，ｎzi）によ
って表されるようなそのデータ値１０に対応する勾配法
線が含まれている。
【００２３】したがって、ディスプレイリスト１８には
表面１６にまたがるボクセル１２上のデータ値１０の座
標およびスケーリングされた勾配法線１４だけが含まれ
ている。このときディスプレイリスト１８がＣＴまたは
ＭＲイメージングシーケンスで収集されたデータの限定
された部分集合となる。イメージング対象の関連構造の
表面１６だけが記述されており、完全に表面１６の内側
または外側にあるボクセル１２からのそれらのデータ値
１０は記述されない。
表面１６にまたがるボクセル１２上のデータ値１０の座
標およびスケーリングされた勾配法線１４だけが含まれ
ている。このときディスプレイリスト１８がＣＴまたは
ＭＲイメージングシーケンスで収集されたデータの限定
された部分集合となる。イメージング対象の関連構造の
表面１６だけが記述されており、完全に表面１６の内側
または外側にあるボクセル１２からのそれらのデータ値
１０は記述されない。
【００２４】画像平面へのディスプレイリストの写像
図３に示すように、画像２２はディスプレイリスト１８
から、ディスプレイリスト１８で表される三次元表面１
６の二次元画像平面２４への数学的投影すなわち写像に
よって作られる。画像平面２４は画素２６の長方形アレ
ーで構成される。各画素は二次元画像平面座標（ｘ’，
ｙ’）をそなえている。この二次元画像平面座標はディ
スプレイリスト１８の点ベクトル２０（図４）の三次元
空間座標（ｘ，ｙ，ｚ）から最初に識別すべきものであ
る。
から、ディスプレイリスト１８で表される三次元表面１
６の二次元画像平面２４への数学的投影すなわち写像に
よって作られる。画像平面２４は画素２６の長方形アレ
ーで構成される。各画素は二次元画像平面座標（ｘ’，
ｙ’）をそなえている。この二次元画像平面座標はディ
スプレイリスト１８の点ベクトル２０（図４）の三次元
空間座標（ｘ，ｙ，ｚ）から最初に識別すべきものであ
る。
【００２５】写像プロセスでは、表面１６が画像２２と
して観測される見掛けの方向を示すように図３に示す観
測軸２８を選択しなければならない。便宜上、この観測
軸２８はディスプレイリスト１８のｚ軸に平行とされ
る。ディスプレイリスト１８の中の各点ベクトル２０
（図４）の空間座標および勾配法線成分を後で簡単に述
べるように適当に変換することにより画像２２の回転が
行われる。観測軸２８（図３）をこのように固定するこ
とにより、ディスプレイリスト１８の点ベクトル２０
（図４）のｘ座標およびｙ座標を画像平面２４（図３）
のｘ’座標およびｙ’座標に直接、写像することができ
る。画像平面２４のスケーリングおよび平行移動は再
び、ディスプレイリスト１８の点ベクトル２０（図４）
を変換することによって行われる。
して観測される見掛けの方向を示すように図３に示す観
測軸２８を選択しなければならない。便宜上、この観測
軸２８はディスプレイリスト１８のｚ軸に平行とされ
る。ディスプレイリスト１８の中の各点ベクトル２０
（図４）の空間座標および勾配法線成分を後で簡単に述
べるように適当に変換することにより画像２２の回転が
行われる。観測軸２８（図３）をこのように固定するこ
とにより、ディスプレイリスト１８の点ベクトル２０
（図４）のｘ座標およびｙ座標を画像平面２４（図３）
のｘ’座標およびｙ’座標に直接、写像することができ
る。画像平面２４のスケーリングおよび平行移動は再
び、ディスプレイリスト１８の点ベクトル２０（図４）
を変換することによって行われる。
【００２６】回転、平行移動、およびスケーリングのこ
れらの変換プロセスの各々では、点ベクトル２０の各ア
ドレスの空間座標および法線成分に４×４の変換行列を
乗算することにより座標および法線値を適当に変換した
新しいベクトルを作成する。これらの変換技術は当業者
には周知のものであり、ジェームス・ディー・フォリー
（Ｊａｍｅｓ Ｄ． Ｆｏｌｅｙ）他著のテキスト、
「対話型コンピュータグラフィックスの基礎」（Ｆｕｎ
ｄａｍｅｎｔａｌｓ ｏｆ Ｉｎｔｅｒａｃｔｉｖｅ
Ｃｏｍｐｕｔｅｒ Ｇｒａｐｈｉｃｓ）（１９８４）に
詳細に述べられている。
れらの変換プロセスの各々では、点ベクトル２０の各ア
ドレスの空間座標および法線成分に４×４の変換行列を
乗算することにより座標および法線値を適当に変換した
新しいベクトルを作成する。これらの変換技術は当業者
には周知のものであり、ジェームス・ディー・フォリー
（Ｊａｍｅｓ Ｄ． Ｆｏｌｅｙ）他著のテキスト、
「対話型コンピュータグラフィックスの基礎」（Ｆｕｎ
ｄａｍｅｎｔａｌｓ ｏｆ Ｉｎｔｅｒａｃｔｉｖｅ
Ｃｏｍｐｕｔｅｒ Ｇｒａｐｈｉｃｓ）（１９８４）に
詳細に述べられている。
【００２７】ディスプレイリストの「クリッピング」も
カット平面をその空間座標（ｘ，ｙ，ｚ）について設定
し、カット平面の片側に入るすべての点ベクトル２０を
捨てることにより行うことができる。カット平面は一般
に次式で記述することができる。 Ａｘ＋Ｂｙ＋Ｃｚ＋Ｄ＝０ （１） 但し、Ａ、Ｂ、Ｃ、およびＤはユーザが設定する定数で
あり、ｘ、ｙ、およびｚは点ベクトル２０の空間座標で
ある。クリッピングは各点ベクトル２０の空間座標の値
を式（１）に代入し、例えば式の右辺が０より大きい点
ベクトルを捨てることにより行われる。
カット平面をその空間座標（ｘ，ｙ，ｚ）について設定
し、カット平面の片側に入るすべての点ベクトル２０を
捨てることにより行うことができる。カット平面は一般
に次式で記述することができる。 Ａｘ＋Ｂｙ＋Ｃｚ＋Ｄ＝０ （１） 但し、Ａ、Ｂ、Ｃ、およびＤはユーザが設定する定数で
あり、ｘ、ｙ、およびｚは点ベクトル２０の空間座標で
ある。クリッピングは各点ベクトル２０の空間座標の値
を式（１）に代入し、例えば式の右辺が０より大きい点
ベクトルを捨てることにより行われる。
【００２８】図９に示すように、この変換プロセスはデ
ィジタルコンピュータ上で動作する回転／平行移動／ス
ケール（ＲＴＳ−ｒｏｔａｔｅ／ｔｒａｎｓｌａｔｅ／
ｓｃａｌｅ）ソフトウェアルーチンにより、またはこれ
らの変換に必要な計算専用の特殊ハードウェアによって
遂行することができる。この両方とも一般的にＲＴＳブ
ロック３０で表されている。図９の装置は通常程度の当
業者には明らかなように完全にソフトウェアにより実現
してもよく、あるいは専用回路で具体化してもよい。
ィジタルコンピュータ上で動作する回転／平行移動／ス
ケール（ＲＴＳ−ｒｏｔａｔｅ／ｔｒａｎｓｌａｔｅ／
ｓｃａｌｅ）ソフトウェアルーチンにより、またはこれ
らの変換に必要な計算専用の特殊ハードウェアによって
遂行することができる。この両方とも一般的にＲＴＳブ
ロック３０で表されている。図９の装置は通常程度の当
業者には明らかなように完全にソフトウェアにより実現
してもよく、あるいは専用回路で具体化してもよい。
【００２９】通常、変換プロセスに固有の近似誤差の結
果として、ＲＴＳブロック３０により第二の変換された
ディスプレイリスト１８’が作成され、原始ディスプレ
イリスト１８が保存される。後続の変換は原始ディスプ
レイリスト１８を使用する。ディスプレイリスト１８お
よび１８’のアドレスｉは変換プロセスによって変化し
ない。したがって、どの変換でも、各リスト１８と１
８’の間の対応する点ベクトル２０（図４）は容易に識
別することができる。
果として、ＲＴＳブロック３０により第二の変換された
ディスプレイリスト１８’が作成され、原始ディスプレ
イリスト１８が保存される。後続の変換は原始ディスプ
レイリスト１８を使用する。ディスプレイリスト１８お
よび１８’のアドレスｉは変換プロセスによって変化し
ない。したがって、どの変換でも、各リスト１８と１
８’の間の対応する点ベクトル２０（図４）は容易に識
別することができる。
【００３０】図８および図９に示すように、変換された
ディスプレイリスト１８’を画像平面２４に写像するプ
ロセスはプロセスブロック３２により変換後のディスプ
レイリスト１８’から第ｉ番目の点ベクトル２０を読み
出すことにより開始される。この第ｉ番目の点ベクトル
２０のｘおよびｙの空間座標ｘi およびｙi を使用して
ｚバッファ３４内のメモリロケーションのアドレッシン
グを行う。
ディスプレイリスト１８’を画像平面２４に写像するプ
ロセスはプロセスブロック３２により変換後のディスプ
レイリスト１８’から第ｉ番目の点ベクトル２０を読み
出すことにより開始される。この第ｉ番目の点ベクトル
２０のｘおよびｙの空間座標ｘi およびｙi を使用して
ｚバッファ３４内のメモリロケーションのアドレッシン
グを行う。
【００３１】図６に示すように、ｚバッファ３４はメモ
リロケーションの二次元アレーである。各メモリロケー
ションアドレスは座標ｘ’，ｙ’を持つ画像平面２４の
特定の画素２６に同定される。写像プロセスでは、第ｉ
番目のベクトル２０上の座標ｘi ，ｙi が画像平面座標
ｘ’，ｙ’となる。ｚバッファ３４の中の各メモリロケ
ーションは最初に、変換されたディスプレイリスト１
８’（図９）の中の点ベクトル２０の最大期待値ｚi よ
り大きい任意に大きい値で充たされる。
リロケーションの二次元アレーである。各メモリロケー
ションアドレスは座標ｘ’，ｙ’を持つ画像平面２４の
特定の画素２６に同定される。写像プロセスでは、第ｉ
番目のベクトル２０上の座標ｘi ，ｙi が画像平面座標
ｘ’，ｙ’となる。ｚバッファ３４の中の各メモリロケ
ーションは最初に、変換されたディスプレイリスト１
８’（図９）の中の点ベクトル２０の最大期待値ｚi よ
り大きい任意に大きい値で充たされる。
【００３２】再び図８および図９を参照して説明する。
図８のプロセスブロック３６に示すように、選択された
第ｉ番目の点ベクトル２０のｚi の値がｚバッファ３４
の中のアドレッシングされたメモリロケーションの中の
値ｚバッファ（ｘi ，ｙi ）と比較される。これは図９
の対応する機能ブロック３７によって行われる。最初の
点ベクトル２０（図４）の場合、ｚバッファ３４が大き
な値に初期設定された結果として値ｚi がｚバッファ
（ｘi ，ｙi ）より小さくなる。この第一の場合、ｚi
すなわち第ｉ番目の点ベクトル２０のｚ空間座標がその
メモリロケーションでｚバッファ３４の値を書き直す。
これにより、図８のプロセスブロック３８に示すように
ｚバッファ（ｘi ，ｙi ）がｚi に等しく設定される。
これは図９のゲート３９によって行われる。
図８のプロセスブロック３６に示すように、選択された
第ｉ番目の点ベクトル２０のｚi の値がｚバッファ３４
の中のアドレッシングされたメモリロケーションの中の
値ｚバッファ（ｘi ，ｙi ）と比較される。これは図９
の対応する機能ブロック３７によって行われる。最初の
点ベクトル２０（図４）の場合、ｚバッファ３４が大き
な値に初期設定された結果として値ｚi がｚバッファ
（ｘi ，ｙi ）より小さくなる。この第一の場合、ｚi
すなわち第ｉ番目の点ベクトル２０のｚ空間座標がその
メモリロケーションでｚバッファ３４の値を書き直す。
これにより、図８のプロセスブロック３８に示すように
ｚバッファ（ｘi ，ｙi ）がｚi に等しく設定される。
これは図９のゲート３９によって行われる。
【００３３】また第一の場合には図８のプロセスブロッ
ク４０および図９の機能ブロック４１で示すように、第
ｉ番目の点ベクトル２０（図４）に対応する勾配法線１
４（図２）を使って画像平面２４の対応する画素２６に
対してｃｏｌｏｒ（ｉ）を計算する。第ｉ番目の点ベク
トル２０の勾配法線１４からｃｏｌｏｒ（ｉ）を計算す
るこのプロセスについては後述する。
ク４０および図９の機能ブロック４１で示すように、第
ｉ番目の点ベクトル２０（図４）に対応する勾配法線１
４（図２）を使って画像平面２４の対応する画素２６に
対してｃｏｌｏｒ（ｉ）を計算する。第ｉ番目の点ベク
トル２０の勾配法線１４からｃｏｌｏｒ（ｉ）を計算す
るこのプロセスについては後述する。
【００３４】このように計算されたｃｏｌｏｒ（ｉ）は
画像バッファ４２に記憶される。図５に示すように、画
像バッファ４２はｚバッファ３４と同じ形式であり、メ
モリロケーションの二次元アレーである。各々のアドレ
スは画像平面２４の画素２６に対応し、その画素２６の
カラーまたはシェージングを保持する。再び図８および
図９を参照して説明する。第一の場合には、第ｉ番目の
点ベクトルに対する計算されたｃｏｌｏｒ（ｉ）が第ｉ
番目の点ベクトル２０（図４）から決定された画像バッ
ファ４２のアドレスｘi 、ｙi で画像バッファ４２に記
憶される。画像バッファ４２の中の予め存在する任意の
値が書き直され、このようにアドレッシングされた画像
バッファの値すなわち画像バッファ（ｘi ，ｙi）がｃ
ｏｌｏｒ（ｉ）に等しくなるように設定される。このよ
うに画像バッファ４２にｃｏｌｏｒ（ｉ）を記憶するこ
とは図８のプロセスブロック４４で表され、図９のゲー
ト４５により行われる。
画像バッファ４２に記憶される。図５に示すように、画
像バッファ４２はｚバッファ３４と同じ形式であり、メ
モリロケーションの二次元アレーである。各々のアドレ
スは画像平面２４の画素２６に対応し、その画素２６の
カラーまたはシェージングを保持する。再び図８および
図９を参照して説明する。第一の場合には、第ｉ番目の
点ベクトルに対する計算されたｃｏｌｏｒ（ｉ）が第ｉ
番目の点ベクトル２０（図４）から決定された画像バッ
ファ４２のアドレスｘi 、ｙi で画像バッファ４２に記
憶される。画像バッファ４２の中の予め存在する任意の
値が書き直され、このようにアドレッシングされた画像
バッファの値すなわち画像バッファ（ｘi ，ｙi）がｃ
ｏｌｏｒ（ｉ）に等しくなるように設定される。このよ
うに画像バッファ４２にｃｏｌｏｒ（ｉ）を記憶するこ
とは図８のプロセスブロック４４で表され、図９のゲー
ト４５により行われる。
【００３５】第二の場合には、図８のプロセスブロック
３６で判定され、図９の機能ブロック３７で遂行される
ように第ｉ番目の点ベクトル２０（図４）に対応する値
ｚiがｚバッファ３４の値ｚバッファ（ｘi ，ｙi ）よ
り大きいとき、ｚバッファも画像バッファ４２も修正さ
れない。後で述べるように、また図８の判定ブロック４
６で判定されるように、ディスプレイリストを切断すべ
きでない場合には、ルーチンが判定ブロック４７にジャ
ンプして、すべての点ベクトル２０（図４）が読み出さ
れたか否かが判定される。すべての点ベクトル２０が読
み出された場合には、プログラムが終了する。すべての
点ベクトル２０が読み出されていない場合には、プログ
ラムがプロセスブロック３２に戻り、ディスプレイリス
ト１８の中の残っている点ベクトル２０に対してプロセ
スブロック３２−４４が反復される。
３６で判定され、図９の機能ブロック３７で遂行される
ように第ｉ番目の点ベクトル２０（図４）に対応する値
ｚiがｚバッファ３４の値ｚバッファ（ｘi ，ｙi ）よ
り大きいとき、ｚバッファも画像バッファ４２も修正さ
れない。後で述べるように、また図８の判定ブロック４
６で判定されるように、ディスプレイリストを切断すべ
きでない場合には、ルーチンが判定ブロック４７にジャ
ンプして、すべての点ベクトル２０（図４）が読み出さ
れたか否かが判定される。すべての点ベクトル２０が読
み出された場合には、プログラムが終了する。すべての
点ベクトル２０が読み出されていない場合には、プログ
ラムがプロセスブロック３２に戻り、ディスプレイリス
ト１８の中の残っている点ベクトル２０に対してプロセ
スブロック３２−４４が反復される。
【００３６】ディスプレイリスト１８の中の各点ベクト
ル２０（図４）を調べた後、画像バッファ４２にはｘi
およびｙi 空間座標が同じで、他の任意の点ベクトル２
０より低い値ｚi をそなえた点ベクトル２０のカラーだ
けが入っている。したがって、観測軸２８（図３）で決
まるような観測者に最も近い値ｚi を持つ点ベクトル２
０だけがディスプレイされる。観測者はより高い空間ｚ
値の方を見て低い空間ｚ値に居るものと仮定する。それ
らの同じｘi およびｙi 座標にある他の点ベクトル２０
は事実上、隠される。
ル２０（図４）を調べた後、画像バッファ４２にはｘi
およびｙi 空間座標が同じで、他の任意の点ベクトル２
０より低い値ｚi をそなえた点ベクトル２０のカラーだ
けが入っている。したがって、観測軸２８（図３）で決
まるような観測者に最も近い値ｚi を持つ点ベクトル２
０だけがディスプレイされる。観測者はより高い空間ｚ
値の方を見て低い空間ｚ値に居るものと仮定する。それ
らの同じｘi およびｙi 座標にある他の点ベクトル２０
は事実上、隠される。
【００３７】画像バッファ４２には最終的に画素平面２
４（図３）の表現が含まれ、これは当業者には知られて
いるように専用のハードウェアによりディスプレイ装置
４８に読み出される。画像２２の対応するディスプレイ
がディスプレイ装置に投影される。画素カラーの計算 図３に示すような画像２２としての表面１６のディスプ
レイでは、表面１６の各要素の向きの人間の目に視覚的
な手掛りを与えることが非常に重要である。これらの視
覚的な手掛りは画像平面２４の画素２６のカラーおよび
シェージングによって与えられる。勾配法線１４の角度
が観測軸２８の角度に近くなるにつれて、（負画像では
なくて少なくとも正画像に対して）印加されるシェージ
ングは明るくなる。観測軸２８からずれた方向の成分を
そなえた勾配法線１４（図２）を示す表面１６の要素は
見えない表面１６の構造を表す。したがって、これらの
勾配法線１４は特定の観測角度に対するビュー（ｖｉｅ
ｗ）からこのような点ベクトル２０（図４）を消去する
ための機構を提供する。相当な成分が観測軸２８に対し
て垂直な方向にある勾配法線１４をそなえた表面１６の
要素はより暗くシェージングされた画素２６によって表
される。
４（図３）の表現が含まれ、これは当業者には知られて
いるように専用のハードウェアによりディスプレイ装置
４８に読み出される。画像２２の対応するディスプレイ
がディスプレイ装置に投影される。画素カラーの計算 図３に示すような画像２２としての表面１６のディスプ
レイでは、表面１６の各要素の向きの人間の目に視覚的
な手掛りを与えることが非常に重要である。これらの視
覚的な手掛りは画像平面２４の画素２６のカラーおよび
シェージングによって与えられる。勾配法線１４の角度
が観測軸２８の角度に近くなるにつれて、（負画像では
なくて少なくとも正画像に対して）印加されるシェージ
ングは明るくなる。観測軸２８からずれた方向の成分を
そなえた勾配法線１４（図２）を示す表面１６の要素は
見えない表面１６の構造を表す。したがって、これらの
勾配法線１４は特定の観測角度に対するビュー（ｖｉｅ
ｗ）からこのような点ベクトル２０（図４）を消去する
ための機構を提供する。相当な成分が観測軸２８に対し
て垂直な方向にある勾配法線１４をそなえた表面１６の
要素はより暗くシェージングされた画素２６によって表
される。
【００３８】第ｉ番目の画素２６（図３）に対する画像
バッファ４２の中のｃｏｌｏｒ（ｉ）の値で表されるそ
の画素２６に対するシェージングは図８のプロセスブロ
ック４０で決定される。これは勾配法線１４（図２）と
ユーザが定義し得る照明軸５１（図３）のドット積を求
めることにより図９の機能ブロック４１により行われ
る。更に、真のカラー情報を探索テーブル（図示しな
い）により付加してもよい。この探索テーブルはデータ
値１０（図２）から明らかにされた組織の型に応じて各
画素２６にカラー（赤、緑、青）を割り当てる役目を果
たす。ここで用いるカラーという用語はディスプレイさ
れる画素２６のシェージングとスペクトルの両方の特性
を含む。
バッファ４２の中のｃｏｌｏｒ（ｉ）の値で表されるそ
の画素２６に対するシェージングは図８のプロセスブロ
ック４０で決定される。これは勾配法線１４（図２）と
ユーザが定義し得る照明軸５１（図３）のドット積を求
めることにより図９の機能ブロック４１により行われ
る。更に、真のカラー情報を探索テーブル（図示しな
い）により付加してもよい。この探索テーブルはデータ
値１０（図２）から明らかにされた組織の型に応じて各
画素２６にカラー（赤、緑、青）を割り当てる役目を果
たす。ここで用いるカラーという用語はディスプレイさ
れる画素２６のシェージングとスペクトルの両方の特性
を含む。
【００３９】カットバッファの作成
図１０に示すように、画像２２を切断する前に画像２２
上で多角形領域５０が追跡される。多角形領域５０の入
力は、当業者には一般に理解されるようにトラックボー
ル等の制御装置によりディスプレイされた画像２２にス
ーパインポーズされるカーソル（図示しない）を動かす
ことにより、また多角形領域５０の頂点を形成する種々
の点５４を示すことにより行われる。入力を容易にする
ため、多角形領域５０は一般に必要最少限の点５４で定
義される。しかし、多角形領域５０は限界内であれば任
意の連続した閉じた表面であってもよい。
上で多角形領域５０が追跡される。多角形領域５０の入
力は、当業者には一般に理解されるようにトラックボー
ル等の制御装置によりディスプレイされた画像２２にス
ーパインポーズされるカーソル（図示しない）を動かす
ことにより、また多角形領域５０の頂点を形成する種々
の点５４を示すことにより行われる。入力を容易にする
ため、多角形領域５０は一般に必要最少限の点５４で定
義される。しかし、多角形領域５０は限界内であれば任
意の連続した閉じた表面であってもよい。
【００４０】図７に示すように、図１０のこれらの点５
４の座標が図９に示すカットバッファ５６に入力され
る。カットバッファ５６はｚバッファ３４および画像バ
ッファ４２と同様の形式になっており、画像平面２４お
よびディスプレイされる画像２２の各画素２６（図３）
に対応するメモリロケーションをそなえている。カット
バッファ５６のメモリロケーションは”１”または”
０”のブール値を保持しており、最初に０に設定され
る。適当な画素座標のメモリロケーションの値を”１”
に設定することにより、多角形領域５０（図１０）の点
５４の座標がカットバッファに記憶される。
４の座標が図９に示すカットバッファ５６に入力され
る。カットバッファ５６はｚバッファ３４および画像バ
ッファ４２と同様の形式になっており、画像平面２４お
よびディスプレイされる画像２２の各画素２６（図３）
に対応するメモリロケーションをそなえている。カット
バッファ５６のメモリロケーションは”１”または”
０”のブール値を保持しており、最初に０に設定され
る。適当な画素座標のメモリロケーションの値を”１”
に設定することにより、多角形領域５０（図１０）の点
５４の座標がカットバッファに記憶される。
【００４１】一旦多角形領域５０のすべての頂点５４が
入力されれば、カットバッファ５６内で表される多角形
領域５０の内側の点も周知の多角形充てんアルゴリズム
により”１”に設定される。カットバッファ５６の中の
メモリロケーションの値が１に設定されるとき、それら
と同じ座標で画像バッファに記憶されるカラーが対照的
なカラー値に設定されて、多角形領域５０が明確に識別
できるようにする。
入力されれば、カットバッファ５６内で表される多角形
領域５０の内側の点も周知の多角形充てんアルゴリズム
により”１”に設定される。カットバッファ５６の中の
メモリロケーションの値が１に設定されるとき、それら
と同じ座標で画像バッファに記憶されるカラーが対照的
なカラー値に設定されて、多角形領域５０が明確に識別
できるようにする。
【００４２】ディスプレイリストの切断
再び図８および図９を参照して説明する。ディスプレイ
リスト１８を切断するプロセスは画像２２のディスプレ
イに使用されるのと同じ手順を使用する。これにより、
カットプロセスは前から存在するハードウェアおよびソ
フトウェアを有効に使うことができる。カットプロセス
では、ディスプレイリスト１８の中の各点ベクトル２０
（図３）がプロセスブロック３２により前と同様に読み
出される。ｚバッファ３４および画像バッファ４２は妥
当な場合、プロセスブロック３６、３８、４０、および
４４により修正される。しかし判定ブロック４６では、
プログラムは判定ブロック４７ではなく判定ブロック５
８に進む。
リスト１８を切断するプロセスは画像２２のディスプレ
イに使用されるのと同じ手順を使用する。これにより、
カットプロセスは前から存在するハードウェアおよびソ
フトウェアを有効に使うことができる。カットプロセス
では、ディスプレイリスト１８の中の各点ベクトル２０
（図３）がプロセスブロック３２により前と同様に読み
出される。ｚバッファ３４および画像バッファ４２は妥
当な場合、プロセスブロック３６、３８、４０、および
４４により修正される。しかし判定ブロック４６では、
プログラムは判定ブロック４７ではなく判定ブロック５
８に進む。
【００４３】判定ブロック５８では、第ｉ番目の点ベク
トル２０（図４）の空間座標ｘi およびｙi を使用して
カットバッファ５６のアドレッシングを行う。アドレッ
シングされたメモリロケーションでのカットバッファ５
６の値が１であれば、すなわち第ｉ番目の点ベクトル２
０が多角形領域５０の中のｘi およびｙi 座標をそなえ
ていることが示されれば、図８のプロセスブロック６４
および図９のゲート６５で示されるように原始ディスプ
レイリスト１８からアドレスｉで読み出された点ベクト
ル２０がカットディスプレイリスト６０の次の空きロケ
ーションｊすなわちカットディスプレイリスト（ｊ）に
記憶される。逆に、カットバッファ５６のそのアドレッ
シングされたメモリロケーションの値が０であれば、す
なわち第ｉ番目の点ベクトル２０が多角形領域５０の外
側のｘi およびｙi 座標をそなえていることが示されれ
ば、図８のプロセスブロック６６および図９のゲート６
７で示されるようにディスプレイリスト１８からアドレ
スｉで読み出された点ベクトル２０が残余ディスプレイ
リスト６２の次の空きロケーションｊ’すなわち残余デ
ィスプレイリスト（ｊ’）に記憶される。
トル２０（図４）の空間座標ｘi およびｙi を使用して
カットバッファ５６のアドレッシングを行う。アドレッ
シングされたメモリロケーションでのカットバッファ５
６の値が１であれば、すなわち第ｉ番目の点ベクトル２
０が多角形領域５０の中のｘi およびｙi 座標をそなえ
ていることが示されれば、図８のプロセスブロック６４
および図９のゲート６５で示されるように原始ディスプ
レイリスト１８からアドレスｉで読み出された点ベクト
ル２０がカットディスプレイリスト６０の次の空きロケ
ーションｊすなわちカットディスプレイリスト（ｊ）に
記憶される。逆に、カットバッファ５６のそのアドレッ
シングされたメモリロケーションの値が０であれば、す
なわち第ｉ番目の点ベクトル２０が多角形領域５０の外
側のｘi およびｙi 座標をそなえていることが示されれ
ば、図８のプロセスブロック６６および図９のゲート６
７で示されるようにディスプレイリスト１８からアドレ
スｉで読み出された点ベクトル２０が残余ディスプレイ
リスト６２の次の空きロケーションｊ’すなわち残余デ
ィスプレイリスト（ｊ’）に記憶される。
【００４４】変換されたディスプレイリスト１８’の点
ベクトル２０（図４）を使用してカットバッファ５６の
アドレッシングを行うが、原始ディスプレイリスト１８
からの点ベクトル２０はカットリスト６０および残余リ
スト６２に分離される。ディスプレイリスト１８に対す
るアドレスｉとディスプレイリスト１８’に対するアド
レスｉとの間の対応がこれを簡単なことにしている。カ
ットディスプレイリスト６０および残余ディスプレイリ
スト６２を作成するために原始ディスプレイリスト１８
を使うことにより、変換プロセスで生じたエラーからカ
ット点ベクトル２０が守られる。しかし、このようなエ
ラーが問題でなければ、変換されたディスプレイリスト
１８’からの点ベクトル２０を使ってカットディスプレ
イリスト６０および残余ディスプレイリスト６２を発生
してもよいことは明らかであろう。
ベクトル２０（図４）を使用してカットバッファ５６の
アドレッシングを行うが、原始ディスプレイリスト１８
からの点ベクトル２０はカットリスト６０および残余リ
スト６２に分離される。ディスプレイリスト１８に対す
るアドレスｉとディスプレイリスト１８’に対するアド
レスｉとの間の対応がこれを簡単なことにしている。カ
ットディスプレイリスト６０および残余ディスプレイリ
スト６２を作成するために原始ディスプレイリスト１８
を使うことにより、変換プロセスで生じたエラーからカ
ット点ベクトル２０が守られる。しかし、このようなエ
ラーが問題でなければ、変換されたディスプレイリスト
１８’からの点ベクトル２０を使ってカットディスプレ
イリスト６０および残余ディスプレイリスト６２を発生
してもよいことは明らかであろう。
【００４５】図１０に示すように、カットディスプレイ
リスト６０の中の点ベクトル２０（図４）は多角形領域
５０の中のディスプレイ装置４８で見える画素２６に対
応する点ベクトル２０だけではなくて、見える画素２６
と同じｘおよびｙの空間座標をそなえ、事実上多角形領
域５０の「下」にある「隠れた」点ベクトル２０でもあ
る。したがって、カット体積は多角形領域５０が変換さ
れたディスプレイリスト１８’のｚ軸に沿って動くとき
多角形領域５０によって掃かれる体積である。これによ
り、ディスプレイされた画像に対して直観的にカット体
積を定義することができ、またＲＴＳ機能ブロック３０
により画像を回転してカット軸を定義することができ
る。カット軸はユーザの希望するカットに応じて任意に
選択することができる。カット軸はそれが観測者から離
れる方向に伸びるように画像を回転させる対話型手順に
よって入力される。
リスト６０の中の点ベクトル２０（図４）は多角形領域
５０の中のディスプレイ装置４８で見える画素２６に対
応する点ベクトル２０だけではなくて、見える画素２６
と同じｘおよびｙの空間座標をそなえ、事実上多角形領
域５０の「下」にある「隠れた」点ベクトル２０でもあ
る。したがって、カット体積は多角形領域５０が変換さ
れたディスプレイリスト１８’のｚ軸に沿って動くとき
多角形領域５０によって掃かれる体積である。これによ
り、ディスプレイされた画像に対して直観的にカット体
積を定義することができ、またＲＴＳ機能ブロック３０
により画像を回転してカット軸を定義することができ
る。カット軸はユーザの希望するカットに応じて任意に
選択することができる。カット軸はそれが観測者から離
れる方向に伸びるように画像を回転させる対話型手順に
よって入力される。
【００４６】このようにして、図９に示すようにディス
プレイリスト１８が迅速にカットディスプレイリスト６
０および残余ディスプレイリスト６２に分離される。こ
れらのディスプレイリストは適当な記憶装置に退避させ
てもよい。通常、この記憶装置には随意に読み出せるフ
ロッピディスク、ハードディスク駆動装置、またはテー
プのようなある形式の磁気記憶媒体が含まれる。
プレイリスト１８が迅速にカットディスプレイリスト６
０および残余ディスプレイリスト６２に分離される。こ
れらのディスプレイリストは適当な記憶装置に退避させ
てもよい。通常、この記憶装置には随意に読み出せるフ
ロッピディスク、ハードディスク駆動装置、またはテー
プのようなある形式の磁気記憶媒体が含まれる。
【００４７】各ディスプレイリスト６０、６２は別々に
ディスプレイおよび操作してもよいし、あるいは組み合
わせ画像のため他のデータと組み合わせてもよい。図１
１では、ディスプレイされた画像２２で見えない隠れた
表面６８を明らかにするため、カットディスプレイリス
ト６０からの画像が９０°回転してディスプレイされて
いる。カット体積は多角形領域５０がディスプレイ画像
２２のｚ軸に沿って掃くとき多角形領域５０のへりによ
って形成される破線７０および図９に示すＲＴＳ機能ブ
ロック３０によりディスプレイリスト１８に予め印加さ
れたカット平面７２で表される。
ディスプレイおよび操作してもよいし、あるいは組み合
わせ画像のため他のデータと組み合わせてもよい。図１
１では、ディスプレイされた画像２２で見えない隠れた
表面６８を明らかにするため、カットディスプレイリス
ト６０からの画像が９０°回転してディスプレイされて
いる。カット体積は多角形領域５０がディスプレイ画像
２２のｚ軸に沿って掃くとき多角形領域５０のへりによ
って形成される破線７０および図９に示すＲＴＳ機能ブ
ロック３０によりディスプレイリスト１８に予め印加さ
れたカット平面７２で表される。
【００４８】点ベクトルのアドレスが任意なディスプレ
イリスト１８、１８’、６０、および６２のリストフォ
ーマットにより、数個のディスプレイリストからの複合
画像の作成は一つのディスプレイリストをもう一つのデ
ィスプレイリストの端に付加するだけの簡単なものとな
る。本発明の特定の利点は画像を多数の画像に切断する
などの画像の対話型操作を行うように表面画像データを
迅速に処理することにある。
イリスト１８、１８’、６０、および６２のリストフォ
ーマットにより、数個のディスプレイリストからの複合
画像の作成は一つのディスプレイリストをもう一つのデ
ィスプレイリストの端に付加するだけの簡単なものとな
る。本発明の特定の利点は画像を多数の画像に切断する
などの画像の対話型操作を行うように表面画像データを
迅速に処理することにある。
【００４９】対話型の使用ができない装置は人間の装置
との対話を最適化するため実時間ディスプレイが必要と
なるので著しく不利である。難しい手順を計画している
外科医の場合には、外科医が阻止構造を切り除くため問
題の構造を画面上に三次元形式でディスプレイすること
ができ、また不必要な遅れなしに随意に構造を回転させ
得ることが非常に重要である。
との対話を最適化するため実時間ディスプレイが必要と
なるので著しく不利である。難しい手順を計画している
外科医の場合には、外科医が阻止構造を切り除くため問
題の構造を画面上に三次元形式でディスプレイすること
ができ、また不必要な遅れなしに随意に構造を回転させ
得ることが非常に重要である。
【００５０】ディスプレイリスト構造を直接切断し得る
こと、およびその分離を行う際に画像の写像を使用し得
ることが本発明の迅速で対話型の処理に対する要点であ
る。本発明の装置および方法は断層撮影のＸ線および磁
気共鳴のイメージングデータの分析およびディスプレイ
に最大の用途があることは疑いがないが、本発明の装置
は超音波、ポジトロン放射形断層撮影、放射形コンピュ
ータ断層撮影、および多モードイメージングを用いる装
置にも同等に適用可能である。更に、本発明は特に医用
画像に適用可能であるが、本発明の装置および方法は三
次元ディスプレイリストを発生し得る任意の装置からの
三次元画像を切断することにも適用可能である。また、
上記の説明から明らかなように、それらのディスプレイ
リストの法線は勾配法線である必要はなく、他の手段に
よって得られる表面の方位を示す表面法線であってもよ
い。
こと、およびその分離を行う際に画像の写像を使用し得
ることが本発明の迅速で対話型の処理に対する要点であ
る。本発明の装置および方法は断層撮影のＸ線および磁
気共鳴のイメージングデータの分析およびディスプレイ
に最大の用途があることは疑いがないが、本発明の装置
は超音波、ポジトロン放射形断層撮影、放射形コンピュ
ータ断層撮影、および多モードイメージングを用いる装
置にも同等に適用可能である。更に、本発明は特に医用
画像に適用可能であるが、本発明の装置および方法は三
次元ディスプレイリストを発生し得る任意の装置からの
三次元画像を切断することにも適用可能である。また、
上記の説明から明らかなように、それらのディスプレイ
リストの法線は勾配法線である必要はなく、他の手段に
よって得られる表面の方位を示す表面法線であってもよ
い。
【００５１】ある実施例について本発明を詳細に説明し
てきたが、熟練した当業者は多数の変形および変更を加
えることができる。したがって、本発明の趣旨と範囲に
入るこのような変形および変更をすべて包含するように
特許請求の範囲を記載してある。
てきたが、熟練した当業者は多数の変形および変更を加
えることができる。したがって、本発明の趣旨と範囲に
入るこのような変形および変更をすべて包含するように
特許請求の範囲を記載してある。
【図１】三次元の各点でのイメージング対象の物理的性
質の測定データ値のボクセルを表わす斜視図である。
質の測定データ値のボクセルを表わす斜視図である。
【図２】定義された表面にまたがる図１の多数のボクセ
ルの立面図であり、それらのボクセルの頂点に対する勾
配法線を示す図である。
ルの立面図であり、それらのボクセルの頂点に対する勾
配法線を示す図である。
【図３】図２のボクセルの三次元点ベクトルのディスプ
レイリストを、ディスプレイのために二次元画像平面へ
写像する様式を図式的に表した概略図である。
レイリストを、ディスプレイのために二次元画像平面へ
写像する様式を図式的に表した概略図である。
【図４】図２のボクセルから作成されたディスプレイリ
ストを表す図表であり、その構造をコンピュータメモリ
のリストとして示す。
ストを表す図表であり、その構造をコンピュータメモリ
のリストとして示す。
【図５】図４のディスプレイリストからディスプレイ装
置上に二次元画像を発生するために使用される画像バッ
ファを表す図表であり、その構造をコンピュータメモリ
のアレーとして示す。
置上に二次元画像を発生するために使用される画像バッ
ファを表す図表であり、その構造をコンピュータメモリ
のアレーとして示す。
【図６】図４のディスプレイリストを、図５の画像バッ
ファへ写像するために使用されるｚバッファを表す図表
であり、その構造をコンピュータメモリのアレーとして
示す。
ファへ写像するために使用されるｚバッファを表す図表
であり、その構造をコンピュータメモリのアレーとして
示す。
【図７】本発明により図４のディスプレイリストを切断
して別々のディスプレイリストとするために使用される
カットバッファを表す図表であり、その構造をコンピュ
ータメモリのアレーとして示す。
して別々のディスプレイリストとするために使用される
カットバッファを表す図表であり、その構造をコンピュ
ータメモリのアレーとして示す。
【図８】図４のディスプレイリストで表される三次元表
面を切断するための本発明の方法のフローチャートであ
る。
面を切断するための本発明の方法のフローチャートであ
る。
【図９】本発明の方法および装置を示す概略ブロック図
である。
である。
【図１０】切断された多角形の構造を示すディスプレイ
装置上の画像の概略図である。
装置上の画像の概略図である。
【図１１】図１０のカット多角形で定義される体積の中
に含まれる点ベクトルのカットディスプレイリストから
作成された画像の概略図であり、ｙ’軸を中心として９
０°回転したカットデータを示す。
に含まれる点ベクトルのカットディスプレイリストから
作成された画像の概略図であり、ｙ’軸を中心として９
０°回転したカットデータを示す。
１４ 勾配法線
１６ 表面
１８ 原始ディスプレイリスト
１８’ 変換されたディスプレイリスト
２４ 画像平面
２６ 画素
３０ ＲＴＳ（回転／平行移動／スケール）機能ブロッ
ク ３４ ｚバッファ ３７ 比較ブロック ３９，４５ ゲート ４１ カラー計算器 ４２ 画像バッファ ４８ ディスプレイ装置 ５０ 多角形領域 ５１ 照明軸 ５６ カットバッファ ６０ カットディスプレイリスト ６２ 残余ディスプレイリスト ６５，６７ ゲート ７２ カット平面
ク ３４ ｚバッファ ３７ 比較ブロック ３９，４５ ゲート ４１ カラー計算器 ４２ 画像バッファ ４８ ディスプレイ装置 ５０ 多角形領域 ５１ 照明軸 ５６ カットバッファ ６０ カットディスプレイリスト ６２ 残余ディスプレイリスト ６５，６７ ゲート ７２ カット平面
─────────────────────────────────────────────────────
フロントページの続き
(51)Int.Cl.5 識別記号 庁内整理番号 ＦＩ 技術表示箇所
Ｇ０６Ｆ 15/20 Ｄ 7218−5Ｌ
(72)発明者 ウイリアム・エドワード・ロレンソン
アメリカ合衆国、ニユーヨーク州、ボルス
トン・レイク、ハースサイド・ドライブ、
14番
(72)発明者 シーグウオルト・ルドケ
アメリカ合衆国、ニユーヨーク州、スコテ
イア、ダウン・ドライブ、10番
Claims (11)
- 【請求項１】 対応する法線をそなえた三次元の各点の
第一のディスプレイリストとして表される表面構造に対
するカット軸に沿った手術切断のシミュレーションを行
うための装置に於いて、 各点の第一のディスプレイリストに作用して各点の座標
を回転させる変換手段、 各点を受けて、回転した座標を画像平面内の画素座標を
持つ画素に写像することにより、表面構造の二次元表現
を作成し、カット軸に平行な線に沿った各点が画像平面
内の単一の画素に写像されるようにする写像手段、 写像手段と通信して画像平面上の多角形領域を記述し、
また多角形領域に写像される各点を識別する手段、およ
び 各点の識別に応じていくつかの点を第二のディスプレイ
リストに記憶する手段を含むことを特徴とする手術切断
シミュレーション装置。 - 【請求項２】 画像平面内に多角形を対話的に作成する
データ入力装置、および表面構造と画像平面内の多角形
のスーパインポーズ画像をディスプレイするディスプレ
イ手段が含まれる請求項１記載の手術切断シミュレーシ
ョン装置。 - 【請求項３】 各点の第一のディスプレイリストに作用
して三次元カット平面の片側の各点を除去する手段が含
まれる請求項１記載の手術切断シミュレーション装置。 - 【請求項４】 上記各点の識別に応答して記憶する手段
により、回転した座標が多角形領域に写像される各点の
第二のディスプレイリストと、回転した座標が多角形領
域に写像されない各点の第三のディスプレイリストとが
作成される請求項１記載の手術切断シミュレーション装
置。 - 【請求項５】 上記写像手段と通信して多角形領域を記
述するための手段には、画素座標によってアドレッシン
グされ、画素座標が多角形領域内にあれば第一の値を保
持するカットバッファが含まれる請求項１記載の手術切
断シミュレーション装置。 - 【請求項６】 対応する法線をそなえた三次元の各点の
第一のディスプレイリストとして表される表面構造に対
するカット軸に沿った手術切断のシミュレーションを行
うための装置に於いて、 各点の第一のディスプレイリストを受けて、各点に対す
る画像平面座標を発生して各点を画像平面の画素に写像
することにより表面構造の二次元表現を作成する写像手
段、 写像手段と通信して画像平面上の多角形領域を記述し、
また多角形領域に写像される各点を識別する手段、 画像平面内の多角形領域を記述するカットバッファであ
って、画像平面座標でアドレッシングされ、画像平面座
標が多角形領域内にあれば第一の値を保持するカットバ
ッファ、 カットバッファと通信し、画像平面座標に等しいカット
バッファのアドレスが第一の値を保持しているかチェッ
クすることにより上記各点の画像平面座標を多角形領域
の画像平面座標と比較する比較手段、および上記比較手
段による比較の結果に応じて各点の第一のディスプレイ
リストから各点の第二のディスプレイリストを発生する
分離手段を含むことを特徴とする手術切断シミュレーシ
ョン装置。 - 【請求項７】 上記分離手段はある点の画像平面座標に
等しいカットバッファのアドレスが第一の値を保持して
いれば上記点を各点の第二のディスプレイリストに記憶
させ、そうでない場合には上記点を第三のディスプレイ
リストに記憶させる請求項６記載の手術切断シミュレー
ション装置。 - 【請求項８】 第一のディスプレイリストの中の各点が
カラーに対応しており、当該装置は、第一のディスプレ
イリストの各点に対する画像平面座標を受ける画像バッ
ファ手段であって、画像平面に垂直な次元に沿った点の
ｚ座標が同じ画像平面座標を有する、前に記憶されたど
の点よりも画像平面に近い場合に上記画像平面座標の点
のカラーを記憶し、それらの座標に対する前に記憶され
たカラーを消去する画像バッファ手段、および対応する
画像平面座標の記憶されたカラーをディスプレイするデ
ィスプレイ手段を含んでいる請求項６記載の手術切断シ
ミュレーション装置。 - 【請求項９】 上記カラーに、上記点に対応する法線と
所定の照明軸とのドット積が含まれる請求項８記載の手
術切断シミュレーション装置。 - 【請求項１０】 画像平面内に多角形領域を対話的に作
成するデータ入力装置、および表面構造と画像平面内の
多角形のスーパインポーズ画像をディスプレイするディ
スプレイ手段が含まれる請求項６記載の手術切断シミュ
レーション装置。 - 【請求項１１】 各点の第一のディスプレイリストに作
用して三次元カット平面の片側の各点を除去する手段が
含まれる請求項６記載の手術切断シミュレーション装
置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US66333391A | 1991-03-01 | 1991-03-01 | |
| US663333 | 1991-03-01 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH057554A true JPH057554A (ja) | 1993-01-19 |
| JPH0793926B2 JPH0793926B2 (ja) | 1995-10-11 |
Family
ID=24661376
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3680192A Expired - Lifetime JPH0793926B2 (ja) | 1991-03-01 | 1992-02-25 | ディスプレイリスト表面データを使用する手術シミュレーション装置 |
Country Status (2)
| Country | Link |
|---|---|
| EP (1) | EP0501812A2 (ja) |
| JP (1) | JPH0793926B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010022602A (ja) * | 2008-07-18 | 2010-02-04 | Fujifilm Ri Pharma Co Ltd | 器官表面画像の表示装置及び方法 |
| JP4824034B2 (ja) * | 2004-12-23 | 2011-11-24 | エレクタ、アクチボラグ | 合成記録ベクトルを用いた最適化された画像記録 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7256780B2 (en) * | 2004-03-04 | 2007-08-14 | Siemens Medical Solutions Usa, Inc. | Visualization of volume-rendered data with occluding contour multi-planar-reformats |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60102318A (ja) * | 1983-11-07 | 1985-06-06 | Daifuku Co Ltd | パレツト位置決め装置 |
| JPS6437678A (en) * | 1987-08-03 | 1989-02-08 | Toshiba Corp | Three-dimensional image processor |
| JPH01135339A (ja) * | 1987-11-19 | 1989-05-29 | Toshiba Corp | 三次元画像診断装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4719585A (en) * | 1985-08-28 | 1988-01-12 | General Electric Company | Dividing cubes system and method for the display of surface structures contained within the interior region of a solid body |
| US4985834A (en) * | 1988-11-22 | 1991-01-15 | General Electric Company | System and method employing pipelined parallel circuit architecture for displaying surface structures of the interior region of a solid body |
-
1992
- 1992-02-25 JP JP3680192A patent/JPH0793926B2/ja not_active Expired - Lifetime
- 1992-02-28 EP EP92301694A patent/EP0501812A2/en not_active Withdrawn
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60102318A (ja) * | 1983-11-07 | 1985-06-06 | Daifuku Co Ltd | パレツト位置決め装置 |
| JPS6437678A (en) * | 1987-08-03 | 1989-02-08 | Toshiba Corp | Three-dimensional image processor |
| JPH01135339A (ja) * | 1987-11-19 | 1989-05-29 | Toshiba Corp | 三次元画像診断装置 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4824034B2 (ja) * | 2004-12-23 | 2011-11-24 | エレクタ、アクチボラグ | 合成記録ベクトルを用いた最適化された画像記録 |
| JP2010022602A (ja) * | 2008-07-18 | 2010-02-04 | Fujifilm Ri Pharma Co Ltd | 器官表面画像の表示装置及び方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0793926B2 (ja) | 1995-10-11 |
| EP0501812A2 (en) | 1992-09-02 |
| EP0501812A3 (ja) | 1994-02-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5995108A (en) | 3D image composition/display apparatus and composition method based on front-to-back order of plural 2D projected images | |
| Reynolds et al. | A dynamic screen technique for shaded graphics display of slice-represented objects | |
| US6181348B1 (en) | Method for selective volume visualization via texture mapping | |
| US4821210A (en) | Fast display of three-dimensional images | |
| CN109584349B (zh) | 用于渲染材料属性的方法和设备 | |
| US5113357A (en) | Method and apparatus for rendering of geometric volumes | |
| JP3483929B2 (ja) | ３次元画像生成方法 | |
| US6990231B2 (en) | Method and apparatus for forming and displaying projection image from a plurality of sectional images | |
| Höhne et al. | Rendering tomographic volume data: adequacy of methods for different modalities and organs | |
| JPH0437465B2 (ja) | ||
| EP3879498A1 (en) | Method of rendering a volume and a surface embedded in the volume | |
| Haubner et al. | Virtual reality in medicine-computer graphics and interaction techniques | |
| Englmeier et al. | Hybrid rendering of multidimensional image data | |
| Turlington et al. | New techniques for efficient sliding thin-slab volume visualization | |
| CN100583161C (zh) | 用于显示以立体数据组成像的对象的方法 | |
| US6392646B1 (en) | Iterative determination of the shortest path between two points on a polygonal surface | |
| JPH057554A (ja) | デイスプレイリスト表面データを使用する手術シミユレーシヨン装置 | |
| Dietrich et al. | Real-time interactive visualization and manipulation of the volumetric data using GPU-based methods | |
| Bornik et al. | Interactive editing of segmented volumetric datasets in a hybrid 2D/3D virtual environment | |
| JPH05189543A (ja) | 三次元画像情報提示方法 | |
| Stewart et al. | Rebuilding the visible man | |
| Tiede et al. | Surface rendering | |
| i Bartrolı et al. | Visualization techniques for virtual endoscopy | |
| Löbbert | Visualisation of two-dimensional volumes | |
| JP3326526B2 (ja) | コンピュータグラフィックス表示装置にボリューム物体をレンダリングする方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19960416 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081011 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091011 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091011 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101011 Year of fee payment: 15 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101011 Year of fee payment: 15 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111011 Year of fee payment: 16 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121011 Year of fee payment: 17 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121011 Year of fee payment: 17 |