知財求人 - 知財ポータルサイト「IP Force」
▶ ▲蘇▼州▲潤▼▲邁▼▲徳▼医▲療▼科技有限公司の特許一覧
特許7437077生理パラメータに基づき血管評定パラメータを取得する方法、装置及び記憶媒体
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-02-14
(45)【発行日】2024-02-22
(54)【発明の名称】生理パラメータに基づき血管評定パラメータを取得する方法、装置及び記憶媒体
(51)【国際特許分類】
A61B 5/0275 20060101AFI20240215BHJP
【ＦＩ】
A61B5/0275
【請求項の数】
15
(21)【出願番号】P 2022525373
(86)(22)【出願日】2019-11-08
(65)【公表番号】
(43)【公表日】2022-12-28
(86)【国際出願番号】 CN2019116664
(87)【国際公開番号】W WO2021087967
(87)【国際公開日】2021-05-14
【審査請求日】2022-04-28
(31)【優先権主張番号】201911066155.4
(32)【優先日】2019-11-04
(33)【優先権主張国・地域又は機関】CN
(73)【特許権者】
【識別番号】522036093
【氏名又は名称】▲蘇▼州▲潤▼▲邁▼▲徳▼医▲療▼科技有限公司
(74)【代理人】
【識別番号】110001896
【氏名又は名称】弁理士法人朝日奈特許事務所
(72)【発明者】
【氏名】▲劉▼ ▲廣▼志
(72)【発明者】
【氏名】▲きょう▼ ▲艷▼君
(72)【発明者】
【氏名】李 建平
(72)【発明者】
【氏名】易 ▲鉄▼慈
(72)【発明者】
【氏名】▲鄭▼ 博
【審査官】高松 大
(56)【参考文献】
【文献】中国特許出願公開第１０８５１１０７５（ＣＮ，Ａ）
【文献】中国特許出願公開第１０３３００８２０（ＣＮ，Ａ）
【文献】中国特許出願公開第１０８５５０１８９（ＣＮ，Ａ）
【文献】中国特許出願公開第１０８７４２５８７（ＣＮ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／０２７５
(57)【特許請求の範囲】
【請求項1】
生理パラメータ取得ユニットによって生理パラメータである糖尿病特徴値、高血圧特徴値、及び性別を取得するステップと、血流速度取得ユニットによって血流速度vを取得するステップと、
大動脈圧波形取得ユニットによって時間に伴い変わる大動脈圧波形をリアルタイムに取得するステップと、
冠状動脈血管評定パラメータユニットによって、前記血流速度v、前記大動脈圧波形、前記生理パラメータに基づき、冠状動脈血管評定パラメータを取得するステップであって、前記冠状動脈血管評定パラメータは拡張期の微小血管抵抗指数CAIFMRである、冠状動脈血管評定パラメータを取得するステップと、
前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取るステップと、
前記vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップと、
を含む、生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法が、記憶されたコンピュータ・プログラムがプロセッサによって実行されることにより、実施されることを特徴とするコンピュータ記憶媒体であって、
CAIFMRは、
但し、
ことを特徴とするコンピュータ記憶媒体。
【請求項2】
前記大動脈拡張期の平均圧を取得するステップにおいて、
但し、
ことを特徴とする請求項１に記載のコンピュータ記憶媒体。
【請求項3】
患者に糖尿病がない場合、0.5≦a≦1とされ、患者に糖尿病がある場合、1＜a≦2とされ、患者の血圧値が90mmHg以上である場合、1＜b≦1.5とされ、患者の血圧値が90mmHg未満である場合、0.5≦b≦1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=3～10とされる、
ことを特徴とする請求項2に記載のコンピュータ記憶媒体。
【請求項4】
患者に糖尿病がない場合、a=1とされ、患者に糖尿病がある場合、a=2とされ、患者の血圧値が90mmHg以上である場合、b=1.5とされ、患者の血圧値が90mmHg未満である場合、b=1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=5とされる、
ことを特徴とする請求項3に記載のコンピュータ記憶媒体。
【請求項5】
前記血流速度を取得するステップは、画像読み取りユニットによって少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップと、
血管段抽出ユニットによって前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップと、
中心線抽出ユニットによって前記血管段の中心線を抽出するステップと、
時間差ユニットによって任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記血管段の中心線の差を計算し、その差の値をΔLとするステップと、
血流速度取得ユニットによって前記ΔLと前記Δtの比に基づき、血流速度を求めるステップと、
を含む、
ことを特徴とする請求項2に記載のコンピュータ記憶媒体。
【請求項6】
前記の前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップは、前記冠状動脈二次元造影画像組からNフレーム冠状動脈二次元造影画像を選び取るステップと、
前記冠状動脈二次元造影画像において、注目する前記血管の最初と最後の点をピックアップし、注目する前記血管段を取得するステップと、
を含む、
ことを特徴とする請求項5に記載のコンピュータ記憶媒体。
【請求項7】
前記の前記血管段の中心線を抽出するステップは、前記冠状動脈二次元造影画像から血管骨格を抽出するステップと、
前記血管段の延在方向、および2点間の最短経路を取る原則に基づくステップと、
前記血管骨格に沿って、前記血管段の中心線を抽出するステップと、
を含む、
ことを特徴とする請求項5に記載のコンピュータ記憶媒体。
【請求項8】
前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記血管段の中心線の差を計算し、その差の値をΔLとし、且つ、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、造影剤が冠状動脈入口である前記血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が前記血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップと、
第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップと、
v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度をそれぞれv1,..., vb,..., va,..., vN-1とするステップと、
を含む、
ことを特徴とする請求項6に記載のコンピュータ記憶媒体。
【請求項9】
前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記血管段の中心線の差を計算し、その差の値をΔLとし、且つ、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1から第1フレームの画像までの時間差及び中心線の長さの差を順に求めるステップと、
v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1から第1フレーム画像までの血流速度をそれぞれ得るステップと、
を含む、
ことを特徴とする請求項6に記載のコンピュータ記憶媒体。
【請求項10】
前記生理的パラメータに基づいて冠動脈血管評定パラメータを取得する方法が、
前記の前記血管段の中心線を抽出するステップの後、前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記血管段の中心線の差を計算し、その差の値をΔLとするステップの前にさらに、
少なくとも2つの体位の冠状動脈二次元造影画像組を読み取るステップと、
前記血管段のジオメトリ構造情報を取得するステップと、
注目する前記血管段に対して図形処理を行うステップと、
前記血管段の血管輪郭線を抽出するステップと、
前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するステップと、
を含む、
ことを特徴とする請求項5に記載のコンピュータ記憶媒体。
【請求項11】
前記の前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、前記三次元血管モデルに基づき、三次元血管モデルの中心線を取得し、前記冠状動脈二次元造影画像により抽出した中心線に対して補正を行い、且つ、中心線の差の値ΔLに対して補正を行い、ΔL'を得るステップと、
前記ΔL'と前記Δtの比に基づき、血流速度vを求めるステップと、
を含む、
ことを特徴とする請求項10に記載のコンピュータ記憶媒体。
【請求項12】
血流速度取得ユニット、大動脈圧波形取得ユニット、及び冠状動脈血管評定パラメータユニットを備え、前記冠状動脈血管評定パラメータユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニットに接続され、
前記血流速度取得ユニットは血流速度vを取得するためのものであり、
前記大動脈圧波形取得ユニットは時間に伴い変わる大動脈圧波形をリアルタイムに取得するためのものであり、
前記冠状動脈血管評定パラメータユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニットが送信した前記血流速度v及び大動脈圧波形を受信し、それから生理パラメータである糖尿病特徴値、高血圧特徴値、及び性別に基づき、冠状動脈血管評定パラメータを取得するためのものであり、前記冠状動脈血管評定パラメータは拡張期の微小血管抵抗指数CAIFMRであり、前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取り、前記vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得し、前記最大血流速度vmax及び前記大動脈拡張期の平均圧に基づき、前記拡張期の微小血管抵抗指数CAIFMRを取得し、
CAIFMRは、
但し、
ことを特徴とする請求項1～11のいずれかに記載のコンピュータ記憶媒体に記憶されたコンピュータ・プログラムがプロセッサによって実行されると、実施される生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法に用いられる生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
【請求項13】
請求項12に記載の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置を含むことを特徴とする冠状動脈分析システム。
【請求項14】
１つまたは複数のプロセッサと、請求項１～１１のいずれか１項に記載のコンピュータ記憶媒体を含むことを特徴とする冠状動脈分析システム。
【請求項15】
１つまたは複数のプロセッサと、請求項１～１１のいずれか１項に記載のコンピュータ記憶媒体を含むことを特徴とする電子装置。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は冠状動脈技術分野に関し、特に、生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法、装置、冠状動脈分析システム及びコンピュータ記憶媒体に関する。
【背景技術】
【0002】
世界保健機関の統計によれば、心血管疾患は既に人の健康を奪う「トップキラー」となっている。近年、血流動力学を用いて心血管疾患の生理及び病理行為を分析することは既に心血管疾患診断の一つの非常に重要な手段となっている。
血液流量及び流速は血流動力学の非常に重要なパラメータとされており、どのように血液流量及び流速を正確に簡単に測定するかが多くの研究学者の研究の重点となっている。
冠状動脈血管評定パラメータとしては、冠状動脈拡張期の血流速度CAIFR、および冠状動脈拡張期の微小血管抵抗指数CAIFMR等を含む。しかし個体群により生命兆候は異なるため、正常な数値の評定標準にも若干違いがある。例えば、高齢者の心筋微小循環機能は低く、血流速度は一般的に若者よりも低いため、もし業界の汎用評価標準を用いる場合、採用する血流速度は実際の値より大きめになり、CAIFMR等が過小評価され、CAIFMR等の測定の正確性が低下する。
血液流量及び流速は血流動力学の非常に重要なパラメータとされており、どのように血液流量及び流速を正確に簡単に測定するかが多くの研究学者の研究の重点となっている。
冠状動脈血管評定パラメータとしては、冠状動脈拡張期の血流速度CAIFR、および冠状動脈拡張期の微小血管抵抗指数CAIFMR等を含む。しかし個体群により生命兆候は異なるため、正常な数値の評定標準にも若干違いがある。例えば、高齢者の心筋微小循環機能は低く、血流速度は一般的に若者よりも低いため、もし業界の汎用評価標準を用いる場合、採用する血流速度は実際の値より大きめになり、CAIFMR等が過小評価され、CAIFMR等の測定の正確性が低下する。
【発明の概要】
【0003】
本発明は、個体群により生命兆候が異なるため、正常な数値の評定標準にも若干違いが出て、冠状動脈血管評定パラメータ測定が不正確になる問題を解決するために、生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法、装置、冠状動脈分析システム及びコンピュータ記憶媒体を提供する。
【0004】
上記目的を実現するために、第一の態様において、本出願は以下の方法を提供する。
【0005】
生理パラメータを取得するステップと、
血流速度vを取得するステップと、
時間に伴い変わる大動脈圧波形をリアルタイムに取得するステップと、
前記血流速度v、前記大動脈圧波形、前記生理パラメータに基づき、冠状動脈血管評定パラメータを取得するステップと、
を含む、
生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
血流速度vを取得するステップと、
時間に伴い変わる大動脈圧波形をリアルタイムに取得するステップと、
前記血流速度v、前記大動脈圧波形、前記生理パラメータに基づき、冠状動脈血管評定パラメータを取得するステップと、
を含む、
生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0006】
任意選択的に、
前記冠状動脈評定パラメータは拡張期の微小血管抵抗指数CAIFMRであり、
前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取るステップと、
前記vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップと、
を含み、
CAIFMRは、
であり、
は
であり、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vmaxは拡張期の最大血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、k=a×b、aは糖尿病特徴値を表し、bは高血圧特徴値を表し、cは性別を表す、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記冠状動脈評定パラメータは拡張期の微小血管抵抗指数CAIFMRであり、
前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取るステップと、
前記vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップと、
を含み、
CAIFMRは、
但し、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0007】
任意選択的に、
患者に糖尿病がない場合、0.5≦a≦1とされ、患者に糖尿病がある場合、1＜a≦2とされ、
患者の血圧値が90mmHg以上である場合、1＜b≦1.5とされ、患者の血圧値が90mmHg未満である場合、0.5≦b≦1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=3～10とされる、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
患者に糖尿病がない場合、0.5≦a≦1とされ、患者に糖尿病がある場合、1＜a≦2とされ、
患者の血圧値が90mmHg以上である場合、1＜b≦1.5とされ、患者の血圧値が90mmHg未満である場合、0.5≦b≦1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=3～10とされる、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0008】
任意選択的に、
患者に糖尿病がない場合、a=1とされ、患者に糖尿病がある場合、a=2とされ、
患者の血圧値が90mmHg以上である場合、b=1.5とされ、患者の血圧値が90mmHg未満である場合、b=1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=5とされる、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
患者に糖尿病がない場合、a=1とされ、患者に糖尿病がある場合、a=2とされ、
患者の血圧値が90mmHg以上である場合、b=1.5とされ、患者の血圧値が90mmHg未満である場合、b=1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=5とされる、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0009】
任意選択的に、
前記血流速度を取得するステップは、
少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップと、
前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップと、
前記血管段の中心線を抽出するステップと、
任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとするステップと、
前記ΔLと前記Δtの比に基づき、血流速度を求めるステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記血流速度を取得するステップは、
少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップと、
前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップと、
前記血管段の中心線を抽出するステップと、
任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとするステップと、
前記ΔLと前記Δtの比に基づき、血流速度を求めるステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0010】
任意選択的に、
前記少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップは、
無線または有線方式により造影画像撮影装置または病院プラットフォームから、少なくとも一つの体位の冠状動脈二次元造影画像組を直接読み取るステップ、又は
記憶装置により、少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップは、
無線または有線方式により造影画像撮影装置または病院プラットフォームから、少なくとも一つの体位の冠状動脈二次元造影画像組を直接読み取るステップ、又は
記憶装置により、少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0011】
任意選択的に、
前記の前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップは、
前記冠状動脈二次元造影画像組からNフレーム冠状動脈二次元造影画像を選び取るステップと、
前記冠状動脈二次元造影画像において、注目する前記血管の最初と最後の点をピックアップし、注目する前記血管段を取得するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップは、
前記冠状動脈二次元造影画像組からNフレーム冠状動脈二次元造影画像を選び取るステップと、
前記冠状動脈二次元造影画像において、注目する前記血管の最初と最後の点をピックアップし、注目する前記血管段を取得するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0012】
任意選択的に、
前記の前記血管段の中心線を抽出するステップは、
前記冠状動脈二次元造影画像から血管骨格を抽出するステップと、
前記血管段の延在方向、および2点間の最短経路を取る原則に基づくステップと、
前記血管骨格に沿って、前記血管段の中心線を抽出するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の前記血管段の中心線を抽出するステップは、
前記冠状動脈二次元造影画像から血管骨格を抽出するステップと、
前記血管段の延在方向、および2点間の最短経路を取る原則に基づくステップと、
前記血管骨格に沿って、前記血管段の中心線を抽出するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0013】
任意選択的に、
前記の前記血管骨格に沿って、前記血管段の中心線を抽出するステップはさらに、
注目する前記血管段上に少なくとも一つのシードポイントを加えるステップと
前記最初と最後の点、シードポイントに基づき、前記血管骨格に沿って、新たに前記血管中心線を生成するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の前記血管骨格に沿って、前記血管段の中心線を抽出するステップはさらに、
注目する前記血管段上に少なくとも一つのシードポイントを加えるステップと
前記最初と最後の点、シードポイントに基づき、前記血管骨格に沿って、新たに前記血管中心線を生成するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0014】
任意選択的に、
前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとし、且つ、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
造影剤が冠状動脈入口である前記血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が前記血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップと、
第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップと、
v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度をそれぞれv1,..., vb,..., va,..., vN-1とするステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとし、且つ、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
造影剤が冠状動脈入口である前記血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が前記血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップと、
第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップと、
v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度をそれぞれv1,..., vb,..., va,..., vN-1とするステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0015】
任意選択的に、
前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとし、且つ、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1から第1フレームの画像までの時間差及び中心線の長さの差を順に求めるステップと、
v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1から第1フレーム画像までの血流速度をそれぞれ得るステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとし、且つ、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1から第1フレームの画像までの時間差及び中心線の長さの差を順に求めるステップと、
v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1から第1フレーム画像までの血流速度をそれぞれ得るステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0016】
任意選択的に、
前記の前記血流速度的最大値である拡張期の血流速度を選び取るステップは、
再帰的アルゴリズムまたはバブルソートにより、前記技術方案11または技術方案12から前記血流速度の最大値である拡張期の血流速度を選び取るステップ、又は、
再帰的アルゴリズムまたはバブルソートにより、前記技術方案11または技術方案12から前記血流速度の最大値である拡張期の血流速度を選び取り、且つ、前記血流速度の最小値である収縮期の血流速度を選び取るステップ
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の前記血流速度的最大値である拡張期の血流速度を選び取るステップは、
再帰的アルゴリズムまたはバブルソートにより、前記技術方案11または技術方案12から前記血流速度の最大値である拡張期の血流速度を選び取るステップ、又は、
再帰的アルゴリズムまたはバブルソートにより、前記技術方案11または技術方案12から前記血流速度の最大値である拡張期の血流速度を選び取り、且つ、前記血流速度の最小値である収縮期の血流速度を選び取るステップ
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0017】
任意選択的に、
前記の前記血管段の中心線を抽出するステップの後、前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとするステップの前にさらに、
少なくとも2つの体位の冠状動脈二次元造影画像組を読み取るステップと、
生理パラメータ、画像体位撮影角度を含む、前記血管段のジオメトリ構造情報を取得するステップと、
注目する前記血管段に対して図形処理を行うステップと、
前記血管段の血管輪郭線を抽出するステップと、
前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の前記血管段の中心線を抽出するステップの後、前記の任意の2つのフレーム冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記分段中心線の差を計算し、その差の値をΔLとするステップの前にさらに、
少なくとも2つの体位の冠状動脈二次元造影画像組を読み取るステップと、
生理パラメータ、画像体位撮影角度を含む、前記血管段のジオメトリ構造情報を取得するステップと、
注目する前記血管段に対して図形処理を行うステップと、
前記血管段の血管輪郭線を抽出するステップと、
前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0018】
任意選択的に、
前記の前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
前記三次元血管モデルに基づき、三次元血管モデルの中心線を取得し、前記冠状動脈二次元造影画像により抽出した中心線に対して補正を行い、且つ、中心線の差の値ΔLに対して補正を行い、ΔL'を得るステップと、
前記ΔL'と前記Δtの比に基づき、血流速度vを求めるステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
前記の前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
前記三次元血管モデルに基づき、三次元血管モデルの中心線を取得し、前記冠状動脈二次元造影画像により抽出した中心線に対して補正を行い、且つ、中心線の差の値ΔLに対して補正を行い、ΔL'を得るステップと、
前記ΔL'と前記Δtの比に基づき、血流速度vを求めるステップと、
を含む、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
【0019】
第二の態様において、本出願は以下の装置を提供する。
血流速度取得ユニット、大動脈圧波形取得ユニット、及び冠状動脈血管評定パラメータユニットを備え、
前記冠状動脈血管評定パラメータユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニットに接続され、
前記血流速度取得ユニットは血流速度vを取得するためのものであり、
前記大動脈圧波形取得ユニットは時間に伴い変わる大動脈圧波形をリアルタイムに取得するためのものであり、
前記冠状動脈血管評定パラメータユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニットが送信した前記血流速度v及び大動脈圧波形を受信し、それから生理パラメータに基づき、冠状動脈評定パラメータを取得するためのものである、
上記のいずれかに記載の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法に用いられる生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
血流速度取得ユニット、大動脈圧波形取得ユニット、及び冠状動脈血管評定パラメータユニットを備え、
前記冠状動脈血管評定パラメータユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニットに接続され、
前記血流速度取得ユニットは血流速度vを取得するためのものであり、
前記大動脈圧波形取得ユニットは時間に伴い変わる大動脈圧波形をリアルタイムに取得するためのものであり、
前記冠状動脈血管評定パラメータユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニットが送信した前記血流速度v及び大動脈圧波形を受信し、それから生理パラメータに基づき、冠状動脈評定パラメータを取得するためのものである、
上記のいずれかに記載の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法に用いられる生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
【0020】
任意選択的に、
さらに、
順に接続される画像読み取りユニット、血管段抽出ユニット、中心線抽出ユニット、前記画像読み取りユニットに接続される時間差ユニット、ジオメトリ情報取得ユニット、及び前記時間差ユニット、前記血流速度取得ユニットにそれぞれ接続される中心線差ユニットを備え、
前記中心線差ユニットは前記中心線抽出ユニットに接続され、
前記ジオメトリ情報取得ユニットは前記冠状動脈血管評定パラメータユニットに接続され、
前記画像読み取りユニットは少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るためのものであり、
前記血管段抽出ユニットは前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信し、前記画像中の注目する血管段を抽出するためのものであり、
前記中心線抽出ユニットは前記血管段抽出ユニットが送信した血管段を受信し、前記血管段の中心線を抽出するためのものであり、
前記時間差ユニットは前記画像読み取りユニットが送信した任意の2つのフレームの冠状動脈二次元造影画像を受信し、前記2つのフレームの冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとするものであり、
前記中心線差ユニットは前記中心線抽出ユニットが送信した2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた分段中心線を受信し、前記分段中心線の差を計算し、その差の値をΔLとするものであり、
前記血流速度取得ユニットは血流速度計算モジュール、拡張期の血流速度計算モジュールを備え、
前記血流速度計算モジュールは前記時間差ユニット、前記中心線差ユニットにそれぞれ接続され、
前記拡張期の血流速度計算モジュールは前記血流速度計算モジュールに接続され、
血流速度計算モジュールは前記時間差ユニット及び前記中心線差ユニットが送信した前記ΔLと前記Δtを受信し、前記ΔLと前記Δtの比に基づき、血流速度を求めるためのものであり、
前記拡張期の血流速度計算モジュールは前記血流速度計算モジュールが送信した血流速度値を受信し、前記血流速度の最大値である拡張期の血流速度を選び取るためのものであり、
前記ジオメトリ情報取得ユニットは前記画像読み取りユニットの冠状動脈二次元造影画像を受信し、患者の生理パラメータ及び画像撮影角度を取得し、且つ、生理パラメータ及び画像撮影角度を前記冠状動脈血管評定パラメータユニットに伝達するためのものである、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
さらに、
順に接続される画像読み取りユニット、血管段抽出ユニット、中心線抽出ユニット、前記画像読み取りユニットに接続される時間差ユニット、ジオメトリ情報取得ユニット、及び前記時間差ユニット、前記血流速度取得ユニットにそれぞれ接続される中心線差ユニットを備え、
前記中心線差ユニットは前記中心線抽出ユニットに接続され、
前記ジオメトリ情報取得ユニットは前記冠状動脈血管評定パラメータユニットに接続され、
前記画像読み取りユニットは少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るためのものであり、
前記血管段抽出ユニットは前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信し、前記画像中の注目する血管段を抽出するためのものであり、
前記中心線抽出ユニットは前記血管段抽出ユニットが送信した血管段を受信し、前記血管段の中心線を抽出するためのものであり、
前記時間差ユニットは前記画像読み取りユニットが送信した任意の2つのフレームの冠状動脈二次元造影画像を受信し、前記2つのフレームの冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとするものであり、
前記中心線差ユニットは前記中心線抽出ユニットが送信した2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた分段中心線を受信し、前記分段中心線の差を計算し、その差の値をΔLとするものであり、
前記血流速度取得ユニットは血流速度計算モジュール、拡張期の血流速度計算モジュールを備え、
前記血流速度計算モジュールは前記時間差ユニット、前記中心線差ユニットにそれぞれ接続され、
前記拡張期の血流速度計算モジュールは前記血流速度計算モジュールに接続され、
血流速度計算モジュールは前記時間差ユニット及び前記中心線差ユニットが送信した前記ΔLと前記Δtを受信し、前記ΔLと前記Δtの比に基づき、血流速度を求めるためのものであり、
前記拡張期の血流速度計算モジュールは前記血流速度計算モジュールが送信した血流速度値を受信し、前記血流速度の最大値である拡張期の血流速度を選び取るためのものであり、
前記ジオメトリ情報取得ユニットは前記画像読み取りユニットの冠状動脈二次元造影画像を受信し、患者の生理パラメータ及び画像撮影角度を取得し、且つ、生理パラメータ及び画像撮影角度を前記冠状動脈血管評定パラメータユニットに伝達するためのものである、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
【0021】
任意選択的に、
さらに、
前記画像読み取りユニットにそれぞれ接続される血管骨格抽出ユニット、三次元血管再構成ユニット、及び前記血管骨格抽出ユニットに接続される輪郭線抽出ユニットを備え、
前記三次元血管再構成ユニットは前記ジオメトリ情報取得ユニット、前記中心線抽出ユニット、前記輪郭線抽出ユニットに接続され、
前記血管骨格抽出ユニットは前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信し、前記画像中の血管骨格を抽出するためのものであり、
前記輪郭線抽出ユニットは前記血管骨格抽出ユニットの血管骨格を受信し、前記血管骨格に基づき、注目する前記血管段の輪郭線を抽出するためのものであり、
前記三次元血管再構成ユニットは前記輪郭線抽出ユニット、前記ジオメトリ情報取得ユニット、前記中心線抽出ユニットが送信した輪郭線、ジオメトリ構造情報及び中心線を受信するため、且つ、前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信するため、且つ、前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するためのものであり、
前記中心線抽出ユニットは前記三次元血管再構成ユニットの三次元血管モデルから新たに前記血管段の中心線を抽出し、且つ、新たに前記中心線の長さを取得するためのものである、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
さらに、
前記画像読み取りユニットにそれぞれ接続される血管骨格抽出ユニット、三次元血管再構成ユニット、及び前記血管骨格抽出ユニットに接続される輪郭線抽出ユニットを備え、
前記三次元血管再構成ユニットは前記ジオメトリ情報取得ユニット、前記中心線抽出ユニット、前記輪郭線抽出ユニットに接続され、
前記血管骨格抽出ユニットは前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信し、前記画像中の血管骨格を抽出するためのものであり、
前記輪郭線抽出ユニットは前記血管骨格抽出ユニットの血管骨格を受信し、前記血管骨格に基づき、注目する前記血管段の輪郭線を抽出するためのものであり、
前記三次元血管再構成ユニットは前記輪郭線抽出ユニット、前記ジオメトリ情報取得ユニット、前記中心線抽出ユニットが送信した輪郭線、ジオメトリ構造情報及び中心線を受信するため、且つ、前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信するため、且つ、前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するためのものであり、
前記中心線抽出ユニットは前記三次元血管再構成ユニットの三次元血管モデルから新たに前記血管段の中心線を抽出し、且つ、新たに前記中心線の長さを取得するためのものである、
上記の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置。
【0022】
第三の態様において、本出願は上記のいずれかに記載の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置を含む冠状動脈分析システムを提供する。
【0023】
第四の態様において、本出願はコンピュータプログラムがプロセッサにより実行される時に上記のいずれかに記載の生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法を実現するコンピュータ記憶媒体を提供する。
【0024】
本出願の実施例により提供する方案により少なくとも以下の有益な效果を奏することができる。
本出願は生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法を提供する。この方法によれば、血流速度v、大動脈圧波形、生理パラメータに基づき、冠状動脈血管評定パラメータを取得し、異なる性別、異なる既往歴の患者に対して冠状動脈血管評定パラメータの個別化測定を行うことができ、より高い狙い性を有し、冠状動脈血管評定パラメータ測定の正確性を高めることができる。
本出願は生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法を提供する。この方法によれば、血流速度v、大動脈圧波形、生理パラメータに基づき、冠状動脈血管評定パラメータを取得し、異なる性別、異なる既往歴の患者に対して冠状動脈血管評定パラメータの個別化測定を行うことができ、より高い狙い性を有し、冠状動脈血管評定パラメータ測定の正確性を高めることができる。
【図面の簡単な説明】
【0025】
ここで説明する図面は本発明をさらに理解するためのものであり、本発明の一部分を構成し、本発明の概略的な実施例及びその記載は本発明を説明するためのものであり、本発明に対する不適切な限定を構成するものではない。
【0026】
【発明を実施するための形態】
【0027】
本発明の目的、技術方案及び長所がより明確になるように、以下に本発明の具体的な実施例及び対応する図面を組み合わせて本発明の技術方案を明瞭、完全に記載する。勿論、記載する実施例は本発明の一部の実施例であり、全ての実施例ではない。本発明における実施例に基づき、当業者が創造的な労働を行わずに得られる他の実施例もすべて本発明の保護範囲に入ると理解されるべきである。
【0028】
以下に図面により本発明の複数の実施形態を開示する。明確に説明できるよう、多くの実践的な細部を以下の記載において一緒に説明する。しかし、これら実践的な細部は本発明を限定するために用いられるのではないと理解されるべきである。つまり、本発明の一部の実施形態において、これら実践的な細部は必須事項ではない。この他、図を簡素化するために、一部の慣用的な構造や構成部材については図中で簡単に概略的に示すだけにする。
【0029】
冠状動脈血管評定パラメータとしては、冠状動脈拡張期の血流速度CAIFR、および冠状動脈拡張期の微小血管抵抗指数CAIFMR等を含む。しかし個体群により生命兆候は異なるため、正常な数値の評定標準にも若干違いがある。例えば、高齢者の心筋微小循環機能は低く、血流速度は一般的に若者よりも低いため、もし業界の汎用評価標準を用いる場合、採用する血流速度は実際の値より大きめになり、CAIFMR等が過小評価され、CAIFMR等の測定の正確性が低下する。
【0030】
上記問題を解決するために、図1に示されるように、本出願は以下のステップを含む生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法を提供する。
【0031】
S000 生理パラメータを取得するステップ、
S010 血流速度vを取得するステップ、
S020 時間に伴い変わる大動脈圧波形をリアルタイムに取得するステップ、
S030 血流速度v、大動脈圧波形、生理パラメータに基づき、冠状動脈血管評定パラメータを取得するステップ。
S010 血流速度vを取得するステップ、
S020 時間に伴い変わる大動脈圧波形をリアルタイムに取得するステップ、
S030 血流速度v、大動脈圧波形、生理パラメータに基づき、冠状動脈血管評定パラメータを取得するステップ。
【0032】
本出願は生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法を提供する。この方法によれば、血流速度v、大動脈圧波形、生理パラメータに基づき、冠状動脈血管評定パラメータを取得し、異なる性別、異なる既往歴の患者に対して冠状動脈血管評定パラメータの個別化測定を行うことができ、より高い狙い性を有し、冠状動脈血管評定パラメータ測定の正確性を高めることができる。
【0033】
生理パラメータに基づき冠状動脈血管評定パラメータを取得する内容であれば、いずれも本出願の保護範囲内に入る。以下に冠状動脈評定パラメータが拡張期の微小血管抵抗指数CAIFMRである場合の具体的な取得方式について説明する。
【0034】
［実施例1］
多くの実験により、高血圧の既往歴、糖尿病の既往歴および性別はいずれも冠状動脈血管評定パラメータ計算の正確性に影響があることが証明されている。したがって、S030における冠状動脈評定パラメータが拡張期の微小血管抵抗指数CAIFMRである場合、具体的に以下のステップを含む。
血流速度vの最大値である拡張期の最大血流速度vmax、を選び取るステップ、好ましくは、本出願は、再帰的アルゴリズムまたはバブルソートにより、血流速度vの最大値を選び取る、
vmax対応する期間は拡張期であり、大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップ、
CAIFMRは、
であり、
は
であり、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vmaxは拡張期の最大血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、k=a×b、aは糖尿病特徴値を表し、bは高血圧特徴値を表し、cは性別を表す。
多くの実験により、高血圧の既往歴、糖尿病の既往歴および性別はいずれも冠状動脈血管評定パラメータ計算の正確性に影響があることが証明されている。したがって、S030における冠状動脈評定パラメータが拡張期の微小血管抵抗指数CAIFMRである場合、具体的に以下のステップを含む。
血流速度vの最大値である拡張期の最大血流速度vmax、を選び取るステップ、好ましくは、本出願は、再帰的アルゴリズムまたはバブルソートにより、血流速度vの最大値を選び取る、
vmax対応する期間は拡張期であり、大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップ、
CAIFMRは、
但し、
【0035】
本出願の一つの実施例において、患者に糖尿病がない場合、0.5≦a≦1とされ、好ましくは、a=1とされ、患者に糖尿病がある場合、1＜a≦2とされ、好ましくは、a=2とされ、
患者の血圧値が90mmHg以上である場合、1＜b≦1.5とされ、好ましくは、b=1.5とされ、患者の血圧値が90mmHg未満である場合、0.5≦b≦1とされ、好ましくは、b=1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=3～10とされ、好ましくは、c=5とされる。
患者の血圧値が90mmHg以上である場合、1＜b≦1.5とされ、好ましくは、b=1.5とされ、患者の血圧値が90mmHg未満である場合、0.5≦b≦1とされ、好ましくは、b=1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=3～10とされ、好ましくは、c=5とされる。
【0036】
図2に示されるように、本出願の一つの実施例において、二次元造影画像により血流速度を取得する時、S010は以下を含む。
S100A 少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ、
S200A 冠状動脈二次元造影画像組から、注目する血管段を抽出するステップ、
S300A 血管段の中心線を抽出するステップ、
S400A 任意の2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、分段中心線の差を計算し、その差の値をΔLとし、
但し、Δt=m×fpsであり、各組の冠状動脈二次元造影画像組には連続再生される複数のフレーム冠状動脈二次元造影画像が含まれるため、mは各組の冠状動脈二次元造影画像組において、選び取った2つのフレームの冠状動脈二次元造影画像が存在する所のフレーム数の差の値を表し、fpsは隣接する2つのフレーム画像の間で切り替わる間隔時間を表し、好ましくは、fps=1/15秒とされ、
ΔLとΔtの比に基づき、血流速度vを求めるステップ。
S100A 少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ、
S200A 冠状動脈二次元造影画像組から、注目する血管段を抽出するステップ、
S300A 血管段の中心線を抽出するステップ、
S400A 任意の2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、分段中心線の差を計算し、その差の値をΔLとし、
但し、Δt=m×fpsであり、各組の冠状動脈二次元造影画像組には連続再生される複数のフレーム冠状動脈二次元造影画像が含まれるため、mは各組の冠状動脈二次元造影画像組において、選び取った2つのフレームの冠状動脈二次元造影画像が存在する所のフレーム数の差の値を表し、fpsは隣接する2つのフレーム画像の間で切り替わる間隔時間を表し、好ましくは、fps=1/15秒とされ、
ΔLとΔtの比に基づき、血流速度vを求めるステップ。
【0037】
図3に示されるように、本出願の一つの実施例において、三次元モデリングにより血流速度を取得する時、S010は以下のステップを含む。
S100B 少なくとも2つの体位の冠状動脈二次元造影画像組を読み取るステップ、
S200B 冠状動脈二次元造影画像組から、注目する血管段を抽出するステップ、
S300B 血管段のジオメトリ構造情報を取得し、そして血管段の中心線を抽出するステップ、
S400B 注目する血管段に対して図形処理を行うステップ、
S500B 血管段の血管輪郭線を抽出するステップ、
S600B 血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するステップ、
S700B 任意の2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとし、三次元血管モデルに基づき、三次元血管モデルの中心線を取得し、冠状動脈二次元造影画像により抽出した中心線に対して補正を行い、且つ、補正後の分段中心線の差を計算し、その差の値をΔL'とし、ΔL'とΔtの比に基づき、血流速度vを求めるステップ、
但し、Δtは、Δt=m×fpsであり、各組の冠状動脈二次元造影画像組には連続再生される複数のフレーム冠状動脈二次元造影画像が含まれるため、mは各組の冠状動脈二次元造影画像組において、選び取った2つのフレームの冠状動脈二次元造影画像が存在する所のフレーム数の差の値を表し、fpsは隣接する2つのフレーム画像の間で切り替わる間隔時間を表し、好ましくは、fps=1/15秒とされる。
S100B 少なくとも2つの体位の冠状動脈二次元造影画像組を読み取るステップ、
S200B 冠状動脈二次元造影画像組から、注目する血管段を抽出するステップ、
S300B 血管段のジオメトリ構造情報を取得し、そして血管段の中心線を抽出するステップ、
S400B 注目する血管段に対して図形処理を行うステップ、
S500B 血管段の血管輪郭線を抽出するステップ、
S600B 血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するステップ、
S700B 任意の2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとし、三次元血管モデルに基づき、三次元血管モデルの中心線を取得し、冠状動脈二次元造影画像により抽出した中心線に対して補正を行い、且つ、補正後の分段中心線の差を計算し、その差の値をΔL'とし、ΔL'とΔtの比に基づき、血流速度vを求めるステップ、
但し、Δtは、Δt=m×fpsであり、各組の冠状動脈二次元造影画像組には連続再生される複数のフレーム冠状動脈二次元造影画像が含まれるため、mは各組の冠状動脈二次元造影画像組において、選び取った2つのフレームの冠状動脈二次元造影画像が存在する所のフレーム数の差の値を表し、fpsは隣接する2つのフレーム画像の間で切り替わる間隔時間を表し、好ましくは、fps=1/15秒とされる。
【0038】
本出願の一つの実施例において、S100A又はS100Bは以下のステップを含む。
無線または有線方式により造影画像撮影装置または病院プラットフォームから、少なくとも一つの体位の冠状動脈二次元造影画像組を直接読み取るステップ、又は
記憶装置により、少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ。
無線または有線方式により造影画像撮影装置または病院プラットフォームから、少なくとも一つの体位の冠状動脈二次元造影画像組を直接読み取るステップ、又は
記憶装置により、少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ。
【0039】
図4に示されるように、本出願の一つの実施例において、S200A又はS200Bは以下のステップを含む。
S210 冠状動脈二次元造影画像組から、Nフレーム冠状動脈二次元造影画像を選び取るステップ、
S220 冠状動脈二次元造影画像において、注目する血管の最初と最後の点をピックアップし、注目する血管段を取得するステップ。
好ましくは、S200A又はS200Bさらに以下を含む。
カテーテルが現れる第1フレーム冠状動脈二次元造影画像を参照画像としてていぎし、完全な冠状動脈が現れた第kフレーム冠状動脈二次元造影画像をターゲット画像として定義し、但しkは1よりも大きい正整数であり、
参照画像からターゲット画像を差し引き、カテーテルの特徴点Oを抽出し、好ましくは、一部の静態ノイズを除去し、さらに、平均値フィルタリングを採用して、一部の動態ノイズを除去し、そして、濃度ヒストグラム解析により、閾値を用いてさらにノイズ除去し、ターゲット画像から参照画像を差し引き、冠状動脈が存在する位置の区域画像を抽出し、区域画像はカテーテルの特徴点をシードポイントとして動的成長を行い、注目する血管段画像を取得する。
S210 冠状動脈二次元造影画像組から、Nフレーム冠状動脈二次元造影画像を選び取るステップ、
S220 冠状動脈二次元造影画像において、注目する血管の最初と最後の点をピックアップし、注目する血管段を取得するステップ。
好ましくは、S200A又はS200Bさらに以下を含む。
カテーテルが現れる第1フレーム冠状動脈二次元造影画像を参照画像としてていぎし、完全な冠状動脈が現れた第kフレーム冠状動脈二次元造影画像をターゲット画像として定義し、但しkは1よりも大きい正整数であり、
参照画像からターゲット画像を差し引き、カテーテルの特徴点Oを抽出し、好ましくは、一部の静態ノイズを除去し、さらに、平均値フィルタリングを採用して、一部の動態ノイズを除去し、そして、濃度ヒストグラム解析により、閾値を用いてさらにノイズ除去し、ターゲット画像から参照画像を差し引き、冠状動脈が存在する位置の区域画像を抽出し、区域画像はカテーテルの特徴点をシードポイントとして動的成長を行い、注目する血管段画像を取得する。
【0040】
図5に示されるように、本出願の一つの実施例において、S300A又はS300Bは以下のステップを含む。
S310 冠状動脈二次元造影画像から血管骨格を抽出するステップ、
S320 血管段の延在方向、および2点の間の最短経路を取得する原則に基づくステップ、
S330 血管骨格に沿って、血管段の中心線を抽出するステップ。
S310 冠状動脈二次元造影画像から血管骨格を抽出するステップ、
S320 血管段の延在方向、および2点の間の最短経路を取得する原則に基づくステップ、
S330 血管骨格に沿って、血管段の中心線を抽出するステップ。
【0041】
図6に示されるように、本出願の一つの実施例において、S330はさらに以下のステップを含む。
S331 注目する血管段上に少なくとも一つのシードポイントを加えるステップ、
S332 最初と最後の点、シードポイントに基づき、血管骨格に沿って、新たに血管中心線を生成するステップ。
S331 注目する血管段上に少なくとも一つのシードポイントを加えるステップ、
S332 最初と最後の点、シードポイントに基づき、血管骨格に沿って、新たに血管中心線を生成するステップ。
【0042】
本出願の一つの実施例において、S400A又はS700Bは二種類の取得方法を含み、方法（1）は図7に示されるように、以下のステップを含む。
S410I 造影剤が冠状動脈入口である血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップ、
S420I 第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップ。
S430I v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度をそれぞれv1,..., vb,..., va,..., vN-1とするステップ。
S410I 造影剤が冠状動脈入口である血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップ、
S420I 第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップ。
S430I v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度をそれぞれv1,..., vb,..., va,..., vN-1とするステップ。
【0043】
本出願の一つの実施例において、S400A又はS700Bは二種類の取得方法を含み、方法（2）は図8に示されるように、以下のステップを含む。
S410II 造影剤が冠状動脈入口である血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップ、
S420II 第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1フレームから第1フレームの画像までの時間差及び中心線の長さの差を順に求めるステップ、
S430II v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1フレームから第1フレーム画像までの血流速度をそれぞれ得るステップ。
S410II 造影剤が冠状動脈入口である血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップ、
S420II 第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1フレームから第1フレームの画像までの時間差及び中心線の長さの差を順に求めるステップ、
S430II v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレームから第bフレームまで、第N-1フレームから第b-1フレームまで、...、第N-b-aフレームから第N-aフレームまで、...、第N-b+1フレームから第1フレーム画像までの血流速度をそれぞれ得るステップ。
【0044】
本出願は再帰的アルゴリズムまたはバブルソートにより、血流速度の最小値である収縮期の血流速度を選び取ることもできる。
【0045】
［実施例2］
図9に示されるように、本出願は上記のいずれかの生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法に用いられる、以下の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置を提供する。
血流速度取得ユニット1、大動脈圧波形取得ユニット2、及び冠状動脈血管評定パラメータユニット3を備え、
冠状動脈血管評定パラメータユニット3は血流速度取得ユニット1、大動脈圧波形取得ユニット2に接続され、
血流速度取得ユニット1は血流速度vを取得するためのものであり、
大動脈圧波形取得ユニット2は時間に伴い変わる大動脈圧波形をリアルタイムに取得するためのものであり、
冠状動脈血管評定パラメータユニット3は血流速度取得ユニット1、大動脈圧波形取得ユニット2が送信した血流速度v及び大動脈圧波形を受信し、それから生理パラメータに基づき、冠状動脈評定パラメータを取得するためのものである。
図9に示されるように、本出願は上記のいずれかの生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法に用いられる、以下の生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置を提供する。
血流速度取得ユニット1、大動脈圧波形取得ユニット2、及び冠状動脈血管評定パラメータユニット3を備え、
冠状動脈血管評定パラメータユニット3は血流速度取得ユニット1、大動脈圧波形取得ユニット2に接続され、
血流速度取得ユニット1は血流速度vを取得するためのものであり、
大動脈圧波形取得ユニット2は時間に伴い変わる大動脈圧波形をリアルタイムに取得するためのものであり、
冠状動脈血管評定パラメータユニット3は血流速度取得ユニット1、大動脈圧波形取得ユニット2が送信した血流速度v及び大動脈圧波形を受信し、それから生理パラメータに基づき、冠状動脈評定パラメータを取得するためのものである。
【0046】
図10に示されるように、本出願の一つの実施例において、さらに、
順に接続される画像読み取りユニット4、血管段抽出ユニット5、中心線抽出ユニット6、画像読み取りユニット4に接続される時間差ユニット7、ジオメトリ情報取得ユニット8を備え、
中心線差ユニット9は中心線抽出ユニット6に接続され、
時間差ユニット7、中心線差ユニット9はそれぞれ血流速度取得ユニット1に接続され、
ジオメトリ情報取得ユニット8は冠状動脈血管評定パラメータユニット3に接続され、
画像読み取りユニット4は少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るためのものであり、
血管段抽出ユニット5は画像読み取りユニット4が送信した冠状動脈二次元造影画像を受信し、画像中の注目する血管段を抽出するためのものであり、
中心線抽出ユニット6は血管段抽出ユニット5が送信した血管段を受信し、血管段の中心線を抽出するためのものであり、
時間差ユニット7は画像読み取りユニット4が送信した任意の2つのフレームの冠状動脈二次元造影画像を受信し、2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとするものであり、
中心線差ユニット9は中心線抽出ユニット6が送信した2つのフレームの冠状動脈二次元造影画像中において造影剤が血管段を流れた分段中心線を受信し、分段中心線の差を計算し、その差の値をΔLとするものであり、
血流速度取得ユニット1は血流速度計算モジュール101、拡張期の血流速度計算モジュール102を備え、
血流速度計算モジュール101は時間差ユニット7、中心線差ユニット8にそれぞれ接続され、
拡張期の血流速度計算モジュール102は血流速度計算モジュール101に接続され、
血流速度計算モジュール101は時間差ユニット7及び中心線差ユニット8が送信したΔLとΔtを受信し、ΔLとΔtの比に基づき、血流速度を求めるためのものであり、
拡張期の血流速度計算モジュール102は血流速度計算モジュールが送信した血流速度値を受信し、血流速度の最大値である拡張期の血流速度を選び取るためのものであり、
ジオメトリ情報取得ユニット8は画像読み取りユニット4の冠状動脈二次元造影画像を受信し、患者の生理パラメータ及び画像撮影角度を取得し、且つ、生理パラメータ及び画像撮影角度を冠状動脈血管評定パラメータユニット3に伝達するためのものである。
順に接続される画像読み取りユニット4、血管段抽出ユニット5、中心線抽出ユニット6、画像読み取りユニット4に接続される時間差ユニット7、ジオメトリ情報取得ユニット8を備え、
中心線差ユニット9は中心線抽出ユニット6に接続され、
時間差ユニット7、中心線差ユニット9はそれぞれ血流速度取得ユニット1に接続され、
ジオメトリ情報取得ユニット8は冠状動脈血管評定パラメータユニット3に接続され、
画像読み取りユニット4は少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るためのものであり、
血管段抽出ユニット5は画像読み取りユニット4が送信した冠状動脈二次元造影画像を受信し、画像中の注目する血管段を抽出するためのものであり、
中心線抽出ユニット6は血管段抽出ユニット5が送信した血管段を受信し、血管段の中心線を抽出するためのものであり、
時間差ユニット7は画像読み取りユニット4が送信した任意の2つのフレームの冠状動脈二次元造影画像を受信し、2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとするものであり、
中心線差ユニット9は中心線抽出ユニット6が送信した2つのフレームの冠状動脈二次元造影画像中において造影剤が血管段を流れた分段中心線を受信し、分段中心線の差を計算し、その差の値をΔLとするものであり、
血流速度取得ユニット1は血流速度計算モジュール101、拡張期の血流速度計算モジュール102を備え、
血流速度計算モジュール101は時間差ユニット7、中心線差ユニット8にそれぞれ接続され、
拡張期の血流速度計算モジュール102は血流速度計算モジュール101に接続され、
血流速度計算モジュール101は時間差ユニット7及び中心線差ユニット8が送信したΔLとΔtを受信し、ΔLとΔtの比に基づき、血流速度を求めるためのものであり、
拡張期の血流速度計算モジュール102は血流速度計算モジュールが送信した血流速度値を受信し、血流速度の最大値である拡張期の血流速度を選び取るためのものであり、
ジオメトリ情報取得ユニット8は画像読み取りユニット4の冠状動脈二次元造影画像を受信し、患者の生理パラメータ及び画像撮影角度を取得し、且つ、生理パラメータ及び画像撮影角度を冠状動脈血管評定パラメータユニット3に伝達するためのものである。
【0047】
本出願の一つの実施例において、さらに、
画像読み取りユニット4にそれぞれ接続される血管骨格抽出ユニット10、三次元血管再構成ユニット11、及び血管骨格抽出ユニット10に接続される輪郭線抽出ユニット12を備え、
三次元血管再構成ユニット11はジオメトリ情報取得ユニット8、中心線抽出ユニット6、輪郭線抽出ユニット12に接続され、
血管骨格抽出ユニット10は画像読み取りユニット4が送信した冠状動脈二次元造影画像を受信し、画像中の血管骨格を抽出するためのものであり、
輪郭線抽出ユニット12は血管骨格抽出ユニット10の血管骨格を受信し、血管骨格に基づき、注目する血管段の輪郭線を抽出するためのものであり、
三次元血管再構成ユニット11は輪郭線抽出ユニット12、ジオメトリ情報取得ユニット8、中心線抽出ユニット6が送信した輪郭線、ジオメトリ構造情報及び中心線を受信するため、且つ、画像読み取りユニットが送信した冠状動脈二次元造影画像を受信するため、且つ、血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するためのものであり、
中心線抽出ユニット6は三次元血管再構成ユニット11の三次元血管モデルから新たに血管段の中心線を抽出し、且つ、新たに中心線の長さを取得するためのものである。
画像読み取りユニット4にそれぞれ接続される血管骨格抽出ユニット10、三次元血管再構成ユニット11、及び血管骨格抽出ユニット10に接続される輪郭線抽出ユニット12を備え、
三次元血管再構成ユニット11はジオメトリ情報取得ユニット8、中心線抽出ユニット6、輪郭線抽出ユニット12に接続され、
血管骨格抽出ユニット10は画像読み取りユニット4が送信した冠状動脈二次元造影画像を受信し、画像中の血管骨格を抽出するためのものであり、
輪郭線抽出ユニット12は血管骨格抽出ユニット10の血管骨格を受信し、血管骨格に基づき、注目する血管段の輪郭線を抽出するためのものであり、
三次元血管再構成ユニット11は輪郭線抽出ユニット12、ジオメトリ情報取得ユニット8、中心線抽出ユニット6が送信した輪郭線、ジオメトリ構造情報及び中心線を受信するため、且つ、画像読み取りユニットが送信した冠状動脈二次元造影画像を受信するため、且つ、血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するためのものであり、
中心線抽出ユニット6は三次元血管再構成ユニット11の三次元血管モデルから新たに血管段の中心線を抽出し、且つ、新たに中心線の長さを取得するためのものである。
【0048】
本出願は上記のいずれかの生理パラメータに基づき冠状動脈血管評定パラメータを取得する装置を含む冠状動脈分析システムを提供する。
【0049】
本出願はコンピュータプログラムがプロセッサにより実行される時に上記のいずれかの生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法を実現するコンピュータ記憶媒体を提供する。
【0050】
当業者の認識として、本発明の各々の態様はシステム、方法又はコンピュータプログラム製品として実現することができる。したがって、本発明の各々の態様は具体的に以下の形式として実現できる。即ち、完全なハードウエア実施形態、完全なソフトウエア実施形態(ファームウエア、常駐ソフトウエア、マイクロコード等を含む)、 又はハードウエア及びソフトウエアの態様を組み合わせた実施形態、ここでは総じて“回路”、“モジュール”又は“システム”と称することができる。この他、一部の実施例において、本発明の各々の態様はさらに一つの又は複数のコンピュータ読み取り可能な媒体におけるコンピュータプログラム製品の形式として実現することができ、当該コンピュータ読み取り可能な媒体にはコンピュータ読み取り可能なプログラムコードが含まれる。本発明の実施例の方法及び/又はシステムの実施形態は手動的、自動的又はその組み合わせ方式により選択されたタスクを実行又は完成させることに関することができる。
【0051】
例えば、本発明の実施例に基づき選択されるタスクを実行するためのハードウエアをチップ又は回路として実現することができる。ソフトウエアとして、発明の実施例に基づき選択されるタスクを、コンピュータにより使用される如何なる適切な操作システムが実行する複数のソフトウエアコマンドとして実現できる。本発明の例示的な実施例において、本明細書の方法及び/又はシステムの例示的な実施例に基づく一つの又は複数のタスクをデータプロセッサにより実行できる。例えば、複数のコマンドを実行するための計算プラットフォーム。任意選択的に、当該データプロセッサはコマンド及び/又はデータを記憶するための揮発性メモリ及び/又はコマンド及び/又はデータを記憶するための不揮発性メモリを含む。例えば、磁気ハードディスク及び/又はリムーバブル媒体。任意選択的に、ネットワーク接続も提供する。任意選択的に、ディスプレイ及び/又は例えばキーボードやマウスなどのユーザ入力機器も提供する。
【0052】
一つの又は複数のコンピュータ読み取り可能な如何なる組み合わせも利用できる。コンピュータ読み取り可能な媒体はコンピュータ読み取り可能な信号媒体又はコンピュータ読み取り可能な記憶媒体とすることができる。コンピュータ読み取り可能な記憶媒体は、例えば電気、磁気、光、電磁、赤外線、又は半導体のシステム、装置又はデバイス、若しくは任意の以上の組み合わせとすることができるが、これらに限られない。コンピュータ読み取り可能な記憶媒体の更なる具体的な例(全てを挙げるものではない)は以下を含む。
【0053】
一つの又は複数の導線を有する電気接続、ポータブルコンピュータディスク、ハードディスク、ランダムアクセスメモリ(RAM)、リードオンリーメモリ(ROM)、書き換え可能なリードオンリーメモリ(EPROM又はフラッシュメモリ)、光ファイバー、コンパクトディスクリードオンリーメモリ(CD-ROM)、光記憶装置、磁気記憶装置、若しくは上記の任意の適切な組み合わせ。本明細書において、コンピュータ読み取り可能な記憶媒体はプログラムを含む又は記憶する如何なる有形媒体とすることができ、当該プログラムはコマンド実行システム、装置又はデバイスにより使用される若しくはこれらと組み合わせて使用することができる。
【0054】
コンピュータ読み取り可能な信号媒体は、ベースバンドに含まれる若しくは搬送波の一部として伝播できるデータ信号とすることができ、コンピュータ読み取り可能なプログラムコードを搭載できる。このように伝播されるデータ信号は複数の形式を採用でき、電磁信号、光信号又は上記の任意の適切な組み合わせを含むが、これらに限られない。コンピュータ読み取り可能な信号媒体はさらに、コンピュータ読み取り可能な記憶媒体以外の如何なるコンピュータ読み取り可能な媒体とすることができ、当該コンピュータ読み取り可能な媒体は、コマンド実行システム、装置又はデバイスにより使用される若しくはこれらと組み合わせて使用するプログラムを送信、伝播又は伝送できる。
【0055】
コンピュータ読み取り可能な媒体に含まれるプログラムコードは如何なる適切な媒体を用いて伝送することができ、無線、有線、光ケーブル、RF等、若しくは上記の任意の適切な組み合わせが含まれる（但しこれらに含まれない）。
【0056】
例えば、一つの又は複数のプログラミング言語の如何なる組み合わせでも本発明の各態様に用いられる操作を実行するためのコンピュータプログラムコードをプログラミングすることができ、例えばJava、Smalltalk、C++等のターゲット型プログラミング言語と通常プロセスのプログラミング言語、例えば"C"プログラミング言語又は類似したプログラミング言語を含む。プログラムコードは完全にユーザコンピュータ上で実行することも、一部をユーザコンピュータ上で実行することもでき、一つの独立したソフトウェアパッケージとして実行することも、一部をユーザコンピュータ上で、一部をリモートコンピュータ上で実行することもでき、若しくは完全にリモートコンピュータ又はサーバ上で実行することもできる。リモートコンピュータにかかる場合、リモートコンピュータは、ローカルエリアネットワーク(LAN)又はワイドエリアネットワーク(WAN)を含む任意の種類のネットワークによりユーザコンピュータに接続でき、若しくは、外部コンピュータに接続できる(例えばインターネットサービスプロバイダーによりインターネットを通じで接続する)。
【0057】
フローチャート及び/又はブロック図の各ブロック及びフローチャート及び/又はブロック図中の各ブロックの組み合わせは、いずれもコンピュータプログラムコマンドにより実現できると理解されるべきである。これらコンピュータプログラムコマンドは汎用コンピュータ、専用コンピュータ又は他のプログラム可能なデータ処理装置のプロセッサに提供され、一つの機器として生産することができ、これにより、これらコンピュータプログラムコマンドはコンピュータ又は他のプログラム可能なデータ処理装置のプロセッサの実行時に、フローチャート及び/又はブロック図中の一つの又は複数のブロック中で規定される機能/動作を実現できる装置とすることができる。
【0058】
これらコンピュータプログラムコマンドは、コンピュータ読み取り可能な媒体に記憶することもでき、これらコマンドは、コンピュータ、他のプログラム可能なデータ処理装置、又は他の設備を特定な方式で作動させ、そして、コンピュータ読み取り可能な媒体に記憶されたコマンドはフローチャート及び/又はブロック図中の一つの又は複数のブロック中で規定される機能/動作を実現できるコマンドを含む製品(article of manufacture)とすることができる。
【0059】
さらに、コンピュータ、他のプログラム可能なデータ処理設備又は他の設備で一連の操作ステップを実行させるように、コンピュータ(例えば、冠状動脈分析システム)又は他のプログラム可能なデータ処理設備にコンピュータプログラムコマンドをロードすることができ、これにより、コンピュータ、他のプログラム可能な装置又は他の設備で実行されるコマンドが、フローチャート及び/又は一つの又は複数のブロック図のブロック中で指定される機能/動作を実現するためのプロセスを提供するように、コンピュータが実行するプロセスを生成することができる。
【0060】
本発明の以上の具体的な実例により、本発明の目的、技術方案及び有益な效果について、より詳細に説明した。以上は本発明の具体的な実施例に過ぎず、本発明を限定するために用いられるものではなく、本発明の思想及び原則を逸脱しない限り、行われる如何なる修正、均等差し替え、改良なども、全て本発明の保護範囲に包含されると理解されるべきである。
【符号の説明】
【0061】
血流速度取得ユニット1
血流速度計算モジュール101
拡張期の血流速度計算モジュール102
大動脈圧波形取得ユニット2
冠状動脈血管評定パラメータユニット3
画像読み取りユニット4
血管段抽出ユニット5
中心線抽出ユニット6
時間差ユニット7
ジオメトリ情報取得ユニット8
中心線差ユニット9
血管骨格抽出ユニット10
三次元血管再構成ユニット11
輪郭線抽出ユニット12
血流速度計算モジュール101
拡張期の血流速度計算モジュール102
大動脈圧波形取得ユニット2
冠状動脈血管評定パラメータユニット3
画像読み取りユニット4
血管段抽出ユニット5
中心線抽出ユニット6
時間差ユニット7
ジオメトリ情報取得ユニット8
中心線差ユニット9
血管骨格抽出ユニット10
三次元血管再構成ユニット11
輪郭線抽出ユニット12
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
