JPWO2021087968A5 - - Google Patents
Download PDFInfo
- Publication number
- JPWO2021087968A5 JPWO2021087968A5 JP2022525374A JP2022525374A JPWO2021087968A5 JP WO2021087968 A5 JPWO2021087968 A5 JP WO2021087968A5 JP 2022525374 A JP2022525374 A JP 2022525374A JP 2022525374 A JP2022525374 A JP 2022525374A JP WO2021087968 A5 JPWO2021087968 A5 JP WO2021087968A5
- Authority
- JP
- Japan
- Prior art keywords
- blood flow
- flow velocity
- unit
- difference
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000017531 blood circulation Effects 0.000 claims description 92
- 210000004351 coronary vessel Anatomy 0.000 claims description 29
- 238000000034 method Methods 0.000 claims description 26
- 230000003205 diastolic effect Effects 0.000 claims description 25
- 206010020565 Hyperaemia Diseases 0.000 claims description 20
- 230000002792 vascular Effects 0.000 claims description 20
- 239000002872 contrast media Substances 0.000 claims description 13
- 210000004204 blood vessel Anatomy 0.000 claims description 12
- 230000000544 hyperemic effect Effects 0.000 claims description 12
- 238000002583 angiography Methods 0.000 claims description 11
- 238000002586 coronary angiography Methods 0.000 claims description 10
- 230000036772 blood pressure Effects 0.000 claims description 6
- 230000035487 diastolic blood pressure Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims 23
- 238000004364 calculation method Methods 0.000 claims 12
- 230000001105 regulatory effect Effects 0.000 claims 6
- 206010012601 diabetes mellitus Diseases 0.000 claims 5
- 238000003384 imaging method Methods 0.000 claims 3
- 230000002194 synthesizing effect Effects 0.000 claims 2
- 206010020772 Hypertension Diseases 0.000 claims 1
- 230000000747 cardiac effect Effects 0.000 claims 1
- 238000004590 computer program Methods 0.000 claims 1
- 230000001144 postural effect Effects 0.000 description 1
Description
任意選択的に、
前記血流速v、大動脈圧波形、生理パラメータに基づき、最大充血状態における微小血管抵抗指数iFMRを取得するステップは、
前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取るステップと、
前記vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップと、
を含み、
iFMRは、
であり、
は、
であり、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vhは最大充血状態における血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、k、cは影響パラメータk=1~3、c=0~10を表す、
上記の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。
前記血流速v、大動脈圧波形、生理パラメータに基づき、最大充血状態における微小血管抵抗指数iFMRを取得するステップは、
前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取るステップと、
前記vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップと、
を含み、
iFMRは、
であり、
は、
であり、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vhは最大充血状態における血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、k、cは影響パラメータk=1~3、c=0~10を表す、
上記の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。
任意選択的に、
前記任意の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とするステップと、
を含む、
上記の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。
[実施例2]
本出願は以下のステップを含む微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法を提供する。
S100 血流速度v、大動脈圧波形、生理パラメータに基づき、拡張期の微小血管抵抗指数iFMRを取得するステップ、
当該ステップにおいて、
1)図2に示されるように、二次元造影画像により血流速度を取得する時、血流速度vを取得する方法は以下のステップを含む。
S110 少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ、
S120 冠状動脈二次元造影画像組から、注目する血管段を取得するステップ、
当該ステップの具体的な方法は図3に示されるように、以下のステップを含む。
S121 冠状動脈二次元造影画像組からNフレーム冠状動脈二次元造影画像を選び取るステップ、
S122 冠状動脈二次元造影画像において、注目する血管の最初と最後の点をピックアップし、注目する血管段を取得するステップ、
S130 血管段の中心線を取得するステップ、
当該ステップの具体的方法は図4に示されるように以下のステップを含む。
S131 冠状動脈二次元造影画像から血管骨格を抽出するステップ、
S132 血管段の延在方向、および2点間の最短経路を取る原則に基づくステップ、
S133 血管骨格に沿って、血管段の中心線を取得するステップ、
S140 任意の2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、分段中心線の差を計算し、その差の値をΔLとするステップ、
S150 ΔLとΔtの比に基づき、血流速度vを求めるステップ、
2)大動脈圧波形の取得方法は以下のステップを含む。
侵襲血圧センサ又は非侵襲血圧器によりリアルタイムの大動脈圧Paj(但しjは1以上の正整数である)を取得し、その後、時間に基づき大動脈圧波形を生成するステップ、
3)拡張期の微小血管抵抗指数iFMRを取得するステップは以下のステップを含む。
S160 S150で解を求めて得られた血流速度vにおいて最大値である拡張期の最大血流速度vmaxを選び取るステップ、
S170 vmaxに対応する期間は拡張期であり、大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップ、即ち、
である、
S180 計算式
に基づきiFMR値を取得するステップ、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vhは最大充血状態における血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、好ましくは、本出願は再帰的アルゴリズムまたはバブルソートにより血流速度vの最大値を選び取り、k、cは影響パラメータk=1~3、c=0~10を表す、
S200 拡張期の微小血管抵抗指数がiFMR<Kである場合、調節パラメータをr=1とし、且つ、拡張期の微小血管抵抗指数がiFMR≧Kである場合、調節パラメータを
(但しKは100未満の正数であり、好ましくは、K=50とする)とするステップ。
本出願は以下のステップを含む微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法を提供する。
S100 血流速度v、大動脈圧波形、生理パラメータに基づき、拡張期の微小血管抵抗指数iFMRを取得するステップ、
当該ステップにおいて、
1)図2に示されるように、二次元造影画像により血流速度を取得する時、血流速度vを取得する方法は以下のステップを含む。
S110 少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップ、
S120 冠状動脈二次元造影画像組から、注目する血管段を取得するステップ、
当該ステップの具体的な方法は図3に示されるように、以下のステップを含む。
S121 冠状動脈二次元造影画像組からNフレーム冠状動脈二次元造影画像を選び取るステップ、
S122 冠状動脈二次元造影画像において、注目する血管の最初と最後の点をピックアップし、注目する血管段を取得するステップ、
S130 血管段の中心線を取得するステップ、
当該ステップの具体的方法は図4に示されるように以下のステップを含む。
S131 冠状動脈二次元造影画像から血管骨格を抽出するステップ、
S132 血管段の延在方向、および2点間の最短経路を取る原則に基づくステップ、
S133 血管骨格に沿って、血管段の中心線を取得するステップ、
S140 任意の2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、分段中心線の差を計算し、その差の値をΔLとするステップ、
S150 ΔLとΔtの比に基づき、血流速度vを求めるステップ、
2)大動脈圧波形の取得方法は以下のステップを含む。
侵襲血圧センサ又は非侵襲血圧器によりリアルタイムの大動脈圧Paj(但しjは1以上の正整数である)を取得し、その後、時間に基づき大動脈圧波形を生成するステップ、
3)拡張期の微小血管抵抗指数iFMRを取得するステップは以下のステップを含む。
S160 S150で解を求めて得られた血流速度vにおいて最大値である拡張期の最大血流速度vmaxを選び取るステップ、
S170 vmaxに対応する期間は拡張期であり、大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップ、即ち、
である、
S180 計算式
に基づきiFMR値を取得するステップ、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vhは最大充血状態における血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、好ましくは、本出願は再帰的アルゴリズムまたはバブルソートにより血流速度vの最大値を選び取り、k、cは影響パラメータk=1~3、c=0~10を表す、
S200 拡張期の微小血管抵抗指数がiFMR<Kである場合、調節パラメータをr=1とし、且つ、拡張期の微小血管抵抗指数がiFMR≧Kである場合、調節パラメータを
(但しKは100未満の正数であり、好ましくは、K=50とする)とするステップ。
本出願の一つの実施例において、S140は二種類の取得方法を含み、方法(1)は図5に示されるように、以下のステップを含む。
S141I 造影剤が冠状動脈入口である血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップ、
S142I 第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップ、
S143I v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度はそれぞれv1,..., vb,..., va,..., vN-1とするステップ。
S141I 造影剤が冠状動脈入口である血管段の最初の点まで流れた時の冠動脈造影画像を第1フレーム画像とし、造影剤が血管段の最後の点まで流れた時の冠動脈造影画像を第Nフレーム画像とするステップ、
S142I 第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの時間差及び中心線の長さの差を順に求め、時間差をそれぞれΔt1,..., Δtb,..., Δta,..., ΔtN-1とし、中心線の長さの差をそれぞれΔL1,..., ΔLb,..., ΔLa,..., ΔLN-1とするステップ、
S143I v=ΔL/Δtに基づき、但し、vは血流速度を表し、第Nフレーム画像から第N-1フレーム、...、第N-bフレーム、...、第N-aフレーム、...、第1フレーム画像までの血流速度をそれぞれ得て、血流速度はそれぞれv1,..., vb,..., va,..., vN-1とするステップ。
Claims (19)
- 拡張期の微小血管抵抗指数ユニットによって、血流速度v、大動脈圧波形、生理パラメータに基づき、拡張期の微小血管抵抗指数iFMRを取得するステップと、
調節パラメータユニットによって、前記拡張期の微小血管抵抗指数がiFMR<Kである場合、調節パラメータをr=1とし、且つ、
前記拡張期の微小血管抵抗指数がiFMR≧Kである場合、調節パラメータを
(但しKは100未満の正数である)とするステップと、
血流速度取得ユニットによって、式v'=rvhに基づき、補正後の最大充血状態における血流速度を取得するステップと、
を含み、
但し、v'は補正後の最大充血状態における血流速度を表し、vhは最大充血状態における血流速度を表す、
ことを特徴とする微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - vhは、
であり、
但し、vhは最大充血状態における血流速度を表し、
は一つの心周期区域内における平均血流速度を表し、zの値の範囲は1~3の定数であり、xの値の範囲は50~300の定数であり、K=50とする、
ことを特徴とする請求項1に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記拡張期の微小血管抵抗指数ユニットによって、血流速度v、大動脈圧波形、生理パラメータに基づき、最大充血状態における微小血管抵抗指数iFMRを取得するステップは、
前記血流速度vの最大値である拡張期の最大血流速度vmaxを選び取るステップと、
前記最大血流速度vmaxに対応する期間は拡張期であり、前記大動脈圧波形に基づき、大動脈拡張期の平均圧を取得するステップと、
を含み、
iFMRは、
であり、
は、
であり、
但し、
は大動脈拡張期の平均圧を表し、Pa1、Pa2、Pajはそれぞれ大動脈圧波形上の拡張期内の第1の点、第2の点、第jの点に対応する大動脈圧力値を表し、jは拡張期内の大動脈圧波形上に含まれる圧力点の数を表し、vhは最大充血状態における血流速度を表し、全ての血流速度vから最大値を選び取ることにより取得され、k、cは影響パラメータk=1~3、c=0~10を表す、
ことを特徴とする請求項1に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記影響パラメータはk=a×bであり、
但しaは糖尿病特徴値を表し、bは高血圧特徴値を表し、cは性別を表す、
ことを特徴とする請求項3に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 患者に糖尿病がない場合、0.5≦a≦1とされ、患者に糖尿病がある場合、1<a≦2とされ、
患者の血圧値が90mmHg以上である場合、1<b≦1.5とされ、患者の血圧値が90mmHg未満である場合、0.5≦b≦1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=3~10とされる、
ことを特徴とする請求項4に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 患者に糖尿病がない場合、a=1とされ、患者に糖尿病がある場合、a=2とされ、
患者の血圧値が90mmHg以上である場合、b=1.5とされ、患者の血圧値が90mmHg未満である場合、b=1とされ、
患者が男性である場合、c=0とされ、患者が女性である場合、c=5とされる、
ことを特徴とする請求項5に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記血流速度を取得するステップは、
画像読み取りユニットによって、少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るステップと、
血管段抽出ユニットによって、前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップと、
中心線抽出ユニットによって、前記血管段の中心線を抽出するステップと、
時間差ユニットによって、任意の2つのフレームの冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、中心線差ユニットによって、前記中心線の差を計算し、その差の値をΔLとするステップと、
血流速度計算ユニットによって、前記ΔLと前記Δtの比に基づき、血流速度vを求めるステップと、
を含む、
ことを特徴とする請求項1に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記血管段抽出ユニットによって、前記冠状動脈二次元造影画像組から、注目する血管段を抽出するステップは、
前記冠状動脈二次元造影画像組からNフレーム冠状動脈二次元造影画像を選び取るステップと、
前記冠状動脈二次元造影画像において、注目する前記血管の最初と最後の点をピックアップし、注目する前記血管段を取得するステップと、
を含む、
ことを特徴とする請求項7に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記中心線抽出ユニットによって、前記血管段の中心線を抽出するステップは、
血管骨格抽出ユニットによって、前記冠状動脈二次元造影画像から血管骨格を抽出するステップと、
前記血管段の延在方向、および2点間の最短経路を取る原則に基づくステップと、
前記血管骨格に沿って、前記血管段の中心線を抽出するステップと、
を含む、
ことを特徴とする請求項7に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記時間差ユニットによって、任意の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とするステップと、
を含む、
ことを特徴とする請求項7に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記時間差ユニットによって、任意の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フレーム画像までの血流速度をそれぞれ得るステップと、
を含む、
ことを特徴とする請求項7に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記血管段の中心線を抽出するステップの後、前記時間差ユニットによって、任意の2つのフレームの冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとし、且つ、前記中心線差ユニットによって、前記中心線の差を計算し、その差の値をΔLとするステップの前にさらに、
少なくとも2つの体位の冠状動脈二次元造影画像組を読み取るステップと、
前記血管段のジオメトリ構造情報を取得するステップと、
注目する前記血管段に対して図形処理を行うステップと、
前記血管段の血管輪郭線を抽出するステップと、
前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するステップと、
を含む、
ことを特徴とする請求項7に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 前記血流速度計算ユニットによって、前記ΔLと前記Δtの比に基づき、血流速度を求めるステップは、
前記三次元血管モデルに基づき、三次元血管モデルの中心線を取得し、前記冠状動脈二次元造影画像により抽出した中心線に対して補正を行い、且つ、中心線差分値ΔLに対して補正を行い、ΔL'を得るステップと、
前記ΔL'と前記Δtの比に基づき、血流速度vを求めるステップと、
を含む、
ことを特徴とする請求項12に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法。 - 請求項13に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法を含むことを特徴とする生理パラメータに基づき冠状動脈血管評定パラメータを取得する方法。
- 血流速度取得ユニット、大動脈圧波形取得ユニット、生理パラメータ取得ユニット、拡張期の微小血管抵抗指数ユニット及び調節パラメータユニットを備え、
前記拡張期の微小血管抵抗指数ユニットは前記血流速度取得ユニット、前記大動脈圧波形取得ユニット、前記生理パラメータ取得ユニットに接続され、
前記血流速度取得ユニットは、血流速度vを取得するためのものであり、
前記大動脈圧波形取得ユニットは、時間に伴い変わる大動脈圧波形をリアルタイムに取得するためのものであり、
前記生理パラメータ取得ユニットは、性別及び既往歴を含む患者の生理パラメータを取得するためのものであり、
前記拡張期の微小血管抵抗指数ユニットは、前記血流速度取得ユニット、前記大動脈圧波形取得ユニット、前記生理パラメータ取得ユニットが送信した血流速度v、 大動脈圧波形、生理パラメータを受信し、それから血流速度v、大動脈圧波形、生理パラメータに基づき、拡張期の微小血管抵抗指数iFMRを取得するためのものであり、
前記調節パラメータユニットは、前記拡張期の微小血管抵抗指数ユニットのiFMR値を受信し、前記拡張期の微小血管抵抗指数がiFMR<Kである場合、調節パラメータをr=1とし、且つ、前記拡張期の微小血管抵抗指数がiFMR≧Kである場合、調節パラメータを
(但しKは100未満の正数であり)とするためのものである、
ことを特徴とする前記請求項1~13のいずれかに記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法に用いられる微小血管抵抗指数に基づき最大充血状態における血流速度を調節する装置。 - さらに、
順に接続される画像読み取りユニット、血管段抽出ユニット、中心線抽出ユニット、前記画像読み取りユニットに接続される時間差ユニット、前記生理パラメータ取得ユニット、及び前記時間差ユニット、中心線差ユニットにそれぞれ接続される血流速度取得ユニットを備え、
前記中心線差ユニットは前記中心線抽出ユニットに接続され、
前記画像読み取りユニットは、少なくとも一つの体位の冠状動脈二次元造影画像組を読み取るためのものであり、
前記血管段抽出ユニットは、前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信し、前記画像中の注目する血管段を抽出するためのものであり、
前記中心線抽出ユニット、は前記血管段抽出ユニットが送信した血管段を受信し、前記血管段の中心線を抽出するためのものであり、
前記時間差ユニットは、前記画像読み取りユニットが送信した任意の2つのフレームの冠状動脈二次元造影画像を受信し、前記2つのフレームの冠状動脈二次元造影画像において造影剤が前記血管段を流れた時間の差を計算し、その差の値をΔtとするものであり、
前記中心線差ユニットは、前記中心線抽出ユニットが送信した2つのフレームの冠状動脈二次元造影画像において造影剤が血管段を流れた分段中心線を受信し、前記分段中心線の差を計算し、その差の値をΔLとするものであり、
前記血流速度取得ユニットは、血流速度計算モジュール、拡張期血流速度計算モジュールを備え、
前記血流速度計算モジュールは前記時間差ユニット、前記中心線差ユニットにそれぞれ接続され、
前記拡張期血流速度計算モジュールは前記血流速度計算モジュールに接続され、
血流速度計算モジュールは、前記時間差ユニット及び前記中心線差ユニットが送信した前記ΔLと前記Δtを受信し、前記ΔLと前記Δtの比に基づき、血流速度を求めるためのものであり、
前記拡張期血流速度計算モジュールは、前記血流速度計算モジュールが送信した血流速度値を受信し、前記血流速度の最大値である拡張期の血流速度を選び取るためのものであり、
前記生理パラメータ取得ユニットは、前記画像読み取りユニットの冠状動脈二次元造影画像を受信し、患者の生理パラメータ及び画像撮影角度を取得し、生理パラメータ及び画像撮影角度を前記拡張期の微小血管抵抗指数ユニットに伝達するためのものである、
ことを特徴とする請求項15に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する装置。 - さらに、
前記画像読み取りユニットにそれぞれ接続される血管骨格抽出ユニット、三次元血管再構成ユニット、及び前記血管骨格抽出ユニットに接続される輪郭線抽出ユニットを備え、
前記三次元血管再構成ユニットは前記生理パラメータ取得ユニット、前記中心線抽出ユニット、前記輪郭線抽出ユニットに接続され、
前記血管骨格抽出ユニットは、前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信し、前記冠状動脈二次元造影画像中の血管骨格を抽出するためのものであり、
前記輪郭線抽出ユニットは、前記血管骨格抽出ユニットの血管骨格を受信し、前記血管骨格に基づき、注目する前記血管段の輪郭線を抽出するためのものであり、
前記三次元血管再構成ユニットは、前記輪郭線抽出ユニット、前記生理パラメータ取得ユニット、前記中心線抽出ユニットが送信した輪郭線、画像撮影角度及び中心線を受信するため、且つ、前記画像読み取りユニットが送信した冠状動脈二次元造影画像を受信するため、且つ、前記血管段のジオメトリ構造情報に基づき、少なくとも2つの体位の血管の前記中心線、輪郭線を抽出した冠状動脈二次元造影画像を三次元平面上に投影し、三次元血管モデルを合成するためのものであり、
前記中心線抽出ユニットは、前記三次元血管再構成ユニットの三次元血管モデルから新たに前記血管段の中心線を抽出し、且つ新たに前記中心線の長さを取得するためのものである、
ことを特徴とする請求項16に記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する装置。 - 請求項15~17のいずれかに記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する装置を含むことを特徴とする冠状動脈分析システム。
- コンピュータプログラムがプロセッサにより実行される時に請求項1~13のいずれかに記載の微小血管抵抗指数に基づき最大充血状態における血流速度を調節する方法を実現することを特徴とするコンピュータ記憶媒体。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911065694.6A CN110786841B (zh) | 2019-11-04 | 2019-11-04 | 基于微循环阻力指数调节最大充血状态流速的方法及装置 |
CN201911065694.6 | 2019-11-04 | ||
PCT/CN2019/116674 WO2021087968A1 (zh) | 2019-11-04 | 2019-11-08 | 基于微循环阻力指数调节最大充血状态流速的方法及装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2022554298A JP2022554298A (ja) | 2022-12-28 |
JPWO2021087968A5 true JPWO2021087968A5 (ja) | 2023-07-18 |
JP7385318B2 JP7385318B2 (ja) | 2023-11-22 |
Family
ID=69442547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022525374A Active JP7385318B2 (ja) | 2019-11-04 | 2019-11-08 | 微小血管抵抗指数に基づき最大充血状態流速を調節する方法及び装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220254028A1 (ja) |
EP (1) | EP4056110B1 (ja) |
JP (1) | JP7385318B2 (ja) |
CN (1) | CN110786841B (ja) |
WO (1) | WO2021087968A1 (ja) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111227821B (zh) * | 2018-11-28 | 2022-02-11 | 苏州润迈德医疗科技有限公司 | 基于心肌血流量和ct图像的微循环阻力指数计算方法 |
WO2020201942A1 (en) | 2019-04-01 | 2020-10-08 | Cathworks Ltd. | Methods and apparatus for angiographic image selection |
US12039685B2 (en) | 2019-09-23 | 2024-07-16 | Cathworks Ltd. | Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device |
CN112164020B (zh) * | 2020-03-31 | 2024-01-23 | 苏州润迈德医疗科技有限公司 | 精确提取血管中心线的方法、装置、分析系统和存储介质 |
CN112690814B (zh) * | 2020-11-06 | 2022-10-14 | 杭州阿特瑞科技有限公司 | 一种低误差的冠状动脉血流储备分数测量方法 |
CN112971818B (zh) * | 2021-01-28 | 2022-10-04 | 杭州脉流科技有限公司 | 微循环阻力指数的获取方法、装置、计算机设备和存储介质 |
CN116206162B (zh) * | 2023-04-28 | 2023-08-01 | 杭州脉流科技有限公司 | 基于造影影像的冠脉血流储备获取方法、装置及设备 |
CN118013233B (zh) * | 2024-04-07 | 2024-08-02 | 中汽研汽车检验中心(天津)有限公司 | 一种重型商用车预见性巡航测试工况构建方法 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1559345A (zh) * | 2004-02-20 | 2005-01-05 | 无锡贝尔森影像技术有限公司 | 脑血管血液动力学检测方法和仪器 |
US10758125B2 (en) * | 2016-04-29 | 2020-09-01 | Siemens Healthcare Gmbh | Enhanced personalized evaluation of coronary artery disease using an integration of multiple medical imaging techniques |
CA3037435A1 (en) * | 2016-09-20 | 2018-03-29 | Heartflow, Inc. | Systems and methods for monitoring and updating blood flow calculations with user-specific anatomic and physiologic sensor data |
CN107689032A (zh) * | 2017-07-05 | 2018-02-13 | 北京工业大学 | 一种个性化无创计算患者最大充血状态下冠脉分支血流量的方法 |
CN107730540B (zh) | 2017-10-09 | 2020-11-17 | 全景恒升(北京)科学技术有限公司 | 基于高精度匹配模型的冠脉参数的计算方法 |
CN108305246B (zh) * | 2017-11-15 | 2020-10-09 | 深圳科亚医疗科技有限公司 | 基于医学图像预测血流特征的装置及系统 |
CN107978371B (zh) * | 2017-11-30 | 2021-04-02 | 博动医学影像科技(上海)有限公司 | 快速计算微循环阻力的方法及系统 |
CN108550189A (zh) | 2018-05-03 | 2018-09-18 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
CN108742587B (zh) * | 2018-06-20 | 2021-04-27 | 博动医学影像科技(上海)有限公司 | 基于病史信息获取血流特征值的方法及装置 |
CN109009037B (zh) * | 2018-06-20 | 2021-11-19 | 博动医学影像科技(上海)有限公司 | 基于个体性别获取血管压力差的计算方法及装置 |
CN108742570B (zh) * | 2018-06-20 | 2021-10-01 | 博动医学影像科技(上海)有限公司 | 基于冠脉优势类型获取血管压力差的装置 |
CN110367965B (zh) * | 2018-09-19 | 2022-03-08 | 苏州润迈德医疗科技有限公司 | 便捷测量冠状动脉血管评定参数的方法、装置及系统 |
CN110384493A (zh) | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | 测量微循环阻力指数的系统以及冠脉分析系统 |
CN110384494A (zh) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | 测量微循环阻力指数的方法 |
CN109770930B (zh) * | 2019-01-29 | 2021-03-09 | 浙江大学 | 一种冠状动脉微循环阻力的确定方法和装置 |
-
2019
- 2019-11-04 CN CN201911065694.6A patent/CN110786841B/zh active Active
- 2019-11-08 EP EP19952100.6A patent/EP4056110B1/en active Active
- 2019-11-08 JP JP2022525374A patent/JP7385318B2/ja active Active
- 2019-11-08 WO PCT/CN2019/116674 patent/WO2021087968A1/zh unknown
-
2022
- 2022-04-28 US US17/731,607 patent/US20220254028A1/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110786841B (zh) | 基于微循环阻力指数调节最大充血状态流速的方法及装置 | |
EP3226758B1 (en) | Method, apparatus and computer program for determining a blood pressure value | |
CN108665449B (zh) | 预测血流矢量路径上的血流特征的深度学习模型及装置 | |
CN109770888A (zh) | 基于压力传感器和造影图像计算瞬时无波形比率的方法 | |
WO2012073863A1 (ja) | 超音波撮像装置、超音波撮像方法、超音波撮像プログラム | |
EP2967368B1 (en) | Apparatus and methods for computing cardiac output of a living subject via applanation tonometry | |
EP3302285B1 (en) | Apparatus and method for determining blood flow velocity | |
CN110786840B (zh) | 基于生理参数获取血管评定参数的方法、装置及存储介质 | |
WO2014084367A1 (ja) | 医用画像処理装置 | |
CN103976760B (zh) | 左心房压力测量方法以及左心房压力测量装置 | |
US20230210492A1 (en) | Hemodynamii parameter estimation | |
CN111833343A (zh) | 冠状动脉狭窄度估计方法系统及设备 | |
WO2013057999A1 (ja) | 超音波撮像装置及び方法 | |
JPWO2021087968A5 (ja) | ||
EP4185188B1 (en) | Method and device for determining a coronary microvascular resistance score | |
CN114664455A (zh) | 一种冠状动脉血流储备分数计算方法及装置 | |
CN116524134B (zh) | 基于血管内图像的三维血管建模方法及ffr计算方法和系统 | |
JPWO2021087967A5 (ja) | ||
CN110929604A (zh) | 基于造影图像的流速的筛选方法、装置、系统和存储介质 | |
GB2570131A (en) | Fluid flow analysis | |
EP4131154A1 (en) | Coronary artery narrowing detection based on patient imaging and 3d deep learning | |
EP4378394A1 (en) | Hemodynamic parameter estimation | |
CN118695816A (zh) | 血液动力学参数估计 | |
Guala et al. | P47 Abnormal Flow Pattern in Marfan Patients is Related to Aortic Geometric Features: A 4D Flow MRI Study | |
Pasynkov et al. | A Method for Automatic Calculation of the Diameter of a Pulsating Blood Vessel on Ultrasound Images in Video Stream Mode |