JP6807981B2 - 画像位置合わせ装置および方法並びにプログラム - Google Patents

Description

本発明は、複数の画像の位置合わせを行う位置合わせ装置および方法並びにプログラムに関するものである。

従来、心筋、脳、肝臓および膵臓などの血流量または脳内の血流量を計測する手法として、パフュージョン撮影が行われている。

パフュージョン撮影とは、被検体に造影剤を注入し、ＣＴ（Computed Tomography）装置などによって異なるタイミングで複数回撮影する撮影方法である。フェーズ（撮影のタイミング）によって、異なる解剖構造が造影された画像が得られる。具体的には、心臓の場合、初期のフェーズでは造影剤が静脈から注入されることによって、右心房および右心室が造影剤によって染まり、続いて左心室および左心房の順に造影剤が流れ、これにより心臓画像のコントラストが変化する。

パフュージョン撮影によって得られた画像は、パフュージョン解析に用いられる。パフュージョン撮影および解析の主な目的は、対象臓器の血流を解析し、その結果から病気を診断することである。たとえば心筋パフュージョン解析によって得られる血流量は心筋梗塞の診断に用いられる。

パフュージョン（血流）の定量的解析手法にはMaximum Slope法（たとえば非特許文献１参照)やDeconvolution法(たとえば非特許文献２参照)などがある。これらの定量的解析には「時間濃度曲線（Time Density Curve(以下、TDCという)）」が用いられる。TDCとは、臓器のある位置（領域）における画素値（ＣＴ値）の時間変化（造影剤濃度の時間変化）を表したものである。

精度の高い定量的解析を行うには、パフュージョン撮影した画像から正確にTDCを作成する必要がある。そのためには、複数フェーズの画像から、同じ解剖学的位置(領域)の画素値を参照してTDCを作成する必要がある。これには、フェーズ間での対象臓器の位置を正確に合わせる必要があり、画像位置合わせが行われる。

画像の位置合わせは、一方の画像IF(x)（Fixed Image）の点xと、他方の画像IM(x)（Moving Image）の点T(x)とが、解剖学的に一致するように幾何学的変換Tを求める問題である。幾何学的変換Tをパラメタμで記述した場合には、下式（１）に基づいて、最適解を求める問題となる。幾何学的変換Tには、剛体変換（平行移動および回転）、Affine変換および非剛体変換（B-SplineおよびDiffeomorphism）などがある。（たとえば非特許文献３参照）

なお、上式（１）におけるSは画像IF(x)と画像IM(T(・；x))の一致度を表す評価関数であり、ここでは、Sが小さくなるほど２画像間は一致している。一致度の評価指標としては、Sum of Squared Difference(SSD; ２画像間の同位置における画素値の差の二乗の合計)や相互情報量（画素値の間の統計的な依存関係を定量化する指標の１つ）などの種々の評価指標が知られている。

特開２０１１−１２５４３１号公報

"Measurement of tissue perfusion by dynamic computed tomography.", Miles KA, The British journal of radiology (1991), Vol.64, No. 761, pp409-412 "Deconvolution-based CT and MR brain perfusion measurement: theoretical model revisited and practical implementation details", Fieselmann, Andreas and Kowarschik, Markus and Ganguly, Arundhuti and Hornegger, Joachim and Fahrig, Rebecca, Journal of Biomedical Imaging(2011), Vol. 2011, pp14 "Survey of medical image registration", Mani, VRS and others, Journal of Biomedical Engineering and Technology(2013), Vol. 1, No. 2, pp.8-25

ここで、上述したような画像の位置合わせを行う場合、異なる解剖構造が造影された画像同士の位置合わせが問題となる。たとえば一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合、同じ領域が有する画素値が異なっているため、その一致度が低く評価され、正確に位置合わせすることができない。

このような問題を解決する方法として、画像全体に画素値変換を行う方法が考えられる。具体的には、高い画素値を低い画素値に変換する方法が考えられる。

しかしながら、この手法の場合、位置合わせ対象臓器の内部の構造の情報を失う。たとえば心筋パフュージョンの画像の場合、２つの画像間のうち一方の画像は、右心房および右心室だけが造影剤によって高い画素値を示し（図２Ｉ参照）、他方の画像が右心房および右心室だけでなく左心房および左心室も高い画素値を示す（図２II参照）場合がある。このような場合、画像全体に同じ画素値変換を施すと、一方の画像しか高い画素値を示さない左心房および左心室に関しては上述の問題が解決するので良いが、右心室および右心房はどちらも高い画素値を示しているので、位置合わせに用いることができる形状情報が失われてしまい、臓器内部の位置を合わせることができなくなる。

なお、特許文献１には、ＣＴ画像と超音波画像といった異なるモダリティで撮影された画像同士の位置合わせに関する技術が提案されているが、上述した問題を解決する方法については何も提案されていない。

本発明は、上記事情に鑑み、一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合でも、画像間のコントラスト差を低減することができ、高精度に位置合わせを行うことができる画像位置合わせ装置および方法並びにプログラムを提供することを目的とするものである。

本発明の第１の画像位置合わせ装置は、位置合わせ処理対象を含む複数の画像を取得する画像取得部と、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する分割領域設定部と、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の画素値に基づいて、各分割領域の画素値変換方法を決定する画素値変換方法決定部と、画素値変換方法決定部によって決定された分割領域の画素値変換方法を用いて、分割領域の組の分割領域内の画像に画素値変換を施す画素値変換部と、画素値変換部によって画素値変換の施された複数の画像に対して位置合わせ処理を施す位置合わせ処理部とを備える。

また、上記本発明の第１の画像位置合わせ装置においては、画素値変換方法決定部が、組の各分割領域の画素値に基づいて、各分割領域の画素値変換を行うか否かを決定することができる。

本発明の第２の画像位置合わせ装置は、位置合わせ処理対象を含む複数の画像を取得する画像取得部と、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する分割領域設定部と、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す位置合わせ処理部とを備え、位置合わせ処理部が、分割領域の組毎に、その組の各分割領域の画素値に基づいて評価関数を変更して位置合わせ処理を施す。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、位置合わせ処理対象として心臓画像を含む画像を取得し、分割領域設定部は、分割領域として、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域とを設定することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、分割領域設定部は、複数の画像に含まれる大動脈の領域の画素値の変化に基づいて、第１の分割領域と第２の分割領域を設定することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、分割領域設定部は、複数の画像に含まれる大動脈の領域の画素値の変化に基づいて、複数の画像の中から、左心房および左心室の領域の画素値よりも右心房および右心室の領域の画素値の方が高い第１の画像と右心房および右心室の領域の画素値よりも左心房および左心室の領域の画素値の方が高い第２の画像とを特定し、第１の画像を用いて第１の分割領域を設定し、第２の画像を用いて第２の分割領域を設定することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、位置合わせ処理対象として心臓画像を含む画像を取得し、分割領域設定部は、分割領域として、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域と心筋を含む第３の分割領域とを設定することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、位置合わせ処理対象として心臓画像を含む画像を取得し、分割領域設定部は、分割領域として、右心房を含む第１の分割領域と右心室を含む第２の分割領域と左心房を含む第３の分割領域と左心室を含む第４の分割領域とを設定することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、位置合わせ処理対象として心臓画像を含む画像を取得し、分割領域設定部は、分割領域として、右心房を含む第１の分割領域と右心室を含む第２の分割領域と左心房を含む第３の分割領域と左心室を含む第４の分割領域と心筋を含む第５の分割領域とを設定することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、位置合わせ処理対象として心臓画像を含む画像を取得し、分割領域設定部は、心臓の血流方向に垂直な断面によって心臓画像を分割して複数の分割領域を設定する。

また、上記本発明の第１および第２の画像位置合わせ装置において、分割領域設定部は、分割領域に対応する領域が予め設定された基準画像を用いて、複数の画像のうちの１つの画像に対して分割領域を設定し、分割領域の設定結果を用いて、上記１つの画像以外の画像に対して分割領域を設定することができる。

また、上記本発明の第１の画像位置合わせ装置において、画素値変換部は、分割領域内の画像の画素値を低くする画素値変換を施すことができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、異なる時点で撮影された複数の画像を取得することができる。

また、上記本発明の第１および第２の画像位置合わせ装置において、画像取得部は、位置合わせ処理対象として心臓画像を含む画像を取得し、かつ複数の心拍における同じ心位相のタイミングで撮影された複数の画像を取得することができる。

本発明の第１の画像位置合わせ方法は、位置合わせ処理対象を含む複数の画像を取得し、複数の画像の各位置合わせ処理対象を分割して分割領域を設定し、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の画素値に基づいて、各分割領域の画素値変換方法を決定し、その決定した分割領域の画素値変換方法を用いて、分割領域の組の分割領域内の画像に画素値変換を施し、画素値変換の施された複数の画像に対して位置合わせ処理を施す。

本発明の第２の画像位置合わせ方法は、位置合わせ処理対象を含む複数の画像を取得し、複数の画像の各位置合わせ処理対象を分割して分割領域を設定し、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す際、分割領域の組毎に、その組の各分割領域の画素値に基づいて評価関数を変更して位置合わせ処理を施す。

本発明の第１の画像位置合わせプログラムは、コンピュータを、位置合わせ処理対象を含む複数の画像を取得する画像取得部と、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する分割領域設定部と、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の画素値に基づいて、各分割領域の画素値変換方法を決定する画素値変換方法決定部と、画素値変換方法決定部によって決定された分割領域の画素値変換方法を用いて、分割領域の組の分割領域内の画像に画素値変換を施す画素値変換部と、画素値変換部によって画素値変換の施された複数の画像に対して位置合わせ処理を施す位置合わせ処理部として機能させる。

本発明の第２の画像位置合わせプログラムは、コンピュータを、位置合わせ処理対象を含む複数の画像を取得する画像取得部と、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する分割領域設定部と、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す位置合わせ処理部として機能させる画像位置合わせプログラムであって、位置合わせ処理部が、分割領域の組毎に、その組の各分割領域の画素値に基づいて評価関数を変更して位置合わせ処理を施す。

本発明の第１の画像位置合わせ装置および方法並びにプログラムによれば、位置合わせ処理対象を含む複数の画像を取得し、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する。そして、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の画素値に基づいて、各分割領域の画素値変換方法を決定し、その決定した分割領域の画素値変換方法を用いて、分割領域の組の分割領域内の画像に画素値変換を施し、その画素値変換の施された複数の画像に対して位置合わせ処理を施す。このように分割領域毎に画素値変換方法を決定して画素値変換を行うことによって、一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合でも、画像間のコントラスト差を低減することができ、高精度に位置合わせを行うことができる。

本発明の第２の画像位置合わせ装置および方法並びにプログラムによれば、位置合わせ処理対象を含む複数の画像を取得し、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する。そして、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施し、この際、分割領域の組毎に、その組の各分割領域の画素値に基づいて評価関数を変更して位置合わせ処理を施す。このように分割領域毎に評価関数を変更して位置合わせ処理を行うことによって、一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合でも、画像間のコントラスト差を低減することができ、高精度に位置合わせを行うことができる。

本発明の画像位置合わせ装置の第１の実施形態を用いた医用画像診断支援システムの概略構成を示すブロック図 心臓をパフュージョン撮影した画像の一例を示す図 パフュージョン解析の一例を説明するための図 本発明の画像位置合わせ装置の第１の実施形態を用いた医用画像診断支援システムの作用を説明するためのフローチャート 心臓の血流方向に垂直な断面で心臓画像を分割して設定した分割領域の一例を示す図 図５に示す分割領域毎の画素値に基づいて、画素値変換方法を決定する方法を説明するための図 本発明の画像位置合わせ装置の第２の実施形態を用いた医用画像診断支援システムの概略構成を示すブロック図 本発明の画像位置合わせ装置の第２の実施形態を用いた医用画像診断支援システムの作用を説明するためのフローチャート

以下、本発明の画像位置合わせ装置および方法並びにプログラムの第１の実施形態を用いた医用画像診断支援システムについて、図面を参照しながら詳細に説明する。図１は、本実施形態の医用画像診断支援システムの概略構成を示すブロック図である。

本実施形態の医用画像診断支援システムは、心臓などを異なる時点で撮影した２以上の画像を取得し、これらの画像の位置合わせを行った後に、その２以上の画像を用いてパフュージョン解析を行うものである。

本実施形態の医用画像診断支援システムは、具体的には、図１に示すように、画像位置合わせ装置１と、医用画像保管サーバ２と、表示装置３と、入力装置４とを備えている。

画像位置合わせ装置１は、コンピュータに本実施形態の画像位置合わせプログラムをインストールしたものである。

画像位置合わせ装置１は、中央処理装置（ＣＰＵ（central processing unit））、半導体メモリ、およびハードディスクやＳＳＤ(Solid State Drive)等のストレージデバイスを備えている。ストレージデバイスには、本実施形態の画像位置合わせプログラムがインストールされており、この画像位置合わせプログラムが中央処理装置によって実行されることによって、図１に示す画像取得部１０、分割領域設定部１１、画素値変換方法決定部１２、画素値変換部１３、位置合わせ処理部１４、パフュージョン解析部１５および表示制御部１６が動作する。

画像位置合わせプログラムは、ＤＶＤ（Digital Versatile Disc）およびＣＤ−ＲＯＭ（Compact Disc Read Only Memory）などの記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、画像位置合わせプログラムは、ネットワークに接続されたサーバコンピュータの記憶装置もしくはネットワークストレージに対して外部からアクセス可能な状態で記憶され、要求に応じてコンピュータにダウンロードされ、インストールされる。

画像取得部１０は、予め撮影された３次元画像６を取得するものである。３次元画像６は、たとえばＣＴ装置やＭＲＩ（Magnetic Resonance Imaging）装置などによって患者を撮影したものである。本実施形態の画像取得部１０は、いわゆるパフュージョン撮影によって撮影された３次元画像６を取得するものであり、造影剤が注入された心臓を時系列に撮影した２以上のフェーズの心臓画像を取得する。

より具体的には、画像取得部１０は、撮影間隔が１〜数心拍の２以上の心臓画像を３次元画像６として取得する。２以上の心臓画像は、各心拍の収縮期において撮影されたものであることが望ましい。各心拍の収縮期において撮影された各心臓画像に含まれる心臓の形状に大きな違いはなく、１ｃｍ前後のずれである。この１ｃｍ程度のずれが、後述する位置合わせ処理によって調整される。なお、本実施形態においては収縮期の心臓画像を取得するようにしたが、これに限らず、心臓の形状がほぼ一致するタイミングで取得された心臓画像であればその他のフェーズの心臓画像を取得するようにしてもよい。

本実施形態においては、造影剤が注入された心臓を撮影した３次元画像６を取得するようにしたが、撮影対象は心臓に限らず、造影剤が注入された肝臓、膵臓または脳などを撮影した画像を３次元画像６として取得するようにしてもよい。

また、３次元画像６としては、アキシャル断層画像、サジタル断層画像およびコロナル断層画像などの断層画像からなるボリュームデータを取得してもよいし、断層画像単体を取得するようにしてもよい。

３次元画像６は、医用画像保管サーバ２に患者の識別情報とともに予め保管されており、画像取得部１０は、入力装置４などを用いてユーザによって入力された患者の識別情報に基づいて、その識別情報を有する３次元画像６を医用画像保管サーバ２から読み出して一時記憶するものである。

分割領域設定部１１は、画像取得部１０によって取得された２以上の３次元画像６に含まれる位置合わせ処理対象を分割して分割領域を設定するものである。本実施形態における位置合わせ処理対象は心臓である。

本実施形態の分割領域設定部１１は、分割領域として、各３次元画像６内の心臓領域における右心房および右心室を含む第１の分割領域と、左心房および左心室を含む第１の分割領域とを設定する。

図２Ｉ〜図２IIIは、各心拍の収縮期において撮影された各心臓画像の一例を示す図である。図２Ｉ〜図２IIIにおける領域Ｒ１が右心房の領域であり、領域Ｒ２が右心室の領域であり、領域Ｌ１が左心房の領域であり、領域Ｌ２が左心室の領域である。なお、図２Ｉについては、左心室と心筋との境界が明確でないので、これらの領域の境界は図示省略している。第１の分割領域として、領域Ｒ１と領域Ｒ２とを合わせた領域が設定され、第１の分割領域として、領域Ｌ１と領域Ｌ２とを合わせた領域が設定される。また、領域Ｍは心筋の領域であり、領域Ｖは大動脈の領域であるが、これらの領域については、後で詳述する。

分割領域設定部１１には、心臓のアトラス画像（本発明の基準画像に相当する）が予め設定されている。アトラス画像とは、予め心臓を撮影した標準的な３次元画像と、右心房および右心室を含む領域および左心房および左心室を含む領域が定義されたラベル画像とが対応付けられたものである。

分割領域設定部１１は、上述したアトラス画像を用いて、３次元画像６に対して第１の分割領域と第２の分割領域とを設定する。具体的には、分割領域設定部１１は、まず、２以上の３次元画像６のうちの１つの代表する３次元画像６と上述したアトラス画像との位置合わせ処理を行うことによって、上記代表する３次元画像６に対して第１の分割領域と第２の分割領域とを設定する。なお、位置合わせ処理としては、既に公知な剛体位置合わせ処理または非剛体位置合わせ処理を用いることができる。

そして、代表する３次元画像６の第１の分割領域と第２の分割領域の設定結果を用いて、その他の３次元画像６に対して第１の分割領域および第２の分割領域を設定する。具体的には、代表する３次元画像６に設定された第１の分割領域と第２の分割領域の座標値をその他の３次元画像６に対応付けることによって第１の分割領域および第２の分割領域を設定する。このように１つの代表する３次元画像６の分割領域の設定結果をその他の３次元画像６に適用することによって、分割領域の設定処理の負荷を下げることができ、処理の高速化を図ることができる。

なお、本実施形態においては、代表する３次元画像６の第１の分割領域および第２の分割領域の設定結果を用いてその他の３次元画像６の第１の分割領域および第２の分割領域を設定するようにしたが、これに限らず、代表する３次元画像６だけでなく、その他の３次元画像６についても、アトラス画像を用いて第１の分割領域および第２の分割領域を設定するようにしてもよい。

画素値変換方法決定部１２は、各３次元画像６に含まれる心臓の同じ位置に設定された分割領域の組毎に、その組の各分割領域の画素値に基づいて、各分割領域の画素値変換方法を決定するものである。なお、本明細書における画素値変換方法としては、画素値変換を行わないことも含むものとする。

一般的に、心臓画像における造影された部分の画素値は１００以上に増加するため、画素値変換方法決定部１２は、この数値を用いて各分割領域が造影されているか否かを判定し、その判定結果に基づいて画素値変換方法を決定する。なお、以降、説明を分かりやすくするために、２以上の３次元画像６のうちの２つの３次元画像６の関係に絞って説明する。

具体的には、画素値変換方法決定部１２は、第１の３次元画像６の第１の分割領域の平均画素値と第２の３次元画像６の第１の分割領域の平均画素値とを算出する。そして、第１の３次元画像６の第１の分割領域の平均画素値が１００以上であって、かつ第２の３次元画像６の第１の分割領域の平均画素値が１００未満である場合には、すなわち第１の３次元画像６と第２の３次元画像６のそれぞれの第１の分割領域のコントラストが異なる場合には、第１の３次元画像６の第１の分割領域の画素値に対して画素値変換を行うことを決定する。

また、逆に、第１の３次元画像６の第１の分割領域の平均画素値が１００未満であって、かつ第２の３次元画像６の第１の分割領域の平均画素値が１００以上である場合には、第２の３次元画像６の第１の分割領域の画素値に対して画素値変換を行うことを決定する。

一方、第１の３次元画像６の第１の分割領域の平均画素値と第２の３次元画像６の第１の分割領域の平均画素値との両方が１００以上または１００未満である場合には、すなわち両方の第１の分割領域のコントラストが近い場合には、第１の３次元画像６と第２の３次元画像６の第１の分割領域には、画素値変換を行わないことを決定する。第１の３次元画像６の第２の分割領域の平均画素値と第２の３次元画像６の第２の分割領域の平均画素値との両方が１００以上または１００未満である場合には、第１の分割領域と同様に、第１の３次元画像６と第２の３次元画像６の第２の分割領域には、画素値変換を行わないことを決定する。

画素値変換部１３は、画素値変換方法決定部１２によって決定された画素値変換方法を用いて、第１の３次元画像６および第２の３次元画像６の第１の分割領域に画素値変換を施し、第１の３次元画像６および第２の３次元画像６の第２の分割領域に画素値変換を施す。具体的には、画素値変換方法決定部１２において画素値変換を行うと決定された分割領域内の画像に対して、下式（２）に基づいて画素値変換を行う。なお、下式（２）におけるｘは画素値であり、ｆ（ｘ）は画素値変換関数である。下式（２）に基づいて画素値変換を行うことによって、第１の３次元画像６および第２の３次元画像６の第１の分割領域間または第１の３次元画像６および第２の３次元画像６の第２の分割領域間のコントラスト差を小さくすることができる。

位置合わせ処理部１４は、画素値変換部１３によって画素値変換の施された第１の３次元画像６と第２の３次元画像６に対して位置合わせ処理を施すものである。位置合わせ処理としては、剛体位置合わせ処理または非剛体位置合わせ処理が施される。剛体位置合わせ処理としては、たとえばＩＣＰ（Iterative Closest Point）法を用いた処理などを用いることができるが、その他の公知な手法を用いるようにしてもよい。また、非剛体位置合わせ処理としては、たとえばＦＦＤ（Free-Form Deformation）法を用いた処理やＴＰＳ（Thin-Plate Spline）法を用いた処理などを用いることができるが、その他の公知な手法を用いるようにしてもよい。また、剛体位置合わせ処理を行った後に、非剛体位置合わせ処理を行うようにしてもよい。

なお、上記説明では、第１の３次元画像６および第２の３次元画像６の２つの３次元画像の位置合わせ処理について説明したが、実際には、パフュージョン解析に必要な多数（たとえば３０以上）の３次元画像６に対して分割領域の設定、画素値変換方法の決定、画素値変換および位置合わせ処理が行われる。

パフュージョン解析部１５は、位置合わせ処理部１４によって位置合わせ処理の施された２以上の３次元画像６を用いて、パフュージョン解析を行うものである。パフュージョン解析の手法としては、上述したようなMaximum Slope法またはDeconvolution法が用いられる。具体的には、図３Ｉに示すように、ユーザによって表示装置３に表示された３次元画像６上において関心領域ＩＲが指定される。パフュージョン解析部１５は、位置合わせ処理の施された２以上の３次元画像６のそれぞれの関心領域ＩＲ内の平均濃度値を算出し、図３IIに示すような時間濃度曲線（Time Density Curve）を算出し、この時間濃度曲線に基づいてパフュージョン解析を行う。

表示制御部１６は、液晶ディスプレイなどの表示デバイスを備えたものであり、パフュージョン解析部１５による解析結果を表示装置３に表示させるものである。また、表示制御部１６は、画像取得部１０によって取得された２以上の３次元画像６および位置合わせ処理後の２以上の３次元画像６などを表示装置３に表示させるものである。

入力装置４は、ユーザによる種々の設定入力を受け付けるものであり、キーボードやマウスなどの入力デバイスを備えたものである。入力装置４は、たとえば患者の識別情報の設定入力および上述した関心領域ＩＲの設定入力などを受け付けるものである。

なお、タッチパネルを用いることによって表示装置３と入力装置４と兼用するようにしてもよい。

次に、本実施形態の医用画像診断支援システムの作用について、図４に示すフローチャートを参照しながら説明する。

まず、ユーザによる患者の識別情報などの入力に基づいて、その患者の心臓をパフュージョン撮影した２以上の３次元画像６が画像取得部１０によって取得される（Ｓ１０）。

画像取得部１０によって取得された2以上の３次元画像６は分割領域設定部１１に入力される。分割領域設定部１１は、２以上の３次元画像６に対してそれぞれ第１の分割領域および第２の分割領域を設定する（Ｓ１２）。

次に、画素値変換方法決定部１２は、各３次元画像６の第１の分割領域の平均画素値に基づいて、第１の分割領域に対する画素値変換方法を決定し、各３次元画像６の第２の分割領域の平均画素値に基づいて、第２の分割領域に対する画素値変換方法を決定する（Ｓ１４）。

そして、画素値変換方法決定部１２によって決定された画素値変換方法を用いて、画素値変換部１３において、各３次元画像６の第１の分割領域および第２の分割領域に対して画素値変換が施される（Ｓ１６）。

次いで、画素値変換が施された３次元画像６が位置合わせ処理部１４に入力され、位置合わせ処理部１４によって位置合わせ処理が施される（Ｓ１８）。

位置合わせ処理の施された３次元画像６はパフュージョン解析部１５に入力され、パフュージョン解析部１５において、入力された３次元画像６を用いてパフュージョン解析が行われる（Ｓ２０）。

パフュージョン解析部１５によるパフュージョン解析結果は表示制御部１６に入力され、表示制御部１６は、入力されたパフュージョン解析結果を表示装置３に表示させる（Ｓ２２）。

上記実施形態の医用画像診断支援システムによれば、位置合わせ処理対象を含む複数の画像を取得し、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する。そして、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の画素値に基づいて、各分割領域の画素値変換方法を決定し、その決定した分割領域の画素値変換方法を用いて、分割領域の組の分割領域内の画像に画素値変換を施し、その画素値変換の施された複数の画像に対して位置合わせ処理を施す。このように分割領域毎に画素値変換方法を決定して画素値変換を行うことによって、一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合でも、画像間のコントラスト差を低減することができ、高精度に位置合わせを行うことができる。

なお、上記第１の実施形態においては、分割領域設定部１１において、分割領域として、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域とを設定するようにしたが、分割領域の設定方法としては、これに限らない。たとえば分割領域として、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域と心筋を含む第３の分割領域とを設定するようにしてもよい。

または、分割領域として、右心房を含む第１の分割領域と右心室を含む第２の分割領域と左心房を含む第３の分割領域と左心室を含む第４の分割領域とを設定するようにしてもよい。もしくは、分割領域として、右心房を含む第１の分割領域と右心室を含む第２の分割領域と左心房を含む第３の分割領域と左心室を含む第４の分割領域と心筋を含む第５の分割領域とを設定するようにしてもよい。

また、上述したように心臓の解剖学的な構造に基づいて分割領域を設定するのではなく、画素値を変化させる造影剤が流れる方向を考慮して、心臓の血流方向に垂直な断面で心臓を分割して分割領域を設定するようにしてもよい。図５は、心臓の血流方向に垂直な断面で心臓を分割して分割領域を設定した一例を示す図である。図５に示す矢印Ａ１は、右心房および右心室における血流方向を示すものであり、矢印Ａ２は、左心房および左心室における血流方向を示すものであり、Ｄ１〜Ｄ１０は、心臓の血流方向に垂直な断面で心臓を分割して設定された分割領域を示すものである。以下、図５に示すように分割領域Ｄ１〜Ｄ１０を設定した場合の各分割領域の画素値変換方法の決定方法について説明する。

画素値変換方法決定部１２は、まず、第１の３次元画像６および第２の３次元画像６の双方について、血流方向に沿って各分割領域の平均画素値を調べる。具体的には、図５の例の場合、矢印Ａ１に沿って各分割領域の平均画素値を調べた後、矢印Ａ２に沿って各分割領域の平均画素値を調べる。そして、平均画素値が予め設定された閾値以上となる分割領域を特定する。すなわち造影剤によって画素値が高くなっている分割領域を特定する。

そして、画素値変換方法決定部１２は、第１の３次元画像６と第２の３次元画像６との間で同じ位置の分割領域について、いずれか一方の分割領域のみの平均画素値が閾値以上である場合、画素値変換を行う分割領域として決定する。逆に、第１の３次元画像６と第２の３次元画像６との間で同じ位置の分割領域について、両方の分割領域の平均画素値が閾値以上または閾値未満である場合には、画素値変換を行わない分割領域として決定する。

図６は、上述した画素値変換方法の決定方法を説明するための図である。図６に示すように、たとえば第１の３次元画像６の分割領域Ｄ１〜Ｄ１０のうち分割領域Ｄ３〜Ｄ６の平均画素値が閾値以上であり、第２の３次元画像６の分割領域Ｄ１〜Ｄ１０のうち分割領域Ｄ５〜Ｄ８の平均画素値が閾値以上である場合には、分割領域Ｄ３、分割領域Ｄ４、分割領域Ｄ７および分割領域Ｄ８については、画素値変換を行う分割領域として決定する。一方、分割領域Ｄ１、分割領域Ｄ２、分割領域Ｄ５、分割領域Ｄ６、分割領域Ｄ９および分割領域Ｄ１０については、画素値変換を行わない分割領域として決定する。

なお、画素値変換方法を決定した後の画素値変換および位置合わせ処理については、上記第１の実施形態と同様である。

次に、本発明の第２の実施形態を用いた医用画像診断支援システムについて説明する。図７は、第２の実施形態の医用画像診断支援システムの概略構成を示すブロック図である。第２の実施形態の医用画像診断支援システムは、第１の実施形態の医用画像診断支援システムとは、画像位置合わせ装置５の構成が異なる。その他の構成については、第１の実施形態の構成と同様である。

第２の実施形態の画像位置合わせ装置５は、画像取得部５０と、分割領域設定部５１と、位置合わせ処理部５２と、パフュージョン解析部５３と、表示制御部５４とを備えている。画像取得部５０、分割領域設定部５１、パフュージョン解析部５３および表示制御部５４については、上記第１の実施形態の画像取得部１０、分割領域設定部１１、パフュージョン解析部１５および表示制御部１６と同様であるので説明を省略する。

第２の実施形態の位置合わせ処理部５２は、第１の３次元画像６と第２の３次元画像６の同じ位置に設定された分割領域の組毎に、その組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施すものである。

そして、位置合わせ処理部５２は、各分割領域の一致度を評価関数によって評価して位置合わせ処理を行うという点では、第１の実施形態の位置合わせ処理と同様であるが、分割領域の組毎に、その組の各分割領域の画素値に基づいて評価関数を変更して位置合わせ処理を施す点で第１の実施形態の位置合わせ処理とは異なる。

位置合わせ処理部５２は、具体的には、第１の３次元画像６の第１の分割領域の平均画素値と第２の３次元画像６の第１の分割領域の平均画素値とを算出する。そして、第１の３次元画像６の第１の分割領域の平均画素値が１００以上であって、かつ第２の３次元画像６の第１の分割領域の平均画素値が１００未満である場合には、すなわち第１の３次元画像６と第２の３次元画像６のそれぞれの第１の分割領域のコントラストが異なる場合には、下式（３）の評価関数を用いて一致度を算出して位置合わせ処理を行う。

上式（３）におけるS１は一致度の評価関数であり、画像間の対応する分割領域の画素値が画像間で異なる場合に、各画像の分割領域の位置が一致しているほどS１は小さな値となる。また、Nは画像間(Fixed画像とMoving画像の間)で比較する画素の総数、fは、上述した式（２）の画素値変換関数、xiは座標点、IF(xi)はFixed画像のxi における画素値、IM(xi)はMoving画像のxi における画素値、Tは座標点の幾何学的変換関数、μは幾何学的変換のパラメタである。なお、第１の３次元画像６および第２の３次元画像６のいずれか一方をFixed画像とし、他方をMoving画像とする。

また、逆に、第１の３次元画像６の第１の分割領域の平均画素値が１００未満であって、かつ第２の３次元画像６の第１の分割領域の平均画素値が１００以上である場合にも、上式（３）の評価関数を用いて一致度を算出して位置合わせ処理を行う。

一方、第１の３次元画像６の第１の分割領域の平均画素値と第２の３次元画像６の第１の分割領域の平均画素値との両方が１００以上または１００未満である場合には、すなわち両方の第１の分割領域のコントラストが近い場合には、下式（４）の評価関数を用いて一致度を算出して位置合わせ処理を行う。また、第１の３次元画像６の第２の分割領域の平均画素値と第２の３次元画像６の第２の分割領域の平均画素値との両方が１００以上または１００未満である場合にも、第１の分割領域と同様に、下式（４）の評価関数を用いて一致度を算出して位置合わせ処理を行う。

上式（４）におけるS２は一致度の評価関数であり、画素値が一致しているほど値が小さい。また、Nは画像間(Fixed画像とMoving画像の間)で比較する画素の総数、fは、上述した式（２）の画素値変換関数、xiは座標点、IF(xi)はFixed画像のxi における画素値、IM(xi)はMoving画像のxi における画素値、Tは座標点の幾何学的変換関数、μは幾何学的変換のパラメタである。

なお、位置合わせ処理自体は、たとえば上述した非特許文献３に記載の方法を用いることができる。

次に、第２の実施形態の医用画像診断支援システムの作用について、図８に示すフローチャートを参照しながら説明する。

まず、ユーザによる患者の識別情報などの入力に基づいて、その患者の心臓をパフュージョン撮影した２以上の３次元画像６が画像取得部５０によって取得される（Ｓ３０）。

画像取得部５０によって取得された2以上の３次元画像６は分割領域設定部５１に入力される。分割領域設定部５１は、２以上の３次元画像６に対してそれぞれ第１の分割領域および第２の分割領域を設定する（Ｓ３２）。

次に、位置合わせ処理部５２は、各３次元画像６の第１の分割領域の平均画素値に基づいて、第１の分割領域に対する評価関数を決定し、各３次元画像６の第２の分割領域の平均画素値に基づいて、第２の分割領域に対する評価関数を決定し、第１の分割領域と第２の分割領域とでそれぞれ決定した評価関数を用いて一致度を算出することによって、各３次元画像６に対して位置合わせ処理を施す（Ｓ３４）。

位置合わせ処理の施された３次元画像６はパフュージョン解析部５３に入力され、パフュージョン解析部５３において、入力された３次元画像６を用いてパフュージョン解析が行われる（Ｓ３６）。

パフュージョン解析部５３によるパフュージョン解析結果は表示制御部５４に入力され、表示制御部５４は、入力されたパフュージョン解析結果を表示装置３に表示させる（Ｓ３８）。

第２の実施形態の医用画像診断支援システムによれば、位置合わせ処理対象を含む複数の画像を取得し、複数の画像の各位置合わせ処理対象を分割して分割領域を設定する。そして、各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、その組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施し、この際、分割領域の組毎に、その組の各分割領域の画素値に基づいて評価関数を変更して位置合わせ処理を施す。このように分割領域毎に評価関数を変更して位置合わせ処理を行うことによって、一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合でも、画像間のコントラスト差を低減することができ、高精度に位置合わせを行うことができる。

なお、上記第２の実施形態においては、各分割領域の画素値に基づいて、上式（３）の評価関数と上式（４）の評価関数とを切り替えて使用するようにしたが、評価関数としてはこれに限らず、上式（３）の代わりに下式（５）を用いるようにしてもよい。

上式（５）におけるSは一致度の評価関数であり、画素値が一致しているほど値が小さい。また、Nは画像間(Fixed画像とMoving画像の間)で比較する画素の総数、xiは座標点、IF(xi)はFixed画像のxi における画素値、IM(xi)はMoving画像のxi における画素値、Tは座標点の幾何学的変換関数、μは幾何学的変換のパラメタ、ｋおよびσは定数である。

上式（５）におけるSは上式（３）と同様に、一致度の評価関数であり、画像間の対応する分割領域の画素値が異なる場合に、各画像の分割領域の位置が一致しているほどSは小さな値となる。また、Nは画像間(Fixed画像とMoving画像の間)で比較する画素の総数、xiは座標点、IF(xi)はFixed画像のxi における画素値、IM(xi)はMoving画像のxiにおける画素値、Tは座標点の幾何学的変換関数、μは幾何学的変換のパラメタ、並びにｋおよびσは定数である。上式（５）のSは、画素値の差が定数σによって決まる所定の値に近いほど小さくなり、一致度が高くなるような関数である。

なお、上記第１〜第３の実施形態においては、分割領域設定部１１は、心臓のアトラス画像を用いて分割領域を設定するようにしたが、分割領域の設定方法としては、これに限らない。以下、その他の分割領域の設定方法について説明する。

まず、大動脈は左心室につながっているため、造影剤による画素値の変化が左心房および左心室に近い。このことを利用して、左心房および左心室の画素値が造影剤によって高くなっているフェーズの画像を特定する。

具体的には、まず、２以上のフェーズの３次元画像６からそれぞれ大動脈の領域を抽出する。上述したように図２Ｉ〜図２IIIに示すＶの領域が大動脈の領域である。なお、大動脈の領域の抽出方法については、公知の画像処理を用いることができる。

次に、各３次元画像６の大動脈の領域の平均画素値の変化に基づいて、左心房および左心室の画素値が右心房および右心室の画素値よりも高くなっているフェーズの３次元画像を特定する。具体的には、各３次元画像の大動脈の領域の平均画素値の変化に基づいて、大動脈の領域の平均画素値が最大となるフェーズの３次元画像６を特定する。そして、その特定した３次元画像６における画素値が予め設定された閾値以上である領域を特定することによって左心房および左心室の領域を抽出する。

次いで、大動脈の領域の平均画素値が予め設定された下限値以上である３次元画像のうち、大動脈の領域の平均画素値が最小であるフェーズの３次元画像を、右心房および右心室の画素値が左心房および左心室の画素値よりも高くなっているフェーズの３次元画像として特定する。そして、その特定した３次元画像６における画素値が予め設定された閾値以上である領域を特定することによって右心房および右心室の領域を抽出する。

そして、上述したようにして抽出された右心房および右心室の領域を第１の分割領域として設定し、左心房および右心室の領域を第２の分割領域として設定する。

また、第１の分割領域および第２の分割領域の設定方法としては、上記の方法に限らず、たとえば公知な輪郭抽出などの画像処理手法を用いることによって、３次元画像６から右心房および右心室の領域と左心房および左心室の領域を抽出し、その抽出結果を用いて第１の分割領域および第２の分割領域を設定するようにしてもよい。

１，５ 画像位置合わせ装置
２ 医用画像保管サーバ
３ 表示装置
４ 入力装置
６ ３次元画像
１０ 画像取得部
１１ 分割領域設定部
１２ 画素値変換方法決定部
１３ 画素値変換部
１４ 位置合わせ処理部
１５ パフュージョン解析部
１６ 表示制御部
５０ 画像取得部
５１ 分割領域設定部
５２ 位置合わせ処理部
５３ パフュージョン解析部
５４ 表示制御部
Ａ１ 右心房および右心室における血流方向を示す矢印
Ａ２ 左心房および左心室における血流方向を示す矢印
Ｄ１-Ｄ１０ 分割領域
ＩＲ 関心領域
Ｋ 定数
Ｌ１ 左心房の領域
Ｌ２ 左心室の領域
Ｍ 心筋の領域
Ｒ１ 右心房の領域
Ｒ２ 右心室の領域
Ｖ 大動脈の領域

Claims (17)

1. 位置合わせ処理対象を含む複数の画像を取得する画像取得部と、
   前記複数の画像の各位置合わせ処理対象を分割して分割領域を設定する分割領域設定部と、
   前記各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、該組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す位置合わせ処理部とを備え、
   前記位置合わせ処理部が、前記分割領域の組毎に、該組の各分割領域の画素値から前記組に含まれる２つの分割領域のコントラストが異なるか、または、前記２つの分割領域のコントラストが近いかを示す判定結果を得て、前記判定結果がコントラストが異なることを示す場合とコントラストが近いことを示す場合とで、前記評価関数を変更して前記位置合わせ処理を施すことを特徴とする画像位置合わせ装置。
2. 前記位置合わせ処理部は、前記組の各分割領域のうちの２つの分割領域の１つの分割領域の平均画素値が基準値以上であり、かつ、前記２つの分割領域のうち他の分割領域の平均画素値が前記基準値未満である場合の前記評価関数と、前記１つの分割領域と前記他の分割領域の双方の平均画素値が前記基準値以上、または、前記基準値未満である場合の前記評価関数を変更して前記位置合わせ処理を施す請求項１記載の画像位置合わせ装置。
3. 前記１つの分割領域の平均画素値が前記基準値以上であり、かつ、前記他の分割領域の平均画素値が前記基準値未満である場合の前記評価関数は、前記組の分割領域の位置が一致しているか否かを評価する関数であり、
   前記１つの分割領域と前記他の分割領域の双方の平均画素値が前記基準値以上、または、前記基準値未満である場合の前記評価関数は、前記組の分割領域の画素値が一致しているか否かを評価する関数である請求項２記載の画像位置合わせ装置。
4. 前記画像取得部が、前記位置合わせ処理対象として心臓画像を含む画像を取得し、
   前記分割領域設定部が、前記分割領域として、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域とを設定する請求項１〜３のいずれか１項に記載の画像位置合わせ装置。
5. 前記分割領域設定部が、前記複数の画像に含まれる大動脈の領域の画素値の変化に基づいて、前記第１の分割領域と前記第２の分割領域を設定する請求項４に記載の画像位置合わせ装置。
6. 位置合わせ処理対象である心臓画像を含む複数の画像を取得する画像取得部と、
   前記複数の画像の各位置合わせ処理対象の心臓画像を分割して、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域とを設定する際に、前記複数の画像に含まれる大動脈の領域の画素値の変化に基づいて、前記複数の画像の中から、前記左心房および左心室の領域の画素値よりも前記右心房および右心室の領域の画素値の方が高い第１の画像と前記右心房および右心室の領域の画素値よりも前記左心房および左心室の領域の画素値の方が高い第２の画像とを特定し、前記第１の画像を用いて前記第１の分割領域を設定し、前記第２の画像を用いて前記第２の分割領域を設定する分割領域設定部と、
   前記各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、該組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す位置合わせ処理部とを備え、
   前記位置合わせ処理部が、前記分割領域の組毎に、該組の各分割領域の画素値に基づいて前記評価関数を変更して前記位置合わせ処理を施すことを特徴とする画像位置合わせ装置。
7. 前記画像取得部が、前記位置合わせ処理対象として心臓画像を含む画像を取得し、
   前記分割領域設定部が、前記分割領域として、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域と心筋を含む第３の分割領域とを設定する請求項１〜３のいずれか１項に記載の画像位置合わせ装置。
8. 前記画像取得部が、前記位置合わせ処理対象として心臓画像を含む画像を取得し、
   前記分割領域設定部が、前記分割領域として、右心房を含む第１の分割領域と右心室を含む第２の分割領域と左心房を含む第３の分割領域と左心室を含む第４の分割領域とを設定する請求項１〜３のいずれか１項に記載の画像位置合わせ装置。
9. 前記画像取得部が、前記位置合わせ処理対象として心臓画像を含む画像を取得し、
   前記分割領域設定部が、前記分割領域として、右心房を含む第１の分割領域と右心室を含む第２の分割領域と左心房を含む第３の分割領域と左心室を含む第４の分割領域と心筋を含む第５の分割領域とを設定する請求項１〜３のいずれか１項に記載の画像位置合わせ装置。
10. 前記画像取得部が、前記位置合わせ処理対象として心臓画像を含む画像を取得し、
    前記分割領域設定部が、心臓の血流方向に垂直な断面によって前記心臓画像を分割して複数の分割領域を設定する請求項１〜３のいずれか１項に項記載の画像位置合わせ装置。
11. 前記分割領域設定部が、前記分割領域に対応する領域が予め設定された基準画像を用いて、前記複数の画像のうちの１つの画像に対して前記分割領域を設定し、該分割領域の設定結果を用いて、前記１つの画像以外の画像に対して分割領域を設定する請求項１から１０いずれか１項に記載の画像位置合わせ装置。
12. 前記画像取得部が、異なる時点で撮影された複数の画像を取得する請求項１から１１いずれか１項に記載の画像位置合わせ装置。
13. 前記画像取得部が、前記位置合わせ処理対象として心臓画像を含む画像を取得し、かつ複数の心拍における同じ心位相のタイミングで撮影された複数の画像を取得する請求項１から１２いずれか１項に記載の画像位置合わせ装置。
14. 位置合わせ処理対象を含む複数の画像を取得し、
    前記複数の画像の各位置合わせ処理対象を分割して分割領域を設定し、
    前記各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、該組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す際、前記分割領域の組毎に、該組の各分割領域の画素値から前記組に含まれる２つの分割領域のコントラストが異なるか、または、前記２つの分割領域のコントラストが近いかを示す判定結果を得て、前記判定結果がコントラストが異なることを示す場合とコントラストが近いことを示す場合とで、前記評価関数を変更して前記位置合わせ処理を施すことを特徴とする画像位置合わせ方法。
15. コンピュータを、位置合わせ処理対象を含む複数の画像を取得する画像取得部と、
    前記複数の画像の各位置合わせ処理対象を分割して分割領域を設定する分割領域設定部と、
    前記各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、該組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す位置合わせ処理部として機能させる画像位置合わせプログラムであって、
    前記位置合わせ処理部が、前記分割領域の組毎に、該組の各分割領域の画素値から前記組に含まれる２つの分割領域のコントラストが異なるか、または、前記２つの分割領域のコントラストが近いかを示す判定結果を得て、前記判定結果がコントラストが異なることを示す場合とコントラストが近いことを示す場合とで、前記評価関数を変更して前記位置合わせ処理を施すことを特徴とする画像位置合わせプログラム。
16. 位置合わせ処理対象である心臓画像を含む複数の画像を取得し、
    前記複数の画像の各位置合わせ処理対象の心臓画像を分割して、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域とを設定する際に、前記複数の画像に含まれる大動脈の領域の画素値の変化に基づいて、前記複数の画像の中から、前記左心房および左心室の領域の画素値よりも前記右心房および右心室の領域の画素値の方が高い第１の画像と前記右心房および右心室の領域の画素値よりも前記左心房および左心室の領域の画素値の方が高い第２の画像とを特定し、前記第１の画像を用いて前記第１の分割領域を設定し、前記第２の画像を用いて前記第２の分割領域を設定し、
    前記各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、該組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す際に、前記分割領域の組毎に、該組の各分割領域の画素値に基づいて前記評価関数を変更して前記位置合わせ処理を施すことを特徴とする画像位置合わせ方法。
17. コンピュータを、位置合わせ処理対象である心臓画像を含む複数の画像を取得する画像取得部と、
    前記複数の画像の各位置合わせ処理対象の心臓画像を分割して、右心房および右心室を含む第１の分割領域と左心房および左心室を含む第２の分割領域とを設定する際に、前記複数の画像に含まれる大動脈の領域の画素値の変化に基づいて、前記複数の画像の中から、前記左心房および左心室の領域の画素値よりも前記右心房および右心室の領域の画素値の方が高い第１の画像と前記右心房および右心室の領域の画素値よりも前記左心房および左心室の領域の画素値の方が高い第２の画像とを特定し、前記第１の画像を用いて前記第１の分割領域を設定し、前記第２の画像を用いて前記第２の分割領域を設定する分割領域設定部と、
    前記各位置合わせ処理対象の同じ位置に設定された分割領域の組毎に、該組の各分割領域の一致度を評価関数によって評価して位置合わせ処理を施す位置合わせ処理部として機能させる画像位置合わせプログラムであって、
    前記位置合わせ処理部が、前記分割領域の組毎に、該組の各分割領域の画素値に基づいて前記評価関数を変更して前記位置合わせ処理を施す画像位置合わせプログラム。