JP7286763B2 - 血管圧力差の取得装置の作動方法及び血管圧力差の取得装置 - Google Patents

Description

本発明は医療器械分野に関し、特に、血管断面関数、血管圧力差と血管応力の作成方法及び装置に関する。

人体血液中の脂類及び糖類物質の血管壁上での堆積によって、血管壁にプラークを形成し、続いて血管が狭くなってしまう。特に、心臓の冠動脈付近に発生した血管の狭窄化は、心筋への血液供給不足を招き、冠動脈性心臓病、狭心症等の病気を誘発させ、人類の健康に深刻な脅威を与えることになる。

多くの冠動脈生理機能評価技術において、冠血流予備量比（ＦＦＲ）は、現在最も正確な機能的な評価指標として認められている。

冠血流予備量比（ＦＦＲ）は、対象領域血管近端の血流圧力Ｐ_Ａ及び対象領域血管近端と遠端との間の血流圧力の差値△ｐから算出することで得たものである。

しかしながら、プラークの大きさ、プラークの長さ、プラーク形成の角度、プラークの形状及び形状の変化のいずれも血流圧力の差値△ｐの計算結果に影響を与えることになる。従来技術では、血流圧力の差値△ｐを計算する時、いずれもプラークの上記要素を考慮していないため、計算結果と実際値との誤差が大きい。

上記の課題を解決するために、対象領域血管の各位置における断面を取得し、内腔の異なる位置における形態を取得し、さらに、プラークの長さ、プラークの形状等の要素を考慮することが提案されている。

しかし、同一の位置における血管については、異なる時点での内腔の形態も全部が同じではない。したがって、計算結果には依然として一定の誤差が存在する。

上記の課題を解決するために、本発明は誤差が小さい血管断面関数の作成方法を提供することを目的とする。

本発明は、血管断面関数の作成方法であって、
少なくとも１つの心拍周期を含んだ画像データを取得するステップと、
１つの心拍周期内に、若干の特徴時刻を選択するステップと、
画像データに基づき、特徴時刻毎に対応する対象領域血管の空間モデルを生成するステップと、
空間モデル毎に基づき、対象領域血管の軸方向に沿って、対象領域血管の各位置における第一断面モデルを作成するステップと、
第一断面モデル毎に基づき、相応の第一断面関数を作成するステップと、を含む、
血管断面関数の作成方法を提供する。

さらに、
時系列に沿って、若干の特徴時刻は順次にｔ_１，ｔ_２，・・・，ｔ_ｑとなり、対象領域血管の軸方向に沿って、対象領域血管の近端から遠端までの各位置１，２，・・・，ｒの第一断面モデルに対応する第一断面関数が順次にｄ_１，ｄ_２，・・・，ｄ_ｒとなるようにするステップと、
第一断面関数と特徴時刻、位置との対応関係に基づき、第一断面マトリックスＡを作成するステップを含み、

式中、A_rqは特徴時刻t_q、対象領域血管の位置ｒであり、第一断面モデルに対応する第一断面関数はd_rとなり、

であり、かつ、ｑ、ｒは正の整数である。

さらに、第一断面マトリックスＡに基づき、列方向に沿って、特徴時刻毎の対象領域血管の二つの隣接する位置に対応する第一断面関数間の第一差異関数を作成するステップを含む。

さらに、第一差異関数に基づき、特徴時刻毎の対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第一変化関数ｆ(ｘ)を作成するステップを含む。

さらに、第一断面モデルは異なる尺度での第二断面モデルを含み、方法は、
尺度毎での第二断面モデルに基づき、相応の第二断面関数を作成するステップと、
第二断面関数と特徴時刻、位置との対応関係に基づき、異なる尺度での第二断面マトリックスを作成するステップと、を含み、
ここで、尺度が二つの隣接する位置間の距離である。

さらに、尺度は第一尺度、第二尺度、・・・、第ｎ尺度を含み、方法は、
尺度毎での第二断面マトリックスに基づき、列方向に沿って、異なる尺度での対象領域血管の二つの隣接する位置に対応する第二断面関数間の第二差異関数を作成するステップと、
第一尺度、第二尺度、・・・、第ｎ尺度での第二差異関数に基づき、対応する対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)を作成するステップと、を含み、
ここで、ｎは１より大きい正の整数である。

さらに、第一断面マトリックスＡに基づき、行方向に沿って、位置毎の対象領域血管の二つの隣接する特徴時刻に対応する第一断面関数間の第一差異関数を作成するステップを含む。

さらに、第一断面関数は面積関数、直径関数、周長関数及びエッジ距離関数のうちの１種又は複数種を含む。

本発明は、血管圧力差の取得方法であって、
上記血管断面関数の作成方法と、
画像データに基づき、対象領域血管の血流モデルを取得するステップと、
血流モデルに基づき、対象領域血管の血流速度Ｖを取得するステップと、
血流速度Ｖと第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)に基づき、特徴時刻毎の第一血管圧力差△ｐ_１を取得するステップと、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、第一血管圧力差△ｐ_１を重み付け加算し、第二血管圧力差△ｐ_２を取得するステップと、を含む、
血管圧力差の取得方法をさらに提供する。

さらに、第一血管圧力差△p_１の異なる尺度での計算式は、

であり、式中、
c₁，c₂，・・・，c_mはそれぞれ血流速度Ｖの参照係数であり、
α₁，α₂，・・・，α_nはそれぞれ異なる尺度での第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)の重み付け係数であり、
ｍ、ｎは正の整数である。

さらに、空間モデル毎に基づき、病変の冠動脈系における領域情報を取得するステップと、
領域情報に基づき、第一血管圧力差△ｐ_１を修正して、修正後の第三血管圧力差△ｐ_３を取得するステップと、を含み、
ここで、第三血管圧力差△ｐ_３と第一血管圧力差△ｐ_１は、
△ｐ_３＝ω×△ｐ_１の関係式を満たし、ωはずれ補正パラメータであり、かつ、

である。

さらに、領域情報は左冠動脈主幹部、左前下行枝、左回旋枝、右冠動脈及び分岐血管のうちの１種又は複数種を含む。

さらに、ずれ補正パラメータωは、以下の１種又は複数種を満たしている：
病変が左冠動脈主幹部にある場合、ずれ補正パラメータω＝１であり；
病変が左前下行枝にある場合、ずれ補正パラメータω＝０．９～１．０であり；
病変が左回旋枝にある場合、ずれ補正パラメータω＝０．６５～０．８５であり；
病変が右冠動脈にある場合、ずれ補正パラメータω＝０．７５～０．９であり；
病変が分岐血管にある場合、ずれ補正パラメータω＝０．５～０．８５である。

本発明は、血管応力の取得方法であって、
上記方法で血管断面関数を作成するステップと、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、同一の位置の二つの隣接する特徴時刻に対応する第一差異関数の第一絶対値関数を重み付け加算して、対象領域血管の各位置における変形を取得するステップと、を含む、
血管応力の取得方法をさらに提供する。

さらに、対象領域血管の各位置における変形に基づき、対象領域血管の二つの隣接する位置に対応する変形間の変形差異関数を取得するステップを含む。

さらに、変形差異関数の第二絶対値関数に基づき、累積的に加算して、対象領域血管の変形不均一さを取得するステップを含む。

本発明は、血管断面関数の作成装置であって、
少なくとも１つの心拍周期を含んだ画像データを取得するための画像取得モジュールと、
１つの心拍周期内に、若干の特徴時刻を選択するための時間選択モジュールと、
画像データに基づき、特徴時刻毎に対応する対象領域血管の空間モデルを生成するための空間モデル生成モジュールと、
空間モデル毎に基づき、対象領域血管の軸方向に沿って、対象領域血管の各位置における第一断面モデルを作成するための第一断面モデル作成モジュールと、
第一断面モデル毎に基づき、相応の第一断面関数を作成するための第一断面関数作成モジュールと、を備える、
血管断面関数の作成装置をさらに提供する。

さらに、時系列に沿って、若干の特徴時刻が順次にｔ_１，ｔ_２，・・・，ｔ_ｑとなり、対象領域血管の軸方向に沿って、対象領域血管の近端から遠端までの各位置１，２，・・・，ｒの第一断面モデルに対応する第一断面関数が順次にｄ_１，ｄ_２，・・・，ｄ_ｒとなることを備え、装置は、
第一断面関数と特徴時刻、位置との対応関係に基づき、第一断面マトリックスＡを作成するための第一断面マトリックス作成モジュールをさらに備え、

式中、A_rqは特徴時刻t_q、対象領域血管の位置ｒであり、第一断面モデルに対応する第一断面関数はd_rとなり、

であり、かつ、ｑ、ｒは正の整数である。

さらに、装置は、
第一断面マトリックスＡに基づき、列方向に沿って、特徴時刻毎の対象領域血管の二つの隣接する位置に対応する第一断面関数間の第一差異関数を作成するための第一差異関数作成モジュールをさらに備える。

さらに、装置は、
第一差異関数に基づき、特徴時刻毎の対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第一変化関数ｆ(ｘ)を作成するための第一変化関数作成モジュールをさらに備える。

さらに、第一断面モデルが異なる尺度での第二断面モデルを備え、装置は、
尺度毎での第二断面モデルに基づき、相応の第二断面関数を作成するための第二断面関数作成モジュールと、
第二断面関数と特徴時刻、位置との対応関係に基づき、異なる尺度での第二断面マトリックスを作成するための第二断面マトリックス作成モジュールと、をさらに備え、
ここで、尺度が二つの隣接する位置間の距離である。

さらに、尺度は第一尺度、第二尺度、・・・、第ｎ尺度を含み、装置は、
尺度毎での第二断面マトリックスに基づき、列方向に沿って、異なる尺度での対象領域血管の二つの隣接する位置に対応する第二断面関数間の第二差異関数を作成するための第二差異関数作成モジュールと、
第一尺度、第二尺度、・・・、第ｎ尺度での第二差異関数に基づき、対応する対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)を作成するための第二変化関数作成モジュールと、をさらに備え、
ここで、ｎは１より大きい正の整数である。

さらに、装置は、
第一断面マトリックスＡに基づき、行方向に沿って、位置毎の対象領域血管の二つの隣接する特徴時刻に対応する第一断面関数間の第一差異関数を作成するための第一差異関数作成モジュールをさらに備える。

さらに、第一断面関数は面積関数、直径関数、周長関数及びエッジ距離関数のうちの１種又は複数種を含む。

本発明は、血管圧力差の取得装置であって、
上記血管断面関数の作成装置と、
画像データに基づき、対象領域血管の血流モデルを取得するための血流モデル取得モジュールと、
血流モデルに基づき、対象領域血管の血流速度Ｖを取得するための血流速度取得モジュールと、
血流速度Ｖと第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)に基づき、取得特徴時刻毎の第一血管圧力差△ｐ_１を取得するための第一血管圧力差取得モジュールと、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、第一血管圧力差△ｐ_１を重み付け加算し、第二血管圧力差△ｐ_２を取得するための第二血管圧力差取得モジュールと、を備える、
血管圧力差の取得装置をさらに提供する。

さらに、第一血管圧力差△p_１の異なる尺度での計算式は、

であり、式中、
c₁，c₂，・・・，c_mはそれぞれ血流速度Ｖの参照係数であり、
α₁，α₂，・・・，α_nはそれぞれ異なる尺度での第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)の重み付け係数であり、
ｍ、ｎは正の整数である。

さらに、装置は、
空間モデル毎に基づき、病変の冠動脈系における領域情報を取得するための領域情報取得モジュールと、
領域情報に基づき、第一血管圧力差△ｐ_１を修正して、修正後の第三血管圧力差△ｐ_３を取得するための第三血管圧力差取得モジュールと、をさらに備え、
ここで、第三血管圧力差△ｐ_３と第一血管圧力差△ｐ_１は、
△ｐ_３＝ω＊△ｐ_１の関係式を満たし、ωはずれ補正パラメータであり、かつ、

である。

さらに、領域情報は左冠動脈主幹部、左前下行枝、左回旋枝、右冠動脈及び分岐血管のうちの１種又は複数種を含む。

さらに、ずれ補正パラメータωは、以下の１種又は複数種を満たしている：
病変が左冠動脈主幹部にある場合、ずれ補正パラメータω＝１であり；
病変が左前下行枝にある場合、ずれ補正パラメータω＝０．９～１．０であり；
病変が左回旋枝にある場合、ずれ補正パラメータω＝０．６５～０．８５であり；
病変が右冠動脈にある場合、ずれ補正パラメータω＝０．７５～０．９であり；
病変が分岐血管にある場合、ずれ補正パラメータω＝０．５～０．８５である。

本発明は、血管応力の取得装置であって、
上記血管断面関数の作成装置と、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、同一の位置の異なる特徴時刻の第一差異関数の第一絶対値関数を重み付け加算し、対象領域血管の各位置における変形を取得するための変形取得モジュールと、を備える、
血管応力の取得装置をさらに提供する。

さらに、装置は、
対象領域血管の各位置における変形に基づき、対象領域血管の二つの隣接する位置に対応する変形間の変形差異関数を取得するための変形差異関数取得モジュールをさらに備える。

さらに、装置は、
変形差異関数の第二絶対値関数に基づき、累積的に加算し、対象領域血管の変形不均一さを取得するための変形不均一さ取得モジュールをさらに備える。

したがって、本発明に提供される血管断面関数の作成方法は、異なる位置と異なる時刻に基づいて作成された血管断面関数であり、従来技術において空間のみに基づいて、又は時間のみに基づいて作成された血管断面関数と比べると、血管の実情をより精確に反映でき、その後の解析演算に誤差が小さい中間変数を提供し、後期の計算値を実際値にさらに接近させる。

本発明の上記内容をより明らかに分かり易くさせるために、ここに、好ましい実施例を挙げ、図面を参照しながら、次のように詳しく説明する。

以下、図面を参照しながら、本発明の具体的な実施形態について、さらに詳しく説明する。

図１は、本発明実施例中の血管断面関数の作成方法のフローを示す図である。 図２は、本発明実施例中の対象血管の一形態における空間モデルを示す図である。 図３は、図２中の第一断面モデルＤ１の構造を示す図である。 図４は、図２中の第一断面モデルＤ２の構造を示す図である。 図５は、図３と図４中の第一断面モデルＤ１と第一断面モデルＤ２の対応点が互いに重ね合わされた後の構造を示す図である。 図６は、本発明実施例中の対象血管の別の形態における空間モデルを示す図である。 図７は、図６中の第一断面モデルＤ３の構造を示す図である。 図８は、図６中の第一断面モデルＤ４の構造を示す図である。 図９は、図７と図８中の第一断面モデルＤ３と第一断面モデルＤ４の対応点が互いに重ね合わされた後の構造を示す図である。 図１０は、本発明実施例中の血管断面関数の作成装置の構造を示す図である。 図１１は、本発明実施例中の血管圧力差の取得装置の構造を示す図である。 図１２は、本発明実施例中の血管応力の取得装置の構造を示す図である。

以下、特定の具体的な実施例で本発明の実施形態を説明し、当業者は本明細書に開示された内容から本発明のその他の利点及び効果を容易に理解できる。本発明の内容は好ましい実施例を挙げて説明するが、これはこの発明の特徴が該実施形態のみに限定されることを表したものではない。まさしく逆であり、実施形態に基づいて発明を説明する目的は、本発明の特許請求の範囲に基づいて拡張可能なその他の選択又は改造をカバーするためである。本発明への深く理解を提供するために、以下の記載には多くの具体的な細部が含まれる。本発明もこれらの細部を使用せずに実施可能である。この他、本発明の重点を混乱や曖昧にさせるのを避けるために、ある具体的な細部は記載において省略される。説明が必要なのは、コンフリクトのない場合には、本発明中の実施例及び実施例中の特徴を互いに組み合わせることができることである。

図１を参照すると、本発明は血管断面関数の作成方法であって、
少なくとも１つの心拍周期を含んだ画像データを取得するステップと、
１つの心拍周期内に、若干の特徴時刻を選択するステップと、
画像データに基づき、特徴時刻毎に対応する対象領域血管の空間モデルを生成するステップと、
空間モデル毎に基づき、対象領域血管の軸方向に沿って、対象領域血管の各位置における第一断面モデルを作成するステップと、
各第一断面モデルに基づき、相応の第一断面関数を作成するステップと、を含む、
血管断面関数の作成方法を提供する。

説明が必要なのは、「特徴時刻」は、１つの心拍周期の画像データのうち、比較的にはっきりしていて、且つ解析に有利である各画像データに対応する各時刻点でもよいことである。病変血管があった場合には、「対象領域血管」は、少なくとも病変血管を含んだある領域の血管でもよい。「各位置」は、取得した画像データの解像度に基づき、対象領域血管が軸方向に沿って該解像度に対応する最小間隔で連続的に取得した各位置でもよい。「第一断面関数」は、対象領域血管断面の形態（内腔断面の形態）を反映でき、且つ取得した全ての物理量を算出できる関数でもよい。

実施例では、画像データは、ＣＴ、ＯＣＴ、ＩＶＵＳやＸ線などの検出手段により取得することができる。空間モデルは、少なくとも対象領域血管の形状、長さ、直径及び面積のようなパラメータを含む。空間モデルは、二次元、三次元等の多次元でもよいし、その多くとも到達可能な次元数は上記画像データの取得に用いられる検出手段に依存する。

血管にとって、同一の時刻点に、その各位置における収縮と弛緩の程度が一様でないため、内腔の形態がまちまちになってしまう。特に、病変が発生した血管、例えば、プラークを有する血管について、そのプラーク位置における血管とプラークのない位置における血管とでは、同一の時刻点で、その内腔の形態が全て同じではない。たとえ同じく一段のプラークを有する血管内でも、同一の時刻点で、その各位置における内腔の形態も各々同じではない。同様に、同一の位置での血管について、その異なる時刻点での内腔の形態が全部同じでもない。

したがって、本発明に提供される血管断面関数の作成方法は、異なる位置と異なる時刻に基づいて作成された血管断面関数であり、従来技術において空間のみに基づいて、又は時間のみに基づいて作成された血管断面関数と比べると、血管の実情をより精確に反映でき、その後の解析演算に誤差が小さい中間変数を提供し、後期の計算値を実際値にさらに接近させる。

本実施例では、第一断面モデルは、対象領域血管の近端の断面を参照面とし、参照面と対象領域血管の空間モデルの中心軸線との交点を参照点とするステップと、参照点を原点とし、中心軸線に沿って延びる方向をｘ軸とし、座標系を作成するステップと、中心軸線に垂直な方向に沿って対象領域血管の各位置に断面を作り、各断面の内外エッジを座標系に投影させ、対象領域血管の各位置での内腔断面の平面幾何図形を取得し、即ち、第一断面モデルを作成完了するステップによって作成することができる。その他の実施例では、第一断面モデルはその他の座標の作成と投影の方法で作成することができる。

説明が必要なのは、本願の「近端」は、対象領域血管の血流が先に流れる一端であり、「遠端」は、対象領域血管の血流が後に流れる一端であることである。

さらに、該血管断面関数の作成方法は、
時系列に沿って、若干の特徴時刻が順次にｔ_１，ｔ_２，・・・，ｔ_ｑとなり、対象領域血管の軸方向に沿って、対象領域血管の近端から遠端までの各位置１，２，・・・，ｒでの第一断面モデルに対応する第一断面関数が順次にｄ_１，ｄ_２，・・・，ｄ_ｒとなるようにするステップと、
第一断面関数と特徴時刻、位置との対応関係に基づき、第一断面マトリックスＡを作成するステップと、

をさらに含み、
式中、A_rqは特徴時刻t_q、対象領域血管の位置ｒであり、第一断面モデルに対応する第一断面関数はd_rとなり、

であり、かつ、ｑ、ｒは正の整数である。

本実施例では、第一断面関数は、面積関数、直径関数、周長関数及びエッジ距離関数のうちの１種又は複数種を含む。説明が必要なのは、面積関数は、内腔断面面積を表す関数でもよいことである。周長関数は、内腔断面周長を表す関数でもよい。直径関数は、内腔断面直径を表す関数でもよく、その中で、直径は、内腔断面が厳格な幾何学的な円形の直径でもよく、内腔断面が厳格な幾何学的な円形ではない場合の等価直径でもよい。

エッジ距離関数は、同一の特徴時刻に、ある位置における断面エッジ上の各点を別の隣接する位置における断面エッジ上の各点に１対１に対応させ、２つの断面エッジ上の対応点の間の直接距離又は間接距離を計算して取得し、更に加算又は平均することで得る、或いは、同一の位置に、ある特徴時刻に対応する断面エッジ上の各点を別の隣接する特徴時刻に対応する断面エッジ上の各点に１対１に対応させ、２つの断面エッジ上の対応点の間の直接距離又は間接距離を計算して取得し、更に加算又は平均することで得る、と理解できる。

説明が必要なのは、直接距離については、２つの断面を互いに重ね合わせ、２つの断面間の軸方向での距離を省略することで、２つの断面エッジ上の対応点の間の距離を直接計算して取得したものである、と理解できることである。もし２つの断面上の対応点の間の直接距離が共にゼロ、即ち２つの断面のエッジ距離関数が共にゼロであれば、該２つの断面で表される内腔の形態が完全一致であり、さもなくば、該２つの断面で表される内腔の形態が完全一致でない。内腔の形態が完全一致であることは、２つの内腔の対応点における内外輪郭が完全に重ね合わせ可能である、と理解できる。

間接距離については、参照位置を導入することで、２つの断面エッジ上の各点の該参照位置からの距離を計算して取得したものである、と理解できる。例えば、参照位置は、対象領域血管の空間モデルの中心軸線の位置でもよい。もし２つの断面エッジ上の各点の中心軸線からの間接距離が共にある固定値、即ち２つの断面のエッジ距離関数がともにある固定値であれば、該２つの断面で表される内腔の形態が完全一致であり、さもなくば、該２つの断面で表される内腔の形態が完全一致ではない。

さらに、該血管断面関数の作成方法は、
第一断面マトリックスＡに基づき、列方向に沿って、特徴時刻毎の対象領域血管の二つの隣接する位置に対応する第一断面関数間の第一差異関数を作成するステップをさらに含む。

さらに、該血管断面関数の作成方法は、
第一差異関数に基づき、特徴時刻毎の対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第一変化関数ｆ(ｘ)を作成するステップをさらに含む。

ここで、参照点は、即ち上記の対象領域血管の近端の断面を参照面とし、参照面と対象領域血管の空間モデルの中心軸線との交点である。

第一差異関数は、二つの隣接する位置での断面で表される内腔の形態差異を示す関数でもよい。本実施例では、第一断面関数は定数関数である。第一断面関数によって算出して得た第一差異関数も定数関数である。第一変化関数ｆ(ｘ)は、第一差異関数と対応の距離ｘからなる若干対の離散点に対してデータをフィッティングさせることによって得ることができる。

具体的には、図２を参照すると、対象領域血管内の血液が血流方向Ｙに沿って流れる。第一断面関数が面積関数となる時、対象領域血管の位置１と位置２での断面エッジ上の各点を１対１に対応させる。図３を参照し、対象領域血管の位置１での第一断面モデルＤ_１のプラークを有する領域をＡ_１とし、対応の面積をＳ_１とする。図４を参照し、対象領域血管の位置２の第一断面モデルＤ_２のプラークを有する領域をＡ_２とし、対応の面積をＳ_２とする。図５を参照し、Ｄ_１とＤ_２の内腔（プラーク）が重ね合わされないため、血流がＤ_１を経てＤ_２へ流れる時、血流圧力がそれに従って変化することになる。この時、第一差異関数は、２つの断面中の非重畳領域（Ｓ_１，Ｓ_２）と重畳領域面積（Ｓ_３）との割合、又は、非重畳領域（Ｓ_１，Ｓ_２）と総面積（Ｓ_１，Ｓ_２，Ｓ_３）との割合であり、かつ、第一変化関数ｆ(ｘ)>０、即ち、Ｄ_１とＤ_２との間には圧力差が存在する。

図６を参照すると、対象領域血管内の血液が血流方向Ｙに沿って流れる。図７を参照し、対象領域血管の位置３での第一断面モデルＤ_３のプラークを有する領域をＡ_３とし、対応の面積をＳ_３とする。図８を参照し、対象領域血管の位置４での第一断面モデルＤ_４のプラークを有する領域をＡ_４とし、対応の面積をＳ_４とする。ここで、Ｓ_４＝Ｓ_３である。図９を参照し、プラークを有する領域Ａ_３とＡ_４が完全に重ね合わせ可能であり、即ち、非重畳領域（Ｓ_３，Ｓ_４）＝０であり、第一差異関数がゼロであり、かつ、第一変化関数ｆ(ｘ)＝０であり、即ち、Ｄ_３とＤ_４との間には圧力差が存在しない。

同様に、第一断面関数がエッジ距離関数である場合、Ｄ_１とＤ_２との間の第一変化関数

であり、即ち、２つの内腔の形態が完全に重なり合っておらず、Ｄ_１とＤ_２との間には圧力差が存在しており、Ｄ_３とＤ_４との間の第一変化関数ｆ(ｘ)＝０であり、即ち、２つの内腔の形態が完全に重なり合っており、Ｄ_３とＤ_４との間には圧力差が存在しない。

第一差異関数は、第一断面マトリックスＡに基づき、列方向に沿って、特徴時刻毎の対象領域血管の二つの隣接する位置に対応する第一断面関数間の関数を作成し、又は、第一断面マトリックスＡに基づき、行方向に沿って、位置毎での対象領域血管の二つの隣接する特徴時刻に対応する第一断面関数間の関数を作成したものである。

つまり、第一差異関数は、（１）列方向に沿って、（２）行方向に沿ってという２つの異なるパスを介して計算して取得することができる。別の実施例では、第一差異関数は、別のパスを介して、例えば、対角線方向に沿って計算して取得することもできる。異なるパスを介して計算して取得した第一差異関数は同一になり得るし、異なり得る。具体的には、最終的に取得を必要とする計算量、例えば、血管圧力差、血流予備量比、血管応力等に応じて適切なパスを選択できる。

本実施例では、最終的に取得を必要とする計算量が血管圧力差又は血流予備量比である場合、第一差異関数は、第一断面マトリックスＡに、列方向に沿って、即ち特徴時刻が同一で、位置が異なる方向に沿って計算することで取得できる。最終的に取得を必要とする計算量が血管応力である場合、第一差異関数は、第一断面マトリックスＡに、行方向に沿って、即ち位置が同一で、特徴時刻が異なる方向に沿って計算することで取得できる。

さらに、第一断面モデルは異なる尺度での第二断面モデルを備える。

該血管断面関数の作成方法は、
尺度毎での第二断面モデルに基づき、相応の第二断面関数を作成するステップと、
第二断面関数と特徴時刻、位置との対応関係に基づき、異なる尺度での第二断面マトリックスを作成するステップと、をさらに含み、
ここで、尺度は二つの隣接する位置の間の距離である。

さらに、尺度は第一尺度、第二尺度、・・・、第ｎ尺度を含む。

該血管断面関数の作成方法は、
尺度毎での第二断面マトリックスに基づき、列方向に沿って、異なる尺度での対象領域血管の二つの隣接する位置に対応する第二断面関数間の第二差異関数を作成するステップと、
第一尺度、第二尺度、・・・、第ｎ尺度での第二差異関数に基づき、対応する対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)を作成するステップと、をさらに含み、
ここで、ｎは１より大きい正の整数である。

第一尺度、第二尺度、・・・、第ｎ尺度は、異なる病変特徴の検出に用いるものでもよい。例えば、第一種の病変特徴は比較的大きい尺度の第一尺度で検出する必要があり、第ｎ種の病変特徴は比較的小さい尺度の第ｎ尺度で検出する必要がある。異なる尺度により、より精確に各種の病変を検出でき、計算結果と実際値との誤差をさらに減少できる。

本発明は、血管圧力差の取得方法であって、
上記血管断面関数の作成方法と、
画像データに基づき、対象領域血管の血流モデルを取得するステップと、
血流モデルに基づき、対象領域血管の血流速度Ｖを取得するステップと、
血流速度Ｖと第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)に基づき、特徴時刻毎での第一血管圧力差△ｐ_１を取得するステップと、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、第一血管圧力差△ｐ_１を重み付け加算し、第二血管圧力差△ｐ_２を取得するステップと、を含む、
血管圧力差の取得方法をさらに提供する。

ここで、各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合は、例えば、１つの心拍周期内に４つの特徴時刻を選択し、時系列に沿って順次にt₁、t₂、t₃、t₄とし、１つの心拍周期をＴとすると、t₁の１つの心拍周期に占める割合＝(t₂-t₁)/Tであり、t₂の１つの心拍周期に占める割合＝(t₃-t₂)/Tであり、その他はこれによって類推することができる。又は、t₂の１つの心拍周期に占める割合＝(t₂-t₁)/Tであり、t₃の１つの心拍周期に占める割合＝(t₃-t₂)/Tであり、その他はこれによって類推することができる。若しくは、t₁の１つの心拍周期に占める割合＝{1/2(t₂-t₁)+1/2(t₃-t₂)}/Tであり、t₂の１つの心拍周期に占める割合＝{1/2(t₃-t₂)+1/2(t₄-t₃)}/Tであり、その他はこれによって類推することができる、と理解できる。

血流モデルは、固定血流モデルと個性化血流モデルを含む。固定血流モデルは、経験値血流モデルであり、実際の臨床経験に基づき、ビッグデータによる採集及び模擬の方法によって直接作成されたものである。血流速度Ｖは、固定血流モデルから直接取得することができる。血流速度Ｖは固定の数値でもよい。

個性化血流モデルは、安静時血流モデルと負荷時血流モデルをさらに含む。個性化血流モデルが安静時血流モデルである場合、血流速度Ｖは流体が充満状態での速度を計算することで得ることができる。

一実施例では、安静時血流モデルは造影剤血流モデルである。血流速度Ｖはグレースケール時間フィッティング関数を利用して、対象領域血管の造影プロセスにおける造影剤の平均血流速度Ｖ_ｑｃを得ることができる。又は、ＴＩＭＩフレームカウントにより計算し、対象領域血管の造影プロセスにおける造影剤の平均血流速度Ｖ_ｑｃを得ることができる。

別の実施例では、安静時血流モデルはＣＴ血流モデルである。血流速度Ｖは血管樹の形態によって計算して得ることができる。血管樹の形態は少なくとも血管樹の面積、体積及び血管樹中の血管セグメントの内腔直径のうちの１種又は複数種を含む。空間モデルはさらに少なくとも灌注面積と分岐血管角度というパラメータのうちの１種又は複数種を含まなければならない。

個性化血流モデルが負荷状態血流モデルである場合、血流速度Ｖはアデノシルを注射することで、血管が十分に拡張された後の最大血流速度Ｖ_ｍａｘである。

一般的に言うと、対象領域血管が冠動脈領域に位置する場合、血流速度Ｖは最大血流速度Ｖ_ｍａｘとなり、対象領域血管が外周血管系に位置する場合、血流速度Ｖは安静時の平均血流速度Ｖ_ｑｃとなる。

ここで、第一血管圧力差△p_１の異なる尺度での計算式は、

であり、式中、
c₁，c₂，・・・，c_mはそれぞれ血流速度Ｖの参照係数であり、
α₁，α₂，・・・，α_nはそれぞれ異なる尺度での第二変化関数f₁(x)，f₂(x)，・・・，f_n(x)の重み付け係数であり、
ｍ、ｎは正の整数である。

ここで、参照係数c₁，c₂，・・・，cmは血液の粘度による影響要素、血液の乱流による影響要素及び粘性係数等の複数のパラメータ係数を含む。本実施では、ｍ＝２であり、かつ、c₁は血液の流動摩擦によるパラメータ係数であり、c₂は血液の乱流によるパラメータ係数である。

さらに、血管圧力差の取得方法は、
空間モデル毎に基づき、病変が冠動脈系にある領域情報を取得するステップと、
領域情報に基づき、第一血管圧力差△ｐ_１を修正し、修正後の第三血管圧力差△ｐ_３を取得するステップと、をさらに含み、
ここで、第三血管圧力差△ｐ_３と第一血管圧力差△ｐ_１は、
△ｐ_３＝ω＊△ｐ_１の関係式を満たし、ωはずれ補正パラメータであり、かつ、

である。

病変（プラーク）のある位置が異なるため、対象領域血管の供給を受ける心筋体積の領域が異なるから、最終的に計算して得られた対象領域血管圧力差の数値のずれが発生してしまうことになる。したがって、本発明はずれ補正パラメータωを導入することで、病変（プラーク）のある箇所の位置要素の圧力差の計算への影響をさらには減らし、計算結果の正しさを高める。

本実施では、領域情報は、左冠動脈主幹部、左前下行枝、左回旋枝、右冠動脈及び分岐血管のうちの１種又は複数種を含む。

ずれ補正パラメータωは、以下の１種又は複数種を満たしており：
病変が左冠動脈主幹部にある場合、ずれ補正パラメータω＝１であり；
病変が左前下行枝にある場合、ずれ補正パラメータω＝０．９～１．０であり；
病変が左回旋枝にある場合、ずれ補正パラメータω＝０．６５～０．８５であり；
病変が右冠動脈にある場合、ずれ補正パラメータω＝０．７５～０．９であり
病変が分岐血管にある場合、ずれ補正パラメータω＝０．５～０．８５である。

説明が必要なのは、上記病変（プラーク）のあるところの位置要素が対象領域血管圧力差の計算結果に影響を及ぼすことになる以外、その他、例えば、病歴情報、生理情報も圧力差の計算結果に程度が異なる影響を及ぼすことになることである。その他の実施例では、病歴情報や生理情報に対する偏差除去処理を行うこともできる。そのうち、病歴情報は、血流の流速又は血液の粘度に影響を与える循環系の病気、呼吸器の病気、神経系の病気、骨格の病気、消化器の病気、代謝の病気や家族史などを含む。生理情報は、年齢、性別、血圧、体質指数や冠動脈の優勢類型など直接又は間接的に取得可能な情報を含む。

上記方法で取得した圧力差に基づき、冠血流予備量比ＦＦＲをさらに計算して取得することができる。

冠血流予備量比ＦＦＲの取得方法は、少なくとも、
対象領域血管の近端圧力値ｐ_ａを取得するステップと、
近端圧力値ｐ_ａと第二血管圧力差△ｐ_２に基づき、冠血流予備量比ＦＦＲを取得するステップと、を含み、
ここで、冠血流予備量比ＦＦＲは、下記の式を満たしている。
FFR=(p_a-△ｐ_２)/p_a

本発明は、血管応力の取得方法であって、
上記方法に基づき、血管断面関数を作成するステップと、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、同一の位置における二つの隣接する特徴時刻に対応する第一差異関数の第一絶対値関数を重み付け加算し、対象領域血管の各位置における変形を取得するステップと、を含む、
血管応力の取得方法をさらに提供する。

説明が必要なのは、第一差異関数が第一断面マトリックスＡに、行方向に沿って、即ち位置が同一で、特徴時刻が異なる方向に沿って計算することで得られることである。第一絶対値関数は第一差異関数に対して絶対値を取って得られた関数である。

さらに、血管応力の取得方法は、
対象領域血管の各位置における変形に基づき、対象領域血管の二つの隣接する位置に対応する変形間の変形差異関数を取得するステップをさらに含んでもよい。

さらに、血管応力の取得方法は、
変形差異関数の第二絶対値関数に基づき、累積的に加算し、対象領域血管の変形不均一さを取得するステップをさらに含んでもよい。

説明が必要なのは、第二絶対値関数は変形差異関数に対して絶対値を取って得られた関数であることである。

上記方法で得られた応力は、従来技術において空間のみに基づいて、又は時間のみに基づいて計算して得られた応力と比べると、実際値にさらに接近させるので、続いてくる評価及び意思決定の正確さを高めることができる。例えば、正常又は病変の血管壁のリアルタイム変形能力及び力学行為を評価し、正常又は病変の血管壁プラークのひび割れのリスクを評価し、ステント埋め込み手術において、ステントの所在する位置やピンチアウトの程度の心拍周期における病変のプラークの変位、内腔解剖形態への影響、及びステント断裂の潜在的なリスクを評価する。

図１０を参照し、本発明は、血管断面関数の作成装置１０であって、
少なくとも１つの心拍周期を含んだ画像データを取得するための画像取得モジュール１１と、
１つの心拍周期内に、若干の特徴時刻を選択するための時間選択モジュール１２と、
画像データに基づき、特徴時刻毎に対応する対象領域血管の空間モデルを生成するための空間モデル生成モジュール１３と、
空間モデル毎に基づき、対象領域血管の軸方向に沿って、対象領域血管の各位置における第一断面モデルを作成するための第一断面モデル作成モジュール１４と、
第一断面モデル毎に基づき、相応の第一断面関数を作成するための第一断面関数作成モジュール１５と、を備える、
血管断面関数の作成装置１０をさらに提供する。

さらに、該血管断面関数の作成装置１０は、
時系列に沿って、若干の特徴時刻が順次にｔ_１，ｔ_２，・・・，ｔ_ｑとなり、対象領域血管の軸方向に沿って、対象領域血管の近端から遠端までの各位置１，２，・・・，ｒでの第一断面モデルに対応する第一断面関数が順次にｄ_１，ｄ_２，・・・，ｄ_ｒとなり、
第一断面関数と特徴時刻、位置との対応関係に基づき、第一断面マトリックスＡを作成するための第一断面マトリックス作成モジュールをさらに備え、

式中、A_rqは特徴時刻t_q、対象領域血管の位置ｒであり、第一断面モデルに対応する第一断面関数はd_rとなり、

であり、かつ、ｑ、ｒは正の整数である。

さらに、血管断面関数の作成装置１０は、
第一断面マトリックスＡに基づき、列方向に沿って、特徴時刻毎の対象領域血管の二つの隣接する位置に対応する第一断面関数間の第一差異関数を作成する、又は、行方向に沿って、位置毎の対象領域血管の二つの隣接する特徴時刻に対応する第一断面関数間の第一差異関数を作成するための第一差異関数作成モジュールをさらに備える。

さらに、血管断面関数の作成装置１０は、
第一差異関数に基づき、特徴時刻毎の対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第一変化関数ｆ(ｘ)を作成するための第一変化関数作成モジュールをさらに備える。

さらに、第一断面モデルは異なる尺度での第二断面モデルを備える。

該血管断面関数の作成装置１０は、
尺度毎での第二断面モデルに基づき、相応の第二断面関数を作成するための第二断面関数作成モジュールと、
第二断面関数と特徴時刻、位置との対応関係に基づき、異なる尺度での第二断面マトリックスを作成するための第二断面マトリックス作成モジュールと、を更に備え、
ここで、尺度は二つの隣接する位置の間の距離である。

さらに、尺度は第一尺度、第二尺度、・・・、第ｎ尺度を含む。

該血管断面関数の作成装置１０は、
尺度毎での第二断面マトリックスに基づき、列方向に沿って、異なる尺度での対象領域血管の二つの隣接する位置に対応する第二断面関数間の第二差異関数を作成するための第二差異関数作成モジュールと、
第一尺度、第二尺度、・・・、第ｎ尺度での第二差異関数に基づき、対応する対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)を作成するための第二変化関数作成モジュールと、をさらに備え、
ここで、ｎは１より大きい正の整数である。

図１１を参照し、本発明は、血管圧力差の取得装置であって、
上記血管断面関数の作成装置１０と、
画像データに基づき、対象領域血管の血流モデルを取得するための血流モデル取得モジュール２０と、
血流モデルに基づき、対象領域血管の血流速度Ｖを取得するための血流速度取得モジュール３０と、
血流速度Ｖと第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)に基づき、特徴時刻毎での第一血管圧力差△ｐ_１を取得するための第一血管圧力差取得モジュール４０と、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、第一血管圧力差△ｐ_１を重み付け加算し、第二血管圧力差△ｐ_２を取得するための第二血管圧力差取得モジュール５０と、を備える、
血管圧力差の取得装置をさらに提供する。

さらに、血管圧力差の取得装置は、
空間モデル毎に基づき、病変が冠動脈系にある領域情報を取得するための領域情報取得モジュールと、
領域情報に基づき、第一血管圧力差△ｐ_１を修正し、修正後の第三血管圧力差△ｐ_３を取得するための第三血管圧力差取得モジュールと、をさらに備え、
ここで、第三血管圧力差△ｐ_３と第一血管圧力差△ｐ_１は、
△ｐ_３＝ω＊△ｐ_１の関係式を満たし、ωはずれ補正パラメータであり、

である。

図１２を参照し、本発明は、血管応力の取得装置であって、
上記血管断面関数の作成装置１０と、
各特徴時刻と隣接する特徴時刻との時間間隔が１つの心拍周期に占める割合に基づき、同一の位置に異なる特徴時刻の第一差異関数の第一絶対値関数を重み付け加算し、対象領域血管の各位置における変形を取得するための変形取得モジュール６０と、を備える、
血管応力の取得装置をさらに提供する。

さらに、血管応力の取得装置は、
対象領域血管の各位置における変形に基づき、対象領域血管の二つの隣接する位置に対応する変形間の変形差異関数を取得するための変形差異関数取得モジュール７０をさらに備える。

さらに、血管応力の取得装置は、
変形差異関数の第二絶対値関数に基づき、累積的に加算し、対象領域血管の変形不均一さを取得するための変形不均一さ取得モジュール８０をさらに備える。

各装置の具体的な実施について、上記相応の方法を参照されたい。ここでは、詳細を再度記載しない。

上記のように、本発明が提供した上記実施例はただ例示的に本発明の原理及びその効果を説明しただけであり、本発明を制限するためのものではない。当業者は本発明の主旨や範囲に違反せず、上記実施例を修飾又は変更することができる。したがって、本発明に開示された主旨及び技術思想から逸脱せずに当業者によってなされる全ての等価の修飾又は変更が、やはり本発明の請求の範囲に包括されるべきである。

Claims (10)

1. 血管圧力差の取得装置の作動方法において、
   少なくとも１つの心拍周期を含んだ画像データを取得するステップと、
   前記１つの心拍周期内に、若干の特徴時刻を選択するステップと、
   前記画像データに基づき、前記特徴時刻毎に対応する対象領域血管の空間モデルを生成するステップと、
   前記空間モデル毎に基づき、前記対象領域血管の軸方向に沿って、前記対象領域血管の各位置における第一断面モデルを作成するステップと、
   前記第一断面モデル毎に基づき、相応の第一断面関数を作成するステップと、
   時系列に沿って、前記若干の特徴時刻は順次にｔ_１，ｔ_２，・・・，ｔ_ｑとなり、前記対象領域血管の前記軸方向に沿って、前記対象領域血管の近端から遠端までの各前記位置１，２，・・・，ｒの前記第一断面モデルに対応する第一断面関数が順次にｄ_１，ｄ_２，・・・，ｄ_ｒとなるようにするステップと、
   前記第一断面関数と前記特徴時刻、前記位置との対応関係に基づき、第一断面マトリックスＡを作成するステップと、
   

   

   
   を含み、
   式中、A_rqは、特徴時刻がt_qであり、かつ、前記対象領域血管の前記位置がｒであると
   きの前記第一断面モデルの第一断面関数d_rであり、
   

   

   
   であり、かつ、ｑ、ｒは正の整数であり、
   前記第一断面モデルは異なる尺度での第二断面モデルを備え、尺度毎での前記第二断面モデルに基づき、相応の第二断面関数を作成し、前記第二断面関数と前記特徴時刻、前記位置との対応関係に基づき、異なる尺度での第二断面マトリックスを作成し、ここで、前記尺度は二つの隣接する前記位置間の距離であり、前記尺度は、第一尺度、第二尺度、・・・、第ｎ尺度を備え、
   尺度毎での前記第二断面マトリックスに基づき、列方向に沿って、前記異なる尺度での前記対象領域血管の二つの隣接する前記位置に対応する前記第二断面関数間の第二差異関数を作成し、
   前記第一尺度、第二尺度、・・・、第ｎ尺度での前記第二差異関数に基づき、対応する前記対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)を作成し、
   ここで、前記ｎは１より大きい正の整数であり、
   前記画像データに基づき、前記対象領域血管の血流モデルを取得し、
   前記血流モデルに基づき、前記対象領域血管の血流速度Ｖを取得し、
   前記血流速度Ｖと前記第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)に基づき、前記特徴時刻毎の第一血管圧力差△ｐ_１を取得し、
   各前記特徴時刻と隣接する前記特徴時刻との時間間隔が前記１つの心拍周期に占める割合に基づき、前記第一血管圧力差△ｐ_１を重み付け加算し、第二血管圧力差△ｐ_２を取得することを特徴とする、
   血管圧力差の取得装置の作動方法。
2. 前記第一血管圧力差△p_１の前記異なる尺度での計算式は、
   

   

   
   であり、式中、
   c₁，c₂，・・・，c_mはそれぞれ前記血流速度Vの参照係数であり、
   α₁，α₂，・・・，α_nはそれぞれ前記異なる尺度での前記第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)の重み付け係数であり、
   m、nは正の整数であることを特徴とする、
   請求項１に記載の血管圧力差の取得装置の作動方法。
3. 前記空間モデル毎に基づき、処理が必要である冠動脈系における領域情報を取得するステップと、
   前記領域情報に基づき、第一血管圧力差△ｐ_１を修正し、修正後の第三血管圧力差△ｐ_３を取得するステップと、をさらに含み、
   ここで、前記第三血管圧力差△ｐ_３と前記第一血管圧力差△ｐ_１は、
   △ｐ_３＝ω＊△ｐ_１の関係式を満たし、ωはずれ補正パラメータであり、かつ、
   

   

   
   であることを特徴とする、
   請求項１に記載の血管圧力差の取得装置の作動方法。
4. 前記領域情報は、左冠動脈主幹部、左前下行枝、左回旋枝、右冠動脈及び分岐血管のうちの１種又は複数種を含むことを特徴とする、請求項３に記載の血管圧力差の取得装置の作動方法。
5. 左冠動脈主幹部を処理する必要がある場合、前記ずれ補正パラメータω＝１であり、
   左前下行枝を処理する必要がある場合、前記ずれ補正パラメータω＝０．９～１．０であり、
   左回旋枝を処理する必要がある場合、前記ずれ補正パラメータω＝０．６５～０．８５であり、
   右冠動脈を処理する必要がある場合、前記ずれ補正パラメータω＝０．７５～０．９であり、
   分岐血管を処理する必要がある場合、前記ずれ補正パラメータω＝０．５～０．８５であり、
   前記ずれ補正パラメータωは、以上の１種又は複数種を満たしていることを特徴とする、
   請求項４に記載の血管圧力差の取得装置の作動方法。
6. 血管圧力差の取得装置において、
   少なくとも１つの心拍周期を含んだ画像データを取得するための画像取得モジュールと、
   前記１つの心拍周期内に、若干の特徴時刻を選択するための時間選択モジュールと、
   前記画像データに基づき、前記特徴時刻毎に対応する対象領域血管の空間モデルを生成するための空間モデル生成モジュールと、
   前記空間モデル毎に基づき、前記対象領域血管の軸方向に沿って、前記対象領域血管の各位置における第一断面モデルを作成するための第一断面モデル作成モジュールと、
   前記第一断面モデル毎に基づき、相応の第一断面関数を作成するための第一断面関数作成モジュールと、を備え、時系列に沿って、前記若干の特徴時刻は順次にｔ_１，ｔ_２，・・・，ｔ_ｑとなり、前記対象領域血管の前記軸方向に沿って、前記対象領域血管の近端から遠端までの各前記位置１，２，・・・，ｒでの前記第一断面モデルに対応する第一断面関数は順次にｄ_１，ｄ_２，・・・，ｄ_ｒとなり、
   前記第一断面関数と前記特徴時刻、前記位置との対応関係に基づき、第一断面マトリックスＡを作成するための第一断面マトリックス作成モジュールをさらに備え、
   

   

   
   式中、A_rqは、特徴時刻がt_qであり、かつ、前記対象領域血管の前記位置がｒであると
   きの前記第一断面モデルの第一断面関数d_rであり、
   

   

   
   であり、かつ、ｑ、ｒは正の整数であり、
   前記第一断面モデルは異なる尺度での第二断面モデルを備え、
   尺度毎での前記第二断面モデルに基づき、相応の第二断面関数を作成するための第二断面関数作成モジュールと、
   前記第二断面関数と前記特徴時刻、前記位置との対応関係に基づき、異なる尺度での第二断面マトリックスを作成するための第二断面マトリックス作成モジュールと、をさらに備え、
   ここで、前記尺度は二つの隣接する前記位置間の距離であり、
   前記尺度は第一尺度、第二尺度、・・・、第ｎ尺度を備え、
   尺度毎での前記第二断面マトリックスに基づき、列方向に沿って、前記異なる尺度での前記対象領域血管の二つの隣接する前記位置に対応する前記第二断面関数間の第二差異関数を作成するための第二差異関数作成モジュールと、
   前記第一尺度、第二尺度、・・・、第ｎ尺度での前記第二差異関数に基づき、対応する前記対象領域血管のそのいずれか１つの位置から参照点までの距離ｘの変化に従って変化する第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)を作成するための第二変化関数作成モジュールと、をさらに備え、
   ここで、前記ｎは１より大きい正の整数であり、前記画像データに基づき、前記対象領域血管の血流モデルを取得するための血流モデル取得モジュールと、
   前記血流モデルに基づき、前記対象領域血管の血流速度Ｖを取得するための血流速度取得モジュールと、
   前記血流速度Vと前記第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)に基づき、前記特徴時刻毎の第一血管圧力差△ｐ_１を取得するための第一血管圧力差取得モジュールと、
   各前記特徴時刻と隣接する前記特徴時刻との時間間隔が前記１つの心拍周期に占める割合に基づき、前記第一血管圧力差△ｐ_１を重み付け加算し、第二血管圧力差△ｐ_２を取得するための第二血管圧力差取得モジュールと、を備えることを特徴とする、
   血管圧力差の取得装置。
7. 前記第一血管圧力差△p₁の前記異なる尺度での計算式は、
   

   

   
   であり、式中、
   c₁，c₂，・・・，c_mはそれぞれ前記血流速度Vの参照係数であり、
   α₁，α₂，・・・，α_nはそれぞれ前記異なる尺度での前記第二変化関数f_１(x)，f_２(x)，・・・，f_n(x)の重み付け係数であり、
   m、nは正の整数であることを特徴とする、
   請求項６に記載の血管圧力差の取得装置。
8. 前記空間モデル毎に基づき、病変の冠動脈系における領域情報を取得するための領域情報取得モジュールと、
   前記領域情報に基づき、第一血管圧力差△ｐ_１を修正し、修正後の第三血管圧力差△ｐ_３を取得するための第三血管圧力差取得モジュールと、をさらに備え、
   ここで、前記第三血管圧力差△ｐ_３と前記第一血管圧力差△ｐ_１が、
   △ｐ_３＝ω＊△ｐ_１の関係式を満たし、ωはずれ補正パラメータであり、かつ、
   

   

   
   であることを特徴とする、
   請求項６に記載の血管圧力差の取得装置。
9. 前記領域情報は左冠動脈主幹部、左前下行枝、左回旋枝、右冠動脈及び分岐血管のうちの１種又は複数種を含むことを特徴とする、請求項８に記載の血管圧力差の取得装置。
10. 前記左冠動脈主幹部を処理する必要がある場合、前記ずれ補正パラメータω＝１であり、
    前記左前下行枝を処理する必要がある場合、前記ずれ補正パラメータω＝０．９～１．０であり、
    前記左回旋枝を処理する必要がある場合、前記ずれ補正パラメータω＝０．６５～０．８５であり、
    前記右冠動脈を処理する必要がある場合、前記ずれ補正パラメータω＝０．７５～０．９であり、
    前記分岐血管を処理する必要がある場合、前記ずれ補正パラメータω＝０．５～０．８５であり、
    前記ずれ補正パラメータωは、以上の１種又は複数種を満たしていることを特徴とする、
    請求項９に記載の血管圧力差の取得装置。

Images