(57)【要約】【課題】実際に撮影した画像が３相画像のように少ない画像であっても、組織の血流変化モデルを適用して高精度かつ高速な血流動態解析を行うこと。【解決手段】複数患者の医用画像に基づき造影血流の流路となる血管及び臓器の位置に設定した関心領域の画素値を測定し、所定の血管における基準正規化信号強度曲線と、この曲線と所定の血管との関係に関するパラメータを保存するデータベース４１と、この曲線に基づき定めた複数のシミュレーション点に対応付けて薬物動態解析モデルの係数を保存したデータベース４２と、３相撮影された対象患者の医用画像に基づき、造影血流の流路となる血管及び臓器の位置における関心領域の早期相と後期相の造影強度の組み合わせに対応するシミュレーション点を検索する検索部１３と、シミュレーション点に対応付けられた係数を取得する取得部１４と、係数に基づき所定の血管の時間造影強度曲線を推定する推定部１５とを備える。【選択図】図２
2020168210.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも３相撮影された患者の医用画像に基づいて血流動態解析を行う際に、薬物動態解析モデルの入力となる血管の信号強度曲線を何らかの手法で推定する時間造影強度曲線推定部が動作した場合に、解析対象の臓器に関連する部位の信号強度に基づいて定めた早期相と後期相（または平衡相）を選択する相選択部と、選択された早期相と後期相の組み合わせからなる複数のシミュレーション点に対応付けて、最適化計算に基づいて算出した薬物動態解析モデルの最適係数を算出し、薬物動態解析モデル係数データベースとして保存するシミュレーション部と、少なくとも３相撮影された画像の対応する位置の画素において早期相と後期相の組み合わせに対応する前記シミュレーション点を検索する検索部と、前記検索部により検索した前記シミュレーション点に対応付けられた前記係数を取得する取得部と、前記検索部により検索した前記係数に基づいて、前記対象患者における所定の画素の時間造影強度曲線を推定する時間造影強度曲線推定部と、を備える、ことを特徴とする血流動態解析システム。
【請求項２】
前記時間造影強度曲線推定部は、少なくとも３相撮影された患者の医用画像において、薬物動態モデルの入力となる血管位置に設定した関心領域の画素値から推定されるものである、
ことを特徴とする血流動態解析システム。
【請求項３】
前記時間造影強度曲線推定部は、多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースと、
少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定する対象患者血管信号値測定部と、
前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線と、前記基準正規化信号強度曲線と前記対象患者の血管との関係に関するパラメータと、前記対象患者血管信号値測定部で得られた対象患者の血管信号値とから、前記対象患者の前記所定の血管の造影剤の到達距離の指標と循環速度の指標を計算し、前記所定の血管の時間造影強度曲線を推定する時間造影強度曲線推定部、から推定されるものである、
ことを特徴とする血流動態解析システム。
【請求項４】
前記相選択部において、３相の場合は、撮影時刻に基づきそのまま早期相、後期相を選択する、
ことを特徴とする血流動態解析システム。
【請求項５】
前記相選択部において、多相の場合は、解析対象の臓器が肝臓の場合、関連する部位として脾臓を用いて相選択を行う、
ことを特徴とする血流動態解析システム。
【請求項６】
前記シミュレーション部は、
脾臓のコンパートメントモデル式を用いて、動脈に関する時定数パラメータを決定する第１決定部と、
肝臓部のコンパートメントモデル式を用いて、門脈に関する時定数パラメータを決定する第２決定部と、
前記肝臓のコンパートメントモデル式における残りの係数を得るために、造影剤の各血管からの流入速度および注目組織からの流出速度を計算する計算部と、を備え、
前記複数のシミュレーション点に対する前記薬物動態解析モデルの最適な係数の算出を行う、
ことを特徴とする請求項１に記載の血流動態解析システム。
【請求項７】
前記シミュレーション部は、
前記基準正規化信号強度曲線を、大動脈および門脈の全時間造影曲線に基づいて算出し、前記複数のシミュレーション点を、前記大動脈および門脈の全時間造影曲線に基づいて定める、
ことを特徴とする請求項１または請求項２に記載の血流動態解析システム。
【請求項８】
前記シミュレーション部は、
前記薬物動態解析モデル係数データベースの作成における前記早期相の決定を、脾臓の信号強度を用いて行い、脾臓の信号強度の１秒当たりの増加率が最大の相の次の相を早期相とする、
ことを特徴とする請求項１ないし請求項３のいずれか１項に記載の血流動態解析システム。
【請求項９】
前記時間造影強度曲線推定部により推定した前記対象患者における所定の血管の時間造影強度曲線に基づいて、任意の時刻の仮想ダイナミック像を生成する仮想ダイナミック像生成部を備える、
ことを特徴とする請求項１ないし請求項４のいずれか１項に記載の血流動態解析システム。
【請求項１０】
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースを作成するステップと、
シミュレーション部により、前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出するステップと、
前記シミュレーション部により、前記複数のシミュレーション点と、前記薬物動態解析モデルの最適な係数とを対応付けて薬物動態解析モデル係数データベースを作成するステップと、
対象患者造影強度測定部により、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の早期相と後期相の造影強度を測定するステップと、
検索部により、前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索するステップと、
取得部により、前記検索部にて検索した前記シミュレーション点に対応付けられた前記係数を取得するステップと、
間造影強度曲線推定部により、前記取得部にて取得した前記係数に基づいて、前記所定の血管の時間造影強度曲線を推定するステップと、を備える、
ことを特徴とする血流動態解析方法。
【請求項１１】
血流動態解析装置のプログラムであって、前記プログラムは、コンピュータに、
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースを作成するステップと、
シミュレーション部により、前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出するステップと、
前記シミュレーション部により、前記複数のシミュレーション点と、前記薬物動態解析モデルの最適な係数とを対応付けて薬物動態解析モデル係数データベースを作成するステップと、
対象患者造影強度測定部により、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の早期相と後期相の造影強度を測定するステップと、
検索部により、前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索するステップと、
取得部により、前記検索部にて検索した前記シミュレーション点に対応付けられた前記係数を取得するステップと、
間造影強度曲線推定部により、前記取得部にて取得した前記係数に基づいて、前記所定の血管の時間造影強度曲線を推定するステップと、を実行させる、
ことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、血流動態解析システム、血流動態解析方法、およびプログラムに関する。
【背景技術】
【０００２】
肝がん等の診断においては、造影剤を注入して経過観察するダイナミックＣＴ撮影が利用されている。例えば、特許文献１においては、ダイナミックＣＴ撮影による多相画像の主要な入力血管領域に設定した関心領域上の画素値（代表値）から、患者の主血管の血流動態を推定することにより、診断に役立つ情報の提供を行っている。特許文献１においては、推定した血管の血流動態を入力波形として、ピクセルごとに、あるいは、設定した関心領域で、組織の血流変化モデルを適用して、最適化計算により血流動態を推定している。高精度な血流動態解析を行うためには１０相撮影などの大量の画像が必要とされているが、特許文献２では、データ数が少ないとされる多くの病院の日常の検査で撮影されている３相画像でも、同等精度の高精度な血流動態解析を行うことができる方法を示した。このような血流動態解析を対象画像の全画素に対して行うことで、任意時刻の仮想ダイナミック像や血流動態パラメータの分布図(パラメータマップ)を作成することができる。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−２３０７９３号公報
【特許文献２】特許６４９１３９１
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１や特許文献２のように対象画像の全画素すべてに対し血流動態解析を行い、任意時刻の仮想ダイナミック像や血流動態パラメータの分布図を計算しようとすると、対象画像のサイズによって、計算量が増大してしまい、待ち時間が大きくなる問題があった。
【０００５】
そこで、この発明の課題は、画像に対する血流動態解析及び任意時刻の仮想ダイナミック像やパラメータマップ作成の高速化を実現できる血流動態解析システム、血流動態解析方法、およびプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
本発明の血流動態解析システムの例として、まず特許文献２を拡張した態様を説明する。
本発明の血流動態解析システムの一態様は、上記課題を解決するため、
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースと、
前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と、後期相（または平衡相）の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出し、薬物動態解析モデル係数データベースとして保存するシミュレーション部と、
少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の造影強度を測定する対象患者造影強度測定部と、
前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索する検索部と、
前記検索部により検索した前記シミュレーション点に対応付けられた前記係数を取得する取得部と、
前記検索部により検索した前記係数に基づいて、前記対象患者における所定の血管の時間造影強度曲線(Time Density Curve:ＴＤＣ)を推定する時間造影強度曲線推定部と、を備える。
【０００７】
本態様は、造影血流の流路となる主要血管を臨床的な事前知識として利用する。その主要血管の位置に関心領域を設定し、画素値（代表値）の時間変化を測定することで、複数の流路間の造影血流循環にかかる距離や時間の関係を考慮して、関心領域における造影血流動態を推定する。複数患者血管信号値データベースは、多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存したデータベースである。また、複数患者血管信号値データベースは、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する。シミュレーション部は、複数患者血管信号値データベースに保存された基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出し、薬物動態解析モデル係数データベースとして保存する。対象患者血管信号値測定部は、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定する。検索部は、関心領域の早期相と後期相の造影強度の組み合わせに対応するシミュレーション点を検索する。取得部は、検索部により検索したシミュレーション点に対応付けられた係数を取得する。時間造影強度曲線推定部は、検索部により検索した係数に基づいて、対象患者における所定の血管の時間造影強度曲線(Time Density Curve:ＴＤＣ)を推定する。
【０００８】
本発明の血流動態解析方法の一態様は、上記課題を解決するため、
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースを作成するステップと、
シミュレーション部により、前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出するステップと、
前記シミュレーション部により、前記複数のシミュレーション点と、前記薬物動態解析モデルの最適な係数とを対応付けて薬物動態解析モデル係数データベースを作成するステップと、
対象患者造影強度測定部により、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の早期相と後期相の造影強度を測定するステップと、
検索部により、前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索するステップと、
取得部により、前記検索部にて検索した前記シミュレーション点に対応付けられた前記係数を取得するステップと、
間造影強度曲線推定部により、前記取得部にて取得した前記係数に基づいて、前記所定の血管の時間造影強度曲線を推定するステップと、を備える。
【０００９】
本態様は、造影血流の流路となる主要血管を臨床的な事前知識として利用する。その主要血管の位置に関心領域を設定し、画素値（代表値）の時間変化を測定することで、複数の流路間の造影血流循環にかかる距離や時間の関係を考慮して、関心領域における造影血流動態を推定する。複数患者血管信号値データベースは、多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存したデータベースである。また、複数患者血管信号値データベースは、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する。シミュレーション部は、複数患者血管信号値データベースに保存された基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出し、薬物動態解析モデル係数データベースとして保存する。対象患者血管信号値測定部は、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定する。検索部は、関心領域の早期相と後期相の造影強度の組み合わせに対応するシミュレーション点を検索する。取得部は、検索部により検索したシミュレーション点に対応付けられた係数を取得する。時間造影強度曲線推定部は、検索部により検索した係数に基づいて、対象患者における所定の血管の時間造影強度曲線(Time Density Curve:ＴＤＣ)を推定する。
【００１０】
本発明の血流動態解析装置のプログラムの一態様は、上記課題を解決するため、
血流動態解析装置のプログラムであって、前記プログラムは、コンピュータに、
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースを作成するステップと、
シミュレーション部により、前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出するステップと、
前記シミュレーション部により、前記複数のシミュレーション点と、前記薬物動態解析モデルの最適な係数とを対応付けて薬物動態解析モデル係数データベースを作成するステップと、
対象患者造影強度測定部により、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の早期相と後期相の造影強度を測定するステップと、
検索部により、前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索するステップと、
取得部により、前記検索部にて検索した前記シミュレーション点に対応付けられた前記係数を取得するステップと、
間造影強度曲線推定部により、前記取得部にて取得した前記係数に基づいて、前記所定の血管の時間造影強度曲線を推定するステップと、を実行させる。
【００１１】
本態様は、造影血流の流路となる主要血管を臨床的な事前知識として利用する。その主要血管の位置に関心領域を設定し、画素値（代表値）の時間変化を測定することで、複数の流路間の造影血流循環にかかる距離や時間の関係を考慮して、関心領域における造影血流動態を推定する。複数患者血管信号値データベースは、多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存したデータベースである。また、複数患者血管信号値データベースは、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する。シミュレーション部は、複数患者血管信号値データベースに保存された基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出し、薬物動態解析モデル係数データベースとして保存する。対象患者血管信号値測定部は、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定する。検索部は、関心領域の早期相と後期相の造影強度の組み合わせに対応するシミュレーション点を検索する。取得部は、検索部により検索したシミュレーション点に対応付けられた係数を取得する。時間造影強度曲線推定部は、検索部により検索した係数に基づいて、対象患者における所定の血管の時間造影強度曲線(Time Density Curve:ＴＤＣ)を推定する。
【発明の効果】
【００１２】
以上より明らかなように、本発明の血流動態解析システム、血流動態解析方法、およびプログラムによれば、実際に撮影した画像が３相画像のように少ない画像であっても、血流動態パラメータの計算を行い、解析対象臓器の入力血管の時間造影強度曲線を推定し、組織の血流変化モデルを適用して高精度かつ高速に血流動態解析を行うことができる。
【図面の簡単な説明】
【００１３】
【図１】第１実施形態における肝血流動態解析システムを示す概略図である。
【図２】（Ａ）は肝血流動態解析装置の構成を示すブロック図、（Ｂ）はビューワの構成を示すブロック図である。
【図３】肝血流動態解析システムの処理の概要を説明するための図である。
【図４】３相画像における肝臓の流入血管である腹部大動脈と門脈の造影信号強度の一例を示す図である。
【図５】腹部大動脈と門脈について推定した時間造影強度曲線の一例を示す図である。
【図６】造影剤の血流路における経路を説明するための図である。
【図７】ＬＡ（左心房）に設定した注目領域における１０時相の平均ＣＴ値の一例を示す図である。
【図８】造影強度曲線の一例を示す図である。
【図９】ＡＡＯ（腹部大動脈）における造影血流循環の速さ(体内時計)を正規化したデータの一例を示す図である。
【図１０】ＲＲＶ（右腎静脈）における造影血流循環の速さ(体内時計)を正規化したデータの一例を示す図である。
【図１１】ＤＰＩをプロットした図である。
【図１２】正規化信号強度曲線を示す図である。
【図１３】ＤＰＩを早期相帯と後期相帯とに分離した図である。
【図１４】（Ａ）は、早期相について正規化信号強度をＤＰＩで多項式フィッティングした例を示す図、（Ｂ）は、後期相について正規化信号強度をＤＰＩで多項式フィッティングした例を示す図である。
【図１５】（Ａ）は６０人の患者のうちの５人についてのＤＳＩを示す図、（Ｂ）５人についてのＤＳＩをプロットした図である。
【図１６】時刻ＴとＤＳＩ時刻（Ｔｍ＝ＤＳＩ＊Ｔ）の関係を示す図である。
【図１７】ＤＳＩ時刻（Ｔｍ）を用いてＤＰＩを多項式フィッティングした例を示す図である。
【図１８】ＬＡ（左心房），ＡＡＯ（腹部大動脈），ＲＲＶ（右腎静脈），ＳＰＶ（脾静脈），ＳＭＶ（上腸間膜静脈），ＩＶＣ（下大静脈），ＰＶ（門脈）に０秒、４０秒、１３０秒における画素値の例を示す図である。
【図１９】図１８に示す画素値をプロットした図である。
【図２０】サンプリング時刻ＴとＤＳＩ時刻ｍＴとの関係を示す図である。
【図２１】サンプリング時刻ＴとＤＳＩ時刻ｍＴとの関係を示す図に平衡相開始時刻ｅｑ_ｌｉｍｉｔを示した図である。
【図２２】サンプリング時刻Ｔと平衡相時間帯ｅｑＰｈａｓｅとの関係を示す図である。
【図２３】ＤＰＩの一例を示す図である。
【図２４】図２３のようなＤＰＩが得られた場合にＤＳＩ時刻に対応するＤＰＩを計算した一例を示す図である。
【図２５】時刻と補正後のＤＰＩとの関係を示す図である。
【図２６】図２５にＤＰＩピーク値の学習データ数分の平均値ｐｅａｋＤＰＩの直線を示した図である。
【図２７】図２６におけるＤＰＩの経時的変化を示す曲線と、平均値ｐｅａｋＤＰＩの直線との交点に基づいて、早期相、後期相、および平衡相を分離した例を示す図である。
【図２８】図２７に示す早期相で計算した正規化信号強度曲線の例を示す図である。
【図２９】図２７に示す後期相で計算した正規化信号強度曲線の例を示す図である。
【図３０】図２８に示す正規化信号強度曲線と、図２９に示す正規化信号強度曲線とを連結して作成した早期相、後期相、および平衡相の正規化信号強度曲線を示す図である。
【図３１】平均値ｐｅａｋＤＰＩに基づいて、早期相、後期相、および平衡相を分離した図である。
【図３２】ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅを示す図である。
【図３３】ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２から計算を行った正規化信号強度曲線の例を示す図である。
【図３４】早期相の正規化信号強度曲線と、後期相の正規化信号強度曲線とを連結した正規化信号強度曲線を示す図である。
【図３５】サンプリング時刻Ｔに対する早期相と後期相の正規化信号強度曲線、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２から計算を行った正規化信号強度曲線、および平均値ｐｅａｋＤＰＩに基づいて分離した平衡相の開始時刻ｔ１から、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２までの正規化信号強度を連結した図である。
【図３６】早期相の対象血管の造影信号強度ｔ_ｄＸ（３）の対象血管の造影信号強度を示す図である。
【図３７】患者番号Ｉの患者の９０秒における学習データの対象血管の正規化信号強度を示す図である。
【図３８】全時間的な正規化信号強度曲線に推定した対象血管の時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅを適用して得られた時間造影強度曲線の一例を示す図である。
【図３９】平衡相帯を実際の後期相の観測値を用いて補正した時間造影強度曲線の一例を示す図である。
【図４０】平均値ｐｅａｋＤＰＩに基づいて分離した平衡相の開始時刻ｔ１から、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２までの未計算の部分を線形補間した一例を示す図である。
【図４１】未計算の部分の補間を行った後の時間造影強度曲線を示す図である。
【図４２】４１秒で時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅに達した例を示す図である。
【図４３】平衡相の直線を示す図である。
【図４４】補正を行った後の時間造影強度曲線を示す図である。
【図４５】テストフェイズで計算した腹部大動脈(脾動脈の代替)と、血流動態解析の結果である脾臓の時間造影強度曲線を示す図である。
【図４６】テストフェイズで計算した腹部大動脈(肝動脈の代替)及び門脈と、血流動態解析の結果である肝臓の時間造影強度曲線を示す図である。
【図４７】動脈と門脈との全時間造影強度曲線の一例を示す図である。
【図４８】造影強度パターンのシミュレーション点の一例を示す図である。
【図４９】コンパートメントモデル式における係数の最適化の概略手順についてのフローチャートである。
【図５０】コンパートメントモデル式の各係数にサンプリング点(ｐ点)分の値をランダムに割り当てた一例を示す表である。
【図５１】（Ａ）は図５０に示すサンプリング点０における各係数値で肝臓のコンパートメントモデル式を数値解析で解き、時間造影強度曲線を計算した結果を示す図、（Ｂ）は図５０に示すサンプリング点１における各係数値で、上述した肝臓のコンパートメントモデル式を数値解析で解き、時間造影強度曲線を計算した結果を示す図、（Ｃ）は図５０に示すサンプリング点ｐ−１における各係数値で、上述した肝臓のコンパートメントモデル式を数値解析で解き、時間造影強度曲線を計算した結果を示す図である。
【図５２】（Ａ），（Ｂ），（Ｃ）は、各観測値と推定曲線上の対応値との距離を示図である。
【図５３】薬物動態解析モデル係数データベースの検索処理を説明するための図である。
【図５４】ビューワにおける入力画面の一例を示す図である。
【図５５】選択されたスタックおよび部位ＲＯＩ指定ツールの拡大図である。
【図５６】大動脈のボタンをクリックした後にスタック上で当該部位の位置Ａをクリックした例を示す図である。
【図５７】ビューワにおける出力画面の一例を示す図である。
【図５８】（Ａ），（Ｂ）は選択アイコンの一例を示す図である。
【図５９】ビューワにおける出力画面の一例を示す図である。
【図６０】（Ａ）は選択されたスタックのオーバーレイ表示の一例を示す図、（Ｂ）は縦軸を造影強度、横軸を経過時刻として表した造影強度曲線における大動脈のピーク造影強度の時刻を示す図、（Ｃ）は縦軸を造影強度、横軸を経過時刻として表した造影強度曲線における門脈のピーク造影強度の時刻を示す図である。
【図６１】本発明の第２実施形態における大動脈と門脈の全時間造影強度曲線における早期相および平衡相と、コンパートモデル式のパラメータとの関係を説明するための図である。
【図６２】脾臓の全時間造影強度曲線の一例を示す図である。
【図６３】脾臓の全時間造影強度曲線における傾きを説明する図である。
【発明を実施するための形態】
【００１４】
以下、この発明の実施の形態を、図面を参照しながら詳細に説明する。
【００１５】
（第１実施形態）
まず、本発明の第１実施形態について説明する。
【００１６】
＜肝血流動態解析システムの概要＞
図１は、本発明の第１実施形態における肝血流動態解析システム１の概要を示すである。図１に示すように、本実施形態における肝血流動態解析システム１は、肝血流動態解析装置１００と、ビューワ２００と、画像サーバ３００とを備えている。肝血流動態解析装置１００と、ビューワ２００と、画像サーバ３００とは、ネットワークＮを介して通信可能に接続されている。
【００１７】
肝血流動態解析装置１００およびビューワ２００は、いわゆるコンピュータであり、それぞれ内部に演算制御部（ＣＰＵ：Central Processing Unit）、主記憶部（ＲＡＭ：Random Access Memory）、外部記憶部（ＨＤＤ：Hard Disk Drive）、通信コントローラを備えている。
【００１８】
肝血流動態解析装置１００は、コンピュータにキーボードやマウス等の入力デバイスと、液晶ディスプレイ等のモニタ（表示部）とを接続して構成されている。それぞれの外部記憶部には、ＯＳ（オペレーティングシステム）およびアプリケーションのプログラムが記憶されている。通信コントローラは、ネットワークＮの接続線と接続され、例えばＷＷＷ（ワールド・ワイド・ウェブ）やＤＩＣＯＭ（Digital Imaging and COmmunications in Medicine）通信プロトコル、ＴＣＰ／ＩＰプロトコル等のデータ通信制御によりネットワークＮを介して行なうデータ通信を制御する。
【００１９】
画像サーバ３００は、図示を省略する画像診断装置とネットワークＮを介して接続されており、画像診断装置で作成されたＣＴ画像等の医用画像を画像データベースに登録して保管する。画像診断装置は、患者を撮影して、その撮影結果を反映したＣＴ画像等の医用画像作成し、作成した医用画像を画像サーバ３００に送信する。
【００２０】
ビューワ２００は、ユーザによる操作に基づいて、ＣＴ画像を画像サーバ３００から取得する。また、ビューワ２００は、ユーザによる操作に基づいて、肝血流動態解析装置１００によって生成された仮想ダイナミック像等を取得する。
【００２１】
なお、本実施形態では、肝血流動態解析装置１００とビューワ２００とを別体の装置としているが、肝血流動態解析装置１００とビューワ２００とを兼ねる構成であっても良い。
【００２２】
［肝血流動態解析装置の機能ブロック］
図２（Ａ）は、肝血流動態解析装置１００の機能ブロックを表すブロック図である。図２（Ａ）に示すように、肝血流動態解析装置１００は、制御部１０と、記憶部４０と、通信処理部５０とを備えている。制御部１０は、対象患者造影強度測定部１１と、シミュレーション部１２と、検索部１３と、取得部１４と、時間造影強度推定部１５と、仮想ダイナミック像生成部１６としての機能を備えている。
【００２３】
対象患者造影強度測定部１１は、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の造影強度を測定する。対象患者造影強度測定部１１の機能の詳細については後述する。
【００２４】
シミュレーション部１２は、複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出し、後述する薬物動態解析モデル係数データベース４２として保存する。シミュレーション部１２は、後述する複数患者血管信号値データベース４１に保存された基準正規化信号強度曲線に基づいて早期相と後期相を定め、早期相と後期相の信号強度の組み合わせからなるシミュレーション点を定義する。
【００２５】
また、シミュレーション部１２は、第１決定部１２ａ、第２決定部１２ｂ、および計算部１２ｃとしての機能を備える。第１決定部１２ａ、第２決定部１２ｂ、および計算部１２ｃは、複数のシミュレーション点に対する薬物動態解析モデルの最適な係数の算出を行うために用いられる。第１決定部１２ａは、脾臓のコンパートメントモデル式を用いて、動脈に関する時定数パラメータを決定する。第２決定部１２ｂは、肝臓部のコンパートメントモデル式を用いて、門脈に関する時定数パラメータを決定する。計算部１２ｃは、肝臓のコンパートメントモデル式における残りの係数を得るために、造影剤の各血管からの流入速度および注目組織からの流出速度を計算する。
【００２６】
シミュレーション部１２、第１決定部１２ａ、第２決定部１２ｂ、および計算部１２ｃの機能の詳細については後述する。
【００２７】
検索部１３は、記憶部４０における薬物動態解析モデル係数データベース４２における複数のシミュレーション点のうち、前記関心領域の早期相と後期相の造影強度の組み合わせに対応するシミュレーション点を検索する。検索部１３の機能の詳細については後述する。
【００２８】
取得部１４は、検索部１３により検索したシミュレーション点に対応付けられた、薬物動態解析モデルの最適な係数を取得する。取得部１４の機能の詳細については後述する。
【００２９】
時間造影強度曲線推定部１５は、取得部１４により取得した薬物動態解析モデルの最適な係数に基づいて、対象患者における所定の血管の時間造影強度曲線を推定する。時間造影強度曲線推定部１５の機能の詳細については後述する。
【００３０】
仮想ダイナミック像生成部１６は、対象患者における所定の組織の時間造影強度曲線に基づいて、任意の時刻の仮想ダイナミック像を生成する。また、仮想ダイナミック像生成部１３は、時間変化を保持した仮想ダイナミック動画像データを生成するようにしてもよい。仮想ダイナミック像生成部１６の機能の詳細については後述する。
【００３１】
記憶部４０は、対象患者造影強度測定部１１、シミュレーション部１２、時間造影強度曲線推定部１５、および仮想ダイナミック像生成部１６によって算出、測定、推定、解析、または生成されたデータ等の肝血流動態解析処理に必要なデータを記憶する。また、記憶部４０は、複数患者血管信号値データベース４１、および薬物動態解析データベース４２の記憶部としても用いられる。
【００３２】
複数患者血管信号値データベース４１を作成するには、多相撮影された複数患者のＣＴ画像データを含むデータを利用して、所定の部位における血管の時刻と基準正規化信号強度曲線との関係に関するパラメータを計算する。そして、複数患者血管信号値データベース４１は、このように計算したパラメータを保存する。複数患者血管信号値データベース４１の詳細については後述する。
【００３３】
薬物動態解析データベース４２を作成するには、複数患者血管信号値データベース４１に保存された基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点を決定する。次に、複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を保存する。薬物動態解析データベース４２の作成方法の詳細については後述する。薬物動態解析モデルの一例としては、コンパートメントモデルが挙げられる。
【００３４】
通信処理部５０は、ネットワークＮを介して、ビューワ２００および画像サーバ３００との通信を行う。
【００３５】
［ビューワの機能ブロック］
図２（Ｂ）は、ビューワ２００の機能ブロックを表すブロック図である。図２（Ｂ）に示すように、ビューワ２００は、制御部２１０と、表示部２１１と、通信処理部２１２と、記憶部２１３とを備えている。制御部２１０は、表示部２１１、通信処理部２１２、および記憶部２１３の制御を行う。
【００３６】
通信処理部２１２は、ネットワークＮを介して、肝血流動態解析装置１００および画像サーバ３００との通信を行う。
【００３７】
表示部２１１は、例えば液晶ディスプレイまたは有機ＥＬディスプレイ等により構成され、肝血流動態解析装置１００の仮想ダイナミック像生成部１３により生成された仮想ダイナミック像を表示する。後述するように、表示部２１１は、任意の時刻を指定するユーザインターフェースを備えており、ユーザインターフェースにより指定された任意の時刻における仮想ダイナミック像を表示する。表示部２１１における表示例等については後述する。
【００３８】
＜肝血流動態解析システム１の処理の概要＞
次に、肝血流動態解析システム１の処理の概要について、図面を参照しつつ説明する。図３は、肝血流動態解析システム１の処理の概要を説明するための図である。図４は、３相画像における肝臓の流入血管である腹部大動脈と門脈の造影信号強度の一例を示す図である。図５は、腹部大動脈と門脈について推定した時間造影強度曲線の一例を示す図である。
【００３９】
本実施形態の肝血流動態解析システム１は、図３に示すように、実際に撮影したＣＴ画像が、造影前の画像ＩＭ１、造影早期相の画像ＩＭ２、および造影後期相の画像ＩＭ３の３相画像だけであっても、任意の時刻における仮想ダイナミック像としての仮想ダイナミックＣＴ画像ＩＭｎを作り出すシステムである。
【００４０】
本実施形態においては、任意の時刻における仮想ダイナミックＣＴ画像ＩＭｎを作り出すために、組織の血流動態（造影血流の授受関係）を、血管と組織という２つのコンパートメント間の物質の移行としてモデル化し、組織の時間造影強度曲線を推定する方法を用いる。
【００４１】
モデル化する方法としては、薬物動態解析モデル解析法の一例としてのコンパートメントモデル解析法を用いる。コンパートメントモデル解析法は、血管と組織間の複雑な物質の濃度変化を、よりシンプルなモデル式に変換し、シミュレーションを行って、時間造影強度曲線を推定する方法である。この方法を用いることにより、飛び飛びの測定値から、組織の滑らかな時間造影強度曲線を推定することができる。
【００４２】
コンパートメントモデルの入力には、流入血管の時間造影強度曲線が必要となるが、従来の方法では、３相画像から得られるデータだけでは、時間分解能が不十分である。そこで、本実施形態においては、まず、図４に示すような、３相画像における、肝臓の流入血管である腹部大動脈と門脈の造影信号強度に基づいて、腹部大動脈と門脈の時間造影強度曲線を推定する。図４においては、０秒が造影前、４０秒が造影早期相、および１２０秒が造影後期相のそれぞれの腹部大動脈と門脈の造影信号強度を示している。
【００４３】
本実施形態では、これらの３相の造影信号強度に基づいて補間を行うことにより、図５に示すように、腹部大動脈と門脈の時間造影強度曲線を推定する。補間の詳細については後述する。
【００４４】
以上のように、本実施形態における肝血流動態解析システム１の処理は、流入血管の時間造影強度曲線の推定処理と、コンパートメントモデル解析法による組織の時間造影強度曲線の推定処理と、仮想ダイナミックＣＴ画像の表示処理とに分かれている。以下、これらの３つの処理について、それぞれ詳しく説明する。
【００４５】
［Ａ：流入血管の時間造影強度曲線の推定処理］
本実施形態においては、日常的な検査で広く用いられる、患者の造影前、早期相、後期相の３相画像データから高精度な流入血管の時間造影強度曲線を推定する。そのため、造影剤の体循環に関する臨床的事前知識を利用する。また、多相撮影された複数患者の血流動態データベースを用意する。以下、本実施形態における流入血管の時間造影強度曲線の推定方法について詳しく説明する。
【００４６】
（Ａ−１：学習フェイズ）
本実施形態における流入血管の時間造影強度曲線の推定方法は、多相撮影された複数患者の血流動態データベースを用いて、少ない撮影回数で得られた画像データから、主血管(大動脈,門脈等)の尤もらしい時間造影強度曲線を推定する方法である。この方法は、大きく分けて、学習フェイズと、テストフェイズとを備えている。まず、学習フェイズについて説明する。
【００４７】
造影剤は、腕の静脈から注入された後、様々な血管(血流路)を通り全身を回って拡散し、循環する。血流路のなかには、短いものや長いものが存在する。例えば、図６に示すように、造影剤は、ＲＡ（右心房）からＬＡ（左心房）を通り、ＡＯからＲＲＶ（右腎静脈）への腎臓を通る排出系の血流路は短い時間で循環する。また、ＡＯ（大動脈）からＳＭＶ（上腸間膜静脈）への消化器系を通る消化系の血流路は、中間の時間で循環する。さらに、ＡＯ（大動脈）からＩＶＣ（下大静脈）への末端を通る全身系の血流路は、足の末端も含むため、最も時間を要して循環する。
【００４８】
（被験者６０人分の１０時相における平均ＣＴ値の取得）
そこで、本実施形態における学習フェイズは、ＬＡ（左心房），ＡＡＯ（腹部大動脈），ＲＲＶ（右腎静脈），ＳＰＶ（脾静脈），ＳＭＶ（上腸間膜静脈），ＩＶＣ（下大静脈），ＰＶ（門脈）の７箇所の部位に注目領域を設定し、それぞれの注目領域での１０時相の平均ＣＴ値を測定する。そして、このような注目領域での１０時相の平均ＣＴ値を６０人の被験者について用意して、入力データとして用いる。なお、６０人はあくまで目安であり被験者の数は問わない。なお設定する注目領域は、上で述べた血管や臓器に限定されることなく、大動脈およびその分枝（総頚動脈、内頚動脈、外頚動脈、腕頭動脈、鎖骨下動脈、椎骨動脈、気管支動脈、肋間動脈、副腎動脈、卵巣動脈、精巣動脈、腰動脈）、腹腔動脈およびその分枝（総肝動脈、固有肝動脈、胃十二指腸動脈、上前膵十二指腸動脈、後膵十二指腸動脈、左胃動脈、右胃動脈、右胃大網動脈、左胃大網動脈、背側膵動脈、脾動脈、大膵動脈）、上腸間膜動脈およびその分枝（膵十二指腸動脈、空腸動脈、回結腸動脈、右結腸動脈、中結腸動脈）、下腸間膜動脈およびその分枝、総腸骨動脈、外腸骨動脈およびその分枝、内腸骨動脈およびその分枝、門脈、肺動脈のうち一つが含まれていればよい。
【００４９】
ここで、１０時相の平均ＣＴ値とは、造影剤注入後撮影時刻０秒，２２秒，２８秒，３４秒，４０秒，４６秒，５２秒，５８秒，９０秒，２１０秒における平均ＣＴ値をいう。なお、造影剤注入後撮影時刻０秒とは、造影前画像の撮影時刻としている。この時刻や数はこの限りではない。また、平均ＣＴ値とは、注目領域に相当する所定のボクセル領域の平均ＣＴ値のことである。図７は、ＬＡ（左心房）に設定した注目領域における１０時相の平均ＣＴ値の一例を示している。図７では、被験者６０人分の１０時相の平均ＣＴ値のうち、一部を示している。
【００５０】
（造影信号強度の計算）
以上のような各部位における被験者６０人分の１０時相の平均ＣＴ値を取得した後、被験者毎、部位毎で、造影前画像（０秒）のＡＡＯ (腹部大動脈)の画素値と、０秒を除く９時相の注目領域における画素値との差分値を算出することにより、造影信号強度を計算し、図８に示すように造影強度曲線を求める。
【００５１】
（造影血流循環の速さ(体内時計)の正規化）
次に、造影血流循環の速さ(体内時計)の正規化を行う。この正規化は、被験者毎、部位毎で、造影前を除く造影撮影の全てのタイミング、つまり、０秒を除く９時相における造影剤到達スコアを計算することにより行う。
【００５２】
前記正規化の計算は、以下のようなアルゴリズムに従う。
（ａ）各部位の中で最大の造影信号強度を取得する。
（ｂ）各部位において、最大の造影信号強度に対する割合を０から１の範囲で表す。つまり、ある部位の造影信号強度が最大の造影信号強度である場合には、当該部位に最大値１を割り当てる。しかし、ある部位の造影信号強度が最大の造影信号強度ではない場合には、最大の造影信号強度を１とした場合の当該部位の造影信号強度の割合を計算する。
（ｃ）各部位において造影剤が到達したかを判断する。最大の造影信号強度の部位には、造影剤が多く存在していると考えられるので、最大値１が割り当てられた部位には、造影剤が到達したと判断する。
【００５３】
上記のアルゴリズムを、９時相について、時刻の早い順（０秒に近い順）に行う。ある時相において最大値１が割り当てられ、造影剤が到達したと判断された部位は、残りの時相においても最大値１を割り当てる。
【００５４】
図９は、ＡＡＯ（腹部大動脈）における造影血流循環の速さ(体内時計)を正規化したデータの一例を示している。また、図１０は、ＲＲＶ（右腎静脈）における造影血流循環の速さ(体内時計)を正規化したデータの一例を示している。なお、図９および図１０においては、被験者６０人分のデータのうちの一部のデータを示している。
【００５５】
（造影信号強度の正規化）
次に、造影信号強度の正規化を行う。本実施形態では、被験者毎、部位毎で、造影前を除く造影撮影の全てのタイミング、つまり、０秒を除く９時相における造影信号強度を正規化する。この正規化は、次のようなアルゴリズムに従う。
（ａ）各部位の造影信号強度の最大値を取得する。
（ｂ）各部位において、造影信号強度の最大値に対する割合を［０，１］で表す。つまり、ある部位の造影信号強度が最大値である場合には、当該部位に最大値１を割り当てる。しかし、ある部位の造影信号強度が最大値ではない場合には、当該部位に最小値０を割り当てる。
【００５６】
（対象撮影時相のDynamic phase index(ＤＰＩ)の計算）
次に、対象部位(血管)のDynamic phase index(ＤＰＩ)を計算する。ＤＰＩは、対象部位の基準正規化信号強度曲線を推定するための指標であり、撮像時相における造影剤の到達距離の指標である。ＤＰＩは、以下の式を用いて最適化計算により求められる。
【００５７】
ＤＰＩ＝ｗ_ｘ（１）＊ｐＡＡＯ＋ｗ_ｘ（２）＊ｐＲＲＶ＋ｗ_ｘ（３）＊ｐＳＰＶ＋ｗ_ｘ（４）＊ｐＳＭＶ＋ｗ_ｘ（５）＊ｐＩＶＣ
【００５８】
上記の式において、ｗ_ｘ（ｉ）は最適化の変数ベクトルであり、重み係数を表す。ｉは１から５までの数値をとる。また、ｐが付された指標は、各血管の造影剤到達スコアであり、それぞれ以下の内容となっている。
【００５９】
ｐＡＡＯ・・・腹部大動脈の造影剤到達スコア
ｐＲＲＶ・・・右腎静脈の造影剤到達スコア
ｐＳＰＶ・・・脾静脈の造影剤到達スコア
ｐＳＭＶ・・・上腸間膜静脈の造影剤到達スコア
ｐＩＶＣ・・・下大静脈の造影剤到達スコア
【００６０】
ＤＰＩの計算においては、上記の各血管の造影剤到達スコアの他、対象部位の正規化信号強度ｒＸ、および対象部位の平衡相（２１０秒）での正規化信号強度の平均値ｒＸ_ｅｑを入力データとして用いる。なお、平均値ｒＸ_ｅｑは、対象部位の正規化信号強度ｒＸを学習データ数（６０人）で割って平均化したものである。
【００６１】
ＤＰＩの計算における最適化変数の範囲は、それぞれ０〜１の範囲であり、ピークの前と後で分けたときの早期相および後期相それぞれの多項式近似の次数は、それぞれ４次となっている。
【００６２】
本実施形態では、対象部位の正規化造影強度曲線からピークタイミングを計算し、早期相、後期相タイミングに分割、多項式近似を行った際の推定誤差を評価値とし、評価値を最小化する変数ベクトルｗ_ｘを求める。
【００６３】
ＤＰＩの計算における最適化の処理は、以下の手順で行う。
（ａ）まず、ＤＰＩを計算する。図１１は、計算したＤＰＩをプロットした図を示している。計算に使用するタイミングは、連続撮影されているタイミング（２２秒，２８秒，３４秒，４０秒，４６秒，５２秒，５８秒）を使用する。
（ｂ）次に、図１２に示すような正規化信号強度曲線を用いて、図１３に示すように、ＤＰＩを早期相帯と後期相帯とに分離する。
（ｃ）次に、早期相について正規化信号強度をＤＰＩで多項式フィッティングする。図１４（Ａ）は、早期相について正規化信号強度をＤＰＩで多項式フィッティングした例を示している。
（ｄ）そして、後期相について正規化信号強度をＤＰＩで多項式フィッティングする。図１４（Ｂ）は、後期相について正規化信号強度をＤＰＩで多項式フィッティングした例を示している。
【００６４】
フィッティング誤差や反復回数に関する条件を満足するまで、（ａ）から（ｄ）を繰り返す。その結果、規格化された、対象部位（血管）の曲線を得ることができる。
【００６５】
以上のように、ＤＰＩを計算することにより、標準的な血流循環モデルを作成することができる。
【００６６】
（対象患者のDynamic speed index(ＤＳＩ)を計算）
次に、対象患者のDynamic speed index(ＤＳＩ)を計算する。ＤＳＩは、循環時間速さ(体内時計)を求めるための指標であり、対象患者固有の循環速度を示す指標である。上述したＤＰＩの時間特性は、患者によって様々なので、ＤＳＩを計算することによって循環時間速さ(体内時計)を求め、循環時間速さ(体内時計)により、ＤＰＩを正規化する。その結果、標準的な血流循環モデルを作成することができる。ＤＳＩは、以下の式を用いて最適化計算により求められる。
【００６７】
ＤＳＩ＝ｗ_ｘ（１）＊ｐＡＡＯ(Ｉ，Ｊ)＋ｗ_ｘ（２）＊ｐＲＲＶ(Ｉ，Ｊ)＋ｗ_ｘ（３）＊ｐＳＰＶ(Ｉ，Ｊ)＋ｗ_ｘ（４）＊ｐＳＭＶ(Ｉ，Ｊ)＋ｗ_ｘ（５）＊ｐＩＶＣ(Ｉ，Ｊ)
【００６８】
上記の式において、ｗ_ｘ（ｉ）は最適化の変数ベクトルであり、重み係数を表す。ｉは１から５までの数値をとる。また、ｐが付された指標は、ＤＰＩの計算に用いた指標と同様であり、各血管の造影剤到達スコアを表す。それぞれ以下の内容となっている。
【００６９】
ｐＡＡＯ・・・腹部大動脈の造影剤到達スコア
ｐＲＲＶ・・・右腎静脈の造影剤到達スコア
ｐＳＰＶ・・・脾静脈の造影剤到達スコア
ｐＳＭＶ・・・上腸間膜静脈の造影剤到達スコア
ｐＩＶＣ・・・下大静脈の造影剤到達スコア
【００７０】
上記の式において、各造影剤到達スコアに用いられる変数ＩとＪは、次のような意味を持つ。まず、変数Ｉは、患者番号である。次に、変数Ｊは、撮影時刻を示す指標を表す。撮影時刻は、Ｊ＝１のとき２２秒、Ｊ＝２のとき３４秒、Ｊ＝３のとき４０秒、Ｊ＝４のとき４０秒、Ｊ＝５のとき４６秒、Ｊ＝６のとき５２秒、および、Ｊ＝６のとき５８秒を示す。本実施形態では、６０人の患者について、これらの６相の撮影時刻における各血管の造影剤到達スコアを用い、これらのスコアの線形和を求める。
【００７１】
対象部位のＤＳＩは、上記の式を用いて、最適化計算を行うことにより求める。最適化変数の範囲は、それぞれ０〜１の範囲である。
【００７２】
ＤＳＩの最適化(最小化)計算を行う際には、上記ＤＳＩの多項式近似を行った際の推定誤差を評価値とし、この評価値を最小化する変数ベクトルｗ_ｘを求める。最適化(最小化)のアルゴリズムは以下のようになる。
【００７３】
（ａ）重み係数ｗ_ｘにおけるＤＳＩを計算する。例えば、図１５（Ａ），（Ｂ）のように６０人の患者のうちの５人について、重み係数ｗ_ｘにおけるＤＳＩを計算する。
（ｂ）次に、ＤＳＩ時刻（Ｔｍ＝ＤＳＩ＊Ｔ）を計算する。図１６は、時刻Ｔと、ＤＳＩ時刻（Ｔｍ＝ＤＳＩ＊Ｔ）の関係を示す図である。
（ｃ）次に、ＤＳＩ時刻（Ｔｍ＝ＤＳＩ＊Ｔ）を用いて、ＤＰＩを多項式フィッティングする。図１７にＤＳＩ時刻（Ｔｍ）を用いて、ＤＰＩを多項式フィッティングした例を示す。
（ｄ）そして、フィッティング誤差や反復回数に関する条件を満足するまで、（ａ）から（ｃ）を繰り返す。
【００７４】
以上のようにＤＰＩを求め、その後にＤＳＩを求め、ＤＳＩ時刻（Ｔｍ）を用いてＤＰＩを多項式フィッティングすることにより、対象部位(血管)の時刻と正規化信号強度の関係を求め、患者によって異なるＤＰＩの時間特性を統一化することができる。つまり、ＤＳＩというパラメータを導入することにより、患者間で循環速度のばらつきがあり、実時間Ｔでは比較できなかったＤＰＩをＤＳＩ時刻（Ｔｍ）によって比較することができる。
【００７５】
学習フェィズにおいては、予め、入力データとしたＬＡ（左心房），ＡＡＯ（腹部大動脈），ＲＲＶ（右腎静脈），ＳＰＶ（脾静脈），ＳＭＶ（上腸間膜静脈），ＩＶＣ（下大静脈），ＰＶ（門脈）の７箇所の部位に注目領域に対し１０相にわたって測定された６０人分の平均ＣＴ値を複数患者血管信号値データベース４１として、記憶部４０に記憶させる。
【００７６】
また、計算で求めたＤＰＩ、ＤＳＩ、重み係数、および各造影剤到達スコア等のパラメータを、血流動態パラメータ記憶部に記憶しておいてもよい。予め計算し記憶しておいた血流動態パラメータを使用することも可能である。
【００７７】
（Ａ−２：テストフェイズ）
次に、テストフェイズについて説明する。テストフェイズでは、対象患者における３相画像の観測値から、ＤＳＩとＤＰＩを求め、流入血管の時間造影強度曲線を計算する。以下、テストフェイズにおける処理の手順について説明する。
【００７８】
（Ａ−２−１）各部位における注目領域の画素値と造影前画像の腹部大動脈（ＡＡＯ）の画素値との差分値の計算
まず、対象患者造影強度測定部１１は、対象患者において、各部位における注目領域の画素値と、造影前画像の腹部大動脈（ＡＡＯ）の画素値との差分値ｄＸ（ｊ）を、以下の式により計算する。
【００７９】
ｄＸ（ｊ）＝ｍａｘ（０，（Ｘ（ｊ）−ＡＡＯ（１）））
【００８０】
ここで、ｊは、撮影時刻番号であり、例えば、ｊ＝１のとき０秒、ｊ＝２のとき４０秒、ｊ＝３のとき１３０秒を示す。ＡＡＯ（１）は、造影前画像の腹部大動脈（ＡＡＯ）の画素値である。Ｘ（ｊ）は、注目領域における各撮影時刻の画素値である。また、ｍａｘ（Ａ，Ｂ，．．．）という演算は、ＡやＢ、その他全ての引数を比較して大きい方を結果として取得する演算である。
【００８１】
図１８に、ＬＡ（左心房），ＡＡＯ（腹部大動脈），ＲＲＶ（右腎静脈），ＳＰＶ（脾静脈），ＳＭＶ（上腸間膜静脈），ＩＶＣ（下大静脈），ＰＶ（門脈）に０秒、４０秒、１３０秒における画素値の例を示す。図１９は、図１８に示す画素値をプロットした図である。
【００８２】
（Ａ−２−２）対象患者の造影血流循環の速さ(体内時計)の正規化
次に、時間造影強度曲線推定部１５は、対象患者の造影血流循環の速さ(体内時計)の正規化を行う。この正規化を行うために、まず、対象患者の各部位において、早期相（４０秒）のタイミングでの造影剤到達スコアを計算する。計算は、以下のようなアルゴリズムにより行う。
【００８３】
（ａ）まず、各部位の中で最大信号強度Ｍａｘを取得する。
Ｍａｘ＝ｍａｘ（ｔ_ｄＬＡ，ｔ_ｄＡＡＯ，ｔ_ｄＲＲＶ，ｔ_ｄＳＰＶ，ｔ_ｄＳＭＶ，ｔ_ｄＰＶ，ｔ_ｄＩＶＣ）
（ｂ）次に、各部位Ｘ（ＬＡ，ＡＡＯ，ＲＲＶ，ＳＰＶ，ＳＭＶ，ＩＶＣ，ＰＶ）において早期相（４０秒）のタイミングでの造影剤到達スコアｔ_ｐＸ（ｔ_ｐＬＡ，ｔ_ｐＡＡＯ，ｔ_ｐＲＲＶ，ｔ_ｐＳＰＶ，ｔ_ｐＳＭＶ，ｔ_ｐＩＶＣ，ｔ_ｐＰＶ）を計算する。造影剤到達スコアｔ_ｐＸは、最大信号強度Ｍａｘに対する早期相（４０秒）のタイミングでの各部位の差分値ｔ_ｄＸの割合として表され、以下の式により計算する。
ｔ_ｐＸ＝ｔ_ｄＸ／Ｍａｘ
【００８４】
ここで、ｔ_ｄＸは、早期相（４０秒）のタイミングにおける、ある部位Ｘの画素値と造影前のＡＡＯの画素値との差分値を示す。
【００８５】
（ｃ）流路に従った各部位Ｘの造影剤到達スコアｔ_ｐＸの修正
次に、時間造影強度曲線推定部２２は、流路に従って、各部位Ｘの造影剤到達スコアｔ_ｐＸ（ｔ_ｐＡＡＯ，ｔ_ｐＲＲＶ，ｔ_ｐＳＰＶ，ｔ_ｐＳＭＶ）を、次式により修正する。
ｔ_ｐＡＡＯ＝ｍａｘ（ｔ_ｐＡＡＯ，ｔ_ｐＲＲＶ，ｔ_ｐＳＭＶ，ｔ_ｐＩＶＣ）
ｔ_ｐＲＲＶ＝ｍａｘ（ｔ_ｐＲＲＶ，ｔ_ｐＳＭＶ，ｔ_ｐＩＶＣ）
ｔ_ｐＳＰＶ＝ｍａｘ（ｔ_ｐＳＰＶ，ｔ_ｐＩＶＣ）
ｔ_ｐＳＭＶ＝ｍａｘ（ｔ_ｐＳＭＶ，ｔ_ｐＩＶＣ）
【００８６】
（３）対象患者のＤＳＩの算出
次に、時間造影強度曲線推定部２２は、対象患者のＤＳＩを算出する。対象患者のＤＳＩであるｔ_ＤＳＩの算出には、学習フェイズで算出したＤＳＩ重みパラメータを使用する。
ｔ_ＤＳＩ＝ｄｓｉ_ｗ（１）＊ｔ_ｐＡＡＯ＋ｄｓｉ_ｗ（２）＊ｔ_ｐＲＲＶ＋ｄｓｉ_ｗ（３）＊ｔ_ｐＳＰＶ＋ｄｓｉ_ｗ（４）＊ｔ_ｐＳＭＶ＋ｄｓｉ_ｗ（５）＊ｔ_ｐＩＶＣ
【００８７】
ここにおいて、ｄｓｉ_ｗは、各部位におけるＤＳＩの重みパラメータである。
【００８８】
（４）対象患者のＤＰＩの算出
次に、時間造影強度曲線推定部２２は、対象患者のＤＰＩを算出する。対象患者のＤＰＩの算出には、以下の式を用いて行い、学習フェイズで算出したＤＰＩ重みパラメータを使用する。
【００８９】
ｔ_ＤＰＩ＝ｄｐｉ_ｗ（１）＊ｔ_ｐＡＡＯ＋ｄｐｉ_ｗ（２）＊ｔ_ｐＲＲＶ＋ｄｐｉ_ｗ（３）＊ｔ_ｐＳＰＶ＋ｄｐｉ_ｗ（４）＊ｔ_ｐＳＭＶ＋ｄｐｉ_ｗ（５）＊ｔ_ｐＩＶＣ
【００９０】
ここにおいて、ｄｐｉ_ｗは、各部位のＤＰＩの重みパラメータである。
早期相のタイミングでの造影剤到達スコアｐＸを使用する。
【００９１】
（5）対象患者における対象血管の時間造影強度曲線の推定
次に、時間造影強度曲線推定部１５は、対象患者の対象血管の時間造影強度曲線を推定する。サンプリング時刻は、１〜２１０秒とする。対象血管の時間造影強度曲線の推定処理は、以下のような手順により行う。
（ａ）まず、時間造影強度曲線推定部１５は、ＤＳＩ時刻ｍＴを以下の式により計算する。
ｍＴ＝ｔ_ＤＳＩ＊Ｔ
ここにおいて、ｔ_ＤＳＩは対象患者のＤＳＩであり、Ｔはサンプリング時刻である。図２０は、サンプリング時刻ＴとＤＳＩ時刻ｍＴとの関係を示す図である。
【００９２】
（ｂ）次に、時間造影強度曲線推定部２２は、ＤＳＩ時刻ｍＴから、平衡相時間帯であるかどうかを表す変数ｅｑＰｈａｓｅを以下の式により計算する。
ｅｑＰｈａｓｅ＝ｍＴ＞ｅｑ_ｌｉｍｉｔ ？  １  ：  ０
なお、条件式？Ａ：Ｂという演算は、条件式が真であればＡを、偽であればＢを返す演算である。
【００９３】
ここで、ｅｑ_ｌｉｍｉｔは、各血管での平衡相開始時刻(定数)を示す。
図２１は、サンプリング時刻ＴとＤＳＩ時刻ｍＴとの関係を示す図に平衡相開始時刻ｅｑ_ｌｉｍｉｔを示した図であり、図２２は、サンプリング時刻Ｔと平衡相時間帯ｅｑＰｈａｓｅとの関係を示す図である。
【００９４】
（ｃ）時間造影強度曲線推定部１５は、ＤＳＩ時刻と学習フェイズで求めたＤＰＩ重みパラメータの多項式フィッティングにより、ＤＳＩ時刻に対応するＤＰＩを計算する。この際、ＤＰＩの値がマイナスの値になる場合には、ゼロに変換する。例えば、図２３のようにＤＰＩが得られた場合には、図２４のようにＤＰＩの値を変換する。
【００９５】
（ｄ）時間造影強度曲線推定部１５は、早期相での観測画素値を用いてＤＰＩを補正する。
図２４は、時刻と補正前のＤＰＩとの関係を示す図であり、図２５は、時刻と補正後のＤＰＩとの関係を示す図である。
【００９６】
（ｅ）時間造影強度曲線推定部１５は、学習フェイズで得られたＤＰＩピーク値の学習データ数分の平均値ｐｅａｋＤＰＩを用いて、早期相と後期相を分離する。図２６は、図２５にＤＰＩピーク値の学習データ数分の平均値ｐｅａｋＤＰＩの直線を示した図である。図２７は、図２６におけるＤＰＩの経時的変化を示す曲線と、平均値ｐｅａｋＤＰＩの直線との交点に基づいて、早期相、後期相、および平衡相を分離した例を示す図である。
【００９７】
（ｆ）時間造影強度曲線推定部１４は、上述した処理により分離した早期相と後期相とのそれぞれで、ＤＰＩと学習フェイズで得られた多項式係数とを用いて正規化信号強度曲線を計算する。図２８は、図２７に示す早期相で計算した正規化信号強度曲線の例を示している。また、図２９は、図２７に示す後期相で計算した正規化信号強度曲線の例を示している。そして、図３０は、図２８に示す正規化信号強度曲線と、図２９に示す正規化信号強度曲線とを連結して作成した早期相、後期相、および平衡相の正規化信号強度曲線を示している。
【００９８】
図２８におけるｅＤＰＩとは、早期相における推定ＤＰＩ値を示し、ｅＲｘとは、波形の振幅が０〜１に正規化された基準波形を示す。また、図２９におけるｄＤＰＩとは、後期相における推定ＤＰＩ値を示す。さらに、図３０におけるｍｏｄｉｆｉｅｄ＿ｅｄＰＩ ｖｓ ｅＲｘとは、図２８に示す早期相における推定ＤＰＩ値と基準波形との関係、および図２９に示す後期相における推定ＤＰＩ値と基準波形との関係とを連結したことを意味する。
【００９９】
（ｇ）時間造影強度曲線推定部１５は、平衡相における正規化信号強度曲線を推定する。図３１は、図２７と同様に、平均値ｐｅａｋＤＰＩに基づいて、早期相、後期相、および平衡相を分離した図、図３２は、上記（ｂ）の処理で計算した、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅを示す図である。
【０１００】
図３１と図３２とを比較すると明らかなように、平均値ｐｅａｋＤＰＩに基づいて早期相、後期相、および平衡相を分離した場合の平衡相の開始時刻ｔ１と、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２とは異なっている。本実施形態では、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２から、正規化信号強度曲線の計算を行い、平均値ｐｅａｋＤＰＩに基づいて分離した平衡相の開始時刻ｔ１から、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２までの正規化信号強度はゼロとしている。図３３は、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２から計算を行った正規化信号強度曲線の例を示している。
【０１０１】
図３４は、早期相の正規化信号強度曲線と、後期相の正規化信号強度曲線とを連結した正規化信号強度曲線を示す図である。本実施形態では、この早期相と後期相の正規化信号強度曲線に対して、ＤＳＩを使用して、対応する時刻計算を行い、サンプリング時刻Ｔに対する早期相と後期相の正規化信号強度曲線を求めている。
【０１０２】
図３５は、以上のようにして求めたサンプリング時刻Ｔに対する早期相と後期相の正規化信号強度曲線、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２から計算を行った正規化信号強度曲線、および平均値ｐｅａｋＤＰＩに基づいて分離した平衡相の開始時刻ｔ１から、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２までの正規化信号強度を連結した例を示している。
【０１０３】
（ｈ）時間造影強度曲線推定部１５は、対象血管の時間造影強度曲線の最大値を推定する。
対象血管の時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅは、もし、対象患者のｔ_ＤＰＩが対象血管の上限下限ＤＰＩの範囲外ならば、次式により、ｔ_ＤＰＩの値を使用せずに推定する。
【０１０４】
Ａｍｐｌｉｔｕｄｅ＝ｔ_ｄＸ（３）／ｍｅａｎ（ｒＸ（Ｉ，８））
ここで、ｔ_ｄＸ（ｉ）は、対象血管の造影信号強度であり、ｔ_ｄＸ（１）は造影前の対象血管の造影信号強度、ｔ_ｄＸ（２）は早期相の対象血管の造影信号強度、ｔ_ｄＸ（３）は後期相の対象血管の造影信号強度をそれぞれ示す。また、ｍｅａｎ（Ａ）は、ベクトルまたは行列Ａの平均を求める演算である。一例として、図３６に、ｔ_ｄＸ（３）の対象血管の造影信号強度を示す。
【０１０５】
また、ｒＸ（Ｉ，Ｊ）は、学習データにおける対象血管の正規化信号強度である。変数Ｉは、患者番号であり、Ｊは、撮影時刻を示す指標を表す。Ｊ＝８は、９０秒を示す。つまり、ｒＸ（Ｉ，８）は、患者番号Ｉの患者の９０秒における学習データの対象血管の正規化信号強度である。一例として、図３７に、ｒＸ（Ｉ，８）の学習データにおける対象血管の正規化信号強度を示す。
【０１０６】
一方、対象血管の時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅは、対象患者のｔ_ＤＰＩが対象血管の上限下限ＤＰＩの範囲内ならば、次式により、ｔ_ＤＰＩの値を使用して推定する。
【０１０７】
ｔ_ＤＰＩが、ｐｅａｋＤＰＩ以下ならば、対象血管の時間造影強度曲線の最大値は、次式から推定される。
Ａｍｐｌｉｔｕｄｅ＝ｔ_ｄＸ（２）／ｐｏｌｙｖａｌ（ｅＤＰＩｐ，ｔ_ＤＰＩ））
ここで、ｅＤＰＩｐは、学習フェイズで求めた早期相の多項式係数である。なお、ｐｏｌｙｖａｌ（Ａ，Ｂ）は、要素ベクトルＡと多項式係数Ｂによる多項式の値を得る演算である。
【０１０８】
また、ｔ_ＤＰＩが、下限のＤＰＩ以上ならば、対象血管の時間造影強度曲線の最大値は、次式から推定される。
Ａｍｐｌｉｔｕｄｅ＝ｔ_ｄＸ（２）／ｐｏｌｙｖａｌ（ｄＤＰＩｐ，ｔ_ＤＰＩ））
ここで、ｄＤＰＩｐは、学習フェイズで求めた早期相の多項式係数である。
【０１０９】
（ｉ）時間造影強度曲線推定部１５は、最大値の推定結果を適用する。
図３５で示した全時間的な正規化信号強度曲線に、上述した（ｈ）の処理で推定した対象血管の時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅを適用することにより、図３８に示す時間造影強度曲線が得られる。
【０１１０】
（ｊ）時間造影強度曲線推定部１５は、実際の後期相の観測値を用いて、平衡相帯を補正する。
図３５に示す時間造影強度曲線のうち、平衡相帯を、実際の後期相の観測値を用いて補正すると、図３９に示す時間造影強度曲線が得られる。補正は以下の式により行う。
ｅｄＲｏｉ（ｔ（９）：ｔ（１０），ｒｏｉ）＝ｅｄＲｏｉ（ｔ（９）：ｔ（１０），ｒｏｉ）＊（ｔ_ｄＲｏｉ（３，ｒｏｉ）／ｅｄＲｏｉ（ｔ_Ｔ（３），ｒｏｉ）
図３９は補正後の時間造影強度曲線を示し、図３５に示す時間造影強度曲線に対して、後期相の観測値ｔ_ｄＲｏｉ（３，ｒｏｉ）を掛けて補正している。なお、Ａ：Ｂは、任意の整数Ａ、Ｂにおいて、ＡからＢまで１刻みで増加もしくは減少するベクトルもしくは配列を作成する演算である。
【０１１１】
（ｋ）時間造影強度曲線推定部１５は、未計算の部分の計算(補間)を行う。
図３５の時間造影強度曲線において、平均値ｐｅａｋＤＰＩに基づいて分離した平衡相の開始時刻ｔ１から、ＤＳＩ時刻ｍＴに基づく平衡相時間帯ｅｑＰｈａｓｅの開始時刻ｔ２までの未計算の部分を線形補間すると、図４０に示す直線が得られる。また、図４１は、図３５の時間造影強度曲線に、上記未計算の部分の補間を行った後の時間造影強度曲線を示す図である。
【０１１２】
（ｌ）時間造影強度曲線推定部１５は、推定した対象血管の時間造影強度の最大値に達する時刻を計算する。
次に、図４１に示す線形補間後の時間造影強度曲線において、上述した（ｈ）の処理で推定した対象血管の時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅに達する時刻を計算する。図４２に示す例では、４１秒で時間造影強度曲線の最大値Ａｍｐｌｉｔｕｄｅに達している。
【０１１３】
（ｍ）時間造影強度曲線推定部１５は、最後に、推定した時間造影強度曲線のうち、後期相と平衡相の間を平衡相の直線を使って補正する。
図４３（Ａ）は、推定した時間造影強度曲線上に、図４３（Ｂ）で示した平衡相の直線を重ねて表示した図である。もし、平衡相の直線の値が、後期相の値よりも大きければ、後期相の値を平衡相の直線の値で書き換える。その後は、必要があれば、点の間を線形補間等の補間を行う。但し、殆どの場合、必要ないと考えられる。図４４は、このようにして補正を行った後の時間造影強度曲線を示す。
【０１１４】
以上のように、テストフェイズでは、対象患者における３相画像の観測値から求めたＤＰＩ，ＤＳＩパラメータを用いて、時刻から正規化信号強度を得て、流入血管の時間造影強度曲線を計算する。このように、テストフェイズにおいては、ＤＳＩとＤＰＩを３相画像の観測値から求めることができ、流入血管の時間造影強度曲線を簡単に計算可能である。
【０１１５】
［Ｂ：コンパートメントモデル解析法による組織の時間造影強度曲線の推定処理］
本実施形態では、以上ようにして計算した対象患者の流入血管の時間造影強度曲線と、薬物動態解析モデル係数データベース４２とを用いて、コンパートメントモデル解析法により、対象患者の組織の時間造影強度曲線を推定する。コンパートメントモデル解析法を用いた肝血流動態解析の利点としては、まず一つとして、初期値とコンパートメントモデル式の係数値が決まると、時間造影強度曲線を推定できることが挙げられる。また、コンパートメントモデル式の係数値によって、肝血流動態に関して定量評価ができることが挙げられる。
【０１１６】
本実施形態においては、実際の画素から造影強度の時間変化を見て計算によりコンパートメントモデル式の係数を求めるのではなく、造影強度の時間変化とコンパートメントモデル式の係数を１対１に対応させた薬物動態解析モデル係数データベース４２をシミュレーションにより作成し、この薬物動態解析モデル係数データベース４２から係数を検索する。
【０１１７】
薬物動態解析モデル係数データベース４２は、基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点と、線形および非線形的最適化に基づいて算出したコンパートメントモデル式の最適な係数を対応つけたデータベースである。
【０１１８】
（Ｂ−１）薬物動態解析モデル係数データベース４２の作成処理
まず、本実施形態における薬物動態解析モデル係数データベース４２の作成処理について説明する。
【０１１９】
肝臓には、２つの代表的な流入血管として門脈と肝動脈がある。門脈は栄養、肝動脈は酸素を運ぶ機能を有している。肝実質と肝癌とでは、支配血管が異なり、肝実質では門脈血が７０％で、肝動脈血は３０％程度である。また、腫瘍では、肝動脈などの動脈血が支配的になる
【０１２０】
コンパートメントモデル解析では、注目肝組織が、門脈血と肝動脈血のどちらが支配的かを定量的に評価することができる。また、血流流出の状態も定量的に評価できる。したがって、腫瘍の鑑別に有用である。
【０１２１】
コンパートメントモデル解析においては、肝動脈造影強度と門脈造影強度をどうブレンドをしたら、肝臓の造影具合を再現可能かを血流流出を加味し計算する。言い換えれば、薬物動態解析モデル解析部１２は、肝臓の造影具合から、動脈造影強度曲線と門脈造影強度曲線のブレンドの割合を逆算する。
【０１２２】
門脈血と肝動脈血のどちらが支配的かは、推定後のコンパートメントモデルの係数に現れる。以下に、肝臓部のコンパートメントモデル式を示す。なお肝動脈の造影強度は、腹部大動脈の造影強度で近似してもよい。
【０１２３】
2020168210.tif 000003
【０１２４】
上記の式における各係数の意味は以下のようになっている。
Ｃａ(ｔ)：肝動脈の造影強度（時間ｔ）
Ｃｐ(ｔ)：門脈の造影強度（時間ｔ）
Ｃｅｆｓ(ｔ)：細胞外液腔(extracellular fluid space)の造影強度（時間ｔ）
Ｋ１ａ：肝動脈から細胞外液腔への流入速度定数
Ｋ１ｐ：門脈から細胞外液腔への流入速度定数
ｋ２：細胞外液腔から静脈への流出速度定数
τa：肝動脈から細胞外液腔への到達時間
τp：門脈から細胞外液腔への到達時間
【０１２５】
また、以下に、脾臓のコンパートメントモデル式を示す。なお脾動脈の造影強度も、腹部大動脈の造影強度で近似してもよい。
2020168210.tif 000004
【０１２６】
上記の式における各係数の意味は以下のようになっている。
Ｃａ(ｔ)：脾動脈の造影強度（時間ｔ）
Ｃｅｆｓ(ｔ)：細胞外液腔の造影強度（時間ｔ）
Ｋ１：脾動脈から細胞外液腔への流入速度定数
ｋ２：細胞外液腔から静脈への流出速度定数
τa:脾動脈から細胞外液腔への到達時間
【０１２７】
薬物動態解析モデル係数データベース４２を作成する際に行うコンパートメントモデル解析では、臓器の流入血管との関係を利用し、各血管に関するパラメータをキャリブレ−ションする。また、このコンパートメントモデル解析においては、時定数パラメータは、各臓器組織に対し一定と仮定する。さらに、モデルの自由度を効率的に削減し、計算効率を上げている。
【０１２８】
実際の計算の流れは、以下のような３ステップから構成される。
１．脾臓ステップ …上記の脾臓のコンパートメントモデル式を用いて、動脈に関するパラメータを決定する。
２．肝臓ステップ …上記の肝臓部のコンパートメントモデル式を用いて、門脈に関するパラメータを決定する。
３．最終計算ステップ …上記の肝臓部のコンパートメントモデル式を用いて、対象組織の計算を全画素について行う。
以下、それぞれのステップについて詳しく説明する。
【０１２９】
（Ｂ−１−１）脾臓ステップ
脾臓ステップでは、シミュレーション部１２の第１決定部１２ａは、上記の脾臓のコンパートメントモデル式における大動脈の係数のキャリブレーションを行う。脾臓のコンパートメントモデル式を用いる理由は、脾臓に造影剤を送る血管が大動脈から分岐する脾動脈のみだからである。本実施形態では、第１決定部１２ａは、上述したテストフェイズで計算した脾臓の時間造影強度曲線を説明するように、脾臓のコンパートメントモデル式における大動脈の係数を決定する。具体的には、第１決定部１２ａは、脾動脈の造影強度Ｃａ(ｔ)の立ち上がり時間と、脾動脈から細胞外液腔への到達時間τaを決定する。
【０１３０】
図４５は、上述したテストフェイズで計算した脾臓と脾動脈の時間造影強度曲線を示す図である。脾動脈の造影強度Ｃａ(ｔ)の立ち上がり時間は、図４５に丸印で示す脾動脈の造影強度が０の点から、最初の観測値までの立ち上がり時間から算出する。
【０１３１】
また、脾動脈から細胞外液腔への到達時間τaは、図４５に丸印で示すように、脾臓の時間造影強度曲線において、時刻０から、脾臓の時間造影強度曲線が０を超える値に立ち上がるまでの時間から算出する。
【０１３２】
（Ｂ−１−２）肝臓ステップ
肝臓ステップでは、シミュレーション部１２の第２決定部１２ｂは、上記の肝臓のコンパートメントモデル式における門脈の係数のキャリブレーションを行う。肝臓のコンパートメントモデル式を用いる理由は、肝臓に造影剤を送る血管が大動脈と門脈だからである。本実施形態の肝臓ステップでは、第２決定部１２ｂは、大動脈に係る係数は脾臓ステップで求めたものを使用し、上述したテストフェイズで計算した肝臓の時間造影強度曲線を説明するように、肝臓のコンパートメントモデル式における門脈脈の係数を決定する。具体的には、門脈の造影強度Ｃｐ(ｔ)の立ち上がり時間と、門脈から細胞外液腔への到達時間τpとを決定する。
【０１３３】
図４６は、肝臓における上述したテストフェイズで計算した大動脈(肝動脈の代替)と門脈との時間造影強度曲線を示す図である。門脈の造影強度Ｃａ(ｔ)の立ち上がり時間は、図４６に丸印で示すように、門脈の造影強度が０の点から、最初の観測値までの立ち上がり時間から算出する。
【０１３４】
（Ｂ−１−３）最終計算ステップ
最終計算ステップでは、シミュレーション部１２の計算部１２ｃは、上記の肝臓のコンパートメントモデル式における残りの係数を得るために、全画素について、造影剤の各血管からの流入速度および注目組織からの流出速度を計算する。具体的には、計算部１２ｃは、肝動脈から外腔への流入速度定数Ｋ１ａと、門脈から外腔への流入速度定数Ｋ１ｐと、肝細胞外液腔から静脈への流出速度定数ｋ２とを決定する。
【０１３５】
以上のような３ステップにより、上記コンパートメントモデルを数値解析の１手法であるルンゲクッタ法で数値的に解析する。なお、数値解析法は上記のコンパートメントモデルのような微分方程式を解析的に解くことができればよく、ルンゲクッタ法に限らない。
【０１３６】
本実施形態では、初期値Ｃｅｆｓ(０)はゼロとし、解析時間帯を０秒から２１０秒で数値解析を行っているが、この限りではない。
【０１３７】
（Ｂ−１−４）コンパートメントモデル式の係数の最適化処理
薬物動態解析モデル係数データベース４２におけるコンパートメントモデル式の係数の算出は、遺伝的アルゴリズムと非線形最適化の組み合わせにより行う。非線形最適化は、勾配を利用するため、計算(収束)時間が早いが、初期値によって局所最適解に陥ることがある。遺伝的アルゴリズムは、勾配を利用しないため、大域最適解を得ることができるが、計算(収束)に時間がかかる。そこで、本実施形態では、遺伝的アルゴリズムで当たりの良い初期値を見つけ、そこから高速な非線形最適化を適用することで、計算時間を抑えつつ大域最適解または大域最適解に近い解を求めている。
【０１３８】
薬物動態解析モデル係数データベース４２においては、複数のシミュレーション点における造影強度の時間変化について、下記に説明するコンパートメントモデル式における係数の最適化処理を行い、複数のシミュレーション点とコンパートメントモデル式の係数とを１対１で対応させる。シミュレーション点は、上述した時間造影強度曲線の推定処理のテストフェイズで計算した肝臓の大動脈と門脈の全時間造影強度曲線で決定する。大動脈と門脈の全時間造影強度曲線は、ＲＯＩの平均値を使用しており、画素値ノイズの影響が少なく、また、上述のようにして求めた大動脈と門脈の全時間造影強度曲線は、６０人分の血流動態データを基礎にしており、対象患者の画像データだけで決定するよりも、信頼性が高いためである。また、組織の造影強度ではなく、血管の造影強度を用いるのは、組織の造影強度よりも血管の造影強度の方が大きく、シミュレーション範囲が適切になり、無駄になる計算が少なくなるためである。
【０１３９】
図４７に、大動脈と門脈との全時間造影強度曲線の一例を示す。本実施形態における薬物動態解析モデル係数データベース４２は、一例として、早期相と後期相の２相に着目して、複数のシミュレーション点における造影強度の時間変化を２次元空間にプロットする。そこで、早期相のシミュレーション範囲において所定の複数のシミュレーション点を取り、かつ、後期相のシミュレーション範囲において所定の複数のシミュレーション点を取る。
【０１４０】
早期相のシミュレーション範囲の最大値は、腹部大動脈の全時間造影曲線の最大値を採用し、早期相のシミュレーション範囲の最小値は、ゼロとする。また、後期相のシミュレーション範囲の最大値は、後期相の所定の撮影時刻、例えば、１３０［ｓｅｃ］における、腹部大動脈の全時間造影曲線の最大値または門脈の全時間造影曲線の最大値のいずれか大きい方を採用する。そして、後期相のシミュレーション範囲の最小値は、ゼロとする。ただし、これに限定されるものではなく、例えば早期相のシミュレーション範囲の最大値を早期相での右腎静脈の信号値にするなど工夫をしてもよい。
【０１４１】
本実施形態では、一例として、以上のような早期相のシミュレーション範囲を９等分し、かつ、後期相のシミュレーション範囲を９等分して、早期相のシミュレーション範囲を一方の次元、また、後期相のシミュレーション範囲を他方の次元として、２次元空間に１００個のシミュレーション点をプロットする。図４８にこのようにしてプロットした造影強度パターンのシミュレーション点の一例を示す。本実施形態では、このように配置したシミュレーション点のそれぞれに、コンパートメントモデル式の係数を１対１に対応させる。
【０１４２】
以下、それぞれのシミュレーション点におけるコンパートメントモデル式の係数の最適化処理について説明する。図４９に、薬物動態解析モデル係数データベース４２を作成する際のコンパートメントモデル式における係数の最適化の概略手順についてのフローチャートを示す。
【０１４３】
１．シミュレーション部１２は、サンプリング点(ｐ点)分、コンパートメントモデル式の係数にランダムに値を割り当てる（図４９：Ｓ１）。
図５０は、各係数にサンプリング点(ｐ点)分の値をランダムに割り当てた一例を示す表である。例えば、等高線で表された山の地図を例にとると、係数の最適解は山の頂点であり、ランダムに割り当てた各係数の値は、山の頂点を目指す登山者の位置を地図上にランダムに割り当てるイメージである。
【０１４４】
２．シミュレーション部１２は、それぞれの係数値でコンパートメントモデル式を数値解析で解き、時間造影強度曲線を計算する（図５０：Ｓ２）。
図５１（Ａ），（Ｂ），（Ｃ）は、図５０に示すサンプリング点０，１，ｐ−１における各係数値で、上述した肝臓のコンパートメントモデル式を数値解析で解き、時間造影強度曲線を計算した結果を示す。図５１（Ａ），（Ｂ），（Ｃ）においては、あるシミュレーション点に対応する早期相のＣＴ値と後期相のＣＴ値を丸印で示す。このシミュレーション点は、早期相の造影強度差分値が１０であり、後期相の造影強度差分値が６０のシミュレーション点である。計算結果は、推定曲線として描かれている。なお、図５１（Ａ），（Ｂ），（Ｃ）は、サンプリング点０，１，ｐ−１に対する例示であり、実際は、サンプリング点０〜ｐ−１までのｐ点分の計算結果が得られる。
【０１４５】
３．シミュレーション部１２は、評価値(観測値との誤差)を計算する（図５０：Ｓ３）。
観測値との誤差である評価値Ｅは、以下の式により求められる。
【０１４６】
2020168210.tif 000005
【０１４７】
上記の式における各係数の意味は以下のようになっている。
Ｃｅｆｓ(ｔ)：推定した細胞外液腔(extracellular fluid space)の造影強度（時間ｔ）
Ｃｍｅｓ(ｔ)：観測した実際の細胞外液腔(extracellular fluid space)の造影強度（時間ｔ）
【０１４８】
評価値は、図５２（Ａ），（Ｂ），（Ｃ）に示すように、各観測値と推定曲線上の対応値との距離として表される。
【０１４９】
４．シミュレーション部１２は、評価値が良いものを選択する（図５０：Ｓ４）。
評価値が良いものとは、誤差値が最小のものであり、図５２（Ａ），（Ｂ），（Ｃ）の３つ例示の中では、図５２（Ｂ）の場合に誤差値が最小であるため、サンプリング点１において割り当てた係数値を選択する。
【０１５０】
５．シミュレーション部１２は、遺伝的アルゴリズムの各手法に従って、最適解候補付近に割り当てるサンプル点位置を計算する（図５０：Ｓ５）。
これは、遺伝的アルゴリズムの交叉、突然変異、世代交代モデルに相当する操作である。これらの操作は、様々な方法があり、従来から行われている方法であるため、ここでは詳細な説明は割愛する。
【０１５１】
上述した等高線で表された山の地図を例で説明すると、山の頂点を目指す登山者の位置を、山の頂点を含む円内に割り当てるイメージである。
【０１５２】
６．シミュレーション部１２は、設定した条件 を満たすまで、上記ステップ２〜５を繰り返す（図５０：Ｓ６）。
上記ステップ２〜５を繰り返すことにより、上述した等高線で表された山の地図を例では、山の頂点を含む円の半径が徐々に狭められ、山の頂点を目指す登山者の位置が徐々に山の頂点に近づくイメージである。
【０１５３】
７．シミュレーション部１２は、コンパートメントモデル式の最適係数を得る（図５０：Ｓ７）。
以上のようにして、コンパートメントモデル式の係数の大域最適解または大域最適解に近い解を求めることができる。
【０１５４】
本実施形態では、シミュレーション部１２が、全てのシミュレーション点について上記解析を行い、それぞれのシミュレーション点に対応するコンパートメントモデルの最適係数を取得する。シミュレーション部１２は、シミュレーションを行う際に、大動脈の時定数τａと門脈の時定数τｐの再計算を行わず、速度定数の同定のみをシミュレーション点全てに対して実施する。その上で、図４８に示すように、早期相の造影強度差分値と後期相の造影強度差分値との２次元空間にプロットした各シミュレーション点と、コンパートメントモデル式の最適係数とを１対１に対応させ、薬物動態解析モデル係数データベース４２を作成する。
【０１５５】
（Ｂ−２）薬物動態解析モデル係数データベース４２の検索処理
次に、以上のようにして作成した薬物動態解析モデル係数データベース４２の検索処理について説明する。図５３は、本実施形態における薬物動態解析モデル係数データベース４２の検索処理を説明するための図である。
【０１５６】
図５３においては、対象患者の３相の観測値（差分値）を各シミュレーション点よりも大きな丸印で示している。本実施形態では、一例として、検索部１２は、観測値に最も近いシミュレーション点を薬物動態解析モデル係数データベース４２から検索する。図５３に示す例では、矢印で示すシミュレーション点が観測値に最も近く、検索部１２は、図５３に矢印で示すシミュレーション点を、観測値に最も近いシミュレーション点として検索する。なお、本実施形態では、観測値に最も近いシミュレーション点を検索する態様について説明したが、本発明はこのような態様に限定される訳ではない。例えば、機械学習等により、観測値に応じて最適なシミュレーション点を検索するようにすれば、検索の高速化および高精度化を狙うことも可能である。
【０１５７】
次に、取得部１３は、薬物動態解析モデル係数データベース４２から、検索したシミュレーション点に対応付けられたコンパートメントモデル式の最適係数を取得する。
【０１５８】
［Ｃ：仮想ダイナミックＣＴ画像の表示処理］
次に、仮想ダイナミック像生成部１６により、３相のＣＴ画像から、任意の時刻における仮想ダイナミック像を作り出す処理を、ビューワ２００の表示部２１１における表示例等と共に説明する。
【０１５９】
後述する出力画面において仮想ダイナミック像を表示させるには、仮想ダイナミック像の入力画面において、所定の条件を設定する必要がある。以下、この所定の条件について説明する。
【０１６０】
（１）コンパートメントモデル解析に使用するｎ相画像を選択する。
例えば、本実施形態の場合には、３相画像を選択する。
（２）ｎ相画像において必要な部位に注目領域を設定する。
注目領域の設定は、手動または自動で行うことができる。本実施形態では、手動で設定する例について説明する。
（３）造影剤注入後直後に撮影した画像の注入後経過時刻を入力する。
撮影時刻自体はＤＩＣＯＭタグで取得することが可能だが、造影剤注入時刻が不明であるため、造影剤注入後直後に撮影した画像の注入後経過時刻を入力する。
（４）解析対象臓器に応じたコンパートメントモデル式を用いて薬物動態解析モデル係数データベース４２を作成する。
コンパートメントモデル式は、肝臓のコンパートメントモデル式、または、肝臓のコンパートメントモデル式と脾臓のコンパートメントモデル式、あるいは、肝臓のコンパートメントモデル式と肝臓以外の腹部臓器コンパートメントモデル式等を適宜選択することができる。
（５）薬物動態解析モデル係数データベース４２を用いて、コンパートメントモデル式の最適係数を取得する。
【０１６１】
時間造影強度推定部１５は、以上のようにして取得されたコンパートメントモデル式の最適係数を用いて、時間造影強度曲線を推定する。この際、Ｋ１ｐをゼロとすると、仮想的に門脈の流入を無視した時間造影強度曲線(動脈下造影ＣＴ曲線)が得られ、Ｋ１ａをゼロとすると、仮想的に肝動脈の流入を無視した時間造影強度曲線(門脈下造影ＣＴ曲線)が得られる。仮想ダイナミック像生成部１６は、時間造影強度曲線に基づいて、全画素にわたり指定された任意時刻に対応するＣＴ値を、線形補間などで求めることで、ＣＴ値画像を得ることができる。なお、求める方法線形補間に限らない。
【０１６２】
図５４にビューワ２００における入力画面２０１の一例を示す。入力画面２０１には、図５４に示すように、解析に使う画像のスタックを複数選択可能に表示させる。図５４に示す例では、６個のスタックのうち、上段の真ん中のスタックが解析に使うスタックとして選択された場合を示している。スタックを選択し、メニューから血流動態解析ボタンをクリックすると、選択されたスタックには、図５４に示すように、部位ＲＯＩ指定ツール２０２が表示される。
【０１６３】
図５５に選択されたスタックと部位ＲＯＩ指定ツール２０２の拡大図を示す。部位ＲＯＩ指定ツール２０２には、大動脈、門脈、左心房、右腎静脈、脾静脈、上腸間膜静脈、下大静脈、脾臓、および肝臓（実質）を部位ＲＯＩとして指定する部位ＲＯＩ指定ボタン２０３が設けられている。また、部位ＲＯＩ指定ツール２０２には、表示する画面の撮影時刻を指定する撮影時刻ボタン２０４が設けられている。さらに、部位ＲＯＩ指定ツールに２０２は、使用するコンパートメントモデルを設定するラジオボックス２０５が設けられている。
【０１６４】
部位ＲＯＩ指定は、部位ＲＯＩ指定ボタン２０３のいずれか部位名のボタンをクリックした後、スタック上での部位の位置をクリックすることにより行われる。図５６には、大動脈のボタン２０３をクリックした後に、スタック上で当該部位の位置Ａをクリックした例を示している。ＲＯＩの大きさについては、予めプリセットが作成されており、そのプリセットの大きさのＲＯＩを設置する。修正が必要であればドラッグすることにより、ＲＯＩの大きさを修正することができる。
【０１６５】
次に、出力画面について説明する。図５７は、出力画面２２０の一例を示す図である。図５７示すように、出力画面は６分割されており、上段の３画面には、３相のＣＴ画像が表示されている。出力画面２２０上には、図５８（Ａ）に示す選択アイコン２２１が表示される。選択アイコン２２１は、ＶＤＣＴ，ＶＤＣＴＨＡ，ＶＤＣＴＡＰ，Ｋ１ａ，Ｋ１ｐ，Ｋ２，Ｖｄ，ＴＦの８個のアイコンに分かれている。ＶＤＣＴは仮想ダイナミックＣＴ画像、ＶＤＣＴＨＡは門脈造影を行った仮想ダイナミックＣＴ画像、ＶＤＣＴＡＰは肝動脈から造影した仮想ダイナミックＣＴ画像のそれぞれを表示する際に選択する。また、Ｋ１ａ，Ｋ１ｐ，Ｋ２，Ｖｄ，ＴＦは、それぞれのパラメータマップを表示する際に選択する。なお、パラメータのうち、Ｖｄは分布容積(Distribution volume)であり、ＴＦは総血流量(Total blood flow)、ＡＦは動脈血流比(Arterial blood flow)、ＰＦは門脈血流比(Portal venous blood flow)である。
【０１６６】
また、選択アイコン２２１は、図５８（Ｂ）に示すように、アイコンの数を増やしてもよい。図５８（Ｂ）においては、ＶＤＣＴは、大動脈ピーク時刻におけるＶＤＣＴ、大動脈の信号強度と門脈の信号強度の差分が最大になる時刻におけるＶＤＣＴ、および門脈ピーク時刻におけるＶＤＣＴに分かれている。また、ＶＤＣＴＨＡは、大動脈ピーク時刻におけるＶＤＣＴＨＡ、大動脈の信号強度と門脈の信号強度の差分が最大になる時刻におけるＶＤＣＴＨＡ、および門脈ピーク時刻におけるＶＤＣＴＨＡに分かれている。さらに、ＶＤＣＴＡＰは、大動脈ピーク時刻におけるＶＤＣＴＡＰ、大動脈の信号強度と門脈の信号強度の差分が最大になる時刻におけるＶＤＣＴＡＰ、および門脈ピーク時刻におけるＶＤＣＴＨＡに分かれている。また、選択アイコン２２１は、パラメータマップだけでなく、表示したい仮想ダイナミック像の時刻ごとに作成して表示することもできる。
【０１６７】
図５９は、出力画面２２０の表示例を示す図である。出力画面２２０においては、デフォルトで、代表的なタイミングでのＶＤＣＥ-ＣＴ像（仮想ダイナミックＣＴ画像）が表示される。例えば、大動脈信号強度ピークタイミングにおけるＶＤＣＥ-ＣＴ像が表示される。
【０１６８】
また、図５９に示す例では、ＶＤＣＴＡＰの選択アイコンが選択されており、ＶＤＣＥ-ＣＴ像の右隣に、ＶＤＣＥ-ＣＴＡＰ像（肝動脈から造影した仮想ダイナミックＣＴ画像）が表示される。ＶＤＣＥ-ＣＴ像、ＶＤＣＥ-ＣＴＡＰ像、およびＶＤＣＥ-ＣＴＨＡ像（門脈造影を行った仮想ダイナミックＣＴ画像）が表示される画面には、図５９に示すように、時刻スライダー２２２が表示される。時刻スライダー２２２を移動させることにより、指定時刻でのそれぞれの画像を表示させることができる。また、ショートカットキーによってＶＤＣＥ−ＣＴ像の指定時刻を変化させることができる。
【０１６９】
ＶＤＣＥ-ＣＴ像が表示され、または、ＶＤＣＥ-ＣＴＡＰ像、もしくはＶＤＣＥ-ＣＴＨＡ像が表示されているバックグラウンドでは、時刻スライダーを移動させるもしくはショートカットキーで変化させた指定時刻でのＶＤＣＥ-ＣＴ像、または、ＶＤＣＥ-ＣＴＡＰ像、もしくはＶＤＣＥ-ＣＴＨＡ像のボリュームデータを作成する。そして、時刻スライダー２２２が移動されて指定時刻が変わったタイミングで、ボリュームデータは破棄され、新たな指定時刻のボリュームデータを作り直すようになっている。
【０１７０】
また、図５９においては、Ｋ１ｐのパラメータの選択アイコンが選択されており、ＶＤＣＥ-ＣＴＡＰ像の右隣には、Ｋ１ｐのパラメータマップが表示されている。パラメータマップはカラー(ヒートマップ)表示できるようなっている。
【０１７１】
ＶＤＣＥ-ＣＴ像、または、ＶＤＣＥ-ＣＴＡＰ像、もしくはＶＤＣＥ-ＣＴＨＡ像を表示する場合には、仮想的に作成された画像であることが分かるようにオーバーレイを入れるようにしてもよい。
【０１７２】
例えば、出力画面２２０においては、それぞれの画面を例えばダブルクリックすることにより、図６０（Ａ）に示すように、それぞれの画面をオーバーレイ表示してもよい。また、時刻スライダー２２２で時刻指定を行う場合には、各部位の時間造影強度曲線のピーク造影強度の時刻にジャンプするようにしてもよい。図６０（Ａ）の例では、図６０（Ｂ）に示す大動脈のピーク造影強度の時刻と、図６０Ｃ）に示す門脈のピーク造影強度の時刻にジャンプするようにしている。このように、ピーク造影強度の時刻にジャンプすることにより、ピーク造影強度の時刻での対象部位の造影剤染まり具合の診断に有用である。また、２つの入力血管の差分値が最大となる時刻にジャンプすることもまた、診断に有用である。また、ショートカットキー等を用いて、指定された範囲の時刻をループするアニメーションとして描画してもよい。なお、ここでは、大動脈のピーク造影強度の時刻を例示したが、これに限ることなく、例えば、肝動脈のピーク造影強度の時刻にジャンプすることもできる。
【０１７３】
仮想ダイナミック像生成部１６は、生成したＶＤＣＥ−ＣＴ，ＶＤＣＥ−ＣＴＨＡ，ＶＤＣＥ−ＣＴＡＰ像の多断面再構成像（ＭＰＲ、Multi-planar Reconstruction)を作成してもよい。また、仮想ダイナミック像生成部１６は、生成したＶＤＣＥ−ＣＴ，ＶＤＣＥ−ＣＴＨＡ，ＶＤＣＥ−ＣＴＡＰ像の曲面のＭＰＲ像であるＣＰＲ（Curved Planar Reconstruction）を作成してもよい。さらに、また、仮想ダイナミック像生成部１６は、生成したＶＤＣＥ−ＣＴ，ＶＤＣＥ−ＣＴＨＡ，ＶＤＣＥ−ＣＴＡＰ像の最大値投影像（ＭＩＰ、Maximum Intensity Projection）、最小値投影像(ＭＩＮＩＰ、Minimum Intensity Projection)、または、平均値投影像を作成してもよい。さらに、仮想ダイナミック像生成部１６は、生成したＶＤＣＥ−ＣＴ，ＶＤＣＥ−ＣＴＨＡ，ＶＤＣＥ−ＣＴＡＰ像のボリュームレンダリング像を作成してもよい。以上のような多断面再構成像、ＣＰＲ、最大値投影像、最小値投影像、平均値投影像、またはボリュームレンダリング像において、時刻スライダー２２２により任意の時刻を指定して、任意の時刻の各画像を表示することにより、様々な視点からの診断に寄与できる。
【０１７４】
また、ビューワ２００においては、プロファイルやヒストグラムなどの表示部２１１に表示された画像の画像特徴量を計測する計測部を備えてもよい。このような計測部を備えることで、より詳細な解析が可能である。
【０１７５】
以上のように、本実施形態によれば、コンパートモデル式の係数を画素ごとに計算するのではなく、造影強度の時間変化とコンパートモデル式の係数を１対１に対応させたデータベースを作成し、作成したデータベースを検索することにより、コンパートモデル式の係数を得る。したがって、日常の検査で使用されるダイナミックＣＴ検査から得られる３相画像（造影前、早期相、後期相）から、対象患者の造影血流動態（大動脈と門脈の時間造影強度曲線を高精度かつ高速に推定することができる。
【０１７６】
動態解析を全画素ごとに行う場合には、例えば、縦が５１２画素、横が５１２画素だとすると、５１２×５１２で２６２１４４回の計算する必要がある。しかしながら、本実施形態の場合には、例えばシミュレーション点が上述したように、１０×１０の１００個だとすると、１００回の計算で済むことになる。その結果、全画素ごとに計算を行う場合に比べて、約２６２１倍高速に処理することができる。さらに体軸方向に複数断面の解析を行う場合，さらなる高速化が期待できる．例えば、縦が５１２画素、横が５１２画素，体軸方向に１００断面の解析を行う場合、５１２×５１２×１００で２６２１４４００回の計算する必要がある。しかしながら、本実施形態の場合には、例えばシミュレーション点が上述したように、１０×１０の１００個だとすると、１００回の計算で済むことになる。その結果、全画素ごとに計算を行う場合に比べて、約２６２１００倍高速に処理することができる。
【０１７７】
また、本実施形態によれば、造影強度パターンごとのシミュレーション点を設定し、全シミュレーション点について計算を行う場合でも、十分に滑らかな仮想ダイナミックＣＴ像を得ることができ、かつ、ノイズに対しても頑健になるという利点がある。
【０１７８】
また、本実施形態においては、コンパートメントモデル式を解くにあたり、上述した脾臓ステップ、肝臓ステップ、最終計算ステップからなる３ステップ法を適用することにより、未知変数を減らしつつ計算精度を保つことができる。
【０１７９】
また、本実施形態によれば、造影血流の流路となる主要血管を臨床的な事前知識と、複数患者の血流動態データから得られた血流動態パラメータを使用するため、血管の高精度な造影血流動態から全画素に存在する組織の血流動態および時間造影強度曲線を高精度に推定し、任意の時刻の仮想的なダイナミックＣＴ像を提示することができる。本実施形態では、入力血管の時間造影強度曲線の精度が高いため、組織各点の時間造影強度曲線の精度も高くなる。
【０１８０】
さらに、本実施形態は、血流動態の臨床的な事前知識と、複数患者から構築した血流動態パラメータを利用することで、患者固有の血流動態を高精度に推定することができる。それだけでなく、全画素に存在する組織の血流動態及び時間造影強度曲線を推定に使用する組織の血流変化モデルは、肝臓の大動脈と門脈の２つの血管の流入を考慮したモデルであるため、仮想的な肝動脈造影下のＣＴ像（ＣＴＨＡ）や門脈造影下のＣＴ像（ＣＴＡＰ）を計算し、提示することが可能である。
【０１８１】
本実施形態では、仮想ダイナミックＣＴ像として表示することで、最適タイミングと予想される時刻の仮想的な画像を確認することができ、再検査による患者の被爆量増加を防ぐことができる。
【０１８２】
（第２実施形態）
次に、本発明の第２実施形態について説明する。第１実施形態では、３相撮影された対象患者の医用画像に基づいて、薬物動態解析モデル係数データベース４２を作成する態様について説明したが、本実施形態では、例えば、１０相撮影された対象患者の医用画像を用いる。
【０１８３】
図６１は、大動脈と門脈の全時間造影強度曲線における早期相および平衡相と、コンパートモデル式のパラメータとの関係を説明するための図である。図６１に示すように、早期相のタイミングは、大動脈の造影信号強度が大きく、門脈の造影血流も肝臓へはあまり流入していないタイミングである。したがって、早期相のタイミングにおいては、パラメータＫ１ａを強く反映している。
【０１８４】
一方、平衡相のタイミングは、造影剤が全身に拡散し、大動脈と門脈の造影信号強度が均一になり、また臓器組織から造影剤が流出するタイミングである。したがって、平衡相のタイミングにおいては、分布容積Ｖｄは、（（Ｋ１ａ＋Ｋ１ｐ）／ｋ２）を強く反映している。
【０１８５】
したがって、早期相と平衡相との２相を使った解析は、多相撮影された全ての造影強度を使った解析の良い近似となることがわかる。
【０１８６】
そこで、本実施形態においては、１０相などの多相撮影された対象患者の医用画像を用いる場合であっても、多相全相を使うのではなく、造影前、早期相、後期相（または平衡相）の３相を用いて、薬物動態解析モデル係数データベース４２を作成する。
【０１８７】
また、１０相などの多相撮影された対象患者の医用画像を用いる場合には、時間分解能が十分と考えられるため、大動脈と門脈の全時間造影強度曲線（入力波形）は、撮影時刻の情報をそのまま使用する。したがって、第１実施形態において説明したＤＳＩによってＤＰＩを正規化する処理は行わない。ただし、それ以外に時間造影強度曲線に対する補間的処理を行っても構わない。
【０１８８】
本実施形態においては、１０相などの多相撮影された対象患者の医用画像を用いる場合であっても、造影前、早期相、後期相（または平衡相）の３相を用いて薬物動態解析モデル係数データベース４２を作成するので、処理の高速化を図ることができる。
【０１８９】
次に、１０相などの多相撮影された対象患者の医用画像を用いる場合の、早期相と後期相（平衡相）の決定方法について説明する。
【０１９０】
後期相（平衡相）は、造影剤が全身に行き渡り、全身血流で希釈される状態なので、あまり厳密なタイミングは要求されない。平衡相は、放射線科医の経験的に、１２０［ｓｅｃ］以降と考えられている、本実施形態では、１２０［ｓｅｃ］以降の最も早い相を選択する。
【０１９１】
一方、早期相（動脈優位相）のタイミングは、厳密なタイミングが要求されるため、脾臓の造影効果を用いて決定する。本実施形態では、造影剤が脾動脈から流入し，脾臓毛細血管を満たすことによって脾臓の血流量を最も正確に反映した最大造影効果が得られることに着目する。
【０１９２】
図６２は、脾臓の全時間造影強度曲線の一例を示す図である。本実施形態では、図６２のような脾臓の信号値の１秒当たりの増加率（傾き）が最大の相に続く一定時間内の次の相を早期相とする。また、これは患者の血流動態を鑑みて決めてよく、例えば血流循環が早い患者においては、最大の相としてもよい。
【０１９３】
図６３は、脾臓の全時間造影強度曲線における傾きを説明する図である。図６３において、青矢印で示す相が、傾きが最大となる相である。本実施形態では、傾きが最大となる相の次の相（赤矢印）を早期相として選択する。
【０１９４】
傾きが最大となる相の次の相を早期相として選択する理由は、傾きが最大の相は、脾動脈から脾臓の毛細血管へ最大の造影剤量が流入過程にある相と考えられ、この傾きが最大の相の次の相であれば、脾臓の毛細血管に最も多く造影剤が分布する相を撮影していると考えられるためである。そのため、早期相タイミングを決定する注目臓器は、解析対象と関連する臓器であればよく、例えば膵臓など動脈でのみ灌流される他の部位であってもよい。
【０１９５】
以上のように、本実施形態によれば、１０相などの多相撮影された対象患者の医用画像を用いる場合であっても、多相全相を使って薬物動態解析モデル係数データベース４２を作成するのではなく、造影前、早期相、後期相（または平衡相）の３相を用いて薬物動態解析モデル係数データベース４２を作成する。したがって、動態解析における処理の高速化を図ることができる。
【０１９６】
（変形例）
上述した実施形態においては、３相撮影されたＣＴ画像について本発明を適用した態様について説明したが、本発明はこのような態様に限定される訳ではない。例えば、ＭＲＩ、核医学検査を含めた医用画像にも拡張可能である。
【０１９７】
上述した実施形態においては、早期相と後期相のシミュレーション範囲を９等分し、それぞれの相において１０個のシミュレーション点を配置する態様について説明したが、本発明はこのような態様に限定される訳ではない。人間の目でＣＴ値の違いを判別できるのは、１０〜２０［ＨＵ］の範囲であることが知られており、それぞれの相におけるシミュレーション点の配置も、この範囲内の個数であればよい。
【０１９８】
また、上述した実施形態においては、各シミュレーション点を等間隔に配置する態様について説明したが、本発明はこのような態様に限定される訳ではない。各シミュレーション点の配置間隔は等間隔である必要はなく、例えば、機械学習等により、観測値に対応するシミュレーション点として、所定の範囲のシミュレーション点が特定される傾向がある場合には、当該範囲のシミュレーション点の間隔を他の範囲のシミュレーション点の間隔よりも細かくすることが考えられる。
【０１９９】
上述した実施形態においては、遺伝的アルゴリズム及び非線形最適化を用いて対象患者の所定の組織の薬物動態解析モデルの係数の最適な係数を算出する態様について説明したが、遺伝的アルゴリズムは、進化的戦略アルゴリズム等の最適化に代替可能であり、本発明はこのような態様に限定される訳ではない。線形および非線形的最適化を用いて対象患者の所定の組織の薬物動態解析モデルの係数の最適な係数を算出することもできる。
【０２００】
上述した実施形態においては、薬物動態解析モデルの一例として、コンパーメントモデルを用いた態様について説明したが、本発明はこのような態様に限定される訳ではない。例えば、コンパーメントモデル以外にも、Ｍａｘｉｍｕｍ  Ｓｌｏｐｅ法などの方法を用いることができる。
【０２０１】
上述した実施形態においては、細胞外液性造影剤を用いた態様について説明したが、本発明はこのような態様に限定される訳ではない。本発明は、例えば、組織特異性造影剤や血液プール造影剤など，他の種類の造影剤を用いた薬物動態解析にも応用可能である。
【０２０２】
上述した実施形態では、肝血流の動態解析に本発明を適用した例について説明したが、解析対象は肝臓に限らず、肝臓・脾臓・膵臓・腎臓・副腎・膀胱・前立腺・子宮・卵巣・精巣・消化菅・筋肉・リンパ節などの腹部臓器に適用可能である。また、本発明は、腹部臓器のみならず，腹部脈管が撮像範囲に含まれていれば、甲状腺・肺・心臓・脳といった全身の臓器に応用可能である。
【０２０３】
本発明においては、上述したように、全画素でＣｅｆｓ（ｔ）のベクトルや配列がメモリに保持されている状態であるので、薬物動態解析手法を問わず、指定された任意時刻に対応する値を、線形補間等の補間処理で取得することで、ＣＴ値マップ画像を生成することができる。
【０２０４】
上述したように、ＤＳＩは、循環速度を表す指標である。実際に心不全を評価する基準としては、「循環時間の遅延」という臨床所見があるため、循環時間と比例関係のある循環速度の低下は、心不全の診断として有用な指標になる。これを上述したビューアーで表示すると、心不全の診断に寄与することができる。また、上述したように、ＤＰＩは、ＤＳＩと密接な関係があるため、ＤＰＩとＤＳＩとを併記して表示することで、診断精度向上につながる。
【０２０５】
上述した実施形態においては、肝血流動態解析装置１００とビューワ２００とを別体で構成し、ビューワ２００に表示部２１１を設けた態様について説明した。しかしながら、本発明はこのような態様に限定される訳ではなく、肝血流動態解析装置１００とビューワ２００とを一体に構成し、肝血流動態解析装置１００に表示部を設けるようにしてもよい。
【０２０６】
上述した実施形態においては、肝血流動態解析装置１００の記憶部４０を複数患者血管信号値データベース４１の記憶部とした態様について説明した。しかしながら、本発明はこのような態様に限定される訳ではなく、複数患者血管信号値データベース４１をクラウド上のストレージあるいはサーバに記憶させるようにしてもよい。
【０２０７】
上述した解析モデルにおいて、脾臓モデルの際には入力血管として考えられる脾動脈が、肝臓モデルの際には入力血管として考えられる肝動脈と門脈が、それぞれ時間造影強度曲線の補正の対象血管であったが、入力血管の組み合わせはこれに限定されるものではない。例えば、肺モデルにおいては、入力血管として肺動脈と気管支動脈が考えられるが、腹部脈管が撮像範囲に含まれていれば、本発明を用いて肺動脈や気管支動脈の時間造影強度曲線の補正ができ、仮想ダイナミック像の生成が可能である。したがって、ピーク造影強度の時刻にジャンプする表示の対象血管も脾動脈や肝動脈、門脈に限定されるものではなく、時間造影強度曲線が補正された血管にも適用可能である。したがって、差分値を計算する入力血管の組合せも肝動脈と門脈だけに限定されない。
【０２０８】
以上の実施形態は例示であり、この発明の範囲から離れることなく様々な変形が可能である。上述した複数の実施の形態は、それぞれ単独で成立し得るものであるが、実施の形態同士の組みあわせも可能である。また、異なる実施の形態の中の種々の特徴も、それぞれ単独で成立し得るものであるが、異なる実施の形態の中の特徴同士の組みあわせも可能である。
【符号の説明】
【０２０９】
１        肝血流動態解析システム
１０      制御部
１１      対象患者血管信号値測定部
１２      シミュレーション部
１２ａ    第１決定部
１２ｂ    第２決定部
１２ｃ    計算部
１３      検索部
１４      取得部
１５      時間造影強度推定部
１６      仮想ダイナミック像生成部
２１      
２２      時間造影強度推定部
４０      記憶部
４１      複数患者血管信号値データベース
４２      薬物動態解析モデル係数データベース
５０      通信処理部
１００    肝血流動態解析装置
２００    ビューワ
２０１    入力画面
２０２    部位ＲＯＩ指定ツール
２０３    部位ＲＯＩ指定ボタン
２０４    撮影時刻ボタン
２０５    ラジオボックス
２１１    表示部
２２０    出力画面
２２１    選択アイコン
２２２    時刻スライダー
３００    画像サーバ
【図１】
2020168210.tif 000006
【図２】
2020168210.tif 000007
【図３】
2020168210.tif 000008
【図４】
2020168210.tif 000009
【図５】
2020168210.tif 000010
【図６】
2020168210.tif 000011
【図７】
2020168210.tif 000012
【図８】
2020168210.tif 000013
【図９】
2020168210.tif 000014
【図１０】
2020168210.tif 000015
【図１１】
2020168210.tif 000016
【図１２】
2020168210.tif 000017
【図１３】
2020168210.tif 000018
【図１４】
2020168210.tif 000019
【図１５】
2020168210.tif 000020
【図１６】
2020168210.tif 000021
【図１７】
2020168210.tif 000022
【図１８】
2020168210.tif 000023
【図１９】
2020168210.tif 000024
【図２０】
2020168210.tif 000025
【図２１】
2020168210.tif 000026
【図２２】
2020168210.tif 000027
【図２３】
2020168210.tif 000028
【図２４】
2020168210.tif 000029
【図２５】
2020168210.tif 000030
【図２６】
2020168210.tif 000031
【図２７】
2020168210.tif 000032
【図２８】
2020168210.tif 000033
【図２９】
2020168210.tif 000034
【図３０】
2020168210.tif 000035
【図３１】
2020168210.tif 000036
【図３２】
2020168210.tif 000037
【図３３】
2020168210.tif 000038
【図３４】
2020168210.tif 000039
【図３５】
2020168210.tif 000040
【図３６】
2020168210.tif 000041
【図３７】
2020168210.tif 000042
【図３８】
2020168210.tif 000043
【図３９】
2020168210.tif 000044
【図４０】
2020168210.tif 000045
【図４１】
2020168210.tif 000046
【図４２】
2020168210.tif 000047
【図４３】
2020168210.tif 000048
【図４４】
2020168210.tif 000049
【図４５】
2020168210.tif 000050
【図４６】
2020168210.tif 000051
【図４７】
2020168210.tif 000052
【図４８】
2020168210.tif 000053
【図４９】
2020168210.tif 000054
【図５０】
2020168210.tif 000055
【図５１】
2020168210.tif 000056
【図５２】
2020168210.tif 000057
【図５３】
2020168210.tif 000058
【図５４】
2020168210.tif 000059
【図５５】
2020168210.tif 000060
【図５６】
2020168210.tif 000061
【図５７】
2020168210.tif 000062
【図５８】
2020168210.tif 000063
【図５９】
2020168210.tif 000064
【図６０】
2020168210.tif 000065
【図６１】
2020168210.tif 000066
【図６２】
2020168210.tif 000067
【図６３】
2020168210.tif 000068
【手続補正書】
【提出日】令和1年7月19日(2019.7.19)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
少なくとも３相撮影された患者の医用画像に基づいて血流動態解析を行う際に、薬物動態解析モデルの入力となる血管の信号強度曲線を何らかの手法で推定する時間造影強度曲線推定部が動作した場合に、解析対象の臓器に関連する部位の信号強度に基づいて定めた早期相と後期相（または平衡相）を選択する相選択部と、選択された早期相と後期相の組み合わせからなる複数のシミュレーション点に対応付けて、最適化計算に基づいて算出した薬物動態解析モデルの最適係数を算出し、薬物動態解析モデル係数データベースとして保存するシミュレーション部と、少なくとも３相撮影された画像の対応する位置の画素において早期相と後期相の組み合わせに対応する前記シミュレーション点を検索する検索部と、前記検索部により検索した前記シミュレーション点に対応付けられた前記係数を取得する取得部と、前記検索部により検索した前記係数に基づいて、前記対象患者における所定の画素の時間造影強度曲線を推定する時間造影強度曲線推定部と、を備える、ことを特徴とする血流動態解析システム。
【請求項２】
前記時間造影強度曲線推定部は、少なくとも３相撮影された患者の医用画像において、薬物動態モデルの入力となる血管位置に設定した関心領域の画素値から推定されるものである、
ことを特徴とする請求項１に記載の血流動態解析システム。
【請求項３】
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースと、
少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定する対象患者血管信号値測定部と、を備え、
前記時間造影強度曲線推定部は、
前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線と、前記基準正規化信号強度曲線と前記対象患者の血管との関係に関するパラメータと、前記対象患者血管信号値測定部で得られた対象患者の血管信号値とから、前記対象患者の前記所定の血管の造影剤の到達距離の指標と循環速度の指標を計算し、前記所定の血管の時間造影強度曲線を推定する、
ことを特徴とする請求項１に記載の血流動態解析システム。
【請求項４】
前記相選択部において、３相の場合は、撮影時刻に基づきそのまま早期相、後期相を選択する、
ことを特徴とする請求項１ないし請求項３のいずれか一項に記載の血流動態解析システム。
【請求項５】
前記相選択部において、多相の場合は、解析対象の臓器が肝臓の場合、関連する部位として脾臓を用いて相選択を行う、
ことを特徴とする請求項１ないし請求項３のいずれか一項に記載の血流動態解析システム。
【請求項６】
前記シミュレーション部は、
脾臓のコンパートメントモデル式を用いて、動脈に関する時定数パラメータを決定する第１決定部と、
肝臓部のコンパートメントモデル式を用いて、門脈に関する時定数パラメータを決定する第２決定部と、
前記肝臓のコンパートメントモデル式における残りの係数を得るために、造影剤の各血管からの流入速度および注目組織からの流出速度を計算する計算部と、を備え、
前記複数のシミュレーション点に対する前記薬物動態解析モデルの最適な係数の算出を行う、
ことを特徴とする請求項１に記載の血流動態解析システム。
【請求項７】
前記シミュレーション部は、
前記基準正規化信号強度曲線を、大動脈および門脈の全時間造影曲線に基づいて算出し、前記複数のシミュレーション点を、前記大動脈および門脈の全時間造影曲線に基づいて定める、
ことを特徴とする請求項１または請求項２に記載の血流動態解析システム。
【請求項８】
前記シミュレーション部は、
前記薬物動態解析モデル係数データベースの作成における前記早期相の決定を、脾臓の信号強度を用いて行い、脾臓の信号強度の１秒当たりの増加率が最大の相の次の相を早期相とする、
ことを特徴とする請求項１ないし請求項３のいずれか１項に記載の血流動態解析システム。
【請求項９】
前記時間造影強度曲線推定部により推定した前記対象患者における所定の血管の時間造影強度曲線に基づいて、任意の時刻の仮想ダイナミック像を生成する仮想ダイナミック像生成部を備える、
ことを特徴とする請求項１ないし請求項４のいずれか１項に記載の血流動態解析システム。
【請求項１０】
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースを作成するステップと、
シミュレーション部により、前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出するステップと、
前記シミュレーション部により、前記複数のシミュレーション点と、前記薬物動態解析モデルの最適な係数とを対応付けて薬物動態解析モデル係数データベースを作成するステップと、
対象患者造影強度測定部により、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の早期相と後期相の造影強度を測定するステップと、
検索部により、前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索するステップと、
取得部により、前記検索部にて検索した前記シミュレーション点に対応付けられた前記係数を取得するステップと、
時間造影強度曲線推定部により、前記取得部にて取得した前記係数に基づいて、前記所定の血管の時間造影強度曲線を推定するステップと、を備える、
ことを特徴とする血流動態解析方法。
【請求項１１】
血流動態解析装置のプログラムであって、前記プログラムは、コンピュータに、
多相撮影された複数患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の画素値を測定し保存した複数患者血管信号値データベースであって、所定の血管における基準正規化信号強度曲線と、当該基準正規化信号強度曲線と前記所定の血管との関係に関するパラメータを保存する複数患者血管信号値データベースを作成するステップと、
シミュレーション部により、前記複数患者血管信号値データベースに保存された前記基準正規化信号強度曲線に基づいて定めた早期相と後期相の信号強度の組み合わせからなる複数のシミュレーション点に対応付けて、線形および非線形的最適化に基づいて算出した薬物動態解析モデルの最適な係数を算出するステップと、
前記シミュレーション部により、前記複数のシミュレーション点と、前記薬物動態解析モデルの最適な係数とを対応付けて薬物動態解析モデル係数データベースを作成するステップと、
対象患者造影強度測定部により、少なくとも３相撮影された対象患者の医用画像に基づいて、造影血流動態の臨床的な事前知識から導き出された造影血流の流路となる血管及び臓器の位置に関心領域を設定し、当該関心領域の早期相と後期相の造影強度を測定するステップと、
検索部により、前記関心領域の早期相と後期相の造影強度の組み合わせに対応する前記シミュレーション点を検索するステップと、
取得部により、前記検索部にて検索した前記シミュレーション点に対応付けられた前記係数を取得するステップと、
時間造影強度曲線推定部により、前記取得部にて取得した前記係数に基づいて、前記所定の血管の時間造影強度曲線を推定するステップと、を実行させる、
ことを特徴とするプログラム。