CN110522439A - 测量冠状动脉血管评定参数的简化方法、装置及系统 - Google Patents

测量冠状动脉血管评定参数的简化方法、装置及系统 Download PDF

Info

Publication number
CN110522439A
CN110522439A CN201910834871.6A CN201910834871A CN110522439A CN 110522439 A CN110522439 A CN 110522439A CN 201910834871 A CN201910834871 A CN 201910834871A CN 110522439 A CN110522439 A CN 110522439A
Authority
CN
China
Prior art keywords
coronary artery
coronary
measurement
assessment parameters
contrastographic picture
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
CN201910834871.6A
Other languages
English (en)
Inventor
刘广志
王之元
冯亮
李泽华
霍勇
龚艳君
李建平
易铁慈
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.)
Suzhou Yun Medsphere Co Ltd
Suzhou Rainmed Medical Technology Co Ltd
Original Assignee
Suzhou Yun Medsphere Co Ltd
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
Priority claimed from CN201811093192.XA external-priority patent/CN109363651A/zh
Priority claimed from CN201910206541.2A external-priority patent/CN109770888A/zh
Application filed by Suzhou Yun Medsphere Co Ltd filed Critical Suzhou Yun Medsphere Co Ltd
Priority to JP2022513657A priority Critical patent/JP7457105B2/ja
Priority to PCT/CN2019/115027 priority patent/WO2021042477A1/zh
Priority to EP19943973.8A priority patent/EP4026491A4/en
Publication of CN110522439A publication Critical patent/CN110522439A/zh
Priority to US17/683,490 priority patent/US20220183643A1/en
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/504Clinical applications involving diagnosis of blood vessels, e.g. by angiography
    • 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/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • 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
    • 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
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/507Clinical applications involving determination of haemodynamic parameters, e.g. perfusion CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M5/00Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests
    • A61M5/007Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests for contrast media
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • G06T7/0016Biomedical image inspection using an image reference approach involving temporal comparison
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/215Motion-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • 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/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices 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/5235Devices 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
    • A61B6/5241Devices 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 combining overlapping images of the same imaging modality, e.g. by stitching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20036Morphological image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30021Catheter; Guide wire
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104Vascular flow; Blood flow; Perfusion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30172Centreline of tubular or elongated structure
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Cardiology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Dentistry (AREA)
  • Quality & Reliability (AREA)
  • Hematology (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Multimedia (AREA)
  • Anesthesiology (AREA)
  • Human Computer Interaction (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本申请提供了一种测量冠状动脉血管评定参数的简化方法、装置及系统,测量方法包括:对测量血管进行冠脉造影;通过压力导丝测量冠状动脉狭窄远端的压力Pd;选取测量血管的第一体位造影图像和第二体位造影图像;根据第一体位造影图像和第二体位造影图像三维建模获得冠状动脉三维血管模型;根据冠状动脉三维血管模型,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2;根据Pd、T1、T2测量冠状动脉血管评定参数。本申请降低了扩张药物的注射时间,基于冠脉造影图像测量IMR、CFR等冠状动脉血管评定参数,测量过程简单,测试结果准确。

Description

测量冠状动脉血管评定参数的简化方法、装置及系统
技术领域
本发明涉及冠状动脉医学技术领域,特别是涉及一种测量冠状动脉血管评定参数的简化方法、装置、冠状动脉分析系统及计算机存储介质。
背景技术
心脏属于高耗能器官。静息状态下,心肌代谢的摄氧量即可达血液氧含量的60%~80%。因此,在运动等应激状态下,心脏难以通过提高组织的摄氧能力来满足心肌缺氧量增加的需求,而绝大部分是通过增加心肌血流量来保证心肌代谢的氧需要量。心肌微循环占冠状动脉循环组成的95%,通过局部代谢产物、内皮、神经内分泌、肌源性等各种因素发挥调控心肌血流量的作用。研究显示,冠状动脉微循环功能异常是冠心病患者远期预后不良的重要预测因素。
在2013年指南进行了变迁,指出“对于怀疑存在微血管性心绞痛的患者,如果冠脉造影未见明显异常,可以考虑在造影期间腔内注射乙酰胆碱或腺苷进行多普勒测量,计算内皮依赖或非内皮以来的CFR,明确是否存在微循环/心外膜血管的痉挛”,并将此列为IIB类的推荐。
2019年指南增加了1个IIA类推荐和2个IIB类推荐。提出“对于持续存在症状的患者但冠脉造影正常或中等狭窄且iwfr/FFR值保留的患者,应当考虑使用基于导丝测量的CFR和/或微循环阻力测量”,并列为IIA类的推荐。
冠状动脉微血管功能是通过检测微血管对于血管扩张剂的反应来完成的。指南这两个方面的变迁,也提示冠状动脉微血管功能检查的重要性。冠状动脉微血管功能采用的测量指标是指冠状动脉微血管最大扩张程度即冠脉血流储备(Coronary Flow Reserve)CFR,所使用血管扩张剂主要包括作用于血管平滑肌的非内皮依赖性血管扩张剂以及作用于血管内皮细胞的内皮依赖性血管扩张剂,其中包括腺苷和乙酰胆碱。
对于冠状造影未见明显狭窄,但怀疑冠心病(CAD)的患者,我们既往的检查手段是注射腺苷和乙酰胆碱,来检测微血管对于血管扩张剂的反应。现在的检查方法主要包括冠状动脉血流储备分数(FFR)和微循环阻力指数(IMR),IMR通过软压力导丝同步记录冠脉压力和温度,导丝杆上的两个温度感受器探测到温度变化的时间差就可知道盐水从指引导管到达导丝头端温度感受器运行的平均传导时间(transit mean time,Tmn),根据定义冠脉远端的压力Pd与Tmn的乘积就可得出IMR值。但是总体来说现在评估微循环的方法并不多。现有的检查手段简化了流程、提高了安全性,对结果也进行了优化,所以指南的推荐级别较以前有所提升。除此之外,无创检查包括经胸多普勒超声、放射性核素显像技术、核磁共振成像技术等手段在微循环疾病的诊断中具有价值,但是都有不同程度的不足,未能成为微循环功能评估的推荐方法。
现有的CFR测量方法包括:(1)多普勒导丝测量方法,将多普勒导丝送入冠状动脉血管内(病变远端)直接测量静息和最大充血状态下冠状动脉内血流速度,即可计算出CFR。(2)热稀释曲线测量方法,通过双重感应导丝崁有温度-压力感受器,能直接感受冠状动脉内温度的变化,可获得静息和最大充血状态下的冠状动脉内热稀释曲线,采用血流平均传导时间替代冠状动脉流量速度来计算CFR。
通过压力导丝传感器测量IMR和CFR存在如下问题:(1)压力导丝传感器离冠脉口太近测量的Tmn太小导致IMR结果偏小。太远测量的Tmn太大导致IMR结果偏大;(2)在静息态和最大充血态时共要注射6次生理盐水,压力导丝传感器位置如果有移动会导致每次测量的结果不具备可比性,且测量过程繁琐;(3)每次注射盐水可能得到的Tmn差异较大,如果某次数值与其他2个数值相差超过30%,需要再次注射盐水测量,增加盐水注射次数;(4)压力导丝感受器测量注射盐水温度下降不够迅速会导致无法记录数值,要提高注射速度,提高注射量,用更低温度的盐水。影响因素太多;(5)注射后温度没有足够快恢复原始值也会出错,从开始注射,到温度恢复正常时间太长(>0.6秒);可能是注射太慢,注射速度不均匀,注射量太大等。因此,压力导丝感受器的距离、注射盐水速度、注射量、盐水的温度均会直接影响测量结果,导致结果不准确,测量过程繁琐;且长时间持续注射扩张药物,对病人有较大影响产生严重的不适感。
发明内容
本发明提供了一种测量冠状动脉血管评定参数的简化方法、装置、冠状动脉分析系统及计算机存储介质,以解决现有技术中通过压力导丝测量CFR和IMR时存在长时间持续注射扩张药物,对病人有较大影响产生严重的不适感,以及压力导丝测量过程繁琐,测量结果不准确的问题。
为实现上述目的,第一方面,本申请提供了一种测量冠状动脉血管评定参数的简化方法,包括:
对测量血管进行冠脉造影;
通过压力导丝测量冠状动脉狭窄远端的压力Pd
选取所述测量血管的静息态下的第一体位造影图像和扩张状态下的第二体位造影图像;
选取从冠脉狭窄病变近端到远端的一段血管,根据所述第一体位造影图像和所述第二体位造影图像三维建模获得冠状动脉三维血管模型;
注射造影剂,根据所述造影剂在所述冠状动脉三维血管模型内的流动,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2
根据Pd、T1、T2测量冠状动脉血管评定参数。
可选地,上述的测量冠状动脉血管评定参数的简化方法,所述冠状动脉血管评定参数包括:冠脉血流储备CFR、微循环阻力系数IMR。
可选地,上述的测量冠状动脉血管评定参数的简化方法,所述第一体位图像为静息态下的造影图像,所述第二体位图像为扩张态下的造影图像,即:
所述CFR=T1/T2;和/或
所述IMR=Pd×T2
可选地,上述的测量冠状动脉血管评定参数的简化方法,所述冠状动脉血管评定参数包括:冠状动脉血流储备分数FFR,通过所述压力导丝测量冠状动脉入口压力Pa,所述FFR=Pd/Pa
可选地,上述的测量冠状动脉血管评定参数的简化方法,所述第一体位与所述第二体位的夹角大于30°。
可选地,上述的测量冠状动脉血管评定参数的简化方法,在所述通过压力导丝测量冠状动脉狭窄远端的压力Pd之前,包括:向血管内注射扩张药物。
可选地,上述的测量冠状动脉血管评定参数的简化方法,在所述通过压力导丝测量冠状动脉狭窄远端的压力Pd之后,所述对测量血管进行冠脉造影之前,包括:在注射扩张药物的同时向所述血管内注射造影剂。
可选地,上述的测量冠状动脉血管评定参数的简化方法,所述根据所述第一体位造影图像和所述第二体位造影图像三维建模获得冠状动脉三维血管模型包括:
去除第一体位造影图像和所述第二体位造影图像的干扰血管,得到结果图像;
沿着所述冠状动脉的延伸方向,提取每幅所述结果图像的冠脉中心线和直径;
将每根所述冠脉中心线和直径均投射于三维空间上进行三维建模,获得冠状动脉三维血管模型。
可选地,上述的测量冠状动脉血管评定参数的简化方法,所述时间T1和T2根据心跳周期区域被分成的局部区域图像的帧数与每秒传输帧数的比值计算。
第二方面,本申请提供了一种测量冠状动脉血管评定参数的装置,用于上述的测量冠状动脉血管评定参数的简化方法,包括:压力导丝测量单元、提取冠脉造影单元、三维建模单元和参数测量单元,所述提取冠脉造影单元与三维建模单元连接,所述参数测量单元与所述压力导丝测量单元、所述三维建模单元连接;
所述压力导丝测量单元,用于通过压力导丝测量冠状动脉狭窄远端的压力Pd和冠状动脉入口压力Pa
所述提取冠脉造影单元,用于选取所述测量血管的第一体位造影图像和第二体位造影图像;
所述三维建模单元,用于接收所述提取冠脉造影单元传递的第一体位造影图像和所述第二体位造影图像,三维建模获得冠状动脉三维血管模型;
所述参数测量单元,用于接收所述三维建模单元传递的冠状动脉三维血管模型,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2;以及接收所述压力导丝测量单元传递的压力Pd,根据Pd、T1、T2测量冠状动脉血管评定参数。
可选地,上述的测量冠状动脉血管评定参数的装置,所述参数测量单元包括:冠脉血流储备模块、微循环阻力系数模块,和/或冠状动脉血流储备分数模块;所述冠脉血流储备模块、所述微循环阻力系数模块均与所述三维建模单元连接;所述微循环阻力系数模块、所述冠状动脉血流储备分数模均与所述压力导丝测量单元连接;
所述冠脉血流储备模块,用于测量冠脉血流储备CFR,CFR=T1/T2
所述微循环阻力系数模块,用于测量微循环阻力系数IMR,IMR=Pd×T2
所述冠状动脉血流储备分数模块,用于测量冠状动脉血流储备分数FFR,FFR=Pd/Pa
第三方面,本申请提供了一种冠状动脉分析系统,包括:上述的测量冠状动脉血管评定参数的装置。
第四方面,本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述的测量冠状动脉血管评定参数的简化方法。
本申请实施例提供的方案带来的有益效果至少包括:
本申请提供了测量冠状动脉血管评定参数的简化方法,只在测试冠状动脉狭窄远端的压力Pd时,注射扩张药物,造影过程只需几秒钟即可停止扩张药物注射,降低扩张药物的注射时间,然后通过冠脉造影图像进行三维建模,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2;根据Pd、T1、T2测量IMR、CFR等冠状动脉血管评定参数,测量过程简单,测试结果准确,克服了全部使用压力导丝测量IMR、CFR等冠状动脉血管评定参数时的问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本申请的测量冠状动脉血管评定参数的简化方法的一个实施例的流程图;
图2为本申请的测量冠状动脉血管评定参数的简化方法的另一实施例的流程图;
图3为本申请的步骤S400的流程图;
图4为本申请的测量冠状动脉血管评定参数的装置的结构框图;
图5为本申请的参数测量单元的结构框图;
图6为本申请的三维建模单元的一个实施例的结构框图;
图7为本申请的三维建模单元的另一实施例的结构框图;
图8为本申请的图像处理模块的结构框图;
图9为参考图像;
图10为待分割的一目标图像;
图11为待分割的另一目标图像;
图12为增强后的导管图像;
图13为导管特征点的二值化图像;
图14为增强后的目标图像;
图15为冠状动脉所处位置的区域图像;
图16为结果图像;
图17为两个体位造影图像;
图18为由图17结合体位角度以及冠脉中心线生成的冠状动脉三维血管模型图;
下面对附图标记进行说明:
压力导丝测量单元110,提取冠脉造影单元120,三维建模单元130,图像读取模块131,分割模块132,血管长度测量模块133,三维建模模块134,图像处理模块135,图像去噪模块1350,导管特征点提取模块1351,冠状动脉提取模块1352,冠脉中心线提取模块136,血管直径测量模块137,参数测量单元140,冠脉血流储备模块141,微循环阻力系数模块142,冠状动脉血流储备分数模块143。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下将以图式揭露本发明的多个实施方式,为明确说明起见,许多实务上的细节将在以下叙述中一并说明。然而,应了解到,这些实务上的细节不应用以限制本发明。也就是说,在本发明的部分实施方式中,这些实务上的细节是非必要的。此外,为简化图式起见,一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。
如图1所示,本申请提供了一种测量冠状动脉血管评定参数的简化方法,包括:
S100,对测量血管进行冠脉造影;
S200,通过压力导丝测量冠状动脉狭窄远端的压力Pd
S300,选取所述测量血管的静息态下的第一体位造影图像和扩张状态下的第二体位造影图像;
S400,选取从冠脉狭窄病变近端到远端的一段血管,根据第一体位造影图像和第二体位造影图像三维建模获得冠状动脉三维血管模型;
S500,注射造影剂,根据所述造影剂在所述冠状动脉三维血管模型内的流动,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2
S600,根据Pd、T1、T2测量冠状动脉血管评定参数。
本申请的一个实施例中,S600中的冠状动脉血管评定参数包括:冠脉血流储备CFR、微循环阻力系数IMR。
本申请提供了测量冠状动脉血管评定参数的简化方法,只在测试冠状动脉狭窄远端的压力Pd时,注射扩张药物,造影过程只需几秒钟即可停止扩张药物注射,降低扩张药物的注射时间,然后通过冠脉造影图像进行三维建模,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2;根据Pd、T1、T2测量IMR、CFR等冠状动脉血管评定参数,测量过程简单,测试结果准确,克服了全部使用压力导丝测量IMR、CFR等冠状动脉血管评定参数时的问题。
需要注意的是:注射扩张药物包括:静脉或者冠脉内注射扩张药物,打入方式包括:扩张药物可以与造影剂混合注射进入静脉或者冠脉,也可以分次间隔注射均在本申请的保护范围内;只要能够起到扩张作用的药物包括腺苷、ATP等均在本申请的保护范围内。
本申请的一个实施例中,S300中CFR=T1/T2;和/或
IMR=Pd×T2
本申请的一个实施例中,上述的时间T1和T2根据心跳周期区域被分成的局部区域图像的帧数与每秒传输帧数的比值计算;即:
T=N/fps,N表示心跳周期区域被分成的局部区域图像的帧数,fps表示画面每秒播放帧数,通俗来讲就是指动画或视频的画面数,T表示某一体位造影图像中的造影剂从血管段入口到出口所经过的时间T,因此T1和T2可以根据上述公式计算得出。本申请的一个实施例中,fps=10~30;优选地,fps=15。
由于T1和T2均是基于造影图像获得的冠状动脉三维血管模型测得,因此CFR也是通过冠状动脉三维血管模型测得,无需依赖于压力导丝传感器,因此克服了压力导丝传感器在盐水的冲击下容易移动,测量不准确的问题,且测量基于造影图像时无需注射盐水,因此避免了注射盐水速度、注射量、盐水的温度对CFR测量结果的影响,提高了测量的准确率。
IMR的测量是基于压力导丝和造影图像的,由于压力导丝测量的冠状动脉狭窄远端的压力Pd更加准确,然后再基于造影图像测量时间T2,减少了扩张药物的注射时间,降低了盐水的注射次数,降低了盐水对测试结果的影响,进而提高了IMR测量结果的准确率,测试过程简单。
本申请的一个实施例中,S600中的冠状动脉血管评定参数包括:冠状动脉血流储备分数FFR,通过S200中的压力导丝测量冠状动脉入口压力Pa,FFR=Pd/Pa
本申请的一个实施例中,S300中的第一体位图像中选取的第一体位与第二体位图像中选取的第二体位的夹角大于30°。
如图2所示,本申请的一个实施例中,在S100之后,S200之前,包括:S000,向血管内注射扩张药物。
如图2所示,本申请的一个实施例中,在S100之前,包括:S110,在注射扩张药物的同时向血管内注射造影剂。需要注意的是:在造影结束后即可停止注射扩张药物,降低了扩张药物的注射量和注射次数,更加安全、可靠。
如图3所示,本申请的一个实施例中,S400包括:
S410,去除第一体位造影图像和第二体位造影图像的干扰血管,得到结果图像,具体为:
去除第一体位造影图像和第二体位造影图像的干扰血管;
对冠脉造影图像去噪,包括:静态噪声和动态噪声;
将有导管出现的第一帧分割图像定义为参考图像,将有完整冠状动脉出现的第k帧所述分割图像定义为目标图像,k为大于1的正整数;
将所述参考图像减去所述目标图像,提取所述导管的特征点O,具体方法为:将所述参考图像减去所述目标图像;去噪,包括:静态噪声和动态噪声;对所述去噪后的图像进行图像增强;对增强后的导管图像进行二值化处理,得到具有一组导管特征点O的二值化图像;
将所述目标图像减去所述参考图像,提取所述冠状动脉所处位置的区域图像,具体方法为:将所述目标图像减去所述参考图像;去噪,包括:静态噪声和动态噪声;对去噪后的所述图像进行图像增强;根据增强后的所述目标图像中各区域与所述导管特征点的位置关系,确定并提取冠状动脉的区域,即为所述冠状动脉所处位置的区域图像;
所述区域图像以所述导管的特征点作为种子点进行动态生长,获得所述结果图像,具体方法为:对所述冠状动脉所处位置的区域图像进行二值化处理,获得二值化冠状动脉图像;对所述二值化冠状动脉图像进行形态学运算,以所述导管的特征点作为种子点,所述二值化冠状动脉图像依据所述种子点所处位置进行动态区域生长,获得所述结果图像;
S420,沿着冠状动脉的延伸方向,提取每幅结果图像的冠脉中心线和直径;
S430,将每根冠脉中心线和直径均投射于三维空间上进行三维建模,获得冠状动脉三维血管模型,具体方法为:
获取每幅冠脉造影图像的体位拍摄角度;
将每根所述冠脉中心线结合体位拍摄角度投射于三维空间上,进行投影,生成冠状动脉三维血管模型。
如图4所示,本申请提供了一种测量冠状动脉血管评定参数的装置100,用于上述的测量冠状动脉血管评定参数的简化方法,包括:压力导丝测量单元110、提取冠脉造影单元120、三维建模单元130和参数测量单元140,提取冠脉造影单元120与三维建模单元130连接,参数测量单元140与压力导丝测量单元110、三维建模单元130连接;压力导丝测量单元110,用于通过压力导丝测量冠状动脉狭窄远端的压力Pd和冠状动脉入口压力Pa;提取冠脉造影单元120,用于选取测量血管的第一体位造影图像和第二体位造影图像;三维建模单元130,用于接收提取冠脉造影单元传递的第一体位造影图像和第二体位造影图像,三维建模获得冠状动脉三维血管模型;参数测量单元140,用于接收三维建模单元传递的冠状动脉三维血管模型,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2;以及接收压力导丝测量单元传递的压力Pd,根据Pd、T1、T2测量冠状动脉血管评定参数。
如图5所示,本申请的一个实施例中,参数测量单元140包括:冠脉血流储备模块141、微循环阻力系数模块142,和/或冠状动脉血流储备分数模块143,冠脉血流储备模块141、微循环阻力系数模块142均与三维建模单元130连接;微循环阻力系数模块142、冠状动脉血流储备分数模块143均与压力导丝测量单元110连接;冠脉血流储备模块141用于测量冠脉血流储备CFR,CFR=T1/T2;微循环阻力系数模块142用于测量微循环阻力系数IMR,IMR=Pd×T2;冠状动脉血流储备分数模块143用于测量冠状动脉血流储备分数FFR,FFR=Pd/Pa
参数测量单元140还包括:T1测量模块、T2测量模块和CFR测量模块,均与三维建模单元130连接,且T1测量模块、T2测量模块均与CFR测量模块连接。
如图6所示,本申请的一个实施例中,三维建模单元130包括图像读取模块131、分割模块132、血管长度测量模块133和三维建模模块134,分割模块132与图像读取模块131、血管长度测量模块133、三维建模模块134连接;图像读取模块131用于读取造影图像;分割模块132用于选取冠状动脉造影图像的一个心跳周期区域;血管长度测量模块133用于测量心跳周期区域内的血管的长度L,并将血管的长度L传递给分割模块132;三维建模模块134用于根据分割模块132选取的冠脉造影图像进行三维建模,获得冠状动脉三维血管模型。
如图7所示,本申请的一个实施例中,三维建模单元130还包括:图像处理模块135、冠脉中心线提取模块136和血管直径测量模块137,图像处理模块135与冠脉中心线提取模块136连接,三维建模模块134与冠脉中心线提取模块136、血管直径测量模块137连接。图像处理模块135用于接收分割模块132传递的至少两个体位的冠脉造影图像,并去除冠脉造影图像的干扰血管,得到如图17所示的结果图像;冠脉中心线提取模块136用于沿着冠状动脉的延伸方向,提取每幅如图17所示的结果图像的冠脉中心线;血管直径测量模块137用于测量血管直径D;三维建模模块134用于将每根冠脉中心线和直径均投射于三维空间上进行三维建模,获得冠状动脉三维血管模型。本申请实现了根据冠脉造影图像合成冠状动脉三维血管模型,弥补了行业内的空白,对于医学技术领域具有积极的作用。
本申请的一个实施例中,图像处理模块135内部设置图像去噪模块1350,用于对冠脉造影图像去噪,包括:静态噪声和动态噪声。通过去噪模块1350去除冠脉造影图像中的干扰因素,提高图像处理的质量。
如图8所示,本申请的一个实施例中,图像处理模块135内部设置均与冠脉中心线提取模块136连接的导管特征点提取模块1351和冠状动脉提取模块1352,导管特征点提取模块1351与冠状动脉提取模块1352、图像去噪模块1350连接;导管特征点提取模块1351用于将有导管出现的第一帧分割图像定义为如图9所示的参考图像,将有完整冠状动脉出现的第k帧分割图像定义为如图10和如图11所示的目标图像,k为大于1的正整数,对如图10和11所示的目标图像进行增强,得到如图12和14所示的增强后的图像;将如图9所示的参考图像减去如图10和如图11所示的目标图像,提取如图13所示的导管的特征点O;冠状动脉提取模块1352用于将如图10和如图11所示的目标图像减去如图9所示的参考图像,根据如图14所示的增强后的目标图像中各区域与导管特征点的位置关系,确定并提取冠状动脉的区域,即为如图15所示的冠状动脉所处位置的区域图像;如图15所示的区域图像以如图13所示的导管的特征点作为种子点进行动态生长,获得如图16所示的结果图像。
图像处理模块135内部还设置二值化处理模块,用于对图像进行二值化处理,以获得冠状动脉三维血管模型。
下面结合具体实施例对本申请进行具体说明:
实施例1:
如图17所示,为一位患者拍摄的两个体位的冠脉造影图像;左图为体位角度为右前斜RAO:25°和头位CRA:23°的造影图像;右图为体位角度为右前斜RAO:3°和头位CRA:30°的造影图像;
把压力导丝传感器放到患者的冠脉远端(离导引导管开口>5cm);在测量FFR时注射扩张药物,通过压力导丝测得Pa=95.2Hg,Pd=85.8mmHg;
冠状动脉三维血管模型的血管长度L值=120mm;生成的冠状动脉三维血管模型如图18所示;血管直径D值=2~4mm,T1=N1/fps1=26/15=1.7s;T2=N2/fps2=9/15=0.6s;
CFR=T1/T2=1.7/0.6=2.83;
因此IMR=85.8×0.6=51.48;
FFR=85.8/95.2=0.90。
对比例1:
与实施例1的患者相同,对比例1和实施例1均为同一位患者拍摄的同一张冠脉造影图像;
把压力导丝传感器放到患者的冠脉远端(离导引导管开口>5cm),通过导管向血管中注入3ml生理盐水,如果检测到血液温度回复到正常值,则再次通过导管向血管中注入3ml生理盐水,重复上述过程3次,然后记录T1,T1为1.6s;向血管通入扩张药物,使血管达到并保持扩张状态(保证通入扩张药物前后的压力导丝传感器处于相同位置),通过导管向血管中注入3ml生理盐水,如果检测到血液温度回复到正常值,则再次通过导管向血管中注入3ml生理盐水,重复上述过程3次,然后记录T2,T2为0.58s,测得冠脉远端的压力Pd=88mmHg,冠脉入口压力Pa=94.8mmHg;
CFR=1.6/0.58=2.76;
IMR=Pd×T2=88×0.58=51.04;
FFR=88/94.8=0.928。
通过实施例1和对比例1的比较,差值在0.5以内,可知IMR测量结果基本相同,因此实施例1的测量结果准确,且本申请的实施例采用压力导丝测量远端的压力,但是不使用生理盐水,T1、T2通过三维血管模型测得;且通过造影图像实现了IMR的测量,弥补了行业内的空白,操作更加简单,还实现了FFR的测量,且无需采用生理盐水,解决了压力导丝传感器在生理盐水的冲力作用下,导致的压力导丝的位置很难控制的问题,解决了远端压力测量不准的问题。
本申请提供了一种冠状动脉分析系统,包括:上述的测量冠状动脉血管评定参数的装置100。
本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述的测量冠状动脉血管评定参数的简化方法。
所属技术领域的技术人员知道,本发明的各个方面可以实现为系统、方法或计算机程序产品。因此,本发明的各个方面可以具体实现为以下形式,即:完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等),或硬件和软件方面结合的实施方式,这里可以统称为“电路”、“模块”或“系统”。此外,在一些实施例中,本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式,该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或系统的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。
例如,可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件,可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作系统执行的多个软件指令。在本发明的示例性实施例中,由数据处理器来执行如本文的根据方法和/或系统的示例性实施例的一个或多个任务,诸如用于执行多个指令的计算平台。可选地,该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器,例如,磁硬盘和/或可移动介质。可选地,也提供了一种网络连接。可选地也提供显示器和/或用户输入设备,诸如键盘或鼠标。
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项:
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括(但不限于)无线、有线、光缆、RF等等,或者上述的任意合适的组合。
例如,可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码,包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言,诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。
也可以把这些计算机程序指令存储在计算机可读介质中,这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作,从而,存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。
还可将计算机程序指令加载到计算机(例如,冠状动脉分析系统)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程,使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。
本发明的以上的具体实例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (12)

1.一种测量冠状动脉血管评定参数的简化方法,其特征在于,包括:
对测量血管进行冠脉造影;
通过压力导丝测量冠状动脉狭窄远端的压力Pd
选取所述测量血管的静息态下的第一体位造影图像和扩张状态下的第二体位造影图像;
选取从冠脉狭窄病变近端到远端的一段血管,根据所述第一体位造影图像和所述第二体位造影图像三维建模获得冠状动脉三维血管模型;
注射造影剂,根据所述造影剂在所述冠状动脉三维血管模型内的流动,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2
根据Pd、T1、T2测量冠状动脉血管评定参数。
2.根据权利要求1所述的测量冠状动脉血管评定参数的简化方法,其特征在于,所述冠状动脉血管评定参数包括:冠脉血流储备CFR、微循环阻力系数IMR。
3.根据权利要求2所述的测量冠状动脉血管评定参数的简化方法,其特征在于,所述CFR=T1/T2;和/或
所述IMR=Pd×T2
4.根据权利要求1~3任一项所述的测量冠状动脉血管评定参数的简化方法,其特征在于,所述第一体位与所述第二体位的夹角大于30°。
5.根据权利要求1~3任一项所述的测量冠状动脉血管评定参数的简化方法,其特征在于,在所述通过压力导丝测量冠状动脉狭窄远端的压力Pd之前,包括:向血管内注射扩张药物。
6.根据权利要求1~3任一项所述的测量冠状动脉血管评定参数的简化方法,其特征在于,在所述通过压力导丝测量冠状动脉狭窄远端的压力Pd之后,所述对测量血管进行冠脉造影之前,包括:在注射扩张药物的同时向所述血管内注射造影剂。
7.根据权利要求1~3任一项所述的测量冠状动脉血管评定参数的简化方法,其特征在于,所述根据所述第一体位造影图像和所述第二体位造影图像三维建模获得冠状动脉三维血管模型包括:
去除第一体位造影图像和所述第二体位造影图像的干扰血管,得到结果图像;
沿着所述冠状动脉的延伸方向,提取每幅所述结果图像的冠脉中心线和直径;
将每根所述冠脉中心线和直径均投射于三维空间上进行三维建模,获得冠状动脉三维血管模型。
8.根据权利要求1~3任一项所述的测量冠状动脉血管评定参数的简化方法,其特征在于,所述时间T1和T2根据心跳周期区域被分成的局部区域图像的帧数与每秒传输帧数的比值计算。
9.一种测量冠状动脉血管评定参数的装置,用于权利要求1~8任一项所述的测量冠状动脉血管评定参数的简化方法,其特征在于,包括:压力导丝测量单元、提取冠脉造影单元、三维建模单元和参数测量单元,所述提取冠脉造影单元与三维建模单元连接,所述参数测量单元与所述压力导丝测量单元、所述三维建模单元连接;
所述压力导丝测量单元,用于通过压力导丝测量冠状动脉狭窄远端的压力Pd和冠状动脉入口压力Pa
所述提取冠脉造影单元,用于选取所述测量血管的第一体位造影图像和第二体位造影图像;
所述三维建模单元,用于接收所述提取冠脉造影单元传递的第一体位造影图像和所述第二体位造影图像,三维建模获得冠状动脉三维血管模型;
所述参数测量单元,用于接收所述三维建模单元传递的冠状动脉三维血管模型,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T1,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T2;以及接收所述压力导丝测量单元传递的压力Pd,根据Pd、T1、T2测量冠状动脉血管评定参数。
10.根据权利要求9所述的测量冠状动脉血管评定参数的装置,其特征在于,所述参数测量单元包括:冠脉血流储备模块、微循环阻力系数模块,和/或冠状动脉血流储备分数模块;所述冠脉血流储备模块、所述微循环阻力系数模块均与所述三维建模单元连接;所述微循环阻力系数模块、所述冠状动脉血流储备分数模均与所述压力导丝测量单元连接;
所述冠脉血流储备模块,用于测量冠脉血流储备CFR,CFR=T1/T2
所述微循环阻力系数模块,用于测量微循环阻力系数IMR,IMR=Pd×T2
所述冠状动脉血流储备分数模块,用于测量冠状动脉血流储备分数FFR,FFR=Pd/Pa
11.一种冠状动脉分析系统,其特征在于,包括:权利要求9或10所述的测量冠状动脉血管评定参数的装置。
12.一种计算机存储介质,其特征在于,计算机程序被处理器执行时实现权利要求1~8任一项所述的测量冠状动脉血管评定参数的简化方法。
CN201910834871.6A 2018-09-19 2019-09-05 测量冠状动脉血管评定参数的简化方法、装置及系统 Pending CN110522439A (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2022513657A JP7457105B2 (ja) 2018-09-19 2019-11-01 冠状動脈血管評定パラメータを測定する装置及びシステム
PCT/CN2019/115027 WO2021042477A1 (zh) 2018-09-19 2019-11-01 测量冠状动脉血管评定参数的简化方法、装置及系统
EP19943973.8A EP4026491A4 (en) 2018-09-19 2019-11-01 SIMPLIFIED METHOD, APPARATUS AND SYSTEM FOR MEASURING CORONARY ARTERY VESSEL ASSESSMENT PARAMETERS
US17/683,490 US20220183643A1 (en) 2018-09-19 2022-03-01 Simplified method, apparatus and system for measuring coronary artery vascular evaluation parameters

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN201811093192.XA CN109363651A (zh) 2018-09-19 2018-09-19 一种用于测量血流储备分数时获取主动脉压的装置
CN201811093192X 2018-09-19
CN2019102065412 2019-03-19
CN201910206541.2A CN109770888A (zh) 2019-03-19 2019-03-19 基于压力传感器和造影图像计算瞬时无波形比率的方法

Publications (1)

Publication Number Publication Date
CN110522439A true CN110522439A (zh) 2019-12-03

Family

ID=68261494

Family Applications (5)

Application Number Title Priority Date Filing Date
CN201910704327.XA Pending CN110384493A (zh) 2018-09-19 2019-07-31 测量微循环阻力指数的系统以及冠脉分析系统
CN201910704330.1A Pending CN110384494A (zh) 2018-09-19 2019-07-31 测量微循环阻力指数的方法
CN201910834871.6A Pending CN110522439A (zh) 2018-09-19 2019-09-05 测量冠状动脉血管评定参数的简化方法、装置及系统
CN201910835038.3A Active CN110367965B (zh) 2018-09-19 2019-09-05 便捷测量冠状动脉血管评定参数的方法、装置及系统
CN201910846497.1A Pending CN110393516A (zh) 2018-09-19 2019-09-09 基于影像和压力传感器计算微循环指标的方法装置及系统

Family Applications Before (2)

Application Number Title Priority Date Filing Date
CN201910704327.XA Pending CN110384493A (zh) 2018-09-19 2019-07-31 测量微循环阻力指数的系统以及冠脉分析系统
CN201910704330.1A Pending CN110384494A (zh) 2018-09-19 2019-07-31 测量微循环阻力指数的方法

Family Applications After (2)

Application Number Title Priority Date Filing Date
CN201910835038.3A Active CN110367965B (zh) 2018-09-19 2019-09-05 便捷测量冠状动脉血管评定参数的方法、装置及系统
CN201910846497.1A Pending CN110393516A (zh) 2018-09-19 2019-09-09 基于影像和压力传感器计算微循环指标的方法装置及系统

Country Status (5)

Country Link
US (3) US11779294B2 (zh)
EP (3) EP4026492A4 (zh)
JP (2) JP7093584B2 (zh)
CN (5) CN110384493A (zh)
WO (5) WO2020057323A2 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111067494A (zh) * 2019-12-27 2020-04-28 西北工业大学 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法
CN111067495A (zh) * 2019-12-27 2020-04-28 西北工业大学 基于血流储备分数和造影图像的微循环阻力计算方法
CN111340794A (zh) * 2020-03-09 2020-06-26 中山大学 冠状动脉狭窄的量化方法及装置
CN111449750A (zh) * 2020-04-02 2020-07-28 刘强 一种心内科防止大量出血的介入控制系统及控制方法
CN112686991A (zh) * 2021-01-08 2021-04-20 博动医学影像科技(上海)有限公司 一种采用混合方式重建血管正常管腔形态的方法及系统

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110384493A (zh) * 2018-09-19 2019-10-29 苏州润迈德医疗科技有限公司 测量微循环阻力指数的系统以及冠脉分析系统
CN112155580B (zh) * 2019-11-20 2021-07-13 苏州润迈德医疗科技有限公司 基于造影图像修正血流速度和微循环参数的方法及装置
CN110786841B (zh) * 2019-11-04 2021-05-25 苏州润迈德医疗科技有限公司 基于微循环阻力指数调节最大充血状态流速的方法及装置
CN110786840B (zh) * 2019-11-04 2021-06-08 苏州润迈德医疗科技有限公司 基于生理参数获取血管评定参数的方法、装置及存储介质
CN110786842B (zh) * 2019-11-04 2022-10-04 苏州润迈德医疗科技有限公司 测量舒张期血流速度的方法、装置、系统及存储介质
CN110929604B (zh) * 2019-11-11 2023-08-18 苏州润迈德医疗科技有限公司 基于造影图像的流速的筛选方法、装置、系统和存储介质
CN110889896B (zh) * 2019-11-11 2024-03-22 苏州润迈德医疗科技有限公司 获取血管狭窄病变区间及三维合成方法、装置和系统
CN112132882A (zh) * 2019-11-19 2020-12-25 苏州润迈德医疗科技有限公司 从冠状动脉二维造影图像中提取血管中心线的方法和装置
CN112151180B (zh) * 2019-12-05 2024-03-08 苏州润迈德医疗科技有限公司 具有狭窄病变的血管数学模型的合成方法和装置
CN112116711B (zh) * 2019-12-05 2024-01-23 苏州润迈德医疗科技有限公司 用于流体力学分析的圆台血管数学模型的合成方法和装置
CN212346515U (zh) * 2019-12-31 2021-01-15 深圳北芯生命科技有限公司 血管内压力测量系统
US11950893B2 (en) * 2020-01-17 2024-04-09 The Research Foundation For The State University Of New York Apparatus and method of determining dynamic vascular parameters of blood flow
CN111312374B (zh) * 2020-01-21 2024-03-22 上海联影智能医疗科技有限公司 医学图像处理方法、装置、存储介质及计算机设备
CN112089433B (zh) * 2020-09-24 2022-04-08 杭州阿特瑞科技有限公司 一种基于cta和dsa的冠状动脉血流储备分数测量方法
CN112690814B (zh) * 2020-11-06 2022-10-14 杭州阿特瑞科技有限公司 一种低误差的冠状动脉血流储备分数测量方法
CN112419280A (zh) * 2020-11-25 2021-02-26 苏州润迈德医疗科技有限公司 精确获取狭窄病变区间的方法及存储介质
CN112419484B (zh) * 2020-11-25 2024-03-22 苏州润迈德医疗科技有限公司 三维血管合成方法、系统及冠状动脉分析系统和存储介质
CN112971818B (zh) * 2021-01-28 2022-10-04 杭州脉流科技有限公司 微循环阻力指数的获取方法、装置、计算机设备和存储介质
CN113180614B (zh) * 2021-06-02 2023-08-04 北京阅影科技有限公司 无导丝ffr、无导丝imr和无导丝cfr的检测方法
CN112890843A (zh) * 2021-03-24 2021-06-04 苏州润迈德医疗科技有限公司 一种冠状动脉分析系统及操作方法
CN113706559A (zh) * 2021-09-13 2021-11-26 复旦大学附属中山医院 基于医学图像的血管分段提取方法和装置
CN115690309B (zh) * 2022-09-29 2023-07-18 中国人民解放军总医院第一医学中心 一种冠脉cta自动三维后处理方法和装置
CN116035549A (zh) * 2022-11-23 2023-05-02 杭州脉流科技有限公司 获得冠脉血流储备分数的方法和计算机设备
CN115869003A (zh) * 2022-12-30 2023-03-31 杭州脉流科技有限公司 基于ct图像的冠状动脉微循环阻力指数计算方法和装置
CN116664564B (zh) * 2023-07-28 2023-10-31 杭州脉流科技有限公司 基于颅内医学影像获取血流量的方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140181717A1 (en) * 2012-12-26 2014-06-26 Volcano Corporation Measurement Navigation in a Multi-Modality Medical Imaging System
CN104582572A (zh) * 2012-08-16 2015-04-29 株式会社东芝 图像处理装置、医用图像诊断装置以及血压监视器
CN106214140A (zh) * 2016-10-17 2016-12-14 苏州亘科医疗科技有限公司 血管内血压测量导丝
CN108245178A (zh) * 2018-01-11 2018-07-06 苏州润迈德医疗科技有限公司 一种基于x射线冠脉造影图像的血液流动速度计算方法
CN108550189A (zh) * 2018-05-03 2018-09-18 苏州润迈德医疗科技有限公司 基于造影图像和流体力学模型的微循环阻力指数计算方法

Family Cites Families (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6672172B2 (en) * 2000-01-31 2004-01-06 Radi Medical Systems Ab Triggered flow measurement
US8315812B2 (en) * 2010-08-12 2012-11-20 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9119540B2 (en) * 2010-09-16 2015-09-01 Siemens Aktiengesellschaft Method and system for non-invasive assessment of coronary artery disease
JP5309109B2 (ja) * 2010-10-18 2013-10-09 富士フイルム株式会社 医用画像処理装置および方法、並びにプログラム
CN108294735B (zh) * 2012-03-13 2021-09-07 西门子公司 用于冠状动脉狭窄的非侵入性功能评估的方法和系统
US10373700B2 (en) * 2012-03-13 2019-08-06 Siemens Healthcare Gmbh Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance
EP2825091B1 (en) * 2012-03-15 2016-08-24 Siemens Healthcare GmbH A framework for personalization of coronary flow computations during rest and hyperemia
JP6134789B2 (ja) * 2012-06-26 2017-05-24 シンク−アールエックス,リミティド 管腔器官における流れに関連する画像処理
US9247918B2 (en) * 2012-07-09 2016-02-02 Siemens Aktiengesellschaft Computation of hemodynamic quantities from angiographic data
EP3903672B1 (en) * 2012-08-03 2023-11-01 Philips Image Guided Therapy Corporation Devices, systems, and methods for assessing a vessel
US9814433B2 (en) * 2012-10-24 2017-11-14 Cathworks Ltd. Creating a vascular tree model
EP2934307B1 (en) * 2012-12-21 2020-08-05 Volcano Corporation Functional gain measurement technique and representation
EP2956050A4 (en) * 2013-02-18 2016-10-26 Univ Ramot FROM INTRAVASCULAR PRESSURE LOSS DERIVED ARTERY STIFFNESS AND REDUCTION OF LIGHT PRESSURE IMPACT
US9700219B2 (en) * 2013-10-17 2017-07-11 Siemens Healthcare Gmbh Method and system for machine learning based assessment of fractional flow reserve
JP6262027B2 (ja) * 2014-03-10 2018-01-17 東芝メディカルシステムズ株式会社 医用画像処理装置
US9087147B1 (en) * 2014-03-31 2015-07-21 Heartflow, Inc. Systems and methods for determining blood flow characteristics using flow ratio
US10307070B2 (en) * 2014-04-04 2019-06-04 St. Jude Medical Coordination Center Bvba Intravascular pressure and flow data diagnostic systems, devices, and methods
US10803995B2 (en) * 2014-05-05 2020-10-13 Siemens Healthcare Gmbh Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in diseased and hypothetical normal anatomical models
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
US9754082B2 (en) * 2014-05-30 2017-09-05 Heartflow, Inc. Systems and methods for reporting blood flow characteristics
WO2016001017A1 (en) * 2014-06-30 2016-01-07 Koninklijke Philips N.V. Apparatus for determining a fractional flow reserve value
WO2016007656A1 (en) * 2014-07-08 2016-01-14 Scinovia, Corp. Fluid flow rate determinations using velocity vector maps
US9349178B1 (en) * 2014-11-24 2016-05-24 Siemens Aktiengesellschaft Synthetic data-driven hemodynamic determination in medical imaging
EP3087921A1 (en) * 2015-04-27 2016-11-02 Coronary Technologies SARL Computer-implemented method for identifying zones of stasis and stenosis in blood vessels
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
CN105326486B (zh) * 2015-12-08 2017-08-25 博动医学影像科技(上海)有限公司 血管压力差与血流储备分数的计算方法及系统
CN105559810B (zh) * 2015-12-10 2017-08-08 博动医学影像科技(上海)有限公司 血管单位时间血流量与血流速度的计算方法
US10758125B2 (en) * 2016-04-29 2020-09-01 Siemens Healthcare Gmbh Enhanced personalized evaluation of coronary artery disease using an integration of multiple medical imaging techniques
US10674986B2 (en) * 2016-05-13 2020-06-09 General Electric Company Methods for personalizing blood flow models
CN106023202B (zh) * 2016-05-20 2017-06-09 苏州润心医疗科技有限公司 基于心脏ct图像的冠状动脉血流储备分数计算方法
CN106073894B (zh) * 2016-05-31 2017-08-08 博动医学影像科技(上海)有限公司 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和系统
CN106327487B (zh) * 2016-08-18 2018-01-02 苏州润迈德医疗科技有限公司 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法
CN106264514B (zh) * 2016-09-27 2023-05-05 上海爱声生物医疗科技有限公司 一种无线血流储备分数测量系统
CN106419890B (zh) * 2016-11-14 2024-04-30 佛山科学技术学院 一种基于时空调制的血流速度测量装置及方法
CN106650267B (zh) * 2016-12-28 2020-03-17 北京昆仑医云科技有限公司 计算血流储备分数的系统和设置边界条件的方法
US20180192916A1 (en) * 2017-01-10 2018-07-12 General Electric Company Imaging system for diagnosing patient condition
EP3375364A4 (en) * 2017-01-23 2019-01-23 Shanghai United Imaging Healthcare Co., Ltd. SYSTEM AND METHOD FOR ANALYZING THE STATUS OF BLOOD CIRCULATION
US11219741B2 (en) * 2017-08-09 2022-01-11 Medtronic Vascular, Inc. Collapsible catheter and method for calculating fractional flow reserve
CN107730540B (zh) * 2017-10-09 2020-11-17 全景恒升(北京)科学技术有限公司 基于高精度匹配模型的冠脉参数的计算方法
CN107978371B (zh) * 2017-11-30 2021-04-02 博动医学影像科技(上海)有限公司 快速计算微循环阻力的方法及系统
US10580526B2 (en) * 2018-01-12 2020-03-03 Shenzhen Keya Medical Technology Corporation System and method for calculating vessel flow parameters based on angiography
CN108186038B (zh) * 2018-02-11 2020-11-17 杭州脉流科技有限公司 基于动脉造影影像计算冠脉血流储备分数的系统
CN108992057B (zh) * 2018-06-05 2021-08-10 杭州晟视科技有限公司 一种确定冠状动脉血流储备分数ffr的方法和装置
CN108742587B (zh) * 2018-06-20 2021-04-27 博动医学影像科技(上海)有限公司 基于病史信息获取血流特征值的方法及装置
CN110384493A (zh) * 2018-09-19 2019-10-29 苏州润迈德医疗科技有限公司 测量微循环阻力指数的系统以及冠脉分析系统
CN109770888A (zh) * 2019-03-19 2019-05-21 苏州润迈德医疗科技有限公司 基于压力传感器和造影图像计算瞬时无波形比率的方法
CN109523515A (zh) * 2018-10-18 2019-03-26 深圳市孙逸仙心血管医院(深圳市心血管病研究所) 血流储备分数的计算方法
CN109770930B (zh) * 2019-01-29 2021-03-09 浙江大学 一种冠状动脉微循环阻力的确定方法和装置
CN109805949B (zh) * 2019-03-19 2020-05-22 苏州润迈德医疗科技有限公司 基于压力传感器和造影图像计算血流储备分数的方法
CN110120031B (zh) * 2019-04-02 2022-05-31 四川锦江电子科技有限公司 一种得到血管血流储备分数的方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104582572A (zh) * 2012-08-16 2015-04-29 株式会社东芝 图像处理装置、医用图像诊断装置以及血压监视器
US20140181717A1 (en) * 2012-12-26 2014-06-26 Volcano Corporation Measurement Navigation in a Multi-Modality Medical Imaging System
CN106214140A (zh) * 2016-10-17 2016-12-14 苏州亘科医疗科技有限公司 血管内血压测量导丝
CN108245178A (zh) * 2018-01-11 2018-07-06 苏州润迈德医疗科技有限公司 一种基于x射线冠脉造影图像的血液流动速度计算方法
CN108550189A (zh) * 2018-05-03 2018-09-18 苏州润迈德医疗科技有限公司 基于造影图像和流体力学模型的微循环阻力指数计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郭雨龙,郭丽君: "心肌微循环及微循环阻力指数研究近况", 《中国介入心脏病学杂志》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111067494A (zh) * 2019-12-27 2020-04-28 西北工业大学 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法
CN111067495A (zh) * 2019-12-27 2020-04-28 西北工业大学 基于血流储备分数和造影图像的微循环阻力计算方法
CN111067494B (zh) * 2019-12-27 2022-04-26 西北工业大学 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法
CN111340794A (zh) * 2020-03-09 2020-06-26 中山大学 冠状动脉狭窄的量化方法及装置
CN111340794B (zh) * 2020-03-09 2023-07-04 中山大学 冠状动脉狭窄的量化方法及装置
CN111449750A (zh) * 2020-04-02 2020-07-28 刘强 一种心内科防止大量出血的介入控制系统及控制方法
CN112686991A (zh) * 2021-01-08 2021-04-20 博动医学影像科技(上海)有限公司 一种采用混合方式重建血管正常管腔形态的方法及系统

Also Published As

Publication number Publication date
US20220183644A1 (en) 2022-06-16
EP4026491A4 (en) 2023-09-06
WO2021042479A1 (zh) 2021-03-11
EP4026491A1 (en) 2022-07-13
EP4029438A1 (en) 2022-07-20
CN110384494A (zh) 2019-10-29
CN110367965A (zh) 2019-10-25
CN110367965B (zh) 2022-03-08
CN110393516A (zh) 2019-11-01
US20210145390A1 (en) 2021-05-20
EP4026492A4 (en) 2023-10-11
CN110384493A (zh) 2019-10-29
US11779294B2 (en) 2023-10-10
WO2020057323A2 (zh) 2020-03-26
JP2022542189A (ja) 2022-09-29
EP4029438A4 (en) 2023-10-04
WO2021042477A1 (zh) 2021-03-11
WO2020057324A1 (zh) 2020-03-26
US20220183643A1 (en) 2022-06-16
JP7093584B2 (ja) 2022-06-30
EP4026492A1 (en) 2022-07-13
JP2022511211A (ja) 2022-01-31
JP7457105B2 (ja) 2024-03-27
WO2021046990A1 (zh) 2021-03-18

Similar Documents

Publication Publication Date Title
CN110522439A (zh) 测量冠状动脉血管评定参数的简化方法、装置及系统
US10872698B2 (en) Method and system for enhancing medical image-based blood flow computations using physiological measurements
CN107978371B (zh) 快速计算微循环阻力的方法及系统
US10354744B2 (en) Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance
CN105380598B (zh) 用于针对动脉狭窄的自动治疗规划的方法和系统
Boettler et al. New aspects of the ventricular septum and its function: an echocardiographic study
JP2022547103A (ja) 冠状動脈血管評定パラメータを簡単に測定する方法、装置及びシステム
CN110085321A (zh) 对冠脉循环进行多尺度的解剖学和功能建模的方法和系统
CN108992059A (zh) 用于根据血管几何形状和生理学估计血流特性的系统和方法
EP3937186A1 (en) System and method for vascular tree generation using joint prior information
US20210236000A1 (en) Method, device and system for acquiring blood vessel evaluation parameters based on angiographic image
CN109065170A (zh) 获取血管压力差的方法及装置
WO2021087961A1 (zh) 测量舒张期血流速度的方法、装置、系统及存储介质
CN108717874A (zh) 基于特定的生理参数获取血管压力值的方法及装置
CN110432886A (zh) 获取冠脉出口处的平均血流量、流速的方法、装置和系统
Shi et al. Patient-specific non-invasive estimation of pressure gradient across aortic coarctation using magnetic resonance imaging
US20050085707A1 (en) Method and arrangement for determining indicator dilution curves of an indicator in a bloodstream and cardiac parameters
WO2020083390A1 (zh) 获取心表大动脉的血流量的方法、装置、系统及存储介质
WO2022161239A1 (zh) 获取流量损失模型、损失比、供血能力的方法和系统
CN114664455A (zh) 一种冠状动脉血流储备分数计算方法及装置
CN112690814A (zh) 一种低误差的冠状动脉血流储备分数测量方法
Messina et al. Transpulmonary thermodilution
Zacharoulis et al. Measurement of stroke volume from pulmonary artery pressure record in man.
Ozkan et al. Intracardiac Shunts
Gordon et al. Modeling the dynamics of lung tissue with pulsating blood-flow in pulmonary arteries for bioimpedance simulation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20191203