JP7236768B2 - 心筋血流量及びct画像に基づく冠状動脈冠血流予備量比の計算方法 - Google Patents
心筋血流量及びct画像に基づく冠状動脈冠血流予備量比の計算方法 Download PDFInfo
- Publication number
- JP7236768B2 JP7236768B2 JP2021529763A JP2021529763A JP7236768B2 JP 7236768 B2 JP7236768 B2 JP 7236768B2 JP 2021529763 A JP2021529763 A JP 2021529763A JP 2021529763 A JP2021529763 A JP 2021529763A JP 7236768 B2 JP7236768 B2 JP 7236768B2
- Authority
- JP
- Japan
- Prior art keywords
- coronary
- image
- coronary artery
- flow
- aorta
- 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.)
- Active
Links
- 230000002107 myocardial effect Effects 0.000 title claims description 38
- 230000017531 blood circulation Effects 0.000 title claims description 32
- 238000004364 calculation method Methods 0.000 title description 4
- 210000004351 coronary vessel Anatomy 0.000 claims description 84
- 210000000709 aorta Anatomy 0.000 claims description 32
- 238000000034 method Methods 0.000 claims description 32
- 230000000747 cardiac effect Effects 0.000 claims description 31
- 201000000057 Coronary Stenosis Diseases 0.000 claims description 24
- 206010011089 Coronary artery stenosis Diseases 0.000 claims description 23
- 239000008280 blood Substances 0.000 claims description 14
- 210000004369 blood Anatomy 0.000 claims description 14
- 206010020565 Hyperaemia Diseases 0.000 claims description 13
- 230000002966 stenotic effect Effects 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 238000004335 scaling law Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000000295 complement effect Effects 0.000 claims description 5
- 230000000877 morphologic effect Effects 0.000 claims description 5
- 210000004204 blood vessel Anatomy 0.000 claims description 4
- 230000000544 hyperemic effect Effects 0.000 claims description 4
- 210000004165 myocardium Anatomy 0.000 claims description 4
- 230000010412 perfusion Effects 0.000 claims description 4
- 230000036772 blood pressure Effects 0.000 claims description 3
- 238000002592 echocardiography Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 3
- 238000002600 positron emission tomography Methods 0.000 claims description 3
- 238000002603 single-photon emission computed tomography Methods 0.000 claims description 3
- 238000003325 tomography Methods 0.000 claims description 3
- 208000031481 Pathologic Constriction Diseases 0.000 claims description 2
- 230000010339 dilation Effects 0.000 claims description 2
- 208000037804 stenosis Diseases 0.000 claims description 2
- 230000036262 stenosis Effects 0.000 claims description 2
- 210000001367 artery Anatomy 0.000 claims 1
- 230000001746 atrial effect Effects 0.000 claims 1
- 238000013184 cardiac magnetic resonance imaging Methods 0.000 claims 1
- 230000000284 resting effect Effects 0.000 claims 1
- 239000000243 solution Substances 0.000 description 11
- 239000000284 extract Substances 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 230000011218 segmentation Effects 0.000 description 5
- 230000003902 lesion Effects 0.000 description 4
- OIRDTQYFTABQOQ-KQYNXXCUSA-N adenosine Chemical compound C1=NC=2C(N)=NC=NC=2N1[C@@H]1O[C@H](CO)[C@@H](O)[C@H]1O OIRDTQYFTABQOQ-KQYNXXCUSA-N 0.000 description 2
- 210000002376 aorta thoracic Anatomy 0.000 description 2
- 238000002586 coronary angiography Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- XQYZDYMELSJDRZ-UHFFFAOYSA-N papaverine Chemical compound C1=C(OC)C(OC)=CC=C1CC1=NC=CC2=CC(OC)=C(OC)C=C12 XQYZDYMELSJDRZ-UHFFFAOYSA-N 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 229930008281 A03AD01 - Papaverine Natural products 0.000 description 1
- ZKHQWZAMYRWXGA-KQYNXXCUSA-N Adenosine triphosphate Chemical compound C1=NC=2C(N)=NC=NC=2N1[C@@H]1O[C@H](COP(O)(=O)OP(O)(=O)OP(O)(O)=O)[C@@H](O)[C@H]1O ZKHQWZAMYRWXGA-KQYNXXCUSA-N 0.000 description 1
- ZKHQWZAMYRWXGA-UHFFFAOYSA-N Adenosine triphosphate Natural products C1=NC=2C(N)=NC=NC=2N1C1OC(COP(O)(=O)OP(O)(=O)OP(O)(O)=O)C(O)C1O ZKHQWZAMYRWXGA-UHFFFAOYSA-N 0.000 description 1
- 239000002126 C01EB10 - Adenosine Substances 0.000 description 1
- 229960005305 adenosine Drugs 0.000 description 1
- 230000004872 arterial blood pressure Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000005003 heart tissue Anatomy 0.000 description 1
- 238000002608 intravascular ultrasound Methods 0.000 description 1
- 238000010253 intravenous injection Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 229960001789 papaverine Drugs 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000250 revascularization Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000000779 thoracic wall Anatomy 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/155—Segmentation; Edge detection involving morphological operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20156—Automatic seed setting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/404—Angiography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Physiology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Cardiology (AREA)
- Vascular Medicine (AREA)
- Hematology (AREA)
- Pulmonology (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法は以下のステップを含む:
S01:心臓CT画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得て、心筋画像に基づいて心筋体積を特定し、
S02:大動脈画像を処理して全大動脈相補画像を得て、領域成長を行い、冠状動脈口を含む大動脈画像を得て、冠状動脈口を含む大動脈画像と全大動脈相補画像に基づき、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に基づいて冠状動脈口を特定し、
S03:心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出し、
S04:冠状動脈画像を二値化し、等値面画像を描き、冠状動脈三次元格子画像を得て、
S05:最大充血状態における冠動脈入口箇所の総流量Qtotal=心筋体積×心筋血流量×CFRを計算し、CFRは冠状動脈冠血流予備能(coronary flow reserve:CFR)であり、
S06:充血状態における血流速度V1を計算し、
S07:V1を冠動脈狭窄血管の入口流速とし、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔP、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算し、但し、Paは大動脈平均圧であり、冠血流予備量比FFR=Pd/Paを得る。
S61:流量体積スケーリング則及び心臓CT三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流量Q=Qtotal×(V/Vtotal)3/4を特定し、但し、Vtotalは心臓CT三次元再構成のすべての心臓表面冠状動脈の血管体積の和であり、Vは心臓表面冠状動脈ツリー内の任意の一本の血管及びその下流血管中の血管体積の和であり、
S62:流量体積スケーリング則及び心臓CT三次元再構成の心臓表面冠状動脈ツリーに基づき、ツリー内の任意の一本の血管内の血流速度V1=Q/Dを特定し、但し、Dは当該血管の平均断面積である、ことを含む。
血管三次元格子に対して解を求め、数値解法を用いて連続性及びナビエ-ストークス(Navier-Stokes)方程式の解を求め:
P、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速V1であり、
三次元計算流体力学により各冠動脈狭窄の圧力損失ΔP1、ΔP2、ΔP3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔP=ΣΔPi(i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算する、但し、Paは大動脈平均であることを含む。
CT再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし、二次元軸対称モデルを構築し、二次元格子を分割し、数値解法を用いて連続性及びナビエ-ストークス方程式の解を求め:
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速V1であり、
二次元流体力学により各冠動脈狭窄の圧力損失ΔP1、ΔP2、ΔP3・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔP=ΣΔPi(i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算する、但し、Paは大動脈平均圧である、ことを含む。
心臓CT画像を分割し、形態学操作により心臓画像を得て、当該心臓画像をヒストグラム分析して心室心房画像を得て、心臓画像と心室心房画像の差分により心筋画像を得る。図2参照。
大動脈画像の二値化画像に対して形態学膨張を行い、全大動脈の二値画像を得て、そして画素リバースにより全大動脈相補画像を得る。
大動脈中心線上点の平均グレースケールに基づき領域成長を行い、冠状動脈口を含む大動脈画像を得る。図3参照。
冠状動脈口を含む大動脈画像と全大動脈相補画像により画像乗算を行い、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に対して連結領域分析を行い、異なるグレースケールラベルで各連結領域をマークし、冠状動脈口を特定する。図4参照。
心筋画像において冠状動脈口をシードポイントとし、領域成長により冠状動脈を抽出し、冠状動脈の平均グレースケール及び平均分散を計算し、冠動脈グレースケール分布に基づき、冠状動脈方向に沿って冠動脈ツリーを抽出する。図5参照。
ステップ三により、冠状動脈画像データV1を得て、当該データ中のボクセルは空間上に一つの立方体を構成し、冠状動脈部分に属すボクセル画素値は0ではなく(画素値約-3000~3000の間)、残りのボクセル画素値はすべて0である。
本ステップにおいて、ステップ五におけるFFR計算のために、データを空間三次元格子データV3に変える必要がある。
冠状動脈画像データV1中のボクセルをトラバーサルし、簡単な画素値判断を行い、もし画素A1が0に等しい場合、当該画素値は変わらず、もしA1が0に等しくない場合、A1の画素値を1とする。
最終的に一つの新しい画像データV2が得られ、当該画像において、冠状動脈部分に属するボクセル画素値は1であり、残りの部分は0である。
ボクセルは一つの極小の六面体であると定義され、隣接する上下層の間の四つの画素が立方体上の八つの頂点を構成する。等値面とは空間中のある同じ属性値を有するすべての点の集合である。これは次のように表すことができる:
{(x、y、z)│f(x、y、z)=c}、cは定数である
本方法におけるcは三次元再構成過程において与えられた画素値1である。
(1)原始データを前処理した後、特定の配列に読み込む。
(2)格子データボディから一つのユニットボディを抽出して現在のユニットボディとすると同時に当該ユニットボディのすべての情報を得る。
(3)現在のユニットボディの8個の頂点の函数値と所定の等値面値Cを比較し、当該ユニットボディの状態表を得る。
(4)現在のユニットボディの状態表インデックスから、等値面と交わるユニットボディエッジを探し出し、そして線形補間の方法により各交点の位置座標を計算する。
(5)中心差分法により現在のユニットボディの8個の頂点の法線ベクトルを求め、それから線形補間の方法により三角パッチの各頂点の法線を得る。
(6)各三角パッチの頂点座標及び頂点の法線ベクトルに基づき等値面イメージの作成を行う。
最終的に冠状動脈の三次元格子画像データV3を得る。図6参照。
心筋コントラストエコー法(MCE)または単一光子放射断層撮影法(SPECT)または陽電子放出断層撮影法(PET)または心臓核磁気共鳴画像法(MRI)またはCT灌流等非侵襲的測定により、安静状態心筋血流量及び冠状動脈冠血流予備能(CFR)を特定し、心筋体積、心筋血流量、CFRにより、最大充血状態における冠動脈入口箇所(左冠動脈ツリー及び右冠動脈ツリーの和を含む)の総流量Qtotal=心筋体積×心筋血流量×CFRを計算する。
V1を冠動脈狭窄血管の入り口流速とし、計算流体力学(CFD)方法により各冠動脈狭窄の厚力損失ΔP1、ΔP2、ΔP3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔP=ΣΔPi(i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算し、但し、Paは大動脈平均圧であり、最後に式FFR=Pd/Paにより冠血流予備量比を計算する。
CT再構成のジオメトリ構造に基づき、三次元格子を分割し、数値解法(例えば、有限差分、有限要素法、有限体積法等)により連続性及びナビエ-ストークス(Navier-Stokes)方程式の解を求め:
P、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速V1であり、
式[A1]及び[A2]に基づき、三次元CFDを実行して各冠動脈狭窄の圧力損失ΔP1、ΔP2、ΔP3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔP=ΣΔPi(i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算し、但し、Paは大動脈平均圧である。
CT再構成のジオメトリ構造に基づき、狭窄のある血管を真っすぐに引き伸ばし(二次元軸対称モデル)、二次元格子を分割し、数値解法(例えば、有限差分、有限要素法、有限体積法等)により連続性及びナビエ-ストークス方程式の解を求め:
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速V1であり、
式[A3]~[A5]に基づき、二次元CFDを実行し、各冠動脈狭窄の圧力損失ΔP1、ΔP2、ΔP3等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔP=ΣΔPi(i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算し、但し、Paは大動脈平均圧である。
Claims (7)
- S01:心臓を含むCT画像を分割し、形態学操作により心臓のみを含む心臓画像を得て、前記心臓画像をヒストグラム分析して心室心房のみを含む心室心房画像を得て、前記心臓画像と前記心室心房画像の差分により心筋を含む画像を得て、前記心筋を含む画像に基づいて心筋体積を特定し、
S02:前記心臓画像から心臓の大動脈画像を得て、前記大動脈画像の二値化画像に対して形態学膨張を行い、心臓の全大動脈の二値画像を得て、そして画素リバースにより心臓の全大動脈相補画像を得て、前記大動脈画像において前記大動脈の中心線上点の平均グレースケールに基づき領域成長を行い、冠状動脈口を含む大動脈画像を得て、前記冠状動脈口を含む大動脈画像と前記全大動脈相補画像に基づき、冠状動脈口を含む画像を得て、冠状動脈口を含む画像に基づいて前記大動脈における前記冠状動脈口の位置を特定し、
S03:前記心筋を含む画像において前記冠状動脈口をシードポイントとし、領域成長により前記冠状動脈を抽出し、前記冠状動脈の平均グレースケールを計算し、前記冠状動脈のグレースケール分布及び前記冠状動脈の平均グレースケールを比較し、前記冠状動脈の成長方向に沿って冠状動脈ツリーを抽出し、
S04:前記冠状動脈ツリーの画像を二値化し、二値化したデータに基づいて等値面画像を三次元座標系に描き、最終的に前記冠状動脈の三次元格子画像を得て、
S05:最大充血状態における前記冠状動脈口の総流量Qtotal=心筋体積×安静状態心筋血流量×CFRを計算し、CFRは冠状動脈冠血流予備能(coronary flow reserve:CFR)であり、
S06:最大充血状態における血流速度V1を計算し、
S07:V1を冠状動脈狭窄血管の入口流速とし、冠状動脈入口から冠状動脈狭窄遠位までの圧力損失ΔP、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算し、但し、Paは心臓の大動脈平均圧であり、冠血流予備量比FFR=Pd/Paを得る、
ステップを含み、
前記ステップS06は、
S61:流量体積スケーリング則及び三次元再構成の前記冠状動脈ツリーに基づき、前記冠状動脈ツリー内の任意の一本の血管内の血流量Q=Q total ×(V/V total ) 3/4 を特定し、但し、Q total は前記最大充血状態における前記冠状動脈口の総流量であり、V total は三次元再構成の前記冠状動脈ツリーの中のすべての冠状動脈の血管体積の和であり、Vは前記冠状動脈ツリー内の任意の一本の血管及びその下流血管中の血管体積の和であり、
S62:流量体積スケーリング則及び三次元再構成の前記冠状動脈ツリーに基づき、前記冠状動脈ツリー内の任意の一本の血管内の血流速度V 1 =Q/Dを特定し、但し、Dは当該血管の平均断面積である、
ことを含み、
前記ステップS07は具体的に、
血管三次元格子に対して解を求め、数値解法を用いて連続性及びナビエ-ストークス(Navier-Stokes)方程式の解を求め:
P、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠状動脈狭窄血管の入口流速V 1 であり、
三次元計算流体力学により各冠状動脈狭窄の圧力損失ΔP 1 、ΔP 2 、ΔP 3 ・・・、冠状動脈入口から冠状動脈狭窄遠位までの圧力損失ΔP=ΣΔP i (i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算し、但し、Paは大動脈平均である、ことを含む、
ことを特徴とする心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。 - 前記ステップS02において、冠状動脈口を含む画像を得た後、前記冠状動脈口を含む画像に対して連結領域分析を行い、異なるグレースケールラベルで各連結領域をマークし、前記大動脈における前記冠状動脈口の位置を特定する、ことを特徴とする請求項1に記載の心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。
- 前記ステップS02において、前記心臓画像において、大動脈断面が円状になる特徴を利用し、心臓の上行大動脈及びその中心線を抽出し、心臓の大動脈画像を得る、ことを特徴とする請求項1に記載の心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。
- 前記ステップS04において冠状動脈画像二値化は、
前記冠状動脈画像V1中のボクセルをトラバーサルし、もしボクセル画素が0に等しい場合、当該画素値は変わらず、もし0に等しくない場合、画素値を1とし、これらの画素値は一つの新しいデータV2となることを含む、
ことを特徴とする請求項1に記載の心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。 - 前記ステップS05において心筋コントラストエコー法(MCE)または単一光子放射断層撮影法(SPECT)または陽電子放出断層撮影法(PET)または心臓核磁気共鳴画像法(MRI)またはCT灌流により、安静状態心筋血流量及び冠状動脈冠血流予備能(CFR)を特定する、
ことを特徴とする請求項1に記載の心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。 - 前記ステップS07は、
CT再構成の前記冠状動脈の三次元格子画像に基づき、狭窄のある血管を真っすぐに引き伸ばし、二次元軸対称モデルを構築し、二次元格子を分割し、数値解法を用いて連続性及びナビエ-ストークス方程式の解を求め:
入口境界条件は、最大充血状態における冠状動脈狭窄血管の入口流速V1であり、
二次元流体力学により各冠状動脈狭窄の圧力損失ΔP1、ΔP2、ΔP3・・・、冠状動脈入口から冠状動脈狭窄遠位までの圧力損失ΔP=ΣΔPi(i=1、2、3・・・)、狭窄遠位冠状動脈内平均圧Pd=Pa-ΔPを計算する、但し、Paは大動脈平均圧である、ことを含む、ことを特徴とする請求項1に記載の心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。 - 前記ステップS07はさらに、
前記冠状動脈の異なるタイプの弯曲に対して、前記冠状動脈の三次元格子画像で冠状動脈入口から冠状動脈狭窄遠位までの圧力損失を計算し、前記二次元軸対称モデルに照らして計算し、各タイプの弯曲の二次元軸対称結果に対する修正係数を記憶するためのデータベースを構築し、
圧力損失を得た後、前記データベース中の修正係数に照らし、修正後の冠状動脈入口から冠状動脈狭窄遠位までの圧力損失を得て、その後前記FFRを計算する、
ことを含む、ことを特徴とする請求項6に記載の心筋血流量及びCT画像に基づく冠状動脈冠血流予備量比の計算方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811438744.6 | 2018-11-28 | ||
CN201811438744.6A CN111227822B (zh) | 2018-11-28 | 2018-11-28 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
PCT/CN2019/071203 WO2020107667A1 (zh) | 2018-11-28 | 2019-01-10 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2022508238A JP2022508238A (ja) | 2022-01-19 |
JP7236768B2 true JP7236768B2 (ja) | 2023-03-10 |
Family
ID=70853702
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021529763A Active JP7236768B2 (ja) | 2018-11-28 | 2019-01-10 | 心筋血流量及びct画像に基づく冠状動脈冠血流予備量比の計算方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11896416B2 (ja) |
JP (1) | JP7236768B2 (ja) |
CN (1) | CN111227822B (ja) |
WO (1) | WO2020107667A1 (ja) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109805949B (zh) * | 2019-03-19 | 2020-05-22 | 苏州润迈德医疗科技有限公司 | 基于压力传感器和造影图像计算血流储备分数的方法 |
WO2020210948A1 (zh) * | 2019-04-15 | 2020-10-22 | 博动医学影像科技(上海)有限公司 | 获取冠脉血流量及血流速度的方法和装置 |
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 |
CN111754522B (zh) * | 2020-06-19 | 2021-08-03 | 上海杏脉信息科技有限公司 | 冠状动脉血流动力学数据的获取方法及装置 |
CN111815589B (zh) * | 2020-06-29 | 2022-08-05 | 苏州润迈德医疗科技有限公司 | 基于ct序列图像获取无干扰冠脉树图像的方法和系统 |
CN111815585B (zh) * | 2020-06-29 | 2022-08-05 | 苏州润迈德医疗科技有限公司 | 基于ct序列图像获取冠脉树和冠脉入口点的方法和系统 |
CN111815586B (zh) * | 2020-06-29 | 2022-08-05 | 苏州润迈德医疗科技有限公司 | 基于ct图像获取左心房、左心室的连通域的方法和系统 |
CN111652881A (zh) * | 2020-07-01 | 2020-09-11 | 杭州脉流科技有限公司 | 基于深度学习的冠脉重构和血流储备分数计算方法、装置、设备以及可读存储介质 |
CN112704505B (zh) * | 2020-11-20 | 2022-05-24 | 杭州阿特瑞科技有限公司 | 一种利用cta和dsa测量冠状动脉血流储备分数的方法 |
CN112535466A (zh) * | 2020-12-16 | 2021-03-23 | 成都全景恒升科技有限公司 | 一种基于血管影像的血流储备分数计算方法 |
CN112837306B (zh) * | 2021-02-20 | 2022-11-22 | 薛竟宜 | 基于深度学习和中智理论的冠状动脉病变功能学定量方法 |
KR102460948B1 (ko) * | 2021-03-19 | 2022-10-28 | 제주대학교 산학협력단 | 관상동맥 협착의 중증도 평가에 필요한 정보를 제공하는 방법 |
CN113096141B (zh) * | 2021-04-19 | 2022-01-11 | 推想医疗科技股份有限公司 | 冠状动脉分割方法以及冠状动脉分割装置 |
CN113139959B (zh) * | 2021-05-17 | 2021-10-01 | 北京安德医智科技有限公司 | 获得心肌桥图像的方法及装置、电子设备和存储介质 |
CN114298988B (zh) * | 2021-12-17 | 2022-11-15 | 深圳市阅影科技有限公司 | 计算血管功能学指标的方法、装置、存储介质及处理器 |
CN114882099A (zh) * | 2022-04-22 | 2022-08-09 | 中国人民解放军陆军第九五〇医院 | 冠脉血管供血量分析装置及运行方法 |
CN116051458A (zh) * | 2022-11-23 | 2023-05-02 | 杭州脉流科技有限公司 | 基于术前冠脉造影影像评估术后冠脉血流储备分数的方法和计算机设备 |
CN116188336B (zh) * | 2023-04-17 | 2023-12-22 | 柏意慧心(杭州)网络科技有限公司 | 基于血管造影的心肌病形态学计算方法、装置及存储介质 |
CN117036530B (zh) * | 2023-10-10 | 2023-12-08 | 杭州脉流科技有限公司 | 基于跨模态数据的冠状动脉血流储备分数预测方法及装置 |
CN117197164B (zh) * | 2023-11-08 | 2024-03-08 | 中国医学科学院北京协和医院 | 用于计算心肌血管供血区域的管道流域计算方法及系统 |
CN117476238B (zh) * | 2023-11-08 | 2024-05-28 | 中国医学科学院北京协和医院 | 用于评估心肌缺血的计算流体力学仿真分析方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140341453A1 (en) | 2013-05-17 | 2014-11-20 | Bailing Hsu | Myocardial blood flow quantitation with dynamic spect or spect/ct imaging |
JP2016511649A (ja) | 2013-01-15 | 2016-04-21 | キャスワークス・リミテッドCathWorks Ltd. | 血流予備量比の算出 |
CN106023202A (zh) | 2016-05-20 | 2016-10-12 | 苏州润心医疗科技有限公司 | 基于心脏ct图像的冠状动脉血流储备分数计算方法 |
US20170095292A1 (en) | 2015-10-02 | 2017-04-06 | Heartflow, Inc. | System and method for diagnosis and assessment of cardiovascular disease by comparing arterial supply capacity to end-organ demand |
JP2017070742A (ja) | 2015-10-05 | 2017-04-13 | 東芝メディカルシステムズ株式会社 | 画像処理装置、画像処理方法、及びプログラム |
JP2017512577A (ja) | 2014-03-31 | 2017-05-25 | ハートフロー, インコーポレイテッド | 流量比を用いて血流の特徴を決定するシステム及び方法 |
CN108550189A (zh) | 2018-05-03 | 2018-09-18 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9974508B2 (en) * | 2011-09-01 | 2018-05-22 | Ghassan S. Kassab | Non-invasive systems and methods for determining fractional flow reserve |
US9814433B2 (en) * | 2012-10-24 | 2017-11-14 | Cathworks Ltd. | Creating a vascular tree model |
US9757073B2 (en) * | 2012-11-06 | 2017-09-12 | Koninklijke Philips N.V. | Fractional flow reserve (FFR) index |
US10052031B2 (en) * | 2013-03-04 | 2018-08-21 | Siemens Healthcare Gmbh | Determining functional severity of stenosis |
WO2015164086A1 (en) * | 2014-04-22 | 2015-10-29 | Siemens Aktiengesellschaft | Method and system for hemodynamic computation in coronary arteries |
CN103932694A (zh) * | 2014-05-07 | 2014-07-23 | 霍云龙 | 精确诊断心肌血流储备分数(ffr)的方法和设备 |
CN107427268B (zh) * | 2014-11-14 | 2023-07-28 | 西门子保健有限责任公司 | 用于基于纯几何机器学习的血流储备分数的方法和系统 |
JP6692809B2 (ja) * | 2014-11-14 | 2020-05-13 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 経皮的冠動脈インターベンション計画インタフェース、並びに関連するデバイス、システム、及び方法 |
CN105078425B (zh) * | 2015-09-09 | 2016-06-08 | 苏州润心医疗科技有限公司 | 冠状动脉负荷检测系统及检测方法 |
JP6484760B2 (ja) * | 2015-11-05 | 2019-03-13 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 非侵襲的血流予備量比(ffr)に対する側副血流モデル化 |
CN105326486B (zh) * | 2015-12-08 | 2017-08-25 | 博动医学影像科技(上海)有限公司 | 血管压力差与血流储备分数的计算方法及系统 |
CN105559810B (zh) * | 2015-12-10 | 2017-08-08 | 博动医学影像科技(上海)有限公司 | 血管单位时间血流量与血流速度的计算方法 |
CN106327487B (zh) * | 2016-08-18 | 2018-01-02 | 苏州润迈德医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
CN106473731A (zh) * | 2016-10-25 | 2017-03-08 | 北京工业大学 | 基于个性化冠状动脉分支血流量的ffrct计算方法 |
CN108109698B (zh) * | 2016-12-28 | 2021-04-20 | 北京科亚方舟医疗科技股份有限公司 | 计算血流储备分数的系统和设置边界条件的方法 |
CN107978371B (zh) * | 2017-11-30 | 2021-04-02 | 博动医学影像科技(上海)有限公司 | 快速计算微循环阻力的方法及系统 |
CN108186038B (zh) * | 2018-02-11 | 2020-11-17 | 杭州脉流科技有限公司 | 基于动脉造影影像计算冠脉血流储备分数的系统 |
CN108511075B (zh) * | 2018-03-29 | 2022-10-25 | 杭州脉流科技有限公司 | 一种非侵入式获取血流储备分数的方法和系统 |
CN108735270A (zh) * | 2018-05-25 | 2018-11-02 | 杭州脉流科技有限公司 | 基于降维模型的血流储备分数获取方法、装置、系统和计算机存储介质 |
-
2018
- 2018-11-28 CN CN201811438744.6A patent/CN111227822B/zh active Active
-
2019
- 2019-01-10 WO PCT/CN2019/071203 patent/WO2020107667A1/zh active Application Filing
- 2019-01-10 JP JP2021529763A patent/JP7236768B2/ja active Active
-
2021
- 2021-05-24 US US17/328,550 patent/US11896416B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016511649A (ja) | 2013-01-15 | 2016-04-21 | キャスワークス・リミテッドCathWorks Ltd. | 血流予備量比の算出 |
US20140341453A1 (en) | 2013-05-17 | 2014-11-20 | Bailing Hsu | Myocardial blood flow quantitation with dynamic spect or spect/ct imaging |
JP2017512577A (ja) | 2014-03-31 | 2017-05-25 | ハートフロー, インコーポレイテッド | 流量比を用いて血流の特徴を決定するシステム及び方法 |
US20170095292A1 (en) | 2015-10-02 | 2017-04-06 | Heartflow, Inc. | System and method for diagnosis and assessment of cardiovascular disease by comparing arterial supply capacity to end-organ demand |
JP2017070742A (ja) | 2015-10-05 | 2017-04-13 | 東芝メディカルシステムズ株式会社 | 画像処理装置、画像処理方法、及びプログラム |
CN106023202A (zh) | 2016-05-20 | 2016-10-12 | 苏州润心医疗科技有限公司 | 基于心脏ct图像的冠状动脉血流储备分数计算方法 |
CN108550189A (zh) | 2018-05-03 | 2018-09-18 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2022508238A (ja) | 2022-01-19 |
US11896416B2 (en) | 2024-02-13 |
CN111227822B (zh) | 2022-02-11 |
CN111227822A (zh) | 2020-06-05 |
US20210275124A1 (en) | 2021-09-09 |
WO2020107667A1 (zh) | 2020-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7236768B2 (ja) | 心筋血流量及びct画像に基づく冠状動脈冠血流予備量比の計算方法 | |
JP7236769B2 (ja) | 心筋血流量及びct画像に基づく微小循環抵抗指数の計算方法 | |
WO2019210553A1 (zh) | 基于造影图像和流体力学模型的微循环阻力指数计算方法 | |
US10460204B2 (en) | Method and system for improved hemodynamic computation in coronary arteries | |
CN106023202B (zh) | 基于心脏ct图像的冠状动脉血流储备分数计算方法 | |
JP6611959B2 (ja) | 血管圧力差の計算方法およびシステム | |
JP5953438B2 (ja) | 冠状動脈を通る血流量を決定する方法 | |
US10803995B2 (en) | Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in diseased and hypothetical normal anatomical models | |
JP6484760B2 (ja) | 非侵襲的血流予備量比(ffr)に対する側副血流モデル化 | |
JP6749917B2 (ja) | iFR−CT | |
Glaßer et al. | Combined visualization of wall thickness and wall shear stress for the evaluation of aneurysms | |
CN110706770B (zh) | 心脏数据处理设备及处理方法、计算机可读存储介质 | |
CN112690814B (zh) | 一种低误差的冠状动脉血流储备分数测量方法 | |
US10332255B2 (en) | Method for assessing stenosis severity in a lesion tree through stenosis mapping | |
CN115115735A (zh) | 一种基于多期相冠脉ct造影的内皮动应变的快速计算系统及方法 | |
US20230410307A1 (en) | Method and system for visualization | |
Siogkas et al. | CORONARY ATHEROSCLEROSIS ASSESSMENT: A NEW ANATOMICAL, FUNCTIONAL, MORPHOLOGICAL AND BIO-MECHANICAL APPROACH | |
Maklad et al. | Extraction of liver volumetry based on blood vessel from the portal phase CT dataset |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210525 |
|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7426 Effective date: 20210604 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220526 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220621 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20220921 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20221012 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230131 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230220 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7236768 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |