(57)【特許請求の範囲】
【請求項１】
同一被写体についての２つの画像の位置合わせを行う位置合わせ手段を備えた画像の位置合わせ装置において、
前記２つの画像が３次元透過画像と２次元透過画像であり、
前記位置合わせ手段が、前記３次元透過画像から前記２次元透過画像と相関が高い２次元投影画像を作成し、該作成された２次元投影画像を前記２次元透過画像と位置合わせするものであって、前記２次元投影画像を作成する際に、骨部と軟部との濃度値の比が前記２次元投影画像と前記２次元透過画像とにおいて同等になるように前記２次元投影画像に骨部と軟部との透過率を変える濃度補正を加えるものであることを特徴とする画像の位置合わせ装置。
【請求項２】
同一被写体についての２つの画像の位置合わせを行う位置合わせ手段と、該位置合わせ手段により位置合わせがなされた２つの画像を用いて所定の画像処理を行う画像処理手段とを備えた画像処理装置において、
前記同一被写体についての２つの画像が３次元透過画像と２次元透過画像であり、前記位置合わせ手段が、前記３次元透過画像から前記２次元透過画像と相関が高い２次元投影画像を作成し、該作成された２次元投影画像を前記２次元透過画像と位置合わせするものであって、前記２次元投影画像を作成する際に、骨部と軟部との濃度値の比が前記２次元投影画像と前記２次元透過画像とにおいて同等になるように該２次元投影画像に骨部と軟部との透過率を変える濃度補正を加えるものであることを特徴とする画像処理装置。
【請求項３】
前記所定の画像処理が、サブトラクション処理であることを特徴とする請求項２記載の画像処理装置。
【請求項４】
前記位置合わせ手段が、前記３次元透過画像および前記２次元透過画像を縮小して縮小画像データを作成する画像縮小手段を有するものであり、
前記位置合わせ手段において前記画像縮小手段により作成された前記３次元透過画像を用いて前記２次元投影画像が作成されたときに、該２次元投影画像と前記２次元透過画像とに対しそれぞれテンプレート領域と探索領域とを設定して対応する領域ごとにシフト量を求め、ワーピングを用いて局所的に位置合わせを行うことにより、局所的位置合わせを行う局所的位置合わせ手段をさらに備え、
前記画像処理手段が、局所的位置合わせがなされた前記２次元投影画像と前記２次元透過画像とを入力して画素を対応させた画像間演算を行い、サブトラクション画像を生成するものであることを特徴とする請求項３記載の画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は画像の位置合わせ装置および画像処理装置に関し、詳細には、同一被写体について得られた２つの画像の位置合わせを行う画像の位置合わせ装置、および、同一被写体について得られた２つの画像の位置合わせを行い、該位置合わせがなされた２つの画像を用いて画像処理を行う画像処理装置に関するものである。
【０００２】
【従来の技術】
従来より、同一被写体についての２以上の画像を比較読影して、両画像間の差異を調べ、その差異に基づいて被写体の検査等を行うことが、種々の分野において行われている。
【０００３】
例えば、工業製品の製造分野においては、ある製品について新品の状態のときに撮影された画像と当該製品の耐久試験後に撮影された画像とを比較読影して、両者の差異の大きな部位に注目することにより、製品の耐久性を向上させるべき部位を検討することが行われており、また医療分野においては、ある患者の疾患部位について時系列的に撮影された複数枚の放射線画像を医師が比較読影することにより、当該疾患の進行状況や治癒状況を把握して治療方針を検討することが行われている。
【０００４】
このように、２以上の画像を比較読影することが日常的に各種の分野で行われているが、比較読影の際には、２以上の画像を画像表示装置やプリンタ等により可視画像として出力し、これらの画像を並べて参照するのが一般的である。この際、読影者にとって最も関心があるのはこれらの画像間の差異であるが、差異を発見する能力は読影者の経験や技量により異なるものであり、さらに、その差異が小さければ小さいほど発見は困難である。このため、読影者の技量等に因ることなく、比較読影の性能を向上させることが求められている。
【０００５】
そこで、比較読影の対象とされる２つの画像間で解剖学的特徴位置を対応させた減算（サブトラクション）処理をはじめとする画像間演算を行い、画像間の差異を抽出・強調したサブトラクション画像を得る技術が既に提案されている（特願平11−342900号公報等）。このように画像間の差異のみを抽出・強調することにより、読影者に対して画像間の差異を確実に認識させることができるため、進行または治癒する病変部の見落とし等を防止することができると考えられる。
【０００６】
また、この画像間演算の際に、各画像中に表われた構造物の位置（解剖学的特徴位置）を２つの画像間で対応させる位置合わせを行う必要がある。この位置合わせの技術としては、例えば、２つの画像間で平行移動、回転または拡大・縮小という全体的な位置合わせ（例えばアフィン変換等を用いた線形的な位置合わせ）を行ない、さらに、この全体的な位置合わせがなされた２つの画像のうち一方の画像に多数の小領域である関心領域（テンプレート領域）を設定するとともに、他方の画像に各テンプレート領域にそれぞれ対応する、各テンプレート領域よりも大きい領域の探索領域を設定し、テンプレート領域と対応する探索領域との各組において、テンプレート領域と画像が略一致する、探索領域内の部分領域（対応テンプレート領域）を求め、一方の画像における各テンプレート領域と、他方の画像における各対応テンプレート領域との対応位置関係に基づいて、一方の画像における各テンプレート領域を他方の画像における各対応テンプレート領域に一致させるためのシフト量を求め、求められたシフト量に基づいて、上記全体的な位置合わせがなされた２つの画像に対し、カーブフィッティング（２次元ｎ次多項式，ｎ≧２）による非線形歪変換（ワーピング）を用いた局所的な位置合わせを行う、２段階の位置合わせ技術（特開平7-37074号公報、米国特許No.5,982,915）等が知られている。この２段階の位置合わせ技術によれば、２つの画像を比較的良好に位置合わせすることができる。
【０００７】
【発明が解決しようとする課題】
ところで、同一の被写体を異なる時期に撮影すると、撮影時の被写体の体位が変動する場合がある。特に、人体の胸部放射線画像を撮影するときには、患者（人体）の撮影時の立ち位置や向きが変わったり、患者が前傾または後傾してしまうことがあり、その結果、撮影時点が互いに異なる２つの２次元透過画像を比較すると、この体位変動に応じて肋骨等の骨部と血管や気管等の軟部とが異なる方向に移動している場合がある。従来、骨部と軟部とが異なる方向に移動している２つの２次元透過画像間で原画像（或いは、原画像の縮小画像、原画像のボケ画像等）を用いて上述したような全体的な位置合わせを行なうと、軟部組織について全体的に位置合わせされた画像が得られていた。
【０００８】
しかしながら、軟部に全体的に位置合わせされた画像を用いて上述したような局所的な位置合わせを行なうと、テンプレート領域中に存在している肋骨が対応する探索領域中に存在していない場合があり、その結果、肋骨を位置合わせ指標とする局所的な位置合わせの後に作成されるサブトラクション画像中に非常に大きな肋骨の位置ずれアーチファクトが生じ、２つの画像の差異（病変部等）が見難くなってしまうことがあった。
【０００９】
これに対し、画像中の特定の構造物を重視した画像を用いて、特定の構造物同士が合うように位置合わせを行う技術が本出願人により提案されている（特願2001-44623、特願2001-129554）。この技術を利用すれば、骨部組織を重視した画像を用いることにより骨部組織について全体的に位置合わせされた画像を得ることが可能になり、サブトラクション画像中の肋骨のアーチファクトを減少させて骨部組織に存在している病変部を高精度で抽出することができる。
【００１０】
しかしながら、撮影時の体位が変動している２つの２次元透過画像の位置合わせを骨部組織に合わせて行うと、軟部組織の位置合わせ精度が低くなり、軟部組織に存在している病変部がサブトラクション画像中において見難くなってしまう。つまり、撮影時に被写体の体位変動があったときには、２つの２次元透過画像の位置合わせには限界があり、すべての構造物同士が同時に合うように位置合わせを行うことは困難であった。
【００１１】
一方近年では、撮影技術の向上に伴い、ＣＴ装置やＭＲＩ装置などを利用して３次元の透過画像を得ることが可能になっている。このような３次元透過画像を過去に撮影した２次元透過画像と比較して経過を観察したい場合があるが、撮影時の体位が変動していると比較が容易ではない。また、費用の面から、初回は３次元透過画像の撮影を行い２回目以降は比較的安価な２次元透過画像の撮影を行って経過を観察することもある。
【００１２】
本発明は、上記事情に鑑み、撮影時に被写体の体位変動があった場合でも、３次元透過画像と２次元透過画像の位置合わせを行うことを可能にした画像の位置合わせ装置および画像処理装置を提供することを目的とするものである。
【００１３】
【課題を解決するための手段】
本発明の画像の位置合わせ装置は、同一被写体についての２つの画像の位置合わせを行う位置合わせ手段を備えた画像の位置合わせ装置において、２つの画像が３次元透過画像と２次元透過画像であり、位置合わせ手段が、３次元透過画像から２次元透過画像と相関が高い２次元投影画像を作成し、作成された２次元投影画像を２次元透過画像と位置合わせするものであることを特徴とするものである。
【００１４】
ここで、「３次元透過画像から２次元投影画像を作成する」とは、３次元透過画像に投影変換を施して２次元投影画像を得ることを意味するものである。この投影変換の手法としては透視投影法を利用することが好ましい。
【００１５】
また、３次元透過画像から２次元透過画像と相関が高い２次元投影画像を作成する手法としては、例えば、投影変換に先立って、２次元透過画像と相関が高い２次元投影画像が得られるように３次元透過画像に対して３次元的な角度変換等の変換処理を施し、この変換処理が施された３次元透過画像に対して投影変換を施す手法を用いることができる。この３次元的な変換処理としては、３次元アフィン変換を利用することができる。
【００１６】
３次元透過画像としては、ＣＴ装置やＭＲＩ装置を用いて被写体を撮影することにより得られる３次元画像を利用することが可能であり、この他にも、被写体の奥行き情報を有する３次元的な種々の透過画像を利用することができる。また、２次元透過画像とは、被写体を放射線撮影することにより得られる放射線画像などを意味するものである。
【００１７】
本発明の画像処理装置は、同一被写体についての２つの画像の位置合わせを行う位置合わせ手段と、位置合わせ手段により位置合わせがなされた２つの画像を用いて所定の画像処理を行う画像処理手段とを備えた画像処理装置において、同一被写体についての２つの画像が３次元透過画像と２次元透過画像であり、位置合わせ手段が、３次元透過画像から２次元透過画像と相関が高い２次元投影画像を作成し、作成された２次元投影画像を２次元透過画像と位置合わせするものであることを特徴とするものである。
【００１８】
ここで、所定の画像処理としては、例えばサブトラクション処理を挙げることができる。また、同一被写体についての２つの画像を、時系列的に異なる時期に撮影された画像としてもよい。
【００１９】
なお、本発明の画像処理装置における位置合わせ手段は、上記本発明の画像の位置合わせ装置における位置合わせ手段と同等のものである。
【００２０】
【発明の効果】
本発明の画像の位置合わせ装置によれば、同一被写体についての３次元透過画像と２次元透過画像との位置合わせを、３次元透過画像から２次元透過画像と相関が高い２次元投影画像を作成することにより行うから、撮影時に被写体の体位変動があった場合でも、３次元透過画像と２次元透過画像の位置合わせを行うことができる。
【００２１】
すなわち、例えば３次元透過画像に３次元的な変換処理を施した後に投影変換を行う手法などによって、２次元透過画像と相関が高い２次元投影画像を３次元透過画像から作成することができるから、撮影時に被写体の体位変動があった場合でも、その変動に影響されることなくすべての構造物同士が合うように位置合わせを行うことができる。
【００２２】
本発明の画像処理装置によれば、同一被写体についての３次元透過画像と２次元透過画像との位置合わせを、３次元透過画像から２次元透過画像と相関が高い２次元投影画像を作成することにより行い、作成された２次元投影画像と２次元透過画像とを用いて所定の画像処理を行うから、撮影時に被写体の体位変動があった場合でも、すべての構造物同士が合うように位置合わせがなされた２つの画像を用いて所定の画像処理を行うことが可能になる。
【００２３】
なお、所定の画像処理をサブトラクション処理とすれば、すべての構造物同士が合うように位置合わせがなされた２つの画像を用いてサブトラクション処理を行うことができるので、サブトラクション画像中のアーチファクトを減少させることが可能となり、２つの画像間の差異を、その存在位置に拘わらず、見易い状態で抽出することができる。
【００２４】
【発明の実施の形態】
以下、本発明の画像の位置合わせ装置の具体的な実施の形態を図面に基づいて説明する。図１は本発明の画像の位置合わせ装置の一実施形態を示す概略構成図である。
【００２５】
本実施形態の画像の位置合わせ装置に備えられた位置合わせ手段10は、同一被写体についての３次元過去画像データＰｐおよび２次元現在画像データＰｎを入力し、それぞれの１／１０縮小画像データを作成する画像縮小手段11と、画像縮小手段11により作成された３次元過去縮小画像データＰｐ′および２次元現在縮小画像データＰｎ′を入力し、３次元過去縮小画像データＰｐ′に対して３次元アフィン変換による座標変換を施した後、透視投影による投影変換を行い２次元過去画像データＰｐ″を取得する位置合わせ画像取得手段13とを備えた構成である。
【００２６】
位置合わせ手段10に入力される３次元過去画像データＰｐは、ＣＴ装置により撮影された３次元の胸部透過画像を表す画像データであり、２次元現在画像データＰｎよりも時系列的に早い時期に撮影されたものである。また、２次元現在画像データＰｎは通常の放射線撮影装置により撮影された２次元の胸部透過画像を表す画像データである。
【００２７】
次に、本実施形態の画像の位置合わせ装置の作用を説明する。
【００２８】
位置合わせ手段10に３次元過去画像データＰｐおよび２次元現在画像データＰｎが入力されると、２つの画像データＰｐ，Ｐｎは画像縮小手段11に入力され、それぞれ１／１０画像データに縮小される。
【００２９】
画像縮小手段11により作成された３次元過去縮小画像データＰｐ′および２次元現在縮小画像データＰｎ′は、位置合わせ画像取得手段13に入力される。位置合わせ画像取得手段13では、まず最初に、３次元過去縮小画像データＰｐ′に対して３次元アフィン変換による変換処理が施される。
【００３０】
３次元アフィン変換は、画像データに対して３次元的な拡大・縮小、回転、平行移動等の変換処理を施すものであり、１２個の係数の値を変動させることにより、拡大・縮小の比率、回転角度、平行移動量などを変動させることができる。ここでは、予め人体の基本的な３次元的体位変動（前傾・後傾、側湾、側旋等）に対応する係数の組み合わせを用意しておき、この組み合わせによる変換処理（以下、基準変換処理とする。）を３次元過去縮小画像データＰｐ′に対して施す。
【００３１】
基準変換処理が施された３次元過去縮小画像データＰｐ′は、透視投影による投影変換に供される。投影変換とは、３次元画像データに基づいて３次元画像を２次元平面に投影することにより、３次元画像データを２次元画像データに変換するものである。また、透視投影とは、近くのものは大きく、遠くのものは小さく見え、２本の平行線は１点に収れんして見えるといったような人間の感覚に近い性質を持った投影法である。図２は透視投影による投影変換の概念を示す図である。図２に示すように、３次元過去縮小画像Ｐｐ′（以下簡単のため、画像データとその画像データが担持する画像に同じ符号を付して説明する。）を投影中心Ｏから２次元平面に透視投影することにより、３次元過去縮小画像Ｐｐ′の２次元投影画像Ｐｉが得られる。この際、２次元投影画像Ｐｉにおける各画素の画素値は、３次元過去縮小画像Ｐｐ′の奥行き方向（投影方向）の画素値を総和（または平均など）した値とする。
【００３２】
このような手法により、位置合わせ画像取得手段13では基準変換処理が施された３次元過去縮小画像データＰｐ′に基づいて２次元投影画像データＰ1を得、さらに、１２個の係数の値を変動させながら３次元過去縮小画像データＰｐ′に対して３次元アフィン変換による変換処理を順次施すことにより得られた各変換処理後の３次元過去縮小画像データＰｐ′それぞれに基づいて２次元投影画像データＰ2〜ＰNを得、得られた各２次元投影画像データＰｉ（ｉ＝１〜Ｎ）と、２次元現在縮小画像データＰｎ′との相関値をそれぞれ算出し、２次元現在縮小画像データＰｎ′と最も相関が高い２次元投影画像データＰｉを２次元過去画像データＰｐ″として取得する。
【００３３】
本実施形態の画像の位置合わせ装置によれば、３次元過去縮小画像データＰｐ′に対して１２個の係数を変動させながらアフィン変換を施し、各変換ごとに投影変換を行い、各投影変換により得られた投影画像データＰｉそれぞれと２次元現在縮小画像データＰｎ′とを比較して、２次元現在縮小画像データＰｎ′と最も相関が高い投影画像データＰｉを２次元過去画像データＰｐ″として取得するから、撮影時に被写体の体位変動があった場合でも、その体位変動の変動方向や変動量に影響されることなく、すべての構造物同士が合うように位置合わせを行うことが可能になる。
【００３４】
なお、上記実施形態においては、位置合わせ画像を取得する際に画像縮小手段11により縮小された２つの画像データＰｐ′，Ｐｎ′を用いて一連の処理を行う手法を示したが、アフィン変換の係数を変動させながら相関値を算出し、最も相関が高い時のアフィン変換の係数を求める処理のみを縮小画像データに基づいて行い、求められたアフィン変換の係数に基づいて透視投影による投影変換を行って位置合わせ画像を取得する処理を縮小前の画像データを用いて行う形態でもよい。
【００３５】
また、上記位置合わせ画像取得手段13において透視投影により２次元投影画像データＰｉを得る際に、肋骨部分の透過率を他の部分の透過率と変えて重み付け総和をとり、この重み付け総和によって２次元投影画像Ｐｉを作成する形態としてもよい。すなわち、ＣＴ装置と通常の放射線撮影装置は異なる線質の放射線を使用している可能性が高く、単純にＣＴ画像の奥行き方向の画素値（ＣＴ値）の総和（または平均）を取って透視投影画像を作成すると骨部と軟部の画素値（濃度値）の比が通常の放射線画像と異なってしまうことがあるため、上記のように、例えば肋骨部分（骨部）の透過率を変えることにより濃度補正を加えて投影画像を作成することが望ましい。以下、その具体的な手法について説明する。
【００３６】
通常の透視投影法により３次元ＣＴ画像の２次元投影画像を作成する際には、例えば次式、
（ｖ1＋ｖ2＋…＋ｖN）／Ｎ
により、同一視線方向（投影後の２次元画像において同一のｘ，ｙ座標となる線）の画素値（ＣＴ値）の平均をとって２次元投影画像における座標（ｘ，ｙ）の画素値とする。但し、ｖ1〜ｖNは、ＣＴ画像の各奥行き位置（同一視線方向の奥行きをＮ個に分割した位置）におけるＣＴ値である。
【００３７】
一方、重み付け総和により２次元投影画像Ｐｉを作成する際には、まずＭＩＰ法により３次元ＣＴ画像の同一視線方向における最大のＣＴ値を用いて２次元画像を作成し、この２次元画像に対して固定閾値処理を施し低濃度部分を検出することにより画像中に低濃度で表れている肋骨部分の抽出を行う。この際、ＭＩＰ法で２次元画像を作成する際に用いた最大ＣＴ値の奥行き位置を記憶しておくことにより、肋骨が存在している奥行き位置を略特定することができる。なお、特願2000-298100号に記載の方法等により予め胸郭領域を認識し、胸郭領域内のみで肋骨の抽出を行えば、より肋骨部分の抽出精度を向上させることができる。
【００３８】
次に、肋骨が存在している奥行き位置をｊとして、重み付け総和による２次元投影画像の各座標における画素値を次式、
（α1ｖ1＋α2ｖ2＋…＋αjｖj＋…＋αNｖN）／（α1＋α2＋…＋αN）
により得る。但し、α1〜αNは透過率であり、濃度補正の程度に応じてαjの値のみを変動させる。例えば肋骨をより低濃度にするにはαjを他より大きくし、より高濃度にするにはαjを他より小さくすればよい。なお、肋骨の太さは略一定なので、上記奥行き位置ｊは画像サイズとの関係で適応的に変えてもよい。
【００３９】
この手法によれば、肋骨部分の透過率のみを変動させることにより、放射線画像と同等の骨部と軟部の濃度比を有する２次元投影画像を得ることができるので、２次元現在縮小画像データＰｎ′とより相関が高い２次元過去画像データＰｐ″を得ることができる。なお、この手法により採用される３次元画像はＣＴ画像に限定されるものではない。また、上記の実施形態では肋骨部分の透過率を変化させたが、逆に肋骨以外の部分の透過率を変化させてもよい。
【００４０】
なお、２次元現在画像データＰｎと３次元過去画像データＰｐの解像度やサイズ等が異なる場合には、最初に前処理として、この２つの画像データを同じ解像度やサイズに揃えてから位置合わせを行うことが望ましい。
【００４１】
次に、本発明の画像処理装置の具体的な実施の形態を図面に基づいて説明する。図３は本発明の画像処理装置の一実施形態を示す概略構成図である。なお、上記実施の形態による画像の位置合わせ装置と同等の要素については、図１と同じ符号を付してその説明は特に必要のない限り省略する。
【００４２】
本実施形態の画像処理装置は、同一被写体についての３次元過去画像データＰｐおよび２次元現在画像データＰｎを入力し、２次元現在画像データＰｎの１／１０縮小画像データ（２次元現在縮小画像データ）Ｐｎ′と、該２次元現在縮小画像データＰｎ′に位置合わせがなされた２次元過去画像データＰｐ″とを取得する位置合わせ手段10と、位置合わせ手段10から２次元現在縮小画像データＰｎ′と２次元過去画像データＰｐ″を入力し、この２つの画像データを用いてサブトラクション処理を行いサブトラクション画像Ｐｓを取得するサブトラクション処理手段20とにより構成される。
【００４３】
次に、本実施形態の画像処理装置の作用を説明する。
【００４４】
サブトラクション処理手段20は、位置合わせ手段10において位置合わせがなされている２次元現在縮小画像データＰｎ′と２次元過去画像データＰｐ″を入力し、両画像データＰｎ′，Ｐｐ″の画素を対応させて画像間演算を行い、サブトラクション画像データＰｓを生成する。サブトラクション画像データＰｓは、例えば、
Ｐｓ＝ａ１・Ｐｎ′−ａ２・Ｐｐ″  （ただし、ａ１，ａ２は定数）
等の画像間演算により得ることができる。また、下記式で表されるように、両画像間で差異のない部分を中間値（中間濃度、中間輝度等）で再生するため、一定値（１０bitの画像のときは、中間値ｍｉｄ＝５１２）を加算するとともに、差異のコントラストを調整することが好ましい。
【００４５】
Ｐｓ＝（ａ１・Ｐｎ′−ａ２・Ｐｐ″）×cont＋mid
ただし、contはコントラスト係数、midは中間値をそれぞれ表す定数である。
【００４６】
本実施形態による画像処理装置によれば、位置合わせ手段10においてすべての構造物同士が合うように位置合わせがなされた２次元現在縮小画像データＰｎ′と２次元過去画像データＰｐ″を用いてサブトラクション処理を行うので、サブトラクション画像中のアーチファクトを減少させ、２つの画像間の差異を見易くすることができる。
【００４７】
また、位置合わせ手段10を、上記画像の位置合わせ装置の実施形態において記載したような重み付け総和により２次元投影画像Ｐｉを作成する形態としてもよい。こうすれば、位置と濃度についてより相関が高い２つの画像を用いてサブトラクション処理を行うことができるので、サブトラクション画像中のアーチファクトをさらに減少させることができる。
【００４８】
なお、上記実施形態による画像処理装置においては、位置合わせ手段10により作成された２次元過去画像データＰｐ″と２次元現在縮小画像データＰｎ′をそのまま用いてサブトラクション画像データＰｓを生成しているが、位置合わせ手段10により作成された２次元過去画像データＰｐ″と２次元現在縮小画像データＰｎ′に対し、従来の位置合わせ手法における局所的位置合わせを行う（両画像にそれぞれテンプレート領域と探索領域を設定して対応する領域ごとにシフト量を求め、ワーピングを用いて局所的に位置合わせを行う）局所的位置合わせ手段をさらに設け、サブトラクション処理手段20がこの局所的位置合わせがなされた２つの画像を入力して画素を対応させた画像間演算を行い、サブトラクション画像データを生成する形態でもよい。
【００４９】
また、上記各実施形態における位置合わせ手段10により作成された２次元過去画像データＰｐ″を２次元現在縮小画像データＰｎ′（或いは２次元現在画像データＰｎ）とともに表示装置等に例えば並列して出力して表示させ、２次元画像同士を比較読影するようにすることもできる。さらに、サブトラクション画像Ｐｓを並べて表示させてもよい。図４（１）〜（４）は、現在画像Ｐｎ′（或いはＰｎ）、過去画像Ｐｐ″、サブトラクション画像Ｐｓの表示例を示したものである。また、同図（５）〜（８）に示したように、ＣＴスライス画像ＣＴ１，ＣＴ２を同時に並べて表示させてもよい。例えば、サブトラクション画像上で変化が観察された位置をマウス等で指示し、対応する体軸方向のＣＴスライス画像を表示する形態などが考えられる。
【００５０】
なお、上記各実施形態においては、現在画像を２次元の画像とし、現在画像よりも時系列的に早い時期に撮影された過去画像を３次元の画像とした例を示したが、この形態に限定されるものではなく、その逆の形態、すなわち過去画像を２次元の画像とし、現在画像を３次元の画像とする形態でもよい。また、処理時間の短縮のため、位置合わせ手段10において３次元過去画像データＰｐおよび２次元現在画像データＰｎを縮小し、縮小画像に基づいて位置合わせや画像間演算を行う例を示したが、両画像データを縮小しないで各処理を行う形態でもよい。
【図面の簡単な説明】
【図１】本発明の画像の位置合わせ装置の一実施形態を示す図
【図２】透視投影による３次元過去縮小画像Ｐｐ′の投影変換の概念図
【図３】本発明の画像処理装置の一実施形態を示す図
【図４】画像の表示形態の一例を示す図
【符号の説明】
10    位置合わせ手段
11    画像縮小手段
13    位置合わせ画像取得手段
20    サブトラクション処理手段
【図１】
0004104054.tif 000002
【図２】
0004104054.tif 000003
【図３】
0004104054.tif 000004
【図４】
0004104054.tif 000005