JP2023553556A - 血管内イメージングに基づく血流速度、血流予備量比の分析方法 - Google Patents
血管内イメージングに基づく血流速度、血流予備量比の分析方法 Download PDFInfo
- Publication number
- JP2023553556A JP2023553556A JP2023535452A JP2023535452A JP2023553556A JP 2023553556 A JP2023553556 A JP 2023553556A JP 2023535452 A JP2023535452 A JP 2023535452A JP 2023535452 A JP2023535452 A JP 2023535452A JP 2023553556 A JP2023553556 A JP 2023553556A
- Authority
- JP
- Japan
- Prior art keywords
- blood flow
- blood
- flow velocity
- examined
- blood vessel
- 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
Links
- 230000017531 blood circulation Effects 0.000 title claims abstract description 133
- 238000003384 imaging method Methods 0.000 title claims abstract description 47
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 78
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000012014 optical coherence tomography Methods 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000002608 intravascular ultrasound Methods 0.000 claims abstract description 13
- 239000008280 blood Substances 0.000 claims description 22
- 210000004369 blood Anatomy 0.000 claims description 22
- 230000000747 cardiac effect Effects 0.000 claims description 21
- 239000012530 fluid Substances 0.000 claims description 11
- 230000003287 optical effect Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 abstract description 8
- 230000004089 microcirculation Effects 0.000 abstract description 8
- 238000005457 optimization Methods 0.000 abstract 1
- 239000003795 chemical substances by application Substances 0.000 description 13
- 239000002872 contrast media Substances 0.000 description 8
- 210000004351 coronary vessel Anatomy 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 239000007788 liquid Substances 0.000 description 4
- 238000013146 percutaneous coronary intervention Methods 0.000 description 4
- 239000002504 physiological saline solution Substances 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 201000001320 Atherosclerosis Diseases 0.000 description 1
- 208000037260 Atherosclerotic Plaque Diseases 0.000 description 1
- 201000000057 Coronary Stenosis Diseases 0.000 description 1
- 206010011089 Coronary artery stenosis Diseases 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002498 deadly effect Effects 0.000 description 1
- 230000035487 diastolic blood pressure Effects 0.000 description 1
- 230000000916 dilatatory effect Effects 0.000 description 1
- 238000007865 diluting Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000002565 electrocardiography Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 230000002966 stenotic effect Effects 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 231100000216 vascular lesion Toxicity 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
Abstract
本発明は、血管内イメージングに基づく血流速度の分析方法、血流予備量比の分析方法、及び微小循環抵抗分析方法を提供する。血管内イメージングに基づく血流速度の分析方法は、光干渉断層画像又は血管内超音波画像である検査対象の血管の第一腔内画像を取得するステップS1と、前記第一腔内画像に基づいて、血流信号を取得して、前記血流信号に基づいて前記検査対象の血管内の血流速度を計算するステップS2と、を含む。本発明の実施例の血管内イメージングに基づく血流速度の分析方法によれば、血管内のダイナミックな血流値、即ち、血流速度を正確に測定することができ、より多くの血管内の血流情報を得ることができ、かつFFR計算、微小循環抵抗の計算等の最適化に用いることができる。【選択図】図1
Description
本発明は、検査対象となる血管内インターベンション画像の医療器械分野に関し、特に、血管内イメージングに基づく血流速度、血流予備量比、微小循環抵抗の分析方法に関する。
冠状動脈疾病は既に世界で一番患者数の多い致死的な疾病になっている。現在、経皮的冠状動脈インターベンション(percutaneous coronary intervention,PCIと略称する)は冠状動脈疾病の有効な治療方法の1つである。
光干渉断層イメージング(OCT)、血管内超音波(IVUS)等の腔内画像技術は、動脈粥状硬化症病変プラークの識別、ステントの設置及び評価に対して臨床的に極めて重要な意義があり、既にPCI治療の指導に用いられた。その中で、OCTが高解像度(縦方向:10μm、横方向:20~40μm)を有しているので、冠動脈の構造や形態的特徴の識別及び破損し易いプラークの検出において圧倒的な優勢を有しているが、形態と構造において動脈粥状硬化症及び冠動脈狭窄を評価するだけに限定されて、検査対象となる血管形態変化による冠動脈血流と供血等の機能的な変化を正確に反映することができない。
一方、血流予備量比(fractional flow reserve,FFRと略称する)は、冠状動脈の狭窄部の遠位端と近位端との圧力差の値を測定することにより、狭窄病変の検査対象の血管の供血機能への影響を効果的に反映することができ、冠動脈から血液が供給される心筋の虚血を引き起こすか否かを評価することができる。
現在、FFRは、既に臨床においてPCI治療を診断し、指導して評価するゴールドスタンダードになっている。しかしながら、従来のFFRは、圧力ガイドワイヤーを介して血液の圧力を測定する必要があるが、このような検査操作が複雑で、費やす時間が比較的長く、所要の手術消耗品(FFRガイドワイヤー)の値段が高い。また、検査対象の血管への拡張薬物の注射による副作用は患者に気分が悪い反応をさせることになり、ガイドワイヤーのインターベンションプロセスにおいて患者の検査対象の血管への損傷を引き起こし易い。上記の原因によりFFRの普及が制限されている。
また、従来の画像処理学的な計算によってFFRを得る方法は、通常、流体力学血流モデルを計算する境界条件として一定の平均血流圧力を使用することに基づくものである。しかし、実際の血流圧力は何時までも変わらないものではなく、被検査者の個体化の差異及び下流血管の病変程度のいずれも実際の血流状況に影響を与えている。
したがって、FFR計算をさらに最適化すると共に、さらに心周期中の異なる状態間の差圧を把握することや、微小循環抵抗等を正確に計算することなどに用いられるように、血管内のダイナミックな血流値を正確に測定できる方法は差し迫って必要とされている。
これに鑑みて、本発明の目的は、血管内イメージングに基づく血流速度を簡単で、高速で正確に計算できる分析方法を提供することにある。
本発明は、高速で、正確な血流予備量比の分析方法をさらに提供する。
このほか、本発明は、高速で、正確な微小循環抵抗の分析方法をさらに提供する。
上記の目的を達成するために、本発明は下記の技術手段を採用する。
本発明の第一側面の実施例の血管内イメージングに基づく血流速度の分析方法において、光干渉断層画像又は血管内超音波画像である検査対象の血管の第一腔内画像を取得するステップS1と、
前記第一腔内画像に基づいて血流信号を取得し、前記血流信号に基づいて前記検査対象の血管内の血流速度を計算するステップS2と、を含む。
前記第一腔内画像に基づいて血流信号を取得し、前記血流信号に基づいて前記検査対象の血管内の血流速度を計算するステップS2と、を含む。
本発明の実施例の血流速度の分析方法によれば、血管内のダイナミックな血流値、即ち血流速度を正確に測定することができ、より多くの血管内の血流情報を得ることができる。
第一腔内画像としては、光干渉断層画像(OCT画像)又は血管内超音波画像(IVUS画像)であってもよい。その中で、OCTは、高解像度(縦方向:10μm、横方向:20~40μm)を有しているので好ましい。
さらに、前記ステップS2では、前記血流速度はドップラーシフトにより得られる。つまり、画像における同一位置での前記血流信号の位相差を2回測定することで血流速度を計算する。ドップラーシフトによりダイナミックな血流速度を正確に計算することができる。
さらに、前記ステップS1は、
イメージングカテーテルを前記検査対象の血管の所定位置に挿入するステップS11と、
前記イメージングカテーテルを前記所定位置で回転させると共に、1つ又は複数の心周期である所定時間内に所定の時間間隔で走査し、前記第一腔内画像を得るステップS12と、を含み、
ステップS2では、隣り合う2回の走査によって得られた前記血流信号に基づいてドップラーシフトを取得し、前記ドップラーシフトに基づいて前記血流速度を確定する。
イメージングカテーテルを前記検査対象の血管の所定位置に挿入するステップS11と、
前記イメージングカテーテルを前記所定位置で回転させると共に、1つ又は複数の心周期である所定時間内に所定の時間間隔で走査し、前記第一腔内画像を得るステップS12と、を含み、
ステップS2では、隣り合う2回の走査によって得られた前記血流信号に基づいてドップラーシフトを取得し、前記ドップラーシフトに基づいて前記血流速度を確定する。
つまり、先ず、イメージングカテーテルを検査対象の血管の所定位置に挿入し、この後、所定時間内に所定の時間間隔で走査し、隣り合う2回の走査によって得られた前記血流信号を確定することで、ドップラーシフトを取得することができる。これによって血流速度を確定する。その中で、1つの心周期内の異なる時期に、例えば、拡張期血圧又は収縮期血圧状態で、その血流速度が異なるので、1つ又は複数の心周期をモニタリングすることで、より多くの血流情報を得ることができる。
その中で、前記心周期は、例えばX線造影又は心電図により測定される。つまり、X線造影技術又は心電図技術を組み合わせて心周期をモニタリングし、モニタリングして得られた心周期の情報に基づき、測定された血流速度に対応する心周期の状態を確定し、より多くの血管及び血流情報を取得することができる。
さらに、前記ドップラーシフトは同一位置において隣り合う時刻で2回走査することで得られた画像中の前記血流信号により取得され、又は、前記ドップラーシフトは異なるフレーム画像中の同一位置での前記血流信号により取得される。
つまり、ドップラーシフトの測定は、同一位置において近い時刻で2回走査した血流信号によって得られることができ、例えば、イメージングカテーテルが検出される検査対象の血管の遠位端にある場合、遅い速度で回転する状態になり、隣り合う走査Aの間に十分な重なり合いを有することを保証することによって、同一位置における、近い時刻の信号の位相の変化を近似して得ることができる。又は、ドップラーシフトは異なるフレーム画像中の同一位置における血流信号により得られ、異なるフレーム画像の間にはイメージングエリアの重なり合いも存在していることになるので、異なるフレーム画像に対応する同一位置の信号を得ることもできる。
さらに、前記第一腔内画像は光干渉断層画像であり、前記ステップS12では、走査する前に、前記検査対象の血管に第一所定量のリンス剤を注入し、前記リンス剤と血液を混合させる。
遠位端の血流速度を測定する場合、先ずイメージングカテーテルを検査対象の血管の遠位端に挿入し、血液によるOCT光の強い散乱と吸収のため、検査対象の血管中に少量のリンス剤を注入して血液と混合させて一部の血液を押し流すことによって、光信号の侵入深さを高める必要があり、リンス剤としては、例えば、X線造影剤又は生理食塩水であってもよい。リンス剤と検査対象の血管中の血液を混合した後、イメージングカテーテルを回転させると共に、同一位置における複数の血流信号を取得し、複数の血流信号に基づいてドップラーシフトを取得して、ドップラーシフトに基づいて最大血流速度を確定する。
従来の血流速度測定方法には温度希釈法が知られており、既知温度の冷食塩水を一定の速度でカテーテルを介して冠状動脈に注射し、冠動脈内に挿入された温度センサーにより血流温度降下の幅と時間を測定することで、血流速度を推定する。もし該方法を使用して血流速度の検出を行うと、血流速度の検出、腔内画像の取得はステップを分けて行う必要がある。
本実施例の分析方法によれば、OCT画像を用いてドップラーシフトを取得し、走査する前に、検査対象の血管中に少量のリンス剤、例えば、X線造影剤、生理食塩水を注入して検査対象の血管中の血液を希釈することによって、その中の血流信号のドップラーシフトを取得することができ、更に該ドップラーシフトに基づいてダイナミックな血流値を取得する。
具体的には、血液中に適量の標的物(即ちX線造影剤、生理食塩水)を加えることで、心臓の鼓動で該標的物を冠動脈に入らせ、測定位置、即ちOCTレンズ位置に到達させ、先に到達した標的物と血液が混合されたため、この部分の混合液体から戻ったOCT信号はシステムにより検出されることができ、異なる時刻で取得した信号の間に血液の流動によってドップラーシフトが生じ、ドップラーシフトを測定することで、下式によって混合液体の流速(即ち血流速度)を算出することができる。
ここで、Vzは検査対象の血管の方向に沿って流れる血流速度、λcはOCT光源の中心波長、Δφは2回の検査間の位相変化、ΔTは2回の検査間の時間間隔、nは被測定液体の屈折率、θはカテーテルの射出光束と血流方向との夾角である。
計算式から分かるように、中心波長、屈折率、角度が既知であるため、2回の検査の時間間隔、及び該2回の検査に対応するドップラーシフト(即ち位相変化)を確定しさえすれば、対応の血流速度を得ることができる。
システムの位相分解能が一定である場合、測定可能な血流速度の最大値及び精度は時間間隔ΔTによって決まる。つまり、時間間隔を調整することで、より高い精度を得ることができる。また、ΔTは実際に隣接するライン間の走査通過の時間間隔に対応し、もし、ΔTを小さく設定しすぎると、解像度の制限を受けて正確に識別することができず、この時、時間間隔を適切に増加させることでサンプリング密度を減らすことができる。
本発明の第二側面の実施例の血流予備量比の分析方法において、
上記のいずれか1つに記載の血流速度の分析方法によって分析して得られる血流速度であって、検査対象の血管中の1つ又は複数の心周期内の血流速度の平均流速を取得するステップS10と、
前記検査対象の血管の第二腔内画像を取得し、前記第二腔内画像及び前記平均流速に基づいて流体力学分析計算を行い、前記血流予備量比を得るステップS20と、を含み、
前記第二腔内画像は光干渉断層画像又は血管内超音波画像である。
上記のいずれか1つに記載の血流速度の分析方法によって分析して得られる血流速度であって、検査対象の血管中の1つ又は複数の心周期内の血流速度の平均流速を取得するステップS10と、
前記検査対象の血管の第二腔内画像を取得し、前記第二腔内画像及び前記平均流速に基づいて流体力学分析計算を行い、前記血流予備量比を得るステップS20と、を含み、
前記第二腔内画像は光干渉断層画像又は血管内超音波画像である。
腔内画像に基づいて数値流体動力学シミュレーションを直接行う場合、通常、血流速度を1つの定数に設定する。しかし、実際の血流速度は被検査者の個体差異、病変状況、及び心周期等によって変化が発生し、リアルタイムの血流測定結果は検査対象の血管内の血流分析をより正確に測定・分析し、流体力学シミュレーションと血流予備量比の計算をさらに最適化することができる。つまり、本発明の血流予備量比の分析方法によれば、第一腔内画像により心周期全体の血流速度の平均流速を確定すると共に、該検査対象の血管の第二腔内画像を合わせて流体動力学分析計算を行って、血流予備量比(FFR)を得る。従来のFFRが流体力学血流モデルを計算する境界条件として固定の平均血流圧力を用いるのに比べて、検査対象の血管の腔内画像に基づいてそのリアルタイムの血流状況を計算し、腔内画像から得られる検査対象の血管の具体的な構造と組み合わることで、より正確な血流予備量比を算出することができる。
さらに、前記ステップS10では、前記血流速度は前記検査対象の血管の遠位端及び/又は近位端の平均速度を含み、前記ステップS20では、前記第二腔内画像を得た後、前記第二腔内画像に基づいて血管のセグメンテーションを行い、血管の構造形態を三次元的に再構成して、前記平均流速に基づいて流体動力学分析計算を行い、前記血流予備量比を得る。
つまり、遠位端の平均速度を用いてもよく、近位端の平均速度を用いてもよく、さらに、遠位端及び近位端の平均速度を同時に用いて第二腔内画像から得られた血管の構造形態と合わせて、流体動力学分析計算によってFFR値を得ることもできる。その中で、遠位端位置における血管が比較的細くて、イメージングカテーテルはその平均流速に大きな影響を及ぼす可能性があるものの、手順から言うと、遠位端で速度を測定し、後退して第二腔内画像を得ることは、ワンステップで行うことができ、手順がさらに簡素化される。したがって、具体的な条件に合わせて遠位端の平均流速、近位端の平均流速、遠位端と近位端の総合平均流速を選択して使用することができる。
さらに、前記ステップS20は、
イメージングカテーテルを前記検査対象の血管の遠位端に到達させることと、
前記イメージングカテーテルを回転させると共に、前記検査対象の血管の近位端に後退させて、前記検査対象の血管の前記第二腔内画像を得ることと、を含む。
イメージングカテーテルを前記検査対象の血管の遠位端に到達させることと、
前記イメージングカテーテルを回転させると共に、前記検査対象の血管の近位端に後退させて、前記検査対象の血管の前記第二腔内画像を得ることと、を含む。
つまり、第二腔内画像を取得する具体的なフローとして、先ずイメージングカテーテルを前記検査対象の血管の遠位端に置き、この後、イメージングカテーテルを回転させると共に、前記検査対象の血管の近位端に後退させ、後退途中で走査を行って検査対象の血管全体の構造情報を反映する第二腔内画像を得る。
血流流速の測定を組み合わせて言うと、次のように行うことができる。先ずイメージングカテーテルを血管遠位端に到達させ、この後、イメージングカテーテルを回転させて該位置における第一腔内画像を取得し、この後、イメージングカテーテルを回転させると共に、後退して走査を行い、第二腔内画像を取得し、第一腔内画像に基づいて遠位端の血流速度の平均流速を解析し、第二腔内画像に基づいて血管の構造形態を取得し、平均流速及び構造形態に基づいて、流体動力学解析、例えば、CFDシミュレーションによって、FFR値を得ることができる。
これにより、1回の後退で血流速度及び腔内画像を同時に取得することができ、測定のステップと時間を減らして、血管の病変に対する機能性分析の正確性と効率を大いに高めた。
さらに、前記第二腔内画像は光干渉断層画像であり、前記イメージングカテーテルを後退させる前に、前記検査対象の血管内に第二所定量のリンス剤を注入して前記血液を押し流す。
つまり、OCT画像を第二腔内画像として使用する場合、後退する前に、血管内に大量のリンス剤(例えば、X線造影剤、生理食塩水等)を注入して血液を押し流す必要があり、血管内の構造形態もより正確に反映する。
本発明の第三側面の実施例の微小循環抵抗の分析方法において、
検査対象の血管の近位端の圧力を取得するステップと、
上記いずれか1つに記載の血流予備量比の分析方法に基づいて血流予備量比を分析するステップと、
近位端の圧力と前記血流予備量比とに基づいて、遠位端の圧力を得るステップと、
前記遠位端の圧力に基づき、及び上記いずれか1つに記載の血流速度の分析方法によって得られた検査対象の血管の遠位端の1つ又は複数の心周期に相応する血流速度の平均流速に基づいて、前記検査対象の血管の前記微小循環抵抗を確定するステップと、を含む。
検査対象の血管の近位端の圧力を取得するステップと、
上記いずれか1つに記載の血流予備量比の分析方法に基づいて血流予備量比を分析するステップと、
近位端の圧力と前記血流予備量比とに基づいて、遠位端の圧力を得るステップと、
前記遠位端の圧力に基づき、及び上記いずれか1つに記載の血流速度の分析方法によって得られた検査対象の血管の遠位端の1つ又は複数の心周期に相応する血流速度の平均流速に基づいて、前記検査対象の血管の前記微小循環抵抗を確定するステップと、を含む。
つまり、血流速度の平均流速により解析してFFR値を得た後、血管近位端の圧力から遠位端の圧力を確定し、遠位端の圧力及び平均流速によって微小循環抵抗を確定することができる。
本発明の実施例の目的、技術手段及び優れた点をよりはっきりさせるために、以下、本発明の実施例の図面を参照しながら、本発明の実施例の技術手段を明確かつ十分に記載する。記載される実施例は本発明の一部の実施例であり、全部の実施例ではないことは明らかである。記載された本発明の実施例に基づいて、当業者が得られたあらゆるその他の実施例のいずれも、本発明の保護範囲に属する。
本発明に使用された技術用語又は科学用語は、別に定義されない限り、本発明の属する分野における通常の知識を有する者によって理解される通常の意味であるべきである。本発明に使用される「第一」、「第二」及び類似の語句は、如何なる順序、数量又は重要性を表すものではなく、異なる構成部分を区分するためのものだけである。同様に、「1つ」又は「一」等の類似語句も数量を制限することを表すものではなく、少なくとも1つが存在することを表す。「接続」又は「繋がり」等の類似語句は物理的又は機械的接続を限定するものではなく、直接的または間接的にかかわらず、電気的な接続を含んでもよい。「上」、「下」、「左」、「右」等は相対位置関係を表すためのものだけであり、説明される対象の絶対位置が変わった後、該相対位置関係も相応に変わる。
以下、本発明の血流速度の分析方法、血流予備量比の分析方法、及び微小循環抵抗の分析方法について、具体的な実施例によりさらに詳細に説明する。
実施例1:OCT画像のドップラーシフトにより血流速度を分析する。
図1に示すように、本実施例の血流速度の分析は、
1a)イメージングカテーテルを検査対象の血管の所定位置に挿入するステップと、
1b)血管内に少量のX線造影剤を注入するステップと、
1c)X線造影剤と血液との混合液が該所定位置に到達した時に、イメージングカテーテルをその場で回転させて、所定の時間間隔、例えば1s(秒)で走査を行い、
同時に、心電図により心周期をモニタニングし、前記イメージングカテーテルが1つ又は複数の心周期の第一腔内画像を取得することを確定するステップと、
1d)隣り合う2回の走査によって得られた第一腔内画像につき、下式によって血流速度を算出するステップと、を含み、
1a)イメージングカテーテルを検査対象の血管の所定位置に挿入するステップと、
1b)血管内に少量のX線造影剤を注入するステップと、
1c)X線造影剤と血液との混合液が該所定位置に到達した時に、イメージングカテーテルをその場で回転させて、所定の時間間隔、例えば1s(秒)で走査を行い、
同時に、心電図により心周期をモニタニングし、前記イメージングカテーテルが1つ又は複数の心周期の第一腔内画像を取得することを確定するステップと、
1d)隣り合う2回の走査によって得られた第一腔内画像につき、下式によって血流速度を算出するステップと、を含み、
ここで、Vzは検査対象の血管方向に沿って流れる血流速度、λcはOCT光源の中心波長、Δφは2回の検査間の位相変化、ΔTは2回の検査間の時間間隔、nは被測定液体の屈折率、θはカテーテルの射出光束と血流方向との夾角である。
1つの心周期全体における異なる状態にそれぞれ対応する血流速度が得られる。
実施例2:OCT画像により血流予備量比を分析する。
図2に示すように、高速OCTシステムによって血流予備量比を分析する。本実施例では高速OCTシステムを用いて1回で検出することで、血流速度の平均流速及び構造形態を解析するための腔内画像(即ちOCT後退画像)を取得した。
その中で、図2に示すように、高速OCTシステムは、OCTイメージングシステム101、走査装置102、電動制御機構103、及びカテーテル104を含んで構成される。その中で、OCTイメージングシステム101は、OCT光源、及び画像採集システムを含む。
OCT後退画像を取得する操作において、電動制御機構103は、カテーテル104(即ちイメージングレンズ)を検査対象の血管105の遠位端に移動させ、この後、電動制御機構103は、360度の回転走査を行うように走査装置102を制御して、後退により螺旋状3D走査を実現し、OCT画像を得る。血液によるOCT光の強い散乱と吸収を考慮し、カテーテル104が指定の位置、即ち遠位端に到達した後、血管内に大量のリンス剤、例えば、X線造影剤又は生理食塩水を注入して血液を押し流し、リンス剤が遠位端に到達した場合(つまり、OCTイメージングシステム101はOCT信号をはっきり取得できる場合)、電動制御機構103により、360度回転走査を行って後退するように走査装置102を制御する。
また、血流信号のドップラーシフトを取得するために、具体的な操作は次の通りである。電動制御機構103はカテーテル104(つまり、イメージングレンズ)を検査対象の血管105の所定位置、例えば遠位端に移動させ、この後、血管内に少量のリンス剤、例えばX線造影剤又は生理食塩水を注入して該リンス剤と血液を混合させ、混合した血液が遠位端に到達した時、電動制御機構103は所定時間内(例えば、1つ又は複数の心周期の時間に相当)に所定の時間間隔で360度の回転走査を行うように走査装置102を制御することで、異なる時刻の血流信号を取得し、異なる時刻の血流信号からドップラーシフトを取得し、ドップラーシフトに基づいて血流速度を確定する。
具体的には、図3は本実施例の血流予備量比の分析方法のフローを示す図である。
図3に示すように、本実施例に係る血流予備量比の分析方法は、
2A)検査対象の血管の遠位端の血流速度の心周期全体内の平均流速を取得し、
具体的には、実施例1と同様の方法を用いて分析・計算することで取得するステップと、
2B)血管内に大量のX線造影剤を注入して血液を押し流すステップと、
2C)X線造影剤が遠位端に到達した場合、イメージングカテーテルを回転させると共に、近位端まで後退させて、OCT後退画像を得るステップと、
2D)OCT後退画像に対して血管のセグメンテーションを行い、血管の構造形態を三次元的に再構成するステップと、
2E)上記得られた遠位端の平均流速、及び血管の構造形態に基づいて、流体動力学シミュレーションを行い、FFR値を得るステップと、を含む。
2A)検査対象の血管の遠位端の血流速度の心周期全体内の平均流速を取得し、
具体的には、実施例1と同様の方法を用いて分析・計算することで取得するステップと、
2B)血管内に大量のX線造影剤を注入して血液を押し流すステップと、
2C)X線造影剤が遠位端に到達した場合、イメージングカテーテルを回転させると共に、近位端まで後退させて、OCT後退画像を得るステップと、
2D)OCT後退画像に対して血管のセグメンテーションを行い、血管の構造形態を三次元的に再構成するステップと、
2E)上記得られた遠位端の平均流速、及び血管の構造形態に基づいて、流体動力学シミュレーションを行い、FFR値を得るステップと、を含む。
実施例3:血管内超音波画像(IVUS画像)により血流予備量比を分析する。
図4は本実施例の血流予備量比の分析方法のフローを示す図である。
図4に示すように、本実施例の血流予備量比の分析方法は、以下のステップを含む。
3A)検査対象の血管の遠位端の血流速度の心周期全体内の平均流速を取得する。
該ステップでは、実施例1と異なるのは、血管内にリンス剤を注入する必要がない。つまり、ステップ1b)を実行する必要がなく、イメージングカテーテルが遠位端に到達した後、直ちに該イメージングカテーテルを回転させて血管内超音波画像を取得し、この後、該血管内超音波画像に基づいてドップラーシフトを取得して、ドップラーシフトに基づいて血流速度を確定してもよい。具体的なその他の細部は実施例1を参照して行ってもよい。ここでは繰り返して述べない。
3B)イメージングカテーテルを回転させると共に、近位端まで後退させ、IVUS後退画像を得る。
3C)IVUS後退画像に対して血管のセグメンテーションを行い、血管の構造形態を三次元的に再構成する。
3D)上記得られた遠位端の平均流速及び血管の構造形態に基づいて、流体動力学シミュレーションを行い、FFR値を得る。
3C)IVUS後退画像に対して血管のセグメンテーションを行い、血管の構造形態を三次元的に再構成する。
3D)上記得られた遠位端の平均流速及び血管の構造形態に基づいて、流体動力学シミュレーションを行い、FFR値を得る。
実施例4:微小循環抵抗分析
図5は本実施例の微小循環抵抗分析のフローを示す図である。
図5に示すように、本実施例の微小循環抵抗の分析方法は、以下のステップを含む。
4A)平均流速V及びFFR値を確定する。
具体的には、実施例1の方法を参照して平均速度Vを求め、実施例2の方法を参照してFFR値を求めてもよい。
4A)平均流速V及びFFR値を確定する。
具体的には、実施例1の方法を参照して平均速度Vを求め、実施例2の方法を参照してFFR値を求めてもよい。
4B)近位端圧力Paを取得して、遠位端の圧力Pdを計算し、ここで、Pd=Pa*FFRである。
ここで、近位端圧力Paは圧力ガイドワイヤーにより測定で得られたものであってもよいし、経験に基づいて設定されたものなどであってもよい、
4C)平均流速及び遠位端の圧力に基づいて、微小循環抵抗IMRを計算し、ここで、
IMR=Pd/Vである。
ここで、近位端圧力Paは圧力ガイドワイヤーにより測定で得られたものであってもよいし、経験に基づいて設定されたものなどであってもよい、
4C)平均流速及び遠位端の圧力に基づいて、微小循環抵抗IMRを計算し、ここで、
IMR=Pd/Vである。
以上は、本発明の好適な実施形態であり、当業者にとって、本発明の上記原理から逸脱しない前提下で、さらに若干の改良や潤色を行ってもよく、これらの改良や潤色も本発明の保護範囲に入ると見なされるべきであることに注意されたい。
Claims (10)
- 血管内イメージングに基づく血流速度の分析方法において、
光干渉断層画像又は血管内超音波画像である検査対象の血管の第一腔内画像を取得するステップS1と、
前記第一腔内画像に基づいて血流信号を取得し、前記血流信号に基づいて前記検査対象の血管内の血流速度を計算するステップS2と、
を含むことを特徴とする血管内イメージングに基づく血流速度の分析方法。 - 前記ステップS2では、前記血流速度は、ドップラーシフトにより得られることを特徴とする請求項1に記載の分析方法。
- 前記ステップS1は、
イメージングカテーテルを前記検査対象の血管の所定位置に挿入するステップS11と、
前記イメージングカテーテルを前記所定位置で回転させると共に、1つ又は複数の心周期である所定時間内に所定の時間間隔で走査し、前記第一腔内画像を得るステップS12と、を含み、
ステップS2では、隣り合う2回の走査によって得られた前記血流信号に基づいてドップラーシフトを取得し、前記ドップラーシフトに基づいて前記血流速度を確定することを特徴とする請求項2に記載の分析方法。 - 前記心周期は、X線造影或いは心電図により測定され、
又は、前記ドップラーシフトは、同一位置において隣り合う時刻で2回走査した画像中の前記血流信号により得られ、
又は、前記ドップラーシフトは、異なるフレーム画像における同一位置での前記血流信号により得られることを特徴とする請求項3に記載の分析方法。 - 前記第一腔内画像は、光干渉断層画像であり、前記ステップS12では、走査する前に、前記検査対象の血管に第一所定量のリンス剤を注入し、前記リンス剤と血液を混合させることを特徴とする請求項3に記載の分析方法。
- 血流予備量比の分析方法において、
請求項1~5のいずれか1項に記載の血管内イメージングに基づく血流速度の分析方法によって分析して得られる血流速度であって、検査対象の血管中の1つ又は複数の心周期内の血流速度の平均流速を取得するステップS10と、
前記検査対象の血管の第二腔内画像を取得し、前記第二腔内画像及び前記平均流速に基づいて流体力学分析計算を行い、前記血流予備量比を得るステップS20と、を含み、
前記第二腔内画像は光干渉断層画像又は血管内超音波画像であることを特徴とする血流予備量比の分析方法。 - 前記ステップS10では、前記血流速度は、前記検査対象の血管の遠位端及び/又は近位端の平均速度を含み、
前記ステップS20では、前記第二腔内画像を得た後、前記第二腔内画像に基づいて血管のセグメンテーションを行い、血管の構造形態を三次元的に再構成して、前記平均流速に基づいて流体動力学分析計算を行い、前記血流予備量比を得ることを特徴とする請求項6に記載の血流予備量比の分析方法。 - 前記ステップS20は、
イメージングカテーテルを前記検査対象の血管の遠位端に到達させることと、
前記イメージングカテーテルを回転させると共に、前記検査対象の血管の近位端に後退させ、前記検査対象の血管の前記第二腔内画像を得ることと、を含むことを特徴とする請求項7に記載の血流予備量比の分析方法。 - 前記第二腔内画像は、光干渉断層画像であり、前記イメージングカテーテルを後退させる前に、前記検査対象の血管内に第二所定量のリンス剤を注入して前記血液を押し流すことを特徴とする請求項8に記載の血流予備量比の分析方法。
- 微小循環抵抗の分析方法において、
検査対象の血管の近位端の圧力を取得するステップと、
請求項6~9のいずれか1項に記載の血流予備量比の分析方法に基づいて血流予備量比を分析するステップと、
近位端の圧力と前記血流予備量比とに基づいて、遠位端の圧力を得るステップと、
前記遠位端の圧力、及び請求項1~4のいずれか1項に記載の血流速度の分析方法によって得られた検査対象の血管の遠位端の1つ又は複数の心周期に相応する血流速度の平均流速に基づいて、前記検査対象の血管の前記微小循環抵抗を確定するステップと、
を含むことを特徴とする微小循環抵抗の分析方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110474542.2A CN113180631A (zh) | 2021-04-29 | 2021-04-29 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
CN202110474542.2 | 2021-04-29 | ||
PCT/CN2022/089546 WO2022228464A1 (zh) | 2021-04-29 | 2022-04-27 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023553556A true JP2023553556A (ja) | 2023-12-22 |
Family
ID=76980914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023535452A Pending JP2023553556A (ja) | 2021-04-29 | 2022-04-27 | 血管内イメージングに基づく血流速度、血流予備量比の分析方法 |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP2023553556A (ja) |
CN (1) | CN113180631A (ja) |
WO (1) | WO2022228464A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113180631A (zh) * | 2021-04-29 | 2021-07-30 | 博动医学影像科技(上海)有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
CN113616160B (zh) * | 2021-09-14 | 2024-02-06 | 苏州博动戎影医疗科技有限公司 | 基于多模态医学影像的ffr确定方法、装置、设备及介质 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6411459B2 (ja) * | 2013-04-12 | 2018-10-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 冠血流予備量比シミュレーションに対する形状感知超音波プローブ |
US10307056B2 (en) * | 2013-12-05 | 2019-06-04 | Bioptigen, Inc. | Systems and methods for quantitative doppler optical coherence tomography |
US9687147B2 (en) * | 2014-12-02 | 2017-06-27 | Nidek Co., Ltd. | Optical coherence tomography device and control program |
CN105286779B (zh) * | 2015-10-30 | 2017-06-30 | 温州医科大学 | 一种在体视网膜血流动力学的成像与绝对流速测量方法 |
CN105326486B (zh) * | 2015-12-08 | 2017-08-25 | 博动医学影像科技(上海)有限公司 | 血管压力差与血流储备分数的计算方法及系统 |
CN105559810B (zh) * | 2015-12-10 | 2017-08-08 | 博动医学影像科技(上海)有限公司 | 血管单位时间血流量与血流速度的计算方法 |
CN106175700A (zh) * | 2016-09-27 | 2016-12-07 | 广东永士达医疗科技有限公司 | 一种应用于人体开放管腔道的oct探头 |
CN107550519A (zh) * | 2017-08-22 | 2018-01-09 | 深圳先进技术研究院 | 一种多功能血管内超声成像装置 |
CN107730540B (zh) * | 2017-10-09 | 2020-11-17 | 全景恒升(北京)科学技术有限公司 | 基于高精度匹配模型的冠脉参数的计算方法 |
CN107978371B (zh) * | 2017-11-30 | 2021-04-02 | 博动医学影像科技(上海)有限公司 | 快速计算微循环阻力的方法及系统 |
CN109065170B (zh) * | 2018-06-20 | 2021-11-19 | 博动医学影像科技(上海)有限公司 | 获取血管压力差的方法及装置 |
WO2020077381A1 (en) * | 2018-10-17 | 2020-04-23 | The University Of Melbourne | Coronary evaluation system and process |
CN112384138B (zh) * | 2018-10-25 | 2023-07-25 | 苏州润迈德医疗科技有限公司 | 获取心表大动脉的血流量的方法、装置、系统及存储介质 |
CN112155580B (zh) * | 2019-11-20 | 2021-07-13 | 苏州润迈德医疗科技有限公司 | 基于造影图像修正血流速度和微循环参数的方法及装置 |
CN111134651B (zh) * | 2019-12-09 | 2022-03-08 | 杭州脉流科技有限公司 | 基于腔内影像计算血流储备分数的方法、装置、系统以及计算机存储介质 |
CN111134614A (zh) * | 2019-12-20 | 2020-05-12 | 苏州大学 | 基于oct测量眼球血管中血流的绝对速度的方法及系统 |
CN111462117A (zh) * | 2020-05-19 | 2020-07-28 | 深圳北芯生命科技有限公司 | 基于血管图像的数据处理系统及其数据处理方法 |
CN212365043U (zh) * | 2020-05-19 | 2021-01-15 | 深圳北芯生命科技有限公司 | 基于血管图像的数据处理系统 |
CN111754506B (zh) * | 2020-07-01 | 2024-02-06 | 杭州脉流科技有限公司 | 基于腔内影像的冠脉狭窄率计算方法、装置、系统和计算机存储介质 |
CN112450882B (zh) * | 2020-11-20 | 2022-07-29 | 深圳先进技术研究院 | 超声探头、内窥镜、内窥成像系统以及内窥成像方法 |
CN112535466A (zh) * | 2020-12-16 | 2021-03-23 | 成都全景恒升科技有限公司 | 一种基于血管影像的血流储备分数计算方法 |
CN113180631A (zh) * | 2021-04-29 | 2021-07-30 | 博动医学影像科技(上海)有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
-
2021
- 2021-04-29 CN CN202110474542.2A patent/CN113180631A/zh active Pending
-
2022
- 2022-04-27 JP JP2023535452A patent/JP2023553556A/ja active Pending
- 2022-04-27 WO PCT/CN2022/089546 patent/WO2022228464A1/zh active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2022228464A9 (zh) | 2023-12-07 |
CN113180631A (zh) | 2021-07-30 |
WO2022228464A1 (zh) | 2022-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107730540B (zh) | 基于高精度匹配模型的冠脉参数的计算方法 | |
EP2453791B1 (en) | Apparatus for measuring flow and pressure within a vessel | |
EP4026491A1 (en) | Simplified method, apparatus and system for measuring coronary artery vessel evaluation parameters | |
US9526424B2 (en) | Apparatus, systems, and methods of in-vivo blood clearing in a lumen | |
EP2967380B1 (en) | Methods, systems and computer program products for non-invasive determination of blood flow distribution using speckle imaging techniques and hemodynamic modeling | |
WO2017097073A1 (zh) | 血管压力差与血流储备分数的计算方法及系统 | |
US20140276059A1 (en) | Externally imaging a body structure within a patient | |
JP2023553556A (ja) | 血管内イメージングに基づく血流速度、血流予備量比の分析方法 | |
KR20070058523A (ko) | 혈관절편 영상화 방법 및 장치 | |
CN111462117A (zh) | 基于血管图像的数据处理系统及其数据处理方法 | |
US11944778B2 (en) | Methods and systems for automatic pullback trigger | |
Mackle et al. | Fibre optic intravascular measurements of blood flow: A review | |
Kips et al. | Identifying the vulnerable plaque: A review of invasive and non-invasive imaging modalities | |
Taki et al. | Overview of different medical imaging techniques for the identification of coronary atherosclerotic plaques | |
CN212365043U (zh) | 基于血管图像的数据处理系统 | |
Pandian et al. | Intravascular ultrasound estimation of arterial stenosis | |
Li et al. | Temporal averaging for quantification of lumen dimensions in intravascular ultrasound images | |
Sethuraman et al. | Intravascular photoacoustic imaging of atherosclerotic plaques: ex-vivo study using a rabbit model of atherosclerosis | |
Moridani | Detection and Quantification of Coronary Atherosclerotic Plaque Using Different Imaging Modalities | |
JP2024512338A (ja) | 自動引き戻しトリガのシステム及び方法 | |
AU2010298337B2 (en) | Apparatus, systems, and methods of in-vivo blood clearing in a lumen | |
Dimopoulos et al. | Approach to intermediate lesions | |
Porto et al. | Integrated Imaging 11 | |
Fitzgerald et al. | INTRAVASCULAR ULTRASOUND | |
Honda et al. | 19 Intravascular Imagin g Techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230609 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230609 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240221 |