(57)【特許請求の範囲】
【請求項１】
時系列に撮影された複数の心臓画像を含む画像を取得する画像取得部と、
前記各画像の同一位置の画素値の変化を前記各画像の複数位置についてそれぞれ取得する画素値変化取得部と、
前記各画像の複数位置についてそれぞれ取得された前記画素値の変化を少なくとも心筋のクラス、右心房および右心室のクラス、並びに左心室および左心房のクラスにクラスタリングするクラスタリング部と、
前記各画像の各画素の前記クラスの情報および画素値に基づいて、前記各画像の心臓画像を複数の心臓の領域に分割する領域分割部と、
前記各画像の前記複数の心臓の領域の情報に基づいて、該各画像に対して位置合わせ処理を施す位置合わせ処理部とを備えたことを特徴とする画像位置合わせ装置。
【請求項２】
前記領域分割部が、前記右心房および右心室のクラスの画素値の統計値が、予め設定された閾値以下である場合には、該右心房および右心室のクラスを前記心筋のクラスに置き換え、前記右心房および右心室のクラスに対応する領域と前記心筋のクラスに対応する領域を同一領域として分割する請求項１に記載の画像位置合わせ装置。
【請求項３】
前記領域分割部が、前記左心房および左心室のクラスの画素値の統計値が、予め設定された閾値以下である場合には、該左心房および左心室のクラスを前記心筋のクラスに置き換え、前記左心房および左心室のクラスに対応する領域と前記心筋のクラスに対応する領域を同一領域として分割する請求項１または２に記載の画像位置合わせ装置。
【請求項４】
前記領域分割部が、前記各画像の各画素の前記クラスの情報および画素値を用いてグラフカットを行うことにより前記各画像を複数の領域に分割する請求項１から３いずれか１項に記載の画像位置合わせ装置。
【請求項５】
前記位置合わせ処理部が、前記各画像の複数の領域の情報に基づいて、該各画像の各領域の画素値を領域毎に異なる代表値に置き換え、該置き換えられた代表値に基づいて、前記位置合わせ処理を行う請求項１から４いずれか１項に記載の画像位置合わせ装置。
【請求項６】
前記領域分割部が、前記右心房および右心室のクラスの画素値の統計値が、予め設定された閾値以下である場合には、該右心房および右心室のクラスを前記心筋のクラスに置き換え、前記右心房および右心室のクラスに対応する領域と前記心筋のクラスに対応する領域を同一領域として分割し、
前記位置合わせ処理部が、前記位置合わせ処理を施す一方の画像の前記右心房および右心室のクラスに対応する領域が、前記心筋のクラスに対応する領域と同一領域として分割され、他方の画像の前記右心房および右心室のクラスに対応する領域と前記心筋のクラスに対応する領域とが別の領域として分割された場合には、前記他方の画像の前記右心房および右心室のクラスに対応する領域の画素値を前記心筋のクラスに対応する領域の代表値に置き換える請求項５に記載の画像位置合わせ装置。
【請求項７】
前記領域分割部が、前記左心房および左心室のクラスの画素値の統計値が、予め設定された閾値以下である場合には、該左心房および左心室のクラスを前記心筋のクラスに置き換え、前記左心房および左心室のクラスに対応する領域と前記心筋のクラスに対応する領域を同一領域として分割し、
前記位置合わせ処理部が、前記位置合わせ処理を施す一方の画像の前記左心房および左心室のクラスに対応する領域が、前記心筋のクラスに対応する領域と同一領域として分割され、他方の画像の前記左心房および左心室のクラスに対応する領域と前記心筋のクラスに対応する領域とが別の領域として分割された場合には、前記他方の画像の前記左心房および左心室のクラスに対応する領域の画素値を前記心筋のクラスに対応する領域の代表値に置き換える請求項５または６に記載の画像位置合わせ装置。
【請求項８】
前記クラスタリング部が、前記画素値の変化のピーク値および該ピーク値となる時間の少なくとも１つに基づいて、前記クラスタリングを行う請求項１から７いずれか１項に記載の画像位置合わせ装置。
【請求項９】
コンピュータが、
時系列に撮影された複数の心臓画像を含む画像を取得し、
前記各画像の同一位置の画素値の変化を前記各画像の複数位置についてそれぞれ取得し、
前記各画像の複数位置についてそれぞれ取得された前記画素値の変化を少なくとも心筋のクラス、右心房および右心室のクラス、並びに左心室および左心房のクラスにクラスタリングし、
前記各画像の各画素の前記クラスの情報および画素値に基づいて、前記各画像の心臓画像を複数の心臓の領域に分割し、
前記各画像の前記複数の心臓の領域の情報に基づいて、該各画像に対して位置合わせ処理を施すことを特徴とする画像位置合わせ方法。
【請求項１０】
コンピュータを、
時系列に撮影された複数の心臓画像を含む画像を取得する画像取得部と、
前記各画像の同一位置の画素値の変化を前記各画像の複数位置についてそれぞれ取得する画素値変化取得部と、
前記各画像の複数位置についてそれぞれ取得された前記画素値の変化を少なくとも心筋のクラス、右心房および右心室のクラス、並びに左心室および左心房のクラスにクラスタリングするクラスタリング部と、
前記各画像の各画素の前記クラスの情報および画素値に基づいて、前記各画像の心臓画像を複数の心臓の領域に分割する領域分割部と、
前記各画像の前記複数の心臓の領域の情報に基づいて、該各画像に対して位置合わせ処理を施す位置合わせ処理部として機能させることを特徴とする画像位置合わせプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の画像の位置合わせを行う画像位置合わせ装置および方法並びにプログラムに関するものである。
【背景技術】
【０００２】
従来、心筋、脳、肝臓および膵臓などの血流量を計測する手法として、パフュージョン撮影が行われている。
【０００３】
パフュージョン撮影とは、被検体に造影剤を注入し、ＣＴ（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）などがある。（たとえば非特許文献３参照）
【０００８】
【数１】
0006611660.tif 000002
なお、上式（１）における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に示すような心臓のアキシャル断層画像である場合、画素値変化取得部１１によって取得される画素値変化としては、大まかに図４に示すような右心房および右心室の領域の画素値変化ＲＶと、左心房および左心室の領域の画素値変化ＬＶと、心筋の領域の画素値変化ＶＶとに分類される。なお、図２Ｉ〜図２IIIにおける領域Ｒ１が右心房の領域であり、領域Ｒ２が右心室の領域であり、領域Ｌ１が左心房の領域であり、領域Ｌ２が左心室の領域であり、領域Ｍが心筋の領域である。図２Ｉについては、左心室と心筋との境界が明確でないので、これらの領域の境界は図示省略している。
【００４３】
クラスタリング部１２は、各画素位置の画素値変化を心筋の領域のクラスと、右心房および右心室の領域のクラスと、左心房および左心室の領域のクラスと、その他の領域のクラスとにクラスタリングする。クラスタリング部１２は、各画素位置の画素値変化のピーク値とそのピーク値となる時間を用いてクラスタリングを行う。具体的には、マハラノビス距離を用いたk-means法を用いてクラスタリングを行う。なお、クラスタリングの方法としては、k-means法に限らず、その他の公知なクラスタリング手法を用いることができる。
【００４４】
なお、各３次元画像の全ての画素位置の画素値変化を取得することなく、一部の画素位置の画素値変化を取得するようにした場合には、その画素値変化をクラスタリングした結果を用いて、上記一部の画素位置以外の画素に対してクラス情報を付与するようにすればよい。
【００４５】
領域分割部１３は、各３次元画像６の各画素のクラスの情報および画素値に基づいて、各３次元画像６を複数の領域に分割するものである。具体的には、本実施形態の領域分割部１３は、各３次元画像６を右心房および右心室の領域と、左心房および左心室の領域と、心筋の領域とに分割する。領域の分割方法としては、たとえばグラフカット法を用いることができる。
【００４６】
グラフカット法を用いる場合には、たとえば右心房および右心室の領域とその他の領域と分割する場合には、右心房および右心室の領域とその他の領域とにノードＳとノードＴが設定され、隣接する画素間の画素値の差に基づいて隣接する画素を繋ぐｎ−ｌｉｎｋが設定され、クラスの情報に基づいて各画素とノードＳまたはノードＴとをつなぐｔ−ｌｉｎｋが設定される。なお、左心房および左心室の領域とその他の領域とを分割する場合、および心筋の領域とその他の領域とを分割する場合も同様である。また、多クラスのグラフカット法を用いて右心房および右心室の領域と、左心房および左心室の領域と、心筋の領域とを同時に抽出するようにしてもよい。
【００４７】
ここで、たとえば図２Ｉに示す３次元画像６のように左心房および左心室の領域と心筋の領域との境界が明確でない場合がある。このような場合、グラフカット法を用いて領域分割した場合、高精度な領域分割ができない場合がある。そこで、本実施形態の領域分割部１３は、たとえば左心房および左心室のクラスに分類された領域内の平均画素値を算出し、その平均画素値が予め設定された閾値（たとえば１００）以下である場合には、その左心房および左心室のクラスの情報を心筋のクラスの情報に置き換え、左心房および左心室のクラスに対応する領域と心筋のクラスに対応する領域を同一の心筋の領域として分割する。また、右心房および右心室の領域も同様に、右心房および右心室のクラスに分類された領域内の平均画素値を算出し、その平均画素値が予め設定された閾値（たとえば１００）以下である場合には、その右心房および右心室のクラスの情報を心筋のクラスの情報に置き換え、右心房および右心室のクラスに対応する領域と心筋のクラスに対応する領域を同一の心筋の領域として分割する。
【００４８】
上述したように領域分割した場合、右心房および右心室の構造情報または左心房および左心室の構造情報がなくなることになるが、グラフカット法によって誤って領域分割された場合と比較すると、位置合わせ処理の精度を高くすることができる。
【００４９】
位置合わせ処理部１４は、領域分割部１３によって分割された領域の情報に基づいて、各３次元画像６に対して位置合わせ処理を施すものである。具体的には、各３次元画像６の各画素値を、分割された領域毎に予め設定された定数の代表値に置き換える。たとえば各３次元画像６を右心房および右心室の領域と、左心房および左心室の領域と、心筋の領域とに分割した場合には、各３次元画像６の右心房および右心室の領域の画素値を代表値「０」に置き換え、左心房および左心室の領域の画素値を代表値「４００」に置き換え、心筋の領域の画素値を「７００」に置き換え、その他の領域の画素値を代表値「−１００」に置き換える。
【００５０】
そして、上記のように画素値を定数に置き換えた画像間で位置合わせ処理を行う。位置合わせ処理としては、剛体位置合わせ処理または非剛体位置合わせ処理が施される。剛体位置合わせ処理としては、たとえばＩＣＰ（Iterative Closest Point）法を用いた処理などを用いることができるが、その他の公知な手法を用いるようにしてもよい。また、非剛体位置合わせ処理としては、たとえばＦＦＤ（Free-Form Deformation）法を用いた処理やＴＰＳ（Thin-Plate Spline）法を用いた処理などを用いることができるが、その他の公知な手法を用いるようにしてもよい。また、剛体位置合わせ処理を行った後に、非剛体位置合わせ処理を行うようにしてもよい。
【００５１】
そして、画素値を定数に置き換えた画像間の位置合わせ処理の際に算出された変形ベクトルを、画素値を置き換える前の３次元画像６に適用して３次元画像６を変形することによって、位置合わせ処理後の３次元画像６を生成する。
【００５２】
なお、位置合わせ処理を施す一方の画像が、上述したように右心房および右心室と心筋との境界が明確でなく、領域分割部１３において、右心房および右心室の領域と心筋の領域とが同一の領域として分割され、他方の画像が、右心房および右心室の領域と心筋の領域とが別の領域として分割されている場合には、位置合わせ処理部５２は、他方の３次元画像６の右心房および右心室の領域の画素値を心筋の領域の代表値に置き換える。これにより高精度に位置合わせ処理を行うことができる。
【００５３】
また、位置合わせ処理を施す一方の画像が、左心房および左心室と心筋との境界が明確でなく、領域分割部１３において、左心房および左心室の領域と心筋の領域とが同一の領域として分割され、他方の画像が、左心房および左心室の領域と心筋の領域とが別の領域として分割されている場合にも、位置合わせ処理部５２は、他方の３次元画像６の左心房および左心室の領域の画素値を心筋の領域の代表値に置き換える。
【００５４】
パフュージョン解析部１５は、位置合わせ処理部１４によって位置合わせ処理の施された２以上の３次元画像６を用いて、パフュージョン解析を行うものである。パフュージョン解析の手法としては、上述したようなMaximum Slope法またはDeconvolution法が用いられる。具体的には、図５Ｉに示すように、ユーザによって表示装置３に表示された３次元画像６上において関心領域ＩＲが指定される。パフュージョン解析部１５は、位置合わせ処理の施された２以上の３次元画像６のそれぞれの関心領域ＩＲ内の画像の平均濃度値を算出し、図５IIに示すような時間濃度曲線（Time Density Curve）を算出し、この時間濃度曲線に基づいてパフュージョン解析を行う。
【００５５】
表示制御部１６は、液晶ディスプレイなどの表示デバイスを備えたものであり、パフュージョン解析部１５による解析結果を表示装置３に表示させるものである。また、表示制御部１６は、画像取得部１０によって取得された２以上の３次元画像６および位置合わせ処理後の２以上の３次元画像６などを表示装置３に表示させるものである。
【００５６】
入力装置４は、ユーザによる種々の設定入力を受け付けるものであり、キーボードやマウスなどの入力デバイスを備えたものである。入力装置４は、たとえば患者の識別情報の設定入力、上述した関心領域ＩＲの設定入力、および画素値変化を取得する画素位置の設定入力などを受け付けるものである。
【００５７】
なお、タッチパネルを用いることによって表示装置３と入力装置４と兼用するようにしてもよい。
【００５８】
次に、本実施形態の医用画像診断支援システムの作用について、図６に示すフローチャートを参照しながら説明する。
【００５９】
まず、ユーザによる患者の識別情報などの入力に基づいて、その患者の心臓をパフュージョン撮影した２以上の３次元画像６が画像取得部１０によって取得される（Ｓ１０）。
【００６０】
画像取得部１０によって取得された２以上の３次元画像６は画素値変化取得部１１に入力される。画素値変化取得部１１は、各３次元画像６の各画素について、画素値の時間変化を取得する（Ｓ１２）。
【００６１】
次に、各３次元画像６の各画素位置の画素値変化がクラスタリング部１２に入力され、クラスタリング部１２は、各画素位置の画素値変化を心筋の領域のクラスと、右心房および右心室の領域のクラスと、左心房および左心室の領域のクラスと、その他の領域のクラスとにクラスタリングする（Ｓ１４）。
【００６２】
そして、各３次元画像６の各画素のクラスの情報および画素値に基づいて、領域分割部１３によって各３次元画像６が複数の領域に分割される（Ｓ１６）。具体的には、各３次元画像６は、右心房および右心室の領域と、左心房および左心室の領域と、心筋の領域とに分割される。
【００６３】
次いで、位置合わせ処理部５２によって、各３次元画像６の各画素値が、分割された領域毎に予め設定された代表値に置き換えられ、その代表値に置き換えられた画像間で位置合わせ処理が施され、その変形ベクトルを用いて各３次元画像６を変形することによって、位置合わせ処理後の３次元画像６を生成する（Ｓ１８）。
【００６４】
位置合わせ処理の施された３次元画像６はパフュージョン解析部１５に入力され、パフュージョン解析部１５において、入力された３次元画像６を用いてパフュージョン解析が行われる（Ｓ２０）。
【００６５】
パフュージョン解析部１５によるパフュージョン解析結果は表示制御部１６に入力され、表示制御部１６は、入力されたパフュージョン解析結果を表示装置３に表示させる（Ｓ２２）。
【００６６】
上記実施形態の医用画像診断支援システムによれば、時系列に撮影された複数の画像を取得し、各画像の同一位置の画素値の変化を各画像の複数位置についてそれぞれ取得する。そして、各画像の複数位置についてそれぞれ取得された画素値の変化を複数のクラスにクラスタリングし、各画像の各画素のクラスの情報および画素値に基づいて、各画像を複数の領域に分割し、各画像の複数の領域の情報に基づいて、各画像に対して位置合わせ処理を施す。
【００６７】
このように各画像における各画素の画素値の変化をクラスタリングした結果を用いて各画像を複数の領域に分割し、その分割された領域の情報を用いて位置合わせ処理を行うことによって、一方の画像では高い画素値を有する構造が、他方の画像では低い画素値を示している場合でも、そのコントラスト差を低減することができ、高精度に位置合わせを行うことができる。
【００６８】
なお、上記実施形態においては、位置合わせ処理部５２において、各３次元画像６の画素値を領域毎の代表値に置き換えて位置合わせ処理を行うようにしたが、位置合わせ処理の方法としては、必ずしもこの手法に限らない。
【００６９】
たとえば領域分割部１３において、各領域に対してラベルを付与してラベル画像を生成する。具体的には、たとえば各３次元画像６の右心房および右心室の領域に対してラベル「１」を付与し、左心房および左心室の領域に対してラベル「２」を付与し、心筋の領域にラベル「３」を付与し、その他の領域にラベル「４」を付してラベル画像を生成し、そのラベル画像を領域分割部１３から位置合わせ処理部１４に出力させる。
【００７０】
そして、位置合わせ処理部１４は、下式（２）の評価関数を用いて一致度を算出する。
【００７１】
【数２】
0006611660.tif 000003
上式（２）におけるSは一致度の評価関数であり、Nは画像間(Fixed画像とMoving画像の間)で比較する画素の総数、xiは座標点、ＬF(xi)はFixed画像のxi におけるラベル、ＬM(xi)はMoving画像のxi におけるラベル、Tは座標点の幾何学的変換関数、μは幾何学的変換のパラメタである。
【００７２】
また、上式（２）における関数Ｃは、下式（３）である。関数Ｃは、ｌ１とｌ２のラベルが一致している場合には定数Ｋを出力し、不一致の場合には０（ゼロ）を出力する。関数Ｃをこのように定義した場合、Sが大きいほど一致していることになる。
【００７３】
【数３】
0006611660.tif 000004
位置合わせ処理部１４は、一致度Sが高くなるような変形パラメタμを求め、その変形パラメタμを用いて各３次元画像６を変形することによって、位置合わせ処理後の３次元画像６を生成する。
【符号の説明】
【００７４】
１     画像位置合わせ装置
２     医用画像保管サーバ
３     表示装置
４     入力装置
６     ３次元画像
１０   画像取得部
１１   画素値変化取得部
１２   クラスタリング部
１３   領域分割部
１４   位置合わせ処理部
１５   パフュージョン解析部
１６   表示制御部
ＩＲ   関心領域
Ｌ１   左心房の領域
Ｌ２   左心室の領域
ＬＶ   左心房および左心室の領域の画素値変化
Ｍ     心筋の領域
Ｒ１   右心房の領域
Ｒ２   右心室の領域
ＲＶ   右心房および右心室の領域の画素値変化
ＶＶ   心筋の領域の画素値変化
【図１】
0006611660.tif 000005
【図２】
0006611660.tif 000006
【図３】
0006611660.tif 000007
【図４】
0006611660.tif 000008
【図５】
0006611660.tif 000009
【図６】
0006611660.tif 000010