(57)【要約】【課題】構造流体解析結果を、血管の異常診断に容易に適用可能とする。【解決手段】注入プロファイル取得部が、被検体の血管に注入された造影剤の注入プロファイルを取得し、時系列ＣＴ値取得部が、造影剤注入後の血管内腔における時系列ＣＴ値を取得する。また、構造流体解析部が、注入プロファイルおよび時系列ＣＴ値を入力として構造流体解析シミュレーションを行う。第１算出部は、注入プロファイルおよび血管の異常の有無が判明している時系列ＣＴ値を用いて構造流体解析シミュレーションした結果から、構造流体解析シミュレーションと同じ入力と出力の関係を表すノミナルモデルを算出して記憶装置に記憶する。そして、差分モデル算出部が、血管の異常診断を行う被検体の時系列ＣＴ値と記憶装置に記憶されているノミナルモデルとの差分を検出することで、血管の異常の度合いを示す差分モデルを算出し、診断部が差分モデルを評価して血管の異常を診断する。【選択図】図２８
2015097759.tif 000002
【特許請求の範囲】
【請求項１】
被検体の血管に注入された造影剤の注入プロファイルを取得する注入プロファイル取得部と、
前記造影剤注入後の前記血管内腔における時系列ＣＴ値を取得する時系列ＣＴ値取得部と、
前記注入プロファイルおよび前記時系列ＣＴ値を入力として構造流体解析シミュレーションを行う構造流体解析部と、
前記注入プロファイルおよび血管の異常の有無が判明している時系列ＣＴ値を用いて前記構造流体解析シミュレーションした結果から、前記構造流体解析シミュレーションと同じ入力と出力の関係を表すノミナルモデルを算出する第１算出部と、
前記ノミナルモデルを記憶する記憶装置と、
血管の異常診断を行う被検体の時系列ＣＴ値と前記記憶装置に記憶されている前記ノミナルモデルとの差分を検出することで、血管の異常の度合いを示す差分モデルを算出する差分モデル算出部と、
前記差分モデルを評価して血管の異常を診断する診断部と
を有する血管解析装置。
【請求項２】
前記時系列ＣＴ値取得部で取得された、血管の異常診断を行う被検体の時系列ＣＴ値から、目標領域の時系列ＣＴ値を抽出する時系列ＣＴ値抽出部と、
前記時系列ＣＴ値抽出部で抽出された、目標領域の前記時系列ＣＴ値、および前記注入プロファイル取得部により取得された、前記血管の異常診断を行う被検体に対して注入された造影剤の注入プロファイルから、前記ノミナルモデルに対応する第２入出力モデルを算出する第２算出部を有し、
前記差分モデル算出部は、前記ノミナルモデルと前記第２入出力モデルとの差分を示す前記差分モデルを算出すること
を特徴とする請求項１に記載の血管解析装置。
【請求項３】
前記時系列ＣＴ値取得部は、血管の異常発生後の目標領域で撮影されたＣＴ画像から前記時系列ＣＴ値を取得すること
を特徴とする請求項１または請求項２に記載の血管解析装置。
【請求項４】
前記血管の異常部位の前後を含む血管の芯線方向にＣＴ値を取得し、血管の走行方向における造影剤濃度勾配を算出する濃度勾配算出部を有し、
前記時系列ＣＴ値として、算出された前記造影剤濃度勾配の時系列ＣＴ値を用いること
を特徴とする請求項１から請求項３のうち、いずれか一項に記載の血管解析装置。
【請求項５】
前記第１算出部および前記第２算出部のうち、少なくとも一方は、有限個のパラメータからなる線形差分方程式で表現されるモデルを算出すること
を特徴とする請求項１から請求項４のうち、いずれか一項に記載の血管解析装置。
【請求項６】
注入プロファイル取得部が、被検体の血管に注入された造影剤の注入プロファイルを取得する注入プロファイル取得ステップと、
時系列ＣＴ値取得部が、前記造影剤注入後の前記血管内腔における時系列ＣＴ値を取得する時系列ＣＴ値取得ステップと、
構造流体解析部が、前記注入プロファイルおよび前記時系列ＣＴ値を入力として構造流体解析シミュレーションを行う構造流体解析ステップと、
第１算出部が、前記注入プロファイルおよび血管の異常の有無が判明している時系列ＣＴ値を用いて前記構造流体解析シミュレーションした結果から、前記構造流体解析シミュレーションと同じ入力と出力の関係を表すノミナルモデルを算出して記憶装置に記憶する第１算出ステップと、
差分モデル算出部が、血管の異常診断を行う被検体の時系列ＣＴ値と前記記憶装置に記憶されている前記ノミナルモデルとの差分を検出することで、血管の異常の度合いを示す差分モデルを算出する差分モデル算出ステップと、
診断部が、前記差分モデルを評価して血管の異常を診断する診断ステップと
を有する血管解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施の形態は、血管解析装置、および血管解析方法に関する。
【背景技術】
【０００２】
三大疾患の一つである心疾患の原因となる冠動脈の狭窄、脳動脈瘤、またはそれらの予兆となる頸動脈のプラークによる狭窄を、非侵襲に予防、診断するための技術が望まれている。冠動脈の狭窄は虚血性心疾患に至る重大な病変である。心疾患ではステント治療を行うべきかどうか等を判断する必要があるが、冠動脈の狭窄診断としては、カテーテルによる冠動脈造影検査（CAG：Coronary  Angiography）が主流である。
【０００３】
冠動脈の器質的病変の診断指標として、心筋血流予備量比（FFR：Fractional  Flow  Reserve）があり、狭窄非存在下の最大冠血流に対する狭窄存在下の最大冠血流の比率で定義される。この血流比は、狭窄近位部冠内圧に対する狭窄遠位部冠内圧の比率とほぼ同じになり、カテーテル先の圧力センサで圧力比を測定する。
【０００４】
ここで、心臓のＸ線ＣＴ画像（ＣＴ：Ｃomputed  Ｔomography、以下、単にＣＴ画像という）で冠動脈の狭窄解析ができれば、カテーテル手術を必要とせず、低侵襲で患者負担低減となり、医療コスト節約にもつながる。現状はＣＴ画像によるプラークまたは血栓の大きさに基づく指標しか非侵襲な判断材料がない。
【０００５】
ＣＴ画像から血栓前後の圧力差等を構造流体シミュレーションにより解析できれば、血栓（またはプラーク）の及ぼす影響の定量化が期待できる。そして、得られた解析結果から、例えば狭窄の程度（狭窄度）等の血管の異常を診断できれば医師や患者にとって大きなメリットとなる。
【０００６】
しかし、解析結果には、各時刻における流速・流量・圧力データ等の大量のデータを含んでいる。このため、血管の異常診断用の画像を得るには、多大な計算リソースおよび多大な計算時間が必要となる。このようなことから、血管の異常診断に対して、上述の解析結果は、利用し難いものとなる恐れがある。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００８−２４１４３２号公報
【非特許文献】
【０００８】
【非特許文献１】門岡ら（ＩＴＵジャーナル，心臓シミュレータが拓くテーラーメード医療〜世界最先端の心臓シミュレータとその適用例のご紹介〜，Vol.41,No.6）
【非特許文献２】Ｊａｍｅｓ  Ｋ．Ｍｉｎ．，ｅｔ．ａｌ．，“Ｒａｔｉｏｎａｌｅ  ａｎｄ  ｄｅｓｉｇｎ  ｏｆ  ｔｈｅ  ＤｅＦＡＣＴＯ（Ｄｅｔｅｒｍｉｎａｔｉｏｎ  ｏｆ  Ｆｒａｃｔｉｏｎａｌ  Ｆｌｏｗ  Ｒｅｓｅｒｖｅ  ｂｙ  Ａｎａｔｏｍｉｃ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｉｃ  ＡｎｇｉＯｇｒａｐｈｙ）ｓｔｕｄｙ”，Ｊｏｕｒｎａｌ  ｏｆ  Ｃａｒｄｉｏｖａｓｃｕｌａｒ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ（２０１１）５，３０１−９
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明が解決しようとする課題は、構造流体シミュレーションによる解析結果を、血管の異常診断に容易に適用可能とした血管解析装置、および血管解析方法を提供することである。
【課題を解決するための手段】
【００１０】
実施の形態によれば、注入プロファイル取得部が、被検体の血管に注入された造影剤の注入プロファイルを取得し、時系列ＣＴ値取得部が、造影剤注入後の血管内腔における時系列ＣＴ値を取得する。また、構造流体解析部が、注入プロファイルおよび時系列ＣＴ値を入力として構造流体解析シミュレーションを行う。第１算出部は、注入プロファイルおよび血管の異常の有無が判明している時系列ＣＴ値を用いて構造流体解析シミュレーションした結果から、構造流体解析シミュレーションと同じ入力と出力の関係を表すノミナルモデルを算出して記憶装置に記憶する。そして、差分モデル算出部が、血管の異常診断を行う被検体の時系列ＣＴ値と記憶装置に記憶されているノミナルモデルとの差分を検出することで、血管の異常の度合いを示す差分モデルを算出し、診断部が、差分モデルを評価して血管の異常を診断する。
【図面の簡単な説明】
【００１１】
【図１】図１は、実施の形態の医用画像診断装置（Ｘ線コンピュータ断層撮影装置）の概略的なハードウェア構成図である。
【図２】図２は、実施の形態の医用画像診断装置の構造流体解析の対象領域に関する力学モデルの一例を示す図である。
【図３】図３は、実施の形態の医用画像診断装置の構造流体解析処理の流れを示す図である。
【図４】図４は、実施の形態の医用画像診断装置の画像処理装置の機能ブロック図である。
【図５】図５は、血管の芯線に直交する断面を模式的に示す図である。
【図６】図６は、画像追尾処理に利用される血管芯線形態の時間的変化を示す図である。
【図７】図７は、時刻ｔと時刻ｔ＋Δｔとの間における追尾処理の一例を示す図である。
【図８】図８は、画像解析・追尾処理による血管芯線の曲げ変形および回転変位の算出例を示す図である。
【図９】図９は、実施の形態の医用画像診断装置の力学モデル構築部により構築される形状モデルの芯線に直交する断面を示す図である。
【図１０】図１０は、実施の形態の医用画像診断装置の力学モデル構築部による形状モデルへの強制変位履歴の割り当てを説明するための図である。
【図１１】図１１は、実施の形態の医用画像診断装置の力学モデル構築部による形状モデルへの強制変位履歴の他の割り当て方法を示す図であり、血管形状変形指標がねじれの場合における割り当て例を示す図である。
【図１２】図１２は、実施の形態の医用画像診断装置の力学モデル構築部による形状モデルへの強制変位履歴の他の割り当て方法を示す図であり、血管形状変形指標が曲げの場合における割り当て例を示すである。
【図１３】図１３は、実施の形態の医用画像診断装置の統計的同定部による、階層ベイズモデルおよびマルコフ連鎖モンテカルロ法による負荷条件（血管内の平均圧力）に関する事後分布算出および平均内圧の同定を説明するための図である。
【図１４】図１４は、実施の形態の医用画像診断装置の統計的同定部による、階層ベイズモデルおよびマルコフ連鎖モンテカルロ法による材料モデルパラメータに関する事後分布算出および材料モデルパラメータ（血管壁の等価弾性率）の同定を説明するための図である。
【図１５】図１５は、実施の形態の医用画像診断装置における、力学的指標の一つである内圧の空間分布の表示例を示す図である。
【図１６】図１６は、実施の形態の医用画像診断装置における、血液流量指標の一つである流速値の空間分布の表示例を示す図である。
【図１７】図１７は、実施の形態の医用画像診断装置における、左冠動脈起始部の血圧に関するグラフの表示例を示す図である。
【図１８】図１８は、実施の形態の医用画像診断装置における、ＬＣＸとＬＤＡとの分岐点付近の血圧に関するグラフの表示例を示す図である。
【図１９】図１９は、実施の形態の医用画像診断装置における、芯線方向に関する血圧変化に関するグラフの表示例を示す図である。
【図２０】図２０は、実施の形態の医用画像診断装置の力学モデル構築部による強制変位履歴の他の割り当て例を示す図である。
【図２１】図２１は、実施の形態の医用画像診断装置の力学モデル構築部による強制変位履歴の他の割り当て例を示す図である。
【図２２】図２２は、実施の形態の医用画像診断装置の力学モデル構築部による強制変位履歴の他の割り当て例を示す図である。
【図２３】図２３は、実施の形態の医用画像診断装置における繊維グループの変形を示す図である。
【図２４】図２４は、実施の形態の医用画像診断装置の肉厚円筒の力学モデルの直交断面を示す図である。
【図２５】図２５は、実施の形態の医用画像診断装置における、微小扇形要素の拡大図である。
【図２６】図２６は、実施の形態の医用画像診断装置における、狭窄モデル同定を説明するための図である。
【図２７】図２７は、実施の形態の医用画像診断装置における、直接推定結果と狭窄モデルによる推定結果を比較するための図である。
【図２８】図２８は、実施の形態の医用画像診断装置で血管の異常診断を行う場合における機能ブロック図である。
【図２９】図２９は、実施の形態の医用画像診断装置における、造影剤の注入パターンを示す図である。
【図３０】図３０は、狭窄の有無に応じた造影剤の濃度変化を示すグラフである。
【図３１】図３１は、造影剤濃度勾配を説明するためのグラフである。
【発明を実施するための形態】
【００１２】
実施の形態の血管解析装置および血管解析方法は、医用画像診断装置で生成された医用画像に含まれる血管領域を構造流体解析するためのコンピュータ装置に適用できる。このコンピュータ装置は、医用画像診断装置に組み込まれていても良いし、医用画像診断装置とは別体のワークステーション等であっても良い。以下、実施の形態の血管解析装置および血管解析方法が適用されたコンピュータ装置が組み込まれた医用画像診断装置を、図面を参照して詳細に説明する。
【００１３】
（概要）
冠動脈等の動脈の流体解析では、血管またはプラークの硬さの影響、心臓拍動が動脈流へ及ぼす影響、動脈入口および動脈出口の境界条件の影響、を適切に考慮する必要がある。
【００１４】
このため、実施の形態の医用画像診断装置は、以下の構造流体解析を行う。すなわち、医用画像診断装置は、医用３次元画像に時間軸を加えた医用４次元画像を画像解析して得られる３次元変形履歴情報、および画像追尾から得られる変位履歴情報から、動脈の血管壁の形状変形指標と変位拘束条件を算出する。また、医用画像診断装置は、硬さ算定領域では、血管壁表面のみ、血管壁の力学モデルの節点に、強制変位履歴による変位自由度拘束条件を与える。
【００１５】
また、医用画像診断装置は、硬さ算定領域以外の領域では、血管壁の力学モデルにおける表面と内部の節点に強制変位履歴を与え、変位拘束条件を与えた構造流体解析を行う。これにより、硬さ算定領域の材料モデルおよび入口、出口の境界条件を同定する構造流体解析を行うことができ、高精度な血管および血行動態の解析が可能となる。そして、医用画像診断装置は、このような高精度な構造流体解析で得られた正常時および狭窄時のノミナルモデルを作成し、診断対象データとのモデル同定から狭窄度を診断する。
【００１６】
具体的には、実施の形態の医用画像診断装置は、ＣＴ（Ｃomputed  Ｔomography）診断装置から血管の４次元画像を取得し、血管の画像追尾から得られた血管部の各点の移動ベクトルを強制変位履歴として算出する。また、医用画像診断装置は、画像追尾情報から血管内腔の容積変化、直径変化、形状変化、または位相変化等の物理量を算出する。また、医用画像診断装置は、算出した物理量を指標として、指定領域部の材料モデル、または境界条件を同定する。
【００１７】
医用画像診断装置は、このような画像追尾と構造流体解析と逆解析（統計的同定法）に基づく血管解析により、血行動態における血圧、血流、応力、またはひずみ等の狭窄等の血管病変の予防、診断に必要な指標を高精度に解析する。そして、医用画像診断装置は、狭窄の有無が判明している典型例を対象に解析を行い、低次の標準モデル（ノミナルモデル）を予め作成しておき、診断対象データとのモデル同定から狭窄度を診断する。
【００１８】
このような実施の形態の医用画像診断装置は、被検体をスキャンするための撮像機構を装備する如何なる種類の画像診断装置にも適用可能である。実施の形態の医用画像診断装置としては、例えばＸ線コンピュータ断層撮影装置（Ｘ線ＣＴ装置）、磁気共鳴診断装置、超音波診断装置、ＳＰＥＣＴ（Single  Photon  Emission  CT）装置、ＰＥＴ（Positron  Emission  Tomography）装置、および放射線治療装置等に適宜利用可能である。以下、説明を具体的に行うため実施の形態の医用画像診断装置は、Ｘ線コンピュータ断層撮影装置であるものとする。
【００１９】
図１は、実施の形態の医用画像診断装置（Ｘ線コンピュータ断層撮影装置）の概略的なハードウェア構成図である。図１に示すように、Ｘ線コンピュータ断層撮影装置は、ＣＴ架台１０とコンソール２０とを有する。ＣＴ架台１０は、コンソール２０の架台制御部２３の制御に従って、Ｘ線で被検体の撮像部位をスキャンする。撮像部位は、例えば心臓である。
【００２０】
ＣＴ架台１０は、Ｘ線管１１、Ｘ線検出器１３、およびデータ収集装置１５を有している。Ｘ線管１１とＸ線検出器１３とは、回転軸Ｚ回りに回転可能にＣＴ架台１０に装備されている。Ｘ線管１１は、造影剤が注入された被検体にＸ線を照射する。Ｘ線検出器１３は、Ｘ線管１１から発生され被検体を透過したＸ線を検出し、検出されたＸ線の強度に応じた電気信号を発生する。
【００２１】
データ収集装置１５は、Ｘ線検出器１３から電気信号を読み出してデジタルデータに変換する。１ビュー毎のデジタルデータのセットは、生データセットと呼ばれている。複数のスキャン時刻に関する時系列の生データセットは、非接触データ伝送装置（図示しない）によりコンソール２０に伝送される。
【００２２】
コンソール２０は、システム制御部２１を中枢として、架台制御部２３、再構成装置２５、画像処理装置２７、入力機器２９、表示機器３１、および記憶装置３３を有している。
【００２３】
架台制御部２３は、ユーザにより入力機器２９を介して設定されたスキャン条件に応じてコンソール２０内の各装置を制御する。
【００２４】
再構成装置２５は、生データセットに基づいて、被検体に関するＣＴ画像のデータを発生する。具体的には、まず、再構成装置２５は、生データセットに前処理を施して投影データセットを発生する。前処理としては、対数変換や不均一補正、キャリブレーション補正等が含まれる。次に、再構成装置２５は、投影データセットに画像再構成処理を施してＣＴ画像のデータを発生する。画像再構成アルゴリズムとしては、ＦＢＰ（filtered backprojection）法等の解析学的画像再構成法や、ＭＬ−ＥＭ（maximum likelihood expectation maximization）法やＯＳ−ＥＭ（ordered subset expectation maximization）法等の逐次近似画像再構成等の既存のアルゴリズムが採用可能である。
【００２５】
実施の形態において再構成装置２５は、時系列の投影データセットに基づいて、時系列のＣＴ画像のデータを発生する。ＣＴ画像は、造影剤により造影された血管に関する画素領域（以下、血管領域と呼ぶことにする。）を含んでいる。なお、ＣＴ画像は、ＣＴ値の２次元空間分布を表現するスライスデータであっても良いし、ＣＴ値の３次元空間分布を表現するボリュームデータであっても良い。以下、ＣＴ画像はボリュームデータであるとする。時系列のＣＴ画像のデータは、記憶装置３３に記憶される。
【００２６】
画像処理装置２７は、時系列のＣＴ画像に基づいて、力学モデルを構築して構造流体解析を実行する。画像処理装置２７の処理の詳細については後述する。
【００２７】
入力機器２９は、ユーザからの各種指令や情報入力を受け付ける。入力機器２９としては、キーボードやマウス、スイッチ等が利用可能である。
【００２８】
表示機器３１は、ＣＴ画像や構造流体解析結果等の種々の情報を表示する。表示機器３１としては、例えばＣＲＴディスプレイや、液晶ディスプレイ、有機ＥＬディスプレイ、プラズマディスプレイ等が適宜利用可能である。
【００２９】
記憶装置３３は、ハードディスク装置等の種々の記憶媒体により構成される。記憶装置３３は、時系列の投影データや時系列のＣＴ画像データ等の種々のデータを記憶する。例えば、記憶装置３３は、時系列のＣＴ画像データをＤＩＣＯＭ（digital  imaging  and  communications  in  medicine）規格に準拠した医用画像ファイル形式で記憶する。また、記憶装置３３は、外部機器により収集された医用データを時系列のＣＴ画像データに医用画像ファイル内において関連付けて記憶しても良い。
【００３０】
システム制御部２１は、中央演算処理装置（ＣＰＵ：central  processing  unit）や読み出し専用メモリ（ＲＯＭ：read  only  memory）、ランダムアクセスメモリ（ＲＡＭ：random  access  memory）を有する。システム制御部２１は、Ｘ線コンピュータ断層撮影装置の中枢として機能する。システム制御部２１は、ＲＯＭやＲＡＭに記憶されている血管解析プログラムを実行して実施の形態の血管構造解析処理を実行する。
【００３１】
なお、システム制御部２１、画像処理装置２７、入力機器２９、表示機器３１、および記憶装置３３は、血管解析装置５０を構成する。実施の形態のように血管解析装置５０は、医用画像診断装置（Ｘ線コンピュータ断層撮影装置）に組み込まれていても良いし、医用画像診断装置とは別体のコンピュータ装置であっても良い。血管解析装置５０が医用画像診断装置とは別体の場合、血管解析装置５０は、医用画像診断装置やＰＡＣＳ（picturearchiving  and  communication  systems）からネットワークを介して時系列のＣＴ画像等の医用データを収集すれば良い。
【００３２】
次に、実施の形態の動作例を詳細に説明する。なお、実施の形態の医用画像診断装置は、心臓血管、頸動脈、または脳動脈等の人体のあらゆる部位の血管を解析対象とすることができる。以下、一例として、心臓の血管が解析対象であることとして説明を進める。
【００３３】
心臓の血管としては、例えば冠動脈と大動脈とが挙げられる。冠動脈は、大動脈の冠動脈起始部から始まり心筋表面を走行し、心外膜側から内膜側に入り込む。冠動脈は、心筋の内膜において無数の毛細管に分岐する。分岐後、無数の毛細管は、再び統合して大心静脈を形成し、冠静脈洞に接続する。冠血管系は、他の臓器と異なり、心筋の収縮および弛緩という力学的変化のなかで、灌流が保障されなければならないという点で特徴的である。
【００３４】
冠血流の特徴としては、心筋収縮による機械的血流阻害作用で冠動脈起始部の内圧が高くなる収縮期よりも、左心室拡張期に灌流圧が低下したときに多く流れることである。そのため、正常の冠動脈血流速波形は収縮期と拡張期の二峰性であり、拡張期血流が優位である。肥大型心筋症や大動脈弁狭窄症では、収縮期に逆行性波が認められ、大動脈逆流症では、収縮期順行波が大きくなる等、疾患によって特異な血流波形となることが知られている。また、拡張期の順行性波形は左室拡張機能、特に左室弛緩と密接な関係がある。左室弛緩遅延例では拡張期波形のピークが後ろにずれ、また減速脚がゆるやかになる傾向がある。また、このような症例では、頻拍時には拡張期の冠血流は十分に増大できず、心筋虚血を助長すると考えられている。
【００３５】
解剖学的に大動脈起始部から分岐する左右冠動脈に、大動脈圧に等しい冠灌流圧（すなわち、冠動脈が分枝する大動脈起始部の圧力）がかかることにより、冠血流が生じる。冠血流を決定するのは、大動脈圧である駆動圧と共に冠血管抵抗が重要である。１４０〜１８０μｍ以上の太い冠血管には、冠血管低抗の２０％程度が存在するのに対し、１００〜１５０μｍ以下の微小血管には、抵抗成分の残りの多くが存在するといわれる。従って、いわゆる冠狭窄等の無い場合には、抵抗値は冠微小血管の緊張性（トーヌス）に左右される。
【００３６】
血管抵抗因子としては、血管特性、動脈硬化、管狭窄、血液粘性、機械的因子があげられる。冠微小血管のトーヌスは、血管特性、心筋代謝（心筋酸素消費）、神経体液性因子、機械的因子、体液因子としての各種の血管作動性物質、血液粘性に規定され、さらに、心肥大、冠動脈硬化等を含めた様々な病変によっても影響され冠循環障害を起こす。
【００３７】
冠動脈血流拍動は、冠動脈血流の拍動パターン、心筋収縮による心筋内血流の制御、機械的刺激に対する心筋内血管の反応に影響される。心筋収縮が血流を阻害する機序としては、心筋内圧の上昇、心筋内血管容量の変化、心筋内血管の圧迫が挙げられる。心筋拡張期の血流規定因子には、拡張期の冠動脈圧、拡張期の血管外力、心拍数、心周期に占める拡張期の割合、心筋弛緩が存在する。
【００３８】
血管解析装置５０は、時系列のＣＴ画像に基づいて、力学モデルを構築し、力学モデルを利用して心臓の血管についての構造流体解析を実行し、血管内の力学的指標や血管流量指標を高精度に算出する。精度の高い力学的指標や血管流量指標を算出するためには、力学モデルに精度の高い潜在変数を割り当てる必要がある。血管解析装置５０は、力学モデルを構築する際、初期的な力学モデルに逆解析を施して潜在変数を統計的に同定する。これにより血管解析装置５０は、高精度の潜在変数を決定することができる。
【００３９】
力学的指標は、血管壁や血液に関する力学的な指標を意味する。血管壁に関する力学的指標としては、例えば血管壁の変位に関する指標、血管壁に生じる応力やひずみに関する指標、血管内腔に負荷される内圧分布に関する指標、血管の硬さ等をあらわす材料特性に関する指標等に分類される。血管の硬さ等をあらわす材料特性に関する指標は、血管組織の応力とひずみの関係を示す曲線の平均的な傾き等が挙げられる。
【００４０】
血液に関する力学的な指標における血管流量指標は、血管を流れる血液に関する血行動態の指標を意味する。血管流量指標としては、例えば血液の流量、血液の流速、血液の粘性等が挙げられる。なお、潜在変数は、例えば、血管の材料構成式、または血液の材料構成式といった材料モデルのパラメータ（例えばヤング率やポアソン比等）、血管内腔に負荷される内圧分布等の負荷条件パラメータ、構造解析や流体解析の境界条件パラメータ、時系列の形態指標や形状変形指標の不確定性に関連するばらつき分布パラメータの少なくとも一つを含む。
【００４１】
ここで、時系列の形態指標や形状変形指標の不確定性に関連するばらつき分布パラメータとは、医用画像データには、各ＣＴ値のノイズに起因したばらつき分布や、生体組織の境界閾値の曖昧性に起因する確率分布等が存在することを考慮して、血管組織や血液の境界座標および特徴点（血管分岐部や造影剤分散配置等）の空間座標における不確定性、または幾何学的構造パラメータ（芯線に垂直な断面の内腔半径等）の不確定性、または医用画像データ自体（ＣＴ値や境界閾値等）の不確定性を、確率分布として表現したものである。
【００４２】
力学モデルは、血管や血液の挙動を表現するための数値モデルである。力学モデルは、構造流体解析の手法に応じて異なるタイプを有している。例えば、力学モデルは、連続体力学モデルと簡易的力学モデルとに分類される。連続体力学モデルは、例えば、有限要素法（ＦＥＭ：finite element method）や境界要素法に用いられる。簡易的力学モデルは、例えば材料力学に基づく材料力学モデルと流れ学に基づく流体力学モデルとに分類される。
【００４３】
なお、以下の説明において特に言及しない場合、力学モデルのタイプについては特に限定しないものとする。初期的な力学モデルは、潜在変数の確率分布や変数範囲から得られる潜在変数のパラメータに関するサンプリング集合（各パラメータの組み合わせの集合）が割り当てられた力学モデルを意味するものとする。
【００４４】
図２は、構造流体解析の対象領域（以下、解析対象領域と呼ぶ）に関する力学モデルＭ１の一例を示す図である。図２に示すように、力学モデルＭ１は、大動脈領域Ｒ１と右冠動脈領域Ｒ２と左冠動脈領域Ｒ３とを有している。血液は、大動脈から右冠動脈または左冠動脈へ流れる。
【００４５】
図２に示すように、力学モデルＭ１の大動脈起始部側の末端が血流の入口に設定され、右冠動脈領域の末端と左冠動脈領域の末端とが血流の出口に設定される。入口と出口との各々に境界条件が設定される。入口に関する境界条件は、例えば入口における血流の流速または血流による圧力、またはそれらの変化率を含む。出口に関する境界条件は、例えば、出口における血流の流速または血流による圧力、またはそれらの変化率を含む。
【００４６】
大動脈、右冠動脈、および左冠動脈の変形は、血流に起因する血管壁への力学的作用、心臓の拍動による血管壁への力学的作用（外力）、血管断面境界の負荷条件、血管壁の材料モデル、血管の無応力状態、および血管壁の幾何学的形状等の様々な要因に依存する。
【００４７】
ここで、血流に起因する血管壁への力学的作用は、例えば、血流に起因する内圧と血流に起因するせん断力とを含む。血流に起因する内圧により、血管半径方向あるは血管内腔面の垂直方向に変形する。心臓の拍動による血管壁への力学的作用と、血流に起因するせん断力は、血管芯線方向に関する伸縮、ねじり、曲げといった血管壁への力学的作用による変形を生じさせる。
【００４８】
血管芯線方向に関する伸縮、ねじり、曲げといった血管の全体的および局所的な変形挙動は、大動脈領域Ｒ１、右冠動脈領域Ｒ２、および左冠動脈領域Ｒ３、またはその一部、またはその周辺組織に、強制変位（移動ベクトルや回転変位）または荷重ベクトルの時間的変化を与えることで、負荷条件を設定する。また、血流に起因する内圧に基づく血管半径方向または内腔面に垂直な方向への変形は、血管内腔への圧力分布の時間的変化として与えることで、負荷条件を設定する。
【００４９】
大動脈領域Ｒ１、右冠動脈領域Ｒ２、および左冠動脈領域Ｒ３には、構造流体解析において、強制変位による変位拘束条件が割り当てられる。これにより、構造流体解析における血管壁の変形自由度を縮小することができ、計算収束性の安定化や、解析時間の短縮を実現できる。
【００５０】
また、例えば、血管の形状の変形度合は、血管壁の材料に依存する。そのため、大動脈領域Ｒ１、右冠動脈領域Ｒ２、および左冠動脈領域Ｒ３に材料モデルが割り当てられる。また、例えば、血管の形状の変形度合は、血管の無応力状態に依存する。負荷条件の初期値として、血管の残留応力分布を負荷条件として与えても良い。
【００５１】
ここで、血管解析対象の数値計算用力学モデルの空間を離散化した節点集合および節点から構成される要素集合において、材料モデルや境界条件や負荷条件等の解析条件を同定しない領域と、同定する領域にわけ、解析条件を同定しない領域内の節点には強制変位履歴の変位拘束条件を与え、材料モデルを同定する領域では、血管壁表面（外表面）の節点のみ強制変位履歴の変位拘束条件を与え、血管壁内部は変位自由度を確保し変位拘束を与えない。これにより、構造流体解析の変形自由度をおさえることができ、安定的かつ効率的に解析を行うことができる。
【００５２】
ただし、血管壁表面上に緩衝用のダミーの要素集合を設け、その表面の節点に強制変位を与えてもよい。これにより、血管内腔にプラーク等による突起がある場合や、血管分岐部等、内圧による負荷が、芯線方向の断面外の変形に影響を与える場合には、血管内腔だけではなく血管壁の形態指標も参照して血管に負荷される荷重ベクトルと内圧を分離して同定することができる。また、生理学的には壁表面の脂肪層を模擬しており、一方、数値計算上は、血管壁表面に強制変位を与えることで血管壁内部に局所的に実際とは異なる高い応力が発生することを避ける効果もある。
【００５３】
これら材料モデル、境界条件、および負荷条件等の潜在変数に関するパラメータは、後述する力学モデルに基づく逆解析（統計的同定処理）により同定される。逆解析により同定された精確な潜在変数は、力学モデルに割り当てられる。精確な潜在変数が割り当てられた力学モデルにより、解析対象血管領域外の血管や心臓等の外部要因による当該解析対象血管領域への影響を加味した構造流体解析または流体解析または構造解析または画像解析に基づく血行動態解析を実行することができる。
【００５４】
血管解析装置５０は、力学モデルの構築に関し、逆解析による潜在変数の同定により、次の４点の困難を解決することができる。困難１．冠動脈の材料モデルの同定方法。困難２．心臓の形状の変形の冠動脈への影響の組み込み。困難３．冠動脈の境界条件の同定方法。困難４．医用画像データの不確定性に起因したばらつきを有する血管形状による画像解析や構造流体解析。この４点の困難の克服により、血管解析装置５０は、逆解析による潜在変数の同定を行わない従来の血管構造流体解析に比して、解析精度の向上を実現する。
【００５５】
次に、実施の形態の医用画像診断装置における構造流体解析処理の詳細について説明する。図３は、システム制御部２１の制御のもとに行われる構造流体解析処理の典型的な流れを示す図である。図４は、画像処理装置２７の機能ブロック図である。
【００５６】
図３に示すように、構造流体解析処理においては、まず、システム制御部２１により記憶装置３３から処理対象の医用画像ファイルが読み出され、画像処理装置２７に供給される。医用画像ファイルは、時系列のＣＴ画像のデータの他に、当該被検体に関する血管内腔に関する圧力値のデータ、血液流量指標の観測値のデータ、およびプラーク指標を含んでいる。ＣＴ画像以外ではＭＲＩ画像や超音波エコー画像であってもよい。時系列のＣＴ画像のデータは、時系列のＣＴ値の３次元空間分布を表現するデータである。時系列のＣＴ画像は、例えば１心拍で約２０枚、すなわち、約２０心位相分のＣＴ画像を含んでいる。
【００５７】
図３に示すように、システム制御部２１は、画像処理装置２７に領域設定処理を行わせる（ステップＳ１）。ステップＳ１において画像処理装置２７の領域設定部５１は、時系列のＣＴ画像に含まれる血管領域に構造流体解析の解析対象領域を設定する。解析対象領域は、冠動脈に関する血管領域の任意の一部分に設定される。例えば、領域設定部５１は、ユーザによる入力機器２９を介した指示、または、画像処理により血管領域に解析対象領域と同定対象領域を設定する。
【００５８】
ここで、図５を参照しながら、血管の構造について説明する。図５は、血管の芯線に直交する断面（以下、血管断面と呼ぶ）を模式的に示す図である。図５に示すように、血管は、管状の血管壁を有している。血管壁の中心軸は芯線と呼ばれている。血管壁の内側は内腔と呼ばれている。内腔に血液が流れる。内腔と血管壁との境は血管内壁と呼ばれている。血管壁の外側には心筋等の血管周辺組織が分布している。血管壁と血管周辺組織との境は血管外壁と呼ばれている。血管壁内部にプラークが発生することがある。
【００５９】
図５に示すように、プラークは、例えば、石灰化した石灰化プラーク、粥状プラーク等に分類される。粥状プラークは、やわらかく、血管内壁が破裂して血栓として血管内部に染み出す危険性があり、不安定プラークと呼ばれることもある。従って、プラークの性状を把握することは臨床的に有用である。プラークの性状や存在領域は、医用画像ファイルに含まれるプラーク指標により特定可能である。プラーク指標は、例えば骨のＣＴ値を基準に正規化したＣＴ値の大きさ等により相対的に判別することができる。しかし、血管内部のプラークの変形特性や硬さを解析するのは容易ではない。
【００６０】
ステップＳ１が行われるとシステム制御部２１は、画像処理装置２７に画像解析・追尾処理を行わせる（ステップＳ２）。ステップＳ２において画像処理装置２７の画像解析・追尾処理部５３は、時系列のＣＴ画像に画像処理を施して時系列の血管形態指標と時系列の血管形状変形指標とを算出する。具体的には、画像解析・追尾処理部５３は、時系列のＣＴ画像に画像解析処理を施すことにより時系列の血管形態指標を算出し、時系列のＣＴ画像に追尾処理を施すことにより時系列の血管形状変形指標を算出する。
【００６１】
より詳細には、画像解析・追尾処理部５３は、画像解析処理において、各ＣＴ画像から血管領域を抽出し、血管の内腔に関する画素領域（以下、血管内腔領域と呼ぶ）と血管壁に関する画素領域（以下、血管壁領域と呼ぶ）とを特定する。画像解析・追尾処理部５３は、血管形態指標として、血管の芯線に垂直な断面、または血管内腔面に垂直な面が、血管内腔、血管壁、プラーク領域に交わる領域上の複数の画素の３次元座標を特定する。
【００６２】
なお、血管形態指標は、３次元座標だけでなく、芯線に垂直な断面における一定角度ごとの血管内腔の半径や直径および０°の方向ベクトル、または断面における全角度に対する平均面積や平均半径、または、芯線方向に垂直な複数の断面で囲まれた血管内腔容積、または内腔面に垂直な複数断面で囲まれた血管壁容積やプラーク容積等の幾何学的指標でも良い。
【００６３】
追尾処理において、画像解析・追尾処理部５３は、ユーザからの入力機器２９を介した示または画像処理により、血管領域や血液や造影剤やプロトンにおける特徴点や特徴形状、代表点、画素等の複数の追跡点を設定する。例えば、画像解析・追尾処理部５３は、血管分岐部や表面の特徴形状等の追跡点集合を設定する。各時刻（各心位相）における画像解析・追尾処理部５３の追尾処理により得られた追跡点集合の変位データから、力学モデルの血管壁表面または血管壁内部または血管内腔における節点の変位の時間的変化を補間処理等により算出し、強制変位として与える。
【００６４】
また、例えば画像解析・追尾処理部５３は、力学モデルに血管芯線上の節点を定義する。画像解析・追尾処理部５３は、力学モデルの血管壁表面または血管壁内部または血管内腔における節点の変位の時間的変化から、血管の芯線方向に関する伸縮やねじりや曲げに関する変形を抽出し、血管芯線と芯線に垂直な断面における節点の強制変位として与えることで表現しても良い。このように、血管形状変形指標としては、力学モデルにおける各時刻の節点の強制変位データ（強制変位履歴）を特定する。
【００６５】
以下、図６および図７を参照しながら、画像解析・追尾処理を説明する。図６は、血管芯線の形態の時間的変化を示す図である。図６に示すように、例えば時系列の医用画像は、１心拍につき２０枚のＣＴ画像を含んでいるものとする。すなわち、心位相０％から９５％まで５％おきにＣＴ画像が得られることとする。画像解析・追尾処理部５３により血管領域の芯線が抽出される。図６に示すように、芯線の形態は、心位相の経過に従って変化している。
【００６６】
図７は、時刻ｔと時刻ｔ＋Δｔとの間における追尾処理の一例を示す図である。図７に示すように、血管芯線上にＰ１からＰ１０の力学モデルにおける節点が設定されており、各断面上の血管の力学モデルの節点と力学的につながっている。ただし、血液の力学モデルの節点とは独立である。血管の追跡点の変位データをもとに、血管芯線上のＰ１からＰ２０の節点の変位データを補間等の処理により算出し、各節点に強制変位が設定されているものとする。血管形状変形指標と血管形態指標について説明するため、節点Ｐ１３と節点Ｐ１４とにより規定される局所血管領域ＲＡを考える。
【００６７】
時刻ｔにおいて、芯線方向に関する節点Ｐ１３と節点Ｐ１４との間の距離がＬであり、血管領域の半径がｒであるとする。画像解析・追尾処理部５３から、節点Ｐ１３と節点Ｐ１４の血管芯線方向の伸縮、ねじり、または曲げ等の強制変位を抽出し、節点Ｐ１３の強制変位（３次元空間における移動変位と、芯線方向の回転変位）と節点Ｐ１４の強制変位（３次元空間における移動変位と、芯線方向の回転変位）を算出する。
【００６８】
図８は、血管芯線の曲げ変形および回転変位の算出例を示す図である。図８に示すように、例えばねじり角は、線ａおよび線ｂから構成される面の法線方向ベクトルcの変化から算出しても良い。
【００６９】
図７および図８に示すように、画像解析・追尾処理部５３は、追跡点の座標と移動ベクトルとに基づいて、芯線上の各節点の強制変位（３次元空間における移動変位と、芯線方向の回転変位）を算出し、血管形状変形指標を算出する。例えば、画像解析・追尾処理部５３は、隣り合う２つの節点の座標差の時間変化を芯線方向に関する伸縮距離ΔＬとして算出する。
【００７０】
また、画像解析・追尾処理部５３は、各芯線上の節点について、当該節点と当該節点を含む血管領域断面上の他の節点（血管内腔または血管壁またはプラーク領域における節点）との間の距離の時間変化を半径方向に関する伸縮距離Δｒとして算出する。また、画像解析・追尾処理部５３は、各追跡点について、当該追跡点の近傍の複数の追跡点の座標と移動ベクトルとに基づいて、芯線上の当該節点の芯線方向のねじれ角度Δθを算出する。
【００７１】
また、画像解析・追尾処理部５３は、血液領域の造影剤やプロトンの画像追尾により、血液流量指標として、流速または芯線方向断面の平均流速または平均流量を算出してもよい。
【００７２】
血管形状変形指標は、力学モデルにおける強制変位として利用される。以下、時系列の血管形態指標を「形状履歴」と呼び、時系列の血管形状変形指標を「強制変位履歴」と呼ぶことにする。
【００７３】
ステップＳ２が行われるとシステム制御部２１は、画像処理装置２７に構築処理を行わせる（ステップＳ３）。ステップＳ３において画像処理装置２７の力学モデル構築部５５は、形状履歴（時系列の血管形態指標）と、強制変位履歴（時系列の血管形状変形指標）と、時系列の医用画像（ＣＴ画像，ＭＲＩ画像，超音波エコー画像等のＤＩＣＯＭデータ）とに基づいて、解析対象領域に関する力学モデルを暫定的に構築する。力学モデルは、構造流体解析を行うための解析対象領域に関する数値モデルである。
【００７４】
以下、ステップＳ３について詳細に説明する。力学モデル構築部５５は、まず、医用画像と形状履歴とに基づいて、力学モデル（数理モデル）を解くための形状モデルを構築する。形状モデルは、各時刻における血管領域の幾何学的構造を模式的に表現したものである。形状モデルは、例えば複数の離散化領域に区分されている。各離散化領域の頂点は、節点と呼ばれる。
【００７５】
力学モデル構築部５５は、時刻毎の医用画像に含まれる血管領域と血管形態指標とに基づいて、時刻毎の形状モデルを構築しても良いし、特定の時相の医用画像に含まれる血管領域と血管形態指標とに基づいて、時刻毎の形状モデルを構築しても良い。また、例えば、初期の負荷状態として、解析対象領域に対応する血管に残留応力が存在しないと仮定する場合、無応力状態の時相として、解析対象領域に対応する血管が最も収縮した時相を無応力状態であると仮定する。
【００７６】
図９は、形状モデルの芯線に直交する断面を示す図である。図９に示すように、形状モデルは、芯線から外側に向けて血管内腔領域、血管壁領域を有している。プラークが存在する場合、血管壁領域にプラーク領域を設けても良い。また、血管周辺組織による血管への影響を考慮する場合、血管周辺組織のダミー要素を血管壁領域の外側に設けても良い。
【００７７】
形状モデルが構築されると力学モデル構築部５５は、各潜在変数の確率分布や変数範囲から得られる潜在変数のパラメータに関するサンプリング値（例えばマルコフ連鎖モンテカルロ法等による、各パラメータの組み合わせの集合からのサンプリング）を力学モデルに設定する。
【００７８】
例えば、力学モデル構築部５５は、図２に示すように、大動脈領域Ｒ１の大動脈起始部側の末端に入口に関する境界条件の同定対象の領域（以下、境界条件同定領域）を設定し、右冠動脈領域Ｒ２の末端と左冠動脈領域Ｒ３の末端とに出口に関する境界条件同定領域を設定する。力学モデル構築部５５は、各境界条件同定領域に境界条件の確率分布または変数範囲から得られる境界条件のパラメータに関するサンプリング値を割り当てる。
【００７９】
また、力学モデル構築部５５は、大動脈領域Ｒ１、右冠動脈領域Ｒ２、および左冠動脈領域Ｒ３に材料モデルの同定対象の領域（以下、材料モデル同定領域と呼ぶ）および負荷条件の同定対象の領域（以下、負荷条件同定領域と呼ぶ）を設定する。力学モデル構築部５５は、各材料モデル同定領域に材料モデルの確率分布または変数範囲から得られる材料モデルのパラメータに関するサンプリング値を割り当て、各負荷条件同定領域に負荷条件の確率分布や変数範囲から得られる負荷条件のパラメータに関するサンプリング値を割り当てる。
【００８０】
血管は、流量が０でも残留応力が存在すると言われている。例えば、力学モデル構築部５５は、流量が０の場合の残留応力を負荷条件の初期値として解析対象領域に割り当てても良い。また、力学モデル構築部５５は、幾何学的構造に不確定性がある部位に、幾何学的構造の同定対象の領域（以下、幾何学的構造同定領域と呼ぶ）を設定しても良い。なお、幾何学的構造のパラメータは、幾何学的構造の不確定性に関連するばらつき分布パラメータ、または医用画像データに内在するばらつき分布パラメータであり、各ＣＴ値のばらつき分布や生体組織の境界閾値のばらつき分布等であってよい。詳細は、後述するが、力学モデル構築部５５は、プラーク領域に材料モデルを設定しても良い。材料モデルの詳細については後述する。
【００８１】
形状モデルが構築されると力学モデル構築部５５は、ステップＳ２において算出された時系列の血管形状変形指標、すなわち、強制変位履歴を形状モデルに割り当てる。潜在変数および強制変位履歴が割り当てられた形状モデルを「力学モデル」と呼ぶことにする。
【００８２】
図１０は、形状モデルＭ２は血管や血液の力学モデルの一部を示しており、力学モデルの節点への強制変位履歴の割り当てを説明するための図である。図１０は、形状モデルＭ２の一部分を示している。ただし、図１０では芯線がＭ２内に位置する場合を示しているが、芯線がＭ２外に位置する場合でも良い。図１０に示すように、形状モデルＭ２には複数の節点ＰＮ（ＰＮ１，ＰＮ２）が設定されている。芯線上の節点をＰＮ１と称し、血管や血液の力学モデルにおける節点をＰＮ２と称することにする。形状モデルＭ２は、ダミー要素表面、血管外壁、血管内壁、プラーク領域表面、プラーク領域内部、または血液部、に設定される。力学モデル構築部５５は、形状モデルＭ２の各節点ＰＮ１に強制変位、すなわち、血管形状変化指標を時刻毎に割り当てる。
【００８３】
具体的には、力学モデル構築部５５は、芯線上の隣り合う節点ＰＮ１と節点ＰＮ１とをビーム要素（またはリジッド要素）ＥＢで結ぶ。また、力学モデル構築部５５は、節点ＰＮ１と当該節点ＰＮ１を通る直交断面に含まれる他の節点ＰＮ２とをビーム要素ＥＢで結ぶ。力学モデル構築部５５は、節点ＰＮ１およびビーム要素ＥＢに各血管形状変形指標の形状変位方向に関する拘束条件を割り当てる。
【００８４】
材料モデルまたは血管内腔の内圧を同定する領域では、強制変位としては、芯線方向に関する血管壁（またはダミー要素）表面の伸縮、血管壁（またはダミー要素）表面のねじれ、血管壁（またはダミー要素）表面の曲げ変形が挙げられる。例えば、材料モデルや血管内腔の内圧を同定しない領域では、芯線方向に関する強制変位だけではなく、半径方向に関する血管壁の時系列の伸縮（変位）も強制変位履歴として割り当てられる。
【００８５】
また、血管内腔に突起がある場合や血管分岐部等のように、内圧が、芯線方向断面外への変形に寄与する場合は、その領域には強制変位を与えずに、その領域の周辺部（例えば、ダミー要素の表面節点）のみに強制変位履歴を割り当てる。また、力学モデル構築部５５、節点ＰＮ１およびビーム要素ＥＢに時系列の血管形状変形指標を強制変位履歴として割り当てられる。これにより、血管の全体および局所に関する伸縮変形やねじれ変形や曲げ変形が表現される。
【００８６】
なお、強制変位履歴の割当対象は、芯線上の節点およびビーム要素に限定されない。図１１および図１２は、形状モデルへの強制変位履歴の他の割り当て方法を示す図である。図１１は、血管形状変形指標がねじれの場合における割り当て例を示している。また、図１２は、血管形状変形指標が曲げの場合における割り当て例を示している。図１１および図１２に示すように、力学モデル構築部５５は、形状モデルの表面または内部の節点に直接的に強制変位履歴を割り当てても良い。例えば、ステップＳ２において画像解析・追尾処理部５３は、特徴点についての伸縮量やねじり量や曲げ量等の血管形状変形指標を算出する。力学モデル構築部５５は、算出された血管形状変形指標を特徴点の周囲の節点へ補間（内挿や外挿）することにより直接的に割り当てる。
【００８７】
実施の形態の画像処理装置２７は、ステップＳ３において暫定的に構築された力学モデルを用いて逆解析を施し、力学モデルに設定される潜在変数を統計的に同定する。統計的同定処理は、後述のステップＳ６において行われる。ステップＳ４およびＳ５は、それぞれ統計的同定処理に用いる血管形態指標および血液流量指標を算出するために設けられている。
【００８８】
ステップ３が行われるとシステム制御部２１は、画像処理装置に２７に血管応力解析処理を行わせる（ステップＳ４）。ステップＳ４において画像処理装置２７の血管応力解析部５７は、現段階の力学モデルに血管応力解析を施して時系列の血管形態指標の予測値を算出する。血管形態指標は、既述の血管形態指の何れであっても良いが、例えば血管芯線方向に関する内腔領域の断面形状指標や血管壁の断面形状指標が用いられると良い。
【００８９】
具体的には、内腔領域の断面形状指標は、内腔領域の注目画素の座標値、内腔領域の幾何学的構造パラメータ（内腔領域の半径、内腔領域の直径等）の少なくとも一つである。また、血管壁領域の断面形状指標は、具体的には、血管壁領域の注目画素の座標値、血管壁領域の幾何学的構造パラメータ（血管壁領域の半径、壁領域の直径等）の少なくとも一つである。なお、予測値は、力学モデルに血管応力解析を施して算出された血管形態指標の算出値を意味する。
【００９０】
また、ステップ３が行われるとシステム制御部２１は、画像処理装置２７に血液流体解析処理を行わせる（ステップＳ５）。ステップＳ５において画像処理装置２７の血液流体解析部５９は、暫定的に構築された力学モデルに血液流体解析を施して時系列の血液流量指標の予測値を算出する。血液流量指標は、血流量または流速、またはその空間的・時間的な平均値の少なくとも一つである。なお、予測値は、力学モデルに血液流体解析を施して算出された血液流体指標の算出値を意味する。
【００９１】
ステップＳ４およびＳ５が行われるとシステム制御部２１は、画像処理装置２７に同定処理を行わせる（ステップＳ６）。ステップＳ６において画像処理装置２７の統計的同定部６１は、ステップＳ４において算出された血管形態指標の予測値とステップＳ５において算出された血液流量指標の予測値とが、事前に収集された血管形態指標の観測値と血液流量指標の観測値とに整合するように力学的モデルの潜在変数のパラメータを統計的に同定する。
【００９２】
図４に示すように、統計的同定部６１は、第１統計的同定部６１−１と第２統計的同定部６１−２とを有している。第１統計的同定部６１−１は、血管形態指標の予測値が血管形態指標の観測値に整合するように潜在変数のパラメータを統計的に同定する。第２統計的同定部６１−２は、血液流量指標の予測値が血液流量指標の観測値に整合するように潜在変数のパラメータを統計的に同定する。以下、第１統計的同定部６１−１と第２統計的同定部６１−２とを順番に説明する。
【００９３】
具体的には、ステップＳ６において第１統計的同定部６１−１は、ステップＳ４において算出された血管形態指標の予測値と観測値とに基づくデータ分布を設定する。データ分布は、例えば、血管形態指標の予測値と観測値との誤差に関する多変量正規分布関数を示す。
【００９４】
力学モデルにおける各節点または各要素での予測値と観測値の誤差に関する正規分布関数値を算出し、それらの積をデータ分布として統計的同定処理部において設定する。データ分布は、時刻毎に個別に設定されても良いし、複数時刻まとめて設定されても良い。
【００９５】
次に、第１統計的同定部６１−１は、力学モデルの潜在変数に事前分布（事前確率分布）を割り当てる。具体的には、材料モデル、境界条件、負荷条件、および時系列の形態指標や形状変形指標の不確定性に関連するパラメータの各々に事前分布が割り当てられる。例えば、負荷条件のパラメータの一つである血管内腔に関する圧力値に関する事前分布が割り当てられる。圧力値の取り得る値の範囲（想定範囲）は、経験的に予め限定することができる。第１統計的同定部６１−１は、これら想定範囲内に限定して内圧値に関するモンテカルロシミュレーションを実行することにより、各離散化領域について内圧値の確率分布、すなわち、事前分布を算出する。
【００９６】
また、力学モデル構築部５５は、事前分布として、芯線方向の圧力分布は滑らかであること、また、時間経過に伴う圧力変化が滑らかであること、血流の逆流がないことが観測されている場合には芯線方向の平均的な圧力変化の傾きが負であることを、例えば多変量正規分布関数により数学的に表現した確率分布を事前分布として設定しても良い。
【００９７】
これら想定範囲内に限定して、想定した確率分布に従って、負荷条件のパラメータに関するモンテカルロシミュレーションを実行することができ、力学モデルへ設定するための負荷条件（潜在変数）のサンプリング値を得ることができる。
【００９８】
次に、第１統計的同定部６１−１は、各潜在変数について、事前分布とデータ分布とに統計的同定処理を施すことにより事後分布（事後確率分布）を算出する。統計的同定処理は、例えば、階層ベイズモデルやマルコフ連鎖モデルが挙げられる。そして、第１統計的同定部６１−１は、各潜在変数について、事後分布の最頻値や平均値等の統計値から各潜在変数のパラメータを同定する。例えば、上述の例の場合、血管内腔圧力値に関する事後分布が算出され、この事後分布から血管内腔圧力値の同定値が算出される。
【００９９】
図１３は、階層ベイズモデルおよびマルコフ連鎖モンテカルロ法による負荷条件（血管内の平均圧力）に関する事後分布算出と平均内圧の同定とを説明するための図である。図１３に示すように、血管起始部から延びる血管に石灰化プラーク領域と粥状プラーク領域とが設定されているものとする。石灰化プラーク領域は材料モデル同定領域に設定され、粥状プラーク領域に材料モデル同定領域に設定される。血管起始部から血管芯線方向に沿って進むにつれて血管内圧が降下する。血管芯線に沿って複数の節点が設定される。各節点を含む直交断面（節点断面）において内腔内圧の事後分布が算出され、事後分布の最頻値が同定される。
【０１００】
なお、血管形態指標の観測値としては、例えば、ステップＳ２において算出された血管形態指標が用いられる。
【０１０１】
第２統計的同定部６１−２による処理は、データ分布の算出に用いる指標が異なるだけで第１統計的同定部６１−１による処理と同様である。すなわち、第２統計的同定部６１−２は、まず、ステップＳ５において算出された血液流量指標の予測値と観測値とに基づくデータ分布を設定する。
【０１０２】
次に、第２統計的同定部６１−２は、力学モデルの潜在変数に事前分布を割り当てる。例えば、血管に関する材料モデルのパラメータや血液に関する材料モデルのパラメータ、プラークに関する材料モデルのパラメータに関する事前分布が割り当てられる。これら材料モデルのパラメータとしては、例えば、弾性率等の材料モデルパラメータや、血液の構成式における粘性に関するパラメータが挙げられる。材料モデルのパラメータの想定範囲や確率分布は、経験的に予め設定することができる。
【０１０３】
第２統計的同定部６１−２は、各離散化領域について材料モデルのパラメータの確率分布、すなわち、事前分布を設定し、これら想定範囲内に限定して、想定した確率分布に従って、材料モデルのパラメータに関するモンテカルロシミュレーションを実行することができ、力学モデルへ設定するための材料モデルパラメータ（潜在変数）のサンプリング値を得ることができる。
【０１０４】
次に、第２統計的同定部６１−２は、各潜在変数について、事前分布とデータ分布とに統計的同定処理を施すことにより事後分布を算出し、算出された事後分布の統計値から各潜在変数のパラメータを同定する。例えば、上述の例の場合、材料モデルのパラメータに関する事後分布が算出され、この事後分布から材料モデルのパラメータの同定値が算出される。
【０１０５】
図１４は、階層ベイズモデルおよびマルコフ連鎖モンテカルロ法による材料モデルパラメータに関する事後分布算出と材料モデルパラメータ（血管壁の等価弾性率）の同定とを説明するための図である。図１４に示すように、血管モデルは、図１３と同様であるとする。材料モデル同定領域に限定して、血管壁の材料モデルのパラメータ（例えば、等価弾性率）の事後分布が算出され、事後分布の最頻値が同定される。
【０１０６】
なお、血液流量指標の観測値は、例えば、大動脈に送り出される血流量変化であると仮定し、血管形態指標の観測値を、時系列のＣＴ画像から画像処理により計測される左心室の容積変化値（ＣＦＡ）として用いることができる。造影剤の冠動脈内注入後の造影剤の画像追尾により特徴点の移動量の時間的変化を算出することにより、流速や流量を算出してもよい。また、造影剤の血管芯線方向または時間的な特定領域の濃度変化量を取得し、その濃度変化を各領域の芯線方向距離間隔で除した値や、濃度変化の時間的変化率から、流速や流量を算出してもよい。ＭＲＩの場合はプロトンの画像追尾を用い、超音波エコーの場合には、コントラストエコー図法等により流量を算出する。
【０１０７】
また、解析対象領域の各画素の座標値が確定値であることを前提としない場合、すなわち、解析対象領域の幾何学的構造に不確定性があると仮定する場合、潜在変数に幾何学的構造を含めても良い。この場合、統計的同定部６１は、各節点の座標値の芯線方向に関する所定範囲内の変動、または、解析対象領域の径の所定範囲内の変動を表現した正規分布等の確率分布を事前分布として設定すると良い。この場合、解析対象領域の形状が滑らかであること、また、芯線における節点の順番が不変であるという制約も事前分布として設定しても良い。これら想定範囲内に限定して、想定した確率分布に従って、幾何学的構造のパラメータに関するモンテカルロシミュレーションを実行することができ、力学モデルへ設定するための幾何学的構造の不確定性パラメータ（潜在変数）のサンプリング値を得ることができる。
【０１０８】
なお、各ステップＳ６において第１統計的同定部６１−１による統計的同定処理と第２統計的同定部６１−２による統計的同定処理との両方が行われなくても良い。すなわち、各ステップＳ６においては、第１統計的同定部６１−１による統計的同定処理と第２統計的同定部６１−２による統計的同定処理との何れか一方が行われても良い。
【０１０９】
また、上記の例においては、第１統計的同定部６１−１は、血管形態指標の予測値が血管形態指標の観測値に整合するように潜在変数のパラメータを統計的に同定し、第２統計的同定部６１−２は、血液流量指標の予測値が血液流量指標の観測値に整合するように潜在変数のパラメータを統計的に同定するとした。しかし、統計的同定部６１は、構造−流体連成解析に基づいて、血管形態指標の予測値と血液流量指標の予測値とが血管形態指標の観測値と血液流量指標の観測値とに整合するように潜在変数のパラメータを統計的に同定しても良い。統計的同定部６１による統計的同定処理のさらなる詳細については後述する。
【０１１０】
ステップＳ６が行われるとシステム制御部２１は、画像処理装置２７に設定処理を行わせる（ステップＳ７）。ステップＳ７において画像処理装置２７の力学モデル構築部５５は、ステップＳ６において算出された潜在変数のパラメータを力学モデルに設定する。
【０１１１】
ステップＳ７が行われるとシステム制御部２１は、同定終了条件が満たされたか否かを判定する（ステップＳ８）。ステップＳ８において同定終了条件が満たされていないと判定した場合（ステップＳ８：ＮＯ）、システム制御部２１は、ステップＳ４、Ｓ５、Ｓ６、Ｓ７、およびＳ８を繰り返す。
【０１１２】
ここで、同定終了条件は、同定終了を判定するための指標（以下、同定終了指標と呼ぶ）が規定値に達するか否かにより表現される。同定終了指標としては、例えば、血管形態指標の予測値と観測値との差分値が挙げられる。この場合、システム制御部２１は、この差分値が既定値よりも大きい場合、同定終了条件が満たされていないと判定し、差分値が既定値よりも小さい場合、同定終了条件が満たされたと判定する。
【０１１３】
また、同定終了指標は、例えば、モンテカルロ法のサンプリング点の数でも良い。この場合、システム制御部２１は、このサンプリング点の数が既定値よりも小さい場合、同定終了条件が満たされていないと判定し、サンプリング点の数が既定値よりも大きい場合、同定終了条件が満たされたと判定する。同定終了条件が満たされた場合、力学モデル構築部５５は、その時点の最新の力学モデルを最終的な力学モデルに設定する。
【０１１４】
最終的な力学モデルが構築されると、力学モデル構築部５５は、血管形状変形指標の観測値、最終的な力学モデルにおける負荷条件のパラメータ、および材料モデルのパラメータを関連付けたモデル（以下、関連モデルと呼ぶ）を算出する。関連モデルは、記憶装置３３に記憶される。関連モデルは、検索の容易性等のため患者情報や検査情報等に関連付けて記憶されると良い。なお、血管形態指標や血液流量指標の観測値、最終的な力学モデルにおける負荷条件のパラメータ、および材料モデルのパラメータは、必ずしもモデルの形態で関連付けられる必要はなく、例えば、テーブルまたはデータベースであっても良い。
【０１１５】
上記のステップＳ４、Ｓ５、Ｓ６、Ｓ７、およびＳ８は、同一の同定法で反復しても良いし、異なる同定法で反復しても良い。異なる同定法で反復する場合、例えば、まず、簡易的力学モデルを利用して潜在変数を暫定的に同定し、次に、連続体力学モデルを利用して潜在変数を正確に同定しても良い。このように統計的同定処理を異なる手法で２段階に分けて行うことにより、潜在変数のパラメータを短時間で収束させることができる。簡易的力学モデルを利用する方法としては、内圧および外圧を厚肉円筒の材料力学の式と、ハーゲン・ポアズイユ流れおよび修正ベルヌーイの式とを用いる方法が挙げられる。連続体力学モデルを利用する方法としては、ＦＥＭ構造流体解析が挙げられる。簡易的力学モデルを利用する同定法と連続体力学モデルを利用する同定法との詳細については後述する。
【０１１６】
ステップＳ８において同定終了受件が満たされたと判定した場合（ステップＳ８：ＹＥＳ）、システム制御部２１は、画像解析・追尾処理部５３に修正処理を行わせても良い（ステップＳ９）。ステップＳ９において画像解析・追尾処理部５３は、統計的同定法による逆解析で得られた潜在変数のもとで実施した構造流体解析結果（力学的指標の予測値および血液流体指標の予測値）が観測値（力学的指標の観測値および血液流体指標の観測値）に整合するように、時系列の医用画像に含まれる血管領域の形状を修正しても良い。
【０１１７】
表示機器３１は、修正後の時系列の医用画像に基づく診断結果を表示する。これにより、血管解析装置５０は、最終的な力学モデルを考慮した診断結果を表示することができる。または、表示機器３１は、逆解析による同定とその構造流体解析により観測結果とが整合しない血管箇所・領域を画面に表示しても良い。
【０１１８】
例えば、血管の挙動の動きの速い心位相の画像はぼやけることが多く、医用画像をもとにした画像解析により観測した血管形状には誤差が大きい箇所や領域が存在する。比較的、血管の挙動が安定した心位相のデータでは、ノイズが少ない画像から得られる。そのような誤差分布の小さい血管形状データに基づき、誤差が大きな心位相での血管形状を、力学モデルを用いることで、正しく内挿することができ、そのような誤差が大きい血管箇所や領域について、正しく内挿した形状とともに、元データからのばらつき分布とともに、表示することができる。これにより、血管形状表示の安定性を確保できるとともに、形状の不確定性をユーザが認識できる。
【０１１９】
ステップＳ９が行われるとシステム制御部２１は、画像処理装置２７に血管応力解析処理を行わせる（ステップＳ１０）。ステップＳ１０において画像処理装置２７の血管応力解析部５７は、最終的な力学モデルに血管応力解析を施して時系列の力学的指標の予測値の空間分布を算出する。具体的には、離散化領域毎に力学的指標の予測値が算出される。
【０１２０】
また、ステップＳ９が行われるとシステム制御部２１は、画像処理装置２７に血液流体解析処理を行わせる（ステップＳ１１）。ステップＳ１１において画像処理装置２７の血液流体解析部５９は、暫定的に構築された力学モデルに血液流体解析を施して時系列の血液流量指標の予測値の空間分布を算出する。具体的には、離散化領域毎に血液流量指標の予測値が算出される。
【０１２１】
なお、力学的指標または血液流量指標として、ＦＦＲが算出されても良い。
【０１２２】
ステップＳ９およびＳ１０が行われるとシステム制御部２１は、表示機器３１に表示処理を行わせる（ステップＳ１１）。ステップＳ１１において表示機器３１は、ステップＳ９において算出された時系列の力学的指標の予測値とステップＳ１０において算出された時系列の血液流量指標の予測値とを表示する。
【０１２３】
例えば、表示機器３１は、時系列の力学的指標または時系列の血管流量指標を、時系列の力学モデルを当該予測値に応じた色で動画的に表示する。このため、表示機器３１は、各種の予測値とカラー値（例えば、ＲＧＢ）との関係を示すカラーテーブルを保持している。表示機器３１は、予測値に応じたカラー値をカラーテーブルを利用して特定し、特定されたカラー値に応じた色で当該予測値に対応する離散化領域を表示する。
【０１２４】
図１５は、力学的指標の一つである内圧の空間分布の表示例を示す図である。図１５に示すように、表示機器３１は、力学モデルを構成する各離散化領域を当該離散化領域に関する内圧値に応じた色で動画的に表示する。ユーザは、力学モデルを観察することにより、経時的且つ空間的に変化する力学的指標を色で把握することができる。
【０１２５】
図１６は、血液流量指標の一つである流速値の空間分布の表示例を示す図である。図１６に示すように、表示機器３１は、力学モデルを構成する各離散化領域を当該離散化領域に関する流速値に応じた色で動画的に表示する。ユーザは、力学モデルを観察することにより、経時的且つ空間的に変化する血液流量指標を色で把握することができる。
【０１２６】
例えば、血管内が完全に狭窄している場合、狭窄部位の内圧は、非狭窄部位の内圧よりも小さい。力学的指標として内圧が指定された場合、ユーザは、力学モデル上での局所的な色の違いにより狭窄の有無を判断することができる。また、狭窄部位の流量は、非狭窄部位の流用よりも小さい。血液流量指標として流量が指定された場合、ユーザは、力学モデル上での局所的な色の違いにより狭窄の有無を判断することができる。
【０１２７】
また、血管応力解析部５７は、プラーク領域の材料モデルパラメータの同定結果に基づいて、力学的指標として硬さ値の空間分布を算出しても良い。この場合、表示機器３１は、プラーク領域に関する硬さ値の空間分布を力学モデル上において表示しても良い。また、表示機器３１は、プラーク領域周辺の内圧分布や応力分布、ひずみ分布を表示して良い。ユーザは、これらの表示をプラークの性状と破綻しやすさとを推定することに活用することができる。
【０１２８】
力学的指標および血液流量指標の予測値は、力学モデルの離散化領域の色で表現する方法のみに限定されない。例えば、図１７、図１８、および図１９に示すように、グラフで表示しても良い。なお、図１７は、左冠動脈起始部の血圧に関するグラフである。図１７のグラフの縦軸は正規化した血圧に規定され、横軸は心位相［％］に規定される。図１８は、ＬＣＸとＬＤＡとの分岐点付近の血圧に関するグラフである。図１８のグラフの縦軸は正規化した血圧に規定され、横軸は心位相［％］に規定される。図１９は、芯線方向に関する血圧変化に関するグラフである。図１９のグラフの縦軸は血圧比に規定され、横軸は大動脈からの距離［ｍｍ］に規定される。表示機器３１は、力学的指標および血液流量指標の予測値をグラフで表示することにより、これら値をユーザに簡便に把握させることができる。
【０１２９】
ステップＳ１１が行われると構造流体解析処理が終了する。
【０１３０】
なお、図１０において、強制変位履歴は、形状モデルの芯線部と外壁部とに設定されるとしたが、強制変位履歴の設定箇所は、これに限定されない。例えば、強制変位履歴は、芯線部と外壁部との間の血管壁領域に設定されても良い。
【０１３１】
また、強制変位履歴の拘束条件の割り当て対象は、境界条件および材料モデルを同定するか否かに応じて切り分けられても良い。図２０は、強制変位履歴の他の割り当て例を示す図であり、形状モデルの断面を示している。例えば、図２０の（ａ）に示すように、境界条件および材料モデルを同定する場合、形状モデルの外壁部ＯＷ上の節点ＰＮ２のみに強制変位履歴を割り当て、血管壁領域ＲＶの節点ＰＮ３には強制変位履歴を割り当てなければよい。
【０１３２】
また、図２０の（ｂ）に示すように、境界条件および材料モデルを同定しない場合、形状モデルの外壁部の節点ＰＮ２と血管壁領域ＲＶの節点ＰＮ３との両方に強制変位履歴を割り当てればよい。この場合、芯線上の節点ＰＮ１に強制変位履歴が割り当てられる。また、外壁部ＯＷの節点ＰＮ２と節点ＰＮ１とをビーム要素ＥＢで結び、ビーム要素ＥＢ上の節点ＰＮ２およびＰＮ３にも強制変位履歴を割り当ててもよい。このとき、半径方向に関する収縮および膨張は、ビーム要素ＥＢの伸縮変位で表現する。なお、血管内腔領域ＲＩには、強制変位履歴が割り当てられなくて良い。
【０１３３】
図２１は、強制変位履歴の割り当ての他の例を示す図であり、血管周辺組織のダミー要素ＲＤを含む形状モデルの断面を示している。図２１に示すように、ダミー要素ＲＤは、血管壁領域ＲＮの外側に設定される。形状モデルがダミー要素ＲＤを含む場合、血管壁領域ＲＮに加え、ダミー要素ＲＤにも節点ＰＮ４が設定される。節点ＰＮ４にも強制変位履歴が割り当てられる。
【０１３４】
力学モデル構築部５５は、境界条件および材料モデルを同定する場合、血管壁領域ＲＶに含まれる節点ＰＮ３に強制変位履歴を割り当て、境界条件および材料モデルを同定しない場合、強制変位履歴を割り当てなくても良い。節点ＰＮ３に強制変位履歴を割り当てる場合、内腔領域ＲＩに関する形状指標以外にも血管壁領域ＲＶに関する形状指標も参照して材料モデルの同定が行われる。
【０１３５】
図２２は、強制変位履歴の割り当ての他の例を示す図であり、プラーク領域ＲＰを含む形状モデルの断面を示している。図２２に示すように、プラーク領域ＲＰは、血管壁領域ＲＶに含まれる。プラーク領域ＲＰは、材料モデル同定領域に設定される。プラーク領域ＲＰについては、内腔形状指標、血管壁形状指標、およびプラーク指標を考慮して材料モデルが同定される。
【０１３６】
既述のように、プラーク指標は、例えば、超音波診断装置による組織性状診断により得られたプラークの性状に関するデータである。力学モデル構築部５５は、性状に応じてプラーク領域を複数の部分領域に区分し、複数の部分領域に個別に材料モデル同定領域を設定する。各部分領域には、予め、当該部分領域の性状に応じたパラメータ範囲を設定することが好ましい。既述の統計的同定処理により、各部分領域についての材料モデルパラメータが統計的同定部６１により同定される。そして、ステップＳ１０において表示機器３１が力学的指標として血管の材料特性に関する指標を表示することにより、ユーザは、プラークの性状を正確且つ容易に把握することができる。
【０１３７】
次に、潜在変数の一つである材料モデルについて詳細に説明する。血管の材料モデルとしては、弾性モデル、超弾性モデル、異方性超弾性モデル、粘性特性を考慮した超弾性モデル等が適用可能である。異方性超弾性モデルとしては、例えば、Ｙ．Ｃ．Ｆｕｎにより提案された数理モデルや、Ｈｏｌｚａｐｆｅｌ−Ｇａｓｓｅｒ構成式と呼ばれる数理モデルを適用できる。
【０１３８】
単位参照体積あたりのひずみエネルギーは、以下の（１）式で表わされる。（１）式の第１項は、コラーゲンを含まない等方性基礎材料のせん断変形に関するエネルギーを表している。また、（１）式の第２項は、コラーゲンを含まない等方性基礎材料の体積変形に関するエネルギーを表している。また、（１）式の第３項は、コラーゲン繊維各グループの寄与（繊維方向の分散を考慮）を表している。
【０１３９】
【数１】
2015097759.tif 000003
【０１４０】
図２３は、繊維グループの変形を示す図である。図２３に示すように、円筒形上の外膜を仮定する。芯線方向ｚと周方向θとにより規定される面における平均方向Ａでの繊維グループの変形は、下記の（２）式で表される。
【０１４１】
【数２】
2015097759.tif 000004
【０１４２】
（１）式および（２）式中の材料モデルに関するパラメータとしては、以下の表１に示すように、材料パラメータと繊維分散パラメータとがある。材料パラメータとしては、Ｃ１０、Ｄ，Ｋ１、Ｋ２等が用いられ、繊維分散パラメータとしては、Ｋａｐｐａやγ等が用いられる。各パラメータのデフォルト値と制約条件とは表１に示す通りである。
【０１４３】
【表１】
2015097759.tif 000005
【０１４４】
血液の材料モデルは、以下の（３）式のようなＣａｓｓｏｎ構成式や（４）式のようなＨＢ構成式が好適である。
【０１４５】
【数３】
2015097759.tif 000006
【０１４６】
【数４】
2015097759.tif 000007
【０１４７】
これら材料モデルのパラメータは、上述のステップＳ６において統計的同定部６１により、血管形態指標および血液流体指標を用いた統計的同定処理により同定される。
【０１４８】
次に、統計的同定部６１により行われる統計的同定処理の詳細について説明する。
【０１４９】
時系列の医用画像から計測される血管形態指標および血液流量指標のような観測変数は、不確定性を有している。統計的同定部６１は、このような不確定性が存在する状況下における潜在変数の統計的同定法として、階層ベイズモデルとマルコフ連鎖モンテカルロ法とに基づく統計的手法を活用する。
【０１５０】
既述のように、ステップＳ６において統計的同定部６１は、ステップＳ４において算出された血管形態指標または血液流量指標の予測値と観測値とに基づくデータ分布を設定する。データ分布は、例えば血管形態指標または血液流量指標の予測値と観測値との誤差に関する多変量正規分布関数を示す。データ分布は、時刻毎に個別に設定されても良いし、複数時刻まとめて設定されても良い。
【０１５１】
次に、統計的同定部６１は、形状モデルの強制変位と潜在変数とに事前分布を割り当てる。事前分布は、取り得る値の確率分布を示す。次に、統計的同定部６１は、潜在変数に関する数値シミュレーションのパラメータサーベイを実行し、潜在変数と血管形態指標または血液流量指標との関係を表現するモデルを構築する。例えば、材料モデルパラメータと内圧分布パラメータと血管形態指標または血液流量指標との関係がモデルに規定される。なお、血管形態指標または血液流量指標と潜在変数との関係は、モデルという形態ではなく、データベースまたはテーブルにより規定されても良い。これらモデル、データベース、または、テーブルは、記憶装置３３に記憶される。
【０１５２】
統計的同定部６１は、モデル、データベース、または、テーブルを利用して事前分布から血管形態指標または血液流量指標の確率分布を算出する。統計的同定部６１は、階層ベイズモデルとマルコフ連鎖モンテカルロ法とにより得られる事後分布から潜在変数を統計的に同定する。
【０１５３】
具体的には、この同定問題は、次の３条件を満足しない不良設定問題となる。３条件は、（１）解の存在が保証される、（２）解が唯一に定まる、（３）解がデータに対して連続的に変化し測定誤差に対して解が安定している、である。不良設定問題は、正規化理論とその拡張とからなる枠組み内で捉えると扱いやすい。不良設定問題を解くことは、標準正規化理論では不十分である。不良設定問題を解くためには、内部状態の不連続点を検出し、検出された不連続点を内部状態の推定に役立てる方法論が必要となる。この点で、本実施形態の統計的同定処理においてマルコフ確率場理論が有効となる。
【０１５４】
血管形態指標と血液流量指標とに不確定性が存在する環境下において統計的同定部６１は、適切な制約条件の下において潜在変数の確率分布パラメータを同定する。適切な制約条件を決定するためには、解の性質を事前に知っている必要がある。統計的同定部６１は、解空間の制約条件に関するデータベースをシミュレーションと観測値とに基づいて、発生する。
【０１５５】
統計的同定部６１は、発生されたデータベースを利用して、超多自由度大規模問題に対してマルコフ確率場理論と階層ベイズモデルとに基づく統計的同定処理を実行する。制約条件となる事前分布の設定では、多くの数値実験結果に基づいて、これらの要因に関するパラメータの確率分布が並列的に個別に構成される。統計的同定部６１は、複数の確率分布を統合しデータの欠損を補間することで潜在変数のパラメータを同定する。この処理のために、統計的同定部６１は、マルコフ確率場理論を用いたモデルに基づく階層ベイズ法による推定を行う。解析対象とする構造の変形状態の実測結果から、同定した中間変数をもとに、ある負荷条件と境界条件における圧力や流量分布が推定できるという仕組みである。
【０１５６】
冠動脈の構造流体解析における材料モデル、境界条件、および負荷条件の同定問題は、非線形逆解析と位置付けられ、解の一意性および安定性が保証されない場合が多い。生体組織の材料特性や血圧の現実的に取り得る範囲は先験情報として想定できるため、これらを事前分布の確率分布として設定できる。また、圧力や変位は空間的時間的に滑らかであることも想定できるためこの情報も先験情報として事前分布の確率分布として設定できる。
【０１５７】
または、血液の流れに逆流が生じていない事実を考慮できる場合は血管芯線方向の全体的な圧力分布の傾きは負（圧力降下が存在）であることも制約条件として用いて良い。負荷条件（内圧分布等）、境界条件、および材料モデルに対して、時系列のＣＴ画像に基づく血管形状変形指標の観測値と力学モデルに基づく血管形状変形指標の予測値との２乗誤差分布をデータ分布として設定できる。
【０１５８】
観測可能な平均流量に関する２乗誤差分布もデータ分布として追加してもよい。これらの事前分布とデータ分布とに基づいて、階層ベイズモデルとモンテカルロ法とを利用して事後分布を算出することができる。事後分布の発生確率や分散により、潜在変数のパラメータの同定値を得ることができる。発生確率が高く、分散が小さいほうが確信度合の高い同定値といえる。
【０１５９】
事後分布が多峰性分布となる場合でも、複数の同定値のうちの分散が小さい同定値を選択すれば良い。または、複数の同定値が存在し得る場合、それぞれの同定条件で、構造流体解析を実施し、それぞれの可能性を認識して、同定値や解析結果を、診断や予防の指針情報として活用することができる。時系列のＣＴ画像にも誤差が含まれていることから、力学モデルの各節点に関する血管形態指標にも誤差が含まれる。このため、各血管形態指標を、例えば、時系列のＣＴ画像から計測された血管形態指標の予測値を平均値とした正規分布の確率変数として扱い、位置の空間的順序を保つという制約を含めた上で、事前分布を設定してもよい。
【０１６０】
また、潜在変数のパラメータの同定において、一意性がなく、複数の候補が考えられる場合がある。この場合、時系列のＣＴ画像から計測された血管形態指標の不確定性に従う乱数のサンプリング点に対する、潜在変数の同定値のサンプル集合の変動幅をチェックすることで、各同定値の候補のロバスト性（安定性）を判定する。各同定値の候補のロバスト性に基づいて、最終的な同定値を決定しても良い。
【０１６１】
次に、力学モデルの詳細について説明する。既述のように、力学モデル構築部５５は、力学モデルの種類に応じて異なるタイプの力学モデルを構築することができる。連続体力学に基づくＦＥＭを用いる場合、力学モデル構築部５５は、血管壁の応力解析用のための形状モデル（ＦＥＭモデル）と血液の流体解析用のための形状モデル（ＦＥＭモデル）との両方を構築する。
【０１６２】
材料力学に基づく簡易的な同定法を用いる場合、材料力学における内圧を受ける厚肉円筒の式から近似的に圧力と弾性率と変位との関係を求める。この場合、芯線方向に配列された複数の離散化領域の各々について厚肉円筒近似が形状モデルとして用いられる。具体的には、力学モデル構築部５５は、芯線上に離散的に配列された節点を通る断面上の血管内腔形状と血管壁表面形状と断面中心とを特定する。
【０１６３】
次に力学モデル構築部５５は、血管内腔形状と血管壁表面形状とに基づいて、平均面積、内腔の平均半径、および平均壁厚を算出する。そして力学モデル構築部５５は、平均面積、内腔の平均半径、および平均壁厚に基づいて、各離散化領域の血管領域に厚肉円筒近似を施して形状モデルを構築する。
【０１６４】
流れ力学に基づく簡易的な同定法を用いる場合、血流の平均圧力と平均流量とを近似的に求めるため、流体力学における修正ベルヌーイの式、または、ハーゲン・ポアズイユ流れ（Ｈａｇｅｎ−Ｐｏｉｓｅｕｉｌｌｅ  Ｆｌｏｗ）の式を用いる。この場合、複数の離散化領域の各々の血圧差と流量との関係を近似的に求めるための形状モデルが構築される。ここで、血圧差は、入口の血圧と出口の圧力との圧力差であり、流量は、単位時間あたり入口流量（または流速）と出口流量（または流速）とを意味する。ただし、力学モデル構築部５５は、ステップＳ２において算出された芯線方向に関する伸縮距離とねじれ角とに基づいて、各節点に関する移動ベクトルと当該節点に隣接する節点に関する回転変位とを各時刻で対応させると良い。
【０１６５】
次に、連続体力学モデルを用いる構造流体解析について説明する。
【０１６６】
血管および血液（血管と血液を構成する物質）についての運動学は、その運動を引き起こしている力とは無関係である。血管および血液についての運動学の基本概念は、位置、時間、物体、運動、および変形し得る物質の集合についての直感的概念を数学的用語へ抽象化したものである。血管および血液に関する変形および運動の局所的解析を支配する基礎的運動学テンソルは、変形勾配テンソルＦと速度勾配テンソルＬとである。
【０１６７】
変形勾配テンソルＦは、運動する血管および血液の物質要素に生じる大きさおよび形の変化を規定する。変形勾配テンソルＦは、回転テンソル（正格直交テンソル）Ｒとストレッチングテンソル（正値対称テンソル）Ｕ，Ｖとの積で表わされる。ストレッチングテンソルＵ，Ｖは、まず、基本形態での正規直交ベクトルＲによって定められる方向に関するストレッチを課し、次に、正規直交ベクトルＲによって与えられる剛体回転を課すことによりもたらされる。なお、ストレッチを課す順番と剛体回転を課す順番とは逆でも良い。
【０１６８】
速度勾配テンソルＬは、基準形態に依存せず現在形態にのみ依存する。速度勾配テンソルＬは、運動する血管および血液の物質要素に生じる大きさおよび形の変化が生じる速度を規定する。速度勾配テンソルＬは、ひずみ速度テンソルＤ（対称テンソル）とスピンテンソル（反対称テンソル）とに分離できる。ひずみ速度テンソルＤは、物体がその現在形態をちょうど通過するときのストレッチの変化率を表す。スピンテンソルは、物体がその現在形態をちょうど通過するときの回転の変化率を表す。
【０１６９】
血管の力学モデルの外表面節点（または積分点）の一部に関する変形勾配テンソルと速度勾配テンソルとに変位拘束条件（時間的変化を含む）を割り当てる。力学モデルの内腔節点（または積分点）に関する予測値（変形勾配テンソル、速度勾配テンソル、またはそれらの関数値（例えば、変位や面積でも良い））と、観測値（観測データから得られた変形勾配テンソル、速度勾配テンソル、またはそれらの関数値）とが整合するように、材料モデルのパラメータ、負荷条件（力学モデルの内腔における表面力ベクトル）、および境界条件（血管境界における力ベクトル）を同定する。ここで、血管内部の応力の初期状態は、予め仮定しておいてもよいし、同定してもよい。
【０１７０】
連続体力学に基づく力学モデルは、運動中の血管および血液における質量、運動学、運動量、角運動量、およびエネルギーの平衡を表す方程式を基礎としている。質量、力、熱、および内部エネルギーという概念が基本である。平衡則とは、運動量の時間全微分が物体力と接触力との和に等しくなり、また、角運動量の時間全微分が物体トルクと接触トルクとの和に等しくなり、また、運動エネルギーと内部エネルギーとの時間的変化が、仕事率（力学的エネルギー）と単位時間あたりの熱供給と熱流束との和に等しくなる、ということを意味している。
【０１７１】
平衡則、構成式、および跳躍条件から、変形勾配テンソル、速度勾配テンソル、および応力テンソルといった場の方程式を導き、血管および血液の力学モデルを記述できる。また、血管および血液におけるひずみ場は、適合条件を満たす。
【０１７２】
ここで、構成式とは、密度、内部エネルギー、速度ベクトル、応力テンソル、熱流束ベクトル、および温度からなる１０個のスカラ方程式の組の関係を与えるものである。場の釣り合い式の中の１７個のスカラ場、つまり、密度、内部エネルギー、速度ベクトル、応力テンソル、および熱流束ベクトルのうち、場の釣り合いにより８個のスカラ関係が与えられ、温度を加えた残りの未知量の関係を与える。ただし、物体力ｂと熱源ｒとは既知としている。これらのスカラ場を与える方程式のパラメータが材料モデルパラメータである。
【０１７３】
連続体力学に基づく力学モデルは、有限要素法又は境界要素法による数値解析法により変位ベクトル、応力テンソル、ひずみテンソル、および速度ベクトル等の場の近似解を、与えられた境界条件、負荷条件、および材料モデルのもとに算出することができる。
【０１７４】
構造−流体連成解析において、構造および流体の方程式を解く方法は、一体型解法（monolithic method）と分離型解法（partitioned method）とのいずれでも良い。また、構造と流体との間の境界面での連成は、弱連成でも強連成でもよい。また、流体解析においては、ＡＬＥ法に代表されるような境界面追跡型の手法で移動境界を扱ってもよいし、「Ｉｍｍｅｒｓｅｄ  Ｂｏｕｎｄａｒｙ  ｍｅｔｈｏｄ」、「Ｉｍｍｅｒｓｅｄ  Ｆｉｎｉｔｅ  Ｅｌｅｍｅｎｔ  Ｍｅｔｈｏｄ」、または「Ｆｉｃｔｉｔｉｏｕｓ  Ｄｏｍａｉｎ  Ｍｅｔｈｏｄ」等の境界面補足型の手法でもよい。
【０１７５】
次に、簡易的な力学モデルの例として、内圧と外圧とを受ける厚肉円筒の材料力学の式と、ハーゲン・ポアズイユ流れおよび修正ベルヌーイの式とについて詳細に説明する。
【０１７６】
まず、図２４と図２５とを参照しながら、厚肉円筒の材料力学の式について説明する。図２４は、肉厚円筒の力学モデルの直交断面を示す図である。図２５は、図２４の微小扇形要素の拡大図である。内半径ｒａ、外半径ｒｂの厚肉円筒に内圧ｐａと外圧ｐｂとが作用する場合の応力やひずみ、変位等の式を説明する。Ｅおよびνは材料モデルパラメータである。Ｅは弾性率、νはポアソン比を表している。厚肉円筒では，半径応力σｒも考慮し、円周応力σθ  の半径方向分布も考慮する必要がある。以降では、軸方向のひずみεｚは、断面の位置および向きに関して一様とする。
【０１７７】
また円筒断面は、軸対称である場合について説明するが任意形状であっても良い。円筒断面が軸対称である場合、平衡条件は、任意の断面の半径方向についてのみ考えれば良い。任意の断面上で半径ｒとｒ＋ｄｒの同心円筒と中心角ｄθで切り取られた単位厚さ１の微小扇形要素について、半径方向に関する力の平衡を考える。変形も軸対称であるため、ａｂ面およびｂｃ面には、せん断応力が生じないので、垂直応力のみ作用する。従って、半径方向に関する力の平衡は、以下の（５）式のように表現することができる。
【０１７８】
σr rdθ+2σθr sin(dθ/2)-（σr+(dσr/dr)dr）(r+dr)dθ=0    …（５）
【０１７９】
ここで、ｄｒはｒより小さく、ｄσｒ  はσｒより小さいので、（５）式に含まれる高次の微小項を省略し、sin(dθ/2)≒dθ/2とすれば、（５）式は、以下の（６）式のように表現することができる。
【０１８０】
rdσr/dr+σr-σθ=0    …（６）
【０１８１】
半径ｒにおける半径方向の変位をｕとすれば、ｕ＋ｄｒでの同方向の変位はｕ＋（ｄｕ／ｄｒ）ｄｒとなるので、半径ｒにおける半径方向のひずみεｒはεｒ＝ｄｕ／ｄｒとなる。また、半径方向の変位ｕによって、半径ｒの円は半径ｒ＋ｕの円になる。従って、円周ひずみεθは、以下の（７）式のように表現することができる。
【０１８２】
εθ =(2π(r+u)-2πr)/2πr=u/r    …（７）
【０１８３】
また、応力とひずみの関係式から、以下の（８）式または（９）式が得られる。
【０１８４】
d2u/dr2+(1/r)(du/dr)-u/r2=0    …（８）
【０１８５】
d2u/dr2+d(u/r)/dr=0    …（９）
【０１８６】
（８）式または（９）式を積分すると以下の（１０）式を得ることができる。
【０１８７】
u=c1r+c2/r    …（１０）
【０１８８】
これにより以下の（１１）式、（１２）式、（１３）式が得られる。
【０１８９】
σr=(E/((1+ν)(1-2ν)))(c1-(1-2ν)(c2/r2)+νεz)     …（１１）
【０１９０】
σθ=(E/((1+ν)(1-2ν)))(c1+(1-2ν)(c2/r2)+νεz)     …（１２）
【０１９１】
σz=(Eν/((1+ν)(1-2ν)))(2c1+((1-ν)/ν)εz)     …（１３）
【０１９２】
（１１）式、（１２）式、（１３）式における定数ｃ１，ｃ２は、周辺条件、すなわち、円筒の内周ｒ＝ｒａでσｒ＝−ｐａ、外周ｒ＝ｒｂでσｒ＝−ｐｂから定めることができる。この周辺条件により、（１１）式、（１２）式、および（１３）式から以下の（１４）式、（１５）式、（１６）式をそれぞれ得ることができる。また、変位ｕは、以下の（１７）式のように表現することができる。
【０１９３】
σr=(1/(rb2-ra2))(ra2(1-rb2/r2)pa-rb2(1-ra2/r2)pb)     …（１４）
【０１９４】
σθ=(1/(rb2-ra2))(ra2(1+rb2/r2)pa-rb2(1+ra2/r2)pb)     …（１５）
【０１９５】
σz=2ν(ra2pa-rb2pb)/(rb2-ra2)+Eεz=ν(σr+σθ)+Eεz   …（１６）
【０１９６】
u=((1+ν)(1-2ν)/E)((ra2pa-rb2pb)/(rb2-ra2))r+((1+ν)/E)((ra2rb2)/((rb2-ra2)r))(pa-pb)-νεzr    …（１７）
【０１９７】
（１６）式および（１７）式は、εｚの項が含んでいる。従って、（１６）式は、対象とする円筒の境界における拘束条件に応じて異なる。例えば、対象とする円筒の両端が拘束される場合、εｚ  ＝０となり、両端が開放され、σｚ＝０となる。
【０１９８】
ｒθ面内のせん断応力τｒは、以下の（１８）式のように表現することができる。
【０１９９】
τr=(1/2)(σr-σθ)=((ra2rb2)/((rb2-ra2)r2))pb   …（１８）
【０２００】
θｚ面内のせん断応力τｒ’は、σｚ＝０のとき最大になるのでσｚ＝０とすると、以下の（１９）式のように表現することができる。
【０２０１】
τr’=(1/2)|σθ|=(1/2)((rb2)/(rb2-ra2))(1+ra2/r2)pb   …（１９）
【０２０２】
次に、ハーゲン・ポアズイユ流れおよび修正ベルヌーイの式について説明する。流体が円管内へ流入すると、下流に進むにつれて圧力は降下するとともに、流れの速度分布も徐々に変化する。流れが管内へ流入すると管壁から境界層が発達し、下流へ進むにつれて境界層は厚さを増し、ついには管内の流れは境界層に覆われる。このため、速度分布は管入口のほぼ平らな分布から下流の放物形分布へと変化し、それ以降の速度分布は変化しなくなる。この状態を完全に発達した流れと呼び、管摩擦損失による圧力降下も一定の割合になる。流れが管入口から発達した流れに達する区間を助走区間または入口区間といい、その区間の長さを助走距離または入口長さと呼ぶ。完全に発達した流れの速度分布は下流方向へ変化しないから、以下の（２０）式に示すように、管摩擦損失によって生じる圧力損失ΔＰの作用力と流体の粘性によって生じるせん断応力τの摩擦力は釣り合うことになる。円管内流れの場合、レイノルズ数Ｒｅがおよそ２３００以下のとき層流となるといわれている。
【０２０３】
【数５】
2015097759.tif 000008
【０２０４】
速度分布ｕは層流の場合、以下の（２１）式のように軸対称な回転放物面で表わされる。
【０２０５】
【数６】
2015097759.tif 000009
【０２０６】
流量Ｑは、速度分布ｕを管断面全体にわたって積分することにより、以下の（２２）式のように表現することができる。
【０２０７】
【数７】
2015097759.tif 000010
【０２０８】
断面平均流速ｖは、以下の（２３）式のように表現することができる。
【０２０９】
【数８】
2015097759.tif 000011
【０２１０】
圧力勾配は、流れ方向へ一定で、圧力は減少する。圧力勾配は、管長ｌの圧力降下Δｐを用いると、以下の（２４）式のように表現できる。（２４）式に（２２）式を代入すると（２５）式が得られる。（２５）式は、流量Ｑが圧力損失Δｐに比例することを意味する。この関係を満たす流れをハーゲン・ポアズイユ流れという。
【数９】
2015097759.tif 000012
【０２１１】
【数１０】
2015097759.tif 000013
【０２１２】
損失がある場合の修正ベルヌーイの式は、以下の（２６）式、（２７）式のように表現することができる。
【０２１３】
【数１１】
2015097759.tif 000014
【０２１４】
【数１２】
2015097759.tif 000015
【０２１５】
ここで、ｌは管の長さ、ｄは管内径、ｖは管内平均速度、λは管摩擦係数である。管摩擦係数λは、流れが層流の場合にはレイノルズ数Ｒｅによって、乱流の場合にはレイノルズ数Ｒｅと表面粗さとによって定まる値になる。流体の粘性によって摩擦抵抗が必ず作用する。この摩擦抵抗は、流れを駆動する動力またはエネルギーを消費することになるのでエネルギー損失になる。
【０２１６】
ハーゲンポワズイユ流れの式をダルシー・ワイズバッハの式のように変形すると、圧力差と流速および管内径とに関する以下の（２８）式が得られる。
【０２１７】
【数１３】
2015097759.tif 000016
【０２１８】
以上でハーゲン・ポアズイユ流れおよび修正ベルヌーイの式の説明を終了する。
【０２１９】
上述のように、画像処理装置２７は、この材料力学の式とハーゲン・ポアズイユ流れおよび修正ベルヌーイの式とを利用して潜在変数を同定することができる。例えば、力学モデルとして、血管の変形を厚肉管の材料力学の式を用い、管径変化を内圧変化と弾性率とにより表現する場合について考える。
【０２２０】
無応力状態を初期形状（例えば、血管が最も収縮する状態）と仮定した場合、血管壁およびプラークの弾性率をある値に設定すると、血管内腔の平均半径等の血管形状変形指標の観測値の時間的変化量と内圧の変化量との関係式が得られる。血管形状変形指標の観測値は、時系列のＣＴ画像から計測される。この血管形状変形指標の観測値の時間的変化量に合致するように血管の内圧分布の時間的変化が決定される。この内圧分布の下に血液の流体解析を行うことで血管流量指標の予測値が計測される。この血管流量指標の予測値が観測値に一致しない場合、画像処理装置２７は、最初に決めた血管壁またはプラークの弾性率を変更して、さらに同様の解析を行う。
【０２２１】
これを繰り返すことにより、画像処理装置２７は、血管形状変形指標の観測値と血液流量指標の観測値とに整合する血管壁およびプラークの弾性率、内圧分布、流体解析の圧力境界条件等の潜在変数を決定することができる。この決定方法をより効率的かつ安定的に行うために、階層ベイズモデルとマルコフ連鎖モンテカルロ法とによる統計的同定手法を用いてもよい。
【０２２２】
上記のように、血管解析装置５０は、記憶装置３３、領域設定部５１、画像解析・追尾処理部５３、力学モデル構築部５５、および統計的同定部６１を有している。記憶装置３３は、被検体の血管に関する時系列の医用画像のデータを記憶する。領域設定部５１は、時系列の医用画像に含まれる血管領域に解析対象領域を設定する。画像解析・追尾処理部５３は、時系列の医用画像を画像処理して解析対象領域の時系列の形態指標と時系列の形状変形指標とを算出する。
【０２２３】
力学モデル構築部５５は、時系列の形態指標と時系列の形状変形指標と時系列の医用画像とに基づいて、解析対象領域の構造流体解析に関する力学モデルを暫定的に構築する。統計的同定部６１は、暫定的に構築された力学モデルに基づく血管形態指標および血液流量指標の予測値が予め計測された血管形態指標および血液流量指標の観測値に整合するように解析対象領域に関する力学モデルの潜在変数を同定する。
【０２２４】
上記構成により、実施の形態の血管解析装置５０は、材料モデル、境界条件、負荷条件、および幾何学的構造等の潜在変数を血管形状変形指標と血液流量指標とを用いた逆解析により同定することができる。
【０２２５】
血管解析装置５０は、潜在変数を変更しながら逆解析を反復して行うことにより、上述の４点の困難、すなわち、１．冠動脈の材料モデルの同定方法、２．心臓の形状の変形の冠動脈への影響の組み込み、３．冠動脈の境界条件の同定方法、４．不確定性を有する血管形状を利用した材料モデルや負荷条件、境界条件の同定方法、を全て加味した潜在変数を同定することができる。
【０２２６】
従って、血管解析装置５０は、ＣＴ画像に描出されない血管や心臓等の外部要因による影響を加味した構造流体解析を実行することができる。このように、実施の形態の医用画像診断装置によれば、血管の構造流体解析の精度の向上を図ることができる。
【０２２７】
次に、このような血管の力学モデルによる構造解析、流体解析、または構造流体連成解析結果から得られた様々な指標（血管の断面積、応力の変化、血流の流速・流量・圧力の変化等）を用いてノミナルモデルを作成し、診断対象のデータとのモデル同定により血管の異常を診断する方法について説明する。なお、以下、血管の異常診断として、血管の狭窄度の診断を行う例を説明する。
【０２２８】
ある２組の指標を選択し、一方を入力、他方を出力としたとき、入出力関係に因果関係があるならば入出力モデルを仮定し、モデルを同定することが可能となる。すなわち、図２６に示すように、狭窄の有無に依存しない指標を入力ｕ１とし、狭窄の有無で影響を受ける指標を出力とする。また、狭窄なしの出力をｙ１、狭窄ありの出力をｙ２とすると、狭窄なしの場合は入力ｕ１、出力ｙ１としたモデルＧ１が仮定される。同様に、狭窄ありの場合は、入力ｕ１、出力ｙ２としたモデルＧ２が仮定される。
【０２２９】
このとき、モデルＧ１およびモデルＧ２が有限個のパラメータで表現され、モデル同士を比較できるならば、血管の狭窄度を評価できる。モデルの同定精度が高ければ、推定したモデルＧ１ｈａｔに、入力ｕ１を入力した時に出力される出力ｙ１ｈａｔは、観測したｙ１との誤差が小さい。なお、「ｈａｔ」は、推定値を意味する。
【０２３０】
一般的に推定精度を上げるためには、モデルの次数を上げる必要がある。高次のモデル同士の比較は、評価対象となるパラメータが多くなるため、結果が不安定なものになりやすい。そこで、予め同定したモデルを「ノミナルモデル」とし、その出力ｙ１ｈａｔを計算しておく。ｙ１ｈａｔを入力とし、比較対象となる指標ｙ２を出力とする新たな入出力モデルＧ１２を仮定し、同定する。このモデルは狭窄の有無を表現したものであるため、以下、「狭窄モデル」と呼ぶこととする。ｙ１ｈａｔには、高次のモデルの要素が含まれている。このため、狭窄モデルとして、次数を減らした低次モデルを用いることができる。
【０２３１】
一例を図２７に示す。図２７のグラフは、横軸が時間、縦軸がＣＴ値（造影剤の濃度）を示している。この図２７のグラフにおいて、○で示すグラフは、狭窄ありの出力ｙ２のグラフである。実線のグラフは、ノミナルモデルＧ2を同定し、共通の入力ｕ1を入力したときの出力ｙ２ｈａｔのグラフである。これに対し、点線のグラフは、狭窄なしの出力y1から同定した別のノミナルモデルＧ１で推定した出力ｙ１ｈａｔを入力とし、ｙ２を出力とした狭窄モデルＧ１２を同定し、ｙ１ｈａｔを入力した時の推定値のグラフである。ノミナルモデルＧ１は３／６次、狭窄モデルＧ１２は２／２次であり、少ないパラメータでありながら推定精度が高いのが分かる。
【０２３２】
図２８に、実施の形態の医用画像診断装置の、血管の狭窄度の診断を行う場合における、血管解析装置５０の機能ブロック図を示す。この場合、血管解析装置５０は、構造流体解析を行う上述の画像処理装置２７、およびノミナルモデル（後述するＡＲＸモデルの各種パラメータ）が記憶される記憶装置３３と共に、造影剤注入プロファイル取得部７１を有する。また、血管解析装置５０は、第１入出力モデル算出部７２、時系列ＣＴ画像取得部７３、血管芯線抽出部７４、時系列ＣＴ値抽出部７５、第２入出力モデル算出部７６、狭窄モデル算出部７７、および狭窄度診断部７８を有している。
【０２３３】
造影剤注入プロファイル取得部７１は、注入プロファイル取得部の一例である。時系列ＣＴ画像取得部７３は、時系列ＣＴ値取得部の一例である。画像処理装置２７は、構造流体解析部の一例である。第１入出力モデル算出部７２は、第１算出部の一例である。造影剤注入プロファイル取得部７１、第１入出力モデル算出部７２、時系列ＣＴ画像取得部７３、画像処理装置２７、および記憶装置３３は、時系列モデル生成部の一例である。記憶装置３３は、記憶装置の一例である。時系列ＣＴ値抽出部７５は、時系列ＣＴ値抽出部の一例である。狭窄モデル算出部７７は、狭窄モデル算出部の一例である。狭窄度診断部７８は、診断部の一例である。第２入出力モデル算出部７６は、第２算出部の一例である。
【０２３４】
造影剤注入プロファイル取得部７１〜狭窄度診断部７８は、全てをハードウェアで構成してもよいし、全てをソフトウェアで構成してもよい。または、一部をハードウェアで構成し、残りをソフトウェアで構成してもよい。
【０２３５】
まず、血管の狭窄度の診断を行う場合、上述のノミナルモデルを事前に用意する。順を追って説明すると、ＣＴ画像の撮像には、血管造影のために造影剤を用いることが多い。造影剤が、例えば上腕部の静脈より一定量注入されると、造影剤が通過した領域のＣＴ値（ＨＵ値）が上昇し、血管以外の組織との識別が容易になる。ＣＴ架台１０は、時系列でＣＴ画像を撮影するように、システム制御部２１および架台制御部２３により制御される。これにより、ＣＴ値が変化する速度を観測することができる。
【０２３６】
一例として、図２９に造影剤の注入パターンを示す。図２９の実線のグラフ、点線のグラフ、および一点鎖線のグラフは、それぞれ造影剤の注入パターンを示している。図２９の実線のグラフは、造影剤を約７秒間かけてゆっくり注入した注入パターンを示している。図２９の点線のグラフは、造影剤を約４．５秒間かけて高速に注入した注入パターンを示している。図２９の一点鎖線のグラフは、造影剤を約６秒間かけて中速で注入した注入パターンを示している。なお、造影剤の注入は、担当医師等が手動で行ってもよいし、担当医師監督のもと、注入機を用いて行ってもよい。
【０２３７】
図３０は、造影剤注入後の注目領域におけるＣＴ値の変化のグラフである。図２９の実線のグラフ、点線のグラフ、および一点鎖線のグラフに示した注入パターンのうち、どの注入パターンで造影剤を注入しても、大体、図３０と同様のグラフを得られる。この図３０において、○のグラフは、狭窄なしの場合のＣＴ値の変化を示している。また、図３０において、●のグラフは、狭窄ありの場合のＣＴ値の変化を示している。この図３０は、注目領域に狭窄が無い場合、ＣＴ値の変化は、○のグラフに示すように立ち上がりが早く、また、ピークも高くなる。これに対して、注目領域の上流側に狭窄部が存在する場合、ＣＴ値の変化は、●のグラフに示すように立ち上がりが遅く、また、ピークも低くなる。実施の形態の医用画像診断装置では、造影剤の注入パターンを入力、注目領域でのＣＴ値を出力とした入出力モデルを用いている。
【０２３８】
図２８の説明に戻り、ノミナルモデルを生成する場合、造影剤注入プロファイル取得部７１は、図２９を用いて説明した造影剤の注入プロファイルを取得する。取得された造影剤の注入プロファイルは、画像処理装置２７、第１入出力モデル算出部７２、および第２入出力モデル算出部７６に供給される。
【０２３９】
上述の画像処理装置２７は、時系列ＣＴ画像取得部７３で取得された時系列のＣＴ画像（ＭＲＩ画像等でもよい）を、構造流体解析する。なお、時系列ＣＴ画像取得部７３は、記憶装置３３に記憶されている時系列のＣＴ画像を取得してもよいし、ＣＴ架台１０でリアルタイムかつ時系列で撮像されたＣＴ画像を取得してもよい。画像処理装置２７の解析により、流体の一つである造影剤の注入プロファイルに対応する出力（造影剤の濃度変化）をシミュレーションすることが可能となる。構造流体解析部は、図２９を用いて説明した注入パターンで注入された造影剤の濃度変化をシミュレーションする。そして、このシミュレーション結果となる図３０に示す造影剤の濃度変化を示す情報を、第１入出力モデル算出部７２に供給する。
【０２４０】
第１入出力モデル算出部７２（および第２入出力モデル算出部７６）は、有限個のパラメータからなる線形差分方程式で表現されるモデルを算出する。すなわち、第１入出力モデル算出部７２（および第２入出力モデル算出部７６）は、有限個のパラメータで表現する入出力モデルの一例として、いわゆるＡＲＸモデルを算出する。ＡＲＸは、「Ａｕｔｏ  Ｒｅｇｒｅｓｓｉｖｅ  ｗｉｔｈ  ｅＸｏｇｅｎｏｕｓ  ｉｎｐｕｔｓ）、または（Ａｕｔｏ  Ｒｅｇｒｅｓｓｉｖｅ  ｗｉｔｈ  ｅＸｔｅｒｎａｌ  ｉｎｐｕｔｓ）の略記である。ＡＲＸモデルは、入力をｕ（ｔ）、出力をｙ（ｔ）、式誤差ε（ｔ）とすると、以下の（２９）式で表される。
【０２４１】
A(q)y(t) = B(q)u(t) + ε(t)   …（２９）
【０２４２】
ただし、Ａ（ｑ）は、以下の（３０）式、Ｂ（ｑ）は、以下の（３１）式のとおりである。
【０２４３】
A(q) = 1 + a1q-1 + … + amq-m  …（３０）
【０２４４】
B(q) = b1q-1 + … + bnq-n  …（３１）
【０２４５】
（２９）式における「ｔ」は、各時間でサンプリングされた離散的なサンプリングデータ（図３０の○または●のＣＴ値（サンプリングデータ）を参照）を意味している。また、（３０）式および（３１）式における「ｑ」が、時間軸に沿ってシフトするシフトパラメータを意味している。例えば、「ｑ-1」は、時間的に一つ前の値を用いることを意味している。また、この例のモデルは、合計ｍ＋ｎ個のパラメータで表現されている。このモデルのパラメータＡ、Ｂは、式誤差の二乗和を最小化するように最小二乗法を用いて決定される。
【０２４６】
このように算出されたＡＲＸモデルの各パラメータは、ノミナルモデルとして記憶装置３３に記憶される。一例として、記憶装置３３に記憶されるノミナルモデルとしては、狭窄なしの場合のノミナルモデル、狭窄度３０パーセントの場合のノミナルモデル、狭窄度５０パーセントの場合のノミナルモデル等の複数のノミナルモデルが記憶される。なお、この例では、入出力モデルとしてＡＲＸモデルを用いたが、パラメータで表現可能なモデルであれば、どのようなモデルを用いてもよい。
【０２４７】
次に、このようなノミナルモデルが記憶装置３３に記憶されると、血管の狭窄度診断が可能となる。血管の狭窄度診断時となると、造影剤が注入された被検体の時系列ＣＴ画像を、時系列ＣＴ画像取得部７３が取得する。また血管芯線抽出部７４は、取得された時系列ＣＴ画像から、上述のように血管の芯線を抽出する。時系列ＣＴ値抽出部７５は、芯線が抽出された血管の注目領域のＣＴ値を時系列で抽出する。この時系列で抽出された各ＣＴ値が、図２７のグラフの○印で示すＣＴ値に相当する。
【０２４８】
次に、第２入出力モデル算出部７６は、第１入出力モデル算出部７２と同様に、有限個のパラメータで表現する入出力モデルの一例として、いわゆるＡＲＸモデルを算出する。第２入出力モデル算出部７６は、各ＣＴ値を取得した際に被検体に注入した造影剤の注入プロファイルを、造影剤注入プロファイル取得部７１を介して取得する。そして、第２入出力モデル算出部７６は、造影剤の注入プロファイルおよび時系列で抽出された各ＣＴ値から、記憶装置３３に記憶されているノミナルモデルと同様のパラメータセットとなるＡＲＸモデルを、第２入出力モデルとして算出する。
【０２４９】
最後に、狭窄モデル算出部７７が、記憶装置３３に記憶されているノミナルモデル（図２７の実線のグラフに相当するパラメータ）と第２入出力モデルとの各パラメータの差分となる狭窄モデルを算出し、狭窄度診断部７８は、算出された狭窄モデルのパラメータを、血管の狭窄度として診断する。
【０２５０】
すなわち、造影剤注入プロファイルをノミナルモデルに入力した時の出力結果を入力とし、同じく造影剤注入プロファイルを第２入出力モデルに入力したときの出力結果を出力としてＡＲＸモデルで計算を行う。仮に、第２入出力モデルとノミナルモデルとの間に差が無ければ、得られた狭窄モデルの伝達関数は「１」となり、第２入出力モデルとノミナルモデルとが同一であることを示すようになる。
【０２５１】
狭窄度診断部７８は、狭窄モデルをそのまま評価（伝達関数１からどれくらい離れているか）してもよいし、例えば３０パーセントの狭窄、１０パーセントの狭窄等のように、血管の狭窄度を表示してもよい。狭窄度の異なるデータを使って作成した複数のノミナルモデルで評価を行い、「どのノミナルモデルに近いか」ということを表示してもよい。さらには、狭窄度診断部７８は、狭窄が無しの正常時の値を「０」として、スケールで狭窄の度合いを表示してもよい。すなわち、狭窄の度合いが分かれば、どのような表示形態で表示してもよい。
【０２５２】
以上の説明から明らかなように、実施の形態の医用画像診断装置は、画像追尾、構造流体解析および逆解析（統計的同定法）に基づく血管解析により、血行動態における血圧、血流、応力、またはひずみ等の狭窄等の血管病変の予防、診断に必要な指標を高精度に解析する。そして、実施の形態の医用画像診断装置は、狭窄の有無が判明している典型例を対象に解析を行い、低次の標準モデル（ノミナルモデル）を予め作成しておき、診断対象データとのモデル同定から狭窄度を診断する。
【０２５３】
これにより、画像処理装置で計算された解析結果を有効に利用して、血管の異常診断を可能とすることができる。また、低次の標準モデル（ノミナルモデル）と診断対象データとのモデル同定をから狭窄度を診断することで、少ない計算リソースおよび極めて短時間で血管の異常診断を行うことができる。
【０２５４】
また、時系列ＣＴ値抽出部７５で抽出された各ＣＴ値、および造影剤注入プロファイルから第２入出力モデル算出部７６でＡＲＸモデルを算出し、算出したＡＲＸモデルと、記憶装置３３に記憶されているノミナルモデルとを、狭窄モデル算出部７７で比較している。このため、同じＡＲＸモデル同士を比較することができるため、精度よく、狭窄度を診断することができる。
【０２５５】
（変形例１）
上述の実施の形態の説明では、同じＡＲＸモデル同士を比較するために、第２入出力モデル算出部７６でＡＲＸモデルを算出したが、時系列ＣＴ値抽出部７５で抽出された各ＣＴ値と、記憶装置３３に記憶されているノミナルモデルとを狭窄モデル算出部７７で比較して、狭窄度を診断してもよい。この場合、第２入出力モデル算出部７６を省略することができるため、実施の形態の医用画像診断装置の構成の簡略化および低コスト化を図ることができる。
【０２５６】
（変形例２）
時系列のＣＴ値として、血管の狭窄前後を含む芯線方向にＣＴ値を取得し、血管の走行方向における勾配（造影剤濃度勾配）を算出した勾配の時系列データを用いてもよい。図３１は、造影剤濃度勾配のグラフの一例である。図３１中、実線で示すグラフは、狭窄なしの場合の造影剤濃度勾配のグラフである。造影剤濃度勾配は、負の勾配となり、狭窄なしの場合、早い時間で立下り、早い時間で０レベルに復帰する勾配のグラフとなる。これに対して、図３１中、点線で示すグラフは、狭窄ありの場合の造影剤濃度勾配のグラフである。狭窄ありの場合、緩やかに時間をかけて立下り、また、緩やかに時間をかけて０レベルに復帰する勾配のグラフとなる。このような勾配の時系列データを用いても、上述と同様の効果を得ることができる。
【０２５７】
以上、実施の形態および変形例を説明したが、実施の形態および変形例は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施の形態および変形例は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。実施の形態およびその変形は、発明の範囲や要旨に含まれると共に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【０２５８】
１０  ＣＴ架台
１１  Ｘ線管
１３  Ｘ線検出器
１５  データ収集装置
２０  コンソール
２１  システム制御部
２３  架台制御部
２５  再構成装置
２７  画像処理装置
２９  入力機器
３１  表示機器
３３  記憶装置
５０  血管解析装置
５１  領域設定部
５３  画像解析・追尾処理部
５５  力学モデル構築部
５７  血管応力解析部
５９  血液流体解析部
６１  統計的同定部
７１  造影剤注入プロファイル取得部
７２  第１入出力モデル算出部
７３  時系列ＣＴ画像取得部
７４  血管芯線抽出部
７５  時系列ＣＴ値抽出部
７６  第２入出力モデル算出部
７７  狭窄モデル算出部
７８  狭窄度診断部
【図１】
2015097759.tif 000017
【図３】
2015097759.tif 000018
【図４】
2015097759.tif 000019
【図５】
2015097759.tif 000020
【図６】
2015097759.tif 000021
【図７】
2015097759.tif 000022
【図８】
2015097759.tif 000023
【図９】
2015097759.tif 000024
【図１０】
2015097759.tif 000025
【図１１】
2015097759.tif 000026
【図１２】
2015097759.tif 000027
【図１３】
2015097759.tif 000028
【図１４】
2015097759.tif 000029
【図１７】
2015097759.tif 000030
【図１８】
2015097759.tif 000031
【図１９】
2015097759.tif 000032
【図２０】
2015097759.tif 000033
【図２１】
2015097759.tif 000034
【図２２】
2015097759.tif 000035
【図２３】
2015097759.tif 000036
【図２４】
2015097759.tif 000037
【図２５】
2015097759.tif 000038
【図２６】
2015097759.tif 000039
【図２７】
2015097759.tif 000040
【図２８】
2015097759.tif 000041
【図２９】
2015097759.tif 000042
【図３０】
2015097759.tif 000043
【図３１】
2015097759.tif 000044
【図２】
2015097759.tif 000045
【図１５】
2015097759.tif 000046
【図１６】
2015097759.tif 000047