JP2004510515A - 器官のような物体の3次元仮想検査を行うシステムおよび方法 - Google Patents
器官のような物体の3次元仮想検査を行うシステムおよび方法 Download PDFInfo
- Publication number
- JP2004510515A JP2004510515A JP2002533259A JP2002533259A JP2004510515A JP 2004510515 A JP2004510515 A JP 2004510515A JP 2002533259 A JP2002533259 A JP 2002533259A JP 2002533259 A JP2002533259 A JP 2002533259A JP 2004510515 A JP2004510515 A JP 2004510515A
- Authority
- JP
- Japan
- Prior art keywords
- virtual
- scan
- voxels
- image
- voxel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B23/00—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
- G09B23/28—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30028—Colon; Small intestine
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30084—Kidney; Renal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30172—Centreline of tubular or elongated structure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/028—Multiple view windows (top-side-front-sagittal-orthogonal)
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/20—Indexing scheme for editing of 3D models
- G06T2219/2004—Aligning objects, relative positioning of parts
Abstract
Description
【発明の属する技術分野】
本発明は、体積を基礎とする3次元仮想検査を行うシステムおよび方法に関し、より特には、増強された視覚化およびナビゲーション特性を与えるシステムに関する。
【0002】
【従来の技術】
結腸癌は、世界中で主な死因であり続けている。初期において人間の結腸にポリープとして現れる癌の成長の初期の検出は、患者の回復の可能性を大幅に改善することができる。現在、患者の結腸においてポリープまたは他の塊を検出する2つの慣例的な方法がある。第1の方法は、大腸鏡検査であり、大腸鏡と呼ばれる可撓性光ファイバチューブを使用し、結腸を、前記大腸鏡を直腸に物理的に入れることによって、視覚的に検査する。医師は、前記チューブを操作し、結腸におけるなにか異常な成長を探索することができる。前記大腸鏡検査は、信頼できるが、費用および時間の双方が比較的かかり、患者にとって侵略的で不快な痛い手続きである。
【0003】
第2の検出技術は、バリウム浣腸および結腸の2次元X線画像化である。前記バリウム浣腸を使用し、結腸をバリウムで覆い、2次元X線画像を撮って、結腸の画像を獲得する。しかしながら、バリウム浣腸は、必ずしも結腸全体の図を与えるとは限らず、前記動作を行うときにしばしばオペレータ依存であり、患者を過度の放射線に曝し、大腸鏡ほど感度がよくはない。上述した慣例的な手段における不完全さのため、結腸をポリープに関して検査する、より確実で、それほど侵入的でなく、それほど費用がかからない方法が望ましい。肺のような人間の他の器官を、塊に関して、確実で費用効果的な方法において検査する方法も望ましい。
【0004】
米国における癌死の別の主要な原因は、膀胱癌である。1995年において、50000の膀胱癌の新しいケースが報告され、この病気の結果として11000人の死が報告された。膀胱癌に関する最も一般的な試験は、尿“ディップスティック”または慣例的な尿検査の使用である。しかしながら、これらのような試験は、膀胱癌の検出において、そのより後期の発達段階においてのみ有効であり、癌の成長のサイズまたは位置に関するどのような情報も与えない。現在、膀胱の異常を調べる主な方法である膀胱鏡検査は、正確な結果を与え、任意の以上の関連する細部および位置に関する情報を与えることができる。しかしながら、膀胱鏡検査は、医師に限られた視野しか与えず、サイズの客観的な表示が不足している。加えて、膀胱鏡検査は、厳しい尿道の狭窄あるいは活発な膀胱の出血を有する患者には望ましくない。したがって、特に癌の進行の初期段階において、膀胱癌に関して患者を選別する代わりの手順を開発することが望ましい。
【0005】
コンピュータ断層撮影およびMRI(磁気共鳴撮影画像)のような現在利用可能な医療画像化装置を用いる人間の器官の2次元(“2D”)視覚化は、患者の診断に広く使用されてきた。3次元画像を、操作機械から発生された2次元ピクチャ間に堆積または補間することによって形成することができる。器官の画像化と、3次元空間におけるその体積の視覚化とは、その物理的な侵入がないことと、データ操作の容易さとにより、有利である。しかしながら、3次元体積画像の調査は、内部から器官を仮想的に見る利点を完全に利用するために、正確に行わなければならない。
【0006】
環境の3次元(“3D”)体積仮想画像を見るときに、機能的モデルを使用し、仮想空間を調査しなければならない。ある可能なモデルは、ビューワに関して仮想空間を調査する基準点として使用することができる仮想カメラである。一般的な3D仮想環境内のナビゲーションの文脈におけるカメラ制御は、以前に研究されている。仮想空間のナビゲーションに提案された2つの慣例的な形式のカメラ制御がある。第1のものは、オペレータがカメラを異なった位置および向きにおいて操作し、所望のビューを達成する、カメラのオペレータ完全制御である。オペレータは、前記カメラを実際に操縦する。これは、オペレータが、興味のある特定の区分を調査し、他の区分を無視することを可能にする。しかしながら、大きい領域におけるカメラの完全な制御は、面倒で骨が折れ、オペレータは、調査の開始点と終了点との間のすべての重要な特徴を見るとは限らない。
【0007】
カメラ制御の第2の技術は、計画ナビゲーション方法であり、この方法は、カメラに、とるべき予め決められた経路を割り当て、オペレータによって変更することができない。これは、予定された“オートパイロット”を有するのと同じである。これは、オペレータが、見ている仮想空間において集中することを可能にし、調査している環境の壁の中に進むことについての心配をする必要がない。しかしながら、この第2の技術は、ビューワに、コースを変更するか、飛行経路に沿って見られた興味ある領域を調査する柔軟性を与えない。
【0008】
上述した2つのナビゲーション方法を使用し、これらの個々の欠点を最小にしながら双方の技術の利点を実現することが望ましい。柔軟なナビゲーション技術を、無侵入無痛完全検査を行うために仮想3D空間において表された人間または動物の器官の調査に用いることが望ましい。望まれるナビゲーション技術は、さらに、滑らかな経路と器官中および周囲の完全な調査とを保証しながら柔軟性を許す、オペレータによる3D空間における仮想器官の完全な調査を考慮する。加えて、器官を見るのに必要な計算を最小にする技術を使用することによって、実時間設定における器官の調査を表示することができることが望ましい。また、望まれる技術は、任意の仮想物体の調査に等しく適用可能であるべきである。
【0009】
本発明の他の目的は、見ている物体の一部の視覚化をカスタマイズするために、特定の体積要素を変化する程度に透明または半透明にするために、不透明係数を表現における各々の体積要素に割り当てることである。物体の断面を、前記不透明係数を使用して構成することもできる。
【0010】
発明の要約
本発明は、人間の器官のような物体の3次元視覚化画像を、体積視覚化技術を使用して発生し、前記仮想画像を、オペレータが、予め規定された飛行経路に沿って移動し、ポリープ、嚢胞、または、器官における他の異常な特徴を識別するために、前記予め規定された経路から離れて、前記画像における興味のある特定の部分に位置および見る角度の双方を調節することを可能にするガイド付きナビゲーションシステムを使用して仮想画像を調査する。
【0011】
物体の仮想調査を行う方法は、コントラストエージェントの存在によって膨張した物体の少なくとも1つの画像化走査を行うステップを含む。加えて、コントラストエージェントを取り除いた前記物体の少なくとも1つの画像化走査を得る。前記走査を、複数のボクセルで形成された対応する体積データセットに変換する。次に、画像分割を行い、各々の走査のボクセルを複数のカテゴリに分類する。各々の走査の体積データセットを、共通座標系に登録する。次に、表示動作を、前記体積データセットの少なくとも2つにおける対応する画像が実際的に同時に表示される場合、行うことができる。前記体積データセットの1つにおいて行われる仮想ナビゲーション動作の結果、対応するナビゲーション動作が、少なくとも1つの他の体積データセットにおいて起こる。
【0012】
好適には、前記膨張した物体の少なくとも1つの走査は、前記物体の横走査および冠状走査を含む。同様に、前記コントラストエージェントを取り除いた物体の走査は、前記物体の横走査および冠状走査を含むことが望ましい。この手順は、特に、前記物体が膀胱である場合、仮想膀胱鏡検査を行うのによく適合する。この場合において、前記走査は、一般に、磁気共鳴画像化の形態をとり、前記コントラストエージェントを尿とすることができる。
【0013】
本発明による他の方法は、物体の仮想調査を行うものである。この方法において、前記物体に画像化走査を行い、画像走査データを得る。得られた画像走査データを、複数の体積ユニットまたはボクセルに変換する。前記ボクセル間を補間することによって、拡張されたデータセットを発生する。次に、画像分割を行い、前記ボクセルを複数のカテゴリに分類することができる。前記物体内の体積を、前記物体の管腔内の種ボクセルからの領域成長アルゴリズムを使用すること等によって、前記拡張されたデータセットから抽出する。次に、解像度が低下したデータセットを、前記拡張されたデータセットから発生する。前記拡張されたデータセットからのデータを効率的に格納し、呼び出すために、このデータを、ツリーデータ構造において格納する。次に、画像を、前記拡張されたデータセットおよび解像度が低下したデータセットの双方に関して描画することができる。次に、これらの画像のうち一方を、見るために選択する。一般に、前記解像度が低下したデータセットを、ナビゲーションまたは画像相互作用中に選択し、前記拡張されたデータセットを、高解像度静止表示に関して選択する。
【0014】
仮想血管造影法を行う方法も提供する。この方法において、大動脈の少なくとも一部の画像化走査データを得る。前記画像化走査データを、複数のボクセルを含む体積表現に変換する。前記体積表現を分割し、前記ボクセルを複数のカテゴリのうちの1つに分類する。次に、前記分割された体積表現を解析し、大動脈壁における動脈瘤の少なくとも一部のボクセル表示を識別する。識別された前記動脈瘤の部分から、少なくとも1つの閉鎖表面を、動脈瘤の少なくとも一部のボクセル表示の周囲に発生する。前記閉鎖表面は、前記動脈瘤の輪郭の推定を与える。ナビゲーション経路を、前記大動脈の管腔を通じて確立することができ、長さ、直径、体積および配置のような、前記動脈瘤の特徴を決定することができる。
【0015】
前記仮想血管造影法を行う方法を使用し、動脈瘤の進行を検出および監視することができ、前記仮想血管造影法を行う方法を、ステントグラフトを配置する必要がある特徴を決定するのに使用することができる。
【0016】
複数のボクセルで形成された中空物体の3次元画像表現に関するスケルトンを規定する方法も提供する。ルートボクセルを、前記中空物体内で最初に確認する。次に、距離マップを、前記中空物体内のすべてのボクセルに関して発生する。前記距離マップを、重み付きユークリッド距離を有する26接続キュービックプレートを使用して形成する。前記距離マップにおける局所最大値を有するこれらのボクセルを、前記中空物体におけるブランチの終点として識別する。各々の局所最大値ボクセルに関して、前記ルートボクセルに接続された最も短い経路、または、以前に規定された最短経路を決定する。最短経路の集合を、前記物体の粗いスケルトンとする。この技術は、呼吸器系および心臓血管系のような多ブランチ構造に特に好適である。
【0017】
【発明の実施の形態】
本願において記載した方法およびシステムを、調査すべき任意の物体に適用することができるが、説明する好適実施形態は、人体における器官、特に結腸の調査である。結腸は、長く、捩れており、患者から金銭と、物理的プローブの不快感および危険性との双方を軽減する仮想調査に特に好適である。制限なく、調査できる器官は、肺、胃、胃腸系の一部、心臓および血管を含む。
【0018】
図1は、体積視覚化技術を使用する仮想大腸鏡検査を行うのに必要なステップを示す。ステップ101は、必要ならば、医師または特定の走査機器のいずれかによって調査のために見るために走査すべき結腸を準備する。この準備は、前記結腸を、口から摂取され、胃を通過した後、前記結腸に入る“カクテル”または液体で洗浄することを含む。前記カクテルは、患者に前記結腸内に存在する排泄物を放出させる。使用する物質の一例は、ゴリテリーである。加えて、結腸の場合において、前記結腸を走査および調査が容易になるように拡張するために、空気またはCO2を前記結腸内に押し入れることができる。これを、直腸内に配置された小さいチューブによって行い、約1000ccの空気を前記直腸内に押し入れ、前記直腸を膨張させる。使用するスキャナの形式に応じて、前記結腸における排泄物を前記結腸の壁自身から識別するために、患者にバリウムのような対照物質を飲ませ、なんらかの排出されていない大便を覆う必要があるかもしれない。代わりに、前記結腸を仮想的に調査する方法は、本明細書において後に説明するように、前記仮想調査の前またはその間に、仮想排泄物を除去することができる。ステップ101は、図1において点線によって示すように、すべての調査において行う必要があるわけではない。
【0019】
ステップ103は、調査すべき器官を走査する。前記スキャナを、結腸を走査するスパイラルCTスキャナ、または、例えばキセノンガスで標識付けされた肺を走査するゼニタMRI機械のような、当該技術分野においてよく知られた装置とすることができる。前記スキャナは、前記体積視覚化に必要なデータを発生するために、呼吸を止めている間に、ボディの周りの異なった位置から多数の画像を得ることができなければならない。単一CT画像の一例は、5mm幅、1:1ないし2:1ピッチのX線ビームを使用し、40cm視界を、前記結腸の左結腸曲の上部から直腸まで行う。
【0020】
前記物体の個別のデータ表現を、走査以外の他の方法によって発生することができる。参照によってここに含まれる、1988年7月26日に出願され、1991年8月8日に発行された、カウフマンによる“連続的な3次元幾何学的表現を、3次元ボクセルベースシステム内の個別の3次元ボクセルベース表現に変換する方法”というタイトルの米国特許第5038302号において記載されている技術によって、物体を表現するボクセルデータを幾何学的モデルから得ることができる。加えて、本発明によって3次元ボクセルに変換し、調査することができる画像のコンピュータモデルから、データを発生することができる。この形式のデータの一例は、スペースシャトルを取り巻く乱流のコンピュータシミュレーションである。
【0021】
ステップ204は、前記走査された画像を3次元体積要素(ボクセル)に変換する。結腸を調査する好適実施形態において、前記走査データを、1mmまたは2.5mmの増量において5mm厚スライスに再フォーマットし、1mmスライスにおいて再構成し、各々のスライスを512×512画素として表す。こうすることによって、約1mm3のボクセルを形成する。このようにして、多数の2Dスライスを、前記走査の長さに応じて発生する。次に、前記2Dスライスの組を、3Dボクセルに再構成する。前記スキャナからの2D画像の3Dボクセルへの変換処理を、走査機械それ自身か、コンピュータのような別個の機械のいずれかによって、当該技術分野において既知の技術で行うことができる(例えば、参照によってここに含まれる、1988年11月11日に出願され、1991年1月15日に発行された、カウフマン他による“ボクセルベースデータを格納し、アクセスし、処理する方法および装置”というタイトルの米国特許第4985856号を参照)。
【0022】
ステップ105は、オペレータが調査すべき選択された器官の部分を規定することを可能にする。医師は、ポリープが発達しそうな前記結腸の特定の部分に興味を持つかもしれない。見るべき経路の開始点および終了点を、医師/オペレータによって指示することができる。慣例的なコンピュータおよびコンピュータインタフェース(例えば、キーボード、マウスまたはスペースボール)を使用し、調査すべき前記結腸の部分を指定することができる。座標を有するグリッド系をキーボード入力に使用することができ、または、医師/オペレータは、所望の点において“クリック”することができる。前記結腸の全体像も、望むならば見ることができる。
【0023】
ステップ107は、調査している仮想器官の計画されたまたはガイド付きナビゲーション動作を行う。ガイド付きナビゲーション動作を行うことを、任意の時間においてオペレータによって手動で調節することができる、予め規定されたまたは自動的に予め決定された飛行経路に沿って環境を通るナビゲーティングとして規定する。前記走査データを3Dボクセルに変換した後、前記器官の内部を、選択された開始点から選択された終了点まで通過しなければならない。前記仮想調査を、前記仮想空間を通って進む、前記終了点をさすレンズを有する小さいカメラを有することにおいてモデル化する。前記ガイド付きナビゲーション技術は、前記カメラとのあるレベルの相互作用を与え、その結果、前記カメラは、オペレータの相互作用がない場合において仮想空間を通って自動的にナビゲートすることができ、同時に、必要な場合、オペレータが前記カメラを操作することを可能にする。ガイド付きナビゲーションを達成する好適な例は、ポテンシャル場を用いて前記カメラの動きを制御する、図2および3において詳細に示す物理的に基づいたカメラモデルを使用することである。
【0024】
ステップ107と同時に行うことができるステップ109は、前記ガイド付きナビゲーション動作の選択された経路の沿った前記カメラの視点から前記器官の内部を表示する。3次元表示を、マーチングキューブ技術のような当該技術分野において既知の技術を使用して発生することができる。しかしながら、前記結腸の実時間表示を発生するために、前記仮想器官の表示に必要な膨大な数のデータの計算を減らす技術が必要である。図9は、この表示ステップをより詳細に示す。
【0025】
図1に示す方法を、ボディにおける多数の器官を同時に走査することに用いることもできる。例えば、患者を、結腸および肺の双方における癌の成長に関して調査してもよい。図1の方法を、ステップ103において興味あるすべての領域を走査し、ステップ105において調査すべき現在の器官を選択するように変更する。例えば、医師/オペレータは、最初に結腸を選択して仮想調査し、後に肺を調査してもよい。代わりに、異なった専門の二人の異なった医師が、彼らの個々の専門に関係する、異なった走査された器官を仮想調査してもよい。ステップ109に続いて、調査すべき次の器官を選択し、その部分を規定し、調査する。これを、調査が必要なすべての器官が処理されるまで続ける。
【0026】
図1に関連して説明したステップを、体積要素によって表すことができる任意の物体の調査に用いることもできる。例えば、建築構造物または無生物を同様に表し、調査することができる。
【0027】
図2は、ステップ107におけるガイド付きナビゲーション技術を行う“サブマリン”カメラ制御モデルを示す。ガイド付きナビゲーション中にオペレータ制御がない場合、通常設定ナビゲーションは、前記結腸の一方の選択された端から他方までの飛行経路に沿って前記カメラを自動的に向ける計画されたナビゲーションと同様である。前記計画されたナビゲーション段階中、前記カメラは、前記結腸の表面のより良好な視界を得る、前記結腸の中心においてとどまる。興味ある領域に出合った場合、前記ガイド付きナビゲーションを使用する仮想カメラのオペレータは、前記カメラを相互作用的に特定の領域に近づけ、前記カメラの動きおよび角度を管理し、前記直腸の壁と不本意に衝突することなしに前記興味ある領域を詳細に研究することができる。オペレータは、前記カメラを、キーボードやマウスのような標準的なインタフェース装置、または、スペースボールのような非標準的な装置で制御することができる。カメラを仮想環境において完全に操作するために、前記カメラに関して6自由度が必要である。前記カメラは、水平、垂直およびZ方向(軸217)において移動することができ、他の3自由度(軸219)において回転することができ、前記カメラが仮想環境のすべての側面および角度に移動し、操作することができるようにしなければならない。前記ガイド付きナビゲーションに関するカメラモデルは、2つの質点x1203およびx2205を接続する、延びない重さのないロッド201を含み、双方の質点はポテンシャル場215を受ける。前記ポテンシャル場を、前記カメラを前記器官の壁から押しのけるために、前記器官の壁における高さとして規定する。
【0028】
前記質点の位置をx1およびx2によって与え、これらは同じ質量mを有するとする。カメラをサブマリンx1203の頭に取り付け、このカメラの見る方向はx2x1と一致する。前記サブマリンは、前記2個の質点が前に規定したポテンシャル場V(x)からの力と、任意の摩擦力と、任意のシミュレートされた外力とによって影響を受けると、並進運動と、前記モデルの質量xの中心の周囲の回転とを行うことができる。x1、x2と、xとの間の関係は、以下の通りである。
【数1】
ここで、r、θおよびφは、ベクトルxx1の極座標である。前記モデルの運動エネルギーTを、x1およびx2の運動の運動エネルギーの和として規定する。
【数2】
【0029】
次に、前記サブマリンモデルの動きに関する式を、ラグランジュ方程式を使用して得る。
【数3】
ここで、qjsは、前記モデルの一般化座標であり、これを、時間tの変数として考えることができる。
【数4】
ψは我々のカメラシステムのロール角を示し、これを後に説明する。Fisを一般化された力と呼ぶ。前記サブマリンの制御を、シミュレートされた外力をx1に用いることによって行う。
【数5】
x1およびx2の双方は、前記ポテンシャル場からの力と、各々の質点の速度の反対方向において作用する摩擦とによって影響を受けるとする。したがって、発生される力は、以下のように式化される。
【数6】
ここで、kは前記システムの摩擦係数を示す。外力Fextを、図2に示すように、発生された画像における所望の方向207におけるマウスボタンの単純なクリックによって、オペレータによって適用する。このとき、このカメラモデルはその方向において動く。これは、オペレータがマウスボタンの1クリックのみによって前記カメラの少なくとも5自由度を制御することを可能にする。式(2)、(3)および(5)から、我々のサブマリンモデルの5つのパラメータの加速度を以下のように得ることができる。
【数7】
ここで、
【外1】
および
【外2】
は、各々、xの一次導関数および二次導関数を示し、
【数8】
は、点xにおけるポテンシャルの勾配を示す。
【外3】
の
【外4】
および
【外5】
の
【外6】
を、遠心力およびコリオリ力と呼び、これらは、前記サブマリンの角速度の交換に関係がある。前記モデルが、前記サブマリンのロッドに関して規定された慣性モーメントを持たないため、これらの項は、φの数値計算のオーバーフローを生じる傾向にある。幸いにも、これらの項は、前記サブマリンモデルの角速度が重要な場合にのみ重要になり、これは、本質的に、前記カメラが速く動き過ぎることを意味する。器官を適切に見ることができないため、カメラを非常に速く動かせるようにすることは無意味であるため、前記オーバーフロー問題を回避するため、我々のインプリメンテーションにおいてこれらの項を最小にする。
【0030】
式(6)の最初の3つの式から、以下の条件を満たす場合、前記サブマリンを、前記ポテンシャル場に対する外力によって進ませることはできないことが分かる。
【数9】
前記サブマリンの速度と、外力
【外7】
とは、我々のインプリメンテーションにおいて上限を有するため、前記物体の境界において十分高いポテンシャル値を割り当てることによって、前記サブマリンが、前記環境において前記物体または壁に決して衝突しないことを保証することができる。
【0031】
上述したように、前記カメラシステムのロール角Ψを考える必要がある。可能な1つのオプションは、オペレータに角度の完全な制御を与える。しかしながら、オペレータは、前記カメラを前記モデルのロッドの周りで自由に回転させることができるが、彼または彼女は、容易に方角を見失うかもしれない。好適な技術は、前記カメラの上部の方向を、質量m2を有する振り子301に接続し、この振り子は、図3に示すように前記サブマリンのロッドの周りを自由に回転するとする。振り子の方向r2を、
r2=γ2(cosθcosφsinψ
+sinφcosψ,cosθsinφsinψ
−cosφcosψ,−sinθsinψ)
のように表し、この式は、前記サブマリンの運動に伴うこの振り子の正確な運動を計算することができるが、システム方程式を複雑にし過ぎる。したがって、ロール角Ψを除くすべての一般化座標が一定であるとし、したがって、前記振り子システムに関する独立の運動エネルギーを、
【数10】
として規定する。これは、前記モデルを前記ロール角に関して簡単にする。このモデルにおいて、重力
【数11】
が質点m2において作用し、Ψの加速度を、
【数12】
のようなラグランジュ方程式を使用して得ることができる。式(6)および(7)から、一般化座標
【外8】
およびこれらの導関数
【外9】
を、
【数13】
のようなテイラー級数を使用し、前記サブマリンを自由に動かすことによって、漸近的に計算する。前記サブマリンの動きを滑らかにするために、時間ステップhを、前記動きを滑らかにするためにできるだけ小さい、しかし、計算コストを減少するのに必要な程度に大きい平衡値として選択する。
【0032】
ポテンシャル場の定義
図2におけるサブマリンモデルにおけるポテンシャル場は、前記仮想器官における境界(壁または他の物)を、前記サブマリンカメラが前記壁または他の境界と衝突しないことを保証するために、高いポテンシャルを前記境界に割り当てることによって規定する。前記カメラモデルがオペレータによって高いポテンシャル領域に移動するように試みられた場合、前記カメラモデルは、例えばオペレータが前記境界の背後またはポリープの内側の器官を調査することを望まない限り、そうすることが禁止される。仮想大腸鏡検査を行う場合において、ポテンシャル場値を、体積測定結腸データの各々の部分(体積要素)に割り当てる。興味ある特定の領域が図1のステップ105において開始点および終了点によって指定された場合、前記走査された結腸の選択された領域内のボクセルを、慣例的なブロッキング操作によって識別する。次に、ポテンシャル値を、前記選択された体積のすべてのボクセルxに、以下の3つの距離値、すなわち、終了点からの距離dt(x)と、前記結腸表面からの距離ds(x)と、前記結腸空間の中心線からの距離dc(x)とに基づいて割り当てる。dt(x)を、慣例的な成長戦略によって計算する。前記結腸表面からの距離ds(x)を、表面ボクセル内部からの慣例的な成長技術を使用して計算する。dc(x)を決定するために、最初に前記ボクセルから前記結腸の中心線を抽出し、次に、dc(x)を、慣例的な成長戦略を使用し、前記結腸の中心線から計算する。
【0033】
前記ユーザに指定された開始点とユーザに指定された終了点とによって規定された前記選択された結腸領域の中心線を計算するために、ds(x)の最大値を捜し出し、dmaxと示す。次に、興味ある領域内の各々のボクセルに関して、dmax−ds(x)のコスト値を割り当てる。このようにして、前記結腸表面に近いボクセルは、高いコスト値を有し、前記中心線に近いボクセルは、比較的低いコスト値を有する。次に、前記コスト割り当てに基づいて、当該技術分野において既知である単一始点最短路技術を用い、始点から終点までの最小コスト経路を効率的に計算する。この低コストラインは、調査することを望む前記結腸区分の中心線または輪郭を示す。前記中心線を決定するこの技術は、本発明の好適な技術である。
【0034】
興味ある領域の内部のボクセルxに関するポテンシャル値V(x)を計算するために、以下の式を用いる。
【数14】
ここで、C1、C2、μおよびυは、タスクに関して選択された定数である。前記仮想カメラと前記仮想結腸表面との間のどのような衝突も回避するために、十分大きいポテンシャル値を、前記結腸の外部のすべての点に割り当てる。前記ポテンシャル場の勾配は、したがって、前記サブマリンモデルカメラが、移動しているとき、前記結腸壁に決して衝突しないようにするため、重要になる。
【0035】
前記結腸における経路の中心線を決定する他の技術は、“ピールレイヤ”技術と呼ばれ、この技術を図4ないし8において示す。
【0036】
図4は、前記体積測定結腸の2D断面を、示している前記結腸の2つの側壁401および403と共に示す。2つのブロッキング壁を、オペレータによって、調査する興味のある前記結腸の選択を規定するために、選択する。前記ブロッキング壁を超えて見ることはできない。これは、前記仮想表現を表示する場合の計算数を減少するのを助ける。前記ブロッキング壁は、側壁と共に、調査すべき前記結腸の含まれる体積測定形状を明らかにする。
【0037】
図5は、前記仮想調査の飛行経路の2つの終点、開始体積要素501および終了体積要素503を示す。前記開始点および終了点を、図1のステップ105においてオペレータによって選択する。前記開始点及び終了点と結腸側面との間のボクセルを、図6において“x”で示す領域によって示されるように識別し、マークする。前記ボクセルは、画素の3次元表現である。
【0038】
次に、前記ピールレイヤ技術を、図6における識別されマークされたボクセルに用いる。すべてのボクセルの最外層(結腸壁に最も近い)を、残りのボクセルのただ1つの内部層が存在するようになるまで、1つずつ剥ぎ取る。すなわち、中心点から最も遠く離れた各ボクセルを、前記開始ボクセルと終了ボクセルとの間の経路の分離を生じなくなるまで除去する。図7は、前記仮想結腸におけるボクセルの除去の多数の反復が完了した後の中間結果を示す。前記結腸壁に最も近いボクセルを除去した。図8は、すべての剥ぎ取り反復が完了した後の、前記結腸の中心を下った前記カメラモデルの最終飛行経路を示す。これは、本質的に、前記結腸の中心におけるスケルトンを発生し、前記カメラモデルに関する所望の飛行経路になる。
【0039】
Zバッファ可視性支援
図9は、器官の仮想3次元体積表現においてカメラモデルによって見られる仮想画像を表示する実時間可視性技術を説明する。図9は、図1におけるステップ109に対応する変更されたZバッファを使用する表示技術を示す。前記カメラモデルからあるいは見ることができるボクセルの数は、非常に多い。計算し、視覚化しなければならない要素(またはポリゴン)の合計数を、前記走査された環境におけるボクセルの全体の組から減らさなければ、計算の全体数は、大きな内部領域に関して、前記視覚化表示プロセスを非常に遅くする。しかしながら、本発明において、前記結腸表面において可視であるこれらの画像のみを、表示に関して計算する必要がある。前記走査された環境を、より小さい区分またはセルに再分することができる。次に、前記Zバッファ技術は、前記カメラから可視であるセルの部分のみをレンダする。前記Zバッファ技術を、3次元ボクセル表示にも使用する。変更されたZバッファの使用は、計算すべき可視ボクセルの数を減らし、医師または医療技術者による前記仮想結腸の実時間調査を可能にする。
【0040】
前記中心線をステップ107において計算した興味ある領域を、前記表示技術を用いる前に再分する。セルは、可視性単位になるボクセルの集合的グループである。各セルは、他のセルを見ることができる多数の門を含む。前記結腸を、前記選択された開始点において開始し、中心線1001に沿って前記終了点に向かって移動することによって再分する。次に、前記結腸を、前記中心経路に沿った予め規定されたしきい値距離に達した場合、セル(例えば、図10におけるセル1003、1005および1007)に分割する。前記しきい値距離は、前記視覚化技術を行うプラットフォームの特性と、その格納及び処理容量とに基づく。前記セルサイズは、前記プラットフォームによって格納し、処理することができるボクセル数に直接関係する。大きく変化させることができるが、しきい値距離の一例は5cmである。各々のセルは、図10に示すようなセルの外側を見る門として2つの断面を有する。
【0041】
図9におけるステップ901は、前記カメラを現在含む選択された器官内のセルを識別する。この現在セルを、前記カメラの方向を与える可視であるすべての他のセルと同様に表示する。ステップ903は、以下にさらに詳細に説明するように、(規定された門を経て)前記カメラから可能的に可視であるセルの階層データのスタブツリー(ツリー図)を構成する。前記スタブツリーは、前記カメラに可視であるかもしれないすべてのセルに関するノードを含む。1つ以上のセルは単一方向において可視であるため、前記セルのいくつかは、妨害する物体がなにもなければ、透明であるかも知れない。ステップ905は、隣接するセルエッジの交点を含むセルからのボクセルのサブセットを格納し、どのセルが可視であるかをより効率的に決定するために、これらを前記スタブツリーの外部エッジにおいて格納する。
【0042】
ステップ907は、ループノードが前記スタブツリーにおいて存在するかどうかをチェックする。ループノードは、単一セルの2つ以上のエッジが同じ近くのセルに接する場合、生じる。これは、単一セルが他のセルによって囲まれている場合、生じるかもしれない。ループノードが前記スタブツリーにおいて識別された場合、前記方法はステップ909に続く。ループノードがない場合、前記プロセスはステップ911に進む。
【0043】
ステップ909は、前記ループノードを形成する2つのセルを1つの大きいノードに合体させる。次に、前記スタブツリーを適宜修正する。これは、ループノードのため、同じセルを2回見るという問題を除去する。このステップを、すべての識別されたループノードにおいて行う。次に、前記プロセスはステップ911に続く。
【0044】
次に、ステップ911は、前記Zバッファを最も大きいZ値で開始する。このZ値は、前記スケルトン経路に沿って前記カメラからの距離を規定する。次に、前記ツリーを横切り、各ノードにおける交点値を最初にチェックする。ノード交点が覆われている場合、現在の門シーケンスが閉じることを意味し(前記Zバッファテストによって決定される)、前記ツリーの現在のブランチの横断を停止する。ステップ913は、前記ブランチの各々を横断し、前記ノードが覆われていないかどうかをチェックし、覆われていない場合、これらを表示する。
【0045】
ステップ915は、オペレータのスクリーンにおいて表示すべき画像を、ステップ913において識別された可視セル内の体積要素から、合成による体積レンダリングのような当該技術分野において既知の種々の技術の1つを使用して構成する。示されるセルは、可能的に可視であると識別されたものだけである。この技術は、実時間表示を達成するために計算を必要とするセルの数を制限し、相応して、よりよい性能のための表示の速度を増す。この技術は、実際に表示されようとされなかろうとすべての可能な可視データ点を計算する従来技術に対する改良である。
【0046】
図11Aは、ガイド付きナビゲーションによって調査されており、オペレータに表示する必要がある器官の2次元画像表現である。器官1101は、2つの側壁1102と、経路の中心における物体1105とを示す。この器官を、4つのセルA1151、B1153、C1155およびD1157とに分割した。カメラ1103は、セルD1157の方に面しており、コーン状視野を認めることができる視野ベクトル1107、1108によって規定される視野を有する。可能的に見ることができるセルは、セルB1153、C1155およびD1157である。セルC1155は、セルBによって完全に取り囲まれ、したがって、ノードループを構成する。
【0047】
図11Bは、図11Aにおけるセルから構成されたスタブツリーの表現である。前記カメラを含むノードA1109は、前記ツリーのルートにある。妨害されない可視の経路である視線または視野コーンは、ノードB1110へ引かれる。ノードBは、ノードC1112およびノードD1114の双方への直接可視視線を有し、これらを接続する矢印によって示す。前記見るカメラの方向におけるノードC1112の視線は、ノードB1110と結合する。ノードC1112およびノードB1110は、したがって、図11Cに示すような1つの大きいノードB´1122に合体する。
【0048】
図11Cは、ノードB´1122(ノードBおよびノードCの双方を含む)およびノードD1114に隣接する前記カメラを含むノードA1109を示す。ノードA、B´およびDを、少なくとも部分的にオペレータに対して表示する。
【0049】
図12A−12Eは、視線を妨害する物体を含むセルと共に、前記変更されたZバッファの使用を説明する。物体を、前記仮想結腸の一部におけるある老廃物とすることができる。図12Aは、10のポテンシャルセル、A1251、B1253、C1255、D1257、E1259、F1261、G1263、H1265、I1267およびJ1269を有する仮想空間を示す。前記セルのいくつかは、物体を含む。カメラ1201がセルI1267に位置し、視野ベクトル1203によって示されるようにセルF1261の方を向いている場合、スタブツリーを、図9における流れ図によって示す技術によって発生する。図12Bは、図12Aにおいて示すような仮想表現に関して示す交点ノードによって発生されたスタブツリーを示す。図12Bは、カメラ1210を含むため、前記ツリーのルートノードとしてセルI1267を示す。セルFが前記カメラの視線に直接接続するため、(矢印で示すように)ノードI1211はノードF1213をさす。ノードF1213は、ノードB1215およびノードE1219の双方をさす。ノードB1215は、ノードA1217をさす。ノードC1202は、カメラ1201による視線から完全に妨げられ、したがって、前記スタブツリーに含まれない。
【0050】
図12Cは、ノードI1211をオペレータ用のディスプレイにおいてレンダした後のスタブツリーを示す。ノードI1211を、すでに表示されたため、前記スタブツリーから除去し、ノードF1213はルートになる。図12Dは、ノードF1213をここでレンダし、ノードI1211を連結することを示す。次に、矢印によって接続された前記ツリーにおける次のノードをチェックし、これらがすでにカバーされたか(すでに処理されたか)を確かめる。この例において、ノードB515(と、したがって従属するノードAと)を前記ディスプレイにおいてレンダする必要がないため、セルI1267において位置するカメラからの交差したノードのすべてがカバーされている。
【0051】
図12Eは、ノードE515をチェックし、その交点がカバーされたかどうかを決定することを示す。カバーされているため、図12A−12Eのこの例においてレンダされたノードは、ノードIおよびFだけであり、ノードA、BおよびEは不可視であり、これらのセルが表示される準備ができている必要はない。
【0052】
図9に示す変更されたZバッファ技術は、より少ない計算を可能にし、ボクセル、またはポリゴンのような他のデータ要素によって表された物体に用いることができる。
【0053】
図13は、その壁の一方に沿って大きいポリープが存在する結腸の2次元仮想ビューを示す。図13は、さらに調査すべき患者の結腸の選択された区分を示す。このビューは、2つの結腸壁1301および1303を、1305として示す腫瘍と共に示す。層1307、1309および1311は、前記腫瘍の内部層を示す。医師が前記ポリープまたは腫瘍の層を剥がし、任意の癌または他の有害な要素に関して固まりの内部を見ることができることが望ましい。このプロセスは、実際には、前記固まりへ実際に切断することなしに、前記固まりの仮想生体組織検査を行う。前記結腸をボクセルによって仮想的に表現したら、物体の層を剥ぎ取るプロセスは、図4ないし8と共に説明したように、容易に行われる。特定の断面を調査できるように、前記固まりをスライスすることもできる。加えて、ユーザが規定したスライス1319を、前記腫瘍において任意の方法で形成することができる。ボクセル1319を、以下に説明するように、剥ぎ取ることができ、変更することもできる。
【0054】
伝達関数を興味ある領域における各ボクセルに行うことができ、各ボクセルの半透明性を表す係数を変化させることによって、前記物体を透明、半透明または不透明にすることができる。不透明度係数を、各ボクセルにその密度に基づいて割り当てる。次に、マッピング関数は、前記密度値を、その半透明性を表す係数に変換する。高密度走査ボクセルは、壁か、単にオープンスペースの以外の他の密な物質のいずれかを示す。次に、オペレータまたはプログラムルーチンは、ボクセルまたはボクセルのグループの不透明度係数を変化させ、これらが前記サブマリンカメラモデルに対して透明または半透明にすることができる。例えば、オペレータは、腫瘍全体の内部または外部の腫瘍を見ることができる。また、透明ボクセルを形成し、図9の表示ステップに関して存在しないかのようにみせる。前記オブジェクトの区分の合成を、当該区分におけるボクセルの不透明度係数の重み付き平均を使用して形成することができる。
【0055】
医師が、ポリープの種々の層を見て癌の領域を探すことを望む場合、これを、ポリープの外部の層1305を除去し、第1層1307を生じることによって行うことができる。加えて、第1内部層1307を取り除き、第2内部層1309を見ることができる。前記第2内部層を取り除き、第3内部層1311等を見ることができる。医師は、ポリープ1305をスライスし、所望の区分内のこれらのボクセルのみを見ることもできる。前記スライスする領域を、完全にユーザが規定することができる。
【0056】
不透明係数を加えることは、仮想システムの調査を援助する他の方法でも使用されるかもしれない。老廃物が存在し、ある既知の範囲内の他の特性として密度を有する場合、前記老廃物を、調査中その不透明度係数を変更することによって前記仮想カメラに対して透明にすることができる。これは、患者が前記手順の前に腸洗浄剤を摂取することを回避することを可能にし、前記調査をより速くより容易にする。他の物体を、実際の適用に応じて同様に消すことができる。加えて、ポリープのようないくつかの物体を、コントラスト剤とそれに続く適切な伝達関数の使用によって電子的に強調することができる。
【0057】
図14は、本明細書において記載の技術を使用して人間の器官のような物体の仮想調査を行うシステムを示す。患者1401は、プラットフォーム上で横になり、走査装置1405は、調査すべき1つ以上の器官を含む領域を調査する。走査装置1405は、前記患者の画像を実際に取得する走査部分1403と、エレクトロニクス部分1406とを含む。エレクトロニクス部分1406は、インタフェース1407と、中央処理ユニット1409と、走査データを一時的に格納するメモリ1411と、データを仮想ナビゲーションプラットフォームに送る第2インタフェース1413とを具える。インタフェース1407および1413を、単一のインタフェース構成要素に含めることができ、または、同じ構成要素とすることができる。部分1406における構成要素を、慣例的なコネクタによって共に接続することができる。
【0058】
システム1400において、装置1403の走査部分から与えられたデータを、処理するために部分1405に伝送し、メモリ1411に格納する。中央処理ユニット1409は、走査された2Dデータを3Dボクセルデータに変換し、その結果をメモリ1411の他の部分に格納する。代わりに、前記変換されたデータを、仮想ナビゲーション端末1416に伝送するために直接インタフェースユニット1413に送ることもできる。前記2Dデータの変換を、インタフェース1413から伝送した後に、仮想ナビゲーション端末1416において行ってもよい。好適実施形態において、オペレータが前記仮想調査を行うために、前記変換されたデータを、キャリア1414によって仮想ナビゲーション端末1416に伝送する。前記データを、記憶媒体において格納し、端末1416に物理的に伝送するか、衛星伝送を使用することによるような他の慣例的な方法において伝送することもできる。
【0059】
前記走査されたデータをその3D表現に、前記視覚化レンダリングエンジンが前記データを3D形式にすることを要求するまで、変換しなくてもいい。これは、計算ステップとメモリ記憶空間とを節約する。
【0060】
仮想ナビゲーション端末1416は、前記仮想器官または他の走査された画像を見るスクリーンと、エレクトロニクス部分1415と、キーボード、マウスまたはスペースボールのようなインタフェース制御1419とを含む。エレクトロニクス部分1415は、インタフェースポート1421と、中央処理ユニット1423と、前記端末を動作するのに必要な他の構成要素1427と、メモリ1425とを具える。端末1416における構成要素を、慣例的なコネクタによって一緒に接続する。前記変換されたボクセルデータを、インタフェースポート1421によって受け、メモリ1425に格納する。次に、中央処理ユニット1423は、前記3Dボクセルを仮想表現に組み立て、図2および3において説明したようなサブマリンカメラモデルを動作し、前記仮想調査を行う。前記サブマリンカメラが前記仮想器官を通って移動しているとき、図9において説明したような可視化技術を使用し、前記仮想カメラから可視のこれらの領域のみを計算し、これらをスクリーン1417において表示する。グラフィックスアクセラレータを、前記表示の発生において使用することもできる。オペレータは、インタフェース装置1419を使用し、前記走査された肉体のどの部分を調査したいかを示すことができる。さらに、インタフェース装置1419を使用し、図2およびその伴う記述において考察したように望まれるサブマリンカメラを制御し、動かすことができる。端末部分1415を、ニューヨーク州立大学ストーニーブロック校におけるコンピュータ科学部から一般的に利用可能なCube−4専用システムボックスとすることができる。
【0061】
走査装置1405および端末1416、またはこれらの部品を、同じユニットの一部とすることができる。単一のプラットフォームを使用し、前記走査画像データを受け、必要ならこれを3Dボクセルに接続し、前記ガイド付きナビゲーションを行う。
【0062】
システム1400における重要な特徴は、前記仮想器官を、患者の存在なしに、より後の時間において調査することができることである。加えて、前記仮想調査を、患者を走査している間に行うことができる。前記走査データを、2人以上の医師が前記器官の内部を同時に見ることを可能にする多数の端末に送ることもできる。したがって、ニューヨークにおける医師が、カリフォルニアにおける医師と同時に患者の器官の同じ部分を、このケースを討論しながら見ることができる。代わりに、前記データを、異なった時間において見ることができる。2人以上の医師が、異なったケースにおける同じデータの彼ら自身の調査を行うことができる。多数の仮想ナビゲーション端末を使用し、同じ走査データを見ることができる。前記器官をデータの個別の組を有する仮想器官として再生することによって、精度、コストおよび可能なデータ操作のような領域において多数の利点がある。
【0063】
さらに、上述した技術を、画像分割動作が続く変更された腸準備動作を用いる改善された電子結腸洗浄技術の使用によって、仮想大腸鏡検査用途において拡張することができ、コンピュータ断層撮影(CT)または磁気共鳴画像(MRI)走査中に前記結腸において残る流体および大便を検出し、前記仮想大腸鏡画像から除去することができる。これらのような技術の使用により、前記結腸の慣例的な物理的洗浄と、これに関係する不便さおよび不快さとは、最小になる、または完全に回避される。
【0064】
図15を参照し、電子結腸洗浄の第1ステップは、腸準備(ステップ1510)であり、このステップは、CTまたは磁気共鳴画像(MRI)走査の前に行い、前記結腸において残る残りの大便および流体が、ガスで満たされた結腸内部および結腸壁の画像特性と重大に異なる画像特性を与える条件を形成することを目的とする。好例の腸準備動作は、前記CTまたはMRI走査の前日中、ニューヨーク、ウェストバリーのE−Z−M,Inc.によって製造されるような2.1%W/Vの硫酸バリウム懸濁液の250cc3回の服用することを含む。前記3回の服用は、1日の間で分散させるべきであり、各々、3回の食事と共に摂ることができる。前記硫酸バリウムは、前記結腸内に残るどのような大便の画像も強調するように働く。前記硫酸バリウムの摂取に加えて、前記CTまたはMRI走査前日中、液体の摂取を好適には増加させる。水を摂取させることもできるが、クランベリージュースは、腸の流体を増加させることが知られており、好適である。前記CT走査の前夜と、前記CT走査の朝との双方において、セントルイス、ミッソーリのマリンクロットによって製造されるMD−ガストロビューとして商業的に利用可能な60mlのアミドトリゾ酸メグルミンおよびアミドトリゾ酸ナトリウム溶液を摂り、前記結腸流体の画像特性を強調することができる。リン酸ナトリウムを前記溶液に加え、前記結腸内の大便を液化することができ、これは、前記結腸流体および残りの大便のより均一な強調を与える。
【0065】
上述した好例の予備腸準備動作は、CTスキャン前の1ガロンのゴリテリーの摂取を必要とするかもしれない慣例的な結腸洗浄プロトコルの必要を除去することができる。
【0066】
前記CTスキャンを行う直前、インディアナポリス、インディアナのエリーリリーアンドカンパニーによって製造される1mlのグルカゴンの静脈注射を施し、結腸の衰弱を最小にすることができる。次に、前記結腸を、CO2または空気のような約1000ccの圧縮ガスを使用して膨張させ、前記ガスを、直腸チューブを経て入れることができる。この時点において、慣例的なCTスキャンを行い、前記結腸の領域からデータを取得する(ステップ1520)。例えば、データを、5mm、1.5−2.0:1ピッチのヘリカルモードにおいて動作するGE/CTIスパイラルモードスキャナを使用して取得することができ、1mmスライスにおいて再生することができ、前記ピッチを、患者の高さに基づいて既知の方法において調節する。120kVpおよび200−280maのルーチンの画像プロトコルを、この動作に使用することができる。前記データを、患者のサイズに応じて34から40cmまで変化する視野において512×512画素のアレイサイズを有する1mm厚スライス画像として取得し、再生することができる。これらのようなスライスの数は、一般に、前記患者の高さに応じて、300から450までこれらの条件の下で変化する。前記画像データ組を、体積要素またはボクセルに変換する(ステップ1530)。
【0067】
画像分割を、多数の方法において行うことができる。画像分割のある現在の方法において、局所近接技術を使用し、前記画像データのボクセルを、同様の強度値にしたがって分類する。この方法において、取得した画像の各ボクセルを、隣接するボクセルのグループに関して評価する。興味あるボクセルは、中央ボクセルと呼ばれ、関係する強度値を有する。各ボクセルに関する分類指標を、前記中央ボクセルの値をその隣接ボクセルの各々と比較することによって確立する。前記隣接ボクセルが前記中央ボクセルと同じ値を有する場合、前記分類指標の値を増加する。しかしながら、前記隣接ボクセルが前記中央ボクセルと異なった値を有する場合、前記中央ボクセルに関する分類指標を減少する。次に、前記中央ボクセルを、最大の指標値を有するカテゴリに分類し、これは、局所隣接ボクセルの中で最も一様な隣接ボクセルを有する。各々の分類は、特定の強度範囲を示し、この強度範囲は、画像化されている1つ以上の材料形式の表現である。さらに、この方法を、混合確率関数を得られた前記類似性分類に用いることによって拡張することができる。
【0068】
画像分割の代わりのプロセスを、2つの主な動作、すなわち、低レベル処理と高レベル特徴抽出として行う。低レベル処理中、身体の輪郭の外部の領域を更なる処理から除外し、前記身体輪郭の内部のボクセルを、良好に規定されたク強度特徴のラスに従って大まかに分類する。例えば、腹部領域のCTスキャンは、良好に規定された強度分布を示す傾向にあるデータ組を発生する。図16のグラフは、このような強度分布を、強度しきい値にしたがって分類することができる4つの良好に規定されたピーク、1602、1604、1606、1608を有する好例のヒストグラムとして示す。
【0069】
腹部CTデータ組のボクセルを、強度しきい値にしたがって4つのクラスタに大まかに分類する(ステップ1540)。例えば、クラスタ1は、これらの強度が140未満のボクセルを含むことができる。このクラスタは、一般に、前記ガスを満たされた結腸内部の最も密度が低い領域に対応する。クラスタ2は、2200を超える強度値を有するボクセルを含むことができる。これらの強度値は、前記結腸内部の強調された大便および流体と、骨とに対応する。クラスタ3は、約900から約1080の範囲における強度を有するボクセルを含むことができる。この強度範囲は、一般に、おそらく結腸には関係しない、脂肪や筋肉のような柔らかい組織をあらわす。次に、残りのボクセルを共にクラスタ4として分類することができ、これは、前記結腸壁(前記結腸壁の周囲の粘膜および一部の体積混合物を含む)と、肺組織および軟骨に関係すると考えられる。
【0070】
クラスタ1および3は、前記結腸壁の識別において特に価値がなく、したがって、仮想大腸鏡検査に関する画像分割手順中、実際的な処理を受けない。クラスタ2に関係するボクセルは、前記結腸壁から大便および流体を分離するのに重要であり、前記高レベル特徴抽出動作中、さらに処理される。低レベル処理を、結腸組織に対応する最も高い可能性を有する第4のクラスタにおいて集中する(ステップ1550)。
【0071】
前記第4クラスタにおける各ボクセルに関して、強度ベクトルを、それ自身およびその周囲のボクセルを使用して発生する。前記強度ベクトルは、所定のボクセルに最も近い隣接ボクセルにおける強度における変化の表示を与える。前記強度ベクトルを確立するのに使用する隣接ボクセルの数は重要ではなく、処理オーバヘッドと精度との間のトレードオフを含む。例えば、簡単なボクセル強度ベクトルを、興味あるボクセルと、3つの相互に直交する軸においてこのボクセルを取り囲む、その前後の隣接ボクセルと、その左右の隣接ボクセルと、その上下の隣接ボクセルとを含む、7ベクトルで確立することができる。図17は、25ボクセル強度ベクトルモデルの形態における好例の強度ベクトルを説明する斜視図であり、この強度ベクトルは、選択されたボクセル1702と、その第1順、第2順および第3順の隣接ボクセルとを含む。選択されたボクセル1702をこのモデルの中心点とし、固定ボクセルと呼ぶ。前記固定ボクセルと同じ平面における12個の隣接ボクセルを含むボクセルの平面スライスを、固定スライス1704と呼ぶ。前記固定スライスに対する隣接平面において、各々5ボクセルを有する2つの最も近いスライス1706がある。第1の最も近いスライス1706の隣は、各々1つのボクセルを有する2つの第2の最も近いスライス1708である。前記第4クラスタにおける各ボクセルに関する強度ベクトルの集合を、局所ベクトル列と呼ぶ。
【0072】
腹部画像に関するデータ組は、一般に、300を超えるスライス画像を含み、スライス画像の各々が512×512ボクセルアレイを有し、各ボクセルが関係する25のボクセル局所ベクトルを有するため、前記局所ベクトル列において特徴解析を行い(ステップ1570)、計算上の負担を減らすことが望ましい。あるこのような特徴解析は、主成分解析(PCA)であり、これを前記局所ベクトル列に適用し、特徴ベクトル列の大きさと、クラスタ4のボクセルに関する直交変換行列とを決定することができる。
【0073】
前記CT画像強度のヒストグラム(図16)は、同等の準備および走査パラメータを与えれば、特定のスキャナに関して、患者間でかなり一定になる傾向にある。この観察を信頼し、直交変換行列を確立することができ、この行列は、同じスキャナを同様の条件の下で使用して獲得された訓練データのいくつかの組を使用することによって決定された、予め決定された行列である。このデータから、カールフネン−ロエベ(K−L)変換のような変換行列を、既知の方法において発生することができる。前記変換行列を前記局所ベクトル列に用い、特徴ベクトル列を発生することができる。特徴ベクトル空間領域において、ベクトル量子化技術を使用し、前記特徴ベクトル列を分類することができる。
【0074】
解析的自己適合的アルゴリズムを、前記特徴ベクトルの分類に使用することができる。このアルゴリズムの規定において、{Xi∈R4:i=1,2,3,...,N}を前記特徴ベクトルの列とし、Nを特徴ベクトルの数とし、Kはクラスの最大数を示し、Tを前記データ組に適合するしきい値とする。各クラスに関して、代表的な要素を前記アルゴリズムによって発生する。akをクラスkの代表的な要素とし、nkをクラスkにおける特徴ベクトルの数とする。
【0075】
次に、前記アルゴリズムを以下のように概説することができる。
1.n1=1;a1=X1;
【外10】
=1とおく
2.クラス数[外10]およびクラスパラメータ(ak,nk)を得る
(i=1;i<N;i++)と(j=1;j<[外10];j++)とに関して、
dj=dist(Xi,aj)を計算し、
index=ark min djとして、
(dindex<T)jまたは([外10]<K)ならば、
クラスパラメータを更新する
aindex=(1/(nindex+1))×(nindex・aindex+Xi)
nindex=nindex+1
(dindex<T)jまたは([外10]<K)でなければ、
新たなクラスを発生する
【数15】
【数16】
[外10]=[外10]+1
3.各々の特徴ベクトルを、最も近い隣接ルールにしたがって、クラスにラベル付けする
(i=1;i<N;i++)と(j=1;j<[外10];j++)とに関して、
dj=dist(Xi,aj)を計算し、
index=ark min dj
ボクセルiをクラスインデックスにラベル付けする
【0076】
このアルゴリズムにおいて、dist(x,y)は、ベクトルxとyとの間のユークリッド距離であり、arc min djは、djの最小値を実現する整数である。
【0077】
上述したアルゴリズムは、パラメータTおよびKにのみ依存する。しかしながら、各ボクセルクラスタ内のクラス数に関するKの値は重要ではなく、K=18のような一定値に設定することができる。しかしながら、ベクトル類似性しきい値であるTは、分類結果に大きく影響する。Tの選択された値が大きすぎる場合、1つのクラスのみが発生する。他方において、Tの値が小さすぎる場合、結果として生じるクラスは、望ましくない冗長性を示す。Tの値を、前記特徴ベクトル列の最大の構成要素の変化に等しく設定することによって、最大数の別個のクラスが結果として生じる。
【0078】
前記初期分類処理の結果として、前記選択されたクラスタ内の各ボクセルを1つのクラスに割り当てる(ステップ1570)。仮想大腸鏡検査の好例の場合において、クラスタ4内にいくつかのクラスが存在する。したがって、次のタスクは、クラスタ4におけるいくつかのクラスのうちのどれが前記結腸壁に対応するかを決定することである。前記特徴ベクトルの残りの座標は、前記局所隣接ボクセル内の方向強度変化の情報を含む。前記結腸の内部に関する結腸壁ボクセルは、一般に、クラスタ1のガスボクセルの近くにあるため、しきい値間隔を、代表的なCTデータ組の代表的な結腸壁強度から選択されたデータ標本によって決定し、結腸壁ボクセル候補を識別することができる。特定のしきい値を、各々の特定の画像プロトコルおよび装置に関して選択する。次に、このしきい値間隔を、(同じ機械から、同じ画像プロトコルを使用して取得された)すべてのCTデータ組に用いることができる。前記代表的な要素の第1座標が前記しきい値間隔において位置する場合、対応するクラスを前記結腸壁クラスとみなし、このクラスにおけるすべてのボクセルに結腸壁様ボクセルとラベル付けする。
【0079】
各々の結腸壁様ボクセルは、コロン壁ボクセルである候補である。結腸壁に属さない3つの可能な結果が存在する。第1の場合は、前記結腸内部の大便/流体に近いボクセルに関係する。第2の場合は、ボクセルが肺組織領域にある場合に生じる。第3の場合は、粘膜ボクセルを表す。明らかにこのとき、低レベル分類は、ある程度の分類不確実性を伴う。前記低レベル分類不確実性の原因はさまざまである。例えば、2つ以上の材料形式(すなわち、流体および結腸壁)を含むボクセルによる結果として生じる部分的体積効果が、不確実性の第1の場合を導く。不確実性の第2および第3の場合は、前記部分的体積効果と、前記CT画像が低コントラストなこととの双方による。前記不確実性を除去するために、追加の情報が必要である。したがって、高レベル特徴抽出プロセスを本方法において使用し、前記結腸壁に関するボクセルを他の結腸壁様ボクセルから、前記CT画像の先験的な解剖知識に基づいてさらに識別する(ステップ1580)。
【0080】
前記高レベル特徴抽出手順の最初のステップを、前記低レベル分類結果から肺組織の領域を除去することとしてもよい。図18Aは、肺領域1802を明確に説明する好例のスライス画像である。肺領域1802は、一般的に、図18Bにおいて説明するような結腸壁様ボクセルによって囲まれた接触する3次元体積として識別可能である。この特性が与えられれば、この技術における第1ステップは、前記腫瘍の領域内で種ボクセルを見つけることである。好適には、前記CT画像操作を行うオペレータは、画像範囲を設定し、前記CT走査の最も上部のスライスが結腸ボクセルを含まないようにする。肺の内部が空気で満たされるべきとき、前記種を、前記低レベル分類によって、単に空気ボクセルを選択することによって与える。図18Bの肺領域輪郭が決定されたら、前記肺体積を前記画像スライスから除去することができる(図18C)。
【0081】
高レベル特徴抽出を行うことにおける次のステップを、クラスタ2における強調された大便/流体ボクセルから骨ボクセルを分離することとすることができる。骨組織ボクセル1902は、一般に、前記結腸壁から比較的離れており、前記結腸体積の外部に存在する。それと反対に、残りの大便1906および流体1904は、前記結腸体積内に囲まれる。先験的な解剖知識と、前記低レベル分類プロセスから得られた結腸壁情報とを組み合わせ、大まかな結腸壁体積を発生する。予め決められた数(例えば、3)より多いボクセルユニットによって前記結腸壁から分離され、前記結腸体積の外部のどのベクトルにも、骨とラベル付けし、次に、前記画像から除去する。クラスタ2における残りのボクセルは、前記結腸体積内の大便および液体を表すと仮定することができる(図19A−C参照)。
【0082】
大便1906および流体1904として識別される前記結腸体積内のボクセルを、前記画像から除去し、きれいな結腸管腔および結腸壁画像を発生することができる。一般に、2種類の大便/流体領域が存在する。一方の領域形式は、前記結腸壁に付着した大便1906の残りの領域である。他方の領域形式は、流体1904の大きい体積であり、洗面器状の結腸の折りたたみにおいて集まる(図19A−C参照)。
【0083】
付着した残りの大便領域1906を、これらは前記低レベル分類プロセス中に発生された大まかな結腸体積の内部にあるため、識別し、除去することができる。前記洗面器状の結腸の折りたたみにおける流体1906は、通常、重力の影響により水平面1908を有する。前記水平面の上は通常ガス領域であり、このガス領域は前記流体強度に対してきわめて高いコントラストを示す。したがって、前記流体領域の界面を容易に明示することができる。
【0084】
領域成長戦略を使用し、付着した大便領域1906の輪郭を示すことができ、前記結腸壁体積から離れた部分を除去することができる。同様に、流体領域1904の輪郭も示すことができる。水平面1908を除去した後、前記結腸壁輪郭も除去し、きれいな結腸壁を得る。
【0085】
前記結腸壁ボクセルから粘膜ボクセルを識別することは困難である。上記3次元処理がいくつかの粘膜ボクセルを除去できたとしても、すべての粘膜ボクセルを除去するのは困難である。光学大腸鏡検査において、医師は、結腸粘膜を直接検査し、前記粘膜の色および感触に基づいて病変を探す。仮想大腸鏡検査において、前記結腸壁における大部分の粘膜ボクセルを、より多くの情報を保存するために、残すことができる。これは、3次元体積レンダリングに関してきわめて有用であるかもしれない。
【0086】
前記分割された結腸壁体積から、前記結腸の内部表面、外部表面および壁自体を抽出し、仮想物体として見ることができる。これは、前記結腸の外部壁を内部壁と同様に調査することができるという、慣例的な光学大腸鏡検査に対する独自の利点を提供する。さらに、前記結腸壁および結腸管腔を、前記分割から別々に得ることができる。
【0087】
前記結腸を、実質上画像化の前に空にするため、一般に遭遇する問題は、前記結腸管腔がある程度までしぼむことである。空気またはCO2のような圧縮ガスによる前記結腸の膨張は、しぼんだ領域の頻度を減らす。仮想大腸鏡検査を行うことにおいて、前記しぼんだ領域を通じて飛行経路を自動的に保持することが望ましく、前記走査された画像データを使用し、前記しぼんだ領域において前記結腸管腔を少なくとも部分的に再形成することも望ましい。上述した画像分割方法は、前記結腸壁の内部および外部の双方を有効に得るため、この情報を使用し、前記しぼんだ領域を通じて飛行経路の発生を拡張することができる。前記結腸のしぼんだ領域を通じた前記飛行経路の拡張、または、前記結腸のしぼんだ領域の膨張において、第1ステップは、しぼんだ領域を検出することである。前記結腸壁の外側の周囲からの画像データのグレイスケール値は、前記結腸壁それ自体の中や、脂肪、筋肉または他の種類の組織のような他の領域におけるグレイスケール値より劇的により大幅に変化するという前提を使用し、エントロピ解析を使用して、結腸しぼみの領域を検出することができる。
【0088】
例えば中心線に沿ったグレイスケール値の変化の程度を、エントロピ値によって表し、測定することができる。エントロピ値を計算するために、前記結腸の外側表面におけるボクセルを選択する。これらのような点を、上述した画像分割技術から識別する。5×5×5の立方体ウィンドウを、興味ある画素を中心として前記画素に用いる。前記エントロピ値を計算する前に、前記画像データからノイズをろ過するために、より小さい(3×3×3)ウィンドウを興味ある画素に用いることができる。次に、前記画素について選択されたウィンドウのエントロピ値を、以下の式によって決定することができる。
【数17】
ここで、Eはエントロピ値であり、C(i)は前記ウィンドウ内のi(i=0,1,2,...,255)のグレイスケール値を有する点の数である。次に、各々のウィンドウに関して計算されたエントロピ値を、予め決められたしきい値と比較する。空気の領域に関して、前記エントロピ値は、組織の領域と比べると、かなり低くなる。したがって、前記結腸管腔の中心線に沿って、前記エントロピ値が増加し、予め決められたしきい値を超える場合、しぼんだ領域を示す。前記しきい値の正確な値は重要ではなく、画像プロトコルと、画像装置の詳細とに部分的に依存する。
【0089】
しぼんだ領域が検出されたら、前もって決定された中心線飛行経路を、1ボクセル幅ナビゲーションラインでしぼみの中心を貫通させることによって前記領域を通じて拡張することができる。
【0090】
前記仮想カメラの飛行経路を前記結腸管腔を通じて自動的に延長することに加えて、前記結腸しぼみの領域を、物理的モデル化技術を使用して仮想的に広げ、前記しぼんだ領域の特性のうちいくつかを回復させることができる。この技術において、前記結腸壁の物理的特性のモデルを開発する。このモデルに関して、運動パラメータ、質量密度、ダンピング密度、引張係数、および曲げ係数を、ラグランジュ方程式に関して推定する。次に、拡張力モデル(すなわち、前記結腸中に注入された空気のようなガスまたは流体)を、前記ラグランジュ方程式によって規定されるような前記結腸の弾性特性に従って式化し、適用し、前記結腸画像のしぼんだ領域をその自然な形状に復元するようにする。
【0091】
前記結腸をモデル化するために、有限要素モデルを、前記結腸管腔のしぼんでふさがれた領域に用いる。これを、8ボクセルレンガのような一定のグリッドにおいて前記要素を標本化し、次に、慣例的な体積レンダリング技術を用いることによって行うことができる。代わりに、四面体のような不規則体積表現アプローチを前記しぼんだ領域に用いることができる。
【0092】
外力(空気注入)モデルを前記結腸モデルに用いることにおいて、前記外力の大きさを最初に決定し、前記しぼんだ結腸壁領域を適切に分離する。三次元成長モデルを使用し、結腸内壁表面および結腸外壁表面を並列にトレースすることができる。個々の表面を、前記しぼんだ領域における開始点から成長源点までマークし、前記外力モデルを用い、前記表面を同様かつ自然なように膨張させる。前記内部表面および外部表面の間の領域、すなわち、結腸壁を、共有領域として分類する。外部斥力モデルをこれらの共有領域に用い、前記しぼんだ結腸壁部分を自然なように分離し、膨張させる。
【0093】
仮想調査を受ける結腸のような仮想物体の特徴をより明瞭に視覚化するために、前記物体の種々のテクスチャのレンダリングを与えることが有利である。光学大腸鏡検査中に与えられるカラー画像において観察することができるこれらのようなテクスチャは、CT画像データによって与えられる白黒、グレイスケール画像においてしばしば失われる。したがって、仮想調査中にテクスチャ画像化するシステムおよび方法が要求される。
【0094】
図20は、テクスチャ成分を有する仮想画像を発生する本方法を示すフローチャートである。本方法の目的は、例えばビジブルヒューマンからの赤−緑−青(RGB)における光学大腸鏡検査画像によって得られたテクスチャを、仮想物体を発生するのに使用されるグレイスケールモノクロームCT画像データ上にマッピングすることである。前記光学大腸鏡検査画像を、ビデオカメラのようなカメラからアナログ光学画像を受け、前記画像を、インタフェースポート1431(図14)を経てCPU1423に供給することができるディジタルデータに変換するディジタル「フレームグラバー」1429によるような、慣例的なディジタル画像取得技術によって取得する。このプロセスにおける第1ステップは、前記CT画像データを分割することである(ステップ2010)。上述した画像分割技術を用い、前記グレイスケール画像における強度しきい値を選択し、前記CT画像データを、骨、結腸壁組織、空気等のような種々の組織形式に分類することができる。
【0095】
前記CT画像データにおける画像分割を行うことに加えて、前記光学画像のテクスチャ特徴を、前記光学画像データから抽出することが必要である(ステップ2020)。これを行うために、ガウシアンフィルタを前記光学画像データに用いることができ、続いて、前記データをサブ標本化し、多解像度ピラミッドに分解する。ラプラシアンフィルタおよび可変フィルタを前記多解像度ピラミッドに用い、前記データの指向特徴および非指向特徴を得ることもできる。この方法は、前記テクスチャ特徴を抽出し、獲得することにおいて効率的であるが、このアプローチの実施は、大量のメモリおよび処理パワーを必要とする。
【0096】
前記テクスチャ特徴を前記光学画像から抽出する代わりのアプローチは、ウェーブレット変換を使用することである。しかしながら、ウェーブレット変換は、一般に計算上効率的であるが、慣例的なウェーブレット変換は、軸に平行な方向を持つ特徴のみを獲得することに制限され、興味ある領域に直接用いることができない。これらの制限を克服するために、分離不可フィルタを用いることができる。例えば、リフティングスキームを用い、2つのステップ、予測および更新アプローチを使用して、任意の次元におけるウェーブレット変換に関するフィルタバンクを構成することができる。これらのようなフィルタバンクを、多次元多項式補間に関するブーアロム(Boor−Rom)アルゴリズムによって合成することができる。
【0097】
前記テクスチャ特徴を前記光学画像データから抽出した後、モデルを発生し、これらの特徴を描かなくてはならない(ステップ2030)。これを、例えば、自然のテクスチャに起因する非ガウス分布のエントロピを推定し、操作することを基礎とする非パラメトリック多スケール統計モデルを使用することによって行うことができる。
【0098】
テクスチャモデルを前記光学画像データから発生したら、テクスチャマッチングを行い、これらのモデルを前記分割されたCT画像データに関連させなければならない(ステップ2040)。前記テクスチャが連続的な前記CT画像データの領域において、対応するテクスチャのクラスは容易に調和する。しかしながら、2つ以上のテクスチャ領域の間の境界領域において、前記プロセスはより複雑になる。境界領域の周囲の前記CT画像の分割は、しばしば、ぼけたデータを生じ、すなわち、各々の材料または組織からのテクスチャの割合を反映すると共に各々の種々の重み付けに応じて変化する結果を生じる。前記重みの割合を使用し、マッチング基準の重要性を設定することができる。
【0099】
前記非パラメトリック多スケール統計モデルの場合において、クロスエントロピまたはカルバックライブラー発散アルゴリズムを使用し、境界領域における異なったテクスチャの分布を測定することができる。
【0100】
テクスチャマッチング後、テクスチャ合成を前記CT画像データにおいて行う(ステップ2050)。これを、前記光学画像データからのテクスチャを前記CT画像データに融合させることによって行う。骨によって与えられるような等方性「テクスチャパターンに関して、前記テクスチャを前記光学データから前記CT画像データへ直接標本化することができる。結腸粘膜のような異方性テクスチャ領域に関して、多解像度標本化手順が好ましい。このプロセスにおいて、均質領域および不均質領域に関する選択性再標本化を用いる。
【0101】
代わりに、擬似カラーテクスチャを、前記CTデータから直接形成することができる。各ボクセルに関し、局所領域隣接部を具える多数のCT値を評価し、所定のボクセルに関する擬似カラーを決定することができる。例えば、5×5×5ボクセル立方体状領域、または、マンハッタン距離によって測定された3ユニット内のすべてのボクセルを表す二重ピラミッド。次に、このスカラ値のベクトルを評価し、その後の体積レンダリング中にこのボクセルに関して表示すべきカラーにマッピングする。前記値の局所隣接部ベクトルの評価は、これらのようなことを、局所曲率、均質性/不均質性、または他の幾何学的または空間的関数として計算することができる。
【0102】
体積レンダリング
上述した画像分割およびテクスチャマッピングに加えて、体積レンダリング技術を仮想大腸鏡検査手順と組み合わせて使用し、結果として生じる画像の忠実度をさらに増す。図21は、本発明による体積レンダリングに使用することができるパースペクティブ体積光線投射法を説明する。前記結腸管腔内の選択された仮想始点、例えばカメラ位置から、光線を、近似画像画素の各々を経て投影する(ステップ2100)。各々の光線に関して、第1標本化点を、前記光線に沿った現在画像画素として設定する(ステップ2110)。次に、前記現在標本化点と最も近い結腸壁との間の距離(d)を決定する(ステップ2120)。源z内距離(d)を、予め決められた標本化間隔(i)と比較する(ステップ2130)。距離(d)が標本化間隔(i)より大きい場合、標本化を行わず、前記光線に沿った次の標本化点を、前記光線に沿って距離dジャンプすることによって決定する(ステップ2140)。前記距離が標本化間隔(i)以下である場合、慣例的な標本化をこの点において行い(ステップ2150)、次の標本化点を標本化間隔(i)にしたがって選択する(ステップ2160)。例えば、8個の隣接するボクセルの密度値の間のトリリニア補間を行い、前記標本化点における新たな密度値を決定することができる。
【0103】
図21の方法は、空間跳躍技術を使用し、前記画像面の光線に沿って空の空間を越えて前記結腸壁まで迅速にスキップするため、効率的に光線投影を加速させる。この方法において、標本点から最も近い結腸壁までの距離を、各々の光線に沿って決定する。前記距離が予め決められた間隔(i)より大きい場合、前記光線に沿って次の標本化点までのジャンプを行う。前記最も近い距離情報は、仮想カメラ制御に使用される前記ポテンシャル場からすでに利用可能であるため、追加の距離符号化計算は必要ない。この場合において、表面レンダリングもZバッファ計算も必要なく、結果として、予備的処理時間およびメモリ空間が節約される。
【0104】
代わりに、空間跳躍法は、各々の光線に関する距離情報を、対応する表面レンダリング画像のZバッファから得ることができる。前記表面レンダリング画像および体積レンダリング画像の双方を発生する場合、このアプローチは、前記Zバッファ情報を前記表面レンダリングの結果として与えるため、最小の処理オーバヘッド負担を与える。したがって、この形式の空間跳躍方法は、前記画像空間領域から世界空間領域への深さ変換を行う追加の処理のみを必要とする。
【0105】
距離(d)をステップ2140において横切った前記光線に沿ったこれらの領域に関して、前記光線に沿った領域は、オープンスペースであり、オープンスペース伝達関数に従って値に割り当てることができる。代表的に、オープンスペースは、最終画素値において寄与しない。標本化を行う各々の点に関して、1つ以上の規定された伝達関数を割り当て、元の体積データの異なった範囲の標本値を、異なった色および不透明度と、あるいは他の表示化のパラメータとにマッピングする。例えば、4つの独立した伝達関数を使用し、CT密度値の範囲を、赤、緑、青の指定された色と、不透明度とに、各々0から255の範囲においてマッピングすることによって、異なった材料を決定した。
【0106】
仮想生体組織検査
上述した技術は、柔軟かつ非侵入的生体組織検査を実行するために調査されている領域の仮想電子生体組織検査を行うシステムの基礎を形成することもできる。上記で注意したように、体積レンダリング技術は、1つ以上の規定された伝達関数を使用し、元の体積データの異なった範囲の標本値を、異なった色、不透明度、および、ナビゲーションおよび見ることに関する他の表示可能パラメータにマッピングする。ナビゲーション中、選択された伝達関数は、一般に、最大不透明度を前記結腸壁に割り当て、外側表面が容易に見られるようにする。疑わしい領域を仮想調査中に検出したら、医師は、体積レンダリング手順中に割り当てられた伝達関数をインタラクティブに変更し、結果として、見られている外側表面が実質上透明になり、前記領域の情報が合成され、したがって、前記領域の内部構造が見られるようにする。多数の予め決定された伝達関数を使用し、前記疑わしい領域を、前記プロセスを通じて割り当てられた不透明度を変化させることによって、多数の異なった深さにおいて見ることができる。
【0107】
ポリープ検出
本システムおよび方法を使用し、自動化されたポリープ検出を行うことができる。図13の参照と共に、例えば前記結腸内に生じるポリープ1305は、一般に、結腸壁1301から延在する小さな凸面の丘のような構造の形態をとる。このジオメトリは、前記結腸壁の折りたたみと異なる。したがって、異なったジオメトリモデルを使用し、前記結腸壁におけるこれらのようなポリープを検出することができる。
【0108】
前記結腸管腔の表面を、C−2平滑表面モデルを使用するように、3次元ユークリッド空間における連続的な第2の区別可能な表面として表すことができる。このようなモデルは、その全体の参照によってここに含まれる、1994年、スプリンガー−バーラグ(Springer−Verlag)によって出版された、B.A.デュブロビン(Dubrovin)による「現代のジオメトリ方法および用途」に記載されている。このモデルにおいて、前記結腸の表面における各ボクセルは、ガウス曲率フィールドと呼ばれるガウス曲率を有する関連する幾何学的特長を有する。ポリープを示すかもしれない前記表面における凸面の丘は、前記ガウス曲率フィールドにおいてユニークな局部特徴を備える。したがって、前記ガウス曲率フィールドを特定の局部特徴に関して探索することによって、ポリープを検出することができる。検出されたら、疑わしいポリープを強調し、したがって、医師の注意を引かせることができ、医師は前記疑わしいポリープを測定し、上述した仮想生体組織検査方法を使用し、前記疑わしい領域をさらに調査することができる。
【0109】
中心飛行経路発生
仮想大腸鏡検査の場合において、前記結腸管腔を通る適切なナビゲーションライン、または、飛行経路は、上述したシステムおよび方法の重要な特徴である。前記仮想カメラモデルの飛行経路を決定するいくつかの技術を、図4−8に関して考察したが、図22は、前記結腸管腔を通る中心飛行経路を発生する代わりの方法を示す。前記結腸壁を識別した後、上述した画像分割方法等によって、体積減少アルゴリズムを用い、前記結腸管腔の方向を強調し、前記結腸体積内のその後の探索時間を減少することができる(ステップ2310)。
【0110】
図23は、多解像度解析モデルに基づく好例の体積減少アルゴリズムのステップをさらに示す。この手順において、前記3次元体積を、同じ行列サイズを有するバイナリ画像のスタックによってあらわす(ステップ2310)。集合的に、これらの画像はバイナリデータ組を形成する。離散的ウェーブレット変換を前記バイナリデータ組に用いることができ、これは、前記バイナリデータ組の異なった時間−周波数成分を表す多数のサブデータ組を生じる。例えば、離散的ウェーブレット変換は、8つのサブデータ組を生じるかもしれない。前記サブデータ組を、予め決められたしきい値と比較し、最も低い周波数成分を識別する(ステップ2330)。この成分は、多解像度構造において再帰的に用いるその後の離散的ウェーブレット変換およびしきい値処理ステップに関するバイナリデータ組を形成する(ステップ2340)。仮想大腸鏡検査の場合において、離散的ウェーブレット変換および関連するしきい値処理を、前記最も低い周波数成分を表すその後のサブデータ組において3回再帰的に用いる(3レベル多解像度分解)。
【0111】
図22に戻って、前記減少した結腸体積モデルから、距離マップ技術を用い、前記結腸の2つの端の間、例えば、直腸から盲腸までの最短距離経路を発生することができる(ステップ2215)。この結果として生じる経路は、前記結腸管腔の全体的な方向情報を維持するが、局所的折りたたみによって示される方向を無視する。次に、全体的な結腸内における制御点を、前記最短距離経路を元のデータ空間にマッピングし戻すことによって得ることができる(ステップ2220)。例えば、3レベル多解像度分解の場合において、減少した体積は、元の体積の3分の1であり、既知のアフィン変換を使用し、前記減少した体積モデルを正確に原スケールモデルにマッピングし戻すことができる。前記減少した体積の最短距離経路を、前記原スケール体積中に点の列としてマッピングし戻すこともでき、これらの点を前記結腸内の制御点として使用することができる。
【0112】
好適な飛行経路は、前記結腸管腔の中心線上におけるものである。しかしながら、初期制御点は、前記結腸管腔の中心において正確に位置しないかもしれない。したがって、前記初期制御点を、2等分面アルゴリズムの使用等によって、中心におく(ステップ2230)。例えば、各々の選択された制御点において、2等分面を、前記方向に対して垂直で前記結腸管腔を横切って切断する面として規定することができる。次に、最大ディスクアルゴリズムのような中心化アルゴリズムを各々の2等分面において行うことができる。このようなアルゴリズムは、参照によってここに含まれる、1996年、PAMIのIEEE会報、Vol.18、1055−1066ページ、ゲ(Ge)他による論説「別々のユークリッド距離マップからのスケルトンの発生における」において考察されている。
【0113】
前記制御点を中心化したら、前記飛行経路を、これらの点を接続する線を補間することによって決定することができる(ステップ2240)。仮想大腸鏡検査の場合において、前記補間された飛行経路は、前記結腸管腔内においてほぼ中心にある滑らかな曲線の形状をとることが望ましい。1992年、ケンブリッジ大学出版、プレス他による「Cにおける数値レシピ:科学計算の技術」において記載のように、微分幾何学理論におけるセレー−フルネ定理に基づく限定されたキュービックBスプライン補間アルゴリズムを使用し、好適な滑らかにカーブする飛行経路を確立することができる。
【0114】
図24における分割された結腸管腔の絵的表現と、図25のフローチャートとは、本発明による依然として他の代わりの飛行経路発生方法を説明する。この代わりの方法において、結腸管腔2400の表現を、最初に、管腔2400の長さに沿って多数のセグメント2402に分割する(ステップ2500)。各々のセグメント2402内から、代表的な点2404a−gを選択する(ステップ2520)。次に、各々の代表的な点2404a−gを、個々のセグメントの中心にこれらの点を押すのに使用される物理ベース変形モデルの使用等によって、前記結腸壁に関して中心にする(ステップ2530)。前記代表的な点を中心にした後、前記点を順次結合し、前記仮想カメラモデルに関する中心線飛行経路を確立する(ステップ2540)。前記セグメントが長さにおいて十分に小さい場合、前記中心化した点を、直線セグメント2406a−fで接続することができる。しかしながら、リニアカーブフィッティング技術を用い、前記中心化した点を結合する場合、より滑らかで連続的な飛行経路が確立される。
【0115】
上述した方法の各々を、図14に示すようなシステムを使用し、適切なソフトウェアを与え、CPU1409およびCPU1423の動作を制御することによって実現することができる。
【0116】
パーソナルコンピュータにおける配備に好適な代わりのハードウェア実施形態を図26に示す。本システムは、400MHzを超えるクロック速度において動作するペンティアム(登録商標)IIIプロセッサのような高速マルチタスクプロセッサの形態をとるべきプロセッサ2600を含む。プロセッサ2600を、高速並列データ伝送を与える慣例的なバス構造2620に結合する。バス構造2620には、メインメモリ2630と、グラフィックスボード2640と、体積レンダリングボード2650も接続する。グラフィックスボード2640を、好適には、ダイアモンドマルチメディアシステムズによって製造されるダイアモンドバイパーv770ウルトラのような、テクスチャマッピングを行うことができるものとする。堆積レンダリングボード2650は、参照によってここに含まれる米国特許第5760781号および第5847711号に基づく、ミツビシエレクトリックによるボリュームプロボードの形態をとることができる。慣例的なSVGAまたはRGBモニタのような表示装置2645を、前記画像データを表示するグラフィックスボード2640に効率的に結合する。MRIまたはCTスキャナのようなイメージングスキャナからデータを受け、このようなデータをバス構造2620に伝送するスキャナインタフェースボード2660も設ける。スキャナインタフェースボード2660は、選択されたイメージングスキャナに関する用途特定インタフェース製品であってもよく、または、一般的な目的の入力/出力カードの形態をとることもできる。PCベースシステム2600は、一般に、キーボード、ディジタルポインタ(例えば、マウス)等のようなI/O装置2680をプロセッサ2620に結合するI/Oインタフェース2670を含む。代わりに、前記I/Oインタフェースを、プロセッサ2620にバス2620を介して結合することができる。
【0117】
3次元画像化の場合において、テクスチャ合成、体積レンダリング、多数のデータ処理および処理動作を含むことが必要である。前記結腸管腔およびその周囲の領域によって表されるもののような大きいデータセットに関して、このような処理は、きわめて時間を消費し、メモリが増大するかもしれない。しかしながら、図27のフローチャートにおいて説明した処理方法にしたがって図26のトポロジーを使用することによって、これらのような動作を、比較的低コストのパーソナルコンピュータ(PC)において行うことができる。画像化データを、スキャナインタフェースボード2660およびバス構造2620を介してプロセッサ2620から受け、メインメモリ2630に格納する。この画像データ(画素)を体積要素(ボクセル)表現に変換する(ステップ2710)。メインメモリ2630に格納された前記体積表現を、例えば、主要な体積軸に沿うようなスライス、または、画像化されている領域の他の部分に分割する(ステップ2720)。前記体積分割を、次に、前記体積レンダリングボードに伝送し、体積レンダリング動作に関する体積レンダリングメモリ2655に一時的に格納する(ステップ2730)。局所的に存在する体積レンダリングメモリ2655の使用は、全体の体積の各々のスライスのレンダリング中、バス2620を越えてデータを交換する必要がないため、体積レンダリングにおける速度を上昇させる。体積レンダリングが前記スライスに関して完了したら、前記レンダリングされたデータを、シーケンシャルバッファにおいてメインメモリ2630またはグラフィックスボード2640に伝送し戻す(ステップ2740)。興味あるすべてのスライスをレンダリングした後、前記シーケンシャルバッファにおける内容を、表示ユニット2645における表示に関してグラフィックスボード2640によって処理する(ステップ2750)。
【0118】
多走査ベース仮想調査
上記で考察した技術は、一般に、単一の磁気共鳴画像(MRI)またはコンピュータ断層撮影(CT)走査から取得したデータセットに基づいて仮想画像化を行う。しかしながら、上記で考察した技術は、ある領域の多数の走査を使用してある領域の仮想調査を行うのにも有用である。ある領域の多数の走査を使用することによって、異常な領域の改善された画像化を達成することができ、モーションアーティファクトを減少することができる。ある興味深いこのような用途は、仮想大腸鏡検査を行い、あり得る膀胱のポリープまたは癌に関して患者を選別することにある。
【0119】
図28は、多数のMRI走査を行い、仮想大腸鏡検査のようなある物体の仮想調査を行う方法を説明するフローチャートである。CT画像と異なり、膀胱壁を尿と識別することが困難であるかもしれない場合、MRI画像において、尿を自然なコントラスト剤として使用し、前記膀胱内壁の輪郭を描くことができる。この目的のため、前画像走査プロトコルを用いる(ステップ2805)。4回のMRI走査の第1のものの約1/2時間前、患者に膀胱を空にし、次に1カップの水を飲むように要求する。約1/2時間後、患者に前記膀胱領域の4回の走査のうちの第1回目を行う(ステップ2810)。前記膀胱がいっぱいで膨らんだ第1走査は、T1重み付き横画像化に関するプロトコルに従う。例えば、上記で参照したピッカー(Picker)スキャナを使用する場合、256×256行列サイズ、38cm視野(FOV)、1.5mmスライス厚(隙間なし)、3ミリ秒TE、9ミリ秒TR、30度フリップ角、および、1走査平均を使用するクジェルファスター(KJELL FASTER)プロトコルを使用することができる。もちろん、これらのパラメータはスキャナ特定になる傾向にあり、前記パラメータにおける種々の変更を、許容しうる結果に関して使用することができる。
【0120】
前記膀胱がいっぱいのままで、前記患者に第2MRI走査、走査2を行う(ステップ2815)。前記第2走査は、256×256行列サイズ、38cm視野(FOV)、1.5mmスライス厚(隙間なし)、3ミリ秒TE、9ミリ秒TR、30度フリップ角、および、2走査平均を使用するピッカークジェルファスター(Picker KJELL FASTER)プロトコルのような、T1重み付き冠状画像化に関するプロトコルに従う。
【0121】
上述した2回の画像走査を、互いに直交する軸に沿って行う。この利点は、一方の走査における重大なモーションアーティファクトの領域が、一般に、直交する走査における最小のモーションアーティファクトの領域に対応することである。したがって、第1走査を横方向に、第2走査を冠状方向にとることによって、前記画像走査を登録することができ、前記データセットにおけるモーションアーティファクトを識別し、補償することができる。
【0122】
走査2の後、前記患者に、前記膀胱を開放するように頼み、次に、2回の追加のMRI走査を行う。第3走査(ステップ2820)は、前記第1走査(横画像化)と同じ画像化プロトコルに従う。第4走査(ステップ2825)は、前記第2走査(冠状画像化)と同じ画像化プロトコルに従う。
【0123】
前記画像走査を、ピッカー1.5Tエッジ全身スキャナを使用して取得することができる。T2画像化プロトコルを使用することもできるが、T1画像化プロトコルは、脂肪と尿との間の改善された輪郭描写を与え、より短い取得期間しか必要としないため、T1画像化プロトコルが仮想大腸鏡検査に好適である。代わりに、前記画像操作は、それらの画像化技術に好適なコントラスト剤およびプロトコルを使用するコンピュータ断層撮影または超音波画像走査の形式をとってもよい。
【0124】
前記最初の2回の走査(走査1および2)の間、前記膀胱は膨張し、前記膨張壁は比較的薄い。この場合において、腫瘍のような生理学上変化した場所は、変化しない膀胱壁と比べて異なった率において薄くなるかもしれず、これらの状況の下でより明らかになるかもしれない。前記第3および第4走査中、前記膀胱は実質的に空であり、前記膀胱壁はより厚い。より厚い壁により、前記膀胱壁の通常の組織と生理学上変化した組織との間に、より明白な画像コントラストが結果として生じるかもしれない。
【0125】
前記4回の走査を取得した後、4つの対応するデータセットを別々に処理する。最初に、図15と組み合わせるように上述したように(ステップ2830)、各々のデータセットを、好適には、画像分割する。画像分割中、前記4つのデータセットのボクセルを、膀胱壁、尿、脂肪、境界等の多数のカテゴリに分類する。前記分類は、前記ボクセルの局所強度ベクトルに基づく。前記ボクセルを分類したら、前記膀胱管腔の内部を、空気ボクセルまたは尿ボクセルを選択することによるような、前記膀胱体積内から選択された種ボクセルで開始する領域成長アルゴリズムを使用して識別することができる。
【0126】
前記分割された体積データセットの臨床解析の前に、前記4つのデータセットの共通座標系への登録を行う(ステップ2835)。前記膀胱の形状は、走査ごとに変化し、正確なボクセル−ボクセル登録は、実際上の値のものではない。代わりに、柔軟な登録プロセスが好適である。この柔軟な登録プロセスにおいて、興味ある各々の体積(各々の対応する走査に関してレンダリングされた体積)に関して、前記体積の中心を、前記体積におけるすべてのボクセルの3つの座標を平均すること等によって決定する。
【0127】
次に、デカルト座標系を構成し、この座標系の原点を前記体積の中心点において位置させる。次に、前記座標系の軸を、多数の方法において方向付けることができる。方向付けの好適な選択は、自然な人間の身体の方向に対応し、例えば、Z軸が前記身体の高さに沿って(例えば、つま先から頭まで)通り、Y軸が後ろから前まで通り、X軸が横に(例えば、左から右に)通る。この座標系における長さの単位を、1ボクセル長の任意の単位に対して便利に設定することができ、その絶対値は、MRI走査に関する取得特性に基づいて変化する。同じ画素間隔をすべての走査において使用し、4つのデータセットのすべてを獲得する限りは、前記4つのデータセットの各々に関して一様な値が結果として生じる。
【0128】
登録後、前記4つのデータセットからの画像を、別々にまたは同時に見ることができる(ステップ2845)。好例の表示ウィンドウを図29および30に示す。図29を参照し、前記表示を、各々走査1、走査2、走査3および走査4に対応する4つのサブウィンドウ2905、2910、2915、2920に分割する。制御パネル選択2925を、前記表示のある部分において与え、グラフィカルユーザインタフェース(GUI)を確立し、表示およびナビゲーション機能をユーザに与えることもできる。オペレータが前記画像サブウィンドウの1つにおいてビューを拡大するようなナビゲーションをすると、好適には、対応する動作が他のサブウィンドウビューにおいても同様に起こる。ユーザは、単一ウィンドウ表示への拡張に関して前記ビューのうち1つを選択することもできる。
【0129】
同時に処理されるデータの量を減らすために、前記データセットを、8つの部分またはオクタント等に分割することができる(ステップ2840)。これを、多数の方法において行うことができる。例えば、図29に示すデカルト座標系への参照と共に、前記データを前記座標系の8つの領域、(1)X,Y,Z;(2)X,−Y,Z;(3)X,Y,−Z;(4)X,−Y,−Z;(5)−X,Y,Z;(6)−X,−Y,Z;(7)−X,Y,−Z;および(8)−X,−Y,−Zに分割することができる。
【0130】
図29は、前記4回の走査の各々から得られた前記膀胱管腔の外側の4つのビューを示す。図30は、前記4回の走査の各々から得られた前記膀胱管腔の内側の一部の4つのビューを示す。
【0131】
多解像度画像化および仮想喉頭鏡検査
ここで説明したシステムおよび方法を適合させ、用いて、仮想喉頭鏡検査に好適な多解像度画像化を行うことができる。図31は、仮想喉頭鏡検査を行う方法を説明するフローチャートである。第1に、患者の喉頭の領域の画像走査を取得する(ステップ3105)。これを、コンピュータ断層撮影(CT)または磁気共鳴画像(MRI)技術を使用して行うことができる。しかしながら、この領域におけるCT走査は、重大により速い取得時間(MRIに関する7分に対して30秒)とより高い解像度(MRIに関する1mm立方体ボクセルに対して0.3mm立方体ボクセル)を与えるため、CT走査が好適である。CT走査データを取得するために、GE/CTIスパイラル走査CTスキャナを使用することができる。好適な走査プロトコルは、120keV、200ma、512×512行列サイズ、15cm視野および3mm/2.0:1ピッチである。この走査は、約30秒において計算され、結果として、0.3mm厚の351の画像スライスを生じ、結果として、0.3mm立方体ボクセルを生じる。
【0132】
画像分割を使用し、ボクセルを多数のカテゴリに分類することができる(ステップ3110)。この動作において、変形した自己適合オンラインベクトル量子化(SOVQ)アルゴリズムを使用することができる。このような場合において、前記アルゴリズムは、各々のボクセルを3次までの隣接ボクセルに関して解析し、局所密度特徴を決定する。前記取得されたデータセットにおける各ボクセルは、関連する局所密度ボクセルを有する。前記局所密度ベクトルをカールフネン−ロエベ(K−L)変換を使用して変換することによって、前記体積画像におけるボクセルに関する特徴ベクトルを得ることができる。前記特徴ベクトルに基づいて、前記ボクセルを分類し、ラベル付けする。ボクセル分類は、局所ボクセル密度ベクトルと、最大クラスタ数(MCN)と呼ばれる1つのプリセットパラメータとに部分的に依存する。MCNは、前記データセットに用いられるボクセル分類の数を設定する。CT画像の場合において、人間の目は、4つの識別可能な組織/材料形式を識別することができる。5のMCN値がこの場合において好適である。MRI画像に関し、人間の目は、6つの異なった組織形式の間に差異を認めることができ、7のMCN値を使用することができる。
【0133】
前記画像分割プロセスの一部として、拡張されたデータセットを、前記測定されたデータ点の間の補間によって発生する。例えば、SOVQアルゴリズムの前に、1次ラグランジュ補間を、前記データセットにおける各々のスライスに用いることができる。これは、前記データセットの原スライスの256×256行列サイズを512×512行列サイズに拡張する。加えて、スライス内補間を行い、実際のスライス間の前記データセットをさらに拡張することができる。前記補間されたデータセットを、拡大データセットと呼ぶ。拡大データセットの発生に加え、前記補間プロセスは、前記データにおいてローパスフィルタ効果を有するため、ノイズを抑制し、部分的体積効果を減少させる。
【0134】
2次元ビューイングツールを使用し、種ボクセルを前記喉頭管腔内で選択し、成長アルゴリズムを用いて、前記データセットから喉頭体積を抽出する(ステップ3115)。いくつかの接続されない体積領域が存在するかもしれない前記喉頭のこれらの領域において、多数の種点を選択することができる。
【0135】
前記識別された喉頭体積と、画像分割によって分類された前記領域のボクセルとに関して、次のタスクは、前記仮想喉頭管腔の効率的なナビゲーションおよびビューイングを可能にするように前記データを管理することである。この場合において、詳細レベル(LOD)アプローチを、本方法における使用に適合させ、変更する。このLOD方法において、減少したデータセットを前記拡大データセットから発生する。例えば、512×512×256拡大データセットを、3レベルのしきい値処理を伴う多解像度分解を使用することによって、64×64×32減少体積データセットに減少することができる。次に、前記拡大および減少体積データセットの双方において体積画像をレンダリングするのに使用されるポリゴンを抽出することができる。慣例的なマーチングキューブ法を使用し、前記喉頭の表面に適合するポリゴンを抽出することができる。
【0136】
先行技術において遭遇する1つの問題は、必要なポリゴンの最大数を管理し、前記拡大データセットに関する3次元画像を発生することである。この問題を、本方法において、前記拡大データセットをバイナリ空間分割(BSP)ツリーデータ構造において組織化することによって解決する(ステップ3130)。前記原画像体積を、前記ツリーのルートとして選択する。次に、前記空間を、ほぼ等しい数のポリゴンを含む2つのサブ空間に分割する。この再分プロセスを、各々の結果として生じるサブ空間におけるポリゴン数がしきい値未満になるまで繰り返す。前記しきい値は、システム性能および適用必要条件に基づいて変化することができる。最後の結果として生じるサブ空間を、前記ツリーのリーフノードと呼ぶ。前記再分プロセスが完了したら、前記拡大データセットのボクセルのすべてを、前記BSPツリーのリーフノードに格納する。
【0137】
ナビゲーションまたはビューイング中、ポリゴン選別を、完全に視界の外側のこれらのリーフノードを現在の処理動作から最初に除去することによって適用することができる。残りのポリゴンを前記BSPツリーから呼び戻し、選別されていないこれらの空間において配列し、レンダリングする。したがって、前記BSPツリーは、特定のナビゲーションまたは表示動作に関して前記データセットの関連する部分を選択する有効なツールを提供する。
【0138】
前記拡大および減少データセットを、2レベルLODレンダリングモードにおいて協働して使用する。ユーザが前記物体と、前記視界において回転、シフトまたは他の変更を行うような相互作用をする場合(ステップ3135)、前記減少データセット(64サイズ)からのポリゴンをレンダリングする(ステップ3140)。含まれるポリゴン数が重大により少ないため、前記減少データセット体積との相互作用を、より速く、少ない処理オーバヘッドで行うことができる。増加する速度とのトレードオフは、減少した画像解像度である。予め決められた期間後、ユーザからの相互作用がない場合、拡大データセット(512サイズ)のポリゴンを前記BSPツリーから選択し、レンダリングし、現在の視野の高解像度画像を与える(ステップ3145)。
【0139】
仮想血管造影法
仮想画像化およびナビゲーションに関する前記技術を、仮想血管造影法に適合させ、用いることもできる。この技術を、循環系の種々の異常および疾病の検出および測定に使用することができる。
【0140】
仮想血管造影法のこれらのような用途の1つは、一般に大動脈管の小さい拡大として始まり、大動脈瘤のサイズの増加と共に破裂するより大きな危険性を示す腹部の大動脈瘤の検出である。以前は、唯一の有効な治療方法は切開手術であり、大動脈瘤のレベルにおいて大動脈内に移植用片を配置していた。しかしながらこの手順は、高い程度の関連する不健全性および死亡率を有する。最近開発された、皮膚ごとに配置される大動脈ステントグラフト技術は、著しくより低い複雑割合を持っている。仮想血管造影法は、これらのあまり侵害しない手順を計画するのを助ける有効な方法であり、大動脈瘤を検出し、大動脈瘤の成長を追跡し、手術が必要な場合を決定する有効なツールでもあり得る。
【0141】
図32のフローチャートは、本仮想血管造影法の概要を与える、仮想血管造影法を行うことにおいて、大動脈のような血管の画像走査を取得しなければならない(ステップ3205)。コンピュータ断層撮影(CT)、磁気共鳴画像(MRI)および超音波のような種々の画像化技術を使用することができる。しかしながら、大動脈CT走査が、CT画像において結果として生じる、血と、柔らかい組織と、カルシウム堆積物との間のコントラストのため、一般に好適である。
【0142】
画像走査データセットを取得したら、画像分割技術を前記データセットに用い、前記データセットのボクセルを多数のカテゴリに分類する(ステップ3210)。図15と組み合わせたような上述した画像分割技術は、一般的に利用可能である。この場合において、前記ボクセルの種々の特徴ベクトル値を、血、柔らかい組織およびカルシウム堆積物のようなカテゴリにしたがってグループ化する。血液ボクセルを種として使用し、領域成長アルゴリズムを使用して、前記大動脈管腔の体積および大きさを決定することができる。
【0143】
前記CT画像において、大動脈瘤は、隣接する柔らかい組織によく似た画像特徴を有する。結果として、前記大動脈瘤の完全な輪郭を確立することは困難であるかもしれない。しかしながら、カルシウム堆積物を有する領域は、前記CT走査において重大なコントラストを提供し、これらを使用し、前記血管壁における大動脈瘤の終了点のような、前記大動脈瘤の部分を識別することができる(ステップ3215)。
【0144】
大動脈瘤の一部を検出した後、1つ以上の閉表面を発生し、前記大動脈瘤の輪郭の推定を規定することができる(ステップ3220)。凸面閉表面を、非一様、非論理的なBスプラインを使用し、識別された前記大動脈瘤の点を通るか、これらの近くを通る表面を発生することによって確立することができる。
【0145】
前記閉表面を発生した後、前記大動脈瘤の体積を推定することができる(ステップ3225)。前記体積を推定する1つの方法は、前記推定された閉表面によって囲まれたボクセルの数を計数することである。加えて、前記大動脈瘤のボクセル内で、血流の方向に沿った中心線を、距離変換技術を使用することによって決定することができる。次に、連続的な局所座標系をこの中心線に沿って確立し、前記大動脈瘤の直径を決定することができる。仮想ナビゲーションをこの中心線に沿って、前記結腸のような管腔を通るナビゲーションに関して上述したものと一致するように行う。
【0146】
図33A−Cを参照し、上述した仮想血管造影法を使用し、ステントグラフトの発生および配置において補助し、腹部の大動脈瘤をバイパスするのに使用することができる。図33Aは、腎臓部の動脈の下で前記大動脈の二股の上に位置する腹部の大動脈瘤の単純化した図を示す。前記大動脈の特定の解剖と、前記腹部の大動脈瘤のサイズおよび位置とにおける患者間の変化のため、ステントグラフトを使用して大動脈瘤をバイパスするため、前記グラフトを、個々の大動脈セグメントに特に適合するように設計し、構成しなければならない。図33Bに示すように、これは、とりわけ、必要なグラフトの長さと、前記バイパスされる領域の各々の端におけるインタフェースの点における直径と、インタフェースの角度とを知ることを必要とする。図33Bおよび33Cに示すように、前記大動脈瘤が動脈の枝の近くに位置する場合、二股に分かれたステントグラフトの二股に分かれた端のサイズおよび角度も決定しなければならない。
【0147】
現在まで、これらのような測定は、鼠蹊領域のレベルにおいて行われる挿入により前記大動脈中に挿入されるカテーテルを使用する侵入的な較正血管造影像この技術を、本仮想血管造影法を使用して補足し、ひょっとすると取って代わることができ、この仮想血管造影法は、前記大動脈管腔のブランチを通って構成された中心線を使用する仮想ナビゲーションを使用してこれらのような距離および角度を決定することができる。加えて、前記仮想血管造影法を使用し、ステントグラフトを挿入することができる領域の仮想試験切除を行うことができる。これは、手術者が、動脈表面の下を見ることを可能にし、血栓、堆積物、石灰化像、または、ステントグラフト手順の使用に禁忌を示す他の要因に関する領域を検査することを可能にする。
【0148】
仮想血管造影法の他の用途は、脳に血流を供給する頚動脈を通る画像化、調査およびナビゲーションである。仮想内視鏡検査に関して上述した技術は、血管の場合において完全に適用可能である。例えば、興味ある血管を、取得した画像データから画像分割技術を使用して抽出する。次に、ナビゲーション飛行経路を、前記血管を通って確立することができる。好適には、ポテンシャル場を、ナビゲーションにおける使用に関して前記血管内に構成する。結腸のような他の器官に関するように、前記血管の体積レンダリングされたモデルを発生することができる。前記飛行経路およびポテンシャル場を使用して、前記体積レンダリングされた血管管腔の内部を通ってナビゲーションすることによって、血管狭窄およびプラーク形成のような異常を観察することができる。加えて、仮想試験切除に関して考察した技術を、この文脈において用い、血管壁を評価し、プラークのような血管壁表面における構造を特徴付けることができる。
【0149】
仮想内視鏡検査に関するツリーブランチ探索
中空の器官または物体を通る仮想ナビゲーションに関する経路探索は、重要な仕事である。この目的を達成するために、飛行経路発生のような種々の技術が考察されてきた。多ブランチ構造を与えるような、調査されている物体のジオメトリがより複雑になるにつれて、前記経路探索の仕事はさらにより複雑になる。主要な管腔の中心線を決定するだけでなく、前記主要な管腔から延びる任意のブランチを識別し、位置付けることが望ましい。複雑なブランチ構造を有する器官の一般の例は、主な気道および肺と、心臓血管系と、結腸膨起の折りたたみの存在のため結腸とを含む。各々の器官または物体は、一般に、該物体に関する経路またはスケルトンを規定する特定の挑戦を示す。しかしながら、このおうなスケルトンを発生する一般化された技術を、図34のフローチャートにおいて説明した。
【0150】
図34を参照し、コンピュータ断層撮影(CT)または磁気共鳴画像(MRI)走査のような、興味ある領域の画像走査を取得する(ステップ3405)。上記で考察したように、前記画像走査を、前記画像走査のバイナリ画像を重ね、これらの重ねられた画像から3次元体積ユニットまたはボクセルを規定することによって、3次元体積に変換する(ステップ3410)。前記興味ある領域の体積および複雑さに応じて、前記スケルトンを発生する前に、前記3次元体積のデータセットのサイズに減少することが望ましい。この目的のため、以下により詳細に考察する多解像度データ減少プロセスを使用することができる(ステップ3415)。
【0151】
前記スケルトンは、前記3次元体積のサブセットである。好適には、以下の属性、すなわち、1)前記ツリーのホモトピーを保存する、2)26接続されている、3)1ボクセル厚である、4)前記ブランチの中央軸を推定する、5)比較的滑らかである、を有する。ホモトピーの程度は、ある程度用途特定である。例えば、結腸管腔のスケルトンの発生において、前記スケルトンは、一般に、数ボクセル深さであるかもしれない多数の結腸膨起の折りたたみの存在にもかかわらず、端から端まで単一経路である。しかしながら、心臓血管系および肺系において、数ボクセル深さであるルートからの小さい分枝が、前記系における正当なブランチを表すことができる。
【0152】
図34に戻って、興味ある体積において、ルートボクセルを前記体積において識別する(ステップ3420)。仮想内視鏡検査を行うことにおいて、これを、評価している構造のジオメトリの理解に基づいて手動で行うことができる。
【0153】
次に、距離マップを発生し、前記ツリーにおけるブランチと、前記ブランチおよびルートボクセルの終了点間の距離とを識別する(ステップ3425)。この方法において用いられる仮定は、前記ツリーのルートまで最も長い距離を示す各ブランチにおいて1つのユニークな終了点が存在することである。図35は、26接続ボクセルと呼ばれる3×3×3立方体ボクセル配置を説明する図式的な図である。この配置の中心において、ゼロの距離重みが割り当てられた種ボクセル3505がある。種ボクセル3505の周囲に、26個の接続された隣接ボクセルがあり、これらの隣接ボクセルには、各々の隣接ボクセルと種ボクセルとの間の各々のユークリッド距離に基づいた距離重みが割り当てられる。立方体配置において、ユークリッド距離は、1、1.4および1.7にほぼ等しい
【外11】
の正規化値をとることができる。簡単に処理するために、前記ボクセルに10、14および17の整数値重みを割り当て、関連するユークリッド距離を近似することができる。
【0154】
図36は、図35の26接続立方体距離プレートのユークリッド重み付き距離を使用して前記体積におけるボクセルからルートまでの距離マップを決定するアルゴリズムの擬似コード表現である。発生された距離マップから、ブランチを識別し、前記ブランチの終了点を決定する(ステップ3430)。
【0155】
図36を参照し、前記体積のルートに整数値0をラベル付けする。次に、前記体積におけるボクセルに関して処理待ち行列を形成する。次に、前記ボクセルに、先入れ先出し方法において、前記待ち行列の先頭におけるボクセルとルートボクセルとの間のユークリッド距離を加算することによってラベル付けする。このプロセスを、前記体積におけるボクセルのすべてが前記距離マップにおける値を割り当てられるまで繰り返す。
【0156】
前記距離マップにおけるボクセルのラベル付けは、前記待ち行列の順序に部分的に依存するため、結果として生じる距離マップは、ユニークな解を与えない。しかしながら、前記待ち行列順序にかかわらず、各部ランチに関して少なくとも1つの最も遠い点が常に存在する。加えて、各々のボクセルに関して、ルートボクセルのほかに、前記体積において前記ルートに対してより短い距離を有する少なくとも26の接続された隣接ボクセルが存在する。したがって、前記ブランチの終了点は、前記距離マップを局所最大距離(局所最大値)に関して探索することによって、容易に検出可能である(ステップ3430)。この局所最大値という用語は、相対的な用語である。前記体積を局所最大値に関して評価することにおいて、前記体積を、評価されている物体に適切な種々のサブ空間に分割すべきである。期待される特徴サイズ、ブランチ長、ブランチ直径等を、前記サブ空間分割を決定することにおいて一般に考慮する。
【0157】
前記ブランチの終了点を決定したら、前記終了点からルートボクセルまでの最短経路を決定する(ステップ3435)。前記終了点からルートボクセルまでの最短経路は、前記ツリーのブランチの基本的構造を規定し、前記ブランチの中心線を近似する。これを、前記体積のラフスケルトンと呼ぶ。前記最短経路を、好適には、前記ツリーの最も遠い端におけるブランチから、これらのブランチの端から始めて発生する。最も遠いブランチ終了点から、第1ボクセルを選択し、その26の接続された隣接ボクセルを解析し、どのボクセルが前記終了点からルートまでの最短距離経路においてあるかを決定する。このプロセスを、選択されたボクセルがルートに達するまで繰り返す。これは、結果として、前記最も遠い端からルートまでの1ボクセル幅経路を生じる。他のブランチに関する最短経路の探索も同様である。しかしながら、その後のブランチに関して、前記選択プロセスは、現在経路が以前に割り当てられた経路に達する場合、終了することができる(例えば、前記経路は、前記ルートへのすべての道を導く必要はない)。前記相互接続された最短経路のすべての集合は、前記物体の結果として生じるラフスケルトンである。
【0158】
前記結果として生じるラフスケルトンの用途に応じて、前記ラフスケルトンを改良することが望ましいかもしれない。前記スケルトンを改良するあるステップは、前記スケルトンを前記ブランチ内で中心に集めることである。最も長いブランチから開始して、一様な間隔を、一般に4−8ボクセルの範囲において、前記ブランチに沿って選択する。各々の間隔に関して、前記ラフスケルトンにおけるボクセルの接線方向を計算し、前記接線方向に対して垂直なボクセルと交差する面を決定する。前記面およびボクセルの交線によって規定される2次元領域を形成し、この交線の中心を既知の最大ディスク技術を使用して計算することができる。次に、前記交線の中心をバイキュービックスプラインB補間または他の曲線フィッティング方法を使用して接続することができる。残りのブランチに関して、他のブランチまたはルートに達する終了点を最初に調節し、以前に中心化されたスケルトンブランチの位置を調和させる。次に、中心化を、前記最も長いブランチに関して説明したのと同じように続ける。
【0159】
ステップ3415に戻って参照し、大きいデータセットが含まれる場合、前記データセットのサイズを減少して、処理をスピードアップすると共に処理コストを減少することが必要であるか、少なくとも望ましいかもしれない。前記ツリー構造がある範囲のスケール内で保存され得ることに注意することによって、大きい体積を、構造解析のためのより小さいスケール空間に縮小することができる。
【0160】
多解像度解析に基づく縮小方法を使用することができる。入力データは、CTまたはMRI走査の分割結果から得ることができる同じサイズのバイナリ画像のスタックである。x方向をスライス画像幅に沿ってとり、y方向はスライス画像高さに沿い、z方向はステップずつの方向に沿う。前記ツリー体積における前景ボクセルを128(最大)の値に設定し、背景ボクセルを0(最小)の値に設定する。すべての有理係数にダウベッチーズ(Daubechies)の両直交ウェーブレット変換を用いる。この1次元(1D)離散ウェーブレット変換(DWT)を、最初にx方向に沿って1行ずつ用いる。前記DWTの適用から、より低い周波数成分のみが保持され、もたらされる。前記DTMをバイナリ信号に用いることに注意し、前記より低い周波数成分を結果として生じる2種類の非ゼロ係数が存在する。第1は値128であり、この種類の係数は前記体積内部において位置する。第2は128に等しくない値であり、これらの係数は前記体積の境界に位置する。
【0161】
前記第2の種類の係数を、予め決められたしきい値と比較する。その絶対値が、予め設定されたしきい値T1より大きい場合、前記係数の値を128に設定し、そうでなければ0に設定する。これは、結果として、原データセットの半分の行サイズを有するバイナリ画像のスタックを生じる。次に、同様のしきい値処理を、より低い周波数成分に用いる場合、同じDWTを、前記結果として生じたデータセットに、y方向に沿って、1列ずつ用いる。この結果は、再びバイナリ画像のスタックだが、ここでは、原データセットと比較して半分の行および列サイズを有する。最後に、前記DWTを、前記最後の結果に、z方向に沿って用い、より低い周波数成分を保持する。このステップは、第1レベル分解を完了する。
【0162】
前記第1レベル分解の結果として生じるデータセットは、原データセットと比較して3方向すべてにおいて半分のサイズのものである。前記縮小手順がこのレベルにおいて停止する場合、最後のしきい値処理を用いる。非ゼロおよび非128値のこれらの係数を再評価する。この種類の係数の絶対値が予め設定されたしきい値T2より大きい場合、128として再評価し、そうでなければ0として再評価する。更なる縮小が必要な場合、同じしきい値アルゴリズムをしきい値T1によって用いる。上述したような更なる縮小プロセスを、結果として生じる体積が所望の減少したデータ体積に達するまで再帰的に用いることができる。仮想内視鏡検査において、前記スライス画像を512×512画素サイズとする。最大分解レベルは、通常3であり、結果として64×64減少画素サイズを生じる。
【0163】
前記体積を、前記与えた方法によってすべての方向において等方的に縮小する。2つのプリセットしきい値T1およびT2を使用し、縮小の程度を制御する。前記体積が著しく過縮小である場合、前記減少した体積において連結度が失われてしまうかもしれない。より小さい程度に過縮小である場合、前記減少した体積データセットにおいて2つの別個のブランチが1つのブランチに合体するかもしれない。前記2つのしきい値がより大きくなると、より薄く減少した体積になる。これらの2つのしきい値の範囲は、[0,r×128]であり、0<r<1である。好適には、仮想内視鏡検査に関する範囲は、T1に関してr∈(0.08,0.28)、T2に関してr∈(0.7,0.98)である。正確な決定は、特定の用途の特徴サイズに依存し、縮小体積における構造特徴の忠実度を保持しながら、減少を達成するように選択される。
【0164】
原体積の縮小後、前記ツリーブランチ探索手順を、前記より小さい体積に用いることができる(ステップ3420−3440)。結果として生じるスケルトンを、原スケール空間にマッピングし戻すことができる。原空間に対して測った場合、より小さいスケールのスケルトンの画像は、もはや、前記原スケール空間にける接続された経路のままではない。前記画像におけるこれらのボクセルは、最終スケルトンに関する制御点として働く。前記制御点を、前に説明したようなアルゴリズムを使用して中心に集め、次に、これらを補間し、前記オブジェクトに関する最終スケルトンを形成する。
【0165】
コンピュータ支援診断
ここに説明した仮想調査技術は、種々の条件のコンピュータ支援診断(CAD)に関する用途に適する。例えば、上述したように、器官組織のジオメトリをガウス曲線に関して調査することによって、結腸管腔内部のポリープのような異常ジオメトリを有する領域を自動的に識別することができる。この技術を一般化し、テクスチャ特徴と共に使用して、多数の用途に関するCAD機能性を与えることができる。
【0166】
例えば、上述した膀胱の多走査画像化を使用し、膀胱壁における腫瘍の自動検出を行うことができる。この場合において、一般に膀胱壁内における腫瘍侵入の程度を使用し、膀胱癌の段階を規定する。上述した多走査画像化および画像分割技術を使用し、前記膀胱壁の領域の輪郭を容易に描くことができる。通常の膀胱組織の領域は、一般に、実質的に一様なテクスチャ特徴を示す。しかしながら、腫瘍が前記領域において存在する場合、前記一様なテクスチャ特徴は中断される。したがって、テクスチャ解析を使用し、前記膀胱壁を評価することによって、腫瘍を示すかもしれない領域は、前記一様なテクスチャ内の乱れまたは“ノイジー領域”としてそれ自身を示す。
【0167】
領域のテクスチャを、規定された範囲内のボクセル間の強度相関を特徴付ける確率分布関数(PDF)によって表すことができる。2次元PDFを使用し、テクスチャ特徴を表すことができる。このようなPDFは、すべての方向に沿った2つの最も近いボクセル間の相関を特徴付ける。前記PDFを評価するために、興味ある領域における任意の2つの最も近い隣接ボクセルの強度を、前記興味ある領域に関する標本ベクトルとして記録することができる(例えば文脈)。多数のこのような標本ベクトルを使用し、この文脈に関して前記PDFを評価する累積分布関数(CDF)を発生することができる。各ボクセルに関して、その隣接ボクセルのある範囲内の標本ベクトルを使用し、局所CDFを発生することもできる。
【0168】
コルモロゴフ−スミルノフテストのような統計テストを前記CDFに用い、前記文脈のCDFおよび局所CDFが統計的に等価であるか、例えば、予め規定された信頼レベル内であるかどうかを決定することができる。そうであるならば、現在ボクセルの周囲の局所テクスチャ特徴を前記文脈と同一であるとみなす。そうでなければ、現在ボクセルは、前記文脈のテクスチャ特徴と異なったテクスチャ特徴を示し、腫瘍のような潜在的な異常とみなしてもよい。あるボクセルが前記文脈と統計的に等価であるかどうかを決定するのに使用される信頼レベルを変化させ、検出の感度を増減させることができる。
【0169】
テクスチャ解析にPDFまたはCDFを用いる代わりの方法において、各々のCDFまたはPDFを、機能的線形空間における点とみなすことができる。この空間における2つのDCFまたはPDFの間の距離を、スコロホールド(Skorohold)距離等によって測定することができる。この距離は、PDFの類似性の程度の測定を与える。例えば、局所CDFと文脈CDFとの間の距離を計算することができ、結果として生じる距離を1つ以上の距離しきい値と比較することができる。前記距離が大きい場合、前記局所テクスチャを、前記文脈と異なると考え、このことは、このようなボクセルが、潜在的な異常または腫瘍を有する領域に属することを示すことができる。好適には、前記距離しきい値を、統計的に十分に既知のデータセットの評価に基づいて決定する。
【0170】
上記で計算した距離を、ここで説明したような視覚化技術および体積レンダリング技術と共に使用することができる。例えば、原データセットと比較し得るサイズを有する特徴体積データセットを形成することができる。この新たなデータセットにおける各ボクセルに関する強度を、前記局所CDFと文脈のCDFとの間の距離に基づいて割り当てることができる。この3次元体積データセットを体積レンダリング技術によって見る場合、疑わしい腫瘍を含む領域は、周囲の領域よりも高い画像強度を示す。
【0171】
ポリープの自動検出に関して上記で考察したように、管腔の表面を、C−2滑らか表面モデル等を使用することによって、3次元ユークリッド空間における連続的に二次微分可能表面として表すことができる。このようなモデルにおいて、前記結腸の表面における各ボクセルは、ガウス曲率場と呼ばれる、ガウス曲率を有する関連する幾何学的特長を有する。種々の器官に関して、ある予想される局所特徴を、別個の曲率テンプレートによって特徴付けることができる。例えば、前記結腸の文脈において、予想される局所特徴は、滑らかな湾曲表面と、リング折りたたみと、滑らかな面からの凸面の丘と、滑らかな表面からの台地とを含む。これらの少なくとも2つの特徴は、ポリープまたは腫瘍を表すかもしれない。したがって、前記ガウス曲率場を特定の予め決められた局所特徴テンプレートに関して探索することによって、ポリープ、腫瘍、および、興味ある他の異常を自動的に検出することができる。コンピュータ支援診断を行うことへの表面ジオメトリのこの使用を、単独で使用することもでき、または、上述したテクスチャベースCAD技術と共に使用することもできる。
【0172】
器官全体または興味ある領域をレンダリングして見る代わりとして、観察されている表面を小さい領域またはパッチに分割することができ、これらのパッチを、前記局所曲率テンプレートによって規定する。各々のパッチは、共通の幾何学的特長または曲率テンプレートを有するボクセルを含むべきである。次に、前記パッチのすべてのボクセルを観察することを可能にする単一の視点を、前記パッチに関して決定する。次に、前記パッチに、該パッチがポリープまたは他の異常を表す可能性を示す優先度スコアを割り当てる。次に、オペレータが前記内蔵体積全体をナビゲートして疑わしい領域を探し出すのではなく、前記パッチを、優先度順に別々に観察することができる。もちろん、好適な診断システムは、オペレータが、あるパッチを見ることから、前記器官の文脈におけるパッチを見ることへ容易に変更することができるように、ビューの間で切り替える能力を含む。代わりに、これら2つのビューを同時に与えることができる。再び、前記テクスチャベースアプローチを使用し、このアプローチを補うことができる。前記テクスチャ解析の結果を調査している前記パッチにマッピングすることによって、前記テクスチャ情報を観察し、診断に使用することもできる。
【0173】
上記は、本画像化および調査システムおよび方法の原理を単に説明する。したがって、当業者が、ここに明白に図示または説明していないが、本発明の原理を具体化し、したがって、請求項によって規定される本発明の精神および範囲内である多数のシステム、装置および方法を考案することができることは明らかであろう。
【0174】
例えば、ここで説明した方法およびシステムを用い、動物、魚、または無生物を仮想調査することができる。医療分野における指定された使用に加え、前記技術の用途を使用し、開くことができない密封された物体の内容を検出することができる。前記技術を、ビルディングまたは洞窟のような建築上の構造内で使用することもでき、オペレータが前記構造中をナビゲートすることを可能にすることができる。
【図面の簡単な説明】
【図1】本発明による、物体、特に結腸の仮想調査を行うステップのフローチャートである。
【図2】仮想器官におけるガイド付きナビゲーションを行う「サブマリン」カメラモデルの説明である。
【図3】「サブマリン」カメラのピッチおよびロールをモデル化するのに使用される振り子の説明である。
【図4】2つのブロッキング壁を識別する体積測定結腸の二次元断面を説明する図である。
【図5】開始および終了体積要素を選択する体積測定結腸の二次元断面を説明する図である。
【図6】前記ブロッキング壁および結腸表面によって閉じた離散的サブ体積を示す体積測定結腸の二次元断面を説明する図である。
【図7】多数の層を剥ぎ取られた体積測定結腸の二次元断面を説明する図である。
【図8】残りの飛行経路を含む堆積測定結腸の二次元断面を説明する図である。
【図9】走査された結腸の体積視覚化を発生するステップのフローチャートである。
【図10】セルに再分された仮想結腸の説明である。
【図11】Aは仮想調査されている器官のグラフィカル描写であり、BはAにおける器官を描写すると同時に、発生されたスタブツリーのグラフィカル描写であり、CはAにおける器官を描写すると同時に、発生されたスタブツリーの他のグラフィカル描写である。
【図12】Aはレンダリングすべきシーンおよび前記シーンの特定のセル内の物体のグラフィカル描写であり、BはAにおける画像を描写すると同時に、発生されたスタブツリーのグラフィカル描写であり、C−EはAにおける画像を描写すると同時に、発生されたスタブツリーの他のグラフィカル描写である。
【図13】層を除去したポリープを含む仮想結腸の2次元表現である。
【図14】本発明による人間の器官の仮想調査を行うのに使用されるシステムの図である。
【図15】改善された画像分割方法を示すフローチャートである。
【図16】代表的な腹部CTデータセットのボクセル強度対周波数のグラフである。
【図17】興味あるボクセルと、その選択された隣接ボクセルとを含む強度ベクトル構造の斜視図である。
【図18】Aは肺を含む領域を第1に説明する人間の腹部領域のCT走査からの好例の画像スライスであり、BはAの画像スライスにおける肺領域の識別を説明する絵的図であり、CはBにおいて識別された肺体積の除去を説明する絵的図である。
【図19】Aは結腸および骨の一部を含む領域を第1に説明する人間の腹部領域のCT走査からの好例の画像スライスであり、BはAの画像スライスからの結腸および骨領域の識別を説明する絵的図であり、Cは骨の領域を除去したAの画像走査を説明する絵的図である。
【図20】テクスチャをモノクローム画像データに用いる方法を説明するフローチャートである。
【図21】高速パースペクティブ光線投影技術を用いて体積レンダリングする方法を説明するフローチャートである。
【図22】体積縮小技術を用いて結腸管腔を通る中心飛行経路を決定する方法を説明するフローチャートである。
【図23】図22に示す方法において使用する体積縮小技術をさらに説明するフローチャートである。
【図24】発生された中心飛行経路による分割された結腸管腔の3次元絵的表現である。
【図25】分割技術を用いて結腸管腔を通る中心飛行経路を発生する方法を説明するフローチャートである。
【図26】パーソナルコンピュータバスアーキテクチャに基づくシステム実施形態のブロック図である。
【図27】図26のシステムを使用して体積画像化を行う方法を説明するフローチャートである。
【図28】膀胱のような物体の仮想調査(仮想膀胱鏡検査)を行う多走査方法を説明するフローチャートである。
【図29】図28の仮想膀胱鏡検査による画像化結果を与え、膀胱構造の説明的外側ビューを与えるのに好適な表示ウィンドウの絵的表現である。
【図30】図28の仮想膀胱鏡検査による画像化結果を与え、膀胱構造の説明的内側ビューを与えるのに好適な表示ウィンドウの絵的表現である。
【図31】多解像度ビューイングを使用して喉頭のような物体の仮想調査を行う方法のフローチャートである。
【図32】仮想血管造影検査を行う方法のフローチャートである。
【図33】A−Cは腹部大動脈瘤の存在を説明する大動脈の一部の絵的ビューである。
【図34】物体のスケルトン構造を発生する方法を説明するフローチャートである。
【図35】26接続ユークリッド重み付き立方体距離プレートの図式的な図である。
【図36】図34の方法において使用する距離マップを発生するプロセスの擬似コードを説明する図である。
Claims (61)
- 物体の仮想調査を行う方法において、
コントラスト剤の存在によって膨張させた物体の少なくとも1回の走査を行うステップと、
前記コントラスト剤を除去した前記物体の少なくとも1回の走査を行うステップと、
前記走査を、複数のボクセルを具える対応する体積データセットに変換するステップと、
画像分割を行い、各々の走査のボクセルを複数のカテゴリに分類するステップと、
各々の走査の体積データセットを共通座標系に登録するステップと、
前記体積データセットの少なくとも2つを本質的に同時に表示するステップと、
前記体積データセットの一方において仮想ナビゲーション動作を行い、対応するナビゲーション動作を少なくとも1つの他の体積データセットにおいて生じさせるステップとを含むことを特徴とする、仮想調査を行う方法。 - 請求項1に記載の仮想調査を行う方法において、前記膨張させた物体の少なくとも1回の走査が、前記物体の横走査および冠状走査を含むことを特徴とする、仮想調査を行う方法。
- 請求項2に記載の仮想調査を行う方法において、前記コントラスト剤を除去した前記物体の少なくとも1回の走査が、前記物体の横走査および冠状走査を含むことを特徴とする、仮想調査を行う方法。
- 請求項3に記載の仮想調査を行う方法において、前記物体を膀胱としたことを特徴とする、仮想調査を行う方法。
- 請求項4に記載の仮想調査を行う方法において、前記走査をコンピュータ断層撮影走査としたことを特徴とする、仮想調査を行う方法。
- 請求項4に記載の仮想調査を行う方法において、前記走査を超音波画像走査としたことを特徴とする、仮想調査を行う方法。
- 請求項4に記載の仮想調査を行う方法において、前記走査を磁気共鳴画像走査としたことを特徴とする、仮想調査を行う方法。
- 請求項7に記載の仮想調査を行う方法において、前記コントラスト剤を尿としたことを特徴とする、仮想調査を行う方法。
- 請求項1に記載の仮想調査を行う方法において、前記コントラスト剤を除去した前記物体の少なくとも1回の走査が、前記物体の横走査および冠状走査を含むことを特徴とする、仮想調査を行う方法。
- 請求項1に記載の仮想調査を行う方法において、前記物体を膀胱としたことを特徴とする、仮想調査を行う方法。
- 請求項10に記載の仮想調査を行う方法において、前記走査をコンピュータ断層撮影走査としたことを特徴とする、仮想調査を行う方法。
- 請求項10に記載の仮想調査を行う方法において、前記走査を超音波画像走査としたことを特徴とする、仮想調査を行う方法。
- 請求項10に記載の仮想調査を行う方法において、前記走査を磁気共鳴画像走査としたことを特徴とする、仮想調査を行う方法。
- 請求項13に記載の仮想調査を行う方法において、前記コントラスト剤を尿としたことを特徴とする、仮想調査を行う方法。
- 請求項1に記載の仮想調査を行う方法において、前記膨張させた物体の少なくとも1回の走査と、前記コントラスト剤を除去した前記物体の少なくとも1回の走査とを評価し、前記走査の一方においてコントラストがより可視である領域を識別するステップと、興味ある領域においてより高いコントラストを有する前記走査を評価し、前記物体の生理学的特長を決定するステップとをさらに含むことを特徴とする、仮想調査を行う方法。
- 請求項15に記載の仮想調査を行う方法において、前記画像分割のステップが、ボクセルを前記ボクセルの局所強度ベクトルに基づいて分類するステップを含むことを特徴とする、仮想調査を行う方法。
- 請求項16に記載の仮想調査を行う方法において、前記画像分割のステップが、領域成長アルゴリズムを使用し、前記分類されたボクセルに基づいて前記物体の領域を識別するステップをさらに含むことを特徴とする、仮想調査を行う方法。
- 請求項1に記載の仮想調査を行う方法において、前記体積画像データセットを、前記座標系に関する複数の領域に分割するステップをさらに含むことを特徴とする、仮想調査を行う方法。
- 請求項18に記載の仮想調査を行う方法において、前記複数の領域が、3次元座標系において規定された8個の領域を含むことを特徴とする、仮想調査を行う方法。
- 物体の仮想調査を行う方法において、
前記物体の画像走査を行い、画像走査データを取得するステップと、
前記取得された画像走査データを複数のボクセルに変換するステップと、
前記ボクセル間を補間し、拡張データセットを発生するステップと、
画像分割を行い、前記ボクセルを複数のカテゴリに分類するステップと、
前記拡張データセットから前記物体内部の体積を抽出するステップと、
前記拡張データセットから減少解像度データセットを発生するステップと、
前記拡張データセットをツリーデータ構造に格納するステップと、
前記拡張データセットおよび減少解像度データセットに関する画像をレンダリングするステップと、
前記減少解像度データセットおよび拡張データセットのうち少なくとも1つを表示に関して選択するステップとを含むことを特徴とする、仮想調査を行う方法。 - 請求項20に記載の仮想調査を行う方法において、前記選択するステップが、
画像相互作用中前記減少解像度データセットを選択するステップと、
画像相互作用が予め決定された期間中に起こらなかった場合、前記拡張データセットを選択するステップとを含むことを特徴とする、仮想調査を行う方法。 - 請求項20に記載の仮想調査を行う方法において、前記画像走査をコンピュータ断層撮影走査としたことを特徴とする、仮想調査を行う方法。
- 請求項20に記載の仮想調査を行う方法において、前記画像走査を磁気共鳴画像走査としたことを特徴とする、仮想調査を行う方法。
- 請求項20に記載の仮想調査を行う方法において、前記画像走査を超音波画像走査としたことを特徴とする、仮想調査を行う方法。
- 請求項20に記載の仮想調査を行う方法において、前記物体を喉頭としたことを特徴とする、仮想調査を行う方法。
- 請求項20に記載の仮想調査を行う方法において、前記ツリー構造をバイナリ空間分割ツリー構造としたことを特徴とする、仮想調査を行う方法。
- 仮想血管造影検査を行う方法において、
大動脈の少なくとも一部を含む画像走査データを取得するステップと、
前記画像走査データを、複数のボクセルを含む体積表現に変換するステップと、
前記体積表現を分割し、前記ボクセルを複数のカテゴリのうち1つに分類するステップと、
前記分割された体積表現を解析し、大動脈壁における動脈瘤の少なくとも一部を表すボクセルを識別するステップと、
動脈瘤の少なくとも一部を表すボクセルの周囲の少なくとも1つの閉表面を発生し、前記動脈瘤の輪郭を評価するステップとを含むことを特徴とする、仮想血管造影検査を行う方法。 - 請求項27に記載の仮想血管造影検査を行う方法において、前記画像走査をコンピュータ断層撮影走査としたことを特徴とする、仮想血管造影検査を行う方法。
- 請求項27に記載の仮想血管造影検査を行う方法において、前記画像走査を磁気共鳴画像走査としたことを特徴とする、仮想血管造影検査を行う方法。
- 請求項27に記載の仮想血管造影検査を行う方法において、前記分割動作が、ボクセルを、少なくとも、血、組織およびカルシウム堆積物のカテゴリに分類することを特徴とする、仮想血管造影検査を行う方法。
- 請求項27に記載の仮想血管造影検査を行う方法において、前記動脈瘤の体積を、前記発生された閉表面を使用して評価するステップをさらに含むことを特徴とする、仮想血管造影検査を行う方法。
- 請求項27に記載の仮想血管造影検査を行う方法において、前記大動脈管腔を通るナビゲーション経路を発生するステップをさらに含むことを特徴とする、仮想血管造影検査を行う方法。
- 請求項32に記載の仮想血管造影検査を行う方法において、前記動脈瘤の長さを、前記ナビゲーション経路に基づいて評価するステップをさらに含むことを特徴とする、仮想血管造影検査を行う方法。
- 血管の仮想内視鏡検査を行う方法において、
前記血管の少なくとも一部を含む画像走査データを取得するステップと、
前記画像走査データを、複数のボクセルを含む体積表現に変換するステップと、
前記体積表現を分割し、前記ボクセルを、血、組織およびカルシウム堆積物のカテゴリを含む複数のカテゴリのうちの1つに分類するステップと、
前記血管を通るナビゲーション経路を発生するステップとを含むことを特徴とする、血管の仮想内視鏡検査を行う方法。 - 請求項34に記載の血管の仮想内視鏡検査を行う方法において、前記血管を頚動脈としたことを特徴とする、血管の仮想内視鏡検査を行う方法。
- 請求項34に記載の血管の仮想内視鏡検査を行う方法において、前記ナビゲーション経路に沿って前記頚動脈の直径を決定し、狭窄の領域を識別するステップをさらに含むことを特徴とする、血管の仮想内視鏡検査を行う方法。
- 請求項34に記載の血管の仮想内視鏡検査を行う方法において、前記画像走査をコンピュータ断層撮影走査としたことを特徴とする、血管の仮想内視鏡検査を行う方法。
- 請求項34に記載の血管の仮想内視鏡検査を行う方法において、前記画像走査を磁気共鳴画像走査としたことを特徴とする、血管の仮想内視鏡検査を行う方法。
- 仮想血管造影検査を使用してステントグラフトの特性を決定する方法において、
大動脈の少なくとも一部を含む画像走査データを取得するステップと、
前記画像走査データを、複数のボクセルを含む体積表現に変換するステップと、
前記体積表現を分割し、前記ボクセルを複数のカテゴリのうちの1つに分類するステップと、
前記分割された体積表現を解析し、前記大動脈壁における動脈瘤の少なくとも一部を表すボクセルを識別するステップと、
動脈瘤の少なくとも一部を表すボクセルの周囲の少なくとも1つの閉表面を発生し、前記動脈瘤の輪郭を評価するステップと、
前記動脈瘤輪郭の終了点の位置を識別するステップと、
前記動脈瘤輪郭の終了点間の長さを計算し、前記ステントグラフトの長さを決定するステップと、
前記動脈瘤輪郭の終了点における前記大動脈管腔の直径を計算し、前記ステントグラフトの必要な外側直径を決定するステップとを含むことを特徴とする、仮想血管造影検査を使用してステントグラフトの特性を決定する方法。 - 請求項39に記載の仮想血管造影検査を使用してステントグラフトの特性を決定する方法において、前記動脈瘤と通常の大動脈管腔との界面の角度を決定し、前記ステントグラフトの対応する端の角度方向を決定するステップをさらに含むことを特徴とする、仮想血管造影検査を使用してステントグラフトの特性を決定する方法。
- 請求項39に記載の仮想血管造影検査を使用してステントグラフトの特性を決定する方法において、前記動脈瘤に最も近い動脈ブランチの位置を突き止め、前記ステントグラフトの最大長を決定するステップをさらに含むことを特徴とする、仮想血管造影検査を使用してステントグラフトの特性を決定する方法。
- 請求項41に記載の仮想血管造影検査を使用してステントグラフトの特性を決定する方法において、前記動脈瘤に最も近い動脈ブランチが、腎臓および大腿部動脈ブランチの少なくとも一方を含むことを特徴とする、仮想血管造影検査を使用してステントグラフトの特性を決定する方法。
- 請求項39に記載の仮想血管造影検査を使用してステントグラフトの特性を決定する方法において、前記動脈瘤の端に最も近い大動脈領域の仮想生体組織検査を行い、予測されたグラフトインタフェース位置における前記組織の特性を決定するステップをさらに含むことを特徴とする、仮想血管造影検査を使用してステントグラフトの特性を決定する方法。
- 複数のボクセルで形成された中空物体の3次元画像表現に関するスケルトンを規定する方法において、
前記中空物体内におけるルートボクセルを識別するステップと、
前記中空物体内のすべてのボクセルに関する距離マップを発生するステップであって、前記距離マップが、ユークリッド重み付き距離を有する隣接ボクセルの26接続立方体プレートを使用して形成される、ステップと、
前記距離マップにおける局所最大値を有するボクセルを、前記中空物体におけるブランチの終了点として識別するステップと、
各々の局所最大値ボクセルに関して、前記ルートボクセルに対する最短接続経路、または、前もって規定された最短経路を決定するステップとを含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。 - 請求項44に記載の3次元画像表現に関するスケルトンを規定する方法において、多解像度データ減少を前記3次元画像表現に行い、前記発生および識別動作に関する減少データセットを発生するステップをさらに含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。
- 請求項44に記載の3次元画像表現に関するスケルトンを規定する方法において、前記最短経路を、前記物体の個々のブランチ内の中心に配置するステップをさらに含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。
- 請求項44に記載の3次元画像表現に関するスケルトンを規定する方法において、前記物体が少なくとも1つの血管を含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。
- 請求項44に記載の3次元画像表現に関するスケルトンを規定する方法において、前記物体が肺の気道を含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。
- 請求項44に記載の3次元画像表現に関するスケルトンを規定する方法において、前記物体が膀胱を含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。
- 請求項44に記載の3次元画像表現に関するスケルトンを規定する方法において、前記物体が脊椎動物の脊髄を含むことを特徴とする、3次元画像表現に関するスケルトンを規定する方法。
- 興味ある領域のコンピュータ支援診断を行う方法において、
前記興味ある領域の少なくとも一部を含む画像走査データを取得するステップと、
前記画像走査データを、複数のボクセルを含む体積表現に変換するステップであって、前記ボクセルの少なくとも一部が前記興味ある領域の表面を表す、ステップと、
表面を表すボクセルの前記部分を、異常を示す幾何学的特長およびテクスチャ特徴の少なくとも一方に関して解析するステップとを含むことを特徴とする、コンピュータ支援診断を行う方法。 - 請求項51に記載のコンピュータ支援診断を行う方法において、前記テクスチャ特徴を、ボクセルの前記部分の2つのボクセル間の相関を特徴付ける確率密度関数に含めることを特徴とする、コンピュータ支援診断を行う方法。
- 請求項52に記載のコンピュータ支援診断を行う方法において、前記2つのボクセルを隣接するボクセルとしたことを特徴とする、コンピュータ支援診断を行う方法。
- 請求項52に記載のコンピュータ支援診断を行う方法において、ボクセルの前記部分の強度を使用し、前記確率密度関数の評価を発生することを特徴とする、コンピュータ支援診断を行う方法。
- 請求項54に記載のコンピュータ支援診断を行う方法において、複数のボクセル強度を使用し、前記興味ある領域の累積分布関数と、局所累積分布関数とを発生し、前記局所累積分布関数を文脈累積分布関数と比較し、異常な領域を識別することを特徴とする、コンピュータ支援診断を行う方法。
- 請求項55に記載のコンピュータ支援診断を行う方法において、前記局所累積分布関数と文脈累積分布関数との間の距離を決定し、前記距離が異常の基準を与えることを特徴とする、コンピュータ支援診断を行う方法。
- 請求項56に記載のコンピュータ支援診断を行う方法において、前記距離を使用し、強度値を、前記興味ある領域の表面を表すボクセルに割り当て、該方法が、強度における変化が異常な領域を表すように前記ボクセルを表示するステップをさらに含むことを特徴とする、コンピュータ支援診断を行う方法。
- 請求項57に記載のコンピュータ支援診断を行う方法において、前記興味ある領域が結腸を含み、前記異常がポリープを含むことを特徴とする、コンピュータ支援診断を行う方法。
- 請求項51に記載のコンピュータ支援診断を行う方法において、前記興味ある領域が大動脈を含み、前記異常が腹部の大動脈瘤を含むことを特徴とする、コンピュータ支援診断を行う方法。
- 請求項51に記載のコンピュータ支援診断を行う方法において、前記表面を二次微分可能な表面として表し、各々の表面体積ユニットが関連するガウス曲率を有し、前記ガウス曲率が結合して、前記幾何学的特長を形成することを特徴とする、コンピュータ支援診断を行う方法。
- 請求項59に記載のコンピュータ支援診断を行う方法において、複数の予め決定された幾何学的特長テンプレートを規定し、前記表面の幾何学的特長を前記テンプレートと比較し、幾何学的特長分類を決定することを特徴とする、コンピュータ支援診断を行う方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US23766500P | 2000-10-03 | 2000-10-03 | |
US09/777,120 US7194117B2 (en) | 1999-06-29 | 2001-02-05 | System and method for performing a three-dimensional virtual examination of objects, such as internal organs |
PCT/US2001/030704 WO2002029764A1 (en) | 2000-10-03 | 2001-10-01 | System and method for performing a three-dimensional virtual examination of objects, such as internal organs |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004510515A true JP2004510515A (ja) | 2004-04-08 |
Family
ID=26930898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002533259A Pending JP2004510515A (ja) | 2000-10-03 | 2001-10-01 | 器官のような物体の3次元仮想検査を行うシステムおよび方法 |
Country Status (5)
Country | Link |
---|---|
US (2) | US7194117B2 (ja) |
EP (1) | EP1342225A4 (ja) |
JP (1) | JP2004510515A (ja) |
AU (1) | AU2002211331A1 (ja) |
WO (1) | WO2002029764A1 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006192254A (ja) * | 2005-01-12 | 2006-07-27 | Medison Co Ltd | 超音波診断システム及び方法 |
JP2007229464A (ja) * | 2006-02-28 | 2007-09-13 | Siemens Ag | 身体構成要素内の物質の濃度の決定方法および装置 |
JP2008503302A (ja) * | 2004-06-23 | 2008-02-07 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 3d管状オブジェクトのパラメータに関する情報を表示するための画像処理システム |
JP2010504794A (ja) * | 2006-09-29 | 2010-02-18 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 突起検出方法、システム及びコンピュータプログラム |
JP2012050606A (ja) * | 2010-08-31 | 2012-03-15 | Fujifilm Corp | 内視鏡画像処理装置および方法、並びに、プログラム |
JP2021002402A (ja) * | 2020-10-06 | 2021-01-07 | 株式会社ソニー・インタラクティブエンタテインメント | 情報処理装置 |
Families Citing this family (306)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
US6814752B1 (en) * | 2000-03-03 | 2004-11-09 | Endovascular Technologies, Inc. | Modular grafting system and method |
JP3254451B2 (ja) * | 2000-03-06 | 2002-02-04 | 経済産業省産業技術総合研究所長 | 多チャンネルmri画像処理によるカラー化方法及び装置 |
DE10100572A1 (de) * | 2001-01-09 | 2002-07-11 | Philips Corp Intellectual Pty | Verfahren zur Darstellung des Blutflusses in einem Gefäßbaum |
JP2002297678A (ja) * | 2001-03-30 | 2002-10-11 | Fujitsu Nagano Systems Engineering Ltd | Voxelモデルの座標決定処理方法,座標決定処理プログラムおよび座標決定処理プログラム記録媒体 |
WO2002093496A1 (en) * | 2001-05-15 | 2002-11-21 | Koninklijke Philips Electronics N.V. | Analysis of a multi-dimensional data set |
US6826297B2 (en) * | 2001-05-18 | 2004-11-30 | Terarecon, Inc. | Displaying three-dimensional medical images |
EP1411824A2 (en) * | 2001-07-16 | 2004-04-28 | Art Advanced Research Technologies Inc. | Multi-wavelength imaging of highly turbid media |
WO2003021532A2 (en) * | 2001-09-06 | 2003-03-13 | Koninklijke Philips Electronics N.V. | Method and apparatus for segmentation of an object |
ATE481694T1 (de) * | 2001-10-16 | 2010-10-15 | Koninkl Philips Electronics Nv | Verfahren zur auswahl von ästen für das ausrichten von sonden |
JP4248399B2 (ja) * | 2001-10-16 | 2009-04-02 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 自動枝ラベリング方法 |
EP1436771B1 (en) * | 2001-10-16 | 2011-06-22 | The University of Chicago | Computer-aided detection of three-dimensional lesions |
EP1306805A1 (en) * | 2001-10-25 | 2003-05-02 | Mitsubishi Electric Information Technology Centre Europe B.V. | Image Analysis |
US20110206247A1 (en) * | 2001-11-21 | 2011-08-25 | Dachille Frank C | Imaging system and methods for cardiac analysis |
US7359538B2 (en) | 2001-11-23 | 2008-04-15 | R2 Technology | Detection and analysis of lesions in contact with a structural boundary |
US6711231B2 (en) | 2001-12-10 | 2004-03-23 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus to assist and facilitate vessel analysis |
US7570802B2 (en) * | 2001-12-27 | 2009-08-04 | The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services | Automated centerline detection algorithm for colon-like 3D surfaces |
NO20020780D0 (no) * | 2002-02-15 | 2002-02-15 | Amersham Health As | Bildeanalyse |
JP2003296383A (ja) * | 2002-04-05 | 2003-10-17 | Denso Corp | 3次元モデリングシステム |
US20030210262A1 (en) * | 2002-05-10 | 2003-11-13 | Tripath Imaging, Inc. | Video microscopy system and multi-view virtual slide viewer capable of simultaneously acquiring and displaying various digital views of an area of interest located on a microscopic slide |
US7819806B2 (en) | 2002-06-07 | 2010-10-26 | Verathon Inc. | System and method to identify and measure organ wall boundaries |
US8221322B2 (en) | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods to improve clarity in ultrasound images |
US8221321B2 (en) | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods for quantification and classification of fluids in human cavities in ultrasound images |
GB2391625A (en) | 2002-08-09 | 2004-02-11 | Diagnostic Ultrasound Europ B | Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams |
DE10231061A1 (de) * | 2002-07-10 | 2004-01-22 | Philips Intellectual Property & Standards Gmbh | Verfahren und System zur Verbesserung des Informationsgehaltes in einem Bild |
US6905468B2 (en) * | 2002-09-18 | 2005-06-14 | Diagnostic Ultrasound Corporation | Three-dimensional system for abdominal aortic aneurysm evaluation |
US7873223B2 (en) * | 2002-10-15 | 2011-01-18 | Definiens Ag | Cognition integrator and language |
US7801361B2 (en) * | 2002-10-15 | 2010-09-21 | Definiens Ag | Analyzing pixel data using image, thematic and object layers of a computer-implemented network structure |
US8594410B2 (en) | 2006-08-28 | 2013-11-26 | Definiens Ag | Context driven image mining to generate image-based biomarkers |
JP4346297B2 (ja) * | 2002-10-22 | 2009-10-21 | 株式会社東芝 | X線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 |
JP2004141514A (ja) * | 2002-10-28 | 2004-05-20 | Toshiba Corp | 画像処理装置及び超音波診断装置 |
US8521260B2 (en) * | 2002-12-02 | 2013-08-27 | Yeda Research And Development Co. Ltd. | Characterization of arteriosclerosis by optical imaging |
US7551758B2 (en) * | 2002-12-04 | 2009-06-23 | Koninklijke Philips Electronics N.V. | Medical viewing system and method for detecting borders of an object of interest in noisy images |
JP4409166B2 (ja) * | 2002-12-05 | 2010-02-03 | オリンパス株式会社 | 画像処理装置 |
US7081088B2 (en) * | 2003-01-30 | 2006-07-25 | Siemens Corporate Research, Inc. | Method and apparatus for automatic local path planning for virtual colonoscopy |
ATE357706T1 (de) * | 2003-02-18 | 2007-04-15 | Koninkl Philips Electronics Nv | Bildsegmentierung durch zuweisen von klassen zu adaptivenmesh-primitiven |
US8570323B2 (en) * | 2003-02-18 | 2013-10-29 | Koninklijke Philips N.V. | Volume visualization using tissue mix |
US7304644B2 (en) * | 2003-03-12 | 2007-12-04 | Siemens Medical Solutions Usa, Inc. | System and method for performing a virtual endoscopy |
JP3802508B2 (ja) * | 2003-04-21 | 2006-07-26 | アロカ株式会社 | 超音波診断装置 |
US7457444B2 (en) * | 2003-05-14 | 2008-11-25 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for fast automatic centerline extraction for virtual endoscopy |
US7577298B2 (en) * | 2003-06-20 | 2009-08-18 | Siemens Medical Solutions Usa, Inc. | Method and system for response image feature collection and candidate summit, surface, and core estimation |
US7349563B2 (en) * | 2003-06-25 | 2008-03-25 | Siemens Medical Solutions Usa, Inc. | System and method for polyp visualization |
US7499572B2 (en) * | 2003-06-30 | 2009-03-03 | The Salk Institute For Biological Studies | Surface segmentation from luminance and color differences |
ITMO20030195A1 (it) * | 2003-06-30 | 2005-01-01 | Associazione Angela Serra Per La Ri Cerca Sul Canc | Dispositivo per la raccolta, l'elaborazione |
US7822461B2 (en) * | 2003-07-11 | 2010-10-26 | Siemens Medical Solutions Usa, Inc. | System and method for endoscopic path planning |
WO2005017830A1 (en) * | 2003-08-04 | 2005-02-24 | Siemens Corporate Research, Inc. | Virtual organ unfolding for visualization |
US7397936B2 (en) * | 2003-08-13 | 2008-07-08 | Siemens Medical Solutions Usa, Inc. | Method and system for wavelet based detection of colon polyps |
DE10339979B4 (de) * | 2003-08-29 | 2005-11-17 | Tomtec Imaging Systems Gmbh | Verfahren und Vorrichtung zur Darstellung eines vorbestimmbaren Bereichs in mehrdimensionalen Datensätzen |
JP4306380B2 (ja) * | 2003-09-10 | 2009-07-29 | 株式会社日立メディコ | 医用画像表示方法及び装置 |
US8276091B2 (en) * | 2003-09-16 | 2012-09-25 | Ram Consulting | Haptic response system and method of use |
US7643025B2 (en) * | 2003-09-30 | 2010-01-05 | Eric Belk Lange | Method and apparatus for applying stereoscopic imagery to three-dimensionally defined substrates |
US20050074150A1 (en) * | 2003-10-03 | 2005-04-07 | Andrew Bruss | Systems and methods for emulating an angiogram using three-dimensional image data |
US20050152588A1 (en) * | 2003-10-28 | 2005-07-14 | University Of Chicago | Method for virtual endoscopic visualization of the colon by shape-scale signatures, centerlining, and computerized detection of masses |
EP2189107A3 (en) * | 2003-10-31 | 2010-09-15 | Olympus Corporation | Insertion support system |
US20050094262A1 (en) * | 2003-11-05 | 2005-05-05 | Visx, Incorporated | Microscope magnification sensor |
WO2005055137A2 (en) * | 2003-11-26 | 2005-06-16 | Viatronix Incorporated | Vessel segmentation using vesselness and edgeness |
US20050110791A1 (en) * | 2003-11-26 | 2005-05-26 | Prabhu Krishnamoorthy | Systems and methods for segmenting and displaying tubular vessels in volumetric imaging data |
US7515743B2 (en) * | 2004-01-08 | 2009-04-07 | Siemens Medical Solutions Usa, Inc. | System and method for filtering a medical image |
JP4206044B2 (ja) * | 2004-01-20 | 2009-01-07 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | カルシウム・スコア測定方法および装置 |
DE102004008979B4 (de) * | 2004-02-24 | 2006-12-28 | Siemens Ag | Verfahren zur Filterung tomographischer 3D-Darstellungen nach erfolgter Rekonstruktion von Volumendaten |
US20050197558A1 (en) * | 2004-03-04 | 2005-09-08 | Williams James P. | System and method for performing a virtual endoscopy in a branching structure |
EP1743279A4 (en) * | 2004-04-08 | 2010-09-01 | Yeda Res & Dev | THREE-DAY LUNG CANCER DETECTION, DIAGNOSIS AND EVALUATION OF FORECAST |
US7447344B2 (en) * | 2004-04-16 | 2008-11-04 | Siemens Medical Solutions Usa, Inc. | System and method for visualization of pulmonary emboli from high-resolution computed tomography images |
EP1755446A4 (en) * | 2004-05-14 | 2007-10-31 | James V Manzione | COMBINATION OF IMAGING TECHNOLOGIES WITH MULTIPLE MODALITIES |
EP1769390B1 (en) * | 2004-06-04 | 2014-12-03 | Stereotaxis, Inc. | User interface for remote control of medical devices |
DE102004027709B4 (de) * | 2004-06-07 | 2006-08-10 | Siemens Ag | Verfahren der virtuellen Endoskopie zur medizinischen 3D-Bilddarstellung und -verarbeitung, Computertomografiegerät, Arbeitsstation und Computerprogrammprodukt |
DE102004027708B4 (de) * | 2004-06-07 | 2006-07-27 | Siemens Ag | Verfahren zur medizinischen 3D-Bilddarstellung und -verarbeitung, Computertomografiegerät, Arbeitsstation und Computerprogrammprodukt |
SE528068C2 (sv) * | 2004-08-19 | 2006-08-22 | Jan Erik Solem Med Jsolutions | Igenkänning av 3D föremål |
DE102004043694B4 (de) * | 2004-09-09 | 2006-09-28 | Siemens Ag | Verfahren zur Segmentierung anatomischer Strukturen aus 3D-Bilddaten unter Nutzung topologischer Information |
US7583829B2 (en) * | 2004-09-13 | 2009-09-01 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for embolism analysis |
US7764838B2 (en) * | 2004-09-14 | 2010-07-27 | Siemens Medical Solutions Usa, Inc. | System and method for extracting an object of interest from an image using a robust active shape model |
US7734119B2 (en) * | 2004-09-21 | 2010-06-08 | General Electric Company | Method and system for progressive multi-resolution three-dimensional image reconstruction using region of interest information |
US20120182403A1 (en) * | 2004-09-30 | 2012-07-19 | Eric Belk Lange | Stereoscopic imaging |
US7127095B2 (en) * | 2004-10-15 | 2006-10-24 | The Brigham And Women's Hospital, Inc. | Factor analysis in medical imaging |
US7702137B2 (en) * | 2004-11-10 | 2010-04-20 | M2S, Inc. | Anatomical visualization and measurement system |
EP1820156B1 (en) * | 2004-11-29 | 2017-09-13 | Koninklijke Philips N.V. | Multi-component vessel segmentation |
US10726741B2 (en) | 2004-11-30 | 2020-07-28 | The Regents Of The University Of California | System and method for converting handheld diagnostic ultrasound systems into ultrasound training systems |
US11627944B2 (en) | 2004-11-30 | 2023-04-18 | The Regents Of The University Of California | Ultrasound case builder system and method |
US8480404B2 (en) * | 2004-11-30 | 2013-07-09 | Eric A. Savitsky | Multimodal ultrasound training system |
US10026338B2 (en) | 2004-11-30 | 2018-07-17 | The Regents Of The University Of California | Embedded motion sensing technology for integration within commercial ultrasound probes |
TWI235041B (en) * | 2004-12-09 | 2005-07-01 | Univ Tsinghua | Characteristic points automatically identification method for three-dimensional space scanning data of human body |
DE102004060931A1 (de) * | 2004-12-17 | 2006-07-06 | Siemens Ag | Verfahren zur Vorbereitung der Begutachtung von tomographischen Colon-Aufnahmen |
US8066759B2 (en) | 2005-02-04 | 2011-11-29 | Boston Scientific Scimed, Inc. | Resonator for medical device |
CN101116110B (zh) * | 2005-02-08 | 2013-03-27 | 皇家飞利浦电子股份有限公司 | 医学图像浏览协议 |
DE102005005687A1 (de) * | 2005-02-08 | 2006-08-17 | Siemens Ag | Bildgebendes Diagnosesystem |
US8094906B2 (en) * | 2005-03-17 | 2012-01-10 | Algotec Systems Ltd. | Bone segmentation |
WO2006104468A1 (en) * | 2005-03-31 | 2006-10-05 | Agency For Science, Technology And Research | Method and apparatus for image segmentation |
US20060228003A1 (en) * | 2005-04-06 | 2006-10-12 | Silverstein D A | Method and apparatus for detection of optical elements |
EP1881788A4 (en) * | 2005-05-06 | 2009-11-04 | Yeda Res & Dev | FIGURE AND MOTION ANALYSIS OF ERYTHROCYTES IN BLOOD VESSELS RELATING TO THE CARDIAC CYCLE |
US7595469B2 (en) * | 2005-05-24 | 2009-09-29 | Boston Scientific Scimed, Inc. | Resonator for medical device |
US7532214B2 (en) * | 2005-05-25 | 2009-05-12 | Spectra Ab | Automated medical image visualization using volume rendering with local histograms |
US7889897B2 (en) * | 2005-05-26 | 2011-02-15 | Siemens Medical Solutions Usa, Inc. | Method and system for displaying unseen areas in guided two dimensional colon screening |
US8249687B2 (en) * | 2005-06-02 | 2012-08-21 | Vital Images, Inc. | Systems and methods for virtual identification of polyps |
EP1903943A1 (en) * | 2005-06-21 | 2008-04-02 | Koninklijke Philips Electronics N.V. | Method and device for imaging a blood vessel |
US20070019778A1 (en) * | 2005-07-22 | 2007-01-25 | Clouse Melvin E | Voxel histogram analysis for measurement of plaque |
US7279664B2 (en) | 2005-07-26 | 2007-10-09 | Boston Scientific Scimed, Inc. | Resonator for medical device |
US7379062B2 (en) * | 2005-08-01 | 2008-05-27 | Barco Nv | Method for determining a path along a biological object with a lumen |
US8583220B2 (en) * | 2005-08-02 | 2013-11-12 | Biosense Webster, Inc. | Standardization of catheter-based treatment for atrial fibrillation |
US7877128B2 (en) * | 2005-08-02 | 2011-01-25 | Biosense Webster, Inc. | Simulation of invasive procedures |
US7304277B2 (en) | 2005-08-23 | 2007-12-04 | Boston Scientific Scimed, Inc | Resonator with adjustable capacitor for medical device |
WO2007025081A2 (en) * | 2005-08-24 | 2007-03-01 | Traxtal Inc. | System, method and devices for navigated flexible endoscopy |
US7524282B2 (en) | 2005-08-29 | 2009-04-28 | Boston Scientific Scimed, Inc. | Cardiac sleeve apparatus, system and method of use |
US7773789B2 (en) * | 2005-08-30 | 2010-08-10 | Siemens Medical Solutions Usa, Inc. | Probabilistic minimal path for automated esophagus segmentation |
US7668342B2 (en) | 2005-09-09 | 2010-02-23 | Carl Zeiss Meditec, Inc. | Method of bioimage data processing for revealing more meaningful anatomic features of diseased tissues |
EP1942800B1 (en) * | 2005-09-09 | 2011-09-07 | Given Imaging Ltd. | Concurrent transfer and processing and real time viewing of in-vivo images |
WO2007033352A2 (en) * | 2005-09-12 | 2007-03-22 | 3Vr Security, Inc. | Storage of video analysis data for real-time alerting and forensic analysis |
US7702153B2 (en) * | 2005-10-07 | 2010-04-20 | Siemens Medical Solutions Usa, Inc. | Systems and methods for segmenting object of interest from medical image |
US7680314B2 (en) * | 2005-10-17 | 2010-03-16 | Siemens Medical Solutions Usa, Inc. | Devices, systems, and methods for improving image consistency |
ATE540616T1 (de) * | 2005-10-18 | 2012-01-15 | Koninkl Philips Electronics Nv | Patientienten-scanzeitoptimierung zur pet/spect- bildgebung |
US7423496B2 (en) | 2005-11-09 | 2008-09-09 | Boston Scientific Scimed, Inc. | Resonator with adjustable capacitance for medical device |
CN101454800A (zh) * | 2005-11-09 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | 用于医学成像的自动粪便删除方法 |
US7983462B2 (en) * | 2005-11-22 | 2011-07-19 | Purdue Research Foundation | Methods and systems for improving quality of an image |
US7590272B2 (en) * | 2005-11-23 | 2009-09-15 | Vital Images, Inc. | Colon characteristic path registration |
US8452061B2 (en) * | 2005-11-30 | 2013-05-28 | The Research Foundation Of State University Of New York | Electronic colon cleansing method for virtual colonoscopy |
US8050470B2 (en) * | 2005-12-07 | 2011-11-01 | Siemens Medical Solutions Usa, Inc. | Branch extension method for airway segmentation |
US8081180B2 (en) * | 2006-11-17 | 2011-12-20 | University Of Washington | Function-based representation of N-dimensional structures |
US8401264B2 (en) | 2005-12-08 | 2013-03-19 | University Of Washington | Solid modeling based on volumetric scans |
US7650179B2 (en) * | 2005-12-09 | 2010-01-19 | Siemens Aktiengesellschaft | Computerized workflow method for stent planning and stenting procedure |
US20070163114A1 (en) * | 2006-01-13 | 2007-07-19 | General Electric Company | Methods for fabricating components |
DE102006002037A1 (de) * | 2006-01-16 | 2007-07-19 | Siemens Ag | Verfahren zur Bearbeitung diagnostischer Bilddaten |
EP1825811B1 (en) * | 2006-02-27 | 2015-08-26 | Kabushiki Kaisha Toshiba | Image display apparatus and x-ray computed tomography apparatus |
US7768652B2 (en) | 2006-03-16 | 2010-08-03 | Carl Zeiss Meditec, Inc. | Methods for mapping tissue with optical coherence tomography data |
ATE499895T1 (de) * | 2006-04-12 | 2011-03-15 | Nassir Navab | Virtuelle penetrierende spiegelvorrichtung zur visualisierung von virtuellen objekten in angiografischen anwendungen |
CN101443815A (zh) * | 2006-05-11 | 2009-05-27 | 皇家飞利浦电子股份有限公司 | 重建图像的方法和装置 |
US8041129B2 (en) * | 2006-05-16 | 2011-10-18 | Sectra Ab | Image data set compression based on viewing parameters for storing medical image data from multidimensional data sets, related systems, methods and computer products |
US7825937B1 (en) | 2006-06-16 | 2010-11-02 | Nvidia Corporation | Multi-pass cylindrical cube map blur |
US20080009712A1 (en) * | 2006-06-16 | 2008-01-10 | Adams Mark L | Apparatus and Methods for Maneuvering a Therapeutic Tool Within a Body Lumen |
US8023703B2 (en) * | 2006-07-06 | 2011-09-20 | The United States of America as represented by the Secretary of the Department of Health and Human Services, National Institues of Health | Hybrid segmentation of anatomical structure |
US20080027917A1 (en) * | 2006-07-31 | 2008-01-31 | Siemens Corporate Research, Inc. | Scalable Semantic Image Search |
WO2008017051A2 (en) | 2006-08-02 | 2008-02-07 | Inneroptic Technology Inc. | System and method of providing real-time dynamic imagery of a medical procedure site using multiple modalities |
US7929741B2 (en) * | 2006-08-14 | 2011-04-19 | Siemens Medical Solutions Usa, Inc. | System and method for automated detection of mucus plugs within bronchial tree in MSCT images |
US7756310B2 (en) * | 2006-09-14 | 2010-07-13 | General Electric Company | System and method for segmentation |
EP1914681B1 (en) * | 2006-10-19 | 2010-09-15 | BrainLAB AG | Smooth gray-level based surface interpolation for anisotropic data sets |
US20080144907A1 (en) * | 2006-10-30 | 2008-06-19 | Ein-Yiao Shen | Devices and Methods for Reconstructing Three Dimensional Images |
US7831096B2 (en) * | 2006-11-17 | 2010-11-09 | General Electric Company | Medical navigation system with tool and/or implant integration into fluoroscopic image projections and method of use |
US8077939B2 (en) * | 2006-11-22 | 2011-12-13 | General Electric Company | Methods and systems for enhanced plaque visualization |
US7830381B2 (en) * | 2006-12-21 | 2010-11-09 | Sectra Ab | Systems for visualizing images using explicit quality prioritization of a feature(s) in multidimensional image data sets, related methods and computer products |
US7941213B2 (en) * | 2006-12-28 | 2011-05-10 | Medtronic, Inc. | System and method to evaluate electrode position and spacing |
KR20090107536A (ko) | 2007-01-22 | 2009-10-13 | 캘리포니아 인스티튜트 오브 테크놀로지 | 정량적 3-d 이미징을 위한 방법 |
WO2008091583A2 (en) * | 2007-01-23 | 2008-07-31 | Dtherapeutics, Llc | Image-based extraction for vascular trees |
US8023710B2 (en) * | 2007-02-12 | 2011-09-20 | The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services | Virtual colonoscopy via wavelets |
US20100087752A1 (en) * | 2007-03-22 | 2010-04-08 | Luya Li | Assessment of urinary system function by pattern matching |
CN101681433B (zh) * | 2007-03-26 | 2014-01-08 | 汤姆森特许公司 | 通过彩色图像分割和形状分析检测足球视频中感兴趣对象的方法和装置 |
WO2008130906A1 (en) * | 2007-04-17 | 2008-10-30 | Mikos, Ltd. | System and method for using three dimensional infrared imaging to provide psychological profiles of individuals |
CN101711125B (zh) | 2007-04-18 | 2016-03-16 | 美敦力公司 | 针对非荧光镜植入的长期植入性有源固定医疗电子导联 |
US10950338B2 (en) * | 2007-04-20 | 2021-03-16 | Robert Edwin Douglas | Method and apparatus for generating an artificial intelligence 3D dataset and performing interactive manipulation and rendering of the dataset |
WO2008133959A1 (en) | 2007-04-23 | 2008-11-06 | California Institute Of Technology | Single-lens 3-d imaging device using a polarization-coded aperture maks combined with a polarization-sensitive sensor |
US20080275467A1 (en) * | 2007-05-02 | 2008-11-06 | Siemens Corporate Research, Inc. | Intraoperative guidance for endovascular interventions via three-dimensional path planning, x-ray fluoroscopy, and image overlay |
US8167803B2 (en) | 2007-05-16 | 2012-05-01 | Verathon Inc. | System and method for bladder detection using harmonic imaging |
US8989468B2 (en) * | 2007-05-25 | 2015-03-24 | Definiens Ag | Generating an anatomical model using a rule-based segmentation and classification process |
JP4563421B2 (ja) * | 2007-05-28 | 2010-10-13 | ザイオソフト株式会社 | 画像処理方法及び画像処理プログラム |
US20080300994A1 (en) * | 2007-05-31 | 2008-12-04 | Caterpillar Inc. | System and method for conveying information |
US8184885B2 (en) * | 2007-07-24 | 2012-05-22 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | System and method for visualizing a structure of interest |
US9135401B2 (en) * | 2007-08-02 | 2015-09-15 | Siemens Medical Solutions Usa, Inc. | Information-theoretic view of the scheduling problem in whole-body computer aided detection/diagnosis (CAD) |
US8126226B2 (en) * | 2007-09-20 | 2012-02-28 | General Electric Company | System and method to generate a selected visualization of a radiological image of an imaged subject |
US8135189B2 (en) * | 2007-10-03 | 2012-03-13 | Siemens Medical Solutions Usa, Inc. | System and method for organ segmentation using surface patch classification in 2D and 3D images |
JP5543922B2 (ja) | 2007-11-02 | 2014-07-09 | コーニンクレッカ フィリップス エヌ ヴェ | 強調された冠状動脈表示 |
WO2009064715A1 (en) * | 2007-11-14 | 2009-05-22 | Auckland Uniservices Limited | Method for multi-scale meshing of branching biological structures |
US8144949B2 (en) * | 2007-11-15 | 2012-03-27 | Carestream Health, Inc. | Method for segmentation of lesions |
US11264139B2 (en) | 2007-11-21 | 2022-03-01 | Edda Technology, Inc. | Method and system for adjusting interactive 3D treatment zone for percutaneous treatment |
CN101877996B (zh) * | 2007-11-21 | 2014-10-15 | 美国医软科技公司 | 通过计算机对三维体图像进行处理的方法和用于交互式经皮术前手术规划的系统 |
US9366742B2 (en) * | 2007-12-11 | 2016-06-14 | Koninklijke Philips N.V. | Reducing motion artefacts in MRI |
CN109584996A (zh) * | 2007-12-13 | 2019-04-05 | 皇家飞利浦电子股份有限公司 | 一系列图像中的导航 |
US20090287076A1 (en) * | 2007-12-18 | 2009-11-19 | Boyden Edward S | System, devices, and methods for detecting occlusions in a biological subject |
US8280484B2 (en) * | 2007-12-18 | 2012-10-02 | The Invention Science Fund I, Llc | System, devices, and methods for detecting occlusions in a biological subject |
US9717896B2 (en) | 2007-12-18 | 2017-08-01 | Gearbox, Llc | Treatment indications informed by a priori implant information |
US20090287120A1 (en) | 2007-12-18 | 2009-11-19 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Circulatory monitoring systems and methods |
US8636670B2 (en) | 2008-05-13 | 2014-01-28 | The Invention Science Fund I, Llc | Circulatory monitoring systems and methods |
US9672471B2 (en) | 2007-12-18 | 2017-06-06 | Gearbox Llc | Systems, devices, and methods for detecting occlusions in a biological subject including spectral learning |
US8403862B2 (en) * | 2007-12-20 | 2013-03-26 | Yeda Research And Development Co. Ltd. | Time-based imaging |
US20090198124A1 (en) * | 2008-01-31 | 2009-08-06 | Ralf Adamus | Workflow to enhance a transjugular intrahepatic portosystemic shunt procedure |
US8663120B2 (en) | 2008-04-18 | 2014-03-04 | Regents Of The University Of Minnesota | Method and apparatus for mapping a structure |
US8494608B2 (en) | 2008-04-18 | 2013-07-23 | Medtronic, Inc. | Method and apparatus for mapping a structure |
US8340751B2 (en) | 2008-04-18 | 2012-12-25 | Medtronic, Inc. | Method and apparatus for determining tracking a virtual point defined relative to a tracked member |
US8260395B2 (en) | 2008-04-18 | 2012-09-04 | Medtronic, Inc. | Method and apparatus for mapping a structure |
US8532734B2 (en) | 2008-04-18 | 2013-09-10 | Regents Of The University Of Minnesota | Method and apparatus for mapping a structure |
US8839798B2 (en) | 2008-04-18 | 2014-09-23 | Medtronic, Inc. | System and method for determining sheath location |
CA2665215C (en) * | 2008-05-06 | 2015-01-06 | Intertape Polymer Corp. | Edge coatings for tapes |
US8077948B2 (en) * | 2008-05-06 | 2011-12-13 | Carestream Health, Inc. | Method for editing 3D image segmentation maps |
WO2009147671A1 (en) | 2008-06-03 | 2009-12-10 | Superdimension Ltd. | Feature-based registration method |
US8041095B2 (en) * | 2008-06-11 | 2011-10-18 | Siemens Aktiengesellschaft | Method and apparatus for pretreatment planning of endovascular coil placement |
EP2313868B1 (en) * | 2008-07-18 | 2016-03-16 | Vorum Research Corporation | Method, apparatus, signals and media for producing a computer representation of a three-dimensional surface of an appliance for a living body |
US8200466B2 (en) | 2008-07-21 | 2012-06-12 | The Board Of Trustees Of The Leland Stanford Junior University | Method for tuning patient-specific cardiovascular simulations |
CA2957778C (en) * | 2008-08-07 | 2022-12-13 | Jongtae Yuk | Device, system, and method to measure abdominal aortic aneurysm diameter |
WO2010027391A2 (en) | 2008-08-27 | 2010-03-11 | California Institute Of Technology | Method and device for high-resolution three-dimensional imaging which obtains camera pose using defocusing |
EP2350999A4 (en) * | 2008-09-25 | 2017-04-05 | CAE Healthcare Canada Inc. | Simulation of medical imaging |
US8301224B2 (en) * | 2008-10-09 | 2012-10-30 | Siemens Aktiengesellschaft | System and method for automatic, non-invasive diagnosis of pulmonary hypertension and measurement of mean pulmonary arterial pressure |
US8175681B2 (en) | 2008-12-16 | 2012-05-08 | Medtronic Navigation Inc. | Combination of electromagnetic and electropotential localization |
US8391600B2 (en) * | 2009-01-21 | 2013-03-05 | Analogic Corporation | Histogram-based compound object separation |
DE102009006416B3 (de) * | 2009-01-28 | 2010-08-26 | Siemens Aktiengesellschaft | Verfahren und Verzweigungs-Ermittlungseinrichtung zur Ermittlung einer Verzweigungsstelle sowie einer Mittellinie innerhalb eines Hohlorgans nebst zugehörigem Computerprogramm |
US11464578B2 (en) | 2009-02-17 | 2022-10-11 | Inneroptic Technology, Inc. | Systems, methods, apparatuses, and computer-readable media for image management in image-guided medical procedures |
US8690776B2 (en) | 2009-02-17 | 2014-04-08 | Inneroptic Technology, Inc. | Systems, methods, apparatuses, and computer-readable media for image guided surgery |
US9405886B2 (en) | 2009-03-17 | 2016-08-02 | The Board Of Trustees Of The Leland Stanford Junior University | Method for determining cardiovascular information |
EP3427687A1 (en) * | 2009-05-14 | 2019-01-16 | Covidien LP | Automatic registration technique |
US8773507B2 (en) * | 2009-08-11 | 2014-07-08 | California Institute Of Technology | Defocusing feature matching system to measure camera pose with interchangeable lens cameras |
US8331669B2 (en) * | 2009-08-25 | 2012-12-11 | Siemens Aktiengesellschaft | Method and system for interactive segmentation using texture and intensity cues |
US8494613B2 (en) | 2009-08-31 | 2013-07-23 | Medtronic, Inc. | Combination localization system |
US8494614B2 (en) | 2009-08-31 | 2013-07-23 | Regents Of The University Of Minnesota | Combination localization system |
US8446934B2 (en) * | 2009-08-31 | 2013-05-21 | Texas Instruments Incorporated | Frequency diversity and phase rotation |
US8355774B2 (en) | 2009-10-30 | 2013-01-15 | Medtronic, Inc. | System and method to evaluate electrode position and spacing |
US8819591B2 (en) * | 2009-10-30 | 2014-08-26 | Accuray Incorporated | Treatment planning in a virtual environment |
BR112012012818A2 (pt) * | 2009-11-27 | 2016-08-16 | Dog Microsystems Inc | método e sistema para determinar uma estimativa do apoio topológico de uma estrutura tubular e sua utilização em endoscopia virtual |
US20110188720A1 (en) * | 2010-02-02 | 2011-08-04 | General Electric Company | Method and system for automated volume of interest segmentation |
US8977019B2 (en) * | 2010-02-11 | 2015-03-10 | The Regents Of The University Of Michigan | Methods for microcalcification detection of breast cancer on digital tomosynthesis mammograms |
US20110213260A1 (en) * | 2010-02-26 | 2011-09-01 | Pacesetter, Inc. | Crt lead placement based on optimal branch selection and optimal site selection |
US9060673B2 (en) | 2010-04-28 | 2015-06-23 | Given Imaging Ltd. | System and method for displaying portions of in-vivo images |
US8908940B1 (en) * | 2010-04-29 | 2014-12-09 | Mim Software, Inc. | System and method of applying an arbitrary angle to reformat medical images |
US8600143B1 (en) * | 2010-05-20 | 2013-12-03 | Kla-Tencor Corporation | Method and system for hierarchical tissue analysis and classification |
JP2012016575A (ja) * | 2010-06-07 | 2012-01-26 | Toshiba Corp | 画像処理装置及び医用画像診断装置 |
MY152058A (en) * | 2010-06-21 | 2014-08-15 | Univ Putra Malaysia | A method of constructing at least one 3 dimensional image |
US20120022366A1 (en) * | 2010-07-21 | 2012-01-26 | Marcus Pfister | Registration of aorta to patient via two 2d images for placement of a stent |
US8315812B2 (en) | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
US8157742B2 (en) | 2010-08-12 | 2012-04-17 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
CN103221975B (zh) | 2010-09-03 | 2017-04-19 | 加州理工学院 | 三维成像系统 |
US9589204B2 (en) * | 2010-09-20 | 2017-03-07 | Koninklijke Philips N.V. | Quantification of a characteristic of a lumen of a tubular structure |
US20130170726A1 (en) * | 2010-09-24 | 2013-07-04 | The Research Foundation Of State University Of New York | Registration of scanned objects obtained from different orientations |
US8526700B2 (en) | 2010-10-06 | 2013-09-03 | Robert E. Isaacs | Imaging system and method for surgical and interventional medical procedures |
US9785246B2 (en) | 2010-10-06 | 2017-10-10 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US11231787B2 (en) | 2010-10-06 | 2022-01-25 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
KR101715780B1 (ko) * | 2010-10-11 | 2017-03-13 | 삼성전자주식회사 | 복셀 맵 생성기 및 그 방법 |
US8903144B2 (en) * | 2010-12-01 | 2014-12-02 | Olympus Corporation | Endoscope apparatus and method of measuring object |
US20140125663A1 (en) * | 2010-12-03 | 2014-05-08 | Institute of Automation, Chinese Academy of Scienc | 3d model shape analysis method based on perception information |
DE102011003064A1 (de) * | 2010-12-29 | 2012-07-05 | Robert Bosch Gmbh | Verfahren zur Bearbeitung einer Fläche mittels eines Roboter-Fahrzeugs |
CA2863675C (en) | 2011-02-11 | 2019-04-23 | E4 Endeavors, Inc. | System and method for modeling a biopsy specimen |
US20120259224A1 (en) * | 2011-04-08 | 2012-10-11 | Mon-Ju Wu | Ultrasound Machine for Improved Longitudinal Tissue Analysis |
US9466116B2 (en) * | 2011-04-15 | 2016-10-11 | Siemens Aktiengesellschaft | Method and system for separating tissue classes in magnetic resonance images |
US9226654B2 (en) | 2011-04-29 | 2016-01-05 | Carl Zeiss Meditec, Inc. | Systems and methods for automated classification of abnormalities in optical coherence tomography images of the eye |
KR101805624B1 (ko) * | 2011-08-29 | 2017-12-08 | 삼성전자주식회사 | 장기 모델 영상 생성 방법 및 장치 |
US8417005B1 (en) * | 2011-09-15 | 2013-04-09 | Sunnybrook Health Sciences Centre | Method for automatic three-dimensional segmentation of magnetic resonance images |
US20130100118A1 (en) * | 2011-10-21 | 2013-04-25 | Digital Artforms, Inc. | Systems and methods for human-computer interaction using a two handed interface |
US20130104084A1 (en) * | 2011-10-21 | 2013-04-25 | Digital Artforms, Inc. | Systems and methods for human-computer interaction using a two handed interface |
CN103098090B (zh) * | 2011-12-21 | 2015-01-21 | 中国科学院自动化研究所 | 多参数三维磁共振图像脑肿瘤分割方法 |
US9070203B2 (en) * | 2012-02-08 | 2015-06-30 | Mrl Materials Resources Llc | Identification and quantification of microtextured regions in materials with ordered crystal structure |
US8983232B2 (en) * | 2012-03-29 | 2015-03-17 | Dcg Systems, Inc. | Method for evaluating the centerline of an arbitrarily shaped object |
US8548778B1 (en) | 2012-05-14 | 2013-10-01 | Heartflow, Inc. | Method and system for providing information from a patient-specific model of blood flow |
US11631342B1 (en) | 2012-05-25 | 2023-04-18 | The Regents Of University Of California | Embedded motion sensing technology for integration within commercial ultrasound probes |
DE102012214351B4 (de) * | 2012-08-13 | 2014-11-06 | Siemens Aktiengesellschaft | Angiographisches Untersuchungsverfahren eines Gefäßsystems |
JP5777070B2 (ja) * | 2012-09-14 | 2015-09-09 | 富士フイルム株式会社 | 領域抽出装置、領域抽出方法および領域抽出プログラム |
US9858668B2 (en) | 2012-10-05 | 2018-01-02 | Volcano Corporation | Guidewire artifact removal in images |
JP2015535723A (ja) * | 2012-10-05 | 2015-12-17 | ナサニエル ジェイ. ケンプ, | パラメータ確立、再生、およびアーチファクト除去3次元撮像のための方法およびシステム |
US10070827B2 (en) | 2012-10-05 | 2018-09-11 | Volcano Corporation | Automatic image playback |
JP6219396B2 (ja) * | 2012-10-12 | 2017-10-25 | インテュイティブ サージカル オペレーションズ, インコーポレイテッド | 分岐した解剖学的構造における医療デバイスの位置決定 |
ITGE20130032A1 (it) * | 2013-03-19 | 2014-09-20 | Esaote Spa | Metodo e dispositivo di imaging del sistema cardiovascolare |
JP6026932B2 (ja) * | 2013-03-22 | 2016-11-16 | 富士フイルム株式会社 | 医用画像表示制御装置および方法並びにプログラム |
JP6188379B2 (ja) * | 2013-03-29 | 2017-08-30 | 東芝メディカルシステムズ株式会社 | 医用画像処理装置、医用画像処理システム、医用画像診断装置 |
WO2014168350A1 (ko) * | 2013-04-10 | 2014-10-16 | 재단법인 아산사회복지재단 | 폐동맥과 폐정맥을 구분하는 방법 및 이를 이용한 혈관의 정량화 방법 |
GB201307590D0 (en) | 2013-04-26 | 2013-06-12 | St Georges Hosp Medical School | Processing imaging data to obtain tissue type information |
CN103324783B (zh) * | 2013-05-30 | 2015-12-23 | 国家电网公司 | 基于边折叠的lod模型实时生成方法 |
JP6026357B2 (ja) * | 2013-06-13 | 2016-11-16 | 富士フイルム株式会社 | 仮想内視鏡画像生成装置および方法並びにプログラム |
KR101851221B1 (ko) | 2013-07-05 | 2018-04-25 | 삼성전자주식회사 | 초음파 영상 장치 및 그 제어 방법 |
US9626757B2 (en) | 2013-08-08 | 2017-04-18 | Washington University | System and method for the validation and quality assurance of computerized contours of human anatomy |
US10376715B2 (en) | 2013-08-08 | 2019-08-13 | Washington University | System and method for the validation and quality assurance of computerized contours of human anatomy |
US10380919B2 (en) | 2013-11-21 | 2019-08-13 | SonoSim, Inc. | System and method for extended spectrum ultrasound training using animate and inanimate training objects |
US9996954B2 (en) * | 2013-10-03 | 2018-06-12 | Covidien Lp | Methods and systems for dynamic display of a trace of a physiological parameter |
BR112016011777A2 (pt) | 2013-11-27 | 2017-08-08 | Gen Electric | Aparelhos de bocal de combustível |
US10451282B2 (en) | 2013-12-23 | 2019-10-22 | General Electric Company | Fuel nozzle structure for air assist injection |
JP6695801B2 (ja) | 2013-12-23 | 2020-05-20 | ゼネラル・エレクトリック・カンパニイ | 可撓性支持構造体を備えた燃料ノズル |
WO2015109251A1 (en) * | 2014-01-17 | 2015-07-23 | Truinject Medical Corp. | Injection site training system |
GB2522663B (en) * | 2014-01-31 | 2020-02-12 | Apical Ltd | A method of selecting a region of interest |
CN105095921B (zh) * | 2014-04-30 | 2019-04-30 | 西门子医疗保健诊断公司 | 用于处理尿液沉渣图像的待处理区块的方法和装置 |
US9962136B2 (en) * | 2014-06-04 | 2018-05-08 | The Board Of Trustees Of The Leland Stanford Junior University | Reconstructing time-varying position of individual radioactive cells directly from positron emission tomography (PET) measurements |
KR102238693B1 (ko) * | 2014-06-20 | 2021-04-09 | 삼성전자주식회사 | 포인트 클라우드에서 특징 영역을 추출하는 방법 및 장치 |
AU2015299765A1 (en) * | 2014-08-06 | 2017-02-16 | Commonwealth Scientific And Industrial Research Organisation | Representing an interior of a volume |
CN106796661B (zh) * | 2014-08-12 | 2020-08-25 | 曼蒂斯影像有限公司 | 投影光图案的系统、方法和计算机程序产品 |
US9901406B2 (en) | 2014-10-02 | 2018-02-27 | Inneroptic Technology, Inc. | Affected region display associated with a medical device |
US10188467B2 (en) | 2014-12-12 | 2019-01-29 | Inneroptic Technology, Inc. | Surgical guidance intersection display |
KR102364490B1 (ko) * | 2014-12-15 | 2022-02-18 | 삼성메디슨 주식회사 | 초음파 진단장치, 그에 따른 초음파 진단 방법 및 그에 따른 컴퓨터 판독 가능한 저장매체 |
US10242488B1 (en) * | 2015-03-02 | 2019-03-26 | Kentucky Imaging Technologies, LLC | One-sided transparency: a novel visualization for tubular objects |
JP6470604B2 (ja) * | 2015-03-25 | 2019-02-13 | キヤノン株式会社 | 画像処理装置及び画像処理方法及びプログラム |
US9741104B2 (en) * | 2015-05-18 | 2017-08-22 | Toshiba Medical Systems Corporation | Apparatus, method, and computer-readable medium for quad reconstruction using hybrid filter convolution and high dynamic range tone-mapping |
RU2600282C2 (ru) * | 2015-06-23 | 2016-10-20 | Федеральное государственное бюджетное учреждение "Национальный медицинский исследовательский радиологический центр" Министерства здравоохранения Российской Федерации (ФГБУ "НМИРЦ" Минздрава России) | Способ планирования анатомических сублобарных резекций легких у больных с периферическими объемными образованиями на основе кт-ангиопульмонографии |
US11600201B1 (en) | 2015-06-30 | 2023-03-07 | The Regents Of The University Of California | System and method for converting handheld diagnostic ultrasound systems into ultrasound training systems |
US9949700B2 (en) * | 2015-07-22 | 2018-04-24 | Inneroptic Technology, Inc. | Medical device approaches |
JP6632830B2 (ja) * | 2015-08-12 | 2020-01-22 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
WO2017031175A1 (en) * | 2015-08-17 | 2017-02-23 | University Of Maryland, Baltimore | Automated surgeon performance evaluation |
CN105231978B (zh) * | 2015-09-14 | 2017-03-22 | 袁非牛 | 一种引导式虚拟内窥镜导航方法 |
WO2017055403A1 (en) * | 2015-09-29 | 2017-04-06 | Institut National De La Sante Et De La Recherche Medicale (Inserm) | Device and system for generating ultrasonic waves in a target region of a soft solid and method for locally treating a tissue |
JP6615603B2 (ja) * | 2015-12-24 | 2019-12-04 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置および医用画像診断プログラム |
US11406264B2 (en) | 2016-01-25 | 2022-08-09 | California Institute Of Technology | Non-invasive measurement of intraocular pressure |
JP6639935B2 (ja) * | 2016-02-09 | 2020-02-05 | キヤノンメディカルシステムズ株式会社 | 画像処理装置およびmri装置 |
US9675319B1 (en) | 2016-02-17 | 2017-06-13 | Inneroptic Technology, Inc. | Loupe display |
US10849688B2 (en) | 2016-03-02 | 2020-12-01 | Truinject Corp. | Sensory enhanced environments for injection aid and social training |
US10278778B2 (en) | 2016-10-27 | 2019-05-07 | Inneroptic Technology, Inc. | Medical device navigation using a virtual 3D space |
WO2018136901A1 (en) | 2017-01-23 | 2018-07-26 | Truinject Corp. | Syringe dose and position measuring apparatus |
US10896628B2 (en) | 2017-01-26 | 2021-01-19 | SonoSim, Inc. | System and method for multisensory psychomotor skill training |
US10251013B2 (en) | 2017-06-08 | 2019-04-02 | Microsoft Technology Licensing, Llc | Audio propagation in a virtual environment |
US11259879B2 (en) | 2017-08-01 | 2022-03-01 | Inneroptic Technology, Inc. | Selective transparency to assist medical device navigation |
WO2019037558A1 (zh) | 2017-08-22 | 2019-02-28 | 优酷网络技术(北京)有限公司 | 图像处理方法和装置 |
GB2567636B (en) * | 2017-10-17 | 2021-11-10 | Perspectum Diagnostics Ltd | Method and apparatus for imaging an organ |
US10426424B2 (en) | 2017-11-21 | 2019-10-01 | General Electric Company | System and method for generating and performing imaging protocol simulations |
US11006852B2 (en) | 2017-12-11 | 2021-05-18 | Covidien Lp | Systems, methods, and computer-readable media of estimating thoracic cavity movement during respiration |
US11484365B2 (en) | 2018-01-23 | 2022-11-01 | Inneroptic Technology, Inc. | Medical image guidance |
US11100633B2 (en) * | 2018-06-13 | 2021-08-24 | Cosmo Artificial Intelligence—Al Limited | Systems and methods for processing real-time video from a medical image device and detecting objects in the video |
CN110738890A (zh) * | 2018-07-19 | 2020-01-31 | 苏州敏行医学信息技术有限公司 | 灌肠术智能训练方法及系统 |
CN109508087A (zh) * | 2018-09-25 | 2019-03-22 | 易念科技(深圳)有限公司 | 脑纹信号识别方法及终端设备 |
US11810473B2 (en) | 2019-01-29 | 2023-11-07 | The Regents Of The University Of California | Optical surface tracking for medical simulation |
US11495142B2 (en) | 2019-01-30 | 2022-11-08 | The Regents Of The University Of California | Ultrasound trainer with internal optical tracking |
US11151789B1 (en) | 2019-03-25 | 2021-10-19 | Kentucky Imaging Technologies | Fly-in visualization for virtual colonoscopy |
US11521316B1 (en) | 2019-04-03 | 2022-12-06 | Kentucky Imaging Technologies | Automatic extraction of interdental gingiva regions |
US11044400B1 (en) | 2019-04-03 | 2021-06-22 | Kentucky Imaging Technologies | Frame stitching in human oral cavity environment using intraoral camera |
DE102019212103A1 (de) * | 2019-08-13 | 2021-02-18 | Siemens Healthcare Gmbh | Surrogatmarker basierend auf medizinischen Bilddaten |
CN110853094B (zh) * | 2019-10-11 | 2020-12-29 | 北京文香信息技术有限公司 | 一种目标位置的阈值确定方法、装置、设备及存储介质 |
CN112998630B (zh) * | 2021-03-17 | 2022-07-29 | 安翰科技(武汉)股份有限公司 | 胶囊内窥镜的完备性自检方法、电子设备及可读存储介质 |
CN113017544B (zh) * | 2021-03-17 | 2022-07-29 | 安翰科技(武汉)股份有限公司 | 胶囊内窥镜的分区完备性自检方法、设备及可读存储介质 |
CN113610825B (zh) * | 2021-08-13 | 2022-03-29 | 推想医疗科技股份有限公司 | 术中影像的肋骨识别方法及系统 |
WO2023023170A1 (en) * | 2021-08-17 | 2023-02-23 | University Of Virginia Patent Foundation | Systems, methods, and computer readable media for parametric fdg pet quantification, segmentation and classification of abnormalities |
Family Cites Families (66)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DK273280A (da) * | 1979-06-28 | 1980-12-29 | Schering Ag | Trijoderede 5-aminoisophthalsyrederivater |
US4391280A (en) * | 1980-04-04 | 1983-07-05 | Miller Roscoe E | Enema apparata improvements relating to double contrast studies |
US4630203A (en) * | 1983-12-27 | 1986-12-16 | Thomas Szirtes | Contour radiography: a system for determining 3-dimensional contours of an object from its 2-dimensional images |
US4737921A (en) * | 1985-06-03 | 1988-04-12 | Dynamic Digital Displays, Inc. | Three dimensional medical image display system |
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 |
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 |
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 |
DE3611018C1 (de) * | 1986-03-27 | 1987-06-19 | Wiest Peter P | Vorrichtung zum Insufflieren von Gas |
US4751643A (en) * | 1986-08-04 | 1988-06-14 | General Electric Company | Method and apparatus for determining connected substructures within a body |
US4791567A (en) * | 1986-09-15 | 1988-12-13 | General Electric Company | Three dimensional connectivity system employing an equivalence schema for determining connected substructures within a body |
US4879668A (en) * | 1986-12-19 | 1989-11-07 | General Electric Company | Method of displaying internal surfaces of three-dimensional medical images |
US4823129A (en) * | 1987-02-24 | 1989-04-18 | Bison Instruments, Inc. | Analog-to-digital converter |
FR2613509B1 (fr) * | 1987-04-03 | 1989-06-09 | Thomson Cgr | Procede de calcul et de representation d'images de vues d'un objet |
FR2614163B1 (fr) * | 1987-04-17 | 1989-06-09 | Thomson Cgr | Procede de representation d'images de vues d'un objet |
US4831528A (en) * | 1987-11-09 | 1989-05-16 | General Electric Company | Apparatus and method for improvement of 3D images derived from tomographic data |
US5170347A (en) * | 1987-11-27 | 1992-12-08 | Picker International, Inc. | System to reformat images for three-dimensional display using unique spatial encoding and non-planar bisectioning |
US5038302A (en) * | 1988-07-26 | 1991-08-06 | The Research Foundation Of State University Of New York | Method of converting continuous three-dimensional geometrical representations into discrete three-dimensional voxel-based representations within a three-dimensional voxel-based system |
US5023072A (en) * | 1988-08-10 | 1991-06-11 | University Of New Mexico | Paramagnetic/superparamagnetic/ferromagnetic sucrose sulfate compositions for magnetic resonance imaging of the gastrointestinal tract |
US4991092A (en) * | 1988-08-12 | 1991-02-05 | The Regents Of The University Of California | Image processor for enhancing contrast between subregions of a region of interest |
US4993415A (en) * | 1988-08-19 | 1991-02-19 | Alliance Pharmaceutical Corp. | Magnetic resonance imaging with perfluorocarbon hydrides |
US4987554A (en) * | 1988-08-24 | 1991-01-22 | The Research Foundation Of State University Of New York | Method of converting continuous three-dimensional geometrical representations of polygonal objects into discrete three-dimensional voxel-based representations thereof within a three-dimensional voxel-based system |
FR2636752B1 (fr) * | 1988-09-16 | 1990-10-26 | Gen Electric Cgr | Procede et systeme de correction des defauts d'images d'un scanner dus aux deplacements de ce dernier |
US4984157A (en) * | 1988-09-21 | 1991-01-08 | General Electric Company | System and method for displaying oblique planar cross sections of a solid body using tri-linear interpolation to determine pixel position dataes |
US4985856A (en) * | 1988-11-10 | 1991-01-15 | The Research Foundation Of State University Of New York | Method and apparatus for storing, accessing, and processing voxel-based data |
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 |
US5101475A (en) * | 1989-04-17 | 1992-03-31 | The Research Foundation Of State University Of New York | Method and apparatus for generating arbitrary projections of three-dimensional voxel-based data |
FR2648304B1 (fr) * | 1989-06-12 | 1991-08-30 | Commissariat Energie Atomique | Procede de determination d'un espace a partir d'un espace discret connu pour la reconstruction d'images bi ou tridimensionnelles, dispositif de mise en oeuvre et application du procede |
JP2714164B2 (ja) * | 1989-07-31 | 1998-02-16 | 株式会社東芝 | 三次元画像表示装置 |
US5006109A (en) * | 1989-09-12 | 1991-04-09 | Donald D. Douglas | Method and device for controlling pressure, volumetric flow rate and temperature during gas insuffication procedures |
US5187658A (en) * | 1990-01-17 | 1993-02-16 | General Electric Company | System and method for segmenting internal structures contained within the interior region of a solid object |
US5086401A (en) * | 1990-05-11 | 1992-02-04 | International Business Machines Corporation | Image-directed robotic system for precise robotic surgery including redundant consistency checking |
US5047772A (en) * | 1990-06-04 | 1991-09-10 | General Electric Company | Digital error correction system for subranging analog-to-digital converters |
US5127037A (en) * | 1990-08-15 | 1992-06-30 | Bynum David K | Apparatus for forming a three-dimensional reproduction of an object from laminations |
US5204625A (en) * | 1990-12-20 | 1993-04-20 | General Electric Company | Segmentation of stationary and vascular surfaces in magnetic resonance imaging |
US5270926A (en) * | 1990-12-21 | 1993-12-14 | General Electric Company | Method and apparatus for reconstructing a three-dimensional computerized tomography (CT) image of an object from incomplete cone beam projection data |
US5166876A (en) * | 1991-01-16 | 1992-11-24 | General Electric Company | System and method for detecting internal structures contained within the interior region of a solid object |
DE4117117A1 (de) * | 1991-05-25 | 1992-11-26 | Hoehne Karl Heinz Prof Dr | Dreidimensionale darstellung von raeumlichen strukturen |
US5345490A (en) * | 1991-06-28 | 1994-09-06 | General Electric Company | Method and apparatus for converting computed tomography (CT) data into finite element models |
US5261404A (en) * | 1991-07-08 | 1993-11-16 | Mick Peter R | Three-dimensional mammal anatomy imaging system and method |
US5283837A (en) * | 1991-08-27 | 1994-02-01 | Picker International, Inc. | Accurate estimation of surface normals in 3-D data sets |
US5734384A (en) * | 1991-11-29 | 1998-03-31 | Picker International, Inc. | Cross-referenced sectioning and reprojection of diagnostic image volumes |
US5371778A (en) * | 1991-11-29 | 1994-12-06 | Picker International, Inc. | Concurrent display and adjustment of 3D projection, coronal slice, sagittal slice, and transverse slice images |
US5442733A (en) * | 1992-03-20 | 1995-08-15 | The Research Foundation Of State University Of New York | Method and apparatus for generating realistic images using a discrete representation |
US5295488A (en) * | 1992-08-05 | 1994-03-22 | General Electric Company | Method and apparatus for projecting diagnostic images from volumed diagnostic data |
US5322070A (en) * | 1992-08-21 | 1994-06-21 | E-Z-Em, Inc. | Barium enema insufflation system |
US5319549A (en) * | 1992-11-25 | 1994-06-07 | Arch Development Corporation | Method and system for determining geometric pattern features of interstitial infiltrates in chest images |
US5361763A (en) * | 1993-03-02 | 1994-11-08 | Wisconsin Alumni Research Foundation | Method for segmenting features in an image |
US5365927A (en) * | 1993-11-02 | 1994-11-22 | General Electric Company | Magnetic resonance imaging system with pointing device |
JP3483929B2 (ja) * | 1994-04-05 | 2004-01-06 | 株式会社日立製作所 | 3次元画像生成方法 |
US5458111A (en) * | 1994-09-06 | 1995-10-17 | William C. Bond | Computed tomographic colonoscopy |
US5920319A (en) | 1994-10-27 | 1999-07-06 | Wake Forest University | Automatic analysis in virtual endoscopy |
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 |
US5611025A (en) * | 1994-11-23 | 1997-03-11 | General Electric Company | Virtual internal cavity inspection system |
US6125194A (en) * | 1996-02-06 | 2000-09-26 | Caelum Research Corporation | Method and system for re-screening nodules in radiological images using multi-resolution processing, neural network, and image processing |
US5699799A (en) * | 1996-03-26 | 1997-12-23 | Siemens Corporate Research, Inc. | Automatic determination of the curved axis of a 3-D tube-shaped object in image volume |
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 |
US5971767A (en) | 1996-09-16 | 1999-10-26 | The Research Foundation Of State University Of New York | System and method for performing a three-dimensional virtual examination |
US6343936B1 (en) | 1996-09-16 | 2002-02-05 | The Research Foundation Of State University Of New York | System and method for performing a three-dimensional virtual examination, navigation and visualization |
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 |
US5986662A (en) * | 1996-10-16 | 1999-11-16 | Vital Images, Inc. | Advanced diagnostic viewer employing automated protocol selection for volume-rendered imaging |
US6130671A (en) * | 1997-11-26 | 2000-10-10 | Vital Images, Inc. | Volume rendering lighting using dot product methodology |
US6771262B2 (en) | 1998-11-25 | 2004-08-03 | Siemens Corporate Research, Inc. | System and method for volume rendering-based segmentation |
US20020164061A1 (en) | 2001-05-04 | 2002-11-07 | Paik David S. | Method for detecting shapes in medical images |
CN100405973C (zh) * | 2002-03-14 | 2008-07-30 | Netkisr有限公司 | 分析和显示计算机体层摄影术数据的方法 |
US20050152588A1 (en) * | 2003-10-28 | 2005-07-14 | University Of Chicago | Method for virtual endoscopic visualization of the colon by shape-scale signatures, centerlining, and computerized detection of masses |
US7515743B2 (en) * | 2004-01-08 | 2009-04-07 | Siemens Medical Solutions Usa, Inc. | System and method for filtering a medical image |
-
2001
- 2001-02-05 US US09/777,120 patent/US7194117B2/en not_active Expired - Lifetime
- 2001-10-01 WO PCT/US2001/030704 patent/WO2002029764A1/en active Application Filing
- 2001-10-01 AU AU2002211331A patent/AU2002211331A1/en not_active Abandoned
- 2001-10-01 EP EP01979359A patent/EP1342225A4/en not_active Withdrawn
- 2001-10-01 JP JP2002533259A patent/JP2004510515A/ja active Pending
-
2006
- 2006-12-20 US US11/613,283 patent/US7474776B2/en not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008503302A (ja) * | 2004-06-23 | 2008-02-07 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 3d管状オブジェクトのパラメータに関する情報を表示するための画像処理システム |
JP2006192254A (ja) * | 2005-01-12 | 2006-07-27 | Medison Co Ltd | 超音波診断システム及び方法 |
US7940995B2 (en) | 2005-01-12 | 2011-05-10 | Medison Co., Ltd. | Ultrasound diagnostic system for automatically detecting a boundary of a target object and method of using the same |
JP2007229464A (ja) * | 2006-02-28 | 2007-09-13 | Siemens Ag | 身体構成要素内の物質の濃度の決定方法および装置 |
JP2010504794A (ja) * | 2006-09-29 | 2010-02-18 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 突起検出方法、システム及びコンピュータプログラム |
JP2012050606A (ja) * | 2010-08-31 | 2012-03-15 | Fujifilm Corp | 内視鏡画像処理装置および方法、並びに、プログラム |
JP2021002402A (ja) * | 2020-10-06 | 2021-01-07 | 株式会社ソニー・インタラクティブエンタテインメント | 情報処理装置 |
JP7044846B2 (ja) | 2020-10-06 | 2022-03-30 | 株式会社ソニー・インタラクティブエンタテインメント | 情報処理装置 |
Also Published As
Publication number | Publication date |
---|---|
US20010031920A1 (en) | 2001-10-18 |
EP1342225A1 (en) | 2003-09-10 |
EP1342225A4 (en) | 2008-06-11 |
US20070167718A1 (en) | 2007-07-19 |
US7194117B2 (en) | 2007-03-20 |
US7474776B2 (en) | 2009-01-06 |
WO2002029764A1 (en) | 2002-04-11 |
WO2002029764A9 (en) | 2003-01-23 |
AU2002211331A1 (en) | 2002-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2004510515A (ja) | 器官のような物体の3次元仮想検査を行うシステムおよび方法 | |
US7486811B2 (en) | System and method for performing a three-dimensional virtual examination of objects, such as internal organs | |
JP4435430B2 (ja) | 3次元仮想細分化および検査を実施するシステムおよび方法 | |
US6343936B1 (en) | System and method for performing a three-dimensional virtual examination, navigation and visualization | |
US8682045B2 (en) | Virtual endoscopy with improved image segmentation and lesion detection | |
AU759501B2 (en) | Virtual endoscopy with improved image segmentation and lesion detection | |
US8600125B2 (en) | System and method for computer aided polyp detection | |
WO2000032106A1 (en) | Virtual endoscopy with improved image segmentation and lesion detection | |
JP2009056326A (ja) | 直接検出できない内視鏡検査における自動分析 | |
IL178768A (en) | System and method for optical mapping of texture properties from at least one optical image to a monochromatic array of information | |
Hong | Computer aided diagnosis for virtual endoscopy | |
MXPA01009388A (en) | System and method for performing a three-dimensional virtual segmentation and examination | |
MXPA01009387A (en) | System and method for performing a three-dimensional virtual examination, navigation and visualization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041001 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080115 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20080415 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20080422 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080715 |