(57)【特許請求の範囲】
【請求項１】
撮像対象物（３）の第１のボリュームデータセット（２１）と、この第１のボリュームデータセット（２１）とは異なる時点で得られた前記撮像対象物（３）の第２のボリュームデータセット（２２）とを用意し、
前記第１のボリュームデータセット（２１）のボリュームエレメントを前記第２のボリュームデータセット（２２）のボリュームエレメントへ写像する、変形ベクトルの第１のセットを決定し、
第２のボリュームデータセット（２２）のボリュームエレメントを第１のボリュームデータセット（２１）のボリュームエレメントへ写像する、変形ベクトルの第２のセットを決定し、
中間ボリュームデータセット（２３ａ，２４ａ）の第１のセットを、中間ボリュームデータセットの第１のセットの各中間ボリュームデータセット（２３ａ，２４ａ）について、第１のボリュームデータセット（２１）のボリュームエレメントに変形ベクトルの第１のセットの減弱された変形ベクトルを掛けることによって作成し、中間ボリュームデータセット（２３ａ，２４ａ）の１つに属する減弱された変形ベクトルが、変形係数Ａiを用いて変形ベクトルの第１のセットの変形ベクトルの重み付けを行うことによって算出され、中間ボリュームデータセットの第１のセットがｉ個の中間ボリュームデータセット（２３ａ，２４ａ）を有し、０＜Ａi＜１であり、
中間ボリュームデータセット（２３ｂ、２４ｂ）の第２のセットを、中間ボリュームデータセットの第２のセットの各中間ボリュームデータセット（２３ｂ、２４ｂ）について、第２のボリュームデータセット（２２）のボリュームエレメントに変形ベクトルの第２のセットの減弱された変形ベクトルを掛けることによって作成し、中間ボリュームデータセット（２３ｂ、２４ｂ）の１つに属する減弱された変形ベクトルが、変形係数Ｂiを用いて変形ベクトルの第２のセットの変形ベクトルの重み付けを行うことによって算出され、中間ボリュームデータセットの第２のセットがｉ個の中間ボリュームデータセットを有し、０＜Ｂi＜１であり、
画像変形ボリュームデータセット（２３，２４）を、中間ボリュームデータセットの第１のセットの中間ボリュームデータセット（２３ａ，２４ａ）のグレー値に、（１−変形係数Ａi）に相当する係数Ｋiを掛け、中間ボリュームデータセットの第２のセットの中間ボリュームデータセット（２３ｂ、２４ｂ）のグレー値に、（１−係数Ｋi）に相当する係数Ｋi´を掛けることによって作成し、その場合、対応する中間ボリュームデータセットの変形係数Ａi，Ｂiの合計が１に等しく、対応する変形係数ＡiないしはＢiを掛けられた中間ボリュームデータセットのグレー値が加算され、画像変形ボリュームデータセット（２３，２４）のセットの個々の画像変形ボリュームデータセットが増大する変形係数Ａiに従って分類され、
画像変形ボリュームデータセット（２３，２４）をシーケンスとして（２０）表示する、
ことを特徴とするシーケンスの作成方法。
【請求項２】
変形ベクトルの第１のセットが、
第１のボリュームデータセット（２１）および第２のボリュームデータセット（２２）の解剖学的に対応する各１個のボリュームエレメントを含む個々のボリュームエレメント対を手動で決定し、その場合、第１のボリュームデータセット（２１）のそれぞれのボリュームエレメントが第２のボリュームデータセット（２２）の対応するボリュームエレメントへ写像され、
手動で求められたボリュームエレメント対に属する変形ベクトルを決定し、
変形ベクトルの第１のセットの残りの変形ベクトルを、手動で求められたボリュームエレメント対に基づいて決定される変形ベクトルを適正に補間することによって決定する、
ことによって決定されることを特徴とする請求項１記載の方法。
【請求項３】
変形ベクトルの第１のセットが、
パターン認識手段を用いて第１のボリュームデータセット（２１）および第２のボリュームデータセット（２２）を分析し、
この分析に基づいて、第１のボリュームデータセット（２１）および第２のボリュームデータセット（２２）の解剖学的に対応する各１個のボリュームエレメントを含むボリュームエレメント対を求め、その場合、第１のボリュームデータセット（２１）のそれぞれのボリュームエレメントが第２のボリュームデータセット（２２）の対応するボリュームエレメントへ写像され、
分析に基づいて変形ベクトルの第１のセットの変形ベクトルを決定する、
ことによって決定されることを特徴とする請求項１又は２記載の方法。
【請求項４】
変形ベクトルの第１のセットが、
第１のボリュームデータセット上に置かれる第１の格子を決定し、
第１の格子の個々の格子点を、第１のボリュームデータセットが第２のボリュームデータセット上にできるだけ良好に写像されるまで移動させ、
この移動に基づいて変形ベクトルの第１のセットの変形ベクトルを決定する、
ことによって決定されることを特徴とする請求項１乃至３の１つに記載の方法。
【請求項５】
変形ベクトルの第２のセットが、
第１のボリュームデータセット（２１）および第２のボリュームデータセット（２２）の解剖学的に対応する各１個のボリュームエレメントを含む個々のボリュームエレメント対を手動で決定し、その場合、第２のボリュームデータセット（２２）のそれぞれのボリュームエレメントが第１のボリュームデータセット（２１）の対応するボリュームエレメントへ写像され、
手動で求められたボリュームエレメント対に属する変形ベクトルを決定し、
手動で求められたボリュームエレメント対に基づいて決定される変形ベクトルを適正に補間することにより変形ベクトルの第２のセットの残りの変形ベクトルを決定する、
ことによって決定されることを特徴とする請求項１乃至４の１つに記載の方法。
【請求項６】
変形ベクトルの第２のセットが、
パターン認識手段を用いて第１のボリュームデータセットおよび第２のボリュームデータセットを分析し、
この分析に基づいて、第１のボリュームデータセット（２１）および第２のボリュームデータセット（２２）の解剖学的に対応する各１個のボリュームエレメントを含むボリュームエレメント対を求め、その場合、第２のボリュームデータセット（２２）のそれぞれのボリュームエレメントが第１のボリュームデータセット（２１）の対応するボリュームエレメントへ写像され、
分析に基づいて変形ベクトルの第２のセットの変形ベクトルを決定する、
ことによって決定されることを特徴とする請求項１乃至５の１つに記載の方法。
【請求項７】
変形ベクトルの第２のセットが、
第２のボリュームデータセット上に置かれる第２の格子を決定し、
第２の格子の個々の格子点を、第２のボリュームデータセットが第１のボリュームデータセット上にできるだけ良好に写像されるまで移動させ、
この移動に基づいて変形ベクトルの第２のセットの変形ベクトルを決定する、
ことによって決定されることを特徴とする請求項１乃至６の１つに記載の方法。
【請求項８】
撮像対象物が生物（３）であることを特徴とする請求項１乃至７の１つに記載の方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、２つのボリュームデータセットに基づいてシーケンスを作成する方法に関する。
【０００２】
【従来の技術】
撮像対象物の２つの２次元画像から、画像変形、例えば"モーフィング"により１つの画像シーケンスを作成することが公知である（例えば非特許文献１参照）。画像シーケンスは画像の１つである出発画像で開始し、別の画像である目標画像で終了する。画像シーケンスのために必要な中間画像は補間により出発画像および目標画像から作成される。まず、出発画像および目標画像の互いに対応する一対の点が決定される。これらの対応する一対の点により、目標画像に対する出発画像の変形および出発画像に対する目標画像の変形が算出される。中間画像は減弱された変形によって算出される。出発画像ないし目標画像の画像点は、変形により決定された区間の一部だけ移動させられる。これにより、出発画像から出発し、変形が徐々に増強された画像シリーズおよび変形が徐々に減弱されて目標画像へ至る画像シリーズが得られる。続いてこれら２つのシリーズは画像変形され、それにより出発画像から開始し、目標画像で終了する２次元画像シリーズが形成される。
【０００３】
医学的フォローアップ研究においては、例えば患者の組織構造の大きさ、位置および形状変化を識別するために例えば異なった時点における患者の状況の撮像が行われる。フォローアップ研究の例は患者の腫瘍治療の観察である。腫瘍治療前および治療中に、例えば患者を治療する医師はコンピュータトモグラフを用いて腫瘍の断層画像を撮像し、それらを比較する。断層画像から立体的印象が生じないので、断層画像の評価が困難である。
【０００４】
心臓カテーテルの挿入前に、心臓カテーテルが挿入される心臓のＣＴ画像を作成することができる。ＣＴ画像は、心臓カテーテルの位置をコントロールする、あるいは心臓カテーテルを的確に心臓の解剖学的構造に対して相対的に制御するために、心臓カテーテルの挿入中にＥＣＧでトリガされた心臓のＸ線画像を用いて画像変形される。心筋の収縮ないし弛緩によりカテーテルの位置は心臓の解剖学的構造に対して相対的に著しく変化する可能性があるので、手術前に撮像されたＣＴ画像に関するカテーテルの位置を何時の時点でも視覚化するために、ＣＴ画像の離散的時点とは異なる時点でも画像変形を表示することが望ましい。これを実現するには、侵襲の前に心臓リズムのほぼ全期間をカバーするＣＴ画像のシーケンスを取得する必要がある。これは患者に対するＸ線負荷が大きすぎるので実現不可能である。
【０００５】
【非特許文献１】
Ｄｅｔｌｅｆ  Ｒｕｐｒｅｃｈｔ著「図形データ処理における道具としての幾何学的変形」ドルトムント大学学位請求論文、１９９４年発行、ｐ１４６〜１４９（Ｓｈａｋｅｒ出版社、アーヘン、１９９５年公刊）
【０００６】
【発明が解決しようとする課題】
このため、本発明の課題は、２つの異なった時点で１つの撮像対象物から作成された２つのボリュームデータセットからボリュームデータセットのシーケンスを作成する方法を提供することにある。
【０００７】
【課題を解決するための手段】
本発明の課題は、本発明によれば、
撮像対象物の第１のボリュームデータセットと、この第１のボリュームデータセットとは異なる時点で得られた前記撮像対象物の第２のボリュームデータセットとを用意し、
撮像対象物の第１のボリュームデータセットのボリュームエレメントを撮像対象物の第２のボリュームデータセットのボリュームエレメントへ写像する、変形ベクトルの第１のセットを決定し、
第２のボリュームデータセットのボリュームエレメントを第１のボリュームデータセットのボリュームエレメントへ写像する、変形ベクトルの第２のセットを決定し、
ｉ個の中間ボリュームデータセットを有する中間ボリュームデータセットの第１のセットを、変形ベクトルの第１のセットの変形ベクトルのｉ番目の減弱されたセットを有する第１のボリュームデータセットが中間ボリュームデータセットの第１のセットのｉ番目の中間ボリュームデータセットへ写像されることによって作成し、その場合、変形ベクトルの第１のセットの減弱された変形ベクトルのｉ番目のセットが変形ベクトルの第１のセットの変形ベクトルに変形係数Ａiを掛けることによって算出され、０＜Ａi＜１であり、
ｉ個の中間ボリュームデータセットを有する中間ボリュームデータセットの第２のセットを、変形ベクトルの第２のセットの変形ベクトルのｉ番目の減弱されたセットを有する第２のボリュームデータセットが中間ボリュームデータセットの第２のセットのｉ番目の中間ボリュームデータセットへ写像されることによって作成し、その場合、変形ベクトルの第２のセットの減弱された変形ベクトルのｉ番目のセットが変形ベクトルの第２のセットの変形ベクトルに変形係数Ｂiを掛けることによって算出され、０＜Ｂi＜１であり、
画像変形ボリュームデータセットのセットを、中間ボリュームデータセットの第１のセットの中間ボリュームデータセットのグレー値に、（１−変形係数Ａi）に相当する係数Ｋiを掛け、中間ボリュームデータセットの第２のセットの中間ボリュームデータセットのグレー値に、（１−係数Ｋi）に相当する係数Ｋi´を掛けることによって作成し、その場合、対応する中間ボリュームデータセットの変形係数Ａi，Ｂiの合計が１に等しく、対応する変形係数ＡiないしはＢiを掛けられた中間ボリュームデータセットのグレー値が加算され、画像変形ボリュームデータセットのセットの個々の画像変形ボリュームデータセットが増大する変形係数Ａiに基づいて分類され、
画像変形ボリュームデータセットをシーケンスとして表示する、
ことによって解決される。
【０００８】
本発明による方法によれば、先ず、変形ベクトルの第１のセットおよび第２のセットが、第１のボリュームデータセットのボリュームエレメントが第２のボリュームデータセットのボリュームエレメントへ写像されることによって決定される。これらの変形ベクトルセットを用いて、どのようにして２つのボリュームデータセットにより写像された画像情報が互いに相対的に移動されているかが示される。次に、中間ボリュームデータセットの２つのセットが作成される。変形係数Ａi，Ｂiの値に依存して、個々の中間ボリュームセットの画像情報が相応のボリュームデータセットと関連して強く移動させられる。続いて、画像変形ボリュームデータセットのセットが作成される。このセットは増大する変形係数Ａiに基づいて分類されているので、第１のボリュームデータセットで開始し、その後増大するＡiに基づいて分類されている画像変形ボリュームデータセットを有し、第２のボリュームデータセットで終了するボリュームデータセットのシーケンスが生じる。
【０００９】
第１のボリュームデータセットおよび第２のボリュームデータセットの各１個のボリュームエレメントを含む個々のボリュームエレメント対を手動で決定し、その場合、第１のボリュームデータセットのそれぞれのボリュームエレメントが第２のボリュームデータセットの相応のボリュームエレメントへ写像され、
手動で求められたボリュームエレメント対に属する変形ベクトルを決定し、
変形ベクトルの第１のセットの残りの変形ベクトルを、手動で求められたボリュームエレメント対に基づいて決定される変形ベクトルを適正に補間することによって決定する、
ことによって決定される。
【００１０】
２つのボリュームデータセットのボリュームエレメント対は、例えば２つのボリュームデータセットに含まれている画像情報の対応する輪郭、表面もしくは目立った点である。
【００１１】
本発明の望ましい変形によれば、個々のボリュームエレメントを手動で決定する代わりに、変形ベクトルの第１のセットは、
パターン認識手段を用いて第１のボリュームデータセットおよび第２のボリュームデータセットを分析し、
この分析に基づいて、第１のボリュームデータセットおよび第２のボリュームデータセットの各１個のボリュームエレメントを含むボリュームエレメント対を求め、その場合、第１のボリュームデータセットのそれぞれのボリュームエレメントが第２のボリュームデータセットの相応のボリュームエレメントへ写像され、
分析に基づいて変形ベクトルの第１のセットの変形ベクトルを決定する、
ことによって決定される。
【００１２】
パターン認識手段を用いた分析に基づいて、変形ベクトルの第１のセットを自動的に求めることができる。また、個々の変形ベクトルの手動で求めることと、２つのボリュームデータセットの分析とを混ぜ合わせることも可能である。
【００１３】
本発明の別の望ましい実施態様によれば、変形ベクトルの第１のセットは、
第１のボリュームデータセット上に置かれる第１の格子を決定し、
第１の格子の個々の格子点を、第１のボリュームデータセットが第２のボリュームデータセット上にできるだけ良好に写像されるまで移動させ、
この移動に基づいて変形ベクトルの第１のセットの変形ベクトルを決定する、
ことによって決定される。
【００１４】
格子の決定に基づいて変形ベクトルの第１のセットの変形ベクトルがベクトル界において自動的に求められる。
【００１５】
本発明の実施態様によれば、変形ベクトルの第２のセットは変形ベクトルの第１のセットと類似して決定される。
【００１６】
本発明による方法は、本発明による実施態様により撮像対象物が生物である場合には特に医療技術において用いられる。例えば、本発明の望ましい実施態様に基づいて考えられているように、第１のボリュームデータセットおよび第２のボリュームデータセットは患者の心拍動の２つの異なった時相中に取得可能である。冒頭に記載したように心臓カテーテルの設置前に心臓カテーテルが挿入される心臓のＣＴ画像が作成可能である。このＣＴ画像は心臓カテーテルの挿入中にＥＣＧでトリガされた心臓のＸ線画像を用いて画像変形される。心臓リズムのほぼ全期間をカバーするＣＴ画像のシーケンスの取得は患者にとってＸ線照射負荷が大きすぎるため実現不可能であるので、従来の技術では侵襲前に心拍動リズムの決められた時点に１回のＣＴ画像しか作成されない。心筋の収縮ないし弛緩によりカテーテルの位置が心臓の解剖学的構造に対して相対的に著しく変化するので、画像変形は通常、ＣＴ画像の離散的時点でのみ表示される。しかし、本発明による方法は全心拍動リズム中の心拍動の連続的なシーケンスを僅かなＣＴ画像に基づいて獲得することを可能にする。侵襲前に、例えば心臓の３つのボリュームデータセットが心拍動リズムの異なった時点で例えばコンピュータトモグラフにより取得される。通常、均等に拍動する心臓の連続的なシーケンスを獲得するために画像変形ボリュームデータセットの３つのセットが作成される。例えば画像変形ボリュームデータセットのセットは、心拍動リズムの開始時に作成された第１のボリュームデータセットから心拍動リズムに関して後で取得された第２のボリュームデータセットへの連続的な移行を可能にし、画像変形ボリュームデータセットの別のセットは、第２のボリュームデータセットから第２のボリュームデータセットに関して後で取得された第３のボリュームデータセットへの移行を可能にする。最後に、画像変形ボリュームデータセットの別のセットは第３のボリュームデータセットから第１のボリュームデータセットへの移行を可能にする。このようにして、心拍動は取得された僅かなボリュームデータセットによりシーケンスとして視覚的に表示できる。
【００１７】
本発明の別の変形によれば、撮像対象物が腫瘍である。従って、医師は腫瘍の継時的増殖を２回もしくはそれより少ない回数で取得されたボリュームデータセットに基づいてより良好に観察することができる。この結果、実施された治療の成果もはるかに容易に判定可能となる。腫瘍の変化の視覚化は断層画像に限定されない。というのは、より多くの情報が含まれている３次元ボリュームが観察されるからである。表面陰影表示（Ｓｕｒｆａｃｅ  Ｓｈａｄｅｄ  Ｄｉｓｐｌａｙ（ＳＳＤ））を用いて、区分された表面の写像を視覚化でき、バーチャルレンダリング技法（Ｖｉｒｔｕａｌ  Ｒｅｎｄｅｒｉｎｇ  Ｔｅｃｈｎｉｑｕｅ（ＶＲＴ））を用いて身体の透視画像も視覚化できる。この立体表示により、医師は腫瘍の変化をより容易に発見することが可能となる。
【００１８】
本発明の別の変形によれば、第１のボリュームデータセットおよび第２のボリュームデータセットは同一の型の撮像装置もしくは同一の撮像装置を用いて取得される。
【００１９】
本発明による方法は特に医療技術に適用可能であるので、別の変形によれば撮像装置は医療技術装置である。
【００２０】
【発明の実施の形態】
実施例は添付されている概略図に例として示されている。
【００２１】
図１は、縁部Ｘ線が一点鎖線で示されている角錐状のＸ線束２を放射するＸ線源１を有するコンピュータトモグラフＣＴの概略図を示している。Ｘ線束２は患者３を貫通し、Ｘ線検出器４に当たっている。この実施例ではＸ線源１およびＸ線検出器４は環状のガントリ５に互いに向かい合って配置されている。ガントリ５は環状のガントリ５の中心点を通って延びるシステム軸線６に対して図１には示されていない支持装置に回転可能に支承されている（矢印ａを参照）。
【００２２】
この実施例では患者３は、図１にはやはり示されていない支持装置によりシステム軸線６に沿って移動可能（矢印ｂを参照）に支承されているＸ線透過性テーブル７上に存在する。
【００２３】
従って、Ｘ線源１およびＸ線検出器４は、システム軸線６を中心にして回転可能であり、システム軸線６に沿って患者３に対して相対的に移動可能である測定装置を形成しており、このため患者３はシステム軸線６に対して種々の投射角および種々の位置で投射可能である。この測定装置はスパイラル走査中にシステム軸線６を中心にして回転し、テーブル７は矢印ｂ方向へ連続的に移動する、すなわちＸ線源１およびＸ線検出器４を含む測定装置は患者３に対して相対的にスパイラル軌道Ｃ上で連続的に患者３の関心領域が完全に捕捉されるまで移動する。この際に発生するＸ線検出器４の出力信号から、データ検出装置９がコンピュータ１１へ供給される測定値を形成する。コンピュータ１１上では専門化には周知のコンピュータプログラムがランし、このプログラムが測定値から患者３の関心領域のボリュームデータセットを算出する。このボリュームデータセットはコンピュータ１１のメモリ１ｂ内に記憶され、ボリュームデータセットに属する画像がコンピュータ１１のモニタ１１ａ上に再現される。コンピュータのマウス１１により、ボリュームデータセットに属する画像は例えば回転され、その結果種々の注視方向から観察される。この実施例ではデータ検出装置９は、図示されていない例えばスリップリング装置もしくは無線の伝送区間を含む電気配線８によりＸ線検出器４に接続されている。コンピュータ１１とデータ検出装置９は電気配線１０に接続されている。
【００２４】
この実施例では患者は腫瘍を有している。腫瘍の継時的変化を観察するために、患者３を治療する医師（図１には示されていない）が患者３の腫瘍を含む２つのボリュームデータセットを作成する。医師は異なった２日間で２つのボリュームデータセットを取得し、それらをコンピュータ１１のメモリ１１ｂに記憶させる。この実施例では医師は第１のボリュームデータセットの作成から２週間後に第２のボリュームデータセットを取得する。
【００２５】
腫瘍の変化を識別するために、医師は２つのボリュームデータセットに属する画像をモニタ１１ｂにより並列に表示することができる。ただし、医師は腫瘍の変化についてより良好な印象を得るために２つのボリュームデータセットから、腫瘍の段階的な継時的変化を識別できるシーケンスを作成して得ることができる。このため、医師は、モニタ１１aに表示されコンピュータ１１のマウス１１ｃにより移動可能であるマーク（図１には示されていない）を用いて、複数のボリュームエレメント対をマーキングする。ボリュームエレメント対はそれぞれ第１のボリュームデータセットの１個のボリュームエレメントおよび第２のボリュームデータセットの１個のボリュームエレメントを含む。ボリュームエレメント対は、どのようにして第１の第１のボリュームデータセットのボリュームエレメントが第２のボリュームデータセットの対応するボリュームエレメントへ写像されるかを示す。この実施例では、マーキングされたボリュームエレメントは変化した腫瘍の表面に割付けられている。
【００２６】
ボリュームエレメント対に基づいて、コンピュータ１１上でランするコンピュータプログラムが変形ベクトルの第１のセットを算出し、このセットが第１のボリュームデータセットの画像内容を第２のボリュームデータセットの画像内容へ変形する。コンピュータ１１のコンピュータプログラムはボリュームエレメント対から直接に変形ベクトルの第１のセットの変形ベクトルの一部を算出する。ただし、医師は僅かなボリュームエレメント対のみをマーキングするにすぎないので、コンピュータプログラムは、変形ベクトルの第１のセットの、変形に必要な別の変形ベクトルを、ボリュームエレメント対に基づいて算出された変形ベクトルを補間することによって算出する。
【００２７】
続いて、コンピュータ１１のコンピュータプログラムは、第２のボリュームデータセットの画像内容を第１のボリュームデータセットの画像内容に変形する、変形ベクトルの第２のセットを算出する。コンピュータプログラムはボリュームエレメント対から直接に変形ベクトルの第２のセットの変形ベクトルの一部を算出する。コンピュータプログラムはボリュームエレメント対に基づいて算出された変形ベクトルを補間することによって、変形ベクトルの第２のセットの、変形に必要な別な変形ベクトルを算出する。
【００２８】
次に、コンピュータ１１上でランするコンピュータプログラムはこの実施例ではそれぞれ２個の中間ボリュームデータセットを有する中間ボリュームデータセットの２つのセットを算出する。
【００２９】
中間ボリュームデータセットの第１のセットの２つの各中間ボリュームデータセットは、変形ベクトルの第１のセットの変形ベクトルの各１個の減弱されたセットを有する第１のボリュームデータセットが中間ボリュームデータセットの第１のセットの対応する中間ボリュームデータセットへ写像されることによって算出される。変形ベクトルの第１のセットの減弱された変形ベクトルは、変形ベクトルの第１のセットの変形ベクトルに変形係数Ａi（ｉ＝１；２および０＜Ａi＜１）を掛けることによって算出される。この実施例では中間ボリュームデータセットの第１のセットの第１の中間ボリュームデータセット用の変形係数Ａ1は０．３３であり、中間ボリュームデータセットの第１のセットの第２の中間ボリュームデータセット用の変形係数Ａ2は０．６７である。
【００３０】
中間ボリュームデータセットの第２のセットの２つの中間ボリュームデータセットも同様に算出される。すなわち、第２のボリュームデータセットは、変形ベクトルの第２のセットの変形ベクトルの各１個の減弱されたセットを用いて、中間ボリュームデータセットの第２のセットの対応する中間ボリュームデータセットへ写像される。変形ベクトルの第２のセットの減弱された変形ベクトルは、変形ベクトルの第２のセットの変形ベクトルに変形係数Ｂi（ｉ＝１；２および０＜Ｂi＜１）を掛けることによって算出される。この実施例では中間ボリュームデータセットの第２のセットの第１の中間ボリュームデータセット用の変形係数Ｂ1は０．６７であり、中間ボリュームデータセットの第２のセットの第２の中間ボリュームデータセット用の変形係数Ｂ2は０．３３である。
【００３１】
続いて、コンピュータ１１上でランするコンピュータプログラムがこの実施例では２つの画像変形ボリュームデータセットを含む画像変形ボリュームデータセットのセットを算出する。第１の画像変形ボリュームデータセットは、中間ボリュームデータセットの第１のセットの第１の中間ボリュームデータセットのグレー値に（１−変形係数Ａ1）すなわち０．６７を掛け、中間ボリュームデータセットの第２のセットの第１の中間ボリュームデータセットのグレー値に（１−変形係数Ｂ1）すなわち０．３３を掛けることによって算出される。続いて、対応する中間ボリュームデータセットの個々のボリュームエレメントのグレー値が加算され、それにより第１の画像変形ボリュームデータセットが生成する。
【００３２】
第２の画像変形ボリュームデータセットは、中間ボリュームデータセットの第１のセットの第２の中間ボリュームデータセットに変形係数Ａ2＝０．６７を掛け、第２の中間ボリュームデータセットの第２のボリュームデータセットに変形係数Ｂ2＝０．３３を掛けることによって算出される。続いて、対応する中間ボリュームデータセットのグレー値が加算される。
【００３３】
次に、コンピュータ１１がモニタ１１ａ上に、第１のボリュームデータセットで開始し、第２のボリュームデータセットで終了するシーケンスを表示する。第１のボリュームデータセットと第２のボリュームデータセットの間に番号の増大する順に第１の画像変形ボリュームデータセットおよび第２の画像変形ボリュームデータセットが表示される。このシーケンスは時間的に次々にモニタ１１ａ上に表示され、その結果医師は容易に腫瘍の変化を観察することができる。
【００３４】
図２は、ボリューム２１で開始しボリューム２２で終了するシーケンス２０を示している。ボリューム２１に属するボリュームデータセットは、例えば上述した腫瘍である撮像対象物に関して時間的に、ボリューム２２に属するボリュームデータセットの前に撮像されている。このシーケンスはさらにボリューム２３，２４も含む。ボリューム２３，２４に属する画像変形ボリュームデータセットはそれぞれ２個の中間ボリュームデータセットから作成される。中間ボリュームデータセットに属するボリューム２３ａ〜２４ｂは図３に示されている。
【００３５】
ボリューム２３，２４に属する画像変形ボリュームデータセットを作成するために、先ずボリューム２１，２２に属するボリュームデータセットのボリュームエレメント対が求められ、その結果からボリューム２１の画像内容をボリューム２２の画像内容へ写像する最適の変形ベクトルのセットを獲得することができる。変形ベクトルの第２のセットも同様に求められる。すなわち、変形ベクトルの第２のセットはボリューム２２の画像内容をボリューム２１の画像内容へ写像する。
【００３６】
次に、中間ボリュームデータセットの２つのセットが求められる。この実施例では中間ボリュームデータセットの２つのセットのいずれもが２つの中間ボリュームデータセットを含む。
【００３７】
中間ボリュームデータセットの第１のセットの中間ボリュームデータセットは、ボリューム２１に属するボリュームデータセットが、ボリューム２１の画像内容をボリューム２２へ写像する変形ベクトルの各１個の減弱されたセットを用いて写像されることによって算出される。この実施例では第１の中間ボリュームデータセット用の変形係数Ａ1は０．３３であり、第２の中間ボリュームデータセット用の変形係数Ａ2は０．６７である。これらの中間ボリュームデータセットに属するボリュームは符号２３ａ，２４ａが付され、ボリューム２３ａには変形係数Ａ1＝０．３３が属し、ボリューム２４ａには変形係数Ａ2＝０．６７が属している。
【００３８】
中間ボリュームデータセットの第２のセットの中間ボリュームデータセットは、ボリューム２２に属するボリュームデータセットが、ボリューム２２の画像内容をボリューム２１へ写像する変形ベクトルの各１個の減弱されたセットを用いて写像されることによって算出される。この実施例では第１の中間ボリュームデータセット用の変形係数Ｂ1は０．６７であり、第２の中間ボリュームデータセット用の変形係数Ｂ2は０．３３である。これらの中間ボリュームデータセットに属するボリュームは符号２３ｂ，２４ｂが付され、ボリューム２３ｂには変形係数Ｂ1＝０．６７が属し、ボリューム２４ｂには変形係数Ｂ2＝０．３３が属している。
【００３９】
ボリューム２３に属する画像変形ボリュームデータセットは、ボリューム２３ａに属する中間ボリュームデータセットのグレー値に（１−変形係数Ａ1）に一致する０．６７を掛けられ、ボリューム２３ｂに属する中間ボリュームデータセットのグレー値に（１−変形係数Ｂ1）に一致する０．３３を掛けられることによって算出される。続いて、対応する中間ボリュームデータセットの個々のボリュームエレメントのグレー値が加算される。
【００４０】
ボリューム２４に属する画像変形ボリュームデータセットは、ボリューム２４ａに属する中間ボリュームデータセットのグレー値に（１−変形係数Ａ2）に一致する０．３３を掛けられ、ボリューム２４ｂに属する中間ボリュームデータセットのグレー値に（１−変形係数Ｂ2）に一致する０．６７を掛けられることによって算出される。続いて、対応する中間ボリュームデータセットの個々のボリュームエレメントのグレー値が加算される。
【００４１】
この実施例では第１のボリュームデータセットおよび第２のボリュームデータセットはコンピュータトモグラフＣＴにより作成されている。本発明による方法の場合、それによりボリュームデータセットが作成されるその他の装置もしくは技術装置（必ずしも医療技術装置ないし技術装置に限定されない）にも適用可能である。第１のボリュームデータセットおよび第２のボリュームデータセットもやはり必ずしも同一の撮像装置により作成される必要はない。
【００４２】
ボリュームエレメント対もやはり必ずしも手動で求められる必要はない。これらのボリュームエレメント対はパターン認識手段によっても発見できる。また、ボリュームエレメント対を手動で求めること及び自動で求めることの組み合わせも考慮できる。
【００４３】
上述した実施例は単に１例を示すにすぎない。
【図面の簡単な説明】
【図１】コンピュータトモグラフを示す概略図
【図２】シーケンスを示す概略図
【図３】図２に示されているシーケンスに属する変形された画像内容を示す概略図
【符号の説明】
１    Ｘ線源
２    Ｘ線束
３    患者
４    Ｘ線検出器
５    ガントリ
６    システム軸線
７    テーブル
８    電気配線
９    データ検出装置
１０  電気配線
１１  コンピュータ
１１ａ  モニタ
１１ｂ  メモリ
１１ｃ  マウス
２０  シーケンス
２１〜２４  ボリューム
Ａi，Ｂi 変形ベクトル
ａ，ｂ    矢印
ｃ      スパイラル軌道
Ｋi，Ｋi´  係数
【図１】
0004424916.tif 000002
【図２】
0004424916.tif 000003
【図３】
0004424916.tif 000004