JP4579990B2 - コンピュータ読み取り可能なプログラム記憶装置 - Google Patents

Description

（関連出願）
本出願は、引例によってその全体がここに組み込まれている、２００４年１０月２２日に出願された米国仮出願通番第６０／６２１，１９８号に対する優先権を主張する。

本発明は、画像解析に関し、特にボリュームサンプリングシステムおよび方法に関する。

［従来技術の説明］
特徴値を決定するための、または新しいデータを生成するためのデータ処理は典型的には、データ、例えば３次元コンピュータ断層撮影（ＣＴ）スキャンの全ボリュームよりむしろサブボリュームに対して行われる。典型的にはサブボリュームは、サブボリュームのサンプリング周波数がオリジナルボリュームのサンプリング周波数と一致しないようにオリジナルボリュームから抽出される。例えば、ある幾つかのアルゴリズムは、サブボリュームに等方性であることを要求するかもしれないが、オリジナルボリュームは等方性でないかも知れず、あるいは異なる分解能で等方性であるかもしれない。このような場合、サブボリュームは、予め定められたサンプリング周波数でオリジナルボリュームをサンプリングすることによってオリジナルボリュームから抽出される。このサンプリングのためにオリジナルボリューム内のボクセルは、サブボリューム内のボクセルに対して一対一対応を持たない。

多数のサブボリュームが１つのオリジナルボリュームから抽出されるとき、サブボリュームは典型的に重複している。図１は重複（オーバラップ）の一例を示す。オリジナルボリュームは境界１０１として示されている。このボリュームから２つのサブボリューム１０２、１０３が抽出される。これらのサブボリュームを処理するアルゴリズムは、関心領域（ＲＯＩ：ｒｅｇｉｏｎ−ｏｆ−ｉｎｔｅｒｅｓｔ）１０４上の特徴を抽出していると想定する。この関心領域１０４は、図１に示すように２つの重複しているサブボリューム１０２、１０３に共通であり得る。理想的には２つ以上のサブボリュームに関して決定された特徴値は、これらのサブボリュームがそれらの内部にＲＯＩの同じ部分を含んでいれば、同じであるはずである。同様に、サブボリューム内のＲＯＩの同じ部分が視覚化されるとき、これらは同じに見えるはずである。

実際には、２つのサブボリュームについて計算された特徴値は異なる可能性があり、２つのサブボリュームの視覚化は異なって見えることがある。相違の理由は下記の通りである。第１のサブボリュームが抽出されるとき、オリジナルボリュームはこのサブボリュームの所望のサンプリング周波数を得るようにサンプリングされる。第２のサブボリュームが抽出されるとき、このサンプリングプロセスが繰り返される。オリジナルボリューム内のこれら２つのサブボリュームのサンプリング位置を観察することによって、サブボリューム両者に関するサンプリング位置はサブボリュームが重複する領域内において一致しないことが決定され得る。例えば、図２を参照すると、第１、第２のサブボリュームに関するサンプリング位置２０２、２０３はそれぞれ等間隔であるが、一致していない。これら２つのサブボリュームの共通領域内のボクセル値は異なる。ボクセル値のこの差は、決定された特徴値の差という結果を生じる。

図２は、２つの重複しているサブボリュームのサンプリング位置を示す。ライン２０１はオリジナルサンプルの位置を表し、ライン２０２は第１のサブボリュームのサンプリング位置を表し、ライン２０３は第２のサブボリュームのサンプリング位置を表す。これら２つのサブボリュームのサンプリング位置２０２、２０３は一致しないことが観察できる。これは、２つのサブボリュームに関して同じでない特徴値につながる。

したがって、サブボリュームが重複する領域内においてサンプリング位置が同じである多数のサブボリュームの同期化されたサンプリングシステムおよび方法が必要である。

本発明によれば、特徴抽出のためにサブボリュームを位置合わせさせるためにコンピュータで実施される方法（以下「コンピュータ実施方法」という）は、データボリュームを受信するステップと、サブボリュームの複数のパラメータを供給するステップと、複数のパラメータにしたがって２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップとを含む。本方法は更に、データボリュームをサンプリングすることによって作り出された仮想ボリュームの格子に対する２つ以上のサブボリュームの位置ずれ量（すなわち位置合わせずれ量）を決定するステップと、各サブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせするために、前記位置ずれ量に関係する量によって各サブボリュームのサンプリング開始点をシフトするステップと、仮想ボリュームの格子上でサンプリングされた値を使用して、位置合わせされた２つ以上のサブボリュームにおける特徴を決定するステップとを含む。

複数のパラメータを供給するステップは、２つ以上のサブボリュームの各サブボリュームの中心点を選択するステップと、２つ以上のサブボリュームのサイズと形状とを選択するステップと、２つ以上のサブボリュームの内のボクセル間の間隔を選択するステップと、を含む。本方法は更に、ボクセル間隔をボリューム（体積）単位に変換するステップを含む。本方法は、各サブボリュームのサンプリング開始点とボリューム単位のボクセル間隔とを使用して、特徴を抽出する前に、ボリュームから２つ以上のサブボリュームを抽出するステップを含む。

２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップは、中心点とサブボリュームサイズとボクセル間隔とにしたがって実行される。

位置ずれ量（位置合わせずれ量）は仮想ボリュームの格子に対して決定され、仮想ボリューム格子の間隔はサブボリュームのボクセル間の間隔であり、また仮想ボリュームの格子の原点は世界座標系の原点と一致する。

本方法は、シフトされたサンプリング開始点を、各サブボリュームの（０，０，０）位置を原点として有するボリューム座標に変換するステップを含む。２つ以上のサブボリュームは異なる分解能を持つことができ、またデータボリュームは異なるモダリティを表す２つ以上のデータボリュームを含む。

各サブボリュームのサンプリング開始点をシフトするステップは世界座標系におけるサブボリュームの１つ又は２つの座標軸に対してだけ実行される。サブボリュームは２次元であって、視覚化ならびに特徴は２次元サブボリュームについて決定される。サンプリング開始点の位置合わせは３次元又はそれ以上の次元について決定される。多次元ボリュームの場合に本方法は１次元又は以上の次元を選択するステップを含み、サンプリング開始点の位置合わせは選択された次元について決定される。マルチモダリティボリュームの場合には、各サブボリュームのサンプリング開始点の位置合わせは予め登録された異なるモダリティに適用される。

本発明によれば、特徴抽出のために１つのデータボリュームをサンプリングすることによって抽出された２つ以上のサブボリュームを位置合わせさせる方法ステップを実行するためにコンピュータによって実行可能な命令のプログラムを記憶したコンピュータ読み取り可能なプログラム記憶装置が提供される。方法ステップは、データボリュームを受信するステップと、サブボリュームの複数のパラメータを供給するステップと、複数のパラメータにしたがって２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップとを含む。本方法は更に、データボリュームをサンプリングすることによって作り出された仮想ボリュームの格子に対する２つ以上のサブボリュームの位置ずれ量を決定するステップと、各サブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせするために、前記位置ずれ量に関係する量によって各サブボリュームのサンプリング開始点をシフトするステップと、仮想ボリュームの格子上でサンプリングされた値を用いて、位置合わせされた２つ以上のサブボリュームにおける特徴を決定するステップとを含む。

本発明によれば、視覚化のためにサブボリュームを位置合わせさせるためのコンピュータ実施方法は、データボリュームを受信するステップと、サブボリュームの複数のパラメータを供給するステップと、複数のパラメータにしたがって２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップとを含む。本方法は、データボリュームをサンプリングすることによって作り出された仮想ボリュームの格子に対する２つ以上のサブボリュームの位置ずれ量を決定するステップと、各サブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせするために、前記位置ずれ量に関係する量によって各サブボリュームのサンプリング開始点をシフトするステップと、位置合わせされたサンプリング開始点を有するサブボリュームを表示するステップとを含み、各サブボリューム内の関心領域は同じボクセル値を有する。

複数のパラメータを供給するステップは、２つ以上のサブボリュームの各サブボリュームの中心点を選択するステップと、２つ以上のサブボリュームのサイズと形状とを選択するステップと、２つ以上のサブボリューム内のボクセル間の間隔を選択するステップとを含む。本方法は、ボクセル間隔をボリューム（体積）単位に変換するステップを含む。本方法は、サンプリング開始点とボリューム単位のボクセル間隔とを使用して、サブボリュームを表示するステップの前に、ボリュームから２つ以上のサブボリュームを抽出するステップを含む。

本発明によれば、視覚化のために１つのデータボリュームをサンプリングすることによって抽出された２つ以上のサブボリュームを位置合わせさせる方法ステップを実行するためにコンピュータによって実行可能な命令のプログラムを記憶したコンピュータ読み取り可能なプログラム記憶装置が提供される。視覚化のために１つのデータボリュームをサンプリングすることによって抽出された２つ以上のサブボリュームを位置合わせさせる方法ステップは、データボリュームを受信するステップと、サブボリュームの複数のパラメータを供給するステップと、複数のパラメータにしたがって２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップとを含む。本方法は、データボリュームをサンプリングすることによって作り出された仮想ボリュームの格子に対する２つ以上のサブボリュームの位置ずれ量を決定するステップと、各サブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせするために、前記位置ずれ量に関係する量によって各サブボリュームのサンプリング開始点をシフトするステップと、位置合わせされたサンプリング開始点を有するサブボリュームを表示するステップとを含み、各サブボリューム内の関心領域は同じボクセル値を有する。

本発明の好適な実施形態は図面を参照して以下に詳細に説明されている。

［優れた実施形態の詳細な説明］
本発明の一実施形態によれば、重複領域を有する多数のサブボリュームに関して決定された実質的に同じ特徴値を有するサブボリュームを抽出するためのシステムおよび方法は、重複領域が関心領域（ＲＯＩ）全体、例えば特徴が決定されつつある小結節をカバーすることを想定している。このシステムおよび方法は多数のサブボリュームのサンプリングを同期化するが、サンプリング位置はサブボリュームが重複（オーバラップ）する領域内において実質的に同じである。

図３を参照すると、仮想格子３０３上でオリジナル入力ボリューム３０２をサンプリングすることによって仮想ボリューム３０１が作り出される。仮想格子３０３の原点（０，０，０）は、世界座標系３０４の原点と一致する。更に格子間隔は、サブボリュームの要求された格子間隔と実質的に同じである。要求されたサブボリュームは、仮想ボリューム３０１から抽出される。仮想ボリューム３０１の更なるサンプリングは存在しないので、サブボリュームが重複する領域内のボクセルは実質的に同じである。

仮想ボリュームは、メモリ資源を節約するために、作り出される必要はない。仮想ボリューム格子に位置合わせされた、要求されたサブボリュームは、オリジナルボリュームから直接抽出できる。利用可能な十分なメモリが存在すれば、仮想ボリュームは作り出されることができ、サブボリュームはこの仮想ボリュームから抽出できる。オリジナルボリューム全体を再サンプリングしない仮想ボリュームを使用することによって、本発明の一実施形態による方法は、サブボリュームのための高いサンプリング周波数を考慮している。これは、サブボリュームの処理の効率と精度とを高めることができる。

図４を参照すると、データボリュームと抽出されるサブボリュームを規定する複数のパラメータとを含む入力が供給される（４０１）。これらのパラメータは、世界座標系における中心点（Ｃ_W）と、要求されたサブボリュームサイズ（Ｄ_W）と、サブボリュームのボクセル間隔（ＶＳ_W）とを含み、これらは例えば座標で表される。パラメータは、ユーザによって供給されるか、あるいは予め決定され得る。中心点とサブボリュームサイズとボクセル間隔とを使用して、世界座標系におけるサブボリュームのサンプリング開始点Ｓ_Wが決定される（４０２）。
Ｓ_W＝Ｃ_W−Ｄ_W×ＶＳ_W

サブボリュームのサンプリング開始点は、例えば下記のように表される仮想ボリュームの格子にサブボリュームのサンプリング開始点を位置合わせされるようにシフトされる（４０３）。
Δ＝Ｓ_W ｍｏｄ ＶＳ_W
もし（Δ＞ＶＳ_W／２）であれば、
Ｓ_W＝Ｓ_W＋（ＶＳ_W−Δ）
そうでなければ
Ｓ_W＝Ｓ_W−ＶＳ_W

この決定はサブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせさせる。仮想ボリューム格子の間隔ＶＳ_Wは、要求されたサブボリュームの格子間隔と同じである。Δは、仮想ボリューム格子に対するサブボリュームの位置ずれ量である。サブボリュームのサンプリング開始点のシフト量を減らすために、Δは（ＶＳ_W／２）と比較される。これは、原サンプリング開始点からのシフトを±ＶＳ_W／２に限定する。

サブボリュームのサンプリング開始点は、オリジナルボリュームの（０，０，０）位置を原点として有するボリューム座標に変換される（４０４）。

要求されたボクセル間隔は、ミリメートルからオリジナルボリューム単位に変換される（４０５）。

サブボリュームは、サブボリュームのサンプリング開始点と、ブロック４０４、４０５で決定されたボリューム単位のボクセル間隔とを使用してオリジナルボリュームから抽出される（４０６）。特徴は、抽出された（例えば、位置合わせされた）サブボリュームを使用して、データボリューム上で決定され得る（４０７）。

図５は、２つのサブボリュームが重複する領域内において同じであるサブボリュームに関するサンプリング位置を示す図である。簡単にするために、垂直サンプリングラインだけが示されている。図５には、図２に示す事例に基づいて本発明の一実施形態による結果が示されている。オリジナルボリュームにおける２つのサブボリュームのサンプリング位置２０２、２０３は同じであることが観察できる。その結果、重複領域内のボクセルは、この領域を共有するサブボリューム間で同じである。

位置合わせされたサブボリュームを与えれば、共有されたＲＯＩ１０４について決定された如何なる特徴もサブボリュームに関して同じになる。この特徴不変性は、ＲＯＩの位置に依存しており、また共有されたＲＯＩにないボクセルは特徴の計算と値とに寄与せず、このＲＯＩ全体はサブボリューム間で共有されるという仮定に基づいている。

本発明が種々の形式のハードウエア、ソフトウエア、ファームウエア、専用プロセッサ、あるいはこれらの組合せで実施され得ることは理解されるべきである。一実施形態において、本発明は、プログラム記憶装置に記憶されたアプリケーションプログラムとしてソフトウエアで実施され得る。このアプリケーションプログラムは、任意の適当なアーキテクチャを備えるマシンにアップロードされて、このマシンによって実行され得る。

図６を参照して、本発明の一実施形態によれば、多数のサブボリュームの同期サンプリング方法を実施するためのコンピュータシステム６０１は、とりわけ中央処理装置（ＣＰＵ）６０２と、メモリ６０３と、入力／出力（Ｉ／Ｏ）インタフェース６０４とを含む。コンピュータシステム６０１は、一般にＩ／Ｏインタフェース６０４を介して表示装置６０５とマウス、キーボードのような種々の入力装置６０６とに接続される。表示装置６０５は、仮想ボリュームのビューと、登録された画像とを表示できる。支援回路は、キャッシュメモリ、電源、クロック回路、および通信バスのような回路を含むことができる。メモリ６０３は、ランダムアクセスメモリ（ＲＡＭ）、リードオンリーメモリ（ＲＯＭ）、ディスクドライブ、テープドライブなど、あるいはこれらの組合せを含むことができる。本発明は、メモリ６０３内に記憶され、信号源６０８からの信号を処理するためのＣＰＵ６０２によって実行されるルーチン６０７として実施できる。そのようなものとして、コンピュータシステム６０１は、本発明のルーチン６０７を実行するときに、専用コンピュータシステムになる汎用コンピュータシステムである。

コンピュータプラットホーム６０１はまた、オペレーティングシステムとマイクロ命令コードとを含む。ここに説明された種々の処理と機能は、オペレーティングシステムを介して実行されるマイクロ命令コードの一部またはアプリケーションプログラム（あるいはこれらの組合せ）の一部のいずれかであり得る。更にコンピュータプラットホームには、追加のデータ記憶装置および印刷装置のような種々の他の周辺装置が接続され得る。

更に、図に描かれた、構成システムコンポーネントと方法ステップの一部はソフトウエアで実施され得るので、システムコンポーネント間の（または処理ステップ間の）実際の接続は、本発明がプログラムされる仕方に依存して異なることがある。ここに説明された本発明の教えが与えられれば、従来技術の一般の当業者は、本発明のこれらおよび類似の実施形態または構成を考えることができるであろう。

ここに説明されたシステムおよび方法は、広範囲のアプリケーションに、例えば特徴を決定するためにサブボリュームを処理するアプリケーションに好適である。他の例のアプリケーションは、体の種々の器官に関するコンピュータ支援診断（ＣＡＤ）アプリケーションを含む。

ここに説明されたシステムおよび方法は、２つのアプリケーション、つまり結腸病変の測定と、結腸ＣＡＤにおける特徴検出とに実施されている。

結腸病変の測定において、ユーザは、サブボリュームの抽出のための中心点として使用されるクリックポイントを選択する。このサブボリュームが処理されて、サブボリューム内の病変は自動的に検出される。病変を検出すると、病変のサイズが測定されてユーザに報告される。このアプリケーションに入力されるクリックポイントは、ユーザが病変付近のどこでも自由にクリックできるように変化する。その結果、多数のクリックポイントが多数のサブボリュームを生成する。本発明の一実施形態によれば、多数のサブボリュームは各々、病変に関して同じ測定値を返報する。

同様に結腸ＣＡＤに関しては、結腸内の種々の候補位置についてサブボリュームが抽出される。各サブボリュームに関して特徴が決定される。本発明の一実施形態によれば、決定された特徴値は、候補位置の微妙な変化に伴って変化することはない。これはよりロバストな特徴値を生じる。

種々のモダリティ、すなわちコンピュータ断層撮影（ＣＴ）、磁気共鳴（ＭＲ）、超音波診断（ＵＳ）、単光子放出ＣＴ（ＳＰＥＣＴ）および陽電子放出断層撮影法（ＰＥＴ）を使用して決定されてきたサブボリュームが位置合わせされる。位置合わせされたサブボリュームは、例えば腫瘍検出方法における偽陽性部位を減らして、または特徴検出方法の精度を向上させて、実質的に同じ特徴値を有する特徴を抽出するために実施され得る。本発明の一実施形態によれば、位置合わせされたサブボリュームは、種々のタイプのデータ、例えば異なる分解能を有するデータを分割して処理するために実施され得る。

多数のサブボリュームの同期化されたサンプリングシステムおよび方法に関して実施形態を説明してきたが、前述の教えを参照すれば、当業者によって修正と変形がなされ得ることは留意される。したがって付属の請求項によって規定され本発明の範囲と精神の内にある開示された本発明の特定の実施形態において変更がなされ得ることは理解されるべきである。かくして特許法によって要求される詳細と綿密性とをもって本発明を説明してきたが、特許証によって保護され請求され所望されることは付属の請求項に述べられている。

１つのボリュームから抽出された重複しているサブボリュームを示す図 ２つの重複しているサブボリュームのサンプリング位置の図 本発明の一実施形態による仮想格子上の入力ボリュームの図 本発明の一実施形態による、２つのサブボリュームが重複している領域で同じであるサブボリュームのためのサンプリング位置の図 本発明の一実施形態による方法の流れ図 本発明の一実施形態によるシステムの図

１０１ 境界
１０２ サブボリューム
１０３ サブボリューム
１０４ 関心領域
２０１ ライン
２０２ ライン
２０３ ライン

Claims (19)

1. 特徴抽出のために１つのデータボリュームをサンプリングすることによって抽出された２つ以上のサブボリュームを位置合わせさせる方法ステップを実行するためにコンピュータによって実行可能な命令のプログラムを記憶したコンピュータ読み取り可能なプログラム記憶装置であって、方法ステップは、
   データボリュームを受信するステップと、
   サブボリュームの複数のパラメータを供給するステップと、
   複数のパラメータにしたがって２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップと、
   データボリュームをサンプリングすることによって作り出された仮想ボリュームの格子に対する２つ以上のサブボリュームの位置ずれ量を決定するステップと、
   各サブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせするために、前記位置ずれ量に関係する量によって各サブボリュームのサンプリング開始点をシフトするステップと、
   仮想ボリュームの格子上でサンプリングされた値を用いて、位置合わせされた２つ以上のサブボリュームにおける特徴を決定するステップと
   を含むことを特徴とするコンピュータ読み取り可能なプログラム記憶装置。
2. 複数のパラメータを供給するステップは、
   ２つ以上のサブボリュームの各サブボリュームの中心点を選択するステップと、
   ２つ以上のサブボリュームのサイズと形状とを選択するステップと、
   ２つ以上のサブボリューム内のボクセル間の間隔を選択するステップと
   を含むことを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
3. ボクセル間隔をボリューム単位に変換するステップを更に含むことを特徴とする請求項２に記載のコンピュータ読み取り可能なプログラム記憶装置。
4. 各サブボリュームのサンプリング開始点とボリューム単位に変換されたボクセル間隔とを使用して、特徴を抽出する前に、ボリュームから２つ以上のサブボリュームを抽出するステップを更に含むことを特徴とする請求項３に記載のコンピュータ読み取り可能なプログラム記憶装置。
5. ２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップは、中心点とサブボリュームのサイズとボクセル間隔とにしたがって実行されることを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
6. 位置ずれ量は仮想ボリュームの格子に対して決定され、仮想ボリューム格子の間隔はサブボリュームのボクセル間の間隔であり、また仮想ボリュームの格子の原点は世界座標系の原点と一致することを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
7. シフトされたサンプリング開始点を、各サブボリュームの（０，０，０）位置を原点とするボリューム座標に変換するステップを更に含むことを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
8. ２つ以上のサブボリュームは異なる分解能を有することを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
9. 前記データボリュームは異なるモダリティを表す２つ以上のデータボリュームを含む請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
10. 各サブボリュームのサンプリング開始点をシフトするステップは世界座標系におけるサブボリュームの１つ又は２つの座標軸に対してだけ実行されることを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
11. サブボリュームは２次元であることを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
12. 特徴は２次元サブボリュームについて決定されることを特徴とする請求項１１に記載のコンピュータ読み取り可能なプログラム記憶装置。
13. 各サブボリュームのサンプリング開始点の位置合わせは３次元について決定されることを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
14. １次元以上の次元を選択するステップを更に含み、サンプリング開始点の位置合わせは選択された次元について決定されることを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
15. 各サブボリュームのサンプリング開始点の位置合わせは予め登録された異なるモダリティに亘って適用されることを特徴とする請求項１に記載のコンピュータ読み取り可能なプログラム記憶装置。
16. 視覚化のために１つのデータボリュームをサンプリングすることによって抽出された２つ以上のサブボリュームを位置合わせさせる方法ステップを実行するためにコンピュータによって実行可能な命令のプログラムを記憶したコンピュータ読み取り可能なプログラム記憶装置であって、方法ステップは、
    データボリュームを受信するステップと、
    サブボリュームの複数のパラメータを供給するステップと、
    複数のパラメータにしたがって２つ以上のサブボリュームの各サブボリュームのサンプリング開始点を決定するステップと、
    データボリュームをサンプリングすることによって作り出された仮想ボリュームの格子に対する２つ以上のサブボリュームの位置ずれ量を決定するステップと、
    各サブボリュームのサンプリング開始点を仮想ボリュームの格子に位置合わせするために、前記位置ずれ量に関係する量によって各サブボリュームのサンプリング開始点をシフトするステップと、
    位置合わせされたサンプリング開始点を有するサブボリュームを表示するステップとを含み、
    各サブボリューム内の関心領域は同じボクセル値を有する
    ことを特徴とするコンピュータ読み取り可能なプログラム記憶装置。
17. 複数のパラメータを供給するステップは、
    ２つ以上のサブボリュームの各サブボリュームの中心点を選択するステップと、
    ２つ以上のサブボリュームのサイズと形状とを選択するステップと、
    ２つ以上のサブボリューム内のボクセル間の間隔を選択するステップと
    を含むことを特徴とする請求項１６に記載のコンピュータ読み取り可能なプログラム記憶装置。
18. ボクセル間隔をボリューム単位に変換するステップを更に含むことを特徴とする請求項１７に記載のコンピュータ読み取り可能なプログラム記憶装置。
19. 各サブボリュームのサンプリング開始点とボリューム単位に変換されたボクセル間隔とを使用して、特徴を抽出する前に、ボリュームから２つ以上のサブボリュームを抽出するステップを更に含むことを特徴とする請求項１８に記載のコンピュータ読み取り可能なプログラム記憶装置。

Images