JPH02211587A - 物体内部構造表面の2次元像を表示する装置と方法 - Google Patents

物体内部構造表面の2次元像を表示する装置と方法

Info

Publication number
JPH02211587A
JPH02211587A JP1295187A JP29518789A JPH02211587A JP H02211587 A JPH02211587 A JP H02211587A JP 1295187 A JP1295187 A JP 1295187A JP 29518789 A JP29518789 A JP 29518789A JP H02211587 A JPH02211587 A JP H02211587A
Authority
JP
Japan
Prior art keywords
volume element
volume
elements
numerical value
values
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
Application number
JP1295187A
Other languages
English (en)
Other versions
JP2744490B2 (ja
Inventor
Harvey E Cline
ハーベイ・エリス・クライン
Richard I Hartley
リチャード・イアン・ハートレイ
Siegwalt Ludke
シーグワルト・ルドゲ
Sharbel E Noujaim
シャーベル・エミル・ノージャイム
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of JPH02211587A publication Critical patent/JPH02211587A/ja
Application granted granted Critical
Publication of JP2744490B2 publication Critical patent/JP2744490B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three-dimensional [3D] modelling for computer graphics

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Generation (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 技術分野 この発明は表示装置、更に具体的に云えば、密実な物体
の内部にある規則的な格子位置に於けるその物理的な性
質の数値を利用し得る様な物体の内面構造を速やかに表
示することに関する。
発明の背景 密実な物体の内部にある規則的な格子位置に於ける1種
類又は更に多くの物理的な性質を表わすデータの3次元
配列を求めることは周知である。
こう云うデータは、計算機式軸断層写真法(CAT)X
線走査装置、核磁気共鳴(NMR)作像装置の様な非侵
入形方法、又は超音波又はポジトロン放出断層写真法(
PET)、放出形計算機式断層写真法(ECT)及び多
重モード作像(MMI)の様な他の非侵入形機構によっ
て求めることが出来る。こう云う方式の各々は、密実な
物体の一連のスライスの各々に対する平面状の格子状配
列の数値を発生し、こう云う数値の3次元配列を作る。
密実な物体は典型的には人体又はその一部分であるが、
この方法はその他の天然又は人工の物体にも等しく応用
し得る。CAT走査の場合、物理的な数値はX線吸収係
数である。NMR作像では、物理的な数値はスピン−ス
ピン又はスピン−格子緩和時間である。何れにせよ、測
定された物理的な数値が、その下にある物理的な構造の
組成、密度又は表面の特性の変化を反映している。典型
的には、この様な3次元データ配列は、物体内の立方又
は平行六面体の格子の規則的な位置に分布した3次元(
x、y、z)座標の複数個の組と、夫々の座標位置に伴
う物理的な性質の少なくとも1つの数値(V   )と
で構成されている。立方体yz として隣接する各々の8個1組の位置が「容積要素」と
呼ばれる立方体容積を限定し、容積要素の8個の頂点の
各々に対する物理的な性質の数値が特定される。各々の
容積要素の「近隣」は、この容積要素自体と、面を共有
する直ぐ隣接する6個の容積要素とを含む。即ち容積要
素の近隣は、容積要素の頂点に関連する32個の物理的
な数値を持つ7個の容積要素を含む3次元容積である。
同様に、この様な内部の物理的な数値の3次元配列を利
用して、物体内にある内部構造の可視像を発生すること
も知られている。人体の場合、こうして発生される可視
像は、診断の様な医療用とか或いは外科手順の計画の為
に使うことが出来る。
然し、この様な3次元の内部構造の2次元像を表示する
為には、物理的な数値の配列内にあるこの様な構造の面
の位置を突止めることが必要である。
これを行なう為、配列の数値を、その面に関連する物理
的な性質の数値に対応する1個の閾値又はある範囲の閾
値と比較する。例えば、骨或いはその他の任意の組織は
既知のある範囲の密度の値によって特徴づけることが出
来、配列の数値をそれに対して比較することが出来る。
−旦面の位置が決定されたら、この面に陰影を施して、
2次元表示装置で表示した時、その面の形及び配置の正
しい印象が人間の目に与えられる様にしなければならな
い。この陰影を施す為、面上の各点に於ける面に対する
法線方向又は直交するベクトルの角度方向を観察者の観
察角度と比較する。その後、陰影の強度を調節して、こ
う云う角度の間の差に比例する様にする。この角度差情
報を使って、表示される像に取入れる色をも制御し、こ
うして面の配置を判断する視覚上の別の手掛かりを与え
ることが出来る。観察角度からそれた方向を向く成分を
持つ法線ベクトルは、関連する面が目に見えないから、
無視することが出来る。
内部構造の面を近似する1つの方法が、米国特許節4,
710.876号に記載された所謂「行進キューブ(m
arching cubos) J方法である。この方
法では、容積要素と交差する面のセグメントが、その容
積要素と交差する限られた数の標準化した平面状の多角
形面の内の1つによって近似される。閾値と容積要素の
8個の頂点の数値との間の2進法で表わした差を表わす
ベクトルにより、特定の1つの標準化した面が選択され
る。この様な標準化した多角形面の各々の組に対する、
面と容積要素の交点の座標並びに法線ベクトルを計算す
るか、或いはテーブル・ルックアップ方式によって求め
ることが出来る。標準化した全ての多角形を切片又はタ
イルとして使って、最終的な面をモザイクに組立てる。
法線ベクトルの角度から導き出した適当な強度の値を直
ちに観察する為に表示してもよいし、或いは後で表示す
る為に記憶することが出来る。米国特許節4.729,
098号には、切片状の標準化した多角形の座標を更に
正確に突止める為に、非一次補間を用いた変形の行進キ
ューブ方法が記載されている。
内部構造の面を近似する別の方法が、米国特許節4,7
19,585号に記載された所謂[分割キューブ(di
viding cubes) J方法である。この方法
では、容積要素の頂点に於ける数値を使って、容積要素
内部の規則的な位置にある部分格子位置に於ける数値を
3次元で補間する。こうして補間した部分格子の数値を
使って、面の位置を更に正確に突止めると共に、法線ベ
クトルを更に正確に計算する。この発明は面の位置及び
面の法線を近似する分割キューブ方法の拡張及び改善で
ある。
物理的な性質の数値が同じであるか類似している異なる
内部構造を区別する為、米国特許節4゜751.643
号には、類似した性質の数値を持つ面にラベルをつけ、
関心のある特定の構造内にある「シード(seed) 
J位置に対する隣接性の判断基準を使って、所望の面を
ラベルをつけた全ての面から分離する方法が記載されて
いる。1986年9月15日に出願された係属中の米国
特許出願通し番号節907.333号には、隣接の情報
から連結性を判断することにより、類似する構造を分離
する別の方式が記載されている。1986年12月19
日に出願された係属中の米国特許出願節943,357
号には、内部構造を区別する為の更に別の方式が記載さ
゛れている。この方式では、データ配列を直線的に通っ
て、走査線に沿った全ての異なる構造の位置を、構造の
界面を計数することによって、突止めてラベルをつける
この為、密実な物体の内部にある物理的な性質の数値の
1個の配列を使って、単に同じ1個の配列の数値の操作
により、任意に選ばれた観察角度から見た様な、物体内
の任意に選ばれた内部構造の斜視図を発生することが知
られている。都合の悪いことに、この様にして1個の像
を作るのに必要な処理量は非常に大きく、従来は多くの
進行中の手順と合せる程、像を手早く発生することが困
難であった。更に、外科手順の複雑さ並びに微妙さが極
く最近高まるのにつれて、進行中の医療又は外科手順を
助ける様に、こう云う像を実時間で対話形で発生する必
要も高まりつ−ある。器官の移植、心臓のバイパス手術
及び美容外科は、医学的な知識の境界が拡がるにつれて
、複雑さ並びに精密度を何れも増している。従って、進
行中の手術の途中に、必要に応じて多数の実時間の3次
元像を発生する必要性が高まりつ\ある。
従って、この発明の1つの目的は、任意に選ばれた観察
角度から観察した様な、3次元物体内にある構造の2次
元の斜視像を実時間で発生する方法と装置を提供するこ
とである。
この発明の別の目的は、3次元物体の内部構造の像を比
較的速い速度で発生することが出来る高度に並列のバイ
ブライン形プロセッサを提供することである。
この発明の別の目的は、物体内の3次元構造の2次元の
斜視像を実時間で対話形で発生する方法と装置を提供す
ることである。
発明の要約 この発明の実施例では、内面の2次元の斜視像が、高度
に並列のバイブライン形プロセッサにより、物理的な数
値の3次元の規則的な配列から発生される。このプロセ
ッサが3次元物体の内部構造の像を進行中の外科手順と
合さる位に速やかに発生する。更に具体的に云うと、3
次元データが多重バンク・メモリに記憶されていて、「
容積要素」のデータを並列にアクセスすることが出来る
関心のある1個の容積要素の近隣にある複数個のこの様
な容積要素からのデータが、2つの並列通路によって同
時に処理される。その一方は表示すべき面の位置を設定
し、他方は面上の関連した位置に於ける陰影のレベルを
決定する。関心のある特定の面は、関心のある面に関連
する物理的な定数の数値を表わす閾値定数によって選択
することが出来る。交差する容積要素を突止める為、3
次元データを持っているメモリから、容積要素の頂点の
数値が一度に1組ずつアクセスされる。こうしてアクセ
スされた容積要素の頂点の数値の各組をこの定数と比較
して、関心のある面と交差する容積要素を確認する。交
差しない容積要素は除外し、次の容積要素をアクセスす
る。この様にして見つけられた各々の交差する容積要素
に対し、「近隣」の容積要素に伴うデータの数値をアク
セスして、面の法線の数値を計算することが出来る様に
する。ある容積要素の「近隣」とは、その容積要素を取
囲み且つそれと接する3X3の容積要素の容積内にある
、その容積要素と面を突合せる6個の容積要素を指す。
各々の交差する容積要素に対し、そして交差した容積要
素に対してだけ、面に対して垂直なベクトルの方向(面
の法線の数値)を容積要素の各々の頂点に対して計算す
る。
その後、交差した容積要素を予め選ばれた倍率で部分容
積要素の立方体に分割し、部分容積要素の頂点の数値を
一次補間によって補間して、面と部分容積要素の立方体
との交点を更に精密に突止める。この目的の為、部分容
積要素の頂点の数値を面の定数と比較し、交差する部分
容積要素の立方体を確認する。部分容積要素の立方体の
1個の点(例えば中心)が関心のある面上にあると見な
される。この面上にある関心のある点をx、y、z座標
で表わすが、それを使って容積要素の頂点の法線の間の
3次元補間を実施して、こうして確認された面の点に於
ける法線の角度を求める。その後、普通のグラフ式表示
手段を使って、観察角度及び観察の奥行を選択し、面の
座標及び法線の数値を面の2次元の図に変換し、面を陰
極線管スクリーンに表示する。
この発明の重要な特徴は、プロセッサのうちの、面を突
止めて法線を計算する部分で行なわれる処理が高度に並
列であることである。更に、面との交差部を持つ容積要
素及び部分容積要素だけを処理しさえすればよく、この
為、必要な処理量が大幅に減少する。最後に、容積要素
及び部分容積要素の比較は同時に行なう必要は全くない
から、比較回路をこう云うZFJ類の比較で時分割で用
いても、プロセッサの動作を遅くすることはない。高度
に並列の処理並びにパイプライン・アーキテクチュアを
使うことにより、内面の構造を実時間で発生することが
出来、進行中の手順の助けとして、こう云う面の表示を
使うことが出来る様になる。
−旦関心のある面全体に対する面の点が蓄積されたら、
それを任意の観察する向きに任意に回転させ、面の点の
操作により、データ空間又は観察者の空間で切取ること
が出来る。
最終的な像の分解能は、部分容積要素の倍率を表示装置
の分解能とぴったり合う様に選ぶことによって制御され
、こうして表示装置の分解能を一杯に活用する。
この発明は以下図面について詳しく説明する所から、更
に十分に理解されよう。
詳しい説明 この発明では、一連の容積要素を用いて、関心のある面
と容積要素との交点を突止める。こ\で云う[容積要素
(voxel ) Jは3次元の平行六面体であって、
非侵入形手段によって、密実な物体の内部で測定される
物理的な数値を表わすデータの数値の配列によって限定
される。この発明では、相次ぐNMR又はCAT走査の
スライスからの、データを数値の3次元配列に組立て、
その後それを解析し又は処理して、3次元情報の2次元
像を求めることが出来る。この様な3次元の内部データ
配列を発生することは周知であり、こ\では詳しく説明
しない。この様なデータ配列が、計算機式%式% 鳴(NMR)作像装置、超音波走査、ポジトロン放出断
層写真法(PET)、放出計算機式断層写真法(ECT
)及び多重モード作像(MMI)の様な周知の非侵入形
の方法によって容易に得られることを述べておけば十分
である。こう云う方法により、データの点の平面状の配
列が得られる。
走査する密実な物体を通る規則的な一連の隣接する「ス
ライス」の各々で、この様な1つの平面状の配列が得ら
れる。それを合せると、一連のスライスがデータの数値
の3次元配列を形成する。
この様なデータの数値が配列になっていることが第1図
に見られる。この図は、頂点vl乃至V8を持つ1個の
容積要素10を示している。用途10の様な各々の容積
要素は、データの数値の相次ぐ2つのスライスに跨がっ
ている。容積要素10の各々の頂点Vl乃至■8には、
3次元物体内の対応する空間位置に伴う少なくとも1種
類の物理的な性質の測定値を表わすデータの数値が関係
している。この空間的な位置は、物体内に規則的に隔た
る格子位置を定める規則的なパターンをなす所にある。
格子位置が容積要素10の様な複数個の隣接する容積要
素を定める。
この発明では、密実な物体の内部構造の可視表示を実時
間で、即ち、この内部構造の知識が何よりも重要である
様な進行中の外科又はその他の手順で助けになる位に速
やかに発生する為、データ配列が高度に並列のパイプラ
イン回路アーキテクチュアで記憶され、アクセスされ、
そして処理される。
この発明に関係する手順を簡単にまとめると、次の通り
である。
1、CAT走査又はNMR作像の様な非侵入形手段によ
り、物理的な数値の3次元配列を求めて記憶する。、・ 2、規則的な順序で、−度に1つの容積要素のデータ(
8個のデータの数値)の3次元配列をアクセスする。
36現在の容積要素の頂点の数値を関心のある面を表わ
す定数と比較する。
4、工程(3)の判断により、関心のある面がこの容積
要素と交差しない場合、次の容積要素のデータの数値を
アクセスし、工程(3)を繰返す。
5、関心のある面がこの容積要素と実際に交差すれば、
近隣の数値(24個のデータの数値)を再生し、こう云
う数値を使って、有限の差により、交差した容積要素の
頂点に於ける法線の数値を計算する。
6、分割キューブ方法を用いて、交差した容積要素を、
表示装置の分解能を十分に活用するのに十分な位の小さ
い寸法を持つ複数の部分容積要素(sub−voxel
 )に分割する。
76 この様にして定めた各々の部分容積要素に対し、
容積要素の頂点のデータの数値及び−次補間を用いて、
各々の部分容積要素の頂点に対する物理的なデータの数
値を補間する。
8、−度に1つの部分容積要素ずつ、部分容積要素の頂
点の数値を関心のある面を表わす定数と比較する。
9、関心のある面がこの部分容積要素と交差しない場合
、交差した容積要素にある次の部分容積要素に進み、工
程(8)を繰返す。
10、関心のある面がこの部分容積要素と交差すれば、
その部分容積要素の座標を関心のある面上のある点の座
標として利用する。
11、こう云う部分容積要素の座標を使って、工程(5
)で計算した法線の数値の3次元−次補間を行なう。
12、工程(10)及び(11)で求められた面の座標
及び法線ベクトルの数値を表示リストに記憶し、全ての
部分容積要素を試験するまで、工程(8)に戻る。
13、現在の交差した容積要素の全ての部分容積要素が
工程(8)で試験された時、工程(2)に戻って、配列
の次の容積要素をアクセスする。
全ての容積要素が試験されたら、工程(2)の容積要素
走査手順を終了する。
14、工程(12)で表示リストに記憶された面の位置
の座標を、観察の角度及び奥行によって決定された観察
平面に投影する。
15、切断平面を選んで、関心のある面並びに観察窓の
両方を制限し、それに応じて観察平面の座標の一波作用
をする。
16、観察角度を面のその位置に於ける法線の角度の数
値と比較することにより、こうして投影して)一波作用
にかけた面の位置に於ける表示の強度(又は色)を計算
する。
17、工程(15)で−波作用にかけられた画像の座標
で、かつ投影された面の座標に於ける工程(16)で計
算された強度レベルで、表示装置に表示する。
次にこの手順を他の図面について更に詳しく説明する。
前に述べた様に、第1図は1個の容積要素を示す。然し
、法線の数値を計算する為には、近隣の容積要素の頂点
の数値が必要である。種々の近隣の物理的な数値並びに
容積要素の頂点に対するその座標位置が第2図に更によ
く示されている。
第2図は頂点■1乃至v8を持つ1個の容積要素を示す
。前に述べた様に、各々の容積要素が情報配列の隣接す
る2つのスライスに跨がる。各々の頂点には、その位置
に於ける3次元物体の物理的な性質の測定値を表わす数
が関連している。立方体として隣接する8個の格子位置
が容積要素を構成する。容積要素の8個の頂点Vl乃至
v8の各々には、容積要素の辺の投影に沿った所にある
隣接する3つの格子位置が関連している。この様な「近
い隣り」の格子位置をWl乃至W24で示しであるが、
これが容積要素の近隣を限定する。
近い隣りのデータ点が24個ある。第2図では、容積要
素の頂点Vl乃至v8は塗りつぶした丸で示し、近隣の
データ点W1乃至W24は塗りつぶさない丸で示しであ
る。各々の容積要素に対する関心のあるデータの数値(
Vl乃至v8及びWl乃至W24)が隣接する4つのデ
ータスライスを占める。この発明の装置及び方法では、
こういう32個のデータの数値を使って、容積要素の頂
点に於ける、並びに補間によって、容積要素内にある面
上の点に於ける法線ベクトルを計算する。こう云う法線
ベクトルは有限差分方法を用いて計算され、特に中心差
分方法を用いて計算される。更に具体的に云うと、頂点
v1に於ける法線ベクトルのX成分を計算するには、格
子位置V4及びW5に於けるデータの数値の差を求める
。同様に、頂点V1に伴う法線ベクトルの2成分は、格
子位置V5及びWlに於けるデータの数値の差によって
計算することが出来る。最後に、頂点V1に関連する法
線ベクトルのy成分は、格子位置V2及びWl2に於け
るデータの数値の差として発生される。差分から形成さ
れたベクトルを単位の大きさに倍率を変える為、各々の
ベクトル成分を、全てのベクトル成分の自乗の和の平方
根で除す。
第3図はこの発明の装置の簡略ブロック図である。バル
ク・メモリ3次元データ記憶装置20が、CAT走査又
はNMR走査方法によって得られた物理的な性質の数値
の配列を記憶している。図示の実施例では、記憶装置2
0はRAMメモリの複数個のバンクで構成される。この
為、毎回のアクセスで、複数個の数値が容積要素データ
記憶装置又はレジスタ21に転送される。記憶装置20
のバンクのアドレスは、相次ぐ容積要素の座標に対応す
る。この為、アドレス発生器22にある1個のアドレス
が、ある容積要素全体のデータの数値(第1図に示す様
に8個の数値)を記憶装置20から容積要素レジスタ2
1に転送するのに役立つ。
後で第4図について説明するが、アドレス発生器22は
xrY及び2カウンタで構成することが出来る。そのカ
ウントを合せたものが、メモリ20の容積要素アドレス
を構成する。こう云うカウンタは最初はx=1.y−1
,z−1の容積要素に設定される。発生器22を進める
度に、Xカウントが、データ配列の最大のX寸法まで、
前進する。
データ空間の縁で、Xカウントを1にリセットし、Xカ
ウントを1だけ増数する。Xカウントが容積要素のデー
タ空間の縁に達した時、Xカウントを1にリセットし、
2カウントを1だけ増数する。
こうしてデータ空間全体を、−度に1つの容積要素ずつ
、規則的に走査する。
格子位置vl乃至v8を表わす8個のデータ点が容積要
素記憶装置21に記憶され、オア・ゲート36を介して
8数値比較回路26に印加される。
比較回路26は第5図について詳しく説明する。
比較回路26には導線24を介して面の定数も印加され
る。導線24の面の定数は、表示すべき面に関連する物
理的な性質の数値に対応する。2つ以上の面が同じ又は
類似した面の定数の値を持つ場合、前に述べた米国特許
第4,751.643号又は係属中の米国特許出願通し
番号第907゜333号又は同第943.357号の方
式を用いて、関心のある面を切離すことが出来る。比較
回路26が導線24からの面の定数を、容積要素レジス
タ21に記憶されている容積要素の8個の頂点の数値v
l乃至V8と比較する。容積要素の全ての頂点の数値が
面の定数より大きければ、容積要素は完全に関心のある
面の外側にある。容積要素の全ての頂点の数値が面の定
数より小さければ、この容積要素は完全に関心のある面
の内側にある。
比較によってこの他の結果が出れば、関心のある面はこ
の容積要素と交差する。面の像に寄与するのは交差する
容積要素だけであるから、容積要素が交差しなければ、
導線25を介してアドレス発生器22を前進させ、次の
容積要素からのデータを検査の為に供給する。実際、大
抵の面は、データ配列にある極く小さな百分率の容積要
素だけと交差するので、大抵の容積要素は比較回路26
で試験した後は、それ以上考慮されない。
関心のある面がこの容積要素と交差しない場合、更に処
理するのがよい。その場合、アドレス発生器22に導線
25から信号を加えて、現在の容積要素の最も接近した
6個の隣りの容積要素に対する容積要素アドレスを発生
する。即ち、x、y及び2カウンタを一度に1つずつ、
−時的に1だけ増数並びに減数して、関心のある容積要
素と共通の面を持つ6個の容積要素の頂点の数値をアク
セスする。こう云う隣りの容積要素から求められた24
個の最も近い隣りの容積要素の頂点の数値Wl乃至W2
4(第2図に示す)が、近隣レジスタ23に記憶される
。近隣レジスタ23に全ての近い隣りの頂点の数値Wl
乃至W24が蓄積された時、それを法線ベクトル計算器
27に印加する。
法線ベクトル計算器27が、頂点の数値vl乃至v8に
よって定められた容積要素の8個の頂点の各々に対する
法線ベクトルのx、y及び2成分を計算して正規化する
が、これは第6図について説明する。24個の法線ベク
トルの数値(容積要素の頂点毎に3個)が法線レジスタ
28に記憶される。
法線ベクトル計算器27が頂点の法線の数値を発生して
いるのと同時に、比較回路26によって付能された容積
要素分割器29が、容積要素を導線30から送られた倍
率Sx +  Sy +  Stにより、複数個の部分
キューブに分割する。こう云う倍率は、使う表示装置の
解像度と正確に合う様に選ばれる。各々の次元に対して
別々の倍率を用いてもよいが、3つの次元全部に同じ倍
率を使ってもよい。倍率は、第3図及び第7図の回路に
要求される処理を少なくする為に、2進法(1,2,4
゜8等)にすることが好ましい。後で第7図について説
明するが、容積要素分割器29が、部分容積要素のキュ
ーブの各々の頂点に対する端数頂点座標を発生する。こ
う云う端数頂点座標を容積要素補間器35で用いて、レ
ジスタ21にある容積要素の頂点の数値から部分容積要
素の頂点のデータ数値を発生する。端数頂点座標は、交
差した部分容積要素に対する部分容積要素の法線の数値
を発生する為に、法線補間回路31でも使われる。容積
要素補間回路35で得られた部分容積要素の頂点の数値
が、−度に1つの部分容積要素キューブ毎に、オア・ゲ
ート36を介して、比較回路26で面の一定の値と比較
される。容積要素自体の場合と同じく、部分容積要素キ
ューブの交差は、部分容積要素の全ての頂点の数値が面
の一定の値より全部大きいか或いは全部それより小さい
ことの何れでもないことによって確認される。容積要素
分割器29が、交差する部分容積要素を捜して、容積要
素自体が走査されるのと大体同じ様に、部分容積要素キ
ューブを走査する。走査する部分容積要素が見つかった
時、法線補間器31の出力が付能されて、レジスタ28
にある法線の数値の3次元の補間を補間器31から表示
装置34に転送する。交差した部分容積要素に対する部
分容積要素の座標が、それらの座標にある法線の数値が
法線補間器31から導線33を介して送られるのとと同
期して、分割器29から導線32を介して表示装置34
に同時に供給される。容積要素分割器29は後で第7図
について更に詳しく説明し、容積要素補間回路35は第
9図について詳しく説明し、表示装置34は第10図に
ついて更に詳しく説明する。法線補間器31は第9図に
示す形式の補間回路を3個を用いて構成される。
第3図に示す処理装置によって実施される手順全体は高
度に並列であることが認められよう。並列処理回路で面
の位置及び面の法線が同時に計算される。更に、処理工
程はパイプラインで実施され、パイプラインの各々の素
子がパイプラインの先行する素子からのデータをその入
力に受取り、その出力データをパイプラインの後続の素
子に送出し、その人力にデータが渡されること以外は、
どんな動作をも待つことがない。従って、パイプライン
の種々の点に於ける種々のプロセスが、相異なる次々の
組の数値に対して譬はあるが、同時に実施される。他の
図面の回路に見られるが、個々のプロセスがそれ自体も
高度に並列であり、この為密実な物体の内面を実時間で
可視的に表示するのに十分高速の全体的な処理装置を使
うのが容易になって、外科手術の様な進行中の手順にと
って貴重な助けとなる。
第4図は第3図のアドレス発生器22の詳しい回路図で
あり、アドレス発生器はXカウンタ40、Xカウンタ4
1及び2カウンタ42を有する。各々のカウンタはデー
タ数値の配列の対応する次元に於ける容積要素の合計の
数までカウントアツプする。更に、Xカウンタ40は、
それがサイクルを終る度に、Xカウンタ41を前進させ
、Xカウンタ41は、それがサイクルを終る度に、2カ
ウンタ42を前進させる。カウンタ40乃至42にある
カウントが、カウンタ40乃至42に対して1つずつあ
る部分47,48.49で、容積要素アドレス・レジス
タ45に記憶される。これらの3つのカウンタを合せた
ものが、第3図に示したバルク・メモリ20のアドレス
を形成し、そのアドレスが、レジスタ45の内容によっ
て同定された容積要素に対する8個の頂点の数値vl乃
至V8の全部をアドレスする。バルク・メモリ20は、
例えば何れも各々のアドレスにある容積要素の頂点の数
値を記憶する記憶素子の8個のバンクで構成することが
出来る。メモリ20の異なるバンクの同一のアドレスに
記憶される数値が、勿論同じ容積要素に対する8個の頂
点の数値である。従って、1個のアドレスが容積要素の
8個の頂点の数値を同時にアクセスするのに役立つ。導
線43(第3図の比較回路26から)の前進パルスがカ
ウンタ40乃至42を次の容積要素アドレスへ前進さす
、次の容積要素に対する8個の頂点の数値をアクセスす
る。
関心のある面が、(第3図の比較回路26により)レジ
スタ45にあるアドレスによって同定された容積要素と
交差すると決定されると、導線44が付勢されて、増数
及び減数回路50乃至55を循環的に付能する。カウン
タ40乃至42のカウントを系統的に増数及び減数する
ことにより、レジスタ45にある相次ぐアドレスによっ
て、近隣の容積要素が同定される。こう云う近隣の容積
要素アドレスがバルク・メモリ20から、対応する容積
要素の頂点の数値を検索するのに役立つ。
Wl乃至W24(第2図に示す)の近隣の数値が、この
為第3図の近隣レジスタ23に記憶する為に利用出来る
様になる。具体的に云うと、Xカウンタ41及び2カウ
ンタ42のカウントを一定に保ち、Xカウンタ40のカ
ウントを1だけ増数すると、容積要素Vl−V2−V5
−V6−W5−W6−Wl 3−Wl 4 (第2図)
がアクセスされ、数値W5. W6. Wl 3. W
l 4を(第3図に示す)レジスタ23に記憶すること
が出来る。その後、Xカウンタ20のカウントを1だけ
減数すると、容積要素V3−V4−V7−V8−W9−
Wl 0−Wl 7−Wl gがアクセスされ、近隣の
数値W9.Wl O,Wl 7.Wl 8をレジスタ2
3に記憶することが出来る。同様に、他の対のカウンタ
40乃至42を一定に保ち、3番目のカウンタを増数及
び減数すると、レジスタ23に記憶する為、他の全ての
近隣の数値Wl乃至W24を持つ容積要素をアクセスす
る様に作用する。
カウンタ40乃至42にあるカウントを予め選ばれた範
囲の数値に(容易に考えられる回路の変更によって)任
意に制限して、データ空間の制限された容積だけを走査
することが出来ることは容易に明らかである。関心のあ
る面又は関心のある面の所望の部分が、データ空間の予
定の容積に制限されることが判っている場合、その制限
された空間内にある容積要素だけを処理するのが、処理
時間の節約になる。関心のある容積が前もって判ってい
なければ、−旦所望の容積が判った時に、表示されるデ
ータを制限することも勿論可能である。この様な制限方
式は第1θ図について説明する。
第5図は第3図の8数値比較回路26の詳しい回路図で
ある。第5図の回路を使って、容積要素レジスタ21か
らの容積要素の頂点の数値と補間回路35からの部分容
積要素の頂点の数値との両方を、導線24(第3図)の
面の定数と比較する。
部分容積要素の比較は、容積要素の比較によってその容
積要素と交差することが判った後にだけ、行なわれるの
であるから、この両方の作用に1個の比較器を使うこと
が出来る。その代りに、別個の部分容積要素比較器を使
って、他方の比較器が次の交差する容積要素を捜してい
る間に、交差する容積要素の部分容積要素を捜すことが
出来る。
この目的の為には、第5図の比較回路を単に二重にする
第5図の比較回路は8個の並列比較回路60゜61・・
・・・・62で構成される。各々の回路60乃至62は
多重デイジットの容積要素の頂点の数値(1組の導線6
3乃至65)を導線24に出る面の一定の数値と比較し
て、対応する頂点の数値が面の定数より大きい場合、そ
の場合にだけ、8本の出力導線67乃至69の内の1つ
に出力信号を発生することが出来る。こう云う8個の出
力信号が全部0検出器70及び全部1検出器71に同時
に送出される。検出器70及び71からの出力がオア・
ゲート72で組合され、その出力が導線74を介して、
容積要素及び部分容積要素(容積要素の頂点の比較の為
)を比較する間だけ、第3図のアドレス発生器22を前
進させる様に結合されると共に、(部分容積要素の頂点
の比較の為)第3図の容積要素分割器29に結合される
。オア・ゲート72の出力信号がインバータ73にも印
加され、導線75に出るその出力信号を使って、第3図
の法線ベクトル計算器27を付能する(容積要素の頂点
の比較の為)と共に、(部分容積要素の頂点の比較の為
)第3図の法線補間回路31の出力をゲートする。
第6図は第3図の法線ベクトル計算器27の詳しい回路
図である。第2図に見られる様に、法線(又はN)の数
値が、単に容積要素の頂点に於けるデータの数値の勾配
を決定することによって計算される。即ち、直ぐ隣合う
頂点(各々の軸に沿って)に於ける数値を減算して、必
要な勾配の近似値とする。具体的に云うと N (Vl、)−W5−V4 N (V ly ) =W 12  V 2N (Vl
、)−Wl−V5      (1)N (V2.c)
−W6−V3 N (V2y)−Vl−W7 N (V2.)−W2−VB これは第2図のv3乃至v8に対するものである。
法線ベクトルの大きさではなく、その方向だけに関心が
あるから、こうして得られた値は次の様に正規化する。
と゛して N(Vl、c ) N(Vl、)− N(Vly) N(Vl、y)− N(Vl工) N(Vl□)− り こ\でNは、容積要素の頂点に於ける単位の値を持つ法
線ベクトルである。即ち、Nは法線ベクトルの角度方向
にある単位ベクトルである。こう云う法線の単位ベクト
ルの数値が、第3図に示す法線レジスタ28に記憶され
、容積要素の内部の法線の数値を補間する為に後で使わ
れる。第6図の回路は式(1)及び(2)の手順を直接
的に実行する。然し、第6図の回路は容積要素の1つの
頂点の値だけに対する法線を発生し、従って、容積要素
の8個の頂点の法線の数値を同時に発生する為には、第
6図と同様な8個の回路を設けることが必要である。
容積要素の各々の頂点x、  y及び2に対し、X方向
の直ぐ隣接して跨がる頂点の数値を第6図の減算回路8
0に印加する。同様に、隣接して跨がるyの数値を減算
回路81に印加し、隣接して跨がる2の数値を減算回路
82に印加する。減算回路80乃至82で得られた差が
夫々自乗回路83乃至85に夫々印加される。その出力
信号が加算回路87で加算される。加算回路87の出力
信号が平方根回路88に印加され、その出力信号が除算
回路89乃至91に同時に印加される。除算回路89に
対する他方の入力は減算回路80の出力信号を受取る。
除算回路90.91に対する他方の人力信号は夫々減算
回路81.82から来る。
除算回路89乃至91の出力信号が、第3図の法線レジ
スタ28に記憶すべき単位法線ベクトルである。
第7図は第3図の容積要素分割回路29の詳しい回路図
である。−船釣に、第3図の容積要素分割回路は、関心
のある面と交差する容積要素の中にある容積要素の端数
座標を系統的に発生する。
第7図に示す様に、容積要素分割回路はjカウンタ10
0、jカウンタ101及びにカウンタ102を介する。
これらの3つのカウンタは第4図のカウンタ40乃至4
2と同様に動作し、導線106の前進パルスに応答して
、逐次的に部分容積要素を同定する。倍率s、cl  
Sy+  s、Lが夫々比較回路103,104,10
5に印加され、そこで夫々カウンタ100.101,1
02にあるカウントと比較される。符合すると、対応す
るカウンタがリセットされ、次のカウンタを1だけ前進
させる。この為、カウンタエ00乃至102は夫々基数
S。を計数する。前に述べた様に、倍率は、データ空間
の次元並びに表示装置の解像度に応じて、3つの次元全
部で同じであってもよいし、或いは任意の2つの次元で
同じであってもよい。カウンタ100,101.102
の出力信号が夫々除算回路107,108,109に印
加され、対応する1つの倍率s、、Sy、S、もそれに
対して印加される。除算回路107乃至109の出力信
号は次の式で表わす1組の容積要素端数座標Fx、Fy
、Fzである。
Fx−i/S、   i−0,1,−・−、S、CFy
−j/Sy     j−0,1+  ・・・、 Sン
Fz−に/S2  k−0,1,−=、52容積要素補
間器35(第3図)及び法線補間器31(第3図)が容
積要素端数パラメータFx、Fy、Fzを利用して、容
積要素の頂点のデータ数値の間、並びに容積要素の頂点
の法線の数値の間の3次元の補間を行なう。倍率Sx、
Sy、S2が2進級数の値、即ち2 であれば、除算回
路107乃至109に必要な処理は著しく簡単になり、
2進小数点を1桁だけ左に置換えることを必要とするだ
けである。
第8図は1個の容積要素の斜視図であり、これは第9図
の3次元補間回路を理解するのに役立つ。
第8図に示す容積要素は第1図と同じ様に頂点Vl乃至
v8を持っている。X方向の容積要素端数座標Fxが与
えられると、容積要素の辺の点A1乃至A4が、平面V
l−V2−V5−V6と平面V3−V4−V7−V8の
間のこの端数距離の所に容易に突止められる。同様に、
線B1−82は、線Al−A3及び線A2−A4の間に
y方向の容積要素の端数距MFyの所に突止めることが
出来る。最後に、2方向の容積要素の端数距離Fzが与
えられ\ば、線B1−82に沿って補間された数値V 
 を突止めることが出来る。第9図に示yz す補間回路35はこの3次元補間方式を使って、関心の
ある面と交差する容積要素内にある部分容積要素の頂点
に於けるデータの数値を決定する。
第9図の補間回路は、第3図の説明で述べた様に、関心
のある面と交差すると判った部分容積要素の頂点に於け
る単位法線ベクトルの補間にも用いられる。
第9図の補間回路35は第8図に示した3次元補間を直
接的に実行する。補間回路35によって行なわれる第8
図に示す容積要素の補間は次の様に要約することが出来
る。
AI−V3+Fx (Vl−V3) A2−V4+Fx (V2−V4) A3=V7+Fx (V5−V7) A4−V8+Fx (V6−V8)     (4)B
l−A1+Fy (A2−AI) B2−A3+Fy (A4−A3) V Ijk−B 2 + F z (B I  B 2
 )補間回路35が7つの減算回路112乃至118.
7つの加算回路119乃至125及び6つの乗算回路1
26乃至132を、この式(4)を実現する様な容易に
考えられる形で一緒に接続して構成される。容積要素の
3つの端数値Fx、  Fy、  F2と容積要素の8
個の頂点の数値Vl乃至v8とが、第7図に示す様な容
積要素分割器29から回路35に印加され、部分容積要
素の頂点の数値vIjkを取出す。面と交差するかどう
か、部分容積要素を試験する為には、部分容積要素の8
個の頂点の数値が必要であるから、何れも第9図に示す
様な形式の8個の捕間回路35を設け、各々の回路35
は第7図に示すのと同様であるが、同じ容積要素の別の
頂点を同定する様に、互いにi。
j及びにカウントがずれる相異なる容積要素分割器によ
って駆動される。この代りに、1組の部分容積要素カウ
ンタを使い、ずれは第4図に示すのと同様な増数回路及
び減数回路によって作ってもよい。
m9図の補間回路から得られた補間による部分容積要素
の単位法線ベクトルが、その法線ベクトルを持つ面と構
成した部分容積要素を同定する、容積要素分割器29か
らの1組の座標と同時に、表示装置34(第3図)に送
出される。前に提案シタ様ニ、倍率sx 、Sy 、 
 Sz h’2進法(1) ’l ”の値であれば、部
分容積要素のアドレスは、単に第4図のレジスタ45に
ある容積要素アドレスを最上位デイジットとし、それと
並べて第7図のカウンタ100乃至102にある部分容
積要素アドレスを下位デイジットして並べたものである
。面の単位法線ベクトルの座標は、交差した部分容積要
素の正確に幾何学的な中心にあると想定するのが最も正
確であるが、全ての単位法線ベクトルが一定量だけずれ
ても、その結果得られる像は何等変ることがないから、
既に利用し得る部分容積要素の頂点の座標のどれを用い
ても構わない。
第10図は第3図の表示装置34の詳しいブロック図で
ある。第10図に示す表示装置は、表示リスト記憶回路
140を持ち、その中に関心のある面と交差する全ての
部分容積要素のx−y−z座標と、この各々の座標に関
連する単位法線ベクトルの成分Nx−Ny−Nzとが記
憶される。全ての座標と法線の対を取出すのに要する時
間が、表示の為にこう云う対を処理するのに要する時間
とは無関係であるから、表示リスト140は、座標と法
線の対に対し、先入れ先出しくFIFO)バッファ記憶
装置として作用する。最初の座標と法線の対がリスト1
40に利用出来る時、単位法線ベクトルの成分が角度比
較器142に供給され、この比較器には配向発生器14
1からの観察角度も供給される。比較器142が法線ベ
クトルと観察方向との間の内積に比例する信号を発生し
、この信号を導線143を介して画素強度又は画素色信
号として出力する。
単位法線ベクトルが比較器142に供給されるのと同時
に、交点のx、y及び2座標がデータ空間切取り装置1
44に供給される。データ切取り選択器145が、表示
に使うべきデータ空間の容積を定める6つの切断平面を
選ぶ。関心のある容積はデータ空間全体よりも小さい場
合が大抵であるから、処理されるデータ空間の寸法を制
限することによって、かなりの処理時間を節約すること
が出来る。最も簡単な形では、データ切取り選択器14
5がx、y、z軸上の座標区間を選択し、データ空間切
取り装置144は、この区間に入らない全ての交点を除
外する。これは、全てのデータ座標を全ての切取り座標
と並列に比較することによって、非常に速やかに行なう
ことが、出来る。
データ空間の切取り平面を慎重に選ぶことにより、関心
のある面と交差させることが出来、こうして関心のある
面を通る可視的な断面を作ることが出来る。
3次元データ空間に入っている3次元の面上の面の点が
与えられた場合、こう云う面の点を、2次元表示装置に
表示する為の表示点に変換することが必要である。これ
は、マトリクスの回転及び並進により、データ点を表示
平面に投影することによって達成される。具体的に云う
と、面の点のX、y及び2成分を周知の下記のマトリク
ス操作により、ベクトルxOr  yo l  zQに
よって表示点xd r  Yd +  zdに回転並び
に並進させることが出来る。
X dMMaM 11X + M 12 V + M 
1a Z + X 。
yd−M  X+M22Y+M23Z+1゜I Y d −Mat X + MB2 Y + Mg2 
Z + Z 。
第10図の装置では、式(5)の考えられる全てのマト
リクス係数M1jがマトリクス記憶装置146に記憶さ
れており、所望の表示の配向に対して使うべき特定の係
数が配向発生器111によって選択される。
回転回路147の出力は一連の表示可能な座標xd +
  yd +  zdである。こう云う表示可能な点は
、観察切取り選択器148の制御のもとに、観察空間切
11り回路149によって再び減少することが出来る。
データ切取り選択器145の場合と同じく、選択器14
8によって6つの切取り平面が選択されて、表示される
図は表示データの選ばれた部分、即ち表示窓に制限され
る。その結果、切取り装置149からの座標と線143
を介して送られる比較器142からの画素強度(色)を
含めて、二重に減少した一連の表示可能な点が得られる
。然し、関心のある面は、物理的な空間の面及び選ばれ
た観察角度の畳込み積分をする為に、同じX、及びyd
座標の所で、2つ以上の点が表示平面に回転することが
あることが認められよう。
面を正確に表示する為、観察者に最も近い面の点を、こ
の一番近い面の点より遠くにある(その為にそれによっ
て隠れる)面の点に優先して選ばなければならない。こ
の曖昧さを解決する為、奥行バッファ記憶装置152が
設けられている。
動作につい、て説明すると、奥行バッファ152は、バ
ッファ152の全ての記憶アドレスに考えられる最大の
奥行をロードすることにより、初期設定される。その後
、観察空間切次り装置149からのXd及びy、座標が
アドレス・レジスタ150に印加される。レジスタ15
0のアドレスを使って、奥行バッファ152及びフレー
ム・バッファ154の両方をアドレスする。フレーム・
バッファ154は、表示内の表示可能な各々の点(画素
)に対する表示強度(又は色)情報を、表示装置の対応
する点の座標と対応する記憶アドレスの所に持っている
。切取り装置149からの奥行座標z、1がアンド・ゲ
ート153を介して奥行バッファ152に印加されると
共に、奥行比較器151にも印加される。奥行バッファ
152は、フレーム・バッファ154と同じく、各々の
表示可能な点に対応するアドレスに情報を持っているが
、強度情報ではなく奥行情報を持っている。前に述べた
様に、全てのアドレスに於ける初期値は考えられる最大
の奥行の値、例えば、データ空間の最大寸法である。比
較器151は、前に記憶されたzdの値と同じXd及び
Yd座標の所にある新しい2.の値が、前に記憶された
zdの値より小さい(即ち、表示平面に一層接近してい
る)かどうかを決定する。そうであれば、比較器151
の出力信号によってアンド・ゲート153を付能するこ
とにより、新しい奥行の値が奥行バッファ152に記憶
される。同時に、比較器142からの新しい強度データ
により、アンド・ゲート155を付能することによって
、フレーム・バッファ154の古い強度データが書換え
られる。こうすることにより、表示リスト140の全て
の表示点を使いきった後、フレーム拳バッファ154の
内容は、観察平面に最も接近した表示点だけで構成され
る。フレーム・バッファ154の内容を使って、表示装
置(陰極線間)156の表示される像を更新する。奥行
バッファ152を最大の奥行の値に初期設定することに
より、切取り装置149から受取った最初の奥行の値が
、比較器151を付能することによって、最大の奥行を
書換えることに注意されたい。
第10図では、データ空間切取り選択器145、データ
空間切取り装置144、配向発生器141、角度比較器
142、マトリクス選択器148及びアドレス・レジス
タ150を含む変換が、ビット直列形シリコン・コンパ
イラを使ってVASIの1.25ミクロンのCMOSチ
ップを設計することによって構成された。こうして得ら
れたチップが毎秒937,500個の点の速度で点を変
換して、対話形の3次元表示を行なう。第10図に従っ
て印刷配線板を構成して、サン・ワークステーションを
使って、医学的な物体の品質の高い3次元像が得られた
上に述べた密実な物体の内部構造の表示プロセッサが高
度に並列のパイプライン・アーキテクチュアであること
により、可視像を非常に速やかに発生することが出来る
ことが理解されよう。更に、種々の並列処理回路は大規
模集積回路チップで容易に製造することが出来るので、
この処理回路の寸法及びコストが下げられると同時に、
こう云う回路の動作速度及び信頼性が改善される。
第3図の装置によって実施される方法は、第11図のフ
ローチャートによって表わすことが出来る。開始ボック
ス160から始まって、物理的な数値の3次元配列が高
度に並列のバルク・メモリに記憶され、そこから−度に
1つの容積要素ずつ、並列にアクセスすることが出来る
。ボックス162で、次の容積要素の8個の頂点のデー
タの数値が、バルク・メモリからアクセスされ、判定ボ
ックス163に送られる。関心のある面がこの容積要素
と交差しなければ、再びボックス162に入って、次の
容積要素を検索する。関心のある面がこの容積要素と交
差すれば、同時にボックス164及び166に入る。ボ
ックス164では、近い隣りの容積要素をアクセスし、
隣りの容積要素の頂点の数値をボックス165で使って
、容積要素の各々の頂点に対する単位法線ベクトルを決
定する。ボックス166で、交差した容積要素を分割し
、ボックス167で、分割された容積要素の次の部分容
積要素に対する8個の頂点の数値を、ボックス162で
求められた容積要素の頂点の数値から補間する。判定ボ
ックス168で、関心のある面と交差するかどうか、部
分容積要素の頂点の数値を試験する。交差しなければ、
ボックス、167に再び入って、次の部分容積要素に対
する頂点の数値を決定する。判定ブロック168で部分
容積要素が交差することが検出されると、部分容積要素
の座標を使って、ボックス169で、部分容積要素の座
標に於ける単位法線ベクトルを補間する。ボックス17
0で、部分容積要素の座標及び部分容積要素の単位法線
ベクトルの対を記憶する。
座標がボックス171で処理されて観察の向きを選択し
、データ空間及び観察空間に於けるデータを減少し、減
少した座標を選ばれた観察平面へ回転並びに並進させる
。面の単位法線ベクトルがボックス172で処理されて
、法線ベクトルの角度を観察角度を比較することにより
、色又は強度情報を取出す。ボックス173で、こうし
て得られた画素の色又は画素の強度がボックス171に
よって定められた表示座標の所で表示される。
この発明のこの他の実施例も、この発明の範囲内で当業
者であれば容易に行なうことが出来ることは明らかであ
ろう。
【図面の簡単な説明】
第1図は非侵入形手段によって求められた物理的な数値
の配列内にある8個の格子位置によって限定された1個
の容積要素の斜視図であって、この発明の詳細な説明す
るのに役立つ図、第2図は容積要素の近隣の斜視図であ
って、容積要素の容積を構成する8個の格子位置に対す
る法線の数値を計算するのに必要な近隣の数値を例示す
ると共に、この発明の詳細な説明するのにも役立つ図、
第3図は実時間で3次元の面の像を発生するこの発明の
処理装置の全体的なブロック図、第4図は第3図の処理
装置に役立つ容積要素アドレス発生器の詳しい回路図、
第5図は第3図の処理装置に役立つ8個の数値の比較回
路の詳しい回路図、第6図は第3図の処理装置に役立つ
面の法線の数値を決定する法線ベクトル計算器の詳しい
回路図、第7図は第3図の処理装置に役立つ容積要素分
割回路の詳しいブロック図、第8図は容積要素の斜視図
で、第9図の回路を理解するのに役立つ3次元−次補間
過程を示す図、第9図は第3図の処理装置に役立つ3次
元−次補間回路の詳しいブロック図、第10図は第3図
の処理装置に役立つ表示プロセッサの詳しいブロック図
、第11図はこの発明の装置で行なわれる方法のフロー
チャートである。 [主な符号の説明] vl乃至V8:容積要素の頂点(及びその数値)Wl乃
至W24:近隣の容積要素の頂点(及びその数値) 20:バルク・メモリ 22ニアドレス発生器 26二8数値比較回路 27:法線計算器 29:容積要素分割器 34:表示装置

Claims (1)

  1. 【特許請求の範囲】 1、3次元の密実な物体の内面の2次元像を表示する装
    置に於て、3次元物体内の容積要素を限定する規則的に
    相隔たる格子位置に於ける該物体に関連した少なくとも
    1種類の物理的な性質の数値を表わす3次元信号パター
    ンを記憶する手段と、該記憶する手段から、前記信号パ
    ターンによって表わされる1個の容積要素の数値全体を
    一度に選択的にアクセスする手段と、該容積要素が前記
    3次元物体の予め選ばれた内面と交差するかどうか、頂
    点の数値によって限定される前記容積要素の数値を試験
    する並列回路手段と、交差する容積要素だけを部分容積
    要素に分割する並列回路手段と、各々の部分容積要素が
    前記予め選ばれた面と交差するかどうか、各々の部分容
    積要素の数値を試験する並列回路手段と、前記予め選ば
    れた面と交差する部分容積要素に対する部分容積要素の
    法線の数値を決定する並列回路手段と、部分容積要素の
    座標及び部分容積要素の法線の数値を受取って、該法線
    の数値をスクリーンに表示する表示プロセッサ手段とを
    有し、表示される数値が前記予め選ばれた面の像を表わ
    す様にした装置。 2、アクセスされた容積要素全体と共通の面を持つ容積
    要素に対する複数個の近隣の容積要素の数値を同時にア
    クセスする手段と、該近隣の容積要素の数値を利用して
    、中心の差により、面の法線ベクトルを決定する並列回
    路手段とを更に有する請求項1記載の装置。 3、前記容積要素の数値を試験する並列回路手段が、全
    ての頂点の数値を前記予め選ばれた面を表わす一定の数
    値と同時に比較する手段と、該一定の数値が、前記頂点
    の数値の内の最大及び最小の間にある時を検出する手段
    とで構成されている請求項1記載の装置。 4、前記交差した容積要素だけを部分容積要素に分割す
    る並列回路手段が、容積要素倍率を利用して、容積要素
    内の部分容積要素端数アドレスを発生する並列回路手段
    を有し、該倍率が規則的に相隔たる格子の各々の次元に
    対して異なっている請求項1記載の装置。 5、前記交差した容積要素だけを部分容積要素に分割す
    る並列回路手段が、容積要素倍率を利用して容積要素内
    の部分容積要素端数アドレスを発生する並列回路手段を
    有し、該倍率が規則的に相隔たる格子の少なくとも2つ
    の次元に対して同じである請求項1記載の装置。 6、前記倍率の内の少なくとも1つが、nを整数として
    、2^nで表わされる請求項4記載の装置。 7、前記倍率の内の少なくとも1つが、nを整数として
    、2^nで表わされる請求項5記載の装置。 8、前記表示プロセッサ手段が、部分容積要素の法線の
    数値を観察角度と比較して、前記予め選ばれた面の像に
    対する画素の数値を発生する手段と、前記法線の数値の
    座標を、前記予め選ばれた面の像に対する観察平面上の
    座標に変換する手段と、前記予め選ばれた面の像に含ま
    れる座標を制限する手段と、同じ画素位置を占める表示
    面の画素の曖昧さを解決する手段とを有する請求項1記
    載の装置。 9、前記予め選ばれた面の像に含まれる座標を制限する
    手段が、前記信号パターンの空間に於ける部分容積要素
    の法線の数値を減少する手段で構成されている請求項8
    記載の装置。 10、前記予め選ばれた面の像に含まれる座標を制限す
    る手段が、前記観察平面の空間に於ける部分容積要素の
    法線の数値を減少する手段で構成されている請求項8記
    載の装置。 11、3次元の密実な物体の内面の2次元像を表示する
    方法に於て、3次元物体内の容積要素を限定する規則的
    に相隔たる格子位置に於ける該物体に伴う少なくとも1
    種類の物理的な性質の数値を表わす3次元信号パターン
    を記憶し、該記憶する手段から、信号パターンによって
    表わされる容積要素全体の数値を一度に選択的にアクセ
    スし、該容積要素が前記3次元物体の予め選ばれた内面
    と交差するかどうか、前記頂点の数値によって限定され
    た複数個の容積要素の数値を並列に試験し、交差する容
    積要素だけを複数個の部分容積要素に並列に分割し、各
    々の部分容積要素が前記予め選ばれた面と交差するかど
    うか、各々の部分容積要素の数値を同時に且つ並列に試
    験し、前記予め選ばれた面と交差する部分容積要素に対
    し、部分容積要素の法線の数値を並列に決定し、表示さ
    れる数値が前記予め選ばれた面の像を表わす様に、前記
    法線の数値を表示する工程を含む方法。 12、アクセスされた容積全体と共通の面を持つ容積要
    素に対する複数個の近隣の容積要素の数値を同時にアク
    セスし、該近隣の容積要素の数値を利用して、中心の差
    により、全ての容積要素の頂点の面の法線ベクトルを並
    列に決定する工程を含む請求項11記載の方法。 13、前記試験する工程が、全ての頂点の数値を予め選
    ばれた面を表わす一定の数値と同時に比較し、該一定の
    数値が頂点の数値の最大及び最小の間にある時を決定す
    ることを含む請求項11記載の方法。 14、前記分割する工程が、容積要素の倍率を利用して
    、容積要素内の複数個の部分容積要素の端数アドレスを
    並列に発生することを含み、該倍率が規則的に相隔たる
    格子の各々の次元に対して異なっている請求項11記載
    の方法。 15、前記分割する工程が、容積要素の倍率を利用して
    容積要素内の複数個の部分容積要素の端数アドレスを並
    列に発生することを含み、該倍率は規則的に相隔たる格
    子の少なくとも2つの次元に対して同じである請求項1
    1記載の方法。 16、前記倍率の内の少なくとも1つが、nを整数とし
    て、2^nで表わされる請求項14記載の方法。 17、前記倍率の内の少なくとも1つが、nを整数とし
    て、2^nで表わされる請求項15記載の方法。 18、前記表示する工程が、前記部分容積要素の法線の
    数値を観察方向と比較して、前記予め選ばれた面の像に
    対する画素値を発生し、前記面の点の座標を前記予め選
    ばれた面の像に対する観察平面上の座標に変換し、該予
    め選ばれた面の像に含まれる座標を制限し、同じ画素位
    置を占める表示面の画素の間の曖昧さを解決する工程を
    含む請求項11記載の方法。 19、前記制限する工程が、前記信号パターンの空間に
    於ける前記部分容積要素の法線の数値を減少することを
    含む請求項18記載の方法。 20、前記制限する工程が、前記観察平面の空間に於け
    る部分容積要素の法線の数値を減少することを含む請求
    項18記載の方法。 21、3次元データ配列の2次元像を速やかに表示する
    装置に於て、該配列にあるデータ数値に応答して、前記
    データ配列内の関心のある面上の複数個の点の位置を同
    時に突止める第1のパイプライン回路手段と、やはり前
    記データの数値に応答して、前記複数個の点に於ける表
    示される像の強度を同時に設定する第2のパイプライン
    回路手段とを有する装置。 22、前記第1及び第2のパイプライン回路手段が大規
    模集積回路チップで構成される請求項21記載の装置。 23、3次元データ配列の2次元像を速やかに表示する
    方法に於て、前記配列内の複数個のデータ数値に応答し
    て、前記データ配列内にある関心のある面上の複数個の
    点の位置を突止め、前記データの数値に応答して、前記
    複数個の点に於ける像の強度を同時に設定する工程を含
    む方法。
JP1295187A 1988-11-22 1989-11-15 物体内部構造表面の2次元像を表示する装置と方法 Expired - Fee Related JP2744490B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US07/275,154 US4985834A (en) 1988-11-22 1988-11-22 System and method employing pipelined parallel circuit architecture for displaying surface structures of the interior region of a solid body
US275,154 1988-11-22

Publications (2)

Publication Number Publication Date
JPH02211587A true JPH02211587A (ja) 1990-08-22
JP2744490B2 JP2744490B2 (ja) 1998-04-28

Family

ID=23051088

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1295187A Expired - Fee Related JP2744490B2 (ja) 1988-11-22 1989-11-15 物体内部構造表面の2次元像を表示する装置と方法

Country Status (2)

Country Link
US (1) US4985834A (ja)
JP (1) JP2744490B2 (ja)

Families Citing this family (52)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
FR2656129B1 (fr) * 1989-12-20 1992-03-13 Gen Electric Cgr Procede de reconstruction multi-echelle de l'image de la structure d'un corps.
JPH0793926B2 (ja) * 1991-03-01 1995-10-11 ゼネラル・エレクトリック・カンパニイ ディスプレイリスト表面データを使用する手術シミュレーション装置
DE4117117A1 (de) * 1991-05-25 1992-11-26 Hoehne Karl Heinz Prof Dr Dreidimensionale darstellung von raeumlichen strukturen
US5734384A (en) * 1991-11-29 1998-03-31 Picker International, Inc. Cross-referenced sectioning and reprojection of diagnostic image volumes
EP0549189B1 (en) * 1991-12-23 2003-05-07 General Electric Company Solid model generation by span method using dividing cubes
EP0549185A2 (en) * 1991-12-23 1993-06-30 General Electric Company System for 3D scan conversion of a polygonal model into a point and normal format, displayed utilizing an accelerator circuit
EP0549182A2 (en) * 1991-12-23 1993-06-30 General Electric Company Apparatus and method for displaying surgical cuts in three-dimensional models
DE69233060T2 (de) * 1991-12-23 2004-04-08 General Electric Co. System zur Anzeige von räumlichen Einschnitten für Festkörpermodellflächen
US5303706A (en) * 1992-05-21 1994-04-19 North American Philips Corporation Directional interpolation for magnetic resonance angiography
IL102289A (en) * 1992-06-24 1997-08-14 R Technologies Ltd B V Method and system for processing moving images
US5295488A (en) * 1992-08-05 1994-03-22 General Electric Company Method and apparatus for projecting diagnostic images from volumed diagnostic data
US5361385A (en) * 1992-08-26 1994-11-01 Reuven Bakalash Parallel computing system for volumetric modeling, data processing and visualization
EP0646263B1 (en) * 1993-04-20 2000-05-31 General Electric Company Computer graphic and live video system for enhancing visualisation of body structures during surgery
US5368033A (en) * 1993-04-20 1994-11-29 North American Philips Corporation Magnetic resonance angiography method and apparatus employing an integration projection
US5412563A (en) * 1993-09-16 1995-05-02 General Electric Company Gradient image segmentation method
US5782762A (en) * 1994-10-27 1998-07-21 Wake Forest University Method and system for producing interactive, three-dimensional renderings of selected body organs having hollow lumens to enable simulated movement through the lumen
US5920319A (en) * 1994-10-27 1999-07-06 Wake Forest University Automatic analysis in virtual endoscopy
US6694163B1 (en) 1994-10-27 2004-02-17 Wake Forest University Health Sciences Method and system for producing interactive, three-dimensional renderings of selected body organs having hollow lumens to enable simulated movement through the lumen
JPH08138078A (ja) * 1994-11-09 1996-05-31 Toshiba Medical Eng Co Ltd 画像処理装置
US5485842A (en) * 1994-11-30 1996-01-23 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic scan conversion for three dimensional display processing
US5659629A (en) * 1995-04-07 1997-08-19 General Electric Company Visualization of a multi-dimensional data set acquired with a surface receive coil
US5602891A (en) * 1995-11-13 1997-02-11 Beth Israel Imaging apparatus and method with compensation for object motion
US5900880A (en) * 1996-03-06 1999-05-04 General Electric Company 3-D surfaces generated from a list of cubic elements
US6046746A (en) * 1996-07-01 2000-04-04 Sun Microsystems, Inc. Method and apparatus implementing high resolution rendition of Z-buffered primitives
US6115047A (en) * 1996-07-01 2000-09-05 Sun Microsystems, Inc. Method and apparatus for implementing efficient floating point Z-buffering
US7486811B2 (en) * 1996-09-16 2009-02-03 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination of objects, such as internal organs
US7194117B2 (en) 1999-06-29 2007-03-20 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination of objects, such as internal organs
US6331116B1 (en) * 1996-09-16 2001-12-18 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual segmentation and examination
US8682045B2 (en) * 1997-02-25 2014-03-25 Wake Forest University Health Sciences Virtual endoscopy with improved image segmentation and lesion detection
US6246784B1 (en) 1997-08-19 2001-06-12 The United States Of America As Represented By The Department Of Health And Human Services Method for segmenting medical images and detecting surface anomalies in anatomical structures
US6249286B1 (en) * 1998-10-31 2001-06-19 Hewlett-Packard Company Memory efficient surface normal compression
EP1133257A4 (en) 1998-11-25 2007-12-26 Univ Wake Forest IMPROVED VIRTUAL ENDOSCOPY SYSTEM WITH IMAGE SEGMENTATION AND DETECTION OF LESIONS
US6463168B1 (en) * 1998-12-04 2002-10-08 General Electric Company Apparatus and method for rapid connectivity processing of images
US6753878B1 (en) * 1999-03-08 2004-06-22 Hewlett-Packard Development Company, L.P. Parallel pipelined merge engines
US6385604B1 (en) * 1999-08-04 2002-05-07 Hyperroll, Israel Limited Relational database management system having integrated non-relational multi-dimensional data store of aggregated data elements
US6408292B1 (en) * 1999-08-04 2002-06-18 Hyperroll, Israel, Ltd. Method of and system for managing multi-dimensional databases using modular-arithmetic based address data mapping processes on integer-encoded business dimensions
US20020029207A1 (en) * 2000-02-28 2002-03-07 Hyperroll, Inc. Data aggregation server for managing a multi-dimensional database and database management system having data aggregation server integrated therein
WO2001093745A2 (en) * 2000-06-06 2001-12-13 The Research Foundation Of State University Of New York Computer aided visualization, fusion and treatment planning
WO2002029723A1 (en) 2000-10-02 2002-04-11 The Research Foundation Of State University Of Newyork Enhanced virtual navigation and examination
US7630750B2 (en) * 2001-02-05 2009-12-08 The Research Foundation For The State University Of New York Computer aided treatment planning
US20050154563A1 (en) * 2001-08-27 2005-07-14 Ulf Hassler Device and method for evaluating a characteristic of an object
DE10141807A1 (de) * 2001-08-27 2003-03-27 Fraunhofer Ges Forschung Vorrichtung und Verfahren zum Auswerten einer Beschaffenheit eines Objekts
US7596256B1 (en) 2001-09-14 2009-09-29 The Research Foundation For The State University Of New York Computer assisted detection of lesions in volumetric medical images
US7324104B1 (en) 2001-09-14 2008-01-29 The Research Foundation Of State University Of New York Method of centerline generation in virtual objects
US7260250B2 (en) * 2002-09-30 2007-08-21 The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services Computer-aided classification of anomalies in anatomical structures
TWI268148B (en) * 2004-11-25 2006-12-11 Univ Chung Yuan Christian Image analysis method for vertebral disease which comprises 3D reconstruction method and characteristic identification method of unaligned transversal slices
US20100260390A1 (en) * 2005-11-30 2010-10-14 The Research Foundation Of State University Of New York System and method for reduction of false positives during computer aided polyp detection
US8452061B2 (en) * 2005-11-30 2013-05-28 The Research Foundation Of State University Of New York Electronic colon cleansing method for virtual colonoscopy
US8243334B2 (en) 2008-06-06 2012-08-14 Virginia Venture Industries, Llc Methods and apparatuses for printing three dimensional images
EP2309450B1 (en) * 2009-09-21 2016-04-06 Orcamp AB Method for characterizing a blood vessel
US10952795B2 (en) * 2017-06-30 2021-03-23 Biosense Webster (Israel) Ltd. System and method for glass state view in real-time three-dimensional (3D) cardiac imaging

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6237782A (ja) * 1985-06-05 1987-02-18 ゼネラル・エレクトリツク・カンパニイ 3次元面構造を表示するための装置と方法
JPS63145638A (ja) * 1986-12-10 1988-06-17 株式会社日立製作所 三次元画像処理方式
JPH0275083A (ja) * 1988-09-12 1990-03-14 Nippon Yougiyoushi Kenkyusho:Kk 考古学出土品の輪郭描画装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4622634A (en) * 1983-03-18 1986-11-11 Irex Corporation Parallel processing of simultaneous ultrasound vectors
US4594673A (en) * 1983-06-28 1986-06-10 Gti Corporation Hidden surface processor
JPS6075033A (ja) * 1983-09-30 1985-04-27 株式会社東芝 X線診断装置
IL70213A (en) * 1983-11-13 1988-02-29 Paul Fenster Digital fluorographic image enhancement system
JPS60152942A (ja) * 1984-01-23 1985-08-12 Toshiba Corp Nmr―ctスキャン計画装置
US4709332A (en) * 1984-10-05 1987-11-24 North Shore University Hospital High speed image data processing
US4729098A (en) * 1985-06-05 1988-03-01 General Electric Company System and method employing nonlinear interpolation for the display of surface structures contained within the interior region of a solid body
US4710876A (en) * 1985-06-05 1987-12-01 General Electric Company System and method for the display of surface structures contained within the interior region of a solid body
US4850027A (en) * 1985-07-26 1989-07-18 International Business Machines Corporation Configurable parallel pipeline image processing system
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
US4751643A (en) * 1986-08-04 1988-06-14 General Electric Company Method and apparatus for determining connected substructures within a body
US4821213A (en) * 1986-12-19 1989-04-11 General Electric Co. System for the simultaneous display of two or more internal surfaces within a solid object
US4821210A (en) * 1987-04-02 1989-04-11 General Electric Co. Fast display of three-dimensional images
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6237782A (ja) * 1985-06-05 1987-02-18 ゼネラル・エレクトリツク・カンパニイ 3次元面構造を表示するための装置と方法
JPS63145638A (ja) * 1986-12-10 1988-06-17 株式会社日立製作所 三次元画像処理方式
JPH0275083A (ja) * 1988-09-12 1990-03-14 Nippon Yougiyoushi Kenkyusho:Kk 考古学出土品の輪郭描画装置

Also Published As

Publication number Publication date
US4985834A (en) 1991-01-15
JP2744490B2 (ja) 1998-04-28

Similar Documents

Publication Publication Date Title
JP2744490B2 (ja) 物体内部構造表面の2次元像を表示する装置と方法
US4984157A (en) System and method for displaying oblique planar cross sections of a solid body using tri-linear interpolation to determine pixel position dataes
US4737921A (en) Three dimensional medical image display system
US5166876A (en) System and method for detecting internal structures contained within the interior region of a solid object
US5079699A (en) Quick three-dimensional display
US4989142A (en) Three-dimensional images obtained from tomographic slices with gantry tilt
US5898793A (en) System and method for surface rendering of internal structures within the interior of a solid object
EP0968683B1 (en) Method and apparatus for forming and displaying image from a plurality of sectional images
CA1315902C (en) Minimization of directed points generated in three-dimensional dividing cubes method
JPS6297074A (ja) 3次元の面構造を表示する装置
Goldwasser et al. Real-time display and manipulation of 3-D medical objects: The voxel processor architecture
US4953087A (en) Three-dimensional images obtained from tomographic data having unequally spaced slices
JPH01166267A (ja) 画像化方法および装置
JPH09182747A (ja) 医療技術画像形成装置における画像形成方法
EP0318293B1 (en) Apparatus and method for processing tomographic data
EP0373854B1 (en) Apparatus and method for detecting internal structures contained within the interior region of a solid object
US6191789B1 (en) Ray casting method using hardware
EP0318291B1 (en) Apparatus and method for generating images from tomographic data
Barrett et al. A low-cost PC-based image workstation for dynamic interactive display of three-dimensional anatomy
Flynn et al. Requirements for the display and analysis of three-dimensional medical image data
Delegacz et al. Three-dimensional visualization system as an aid for lung cancer detection
Chang et al. Visualizing the anatomical-functional correlation of the human brain

Legal Events

Date Code Title Description
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: 20090206

Year of fee payment: 11

LAPS Cancellation because of no payment of annual fees