(57)【要約】【課題】血管の狭窄部位または閉塞部位を短時間で精度よく判別する。【解決手段】予め記憶されている被検体の血管の画像データを記憶部から読み出して血管解析モデルを生成し、生成された血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、当該所定位置から各毛細血管までの複数本のルート間の各圧力損失を求め、算出部で求めた複数本の圧力損失と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する。【選択図】図２
2015171486.tif 000002
【特許請求の範囲】
【請求項１】
血管内の血液の流れをシミュレーションによって解析する血流解析システムであって、
撮像によって予め取得した被検体の血管の画像データを記憶する記憶部と、
前記記憶部から画像データを読み出して血管解析モデルを生成する血管解析モデル生成部と、
前記血管解析モデルを表示する表示部と、
前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件を設定する操作部と、
前記設定された各条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、所定位置から各毛細血管までの複数本のルート間の各圧力損失を求める算出部と、
前記算出部で求めた複数本の圧力損失と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部と、
を備えたことを特徴とする血流解析システム。
【請求項２】
血管内の血液の流れをシミュレーションによって解析する血流解析システムであって、
撮像によって予め取得した被検体の血管の画像データを記憶する記憶部と、
前記記憶部から画像データを読み出して血管解析モデルを生成する血管解析モデル生成部と、
前記血管解析モデルを表示する表示部と、
前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件を設定する操作部と、
前記設定された各条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、起点から各毛細血管までの複数本のルート間を流れる血液の各流量を求める算出部と、
前記算出部で求めた複数本の流量と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部と、
を備えたことを特徴とする血流解析システム。
【請求項３】
請求項１または請求項２記載の血流解析システムにおいて、
前記判別部は、予め設定された大動脈圧と複数本のルートの遠位部の内圧の圧力損失、または、予め設定された冠動脈入口部の流量と複数本のルートの遠位部の流量の損失を求め、算出結果と予め決めた基準値を比較し、有意狭窄部位または閉塞部位を有する血管のルートの候補を判別する
ことを特徴とする血流解析システム。
【請求項４】
請求項１ないし請求項３のいずれかに記載の血流解析システムにおいて、
前記各ルートを所定間隔で複数個の血管ブロックに分割する分割処理部を備え、
前記判別部は、血管ブロック同士の継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別する
ことを特徴とする血流解析システム。
【請求項５】
請求項１ないし請求項４のいずれかに記載の血流解析システムにおいて、
前記血管解析モデルは、細線化された一次元要素から生成されている
ことを特徴とする血流解析システム。
【請求項６】
請求項１ないし請求項４のいずれかに記載の血流解析システムにおいて、
前記血管解析モデルは、三次元要素から生成されている
ことを特徴とする血流解析システム。
【請求項７】
請求項６に記載の血流解析システムにおいて、
前記分割処理部によって分割された血管ブロックを細線化し、当該血管ブロックごとに継ぎ目座標、細線化前の継ぎ目部分の直径および血管ブロックの長さを関連付ける細線化処理部を備え、
前記判別部は、細線化された血管ブロックの情報を利用して継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別する
ことを特徴とする血流解析システム。
【請求項８】
請求項４ないし請求項７のいずれかに記載の血流解析システムにおいて、
前記判別部によって判別された狭窄部位または閉塞部位を含む血管ブロックに装着可能なサイズのステントを抽出するステント抽出部を備え、
前記判別部は、抽出されたステントを装着した場合のシミュレーションを実行させて狭窄部位または閉塞部位の血流が所定値に改善されるか否かを判別する
ことを特徴とする血流解析システム。
【請求項９】
血管内の血液の流れをコンピュータによるシミュレーションによって解析する血流解析プログラムであって、
前記コンピュータに予め記憶されている被検体の血管の画像データを記憶部から読み出して血管解析モデルを生成する血管解析モデル生成部と、
設定された前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、当該所定位置から各毛細血管までの複数本のルート間の各圧力損失を求める算出部と、
前記算出部で求めた複数本の圧力損失と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させる
ことを特徴とする血流解析プログラム。
【請求項１０】
血管内の血液の流れをコンピュータによるシミュレーションによって解析する血流解析プログラムであって、
前記コンピュータに予め記憶されている被検体の血管の画像データを記憶部から読み出して血管解析モデルを生成する血管解析モデル生成部と、
設定された前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、当該所定位置から各毛細血管までの複数本のルート間の各流量を求める算出部と、
前記算出部で求めた複数本の流量と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させる
ことを特徴とする血流解析プログラム。
【請求項１１】
請求項９または請求項１０に記載の血流解析プログラムにおいて、
前記予め設定された大動脈圧と複数本のルートの遠位部の内圧の圧力損失、または、予め設定された冠動脈入口部の流量と複数本のルートの遠位部の流量の損失を求め、算出結果と予め決めた基準値を比較し、有意狭窄部位または閉塞部位を有する血管のルートの候補を判別する判別部を機能させる
ことを特徴とする血流解析プログラム。
【請求項１２】
請求項９ないし請求項１０のいずれかに記載の血流解析プログラムにおいて、
前記コンピュータに各ルートを所定間隔で複数個の血管ブロックに分割する分割処理部と、
前記血管ブロック同士の継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させる
ことを特徴とする血流解析プログラム。
【請求項１３】
請求項９ないし請求項１２のいずれかに記載の血流解析プログラムにおいて、
前記血管解析モデルは、細線化された一次元要素から生成されている
ことを特徴とする血流解析プログラム。
【請求項１４】
請求項９ないし請求項１２のいずれかに記載の血流解析プログラムにおいて、
前記血管解析モデルは、三次元要素から生成されている
ことを特徴とする血流解析プログラム。
【請求項１５】
請求項１４に記載の血流解析プログラムにおいて、
前記コンピュータに分割処理部によって分割された血管ブロックを細線化し、当該血管ブロックごとに継ぎ目座標、細線化前の継ぎ目部分の直径および血管ブロックの長さを関連付けて記憶部にさせる細線化処理部と、
細線化された前記血管ブロックの情報を利用して継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させる
ことを特徴とする血流解析プログラム。
【請求項１６】
請求項１２ないし請求項１５のいずれかに記載の血流解析プログラムにおいて、
前記判別部によって判別された狭窄部位または閉塞部位を含む血管ブロックに装着可能なサイズのステントを抽出するステント抽出部と、
前記ステントを装着した場合のシミュレーションを実行させて狭窄部位または閉塞部位の血流が所定値に改善されるか否かを判別する判別部とを機能させる
ことを特徴とする血流解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば冠動脈などの血管を流れる血液の状態を解析し、当該解析結果から血管内の狭窄部位および閉塞部位を判別することを可能にする血流解析システムおよび血流解析プログラムに関する。
【背景技術】
【０００２】
Ｘ線ＣＴ撮像装置などによって予め撮像された脳内血管の画像データを利用して動脈瘤の破裂の危険性をコンピュータのシミュレーションによって判断している。具体的には、血管部位の三次元形状データ、血管の厚み情報および病理情報など基づいて、流体解析および構造解析を行い、血管内の脳動脈瘤に作用する剪断応力ベクトルを求める。当該剪断応力ベクトルの向きおよび血流の状態量から動脈瘤の破裂の危険性を判断している（特許文献１を参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】国際公開公報ＷＯ２０１３−０３１７４１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、上記コンピュータのシミュレーションでは冠動脈などの血管の狭窄部位および閉塞部位の治療が必要な部位を判別することはできない。そこで、現状は次のようにして血管の狭窄部位および閉塞部位を判別している。血管の狭窄部位または閉塞部位は、複数箇所にわたって存在していることがある。そこで、例えば予め撮像した被検体の冠動脈の画像データから作成された冠動脈の三次元画像を、医師が目視してステント治療の必要な狭窄部位および閉塞部位を判別している。あるいは、狭窄部位と診断した箇所の遠位部までガイドワイヤを挿入し、当該ガイドワイヤに備わった圧力センサによって血管の内圧を測定し、主として各主要三分枝の遠位部の血管の内圧が大動脈（すなわち、冠動脈入口部位）に比べて血管拡張条件下でどの程度低下しているかを判断している。すなわち、冠血流予備能比を測定することによって治療適用の有無を判断している。
【０００５】
医師の目視による判断では、治療すべき狭窄部位および閉塞部位を過小評価したり、過大評価したりするおそれがある。また、圧力センサ付きガイドワイヤを利用する場合、大抵の施設で入院による心臓カテーテル検査が必要となり、患者への身体的、時間的および経済的な負担がかかるばかりでなく、検査自体に伴う脳梗塞や冠動脈穿孔などといった重篤な合併症のリスクが生じ、医療経済的にも非効率である。
【０００６】
本発明はこのような事情に鑑みてなされたものであって、冠動脈の治療が必要な狭窄部位および閉塞部位を非侵襲的に短時間で精度よく、数値的かつ客観的に判別することのできる血流解析システムおよび血流解析プログラムを提供することを主たる目的としている。
【課題を解決するための手段】
【０００７】
この発明は、このような目的を達成するために、次のような構成をとる。
【０００８】
すなわち、血管内の血液の流れをシミュレーションによって解析する血流解析システムであって、
撮像によって予め取得した被検体の血管の画像データを記憶する記憶部と、
前記記憶部から画像データを読み出して血管解析モデルを生成する血管解析モデル生成部と、
前記血管解析モデルを表示する表示部と、
前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件を設定する操作部と、
前記設定された各条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、所定位置から各毛細血管までの複数本のルート間の各圧力損失を求める算出部と、
前記算出部で求めた複数本の圧力損失と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部と、
を備えたことを特徴とする。
【０００９】
（作用・効果）  この血液解析システムによれば、生成された血管解析モデルの所定位置から毛細血管まで複数本のルートに、所定条件で血液を流した場合の各圧力損失がシミュレーションによって求められる。したがって、当該圧力損失と予め決めた基準値を比較することにより、狭窄部位または閉塞部位を含む血管のルートを判別することができる。すなわち、従来のように患者に負担をかけることなく、シミュレーションによるコンピュータ解析によって短時間で狭窄部位または閉塞部位を精度よく判別することができる。
【００１０】
また、この発明は、このような目的を達成するために、次のような構成をとる
【００１１】
すなわち、血管内の血液の流れをシミュレーションによって解析する血流解析システムであって、
撮像によって予め取得した被検体の血管の画像データを記憶する記憶部と、
前記記憶部から画像データを読み出して血管解析モデルを生成する血管解析モデル生成部と、
前記血管解析モデルを表示する表示部と、
前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件を設定する操作部と、
前記設定された各条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、起点から各毛細血管までの複数本のルート間を流れる血液の各流量を求める算出部と、
前記算出部で求めた複数本の流量と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部と、
を備えたことを特徴とする。
【００１２】
（作用・効果）  この血流解析システムによれば、生成された血管解析モデルの所定位置から毛細血管まで複数本のルートに、所定条件で血液を流した場合の各ルート間の各流量がシミュレーションによって求められる。したがって、当該流量と予め決めた基準値を比較することにより、狭窄部位または閉塞部位を含む血管のルートを判別することができる。すなわち、従来のように患者に負担をかけることなく、シミュレーションによるコンピュータ解析によって短時間で狭窄部位または閉塞部位を精度よく判別することができる。
【００１３】
なお、上記血流解析システムにおいて、判別部は、予め設定された大動脈圧と複数本のルートの遠位部の内圧の圧力損失、または、予め設定された冠動脈入口部の流量と複数本のルートの遠位部の流量の損失を求め、当該算出結果と予め決めた基準値を比較し、有意狭窄部位または閉塞部位を有する血管のルートの候補を判別することが好ましい。
【００１４】
この血流解析システムによれば、算出された圧力損失または流量の変化から複数本あるルートのうち治療の必要な責任病変（有意狭窄部位または閉塞部位）を含むルートの候補を精度よく判別することができる。例えば、血流などのパラメータを制御することで、血管拡張条件下、いわゆるhyperemia状態においての各主要三分枝の遠位部の血管内圧の大動脈圧に対する比率を求めることができる。圧力センサで測定したときに、治療適応の判断基準として用いられる冠血流予備能比の測定値（FFRmyo）が、例えば０．７５を下回るか否か、すなわち各主要三分枝の遠位部の血管内圧が大動脈圧に比し、例えば２５％以上の圧力損失を伴っているかを判別することができる。
【００１５】
上記各血流解析システムにおいて、各ルートを所定間隔で複数個の血管ブロックに分割する分割処理部を備え、
判別部は、血管ブロック同士の継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別することが好ましい。
【００１６】
この血流解析システムによれば、ルート内のどの部位に狭窄部位または閉塞部位が存在するかを精度よく特定することができる。
【００１７】
なお、上記各血流解析システムにおいて、血管解析モデルは、例えば、三次元要素または細線化された一次元要素のいずれであってもよい。
【００１８】
また、上記血流解析システムにおいて血管解析モデルが三次元要素の場合、分割処理部によって分割された血管ブロックを細線化し、当該血管ブロックごとに継ぎ目座標、細線化前の継ぎ目部分の直径および血管ブロックの長さを関連付ける細線化処理部を備え、
判別部は、細線化された血管ブロックの情報を利用して継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別することが好ましい。
【００１９】
この血流解析システムによれば、三次元の画像データから絞り込んだ特定の数値パラメータのみを利用して圧力損失を求めることができる。その結果、三次元の全画像データを利用する場合に比べて演算処理時間を短縮することができる。
【００２０】
さらに、上記血流解析システムにおいて、例えば、判別部によって判別された狭窄部位または閉塞部位を含む血管ブロックに装着可能なサイズのステントを抽出するステント抽出部を備え、
判別部は、抽出されたステントを装着した場合のシミュレーションを実行させて狭窄部位または閉塞部位の血流が所定値に改善されるか否かを判別することが好ましい。
【００２１】
この血流解析システムによれば、抽出された複数個のステントを狭窄部位または閉塞部位に装着した場合の血流の状態がシミュレーションによって解析されるので、最適なステントを選択することができる。すなわち、径、長さを指定したステントを狭窄部位、または閉塞部位に留置した状態をシミュレーションすることにより、治療後の血流、あるいは血管内圧の予測値を算出することができる。
【００２２】
また、この発明は、このような目的を達成するために、次のような構成をとる
【００２３】
すなわち、血管内の血液の流れをコンピュータによるシミュレーションによって解析する血流解析プログラムであって、
前記コンピュータに予め記憶されている被検体の血管の画像データを記憶部から読み出して血管解析モデルを生成する血管解析モデル生成部と、
設定された前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、当該所定位置から各毛細血管までの複数本のルート間の各圧力損失を求める算出部と、
前記算出部で求めた複数本の圧力損失と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させる
ことを特徴とする。
【００２４】
また、血管内の血液の流れをコンピュータによるシミュレーションによって解析する血流解析プログラムであって、
前記コンピュータに予め記憶されている被検体の血管の画像データを記憶部から読み出して血管解析モデルを生成する血管解析モデル生成部と、
設定された前記血管解析モデルの所定位置から下流にかけて分岐している複数本の毛細血管の位置、毛細血管の流動抵抗および血液の流入部および流出部の圧力、流量、流速のいずれかの境界条件に応じて、血管解析モデルに血液を流した場合のシミュレーションを実行し、当該所定位置から各毛細血管までの複数本のルート間の各流量を求める算出部と、
前記算出部で求めた複数本の流量と予め決めた基準値の比較から狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させる
ことを特徴とする。
【００２５】
（作用・効果）  これら血流解析プログラムによれば、上記血流解析システムを好適に実施させて、血管内の狭窄部位または閉塞部位を精度よく求めることができる。
【００２６】
なお、上記血流解析プログラムは、例えば、コンピュータに予め設定された大動脈圧と複数本のルートの遠位部の内圧の圧力損失、または、予め設定された冠動脈入口部の流量と複数本のルートの遠位部の流量の損失を求め、当該算出結果と予め決めた基準値を比較し、有意狭窄部位または閉塞部位を有する血管のルートの候補を判別させることが好ましい。
【００２７】
この血流解析プログラムによれば、血管の複数本のルートから有意狭窄部位または閉塞部位を含むルートの候補をコンピュータによって精度よく判別させることができる。
【００２８】
また、上記血流解析プログラムにおいて、コンピュータに各ルートを所定間隔で複数個の血管ブロックに分割する分割処理部と、
前記血管ブロック同士の継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させることが好ましい。
【００２９】
なお、上記各血流解析プログラムにおいて利用される血管解析モデルは、例えば、三次元要素または細線化された一次元要素のいずれであってもよい。
【００３０】
さらに、利用する血管解析モデルが三次元要素である場合、コンピュータに分割処理部によって分割された血管ブロックを細線化し、当該血管ブロックごとに継ぎ目座標、細線化前の継ぎ目部分の直径および血管ブロックの長さを関連付けて記憶部にさせる細線化処理部と、
細線化された前記血管ブロックの情報を利用して継ぎ目間の圧力損失を求め、当該圧力損失と基準値を比較して狭窄部位または閉塞部位を有している血管を判別する判別部とを機能させることが、さらに好ましい。
【００３１】
この血流解析プログラムによれば、血管の複数本のルートから判別された狭窄部位または閉塞部位を含むルート内のどの部位に狭窄部位などが存在するかをコンピュータによって精度よく判別させることができる。
【００３２】
なお、上記血流解析プログラムは、例えば、コンピュータに判別部によって判別された狭窄部位または閉塞部位を含む血管ブロックに装着可能なサイズのステントを抽出するステント抽出部と、
前記ステントを装着した場合のシミュレーションを実行させて狭窄部位または閉塞部位の血流が所定値に改善されるか否かを判別する判別部を機能させてもよい。
【００３３】
この血流解析プログラムによれば、狭窄部位などの血流を改善させる最適なステントをコンピュータによって抽出させることができる。
【発明の効果】
【００３４】
本発明の血流解析システムおよび血流解析プログラムによれば、血管の所定位置から毛細血管に分岐している複数本のルート間のどの部位に狭窄部位または閉塞部位が存在するかを短時間で精度よく判別することができる。
【図面の簡単な説明】
【００３５】
【図１】血流解析システムの概略全体構成を示すブロック図である。
【図２】血流解析システムを利用した血流解析の処理を示すフローチャートである。
【図３】三次元要素からなる血管解析モデルの画像を示す図である。
【図４】血管解析モデルの画像の分割処理を示す模式図である。
【図５】分割後の血管ブロックを細線化した画像を示す図である。
【図６】血管ブロックの圧力損失を求める模式図である。
【図７】変形例の血流解析システムの概略全体構成を示すブロック図である。
【図８】変形例の血流解析システムを利用した血流解析の処理を示すフローチャートである。
【図９】変形例の血流解析システムを利用した血流解析の処理を示すフローチャートである。
【発明を実施するための形態】
【００３６】
以下、図面を参照して本発明の一実施例を説明する。
【００３７】
図１に示すように、血流解析システムは、記憶部１、血管解析モデル生成部２、分割処理部３、細線化処理部４、算出部５、判別部６，表示部７および操作部８などから構成されている。以下、各部の構成については、図２に示すフローチャートの動作説明に沿って詳述する。
【００３８】
記憶部１は、被検体の冠動脈などの血管を、例えば、Ｘ線ＣＴ撮像装置、ＭＲＩ（磁気共鳴画像診断装置）およびＤＳＡ装置（デジタル差引血管造影法）などによって予め撮像された画像データを記憶している。また、記憶部１には、後述するように、操作部８から設定入力された各種条件、予め決めたレシピ、算出結果なども記憶する。
【００３９】
血管解析モデル生成部２が、記憶部１から予め記憶されている画像データを読み出し、三次元要素の血管解析モデルを生成する。例えば、本実施例では、図３に示すように、冠動脈の三次元要素の血管解析モデルを生成し、当該血管解析モデルの画像を表示部７に表示する（ステップＳ１）。
【００４０】
キーボード、マウスおよび操作パネルなどの操作部８を操作して各種条件を設定する。例えば、冠動脈の基端側の位置Ａ（起点）、毛細血管の流動抵抗および当該位置Ａから下流にかけて分岐している毛細血管の位置Ａ１、Ａ２…、心臓から冠動脈に送り出される血液の流入部および流出部の圧力、流量および流速のいずれか、および血管をブロック（区画）分割するピッチなどを設定する（ステップＳ２）。
【００４１】
ここで、毛細血管の流動抵抗、血液の圧力、流量および流速は、被検体の平常時および運動レベルに応じた複数のレシピが記憶部１に予め設定登録されており、適宜に設定を変更可能になっている。毛細血管の位置は、画面に表示された血管解析モデルを回転させながら選択してもよいし、指定した冠動脈の基端の血管の直径と表示領域内の下流側の血管の直径とを比較して求めてもよい。すなわち、冠動脈の基端と下流側の血管との直径を比較すると下流側が小さいので、当該直径の相対的な比較によって抽出することもできる。
【００４２】
なお、血管解析モデルには、毛細血管の末端までが生成されていないので、毛細血管の流動抵抗は、一定流量の血液が流れたときの圧力損失が所定の値になるように毛細血管の半径および距離を相対的に変化させた複数パターンのレシピを予め作成してある。
【００４３】
設定状態に基づいて、演算処理用のモデルデータが作成される。例えば、冠動脈の基端から複数箇所の毛細血管の先端までの複数本のルートＲ１〜Ｒｉ（ｉは正の整数）が作成される（ステップＳ３）。
【００４４】
分割処理部３が、図４の点線の楕円で示すように、設定された所定ピッチで各ルートを短小の血管ブロックＤ（Ｄ１、Ｄ２…）に分割する（ステップＳ４）。
【００４５】
分割処理が完了すると、図５に示すように、血管ブロックＤごとに中心線を求めて細線化する（ステップＳ５）。この細線化処理に伴い、隣接する血管ブロックＤ同士の継ぎ目の中心線の座標、継ぎ目の血管ブロックの直径または半径、血管ブロックＤの長さ（ピッチ）を血管ブロックＤごと関連付けて演算処理用のモデルデータを作成する。なお、モデルデータは、記憶部１に記憶される。
【００４６】
モデルデータの作成が完了すると、シミュレーションを実行し、血液の流れを解析する。すなわち、設定した条件で冠動脈の基端側（流入部）から各ルートの末端（流出部）まで血液を流すシミュレーションを実行する。例えば、流体解析法の一種である有限要素法などによってシミュレーションを実行し、血管ブロックＤごとに圧力損失を算出する（ステップＳ６）。
【００４７】
具体的には、算出部５が、各血管ブロックＤのモデルデータを読み出し、図６に示すように、基端側の圧力Ｐ１と終端側の圧力Ｐ２を求める。ここで、血管ブロックＤ内の圧力損失ΔＰは、ΔＰ＝Ｐ１−Ｐ２によって求まる。なお、図６は、狭窄部位を含む血管ブロックＤを例にとって説明している。
【００４８】
血管ブロックＤの圧力損失が求まるごとに、判別部６が、狭窄部位や閉塞部位などの有無を判別する（ステップＳ７）。すなわち、異常がない場合の血管ブロックＤの圧力損失の基準値と求めた圧力損失の算出値を比較する。比較の結果、判別部６は、算出値が基準値を超える場合、狭窄部位または閉塞部位と判別する。なお、狭窄部位と閉塞部位の判別は、予め決めた圧力損失のレベルに応じて判別する。
【００４９】
分割した全ての血管ブロックＤの判別処理が完了するまで、ステップＳ６からステップＳ７の処理が繰り返し実行される（ステップＳ８）。
【００５０】
上記一連の処理が完了し、狭窄部位または閉塞部位を含む血管ブロックＤが検出されると、検出結果が記憶部１に記憶される。当該記憶された検出結果は、血管解析モデルの対応する部位をハイカラーなどのマーキングなどを施して適宜に表示部７に表示する（ステップＳ９）。なお、表示画像は、例えば、三次元の血管解析モデルの画像上に狭窄部位を含む血管ブロックの箇所をマーキングしてもよい。または、細線化した画像上の部位をマーキングしてもよい。以上で、一連の処理が完了する。
【００５１】
なお、上記ステップＳ３からステップＳ４の処理は、コンピュータのプログラムによって実行される。
【００５２】
この構成によれば、三次元要素の血管解析モデルを細線化した血管ブロックＤごとに基端と終端の中心座標、血管の直径およびピッチを関連付けたデータ量の小さいモデルデータを作成し、当該モデルデータを利用して圧力損失を求めるので、三次元要素の全ての画像データを利用する場合に比べて演算処理の時間を大幅に短縮することができる。
【００５３】
また、血管ブロックごとに圧力損失を求めるので、狭窄部位または閉塞部位を精度よく特定することができる。
【００５４】
さらに、狭窄部位などは、Ｘ線ＣＴ撮像装置、ＭＲＩおよびＤＳＡ装置などによって予め撮像された画像データを用いたコンピューターシミュレーションによって判別されるので、患者には一切の負担がかからない。
【００５５】
本発明は上述した実施例のものに限らず、次のように変形実施することもできる。
【００５６】
（１）上記実施例装置は、以下のようなシミュレーションも実行することができる。例えば、図７に示すように、判別結果によって求まった狭窄部位などを含む血管ブロックＤの形状（例えば、直径、長さ）に応じたステントを抽出するステント抽出部９を上述の血流解析システムに付加した構成にする。この構成により、図８に示すフローチャートに沿ってシミュレーションを実行することができる。なお、ステップＳ１からステップＳ９までの処理は、上記実施例と同じであるので説明を省略し、追加されたステップＳ１０移行の処理について詳述する。
【００５７】
判別部６によって狭窄部位を含む血管ブロックＤが求まると、ステント抽出部９が、記憶部１から該当する血管ブロックＤのモデルデータを読み出し、当該血管ブロックＤの直径および長さに応じたステントを記録部１に記録されているステントタイプのデータベースから複数種類の候補を抽出する（ステップＳ１０）。
【００５８】
候補のステントごとに装着した血管の修正解析モデルを再作成する（ステップＳ１１）。修正解析モデルごとの血流状態を確認するために、先の判別時と同じ条件で血液を流したシミュレーションを実行し、装着部位の血管ブロックＤのみの圧力損失を求め（ステップＳ１２）、算出値と基準値との比較によって血流状態の改善レベルを判別する（ステップＳ１３）。全て候補のステントについての処理が繰り返されて（ステップＳ１４）終了すると、圧力損失の結果一覧を表示部７に表示させる（ステップＳ１５）。当該結果一覧から圧力損失が最小となるものなど最適条件となるステントを選択する。
【００５９】
この血管解析システムによれば、狭窄部位または閉塞部位の血流を改善するのに最適なステントを短時間で選択することができる。
【００６０】
（２）上記各実施例システムにおいて、狭窄部位または閉塞部位を特定する前に、これら狭窄部位などを含む血管の候補を複数本のルートから判別するように構成してもよい。例えば、上記各実施例のステップＳ７の判別処理において、次のような処理を実施すればよい。
【００６１】
予め設定された大動脈圧と複数本のルートの遠位部の内圧の圧力損失、または、予め設定された冠動脈入口部の流量と複数本のルートの遠位部の流量の損失を求め、当該算出結果と予め決めた基準値を比較し、有意狭窄部位または閉塞部位を有する血管のルートの候補を判別する。
【００６２】
例えば、血流などの複数個のパラメータを制御することにより、血管拡張条件下、いわゆるhyperemia状態においての各主要三分枝の遠位部の血管内圧の大動脈圧に対する比率を求めることができる。圧力センサ付きガイドワイヤで測定したときに、治療適応の判断基準として用いられる冠血流予備能比の測定値（FFRmyo：Fractional Flow Reserve myo）が、例えば０．７５を下回るか否か、すなわち各主要三分枝の遠位部の血管内圧が大動脈圧に比し、例えば２５％以上の圧力損失を伴っているかを判別することができる。
【００６３】
（３）上記各実施例システムにおいて、狭窄部位などを含む血管ブロックＤが求まった場合、血管ブロックＤのどの部位に含まれるかを、さらに解析するようシミュレーションを実行してもよい。
【００６４】
例えば、抽出された血管ブロックＤを所定ピッチで前後にシフトさせ、シフトごとにシミュレーションを実行し、当該血管ブロックＤのみの圧力損失の変化を求める。算出値が、最大値と最小までの間の偏差の変化量から位置を特定してもよい。あるいは、当該血管ブロックのピットを拡大および縮小したときのシミュレーションを実行して当該血管ブロックＤのみの圧力損失を求めるようにしてもよい。
【００６５】
（４）上記各実施例システムでは、血管ブロックＤごとに圧力損失を求めていたが、流量または流速の変化に基づいて、狭窄部位または閉塞部位を判別してもよい。すなわち、圧力、流量および流速は、流体のシミュレーションにおいて相対的な関係を有するので、これらのパラメータを予め設定することにより、いずれのパラメータを選択して実行しても、狭窄部位または閉塞部位を同じ精度で判別することができる。
【００６６】
すなわち。各血管ブロックの基端から終端を通過する血液の流量または流速をそれぞれ算出し、当該算出値と予め決めた基準値を比較することにより、狭窄部位または閉塞部位を判別することができる。
【００６７】
（５）上記各実施例では、三次元要素の血管解析モデルを利用していたが、一次元要素の血管解析モデルを利用してもよい。
【００６８】
例えば、図９に示すフローチャートに示すように、最初にＸ線ＣＴ画像データで作成される三次元画像を細線化し、一次元要素の血管解析モデルを作成する（ステップＳ１）。その後、上記実施例と同様に、条件設定し（ステップＳ２）、複数本の血管のルートを作成する（ステップＳ３）。ルートが作成されると、所定ピッチでルートを分割する。このとき、分割ブロックの節点ごとに座標および血管の直径または半径がそれぞれ関連付けられ、一次元要素のモデルデータを作成する（ステップＳ４）。以後、当該モデルデータを利用して分割ブロックごとに圧力損失求め（ステップＳ５）、当該算出値と基準値を比較して狭窄部位または閉塞部位を判別する（ステップＳ６）。当該圧力損失は、分割ブロックの全てに対して実施例される（ステップＳ８）。
【００６９】
当該実施例によれば、一次元要素の血管解析モデルのデータ量が、三次元要素の血管解析モデルのデータ量よりも小さいので、演算負荷が小さくなり、ひいては演算処理速度を向上させることができる。
【００７０】
なお、一次元要素の血管解析モデルを利用した場合であっても、上記変形例のように、判別結果によって狭窄部位が求まると、当該狭窄部位に応じたステントの候補を抽出し、血管修正モデルを再作成してシミュレーションを実行し、最適なステントを選択するようにもできる。
【００７１】
（６）上記実施例システムは、個々の病院内に設置したコンピュータのみによって構成してもよいし、個々の病院に設置された複数台のコンピュータと特定の場所に設置されたホストコンピュータをネットワークによって接続して構成してもよい。ネットワークによる構成の場合、個々の病院に設置したコンピュータから画像データを送信し、ホストコンピュータによって血流解析シミュレーションを実行した後に、当該結果を個々の病院のコンピュータにフィードバックするよう構成すればよい。
【符号の説明】
【００７２】
１  …  記憶部
２  …  血管解析モデル生成部
３  …  分割処理部
４  …  細線化処理部
５  …  算出部
６  …  判別部
７  …  表示部
８  …  操作部
９  …  ステント抽出部
Ｄ  …  血管ブロック
Ｒ  …  ルート
【図１】
2015171486.tif 000003
【図２】
2015171486.tif 000004
【図３】
2015171486.tif 000005
【図４】
2015171486.tif 000006
【図５】
2015171486.tif 000007
【図６】
2015171486.tif 000008
【図７】
2015171486.tif 000009
【図８】
2015171486.tif 000010
【図９】
2015171486.tif 000011