JP2024020363A - ハイブリッドな画像-侵襲性-圧力血行動態機能評価 - Google Patents

ハイブリッドな画像-侵襲性-圧力血行動態機能評価 Download PDF

Info

Publication number
JP2024020363A
JP2024020363A JP2023192916A JP2023192916A JP2024020363A JP 2024020363 A JP2024020363 A JP 2024020363A JP 2023192916 A JP2023192916 A JP 2023192916A JP 2023192916 A JP2023192916 A JP 2023192916A JP 2024020363 A JP2024020363 A JP 2024020363A
Authority
JP
Japan
Prior art keywords
blood flow
artery
pressure
medical image
pressure measurement
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.)
Pending
Application number
JP2023192916A
Other languages
English (en)
Inventor
べレヴィル クラウド
Belleville Claude
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Opsens Inc
Original Assignee
Opsens Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Opsens Inc filed Critical Opsens Inc
Publication of JP2024020363A publication Critical patent/JP2024020363A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/0215Measuring pressure in heart or blood vessels by means inserted into the body
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • A61B5/0275Measuring blood flow using tracers, e.g. dye dilution
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4887Locating particular structures in or on the body
    • A61B5/489Blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6846Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive
    • A61B5/6847Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive mounted on an invasive device
    • A61B5/6851Guide wires
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Vascular Medicine (AREA)
  • Hematology (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

【課題】患者特有の血行動態パラメータを計算する方法を記載する。【解決手段】該方法は、血管内圧測定デバイスを用いて動脈における少なくとも1つの圧力測定値を測定するステップと、及び医用撮像機器から動脈の少なくとも1つの医用画像を撮像するステップとを備え、動脈の少なくとも1つの医用画像は少なくとも1つの圧力測定値と同期するものとする。圧力測定値及び医用画像の双方をコンピュータシステムに供給して、少なくとも1つの医用画像から血流を計算し、少なくとも2つの動脈圧力低下及びそれらに対応する血流成分から動脈のパラメータを計算し、また血流及び動脈のパラメータに基づいて1つ又は複数の患者特有の血行動態パラメータを計算する。【選択図】図3

Description

本出願は、参照により全体が本明細書に組み入れられるものとする、2018年6月
27日出願の米国仮出願第62/690,756号の優先権を主張する。
(a) 分野
本明細書で開示する本発明要旨は、概して、血行動態機能を評価する方法に関する。よ
り具体的には、血管内圧測定デバイスを医用画像撮像と組み合わせて使用して血行動態パ
ラメータを決定する方法に関する。
(b) 関連従来技術
冠動脈疾患は、現在では血流に対する冠動脈狭窄の影響を評価することによって診断さ
れる。臨床研究の大多数は、血流予備量比(FFR:Fractional Flow Reserve)の測定
に基づく冠動脈狭窄の治療が患者予後の改善に導くことを示した。FFRは遠位圧力(P
d)対近位圧力(Pa)の比をとることにある。冠血流制限病変がない場合、FFRは1
に等しい。血流制限が増大すると、FFR値は低下する。FFR比は、何ら病変がない同
一動脈と比べた病変がある血流予備量比のパーセンテージで表現する。より具体的には、
0.80より高いFFR値の特定狭窄部位の治療は保留するとともに、FFR値が0.80
以下である場合、その狭窄部位を治療することが患者予後の改善に導くことを示した。
FFRは心外膜狭窄の重症度を測定するが、微小血管機能不全の存在はFFR値を増大
させて、病変の重症度を過小評価するおそれがある。冠血流予備能(CFR:Coronary f
low reserve)は、心外膜抵抗及び微小血管抵抗の双方に依存する。したがって、FFR
とCFRとの間における不一致は、病変部位を治療すべきか否かを決断する際に起こり得
る。このFFRとCFRとの間における不一致は微小血管機能を考慮する正確な診断のた
めの必要性を示している(非特許文献1:Bon-Kwon Koo, JACC, Vol. 10, No. 10, 2017
参照)。急性心筋梗塞の最も多くの場合に正常な心外膜血流を回復するにも係わらず、高
い割合の微小血管障害及び心筋壊死が残ることになる(非特許文献2:Niccoli G et al,
J Am Coll Cardiol 2009; 54:281-292参照)。選択的冠動脈血管形成術は、術中の20
%にも達する心筋梗塞に関連し、微小塞栓術の結果であると考えられる(非特許文献3:
Jaffe R et al., JACC Cardiovasc Interv 2010; 3:695-704参照)。安定狭心症の症状を
有する患者の大部分は正常な冠動脈を呈し、このため微小血管機能の修正は心筋虚血に関
与する要因である(非特許文献4:Schwartz L et al., Arch Intern Med 2001; 161:182
5-1833参照)。肥大性心筋症を有する患者において、微小血管機能不全は、左心室機能不
全及び死亡の強い予測因子である(非特許文献5:Cecchi F et al, N Engl J Med 2003;
349:1027-1035参照)。微小血管抵抗は、ST-セグメント-上昇後の左心室及び臨床転帰
に関連することが示された。虚血時間、ST-セグメント-上昇、血管造影紅潮等級、及び
CFRを含む心筋再かん流有効性における標準臨床的尺度と比較すると、微小血管抵抗は
、重症度分類に関する優れた臨床的価値を有していることを示した、また不全な心筋再か
ん流に対する参考試験とみなすことができる(非特許文献6:Carrick D. et al., Circu
lation 2016, 134:1833-1847参照)。FFRではなく微小血管抵抗は、グレーゾーンFF
R病変部における冠血流の改善に関する最強の予測因子である(非特許文献7:Niida T.
et al, Catheter Cardiovasc Interv. 2018;1-11参照)。
近年、冠血流に対する理解は主に心外膜動脈に集中されてきたが、微小循環は虚血患者
を診断する上で重要な役割を果たすことを示す証拠実体となりつつあり、究極的には最も
適切な心筋再かん流治療プランになり得る。
しかし、微小血管抵抗の測定は、今日のカテーテル検査室(又はカテラボ)において慣
行ではなく、これはすなわち、その方法は巧妙さを要し、また臨床的ワークフローが煩雑
になるからである。微小血管抵抗は、一般的には、遠位圧力(Pd)を動脈内における血
流で除算することによって計算する。遠位圧力は、通常は感圧ガイドワイヤを用いて測定
する。血流を測定するには種々の方法が存在する。血液速度を測定するためには、ドップ
ラー感知ガイドワイヤを使用することができる。この方法は、ワイヤ取扱い特性が悪いこ
とに起因して、またドップラーセンサを有するワイヤ部分のオペレータ依存及び解剖学的
構造依存である指向方向で血液速度測定が変化することに起因して次善策である。他の方
法は、血液が動脈の遠位部分に到達するに要する時間に基づき、この場合、生理食塩水の
ボーラス投入時点と遠位に配置した温度感知ガイドワイヤによる検出時点との時間差を測
定する、いわゆる熱希釈法を行う。熱希釈方法は、可変通過時間測定を招いて複数回の生
理食塩水ボーラス投入を必要とする。さらに、患者間で動脈容積が一定であると仮定し、
この場合最小冠動脈疾患であることをも仮定する。近年、新しい方法は、室温の生理食塩
水の連続既知流量注入後における血液平均温度低下測定に基づく(非特許文献8:Aarnou
dse W. et al., JACC, Vol. 50, N0. 24, 2007参照)。この最近の方法はより正確かつ再
現可能性があると考えられるが、処置時間及びコストを追加させがちな特別カテーテル及
びセットアップを使用する必要がある。ドップラーをベースとする方法を除いて、これら
方法は、平均血流を測定し、またひいては時間依存血流及び時間依存微小血管抵抗を測定
する能力を持たない(すなわち、血流及び抵抗が一定である代わりに時間の関数である場
合に適用不可である)。
近年において、侵襲性圧力測定の必要がなく血流予備量比を決定することを考慮して、
冠動脈の非侵襲性画像をベースとする圧力-血流特性を開発する多数の業績があった。血
管造影、CTスキャン、MRIスキャン、及び従来既知の他の非侵襲性撮像システムを用
いて測定された患者特有の動脈ジオメトリは、これら非侵襲性患者血行動態評価の中心で
ある。これら方法は、とりわけ動脈、微小血管系及び心臓の寄与を含む閉ループモデルを
構築するステップを伴う。心臓及び微小血管系は、抵抗、インダクタンス及びキャパシタ
ンスの集中したシステムでモデル化し、心臓は血流駆動力も含む。他方、動脈は、3次元
(3D)計算モデル、低減次数計算モデル、動脈セグメントの特定形体に基づく異なるタ
イプのモデルを組み合わせたマルチスケールモデルのいずれかを用い、通常周知のナビエ
・ストークス方程式を使用してモデル化する。動脈は、さらに、非侵襲性撮像システムか
ら取得した動脈の患者特有ジオメトリを用いて動脈、又は動脈の所定セグメントにおける
抵抗、インダクタンス及び/又はキャパシタンスを計算する解析若しくは半解析的モデル
又は集中型モデルを使用してモデル化することもできる。
閉ループ系を解明する、ひいては動脈の圧力-血流特性を決定するためには、その系の
境界及び初期条件を知得する必要がある。これら境界及び初期条件は評価対象患者に対し
て不明であるため、これらモデルは、通常すべて母集団の特性からのものである、平均微
小血管抵抗、血管柔軟性、左心室収縮性のような汎用パラメータを使用する。これら仮定
は、動脈の圧力-血流特性を計算する上で大きな誤差を生ずることになり得る。圧力-血流
特性を計算するのに必要なすべての入力パラメータのうち、微小血管抵抗はすべての誤差
の59%にもなる(非特許文献9:Morris et al., JACC. Vol. 2, No. 4, 2017参照)。
患者間における微小血管抵抗変動性は、実際に、閉ループモデルから計算されるとき、臨
床結果に関して上述した通り観測されたように、圧力-血流特性の精度に直接的な影響を
及ぼす。したがって、追加の患者特有測定値を使用することなく患者特有微小血管抵抗を
決定することはできない。
「生理学的測定値を使用して画像ベースの血流計算を向上させる方法及びシステム(Me
thod and system for enhancing image-based blood flow computations using physiolo
gical measurements)」と題した特許文献1(米国公開特許公報第2017/00320
97号)は、微小血管抵抗を計算するため種々の測定値を使用する。第1方法は、血流に
等しいと仮定される投入造影剤が動脈の遠位部分に到達するための総移行時間(Tmn
の測定ステップと、微小血管抵抗を得るため、このTmn値に測定遠位圧力(Pd)を乗
算するステップとを含む。この方法は、患者特有でない一般平均動脈容積を仮定し、した
がって、動脈における患者疾患を考慮していない。さらに、微小血管抵抗に基づく時間を
何ら考慮していない。提案された他の方法は、種々の圧力-血流関係に侵襲性圧力又は血
流測定値を付加すること、例えば、3次元計算流動力学(CFD:Computation Flow Dyn
amic)、低減次数CFD、動脈内の形体に基づく異なる圧力-血流関係を組み合わせるマ
ルチスケール方法、解析若しくは半解析的方法、集中型モデル(抵抗、インダクタンス及
びキャパシタンスモデル)に基づく。圧力又は血流測定値は、モデル精度を改善すること
に充てられる。しかし、これら方法は、冠動脈内における侵襲性デバイスの存在をフルに
活用することはなく、これはすなわち、依然として圧力-血流/動脈ジオメトリ関係に大
きな程度依存するからである。
米国公開特許公報第2017/0032097号明細書
Bon-Kwon Koo, JACC, Vol. 10, No. 10, 2017 Niccoli G et al, J Am Coll Cardiol 2009; 54:281-292 Jaffe R et al., JACC Cardiovasc Interv 2010; 3:695-704 Schwartz L et al., Arch Intern Med 2001; 161:1825-1833 Cecchi F et al, N Engl J Med 2003; 349:1027-1035 Carrick D. et al., Circulation 2016, 134:1833-1847 Niida T. et al, Catheter Cardiovasc Interv. 2018;1-11 Aarnoudse W. et al., JACC, Vol. 50, N0. 24, 2007 Morris et al., JACC. Vol. 2, No. 4, 2017
したがって、微小血管抵抗を信頼性高く測定するとともに、今日の臨床的ワークフロー
に適合するよう運用が迅速かつ容易であるシステムに対する必要性が存在する。
本発明の態様によれば、
・血管内圧力測定デバイスを用いて動脈内の少なくとも1つの圧力測定値を測定するス
テップと、
・医用撮像機器から前記少なくとも1つの圧力測定値と同期する、動脈の少なくとも1
つの医用画像を撮像するステップと、
・前記少なくとも1つの圧力測定値及び前記少なくとも1つの医用画像の双方を計算シス
テムに供給するステップと、
・前記少なくとも1つの医用画像から血流を計算するステップと、
・少なくとも2つの動脈圧力低下及びこれに対応する血流成分から動脈のパラメータを計
算するステップと、及び
・前記血流及び前記動脈のパラメータに基づいて、患者特有血行動態パラメータを計算す
るステップと
を備える、方法を提供する。
一実施形態によれば、前記少なくとも1つの圧力測定値を測定するステップ及び前記少
なくとも1つの医用画像を撮像するステップは、血流の2つの条件で実施し、前記動脈の
パラメータを計算するステップは、前記血流の2つの条件で測定された前記少なくとも2
つの動脈圧力低下及びこれに対応する血流成分を用いて、動脈のパラメータの二次方程式
を解くステップを含む。
一実施形態によれば、前記血流の2つの条件は、より多い血流条件及びより少ない血流
条件を含み、方法は、さらに、より多い血流条件にする薬剤を投与するステップを備え、
この後、より多い血流条件の下で前記少なくとも1つの圧力測定値を測定するステップ及
び前記少なくとも1つの医用画像を撮像するステップを行う。
一実施形態によれば、前記より多い血流条件の誘発は、部分充血又は満充血を誘発する
造影剤又は充血剤を注入することによって行う。
一実施形態によれば、本発明方法は、さらに、前記少なくとも1つの圧力測定値を測定
するステップを微小血管膨張剤の注入と同期調整するために、前記より多い血流条件の誘
発時に音声信号を発生するステップを備える。
一実施形態によれば、本発明方法は、さらに、前記動脈内に放射線吸収造影媒剤を導入
するステップを備え、前記少なくとも1つの医用画像を撮像するステップは、さらに、
前記少なくとも1つの医用画像において、前記動脈の長さに沿って変化する動脈直径を
測定するステップと、
前記動脈における前記放射線吸収造影媒剤の伝播を追跡し、また前記動脈の第1ポイン
トL1から第2ポイントL2に伝播するに要する時間を測定するステップと、
を有し、
前記血流計算ステップは、動脈内の平均血流を計算するため、次式、すなわち、
の容積Vを前記伝播に要する時間で除算するステップを含む。
一実施形態によれば、平面を画定する前記少なくとも1つの医用画像で、前記少なくと
も1つの医用画像の前記平面に直交する動脈直径を測定するステップは、前記放射線吸収
造影媒剤に対する濃度測定法を使用して、前記動脈断面における測定した明暗度に対しラ
ンベルト・ベールの法則を適用することによってxの関数としてDを測定するステップを
含む。
一実施形態によれば、本発明方法は、さらに、前記動脈内への充血剤の冠動脈内又は経
皮的注入を使用するステップと、及び前記動脈のパラメータを使用して前記充血条件に関
連する血流を計算するステップと、を備え、前記患者特有血行動態パラメータを計算する
ステップは、圧力低下及び前記計算した前記充血条件に関連する血流に基づく。
一実施形態によれば、前記患者特有血行動態パラメータは微小血管抵抗を含む。
一実施形態によれば、前記少なくとも1つの圧力測定値を測定するステップ及び前記少
なくとも1つの医用画像を撮像するステップは、安静条件の下で、第1に少なくとも収縮
期の一部分をカバーする収縮期中に、また第2に少なくとも拡張期の一部分をカバーする
拡張期中に実施し、前記動脈のパラメータを計算するステップは、前記収縮期及び前記拡
張期にそれぞれ測定及び計算された前記少なくとも2つの動脈圧力低下及びこれらに対応
する血流成分を使用して前記動脈のパラメータの二次方程式を解くステップを含む。
一実施形態によれば、本発明方法は、さらに、前記動脈内で狭窄部の存在を同定するス
テップと、前記狭窄部の遠位側にある何ら狭窄部がないセグメントを同定するステップと
、及び前記セグメントに沿って複数の圧力測定値を取得するステップとを備える。
一実施形態によれば、前記狭窄部の存在を同定するステップは、前記動脈のパラメータ
のうち1つのパラメータを、前記動脈内で狭窄部の存在を同定するための所定閾値と比較
するステップを含む。
一実施形態によれば、前記狭窄部の存在を同定するステップは、前記少なくとも1つの
医用画像を使用して、50%以上の内腔閉塞の存在を決定することによって前記動脈内で
の前記狭窄部の存在を同定するステップを含む。
一実施形態によれば、前記セグメントに沿って複数の圧力測定値を取得するステップは
、血管内圧力測定デバイスの先端を前記セグメントの最も遠位の位置に送り、また前記血
管内圧力測定デバイスを引き戻すことによって行う。
一実施形態によれば、本発明方法は、さらに、前記セグメントに沿う前記複数の圧力測
定値を使用し、またナビエ・ストークス方程式を数値的に解いてジオメトリベースの血流
を計算するステップと、並びに前記動脈のパラメータ及び前記少なくとも1つの医用画像
から血流を計算するステップに対する補正に適用するため、したがって、前記狭窄部の存
在を考慮する前記患者特有血行動態パラメータに対する補正を適用するため、前記ジオメ
トリベースの血流を使用するステップと、を備える。
一実施形態によれば、前記少なくとも1つの圧力測定値を測定するステップ及び前記少
なくとも1つの医用画像を撮像するステップは、前記少なくとも1つの圧力測定値及び前
記少なくとも1つの医用画像をリアルタイムで同時に取得することによって、又は前記少
なくとも1つの圧力測定値及び前記少なくとも1つの医用画像にタイムスタンプ付けする
ことによって、同期させる。
本発明の他の態様によれば、血管内圧力測定デバイス及び医用撮像機器と通信する計算
システムであり、命令及び前記血管内圧力測定デバイス及び医用撮像機器からのデータを
記憶するメモリと、及び前記命令を実行するプロセッサであって、
・動脈内における前記血管内圧力測定デバイスから少なくとも1つの圧力測定値を受け
取る、
・他の圧力デバイスから少なくとも1つの大動脈圧力を受け取る、
・前記医用撮像機器から動脈の少なくとも1つの医用画像を受け取る、
・前記少なくとも1つの医用画像から血流を計算する、
・少なくとも2つの動脈圧力低下及びこれらに対応する血流成分から前記動脈のパラメ
ータを計算する、及び
・前記血流及び前記動脈のパラメータに基づいて、患者特有血行動態パラメータを計算
する
ための命令を実行する、該プロセッサと、
を備える、システムを提供する。
一実施形態によれば、前記血管内圧力測定デバイスは圧力ガイドワイヤである。
本発明開示の他の特徴及び利点は、添付図面と組み合わせた以下の詳細説明から明らか
になるであろう。
実施形態による、動脈内の圧力測定値を取得する血管内圧測定デバイスを示す概略的断面図である。 実施形態による、モデル化された動脈を示す概略的断面図である。 実施形態による、拡張期血流測定値及び充血性媒剤で血行動態パラメータを決定する方法を示すフローチャートである。 実施形態による、収縮期及び拡張期の血流測定値で、またいかなる充血性媒剤もない状態で血行動態パラメータを決定する方法を示すフローチャートである。 実施形態による、収縮期及び拡張期の血流測定値で、またいかなる充血性媒剤もなく幾何補正した状態で血行動態パラメータを決定する方法を示すフローチャートである。 実施形態による、分岐を有する動脈内の圧力測定値を取得する血管内圧測定デバイスを示す概略的断面図である。 実施形態による、異なるセグメントに異なる狭窄部を離隔させて示すセグメント化した動脈内の圧力測定値を取得する血管内圧測定デバイスを示す概略的断面図である。
添付図面にわたり同様の形体は同様の参照符号で示すことに留意されたい。
実施形態による本発明方法は、血管内圧力測定デバイス(例えば、ガイドワイヤ)のよ
うな侵襲性圧力測定機器と、動脈の医用画像を撮像する医用撮像機器のような非侵襲性機
器との同時(又は同期)組合せを伴う。この方法は、血流の適正測定に依存する他の値を
決定するよう、より精密に血流を決定するためのものである。
最近の従来技術による方法において一般的であるような侵襲性機器に変えて非侵襲性機
器を使用する、またこのことは上述したような欠点を伴う非侵襲性機器使用の代わりに、
本発明による方法は双方ともに使用する。この組合せにより、血管内圧力測定デバイスか
らの精密圧力測定を可能にし、また動脈のジオメトリ及び動脈内の血流データに対してリ
アルデータを有する血行動態パラメータの計算をリアルタイムで個人特定化する医用画像
で補完する。このようにして、提案された非侵襲性機器からジオメトリ又は血流データが
収集されなかった場合よりも、血行動態パラメータは、より高い精度で計算され、また動
脈内における実際のリアルタイム条件をよりよく反映する。
本明細書記載の実施形態の目的は、臨床的ワークフローに適合した方法を提供すること
にあり、すなわち、標的動脈内における血流を測定又は推定できる、また血管内圧力測定
とともに、微小血管抵抗、冠血流予備能(CFR)及び血流予備量比(FFR)のような
血行動態パラメータを計算できる迅速かつ容易な方法を得るにある。さらに、正確な冠血
流を取得するよう侵襲性圧力測定を冠動脈画像に組み合わせることも目的とし、これによ
り、血管内圧力測定と組み合わせるとき血行動態パラメータの計算を可能にする。
図1は、冠動脈疾患(CAD:coronary artery disease)を有する患者の血行動態機
能を評価するのに使用される代表的なセットアップを示す。心外膜動脈1は、動脈開口部
4から遠位端部5に血流を導き、この遠位端部5で血液は酸素交換を生ずる微小血管系に
進入する。微小血管抵抗は、ここに単一抵抗Rで示す。心筋キャパシタンスのような他
の生理学的機能を考慮に入れる微小血管系をモデル化することも可能であることは後で示
す。代表的セットアップにおいて、案内カテーテル3を動脈開口部に配置する。案内カテ
ーテルを使用して、動脈輪郭を明らかにするのに用いられる造影剤を送給し、圧力ガイド
ワイヤ2のような医用デバイスを動脈開口部及び冠動脈遠位に案内し、本明細書で大動脈
圧Paと称される動脈入口の血圧を測定する。大動脈圧は、血液を動脈内に強制的に流す
境界条件として使用することができる。圧力ガイドワイヤによって測定された遠位圧力P
dも境界条件として使用することができる。
図2は、血管系の機能的等価体を示す。開口部における血流及び圧力(近位圧力)はそ
れぞれQ及びPaで表すことができる。心外膜は、圧力低下Δpを生じ、動脈の遠位端部
における遠位圧力Pdをもたらす狭窄部を含むことがあり得る。開口部と遠位端部との間
に大きな分岐がないと仮定する場合、血流Qは一定になる。側方血管(側枝とも称される
)も若干の患者に存在し得る。側枝は、通常圧力ガイドワイヤ測定ポイントPdに対して
遠位位置、したがって、Pdポイントと微小血管系との間における主要血管に接続する。
側枝は微小循環系への付加的血流に寄与し、したがって、狭窄部にわたる血流及び圧力低
下の双方を減少する。側枝からの血流は心外膜血流に付加するが、狭窄部前後の遠位圧力
は変化せず、このことは、FFR測定における側枝血流寄与を考慮に入れることを可能に
する。FFRは以下の関係式1によって定義され、ここでQは狭窄部(通常は動脈)が
ない場合の血流であり、Qは狭窄部がある場合、応力条件の下、すなわち、満充血条件
(通常はアデノシンのような負荷誘導剤によって誘発される)の下で得られる場合におけ
る双方の血流である。血流は圧力の一次関数であり、微小血管抵抗Rはこれら関係量間
における比例パラメータ、この場合、血流は微小血管圧力差及びこれに関連する抵抗R
であると仮定するのは合理的である。微小血管抵抗Rを考慮することは最小限であり、
狭窄部があるか否かに変化せず、また静脈圧力(Pv)を考慮することは無視しうるもの
(Pv≒0)であり、FFRは関係式1の最後の部分から得られる。
圧力ガイドワイヤから得られる遠位圧力測定値は、画像上で位置とともに圧力を同時記
録することによって動脈の遠位位置に関連付けることができる。同時記録手段は、双方の
値を同時に測定及び記録することを意味する。さらに、動脈に沿って1つより多い圧力測
定値を取得し、画像における1つより多い圧力及び関連位置を同時記録することができる
。さらにまた、画像の記録及び動脈に沿って測定した圧力とともに、圧力測定値各々に関
連する画像上の位置を同時記録しながら、圧力ガイドワイヤを引き戻することができる。
圧力センサがどこにあるか、及びおそらく圧力ガイドワイヤ2の先端部に設けられる放射
線不透過性先端部に起因してX線不透過像がどこで急激に変化するかを、圧力ガイドワイ
ヤ2の先端部端縁を検出することによって或る圧力測定値は動脈内で局限することができ
る。
従来技術のシステムは、精度又は容易性のいずれかを改善するよう侵襲性測定を画像と
組み合わせることに傾注していたものであり、血行動態機能の計算は取得した画像のジオ
メトリに強く依存しており、3-D計算血行動態(CFD)、低減次数CFD、0(ゼロ
)-D解析モデル、又は他の集中型モデルから動脈の圧力-血流特性を計算することによっ
て行うものであった。これらすべてのモデルは、圧力-血流特性をモデル化するよう動脈
のジオメトリに基づく。したがって、この実施形態の目的は、測定機器が行う最終的評価
を改善するよう動脈のジオメトリに影響を受けにくい圧力-血流特性を開発することにあ
る。
以下の実施形態は、データを取得、保存、処理、表示、送信するために使用される、又
は使用及び測定したデータを後で示される、少なくとも1つのコンピュータを備えること
が理解される。コンピュータは、処理ユニット、メモリ、必要な通信ポートを備え、信号
上で論理演算を実施できるようにする任意なシステムであると理解されたい。
平均血流測定
第1実施形態は、3次元画像再構築する必要なく、動脈の局所的ジオメトリに対する依
存度が少ない動脈特有血行動態パラメータを測定するよう提案された。放射線吸収造影剤
を動脈開口部に導入する。造影剤の伝播を追跡し、またひいては血液が1つのポイントか
ら他のポイントに伝播するのにかかる時間を決定することができる。
この造影剤伝播を追跡するための様々な従来技術の方法が存在し、これには「デジタル
血管造影像から瞬間及び平均血流量を決定する方法及びシステム(Method and system fo
r determination of instantaneous and average blood flow rates from digital angio
grams)」と題する米国特許第5,150,292号、及び「脈管系内の液体流れを特徴付
けている量を決定するための方法及びデバイス(Method of and device for determining
quantities characterizing the flow of a liquid in a vascular system)」と題する
米国特許第5,150,292号に記載の方法があり、これら双方の特許文献は参照により
本明細書に組み入れられるものとする。参照により本明細書に組み入れられるものとする
特許文献1(米国公開特許公報第2017/0032097号)は、血液が動脈の遠位部
分に到達するのにかかる時間を測定するための類似方法を使用する。その従来技術文献に
おける血流は、このとき、血液が動脈の遠位端部に到達する時間の逆数をとる(Q~1/
min)熱希釈法によって血流を測定するのに使用されるのと同一の関係を用いて計算
される。特定標的動脈の容積は母集団の平均動脈容積と同一であると仮定し、したがって
、それは患者特有又は動脈特有のいずれでもない。温感ガイドワイヤを使用する熱希釈法
は画像ベースのシステムではないため、探査下にある患者動脈の特有容積部へのアクセス
が容易ではない。
本発明実施形態による方法は、むしろ標的動脈容積部を血流計算に結び付ける、及びひ
いては血行動態パラメータ計算に結び付けるステップを備える。動脈容積部は、以下の関
係式2を解くことによって得ることができる。Lは造影剤追跡がスタートする初期位置
、代表的には動脈の開口部であり、またLは圧力ガイドワイヤの感圧部分が配置される
位置である。D(x)は、周知の端縁検出方法の使用により動脈輪郭を検出することによ
って測定される動脈直径である。この直径は動脈の長さに沿って、すなわち、動脈の中心
線に沿って積分する。
体積流量Qmnは、動脈の容積部Vを造影剤がLからLまで移動するに要する時間
で除算することによって計算する。微小血管抵抗Rのような血行動態パラメータは、侵
襲性遠位圧力測定値Pdを体積流量Qmnで除算することによって計算する。冠血流予備
能(CFR)も、充血部で測定した血流を安静時(充血がない時)に測定した血流で除算
することによって計算することができる。
関係式2は、動脈が円形であると仮定するとともに、管腔がほぼ円形ではないことがあ
り得る疾患動脈は想定しない。目視に基づいて、非円形性は動脈断面の過大評価又は過小
評価に至るおそれがある。計算した総動脈容積部に対する影響は、過大評価及び過小評価
は動脈長さにわたり相殺されるという事実によってかなり最小化される。局所的特有ジオ
メトリに極めて鋭敏なジオメトリベースのCFDで導き出した任意な圧力-血流特性に対
抗して、本発明方法は、全体ジオメトリに基づき、またひいては局所的ジオメトリ不正確
さに対してはるかに鋭敏ではない。
この方法の非円形ジオメトリに対する相対的不感性にも係わらず、造影剤密度を用いて
、すなわちデンシトメトリー法の使用により管腔面積精度を改善することができる。造影
剤による放射線吸収は放射線が通過する吸収媒体の深さに依存すること(ランベルト・ベ
ールの法則)は当業者にとって既知である。したがって、断面積で生ずる吸収を測定する
ことによって動脈の垂直方向成分を推定することができる。
吸収媒体を透過する放射線レベルは以下の関係式5によって表現され、Iは入射放射
線強度であり、Iは透過放射線強度であり、thは吸収媒体の厚さ、すなわち見ている
平面に直交する動脈深さであり、またαは造影剤吸収係数に対応する定数である。円形で
ある、例えば、疾患がないと仮定される既知な同一動脈又は他の動脈における断面でI
-calを測定することによってαを推定できる。同様に、αを得るために案内カテーテ
ルを使用することもできる。Ii-calは、IT-calの位置近傍であるが、動脈断
面外側の値をとることによって得るとともに、thcalは画像平面における動脈直径で
ある。αは関係式6を用いて計算する。αを知ることによって、関係式7の使用により動
脈直径を推定することができ、ここIは直交直径を求めるための動脈断面の強度であり
、Iはこの断面近傍であるが動脈外側の強度である。動脈近傍の強度Iを取得する代
わりに、最初にデジタル減算血管造影(DSA)法を適用することによって動脈外側の背
景における任意なポイントを取得することができ、これは、造影剤を有する画像から造影
剤のない画像を減算することによって背景形体を除去することから成る。関係式5~7を
考慮して、ピクセル強度の対数を減算しなければならない。
拡張期血流測定
或る動脈の圧力-血流関係をモデル化するために二次多項式を使用することができ、1
次項は粘性摩擦に関し、2次項は血流分離に関する。従来技術において、多項式の項は、
動脈ジオメトリを解析することによって、また動脈セグメントの各個に適合するジオメト
リベースの解析モデルを適用するために動脈を断片化することによって見出される。した
がって、とくに、狭窄部をモデル化するときには、正確な動脈ジオメトリに極めて影響を
受け易く、複雑な狭窄部は、とくに、断面積の滑らかな変動率及び対称的遷移を有する理
想化された狭窄部としてモデル化される。しかし、このようなジオメトリは現実には観測
されることは決してなく、したがって、これら従来技術は結果に誤差を生じがちにする。
したがって、本発明の他の実施形態は、動脈ジオメトリに関連しないが、その代わりに
機能的特徴に関する動脈特徴付け二次多項式の計算に基づく。
この方法の実施形態を図3で以下に説明し、この場合、画像と圧力測定値との間の同期
性を仮定し、リアルタイムの双方モダリティを同時取得するか、又は双方の取得量にタイ
ムスタンプ付けするかのいずれかを行う。造影剤が動脈の遠位端部に到達するに要する時
Tmnであって、一般的には1回の心拍サイクルよりも短い時間を考慮すると、この方
法は、心拍サイクル内の血流が一定である場合、正確な血流測定をもたらす。しかし、若
干の患者の右冠状動脈(RCA)で収縮期血流が拡張期血流よりも多い2、3の例外を除
いて、拡張期フェーズ中の血流は収縮期フェーズ中の血流よりも相当多い。ステップ10
において、血流が或る程度一定である心臓拡張期中に血流が多く伝播するため、上述の実
施形態で示したように血流は測定されるが、体積流量測定は拡張期又は拡張期の一部分に
限定される。拡張期中の血流がより多いものとして測定するのが好ましいものの、収縮期
又は収縮期の一部分中にも同一測定を実施することができる。安静時の患者では、この拡
張期体積流量Qは、圧力低下Δp =(Pai-Pdiとともに測定し、iは、
所与のレート、代表的には125Hzの圧力測定シーケンスに対応する指数である。負荷
誘導剤(アデノシンのような)又はより一般的には微小血管膨張剤を用いて患者(ストレ
ス付与患者)に充血を誘発させるとともに、患者が充血状態にある間に拡張期中に再び体
積流量Q及び圧力低下Δp を測定する。微小血管膨張剤の使用により、先ず第1血
流条件で測定値を取得し、次いで代表的には適正媒剤で誘発される微小血管膨張を伴う第
2血流条件であって、結果として第1血流条件(通常血流を有し、第1血流条件に比べて
血流が少ない)よりも血流が多くなる第2血流条件で測定値を取得する。
ステップ11において測定した値を使用して、パラメータA及びBを計算することがで
き、A及びBのための関係式8及び9を解くことによって動脈を特徴付けることができる

であり、XはR(安静時)又はS(ストレス付与時)のいずれかを示し、またn及びn
は血流が測定される期間(拡張期又は拡張期の一部分)の境界である。
動脈パラメータA及びBを知ることによって、ステップ12において、任意な圧力低下
(Pa-Pd)に関連する瞬間血流を計算することができる。次のステップでは関係式1
1を用いて安静時又は充血時のいずれかにある患者の全心拍サイクルにわたる各圧力低下
測定値Δp=Pa-Pdのための血流Q を計算する。体積流量は、最終的に関
係式11及び12を使用して関心対象期間n、nにわたり計算することができ、n
及びnは、拡張期若しくは収縮期中のみで、又は他の時間フレーム中で全心拍にわたる
値を得るよう選択される。
微小血管抵抗Rは、次にステップ13において、心拍全体、収縮期若しくは拡張期、
安静時若しくは充血時のいずれかのような関心対象心拍の部分にわたり以下の関係式13
を用いて計算し、ここでP は、Q(関係式14)を計算するのに使用される期間と
同一の関心対象期間にわたる平均遠位圧力である。さらに、充血時と安静時抵抗との間に
おける比のような、微小血管抵抗の任意な組合せを計算することもできる。さらに、とり
わけ以下の関係式15を用いてCFRのような他の指標を計算することもできる。すなわ
ち、
であり、ここでn及びnは全心拍サイクルをカバーするよう選択される。他の境界値
も、安静時血流に対する充血血流の比(すなわち、第1血流条件及び微小血管膨張剤投与
の結果として第1条件におけるよりも多い血流を含む第2血流条件で行った測定値間の比
)を計算するよう選択することができる。
対比誘発充血血流測定
連続充血は、種々の方法によって、最も一般的にはアデノシンの一定静脈内注入を使用
して誘発することができる。静脈内アデノシン又は充血剤の連続注入は、しかし、すべて
のカテーテル検査室、又は「カテラボ(cath lab)」において利用可能ではない。さらに
、アデノシン又は他の薬剤の冠動脈内(IC:intra-coronary)ボーラス投与とともに、
過渡充血を誘発することができる。IC誘発充血の持続時間は極めて短く、オペレータは
充血剤注入と造影剤注入を同期させなければならないため、この方法は課題が多い。他の
実施形態は、上述した実施形態で説明したように、安静条件(又は第1血流条件)におけ
る血流を測定するため1回の造影剤注入した直後に、血流条件を通常から高い状態へと変
化させるよう血流を増加させる微小血管膨張剤として作用する第2造影剤注入を行い、こ
れは造影剤により誘発された充血(又は第2血流条件)と若干同期させる。造影剤は、実
際に満充血のほぼ80%にまで充血を誘発することが知られている。満充血において、動
脈パラメータBの二次項が展開し、また圧力測定が造影剤ベースの血流測定と同期する限
り、第2血流測定をする必要はなく、動脈パラメータA及びBは適正な精度で測定される
ことになる。次に、測定したΔp及びΔpをそれぞれ使用することによって安静条件
及び充血条件に関する血流を計算することができる。
さらに以下に記載するように、血流及び圧力差は、全心拍サイクルにわたり、又は拡張
期若しくは収縮期のみにわたり、測定することができる。微小血管抵抗及び他の血行動態
パラメータは、上述の実施形態で説明したように、血流及び圧力を用い、心拍全体又は拡
張期若しくは収縮期のような期間の一部分を使用して、計算する。
収縮期/拡張期血流測定
図4に示す他の好適な実施形態は、何ら充血剤を必要としない利点を有する。ステップ
20において、体積流量及び光圧力差(Pa-Pd)は、上述の実施形態につき説明した
ように、心拍サイクルの所定期間中に造影剤の波面を追跡するとともに探査下動脈の容積
を決定することによって、測定する。収縮期中及び拡張期中に血流が全く異なることを考
慮すると、体積流量は、収縮期Qsys及び拡張期Qdias中に個別に測定する。これ
ら期間(収縮期対拡張期)は、第1血流条件及び第2血流条件として作用し、安静条件対
充血条件に類似する。1回より多い造影剤注入が収縮期及び拡張期双方をカバーするため
に必要となり得る。さらに、音声信号を使用して、オペレータが適正時点で造影剤を注入
するのに役立てることができる。
安静体積流量Qsys及びQdias、並びに安静平均圧力差Δpsys及びΔpdi
asから、上述の実施形態で説明したようにステップ21で動脈パラメータA及びBを計
算する。
単一の心拍中における圧力差測定は、連続充血を必要とすることなく、またしたがって
、圧力差Δp は、アデノシン等の冠動脈内ボーラス注入に続いて容易に測定される。
ステップ22で示すように、このとき、Δp を使用して、ステップ21でA及びBが
見出される場合に、関係式11の使用により瞬間血流Q を計算する。
体積流量Q (又はQ )及び遠位圧力Pd(又はPd)は関係式12及び1
4を用いて関心対象期間に関して計算する。
微小血管抵抗及び他の血行動態パラメータは上述した実施形態で説明したようにステッ
プ24で計算する。
ジオメトリベースモデルの血流測定向上
動脈内血流を測定するために上述した造影剤波面追跡方法の使用は、この目的のために
適正な精度をもたらすが、結果をさらに向上させる必要があり得る。1-DのCFDモデ
ルは狭窄部を有する動脈には使用できず、これはすなわち、断面積の急激な変化が動脈に
沿う伝播軸線に直交する方向の速度成分であって、1-DのCFDモデルによっては考慮
されない成分を血流に対して発生させるからである。0-D圧力低下モデルも使用できる
が、上述したような圧力低下狭窄部モデルは、臨床的セットアップでは存在し得ない理想
的狭窄ジオメトリを仮定する。
図5に示す本発明の他の実施形態は、上述したように、またこの実施形態のステップ3
0で示すように、造影剤を使用する圧力及び血流の測定で開始し、これに続いて上述した
ように、動脈パラメータの計算を行う。ステップ31において、有意な狭窄部の存在は、
動脈パラメータBが所定閾値を超えるか、又は視覚的に血管造影像で50%超、60%超
、70%超の内腔閉塞割合、又は他の所定閾値内腔閉塞割合を超える局所的狭窄部の存在
を識別するかによって同定する。
ステップ32において、付加的圧力を動脈セグメントに沿って測定し、この測定は、上
述したように、圧力ガイドワイヤをゆっくりと引き戻しながら、セグメントに沿う各圧力
測定に関連する関連画像とともに圧力を同時記録することによって行う。代案として、複
数(2又はそれ以上)の圧力測定を動脈内の既知の位置で記録することができ、最小限に
は、探査下の動脈の近位部分に狭窄部がある場合狭窄部の遠位位置で圧力を測定するか、
又は狭窄部が遠位にある場合狭窄部の近位位置で圧力を測定するか、又は狭窄部に隣接し
て、好適には、狭窄部を含まない動脈の最長セグメントで圧力を測定する。以下に狭窄部
が動脈の近位部分に位置することを仮定する。
ステップ32で他の圧力を測定する動脈セグメントにおける血流は、狭窄部を含んでい
ないので周知の1-DのCFDモデルを用いて計算することができる。1-DのCFDモデ
ルは、以下の関係式16及び17を解くことになるナビエ・ストークス方程式から導き出
すことができる。
ここで、ρは血液濃度であり、αはコリオリ係数であり、またΚは血液粘性及び速度プ
ロファイルに関する抵抗である。動脈が弾性ではなく、分岐が存在しないことを考慮する
と仮定する場合、関係式16の双方の項は、上述したように0に等しい。
動脈セグメントにおける血流を計算するCFDモデル使用の代案としては、ポアズイユ
圧力低下モデル、集中型モデルのような他のジオメトリベースのモデル、又は他の等価な
適切モデルを使用することがある。ポアズイユ法則は関係式18で表現される。血流Qは
、したがって、圧力ガイドワイヤをゆっくりと引き戻しながら圧力測定を行い、またこれ
により動脈の軸線に沿って圧力をマッピングする状況では関係式19の最初の部分を用い
て計算することができ、そうでない場合には関係式19の最後の部分を使用する。
ここで、νは血液粘性、ιは動脈長さ、dは動脈直径である。Qgeomは、全心拍平均
血流、収縮期若しくは拡張期の平均血流のいずれかとすることができ、安静条件又は充血
条件のいずれかに関して計算する。それは、使用されるΔpの部分、及び圧力が安静時
又は充血時のいずれで取得されたかに依存する。
geom が安静時における平均血流であると仮定する。造影剤方法を使用する体積
流量測定は絶対体積流量を決定する際に幾分のバイアス誤差を含み得るが、ステップ30
で得られる2つの血流の比は、バイアス誤差の影響が減少するため、より正確な結果を確
実にもたらす。したがって、先に測定した血流比が正確であると仮定することによって、
ステップ33で充血時血流Qgeom を再計算でき、また補正血流は以下の関係式20
で得られる。
ジオメトリベースの血流Qgeom 及びQgeom 、並びに関連の圧力低下を使用
して、ステップ11又は21の一方により新たな動脈パラメータA及びBを再計算する(
ステップ34)。
動脈のセグメント化
図6に示す有意な分岐7は、分岐の近位容積部が分岐に送給される血流を搬送し、また
したがって、関係式3における動脈容積が比例的に減少するため、血流を過大評価するこ
とになる。血流分割比を決定する1つの方法は、参照により本明細書に組み込まれる、「
循環系輸送時間を決定する方法及びその適用(Methods for the determination of trans
it time circulatory systems and applications of the same)」と題する米国公開特許
公報第2016/0350920号に記載されているような、周知の相対成長法則を使用
することである。その他には、造影剤波面追跡方法を用いて分岐及び動脈における血流を
測定することによって、分割比を決定することができる。関係式3での血流計算に使用さ
れる分岐近位の容積部は、近位容積部に比Qartery/Qbranchを乗算するこ
とによって補正する。
上述の方法は、分岐近位の動脈部分に大きな圧力低下がない場合には、合理的で正確な
結果をもたらす。有意な狭窄がある場合、1つの方法は、分岐6、すなわち、狭窄部の遠
位における圧力を測定するステップから成り、狭窄部の遠位側の動脈セグメントにおける
血流-圧力関係を測定する。他の方法は、2つのセグメントを別個にとり、また2つのセ
ットの動脈パラメータA及びBを計算し、それぞれに対応する血流、次いで微小血管抵抗
、及び他の血行動態パラメータを計算することから成る。
動脈内で圧力ガイドワイヤ2を引き戻しながら、引き戻し移動の経路に沿って圧力測定
を行うことによって、動脈を一連のセグメントに区分し、また動脈特性A及びBを計算し
、次いで関心対象の血行動態パラメータを計算することができる。
側枝付き微小血管抵抗
側枝は動脈を互いに接続する心外膜血管である。図2に示したように、側枝接続は、最
も一般的には圧力ガイドワイヤ2が完全に遠位まで延在するとき、圧力ガイドワイヤのリ
ーチの遠位側である。安静時において、一般的に、側枝は閉じた状態に留まる。充血時(
例えば、充血剤を患者に与えるとき)において、狭窄部の存在は、遠位端部における圧力
を著しく低下させ、供給している動脈の圧力に圧力差を生じ、これにより側枝を開かせる
。側枝は微小循環により血流を増加させるとともに、遠位圧力Pdの増加にも寄与する。
血流及び圧力に対するこの寄与は線形的ではなく、したがって、相殺することはない。
側枝からの不明な血流供給によって引き起こされる誤差を最小化する方法は、楔入圧を
測定し、これを微小血管抵抗の計算に含めることによって側枝の関与を切り離すことから
成る(Martinez et al, Cor. Art. Dis., 2015)。楔入圧は日々の臨床的ワークフローの
一部分ではないため、他の方法は、母集団から導き出したFFRmyoとFFRcor
の間の関係を使用することを開発した。微小血管抵抗の測定における側枝関与をよりよく
考慮する必要性がある。
側枝存在下での収縮期及び拡張期の充血微小血管抵抗は以下の関係式21、すなわち、
および
で表現することができ、ここでQは探査下動脈からの血流であり、またQは側枝血流
である。
上述したように、側枝は充血に応答して開き、またしたがって、安静時微小血管抵抗は
以下の関係式22、すなわち、
および
で表現することができる。
側枝に血流を流し込む圧力を仮定すると、側枝血流は次式、すなわち、
で表現することができ、ここでRは側枝の抵抗であり、Δp d-sys/diasは収
縮期又は拡張期いずれかの充血圧力差Pa-Pbである。側枝抵抗は血流依存ではなく、
またしたがって、RC-sys=RC-diasである。
安静時における拡張期対収縮期の抵抗比Κは以下の関係式、すなわち、
から計算することができる。
充血時における収縮期対拡張期の抵抗比が安静時と同一であると仮定すると、関係式2
1、23及び24は次式、すなわち、
となる。
側枝血流Q は関係式25のRを関係式23に使用して計算することができる。次
に、微小血管抵抗における血流がQ 及びQ の和であることを考慮することによっ
て、血行動態パラメータを計算することができる。
ポストPCIの血行動態評価
安静条件においては、圧力と血流との間の関係は一般的には線形的である。2つ又はそ
れ以上の狭窄部が存在すると、1つの狭窄部にステント挿入による(経皮的冠動脈インタ
ーベンション又はPCIを実施した後の)影響は、排除すべき狭窄部により生ずる圧力低
下を排除することによって容易に決定することができ、Pd/Pa、iFR(瞬間ウェー
ブフリー指数)又はdPR(拡張期圧力比)のような任意な既存安静時指標における遠位
圧力及びポストPCI値を計算することができる。しかし、FFRに関しては、充血時指
数であるため不可能である。1つの狭窄部を排除した後にFFRに対する影響を決定する
方法は楔入圧測定を必要とする。上述したように、このことは臨床的セットアップでは一
般的な慣行ではない。
好適なやり方は、上述した方法のうち1つ、好適には、動脈をセグメント化するステッ
プを含む方法により動脈を特徴付けするステップを伴うものである。セグメント化は、説
明を簡略化するため2つの狭窄部及び動脈の他の部分とともに含むものとするが、セグメ
ントには、開口部から狭窄部の中間までの第1セグメント40及び動脈の中間から遠位端
部までの第2セグメント41があると仮定する(図7参照)。このことは、明らかにびま
ん性疾患がないことを仮定する。上述したように、病変が動脈に沿って拡散している場合
、追加的セグメント化によってモデル化することができる。遠位端部Pdで及び狭窄部P
s相互間での双方において、圧力ガイドワイヤを所定位置に引き戻すことによって、血流
をステップ10又は20により測定し、また圧力をステップ11により測定することを仮
定する。有意な分岐が存在しないと仮定し、このとき各セグメント40及び41に対して
関係式8及び9の2つのセットがある。双方セグメントでの血流は同一であるとともに、
第1セグメントの圧力差はΔps1=p-p及びΔps2=(p-p)-Δp
である。第1セグメント40に関するパラメータAs1及びBs1及び第2セグメント
41に関するパラメータAs2及びBs2は既知である。そこから、平均微小血管抵抗R
を計算する。
びまん性疾患が第1セグメント40に存在しないと仮定すると、一方の狭窄部(狭窄部
1と仮定する)のステント挿入は、このセグメントにおける圧力低下を排除することにつ
ながると仮定することができる。したがって、パラメータの第2セットのみを使用して動
脈をモデル化することができる。しかし、狭窄部1を排除することは、QPost-PC
=PdPost-PCI/Rによって与えられる新たな血流をもたらす。この関係式
を関係式9に挿入すると、次式、すなわち、
にいたり、P を代入することによってP は容易に解ける。このとき、ポストPC
IのFFRはFFR=P /P である。これは一実施例であり、他のセグメント化
を使用して、本明細書記載のびまん性疾患又は他の形体を考慮することもできる。
他のモデル
本明細書の記載は、動脈の圧力-血流関係を特徴付けるため、主に二次モデル、すなわ
ち、関係式8及び9のような線形項及び二次血流項の足し合わせで圧力低下を関係付ける
モデルに基づいているが、他のモデルも使用できる。例えば、収縮期中対拡張期中、又は
安静時中対充血時中のように血流が異なる間に、2つ又はそれ以上の圧力測定値及びこれ
らに対応する(及び同時的)画像を取得するモデルも使用することができる。収縮期対拡
張期は、第1血流条件及び第2血流条件として作用し、また安静時条件対充血時条件も同
様である。非限定的実施例として、他のモデルは、三次項を含むことができる、又は微分
項又は積分項、時間依存項、空間依存項、等々を有することができる。
同時記録(Co-registration)
上述の実施形態が動脈の特定領域に関する血行動態情報を含むことを考慮すると、参照
により全体が本明細書に組み入れられるものとする「」と題する米国公開特許公報第20
06/0052700号に類似するような血管造影図に沿って積分される場合、血行動態
機能の視覚化及び評価を容易にする。
本明細書において微小血管抵抗は純抵抗として考慮するが、微小血管インピーダンスが
抵抗に加えてキャパシタンスを含む(Morris et al., JACC Vol. 8, No. 9 2015)と仮定
することができ、また微小血管抵抗及びキャパシタンスの双方を計算することができる。
本発明の実施形態によれば、計算システムは、適宜上述したような関係式により上述し
た血行動態パラメータを有利に計算することを必要とし、これはすなわち、迅速かつ差し
迫った意思決定に対する介入中に値をリアルタイムで計算することを必要とするからであ
る。計算システムは、命令、とくに、血行動態パラメータを計算するための命令を記憶す
るメモリを備える。このメモリと通信するプロセッサは計算を実施する命令を実行する。
計算システムは、圧力ガイドワイヤ2から圧力測定値を受け取る通信ポート(有線又は
無線とすることができる)を有する。圧力ガイドワイヤ2は、圧力ガイドワイヤから通信
手段(電気的又は光学的な)により信号をレシーバに伝送する圧力センサを有し、レシー
バは、信号を処理して計算システムに供給する、又は計算システムがこのような信号を圧
力ガイドワイヤから受け取るよう構成されている場合には、直接計算システムに供給する
。計算システムは、こうして本明細書記載の実施形態が必要とするような圧力測定値を受
け取る。
撮像機器が撮像する画像も、圧力ガイドワイヤからのデータとともに計算システムに供
給され、これにより圧力測定値は、局所的及び個人付けした動脈ジオメトリのリアル画像
データに同時に補完される。
侵襲性方法に替えて非侵襲性方法を使用する代わりに、侵襲性(圧力ガイドワイヤ)及
び非侵襲性(画像化)の双方を補完物として使用し、より正確な血行動態パラメータを計
算する。計算システムは双方のデータソースを受け取り、血行動態パラメータをリアルタ
イムで計算し、また探査している冠動脈のリアルな局所的条件をよりよく表す血行動態パ
ラメータ値を臨床医に提供し、動脈ジオメトリに関する不正確な一般論を仮定する代わり
に、現実をよりよく反映する血行動態パラメータとしてカテ室でのよりよい意思決定を可
能にする。
好適な実施形態を上述し、また添付図面で図示してきたが、当業者には本開示から逸脱
せず変更を行えることは明らかであろう。このような変更は本発明開示の範囲内に含まれ
る変更例として見なされる。

Claims (18)

  1. ・血管内圧力測定デバイスを用いて動脈内の少なくとも1つの圧力測定値を測定するス
    テップと、
    ・医用撮像機器から前記少なくとも1つの圧力測定値と同期する、動脈の少なくとも1
    つの医用画像を撮像するステップと、
    ・前記少なくとも1つの圧力測定値及び前記少なくとも1つの医用画像の双方を計算シス
    テムに供給するステップと、
    ・前記少なくとも1つの医用画像から血流を計算するステップと、
    ・少なくとも2つの動脈圧力低下及びこれに対応する血流成分から動脈のパラメータを計
    算するステップと、及び
    ・前記血流及び前記動脈のパラメータに基づいて、患者特有血行動態パラメータを計算す
    るステップと
    を備える、方法。
  2. 請求項1記載の方法において、前記少なくとも1つの圧力測定値を測定するステップ及
    び前記少なくとも1つの医用画像を撮像するステップは、血流の2つの条件で実施し、前
    記動脈のパラメータを計算するステップは、前記血流の2つの条件で測定された前記少な
    くとも2つの動脈圧力低下及びこれに対応する血流成分を用いて、動脈のパラメータの二
    次方程式を解くステップを含む、方法。
  3. 請求項2記載の方法において、前記血流の2つの条件は、より多い血流条件及びより少
    ない血流条件を含み、方法は、さらに、より多い血流条件にする薬剤を投与するステップ
    を備え、この後、より多い血流条件の下で前記少なくとも1つの圧力測定値を測定するス
    テップ及び前記少なくとも1つの医用画像を撮像するステップを行う、方法。
  4. 請求項3記載の方法において、前記より多い血流条件の誘発は、部分充血又は満充血を
    誘発する造影剤又は充血剤を注入することによって行う、方法。
  5. 請求項3又は4記載の方法において、さらに、前記少なくとも1つの圧力測定値を測定
    するステップを微小血管膨張剤の注入と同期調整するために、前記より多い血流条件の誘
    発時に音声信号を発生するステップを備える、方法。
  6. 請求項3~5のうちいずれか1項記載の方法において、方法はさらに、前記動脈内に放
    射線吸収造影媒剤を導入するステップを備え、前記少なくとも1つの医用画像を撮像する
    ステップは、さらに、
    前記少なくとも1つの医用画像において、前記動脈の長さに沿って変化する動脈直径を
    測定するステップと、
    前記動脈における前記放射線吸収造影媒剤の伝播を追跡し、また前記動脈の第1ポイン
    トL1から第2ポイントL2に伝播するに要する時間を測定するステップと、
    を有し、
    前記血流計算ステップは動脈内の平均血流を計算するため、次式すなわち、
    の容積Vを前記伝播に要する時間で除算するステップを含む、方法。
  7. 請求項6記載の方法において、平面を画定する前記少なくとも1つの医用画像で、前記
    少なくとも1つの医用画像の前記平面に直交する動脈直径を測定するステップは、前記放
    射線吸収造影媒剤に対する濃度測定法を使用して、前記動脈断面における測定した明暗度
    に対しランベルト・ベールの法則を適用することによってxの関数としてDを測定するス
    テップを含む、方法。
  8. 請求項3~7のうちいずれか1項記載の方法において、さらに、前記動脈内への充血剤
    の冠動脈内又は経皮的注入を使用するステップと、及び前記動脈のパラメータを使用して
    前記充血条件に関連する血流を計算するステップと、を備え、前記患者特有血行動態パラ
    メータを計算するステップは、圧力低下及び前記計算した前記充血条件に関連する血流に
    基づく、方法。
  9. 請求項1~8のうちいずれか1項記載の方法において、前記患者特有血行動態パラメー
    タは微小血管抵抗を含む、方法。
  10. 請求項2記載の方法において、前記少なくとも1つの圧力測定値を測定するステップ及
    び前記少なくとも1つの医用画像を撮像するステップは、安静条件の下で、第1に少なく
    とも収縮期の一部分をカバーする収縮期中に、また第2に少なくとも拡張期の一部分をカ
    バーする拡張期中に実施し、前記動脈のパラメータを計算するステップは、前記収縮期及
    び前記拡張期にそれぞれ測定及び計算された前記少なくとも2つの動脈圧力低下及びこれ
    らに対応する血流成分を使用して前記動脈のパラメータの二次方程式を解くステップを含
    む、方法。
  11. 請求項1~10のうちいずれか1項記載の方法において、さらに、前記動脈内で狭窄部
    の存在を同定するステップと、前記狭窄部の遠位側にある何ら狭窄部がないセグメントを
    同定するステップと、及び前記セグメントに沿って複数の圧力測定値を取得するステップ
    とを備える、方法。
  12. 請求項11項記載の方法において、前記狭窄部の存在を同定するステップは、前記動脈
    のパラメータのうち1つのパラメータを、前記動脈内で狭窄部の存在を同定するための所
    定閾値と比較するステップを含む、方法。
  13. 請求項11項記載の方法において、前記狭窄部の存在を同定するステップは、前記少な
    くとも1つの医用画像を使用して、50%以上の内腔閉塞の存在を決定することによって
    前記動脈内での前記狭窄部の存在を同定するステップを含む、方法。
  14. 請求項11項記載の方法において、前記セグメントに沿って複数の圧力測定値を取得す
    るステップは、血管内圧力測定デバイスの先端を前記セグメントの最も遠位の位置に送り
    、また前記血管内圧力測定デバイスを引き戻すことによって行う、方法。
  15. 請求項11~14のうちいずれか1項記載の方法において、さらに、前記セグメントに
    沿う前記複数の圧力測定値を使用し、またナビエ・ストークス方程式を数値的に解いてジ
    オメトリベースの血流を計算するステップと、並びに前記動脈のパラメータ及び前記少な
    くとも1つの医用画像から血流を計算するステップに対する補正に適用するため、したが
    って、前記狭窄部の存在を考慮する前記患者特有血行動態パラメータに対する補正を適用
    するため、前記ジオメトリベースの血流を使用するステップと、を備える、方法。
  16. 請求項1~15のうちいずれか1項記載の方法において、前記少なくとも1つの圧力測
    定値を測定するステップ及び前記少なくとも1つの医用画像を撮像するステップは、前記
    少なくとも1つの圧力測定値及び前記少なくとも1つの医用画像をリアルタイムで同時に
    取得することによって、又は前記少なくとも1つの圧力測定値及び前記少なくとも1つの
    医用画像にタイムスタンプ付けすることによって、同期させる、方法。
  17. ・血管内圧力測定デバイス及び医用撮像機器と通信する計算システムであり、命令及び前
    記血管内圧力測定デバイス及び医用撮像機器からのデータを記憶するメモリと、及び前記
    命令を実行するプロセッサであって、
    ・動脈内における前記血管内圧力測定デバイスから少なくとも1つの圧力測定値を受け
    取る、
    ・他の圧力デバイスから少なくとも1つの大動脈圧力を受け取る、
    ・前記医用撮像機器から動脈の少なくとも1つの医用画像を受け取る、
    ・前記少なくとも1つの医用画像から血流を計算する、
    ・少なくとも2つの動脈圧力低下及びこれらに対応する血流成分から前記動脈のパラメ
    ータを計算する、及び
    ・前記血流及び前記動脈のパラメータに基づいて、患者特有血行動態パラメータを計算
    する
    ための命令を実行する、該プロセッサと、
    を備える、システム。
  18. 請求項17記載のシステムにおいて、前記血管内圧力測定デバイスは圧力ガイドワイヤ
    である、システム。
JP2023192916A 2018-06-27 2023-11-13 ハイブリッドな画像-侵襲性-圧力血行動態機能評価 Pending JP2024020363A (ja)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US201862690756P 2018-06-27 2018-06-27
US62/690,756 2018-06-27
JP2020568996A JP2021529015A (ja) 2018-06-27 2019-06-27 ハイブリッドな画像−侵襲性−圧力血行動態機能評価
PCT/CA2019/050894 WO2020000102A1 (en) 2018-06-27 2019-06-27 Hybrid image-invasive-pressure hemodynamic function assessment

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2020568996A Division JP2021529015A (ja) 2018-06-27 2019-06-27 ハイブリッドな画像−侵襲性−圧力血行動態機能評価

Publications (1)

Publication Number Publication Date
JP2024020363A true JP2024020363A (ja) 2024-02-14

Family

ID=68985313

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2020568996A Pending JP2021529015A (ja) 2018-06-27 2019-06-27 ハイブリッドな画像−侵襲性−圧力血行動態機能評価
JP2023192916A Pending JP2024020363A (ja) 2018-06-27 2023-11-13 ハイブリッドな画像-侵襲性-圧力血行動態機能評価

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2020568996A Pending JP2021529015A (ja) 2018-06-27 2019-06-27 ハイブリッドな画像−侵襲性−圧力血行動態機能評価

Country Status (5)

Country Link
US (2) US11369277B2 (ja)
EP (1) EP3813654A4 (ja)
JP (2) JP2021529015A (ja)
CA (1) CA3103694A1 (ja)
WO (1) WO2020000102A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230190113A1 (en) * 2020-06-04 2023-06-22 The Regents Of The University Of California Mathematical modeling of blood flow to evaluate hemodynamic significance of peripheral vascular legions
JP2023531566A (ja) * 2020-06-22 2023-07-24 コロヴェンティス リサーチ エービー 微小血管抵抗予備能を決定するためのシステムおよび方法
US11872403B2 (en) 2021-08-06 2024-01-16 Solo Pace Inc. Systems, methods, and apparatus for external cardiac pacing
WO2023014531A1 (en) 2021-08-06 2023-02-09 Solo Pace Inc. Systems, methods, and apparatus for external cardiac pacing
CN116919374B (zh) * 2023-07-19 2024-04-12 西安交通大学 颅内动脉瘤及载瘤动脉内血流动力学参数评估方法及系统

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL8801982A (nl) * 1988-08-09 1990-03-01 Philips Nv Werkwijze en inrichting voor het bepalen van grootheden die kenmerkend zijn voor de stroming van een vloeistof in een vatenstelsel.
US5150292A (en) 1989-10-27 1992-09-22 Arch Development Corporation Method and system for determination of instantaneous and average blood flow rates from digital angiograms
SE0402145D0 (sv) 2004-09-08 2004-09-08 Radi Medical Systems Pressure measurement system
US10217220B2 (en) 2007-01-23 2019-02-26 Dtherapeutics, Llc Methods for the determination of transit time in circulatory systems and applications of the same
CA2926666C (en) * 2009-09-23 2018-11-06 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems, apparatus and methods
US8315812B2 (en) * 2010-08-12 2012-11-20 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
CA2837577C (en) * 2011-05-27 2018-05-29 Lightlab Imaging, Inc. Optical coherence tomography and pressure based systems and methods
US9339348B2 (en) 2011-08-20 2016-05-17 Imperial Colege of Science, Technology and Medicine Devices, systems, and methods for assessing a vessel
JP6181180B2 (ja) * 2012-08-03 2017-08-16 ボルケーノ コーポレイション 脈管を評価するためのシステム
US9595089B2 (en) * 2014-05-09 2017-03-14 Siemens Healthcare Gmbh Method and system for non-invasive computation of hemodynamic indices for coronary artery stenosis
US9349178B1 (en) * 2014-11-24 2016-05-24 Siemens Aktiengesellschaft Synthetic data-driven hemodynamic determination in medical imaging
US10872698B2 (en) * 2015-07-27 2020-12-22 Siemens Healthcare Gmbh Method and system for enhancing medical image-based blood flow computations using physiological measurements
EP3356969B1 (en) 2015-09-29 2020-04-15 Koninklijke Philips N.V. Estimating flow, resistance or pressure from pressure or flow measurements and angiography

Also Published As

Publication number Publication date
US20220248969A1 (en) 2022-08-11
EP3813654A4 (en) 2022-04-06
US20210244293A1 (en) 2021-08-12
JP2021529015A (ja) 2021-10-28
EP3813654A1 (en) 2021-05-05
WO2020000102A1 (en) 2020-01-02
CA3103694A1 (en) 2020-01-02
US11369277B2 (en) 2022-06-28

Similar Documents

Publication Publication Date Title
US20210338186A1 (en) Determining a characteristic of a lumen by measuring velocity of a contrast agent
US10872698B2 (en) Method and system for enhancing medical image-based blood flow computations using physiological measurements
US10803994B2 (en) Vascular flow assessment
US20190110776A1 (en) Methods for Computing Coronary Physiology Indexes Using a High Precision Registration Model
US10984531B2 (en) Determining a luminal-flow-related index using blood velocity determination
JP2024020363A (ja) ハイブリッドな画像-侵襲性-圧力血行動態機能評価
JP7133346B2 (ja) 導管の連続する画像フレームのシーケンスから前記導管を流れる流体を定量的にフロー分析する装置の作動方法および撮像デバイス
CN112543618B (zh) 用于进行定量血液动力学流量分析的方法和装置
Seike et al. Intracoronary optical coherence tomography-derived virtual fractional flow reserve for the assessment of coronary artery disease
CN113711317A (zh) 用于评估冠状动脉疾病的手段和设备
US20210251504A1 (en) System and method for determining pulse wave velocity of vessels
US20230190113A1 (en) Mathematical modeling of blood flow to evaluate hemodynamic significance of peripheral vascular legions
Ágoston et al. The pressure-derived microvascular resistance reserve and its correlation to Doppler MRR measurement—a proof of concept study
EP3643244A1 (en) System and method for determining pulse wave velocity of vessels
Jain et al. Physiological Assessment During Interventional Procedures

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231211

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231211