知財求人 - 知財ポータルサイト「IP Force」
▶ スーチョウ レインメド メディカル テクノロジー カンパニー リミテッドの特許一覧
特表2022-508239心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2022-01-19
(54)【発明の名称】心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法
(51)【国際特許分類】
A61B 6/03 20060101AFI20220112BHJP
【ＦＩ】
A61B6/03 370Z
A61B6/03 360D
【審査請求】有
【予備審査請求】未請求
(21)【出願番号】P 2021529764
(86)(22)【出願日】2019-01-10
(85)【翻訳文提出日】2021-05-25
(86)【国際出願番号】 CN2019071204
(87)【国際公開番号】W WO2020107668
(87)【国際公開日】2020-06-04
(31)【優先権主張番号】201811432014.5
(32)【優先日】2018-11-28
(33)【優先権主張国・地域又は機関】CN
(81)【指定国・地域】
(71)【出願人】
【識別番号】520283048
【氏名又は名称】スーチョウ レインメド メディカル テクノロジー カンパニー リミテッド
(74)【代理人】
【識別番号】110001896
【氏名又は名称】特許業務法人朝日奈特許事務所
(72)【発明者】
【氏名】フー、インフェイ
(72)【発明者】
【氏名】リュウ、クワンジー
(72)【発明者】
【氏名】ウー、シンユン
(72)【発明者】
【氏名】ワン、ツィーユアン
【テーマコード（参考）】
4C093
【Ｆターム（参考）】
4C093AA22
4C093AA26
4C093CA35
4C093CA37
4C093FD03
4C093FD07
4C093FD08
4C093FD11
4C093FF17
4C093FF19
4C093FF20
4C093FF22
4C093FF23
4C093FF28
4C093FF42
4C093FF45
4C093FF47
4C093FG04
4C093FH02
4C093FH06
4C093FH07
(57)【要約】
本発明は、心筋画像を抽出することにより、心筋体積を特定し、冠状動脈口を位置決めし、冠状動脈を正確に分割し、冠状動脈体データのエッジを検出することにより、計算に必要な格子モデルを生成し、非侵襲的測定により安静状態心筋血流量及びＣＦＲを特定し、最大充血状態における冠動脈入口箇所の総流量を計算し、冠動脈ツリーにおける異なる血管内面の最大充血状態の流量を特定し、そして最大充血状態の流速Ｖ1を特定し、最大充血状態における平均移動時間Ｔｍｎを得て、Ｖ1を冠動脈入口流速とし、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、微小循環抵抗指数を計算する、心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法を提供する。本発明は心臓ＣＴ画像に基づき、微小循環抵抗指数を早く、正確に、全自動で得ることができる。
【特許請求の範囲】
【請求項1】
Ｓ０１：心臓ＣＴ画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得て、心筋画像に基づいて心筋体積を特定し、Ｓ０２：大動脈画像を処理して全大動脈相補画像を得て、領域成長を行い、冠状動脈口を含む大動脈画像を得て、冠状動脈口を含む大動脈画像と全大動脈相補画像に基づき、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に基づいて冠状動脈口を特定し、
Ｓ０３：心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出し、
Ｓ０４：冠状動脈画像を二値化し、等値面画像を描き、冠状動脈三次元格子画像を得て、
Ｓ０５：最大充血状態における冠動脈入口箇所の総流量Ｑtotal＝心筋体積×心筋血流量×ＣＦＲを計算し、ＣＦＲは冠状動脈冠血流予備能（ｃｏｒｏｎａｒｙ ｆｌｏｗ ｒｅｓｅｒｖｅ：ＣＦＲ）であり、
Ｓ０６：充血状態における血流速度Ｖ1及び最大充血状態における平均移動時間Ｔｍｎを計算し、
Ｓ０７：Ｖ1を冠動脈狭窄血管の入口流速とし、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧であり、微小循環抵抗指数ＩＭＲ＝Ｐd×Ｔｍｎを得る、
ステップを含む、ことを特徴とする心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項2】
前記ステップＳ０２において、冠状動脈口を含む画像を得た後、冠状動脈口を含む画像に対して連結領域分析を行い、異なるグレースケールラベルで各連結領域をマークし、冠状動脈口を特定する、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項3】
前記ステップＳ０２において、心臓画像において、大動脈断面が円状になる特徴を利用し、上行大動脈及び中心線を抽出し、大動脈画像を得る、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項4】
前記ステップＳ０４において冠状動脈画像二値化は、冠状動脈画像Ｖ１中のボクセルをトラバーサルし、もしボクセル画素が０に等しい場合、当該画素値は変わらず、もし０に等しくない場合、画素値を１とし、一つの新しいデータＶ２を得ることを含む、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項5】
前記ステップＳ０５において心筋コントラストエコー法（ＭＣＥ）または単一光子放射断層撮影法（ＳＰＥＣＴ）または陽電子放出断層撮影法（ＰＥＴ）または心臓核磁気共鳴画像法（ＭＲＩ）またはＣＴ灌流により、安静状態心筋血流量及び冠状動脈冠血流予備能（ＣＦＲ）を特定する、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項6】
前記ステップＳ０６は、Ｓ６１：流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流量Ｑ＝Ｑtotal×（Ｖ／Ｖtotal）3/4を特定し、但し、Ｖtotalは心臓ＣＴ三次元再構成のすべての心臓表面冠状動脈の血液体積の和であり、Ｖは心臓表面冠状動脈ツリー内の任意の一本の血管及びその下流血管中の血液体積の和であり、
Ｓ６２：流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流速度Ｖ1＝Ｑ／Ｄを特定し、但し、Ｄは当該血管の平均直径であり、
Ｓ６３：冠動脈血管長さＬ及び流速Ｖ1に基づき最大充血状態における平均移動時間Ｔｍｎ＝Ｌ／Ｖ1を得ること
を含む、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項7】
前記ステップＳ０７は具体的に、血管三次元格子に対して解を求め、数値解法を用いて連続性及びナビエ－ストークス（Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ）方程式の解を求め：
【数1】
【数2】
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
三次元計算流体力学により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算する、但し、Ｐaは大動脈平均である、
ことを含む、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項8】
前記ステップＳ０７は、ＣＴ再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし、二次元軸対称モデルを構築し、二次元格子を分割し、数値解法を用いて連続性及びナビエ－ストークス方程式の解を求め：
【数3】
【数4】
【数5】
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
二次元流体力学により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算する、但し、Ｐaは大動脈平均圧である、ことを含む、ことを特徴とする請求項１に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【請求項9】
前記ステップＳ０７はさらに、血管の異なるタイプの弯曲に対して、三次元モデルで入口から出口までの圧力差を計算し、二次元軸対称モデルに照らして計算し、各類型の弯曲の二次元軸対称結果に対する修正係数を記憶するためのデータベースを構築し、
圧力を得た後、データベース中の修正係数に照らし、修正後の入口から出口までの圧力差を得て、その後ＩＭＲを計算する、
ことを含む、ことを特徴とする請求項８に記載の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法。
【発明の詳細な説明】
【技術分野】
【0001】
発明は冠状動脈画像学評価分野に関し、特に、心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法に関する。
【背景技術】
【0002】
冠動脈性心疾患は重大な心臓疾患であり、経皮的冠動脈インターベンション（ＰＣＩ）は既に冠動脈性心疾患治療のひとつの重要手段となっており、長年の研究で明らかなように、ＰＣＩ治療による冠動脈性心疾患、特に急性冠症候群（ａｃｕｔｅ ｃｏｒｏｎａｒｙ ｓｙｎｄｒｏｍｅ、ＡＣＳ）の病死率は著しく低下している。しかし既知の研究が示すように、仮に冠動脈インターベンション又は冠動脈バイパス手術の治療を受けても、２５％の患者が心臓表面血管再構成後、組織レベルの心筋再灌流を実現できず、その主な原因は冠動脈微小循環機能障害とされている。冠動脈微小循環機能の損傷は急性心筋梗塞予後を決定する独立的な危険要素であるため、患者の冠動脈微小循環機能状態に対する評価はますます重視されている。
【0003】
現在、冠動脈微小循環の評価は負荷心電図、放射性核種現像技術、核磁気共鳴画像法技術等の方法により行うことができる。しかし以上の各評価技術は検査技術手段自体及び検査員の影響を大きく受けるため、冠動脈微小循環機能を正確に評価することができないとともに、これら指標が反映するのは心臓表面血管及び微小循環の共同作用の結果である。心臓表面血管の影響ではなく、冠動脈微小循環機能状態のみを反映する、より正確で、より安定した指標が必要である。現在のＩＭＲはおそらく上記条件を満たす一つの理想指標である。ＩＭＲは冠状動脈微小血管抵抗を反映する指標であり、冠動脈遠位動脈圧力（Ｐｄ）を最大充血状態における平均移動時間（Ｔｍｎ）の逆数で割ると定義され、言い換えれば、ＰｄとＴｍｎの積であり、単位はｍｍＨｇ・ｓである。ＩＭＲ＜２５は正常、＞３０は異常とされ、微小循環抵抗の上昇を知らせ、２５－３０はグレーゾーンである。
【0004】
従来のＩＭＲ測定方法は０．０１４インチの軟圧ガイドワイヤにより冠動脈圧力及び温度を同時記録し、ガイドワイヤスティック上の2つの温度センサーは温度変化の時間差を探知すると、塩水がガイディングカテーテルからガイドワイヤ先端の温度センサーまで移動した平均移動時間（ｔｒａｎｓｉｔ ｍｅａｎ ｔｉｍｅ、Ｔｍｎ）を知ることができ、ＰdとＴｍｎの積に基づきＩＭＲ値を得ることができる。しかし上記圧力ガイドワイによるＩＭＲ測定は血管末端へのインターベンションを必要とし、手術の難度及びリスクが増大し、また、圧力ガイドワイヤの高い価格もその大規模応用を制限している。
【0005】
冠動脈ＣＴＡは冠動脈狭窄程度を正確に評価でき、且つ血管壁プラーク性質を見分けることができる、一つの非侵襲的で、操作が簡単な冠状動脈病変診断検査方法であり、リスクの高い人をスクリーニングする優先的な方法となる。よって、冠動脈性心疾患患者の血管に対して介入する場合、初期段階で患者冠動脈に対してＣＴＡ評価を行わなければならない。
【0006】
冠動脈ＣＴＡにより非侵襲で得られた微小循環抵抗指数を計算する場合、別途、イメージング検査又は薬は不要であり、根本的に、不必要な冠動脈血管造影と血行再建治療を避けることができる。しかしＣＴＡは充血状態における冠動脈流速を測定できず、直接、微小循環抵抗指数を計算することができない。
【発明の概要】
【0007】
上記技術課題を解決するために、本発明は、非侵襲的測定により安静状態心筋血流量及び冠状動脈冠血流予備能（ＣＦＲ）を特定し、そして冠動脈ツリーにおける異なる血管内面の最大充血状態の流量を特定し、そして最大充血状態の流速Ｖ1を特定することによって、微小循環抵抗指数を早く、正確に、全自動で得ることができる、造影画像及び流体力学モデルにおける微小循環抵抗指数の計算方法を提供することを目的とする。
【0008】
本発明の技術方案は以下のとおり：
心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法は以下のステップを含む：
Ｓ０１：心臓ＣＴ画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得て、心筋画像に基づいて心筋体積を特定し、
Ｓ０２：大動脈画像を処理して全大動脈相補画像を得て、領域成長を行い、冠状動脈口を含む大動脈画像を得て、冠状動脈口を含む大動脈画像と全大動脈相補画像に基づき、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に基づいて冠状動脈口を特定し、
Ｓ０３：心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出し、
Ｓ０４：冠状動脈画像を二値化し、等値面画像を描き、冠状動脈三次元格子画像を得て、
Ｓ０５：最大充血状態における冠動脈入口箇所の総流量Ｑtotal＝心筋体積×心筋血流量×ＣＦＲを計算し、ＣＦＲは冠状動脈冠血流予備能（ｃｏｒｏｎａｒｙ ｆｌｏｗ ｒｅｓｅｒｖｅ：ＣＦＲ）であり、
Ｓ０６：充血状態における血流速度Ｖ1及び最大充血状態における平均移動時間Ｔｍｎを計算し、
Ｓ０７：Ｖ1を冠動脈狭窄血管の入口流速とし、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧であり、微小循環抵抗指数ＩＭＲ＝Ｐd×Ｔｍｎを得る。
心筋血流量及びＣＴ画像に基づく微小循環抵抗指数の計算方法は以下のステップを含む：
Ｓ０１：心臓ＣＴ画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得て、心筋画像に基づいて心筋体積を特定し、
Ｓ０２：大動脈画像を処理して全大動脈相補画像を得て、領域成長を行い、冠状動脈口を含む大動脈画像を得て、冠状動脈口を含む大動脈画像と全大動脈相補画像に基づき、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に基づいて冠状動脈口を特定し、
Ｓ０３：心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出し、
Ｓ０４：冠状動脈画像を二値化し、等値面画像を描き、冠状動脈三次元格子画像を得て、
Ｓ０５：最大充血状態における冠動脈入口箇所の総流量Ｑtotal＝心筋体積×心筋血流量×ＣＦＲを計算し、ＣＦＲは冠状動脈冠血流予備能（ｃｏｒｏｎａｒｙ ｆｌｏｗ ｒｅｓｅｒｖｅ：ＣＦＲ）であり、
Ｓ０６：充血状態における血流速度Ｖ1及び最大充血状態における平均移動時間Ｔｍｎを計算し、
Ｓ０７：Ｖ1を冠動脈狭窄血管の入口流速とし、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧であり、微小循環抵抗指数ＩＭＲ＝Ｐd×Ｔｍｎを得る。
【0009】
好ましい技術方案において、前記ステップＳ０２において、冠状動脈口を含む画像を得た後、冠状動脈口を含む画像に対して連結領域分析を行い、異なるグレースケールラベルで各連結領域をマークし、冠状動脈口を特定する。
【0010】
好ましい技術方案において、前記ステップＳ０２において、心臓画像において、大動脈断面が円状になる特徴を利用し、上行大動脈及び中心線を抽出し、大動脈画像を得る。
【0011】
好ましい技術方案において、前記ステップＳ０４において冠状動脈画像二値化は、冠状動脈画像Ｖ１中のボクセルをトラバーサルし、もしボクセル画素Ａ１が０に等しい場合、当該画素値は変わらず、もしＡ１が０に等しくない場合、Ａ１の画素値を１とし、一つの新しいデータＶ２を得ることを含む。
【0012】
好ましい技術方案において、前記ステップＳ０５において心筋コントラストエコー法（ＭＣＥ）または単一光子放射断層撮影法（ＳＰＥＣＴ）または陽電子放出断層撮影法（ＰＥＴ）または心臓核磁気共鳴画像法（ＭＲＩ）またはＣＴ灌流により、安静状態心筋血流量及び冠状動脈冠血流予備能（ＣＦＲ）を特定する。
【0013】
好ましい技術方案において、前記ステップＳ０６は、
Ｓ６１：流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流量Ｑ＝Ｑtotal×（Ｖ／Ｖtotal）3/4を特定し、但し、Ｖtotalは心臓ＣＴ三次元再構成のすべての心臓表面冠状動脈の血液体積の和であり、Ｖは心臓表面冠状動脈ツリー内の任意の一本の血管及びその下流血管中の血液体積の和であり、
Ｓ６２：流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流速度Ｖ1＝Ｑ／Ｄを特定し、但し、Ｄは当該血管の平均直径であり、
Ｓ６３：冠動脈血管長さＬ及び流速Ｖ1に基づき最大充血状態における平均移動時間Ｔｍｎ＝Ｌ／Ｖ1を得ること、を含む。
Ｓ６１：流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流量Ｑ＝Ｑtotal×（Ｖ／Ｖtotal）3/4を特定し、但し、Ｖtotalは心臓ＣＴ三次元再構成のすべての心臓表面冠状動脈の血液体積の和であり、Ｖは心臓表面冠状動脈ツリー内の任意の一本の血管及びその下流血管中の血液体積の和であり、
Ｓ６２：流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流速度Ｖ1＝Ｑ／Ｄを特定し、但し、Ｄは当該血管の平均直径であり、
Ｓ６３：冠動脈血管長さＬ及び流速Ｖ1に基づき最大充血状態における平均移動時間Ｔｍｎ＝Ｌ／Ｖ1を得ること、を含む。
【0014】
好ましい技術方案において、前記ステップＳ０７は、血管三次元格子に対して解を求め、数値解法を用いて連続性及びナビエ－ストークス（Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ）方程式の解を求め：
但し、
Ｐ、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
三次元計算流体力学により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算する、但し、Ｐaは大動脈平均である、ことを含む。
【数1】
【数2】
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
三次元計算流体力学により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算する、但し、Ｐaは大動脈平均である、ことを含む。
【0015】
好ましい技術方案において、前記ステップＳ０７は、
ＣＴ再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし、二次元軸対称モデルを構築し、二次元格子を分割し、数値解法を用いて連続性及びナビエ－ストークス方程式の解を求め：
但し、ρは血液の密度を表し、ｕz、ｕrはそれぞれｚ方向、ｒ方向の流速を表し、μは血液の動粘度を表し、ｐは血液の圧の強さを表し、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
二次元流体力学により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算する、但し、Ｐaは大動脈平均圧である、ことを含む。
ＣＴ再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし、二次元軸対称モデルを構築し、二次元格子を分割し、数値解法を用いて連続性及びナビエ－ストークス方程式の解を求め：
【数3】
【数4】
【数5】
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
二次元流体力学により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算する、但し、Ｐaは大動脈平均圧である、ことを含む。
【0016】
好ましい技術方案において、前記ステップＳ０７はさらに、血管の異なるタイプの弯曲に対して、三次元モデルで入口から出口までの圧力差を計算し、二次元軸対称モデルに照らして計算し、各類型の弯曲の二次元軸対称結果に対する修正係数を記憶するためのデータベースを構築し、圧力を得た後、データベース中の修正係数に照らし、修正後の入口から出口までの圧力差を得て、その後ＩＭＲを計算する。
【0017】
従来技術に比べ、本発明は以下の長所を有する：
本発明は心筋血流量及び心臓ＣＴ画像により、微小循環抵抗指数を早く、正確に、全自動で得ることができる。非侵襲的測定により、操作が簡単で、手術の難度及びリスクを大幅に低減し、臨床で大規模に応用を進めることができる。
本発明は心筋血流量及び心臓ＣＴ画像により、微小循環抵抗指数を早く、正確に、全自動で得ることができる。非侵襲的測定により、操作が簡単で、手術の難度及びリスクを大幅に低減し、臨床で大規模に応用を進めることができる。
【図面の簡単な説明】
【0018】
以下に図面及び実施例と合わせて本発明を更に説明する。
【0019】
【発明を実施するための形態】
【0020】
本発明の目的、技術方案及び長所をより明確に説明するために、以下に図面を参照しながら発明を実施するための形態と合わせて、本発明について、さらに詳細に説明する。これら説明は単なる例であり、本発明の範囲を限定するものではないと理解されるべきである。この他、以下の説明において、本発明の概念と不要な混同を避けるために、公知の構造及び技術に関する説明は省略する。
【0021】
所定の心臓ＣＴ画像について、リバースメソッドにより、心臓を抽出し、非目標領域である下行大動脈、脊椎、肋骨を対象に処理を行い、一つ一つ胸壁、肺部、椎骨及び下行大動脈等の非心臓組織を取り除いて心臓画像を抽出して得る。得られた心臓画像において、大動脈断面が円状になる特徴を利用し、上行大動脈及び中心線を抽出し、大動脈画像を得る。
【0022】
図１に示されるように、本発明の心筋血流量及びＣＴ画像に基づく微小循環抵抗指数ＩＭＲの計算方法は、心筋画像を抽出し、冠状動脈口を抽出し、冠状動脈を抽出し、冠状動脈格子モデルを生成し、安静状態心筋血流量及び冠状動脈冠血流予備能（ＣＦＲ）を特定し、最大充血状態における冠動脈入口箇所の総流量を計算し、充血状態における血流速度Ｖ1を計算し、微小循環抵抗指数を特定することを含む。具体的に以下のステップを含む：
【0023】
１：心筋画像を抽出する：
心臓ＣＴ画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得る。図２参照。
心臓ＣＴ画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得る。図２参照。
【0024】
２：冠状動脈口を抽出する：
大動脈画像の二値化画像に対して形態学膨張を行い、全大動脈の二値画像を得て、そして画素リバースにより全大動脈相補画像を得る。
大動脈中心線上点の平均グレースケールに基づき領域成長を行い、冠状動脈口を含む大動脈画像を得る。図３参照。
冠状動脈口を含む大動脈画像と全大動脈相補画像により画像乗算を行い、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に対して連結領域分析を行い、異なるグレースケールラベルで各連結領域をマークし、冠状動脈口を特定する。図４参照。
大動脈画像の二値化画像に対して形態学膨張を行い、全大動脈の二値画像を得て、そして画素リバースにより全大動脈相補画像を得る。
大動脈中心線上点の平均グレースケールに基づき領域成長を行い、冠状動脈口を含む大動脈画像を得る。図３参照。
冠状動脈口を含む大動脈画像と全大動脈相補画像により画像乗算を行い、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に対して連結領域分析を行い、異なるグレースケールラベルで各連結領域をマークし、冠状動脈口を特定する。図４参照。
【0025】
３：冠状動脈を抽出する：
心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出する。図５参照。
心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出する。図５参照。
【0026】
４：冠状動脈格子モデルを生成する：
ステップ三により、冠状動脈画像データＶ１を得て、当該データ中のボクセルは空間上に一つの立方体を構成し、冠状動脈部分に属すボクセル画素値は０ではなく（画素値約－３０００～３０００の間）、残りのボクセル画素値はすべて０である。
本ステップにおいて、ステップ五における計算のために、データを空間三次元格子データＶ３に変える必要がある。
ステップ三により、冠状動脈画像データＶ１を得て、当該データ中のボクセルは空間上に一つの立方体を構成し、冠状動脈部分に属すボクセル画素値は０ではなく（画素値約－３０００～３０００の間）、残りのボクセル画素値はすべて０である。
本ステップにおいて、ステップ五における計算のために、データを空間三次元格子データＶ３に変える必要がある。
【0027】
（１）冠状動脈データ二値化
冠状動脈画像データＶ１中のボクセルをトラバーサルし、簡単な画素値判断を行い、もし画素Ａ１が０に等しい場合、当該画素値は変わらず、もしＡ１が０に等しくない場合、Ａ１の画素値を１とする。
最終的に一つの新しい画像データＶ２が得られ、当該画像において、冠状動脈部分に属するボクセル画素値は１であり、残りの部分は０である。
冠状動脈画像データＶ１中のボクセルをトラバーサルし、簡単な画素値判断を行い、もし画素Ａ１が０に等しい場合、当該画素値は変わらず、もしＡ１が０に等しくない場合、Ａ１の画素値を１とする。
最終的に一つの新しい画像データＶ２が得られ、当該画像において、冠状動脈部分に属するボクセル画素値は１であり、残りの部分は０である。
【0028】
（２）等値面生成
ボクセルは一つの極小の六面体であると定義され、隣接する上下層の間の四つの画素が立方体上の八つの頂点を構成する。等値面とは空間中のある同じ属性値を有するすべての点の集合である。これは次のように表すことができる：
｛（ｘ、ｙ、ｚ）│ｆ（ｘ、ｙ、ｚ）＝ｃ｝、ｃは定数である
本方法におけるｃは三次元再構成過程において与えられた画素値１である。
ボクセルは一つの極小の六面体であると定義され、隣接する上下層の間の四つの画素が立方体上の八つの頂点を構成する。等値面とは空間中のある同じ属性値を有するすべての点の集合である。これは次のように表すことができる：
｛（ｘ、ｙ、ｚ）│ｆ（ｘ、ｙ、ｚ）＝ｃ｝、ｃは定数である
本方法におけるｃは三次元再構成過程において与えられた画素値１である。
【0029】
等値面を抽出するフローは以下のとおり：
（１）原始データを前処理した後、特定の配列に読み込む。
（２）格子データボディから一つのユニットボディを抽出して現在のユニットボディとすると同時に当該ユニットボディのすべての情報を得る。
（３）現在のユニットボディの８個の頂点の函数値と所定の等値面値Ｃを比較し、当該ユニットボディの状態表を得る。
（４）現在のユニットボディの状態表インデックスから、等値面と交わるユニットボディエッジを探し出し、そして線形補間の方法により各交点の位置座標を計算する。
（５）中心差分法により現在のユニットボディの８個の頂点の法線ベクトルを求め、それから線形補間の方法により三角パッチの各頂点の法線を得る。
（６）各三角パッチの頂点座標及び頂点の法線ベクトルに基づき等値面イメージの作成を行う。
最終的に冠状動脈の三次元格子画像データＶ３を得る。図６参照。
（１）原始データを前処理した後、特定の配列に読み込む。
（２）格子データボディから一つのユニットボディを抽出して現在のユニットボディとすると同時に当該ユニットボディのすべての情報を得る。
（３）現在のユニットボディの８個の頂点の函数値と所定の等値面値Ｃを比較し、当該ユニットボディの状態表を得る。
（４）現在のユニットボディの状態表インデックスから、等値面と交わるユニットボディエッジを探し出し、そして線形補間の方法により各交点の位置座標を計算する。
（５）中心差分法により現在のユニットボディの８個の頂点の法線ベクトルを求め、それから線形補間の方法により三角パッチの各頂点の法線を得る。
（６）各三角パッチの頂点座標及び頂点の法線ベクトルに基づき等値面イメージの作成を行う。
最終的に冠状動脈の三次元格子画像データＶ３を得る。図６参照。
【0030】
５：充血状態における血流速度Ｖ1を計算する：
心筋コントラストエコー法（ＭＣＥ）または単一光子放射断層撮影法（ＳＰＥＣＴ）または陽電子放出断層撮影法（ＰＥＴ）または心臓核磁気共鳴画像法（ＭＲＩ）またはＣＴ灌流等非侵襲的測定により、安静状態心筋血流量及び冠状動脈冠血流予備能（ＣＦＲ）を特定し、心筋体積、心筋血流量、ＣＦＲにより、最大充血状態における冠動脈入口箇所（左冠動脈ツリー及び右冠動脈ツリーの和を含む）の総流量Ｑtotal＝心筋体積×心筋血流量×ＣＦＲを計算する。
心筋コントラストエコー法（ＭＣＥ）または単一光子放射断層撮影法（ＳＰＥＣＴ）または陽電子放出断層撮影法（ＰＥＴ）または心臓核磁気共鳴画像法（ＭＲＩ）またはＣＴ灌流等非侵襲的測定により、安静状態心筋血流量及び冠状動脈冠血流予備能（ＣＦＲ）を特定し、心筋体積、心筋血流量、ＣＦＲにより、最大充血状態における冠動脈入口箇所（左冠動脈ツリー及び右冠動脈ツリーの和を含む）の総流量Ｑtotal＝心筋体積×心筋血流量×ＣＦＲを計算する。
【0031】
流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流量Ｑ：Ｑ＝Ｑtotal×（Ｖ／Ｖtotal）3/4を特定し、但し、Ｖtotalは心臓ＣＴ三次元再構成のすべての心臓表面冠状動脈（左冠動脈ツリー及び右冠動脈ツリーの和を含む）の血液体積の和であり、Ｖは心臓表面冠状動脈ツリー内の任意の一本の血管及びその下流血管中の血液体積の和である。図７参照。流量体積スケーリング則及び心臓ＣＴ三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流速度Ｖ1：Ｖ1＝Ｑ／Ｄを特定し、但し、Ｄは当該血管の平均直径である（当該血管の血液体積を当該血管の長さで割る）。冠動脈血管長さＬ及び流速Ｖ1に基づき、式Ｔｍｎ＝Ｌ／Ｖ1により、最大充血状態における平均移動時間Ｔｍｎを計算する。
【0032】
６：微小循環抵抗指数の計算：
Ｖ1を冠動脈狭窄血管の入口流速とし、計算流体力学（ＣＦＤ）方法により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧であり、最後に式ＩＭＲ＝Ｐd×Ｔｍｎにより微小循環抵抗指数を計算する。
Ｖ1を冠動脈狭窄血管の入口流速とし、計算流体力学（ＣＦＤ）方法により各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧であり、最後に式ＩＭＲ＝Ｐd×Ｔｍｎにより微小循環抵抗指数を計算する。
【0033】
三次元モデルに対する処理ステップは以下を含む：
ＣＴ再構成のジオメトリ構造に基づき、三次元格子を分割し、数値解法（例えば、有限差分、有限要素法、有限体積法等）により連続性及びナビエ－ストークス（Ｎａｖｉｅｒ‐Ｓｔｏｋｅｓ）方程式の解を求め：
但し、
Ｐ、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
式［Ａ１］及び［Ａ２］に基づき、三次元ＣＦＤを実行して各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧である。
ＣＴ再構成のジオメトリ構造に基づき、三次元格子を分割し、数値解法（例えば、有限差分、有限要素法、有限体積法等）により連続性及びナビエ－ストークス（Ｎａｖｉｅｒ‐Ｓｔｏｋｅｓ）方程式の解を求め：
【数6】
【数7】
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
式［Ａ１］及び［Ａ２］に基づき、三次元ＣＦＤを実行して各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧である。
【0034】
二次元モデルに対して、以下のステップを含む：
ＣＴ再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし（二次元軸対称モデル）、二次元格子を分割し、数値解法（例えば、有限差分、有限要素法、有限体積法等）により連続性及びナビエ－ストークス方程式の解を求め：
但し、ρは血液の密度を表し、ｕz、ｕrはそれぞれｚ方向、ｒ方向の流速を表し、μは血液の動粘度を表し、ｐは血液の圧の強さを表し、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
式［Ａ３］-［Ａ５］に基づき、二次元ＣＦＤを実行し、各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧である。
ＣＴ再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし（二次元軸対称モデル）、二次元格子を分割し、数値解法（例えば、有限差分、有限要素法、有限体積法等）により連続性及びナビエ－ストークス方程式の解を求め：
【数8】
【数9】
【数10】
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ1であり、
式［Ａ３］-［Ａ５］に基づき、二次元ＣＦＤを実行し、各冠動脈狭窄の圧力損失ΔＰ1、ΔＰ2、ΔＰ3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰi（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐd＝Ｐa－ΔＰを計算し、但し、Ｐaは大動脈平均圧である。
【0035】
血管の異なるタイプの弯曲に対して、三次元モデルで入口から出口までの圧力差を計算し、二次元軸対称モデルに照らして計算し、各類型の弯曲の二次元軸対称結果に対する修正係数を記憶するためのデータベースを構築し、圧力算出後データベース中の修正係数に照らし、入口から出口までの圧力差を正確に得て、最後に式ＩＭＲ＝Ｐd×Ｔｍｎにより微小循環抵抗指数を計算する。
【0036】
本発明の上記発明を実施するための形態は本発明の原理を単に例示的に説明又は解釈するためのものであり、本発明を限定するものではないと理解されるべきである。よって、本発明の精神及び範囲を逸脱しない限り、如何なる修正、均等の差し替え、改良等もすべて本発明の保護範囲に含まれる。この他、本発明に付属の請求項は、付属の請求項の範囲及び境界、またはこれらの範囲及び境界の均等形式における全ての変更及び修正例を含むことを意図するものである。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【国際調査報告】