CN105559810A - 血管单位时间血流量与血流速度的计算方法 - Google Patents
血管单位时间血流量与血流速度的计算方法 Download PDFInfo
- Publication number
- CN105559810A CN105559810A CN201510916119.8A CN201510916119A CN105559810A CN 105559810 A CN105559810 A CN 105559810A CN 201510916119 A CN201510916119 A CN 201510916119A CN 105559810 A CN105559810 A CN 105559810A
- Authority
- CN
- China
- Prior art keywords
- blood flow
- curve
- value
- gray scale
- area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000017531 blood circulation Effects 0.000 title claims abstract description 146
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 43
- 238000004364 calculation method Methods 0.000 title claims abstract description 16
- 230000000747 cardiac effect Effects 0.000 claims abstract description 41
- 238000000034 method Methods 0.000 claims description 51
- 239000002872 contrast media Substances 0.000 claims description 31
- 238000002601 radiography Methods 0.000 claims description 17
- 230000002792 vascular Effects 0.000 claims description 17
- 238000002583 angiography Methods 0.000 claims description 14
- 230000000630 rising effect Effects 0.000 claims description 11
- 230000007423 decrease Effects 0.000 claims description 10
- 238000000205 computational method Methods 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 6
- 210000005036 nerve Anatomy 0.000 claims description 6
- 210000002254 renal artery Anatomy 0.000 claims description 6
- 230000002889 sympathetic effect Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 5
- 239000000155 melt Substances 0.000 claims description 5
- 210000001367 artery Anatomy 0.000 claims description 4
- 206010028980 Neoplasm Diseases 0.000 claims description 3
- 230000036770 blood supply Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000010247 heart contraction Effects 0.000 claims description 3
- 230000001225 therapeutic effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 claims description 2
- 238000002844 melting Methods 0.000 claims description 2
- 230000008018 melting Effects 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims description 2
- 238000002679 ablation Methods 0.000 claims 1
- 230000001174 ascending effect Effects 0.000 abstract description 3
- 210000004369 blood Anatomy 0.000 description 11
- 239000008280 blood Substances 0.000 description 11
- 238000002586 coronary angiography Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 4
- PVBALTLWZVEAIO-UHFFFAOYSA-N diodone Chemical compound OC(=O)CN1C=C(I)C(=O)C(I)=C1 PVBALTLWZVEAIO-UHFFFAOYSA-N 0.000 description 4
- 229960001845 diodone Drugs 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 210000004165 myocardium Anatomy 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 239000012895 dilution Substances 0.000 description 3
- 238000010790 dilution Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 210000004351 coronary vessel Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 230000010412 perfusion Effects 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 201000000057 Coronary Stenosis Diseases 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003113 dilution method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 210000005259 peripheral blood Anatomy 0.000 description 1
- 239000011886 peripheral blood Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000002966 stenotic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 210000005166 vasculature Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- 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
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
-
- 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/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/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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Vascular Medicine (AREA)
- Human Computer Interaction (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
Abstract
本发明公开了一种血管单位时间血流量与血流速度的计算方法，包括：选定感兴趣区域并计算与拟合出在一段影像序列中该感兴趣区域灰度随时间变化的拟合曲线，计算灰度拟合曲线下降(或者上升)最快的时间点，在以该时间点为中心的预定时间区间内对灰度变化拟合曲线进行积分得到面积值；获取该面积值对应的单位时间血流量，结合血管管腔面积，可进一步获得血流速度。本发明的创新之处在于，以灰度拟合曲线下降(或者上升)最快位置为中心，计算整一个心动周期时间内的血流量，从而更准确地计算出单位时间血流量与血流速度，有效地避免了不恰当时间区间进行计算所导致的误差。
Description
技术领域
本发明涉及应用于医疗领域，尤其涉及应用在基于影像的单位时间血流量与血流速度的准确、快速、无创计算。
背景技术
斑块引起的血管狭窄影响到心肌血流供应，对人的健康造成威胁。通过冠脉造影，可以显示冠脉狭窄的严重程度，但是不能反映狭窄的功能性意义。血流储备分数(FFR)被评价为诊断冠脉功能学的金标准，其定义为病变冠脉能为心肌提供的最大血流量与该冠脉完全正常时最大供血流量之比，可以通过在冠脉最大充血状态下，狭窄病变远端的压力与狭窄近端压力的比值计算得出。
通过压力传感器对血管进行有创侵入性压力测量不仅工作量大，而且存在着损伤血管的风险。通过三维或者二维定量冠脉造影可以获得冠脉系统的几何模型从而进行计算机流体力学分析，然而当前速度获得方法是通过半自动(人工+计算机)方法获得，重建过程需要精确勾勒出血管边界，历时较久且操作者要具备丰富经验，同时需要人工确定造影剂进入冠状动脉的起始位置(起始帧)以及充盈到目标远端的结束位置(结束帧)，另外，造影剂在微循环充分扩张时的充盈过程较短，通过感兴趣的血管段通常仅占用了一个心跳周期的某个时间段，这样会导致在心跳周期的不同时段所确定的血流速度的差异性非常大。除此之外，通常情况下只对一根狭窄血管做重建，忽略了边支对血流速度的影响，导致该方法最终计算出的血流速度与实际值之间的误差较大，影响了后续的血流储备分数(FFR)求解结果的准确性。
现有技术中，比较典型的计算血流速度的方法有如下几类：
首次分布分析方法：参见对比文件1(WongJT，DucoteJL，TongX，HassaneinMT，SabeeM.AutomatedTechniqueforAngiographicDeterminationofCoronaryBloodFlowandLumenVolume1[J].《AcademicRadiology》，2006，13(2)：186-194)，Wong等人提出了首次分布分析的方法计算血流速度。该方法首先通过碘造影剂校正平台，对半径、长度已知的圆管填充满密度已知的碘造影剂后进行X射线造影，从圆管的造影结果中拟合出整合灰度值随碘质量变化的线性曲线。利用冠脉造影，提取冠脉感兴趣区域的整合灰度值随时间变化曲线，根据拟合的灰度值-碘质量线性曲线，得到感兴趣区域内碘造影剂质量随时间变化曲线，在已知碘造影剂密度的条件下，可以得到感兴趣区域血管内液体体积随时间变化曲线，从体积时间变化曲线中求得血流量和血流速度。该技术的缺陷在于所述的感兴趣区域不会随着造影序列目标血管位置移动而相应移位；且对于平均血流速度计算的时间段选取没有规定，可能是非整数个心动周期。导致该技术方案所求得的血流量和血流速度的精确度较低、数据的稳定性较差。
TIMI记帧法：参见对比文件2(TuS，BarbatoE，KoszegiZetal.Fractionalflowreservecalculationfrom3-dimensionalquantitativecoronaryangiographyandTIMIframecount：afastcomputermodeltoquantifythefunctionalsignificanceofmoderatelyobstructedcoronaryafteries[J].JACCCardiovascularinterventions，2014，(7)：768-777)，该文件2公开了如下方案：播放冠脉造影图像，观察造影剂进入靶血管近端解剖标志点到远端解剖标志点所经历的帧数，根据造影图像的帧频，得到造影剂从近端标志点到远端标志点传输时间。再利用冠脉造影三维重建得到标志点之间的距离，得到造影剂的传输时间，用距离除以传输时间近似为血流速度。
数字跟踪法：参见对比文件3(陈兴新，骆秉锉，杨瑞华，陈莉莉.冠脉造影数字跟踪技术测量冠脉血流速度的临床研究[J].生物医学工程杂志，2007，24(2)：294-298)，该文件3公开如下技术方案：电影环放冠脉造影序列图像，选择靶血管，采用数字跟踪技术测定前后两帧或N帧(标志点)图像移动(径线或曲线)的距离和经过的时间，自动计算血流速度。
Doppler导丝法：参见对比文件1，该文件1还公开了：将Doppler导丝插入到冠脉血管中，利用超声波振源与相对运动的血液间就产生的多普勒效应测量血流速度。
温度稀释法：将压力导丝插入到冠脉血管中，向冠脉中注射室温的生理盐水，其随血液的流动而被稀释，并吸收血液的热量而升高温度。这个温度稀释过程由导管前段的热敏电阻检测，经检测仪记录可得到温度时间稀释曲线。根据热稀释理论，血流速度与指示剂的平均转运时间成反比，根据温度时间稀释曲线可计算血流速度。
上述专利文献尽管从不同角度、不同计算方法中给出了计算单位时间血流量或血流速度的方法，但其仍具有至少以下一种或多个技术缺陷：
(1)血流速度的测量准确度受到血管长度的估计、冠脉造影的短缩重叠现象而限制。
(2)需要人工确定造影剂进入冠状动脉的起始位置(起始帧)以及充盈到目标远端的结束位置(结束帧)，计算准确度受人为因素影响。
(3)血流速度的计算要在造影剂离开感兴趣血管段之前测得，造影剂的充盈速度快，通过感兴趣血管段的时间仅占用了一个心跳周期的某个时间段(通常少于一个心动周期)，这样会导致在心跳周期的不同时段所确定的血流速度通常是非整数心动周期的血流速度均值，与实际平均血流速度存在较大的偏差。
(4)选择大的感兴趣区域，如整个冠状动脉及其下游所灌注的心肌的微血管，进行计算，无法确保造影剂充盈感兴趣区域的过程所对应的图像灰度变化曲线下降或上升时间段为整数个心动周期，这样会导致计算得到的单位时间血流量或血流速度均值与实际单位时间血流量或平均血流速度存在较大的偏差。
(5)直接测量血流速度具有有创性，而且Doppler血流导丝仪器昂贵，病人难以承受，限制了该方法的临床推广应用。
(6)造影感兴趣区域不会随着造影序列中目标血管位置移动而相应移位，导致感兴趣区域内的灰度值变化不仅仅由于造影剂在目标血管中传输导致，还可能由于目标血管移出和边支血管移入导致，从而导致计算的单位时间血流量或血流速度不准确。
发明内容
本发明所要解决的技术问题是提供一种新的血管单位时间血流量与血流速度计算方法，具体方案包括：
一种血管单位时间血流量的计算方法，该方法包括：确定血管的感兴趣区域；计算并拟合出所述感兴趣区域内的灰度拟合曲线；确定一预定时间区间内的最大灰度值曲线或最小灰度值曲线；计算预定时间区间内最大灰度值曲线或最小灰度值曲线与灰度拟合曲线所包围的区域面积值；基于该区域面积值，获取该面积值所对应的单位时间血流量。
一种血管血流速度的计算方法，该方法包括：确定血管的感兴趣区域；计算并拟合出所述感兴趣区域内的灰度拟合曲线；确定一预定时间区间内的最大灰度值曲线或最小灰度值曲线；计算预定时间区间内最大灰度值曲线或最小灰度值曲线与灰度拟合曲线所包围的区域面积值；基于该区域面积值，获取该面积值所对应的单位时间血流量；基于所述单位时间血流量以及血管管腔面积，获得所述血管的血流速度。
优选的，所述感兴趣区域包括注入造影剂的主支血管及其分支。
优选的，通过目标图像跟踪检测感兴趣区域在不同心跳时刻下的位置的变化，从而获得最佳的感兴趣区域。
优选的，所述方法进一步包括：接受血管的X射线造影影像序列，选择感兴趣区域；选择起始时间为造影剂充盈之前，提取每帧造影中感兴趣区域内灰度直方图，通过所述灰度直方图计算出每帧下感兴趣区域内的灰度值，并根据灰度值拟合出灰度随时间变化的灰度拟合曲线。
优选的，所述方法进一步包括：确定第一时间点，以及以第一时间点为中心的一预定时间区间内灰度拟合曲线的最大值和最小值；
优选的，所述最大灰度值曲线是以预定时间区间内灰度拟合曲线的最大值为纵坐标所做的曲线；所述最小灰度值曲线是以预定时间区间内灰度拟合曲线的最小值为纵坐标所做的曲线。
优选的，当所述灰度拟合曲线变化趋势为下降时，所述第一时间点为灰度拟合曲线中灰度值下降最快的时间点；当所述灰度拟合曲线变化趋势为上升时，所述第一时间点为灰度拟合曲线中灰度值上升最快的时间点。
优选的，当灰度拟合曲线变化趋势为下降时，计算灰度拟合曲线上每点的斜率，获取斜率为负值，且斜率绝对值最大的点，所述点即为灰度值下降最快的时间点；当灰度拟合曲线变化趋势为上升时，计算灰度拟合曲线上每点的斜率，获取斜率为正值，所述斜率的值最大的点即为灰度值上升最快的时间点。
优选的，所述区域面积值的计算过程进一步包括：当灰度拟合曲线变化趋势为下降时，获取所述灰度拟合曲线中的第一时间点，在所述以第一时间点为中心的预定时间区间内对灰度拟合曲线进行积分，计算该预定时间区间内最大灰度值曲线与所述灰度拟合曲线所包围的区域面积值；当灰度拟合曲线变化趋势为上升时，获取所述灰度拟合曲线中的第一时间点，在所述以第一时间点为中心的预定时间区间内对灰度拟合曲线进行积分，计算该预定时间区间内最小灰度值曲线与所述灰度拟合曲线所包围的区域面积值。
优选的，所述预定时间区间为整数个心动周期，所述整数大于等于1。
优选的，所述预定时间区间为一个心动周期，包括第一时间点前后的各二分之一心动周期；其中，所述第一时间点前后二分之一心动周期时间区间为造影剂开始充盈血管感兴趣区域后，在没有完全充盈感兴趣区域前的这一时间段。
优选的，通过查对应表，可获得该区域面积值所对应的单位时间血流量，所述对应表为不同面积值与单位时间不同血流量之间的对应表格。
优选的，所述血管管腔面积可由三维定量测量方法获得。
优选的，所述灰度拟合曲线的拟合公式为：
g(t)＝a0+a1t+a2t2+…+antn.其中a0，a1，a2，…an为拟合系数，t为时间。
优选的，根据心电图确定心动周期，或根据拟合前的原始灰度变化曲线计算得到的灰度值峰峰值之间的帧数m，计算心动周期T＝m/f，其中f代表造影的帧频率。
优选的，所述原始灰度变化曲线是由每帧造影影像中感兴趣区域内灰度直方图计算得到的灰度值所作的原始数据曲线。
优选的，通过接受心脏冠脉造影，获得血管感兴趣区域后，所获得的血流速度可用于评价血管狭窄对血流速度的影响，或者后续计算狭窄血管血流储备分数FFR值。
优选的，通过接受肾动脉造影影像，获得血管感兴趣区域后，所获得单位时间血流量或血流速度在肾动脉交感神经消融方面用于评估肾动脉交感神经消融前后的变化，或者用于实时评估消融时单位时间血流量和血流速度的变化。
优选的，通过接受基于有肿瘤区域的血管造影，获得血管感兴趣区域后，所获得的单位时间血流量或血流速度可用于评估肿瘤治疗前跟治疗后供血的变化，以提示治疗效果。
优选的，通过接受动脉血管造影，获得血管感兴趣区域后，所获得的血流速度可用于外周血管计算狭窄血管的压力降或血流储备分数(FFR)值。
本发明的有益效果在于，该技术方案提供了一种新的单位时间血流量与血流速度计算方法，保证了计算得到的单位时间血流量与血流速度是整数个心动周期血流速度均值，从而有效地避免了不恰当时间段的选择导致计算的平均单位时间血流量与血流速度是非整数心动周期内均值造成的计算错误。而利用感兴趣区域图像灰度值随时间的变化求出血流速度，不仅达到了无创诊断，同时能够有选择性地增加或者排除边支血流，以适应不同的应用。
附图说明
图1为冠脉造影图像灰度直方图；
图2A为造影剂充盈前血管灰度变化示意图；
图2B为造影剂充盈后血管灰度变化示意图；
图3为Doppler导丝法测得的不同心动周期内血流速度变化曲线示意图；
图4为感兴趣区域不同心动周期内的原始灰度变化曲线和灰度拟合曲线示意图；
图5为单位时间血流量与血流速度计算原理示意图。
具体实施方式
下面将结合本发明实施例中的附图1-5，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
本发明提供一种计算血管单位时间血流量及血流速度的方法，具体包括如下步骤：首先，确定血管的感兴趣区域(优选的一种方式为，通过接受血管的X射线造影，选择造影中感兴趣区域)；其次，计算并拟合出所述感兴趣区域内灰度拟合曲线；其次，获取预定时间区间内的最大灰度值曲线；其次，计算该预定时间区间内最大灰度值曲线与灰度拟合曲线所包围的区域面积值S；再次，基于该区域面积值S，获取该面积值所对应的单位时间血流量Q；最后，结合血管管腔面积，获得所述血管的血流速度V。
优选的，所述方法还可包括：确定第一时间点，以及以第一时间点为中心的一预定时间区间内灰度拟合曲线的最大值和最小值。优选的，所述第一时间点为灰度拟合曲线中灰度值下降最快的时间点。所述最大灰度值曲线是以预定时间区间内灰度拟合曲线的最大值为纵坐标所做的曲线。
优选的，可以通过目标图像跟踪配准检测感兴趣区域在不同心跳时刻下的位置的变化，从而获得最佳的感兴趣区域。需要强调的是，现有技术中的感兴趣区域不会随着造影序列目标血管位置移动而相应移位，导致感兴趣区域内的灰度值变化不仅仅由于造影剂在目标血管中传输导致，还可能由于目标血管移出和边支血管移入，从而导致计算的血流速度不准确。而造影序列目标血管位置移动十分常见，可由心脏的跳动、病人的呼吸和移动导致。
由于造影剂的充盈速度较快，充盈整个感兴趣区域通常时间较短，因此，优选选择上述预定时间区间内为整一个心动周期(T)，即在所述第一时间点前后的二分之一心动周期的区间内对灰度拟合曲线进行积分，得到面积值S。其中，在第一时间点前后二分之一心动周期时间区间进一步优选为造影剂开始充盈血管感兴趣区域后，在没有完全充盈感兴趣区域前的这一段时间段。在一优选实施例中，该时间段包含了心肌跟微循环灌注，因此，优选的，选定感兴趣区域则为包括该感兴趣血管段灌注的心肌。
通过查对应表，可获得该面积值S所对应的单位时间血流量Q，所述对应表为不同面积值与不同单位时间血流量之间的对应表格，所述表格可以通过多次、重复性的大量常规实验进行获得，并根据后期的实验数据进行更新该表格。
其中，在计算心动周期时，优选的计算方法为，根据心电图数据得到心动周期。或者根据直方图计算得到的原始灰度变化曲线峰峰值之间的帧数m，计算心动周期T＝m/f。其中f代表造影的帧频率。其中，原始灰度变化曲线是利用从每帧造影直方图上直接计算得到的灰度值所做的曲线，是原始数据曲线；灰度拟合曲线是根据原始数据通过拟合手段获得的曲线。
优选所述第一时间点为灰度拟合曲线中灰度值下降最快的时间点。可通过计算灰度拟合曲线上每点的斜率，获得斜率为负值，且绝对值最大的点，所述点则为灰度值下降最快的时间点。
优选的，由三维定量测量的方法，得到血管管腔面积A，平均血流速度V＝Q/A。
需要指出的时，本实施例采用了一般血管造影情况下的灰度值变化为下降趋势的情况进行了分析，即在所得到灰度拟合曲线为下降曲线的情况下，通过选择该曲线的灰度值下降最快的点为第一时间点，并以该第一时间点为中心确定上述预定时间区间，基于上述预定时间区间的灰度拟合曲线的最大值获得最大灰度值曲线，从而计算出以第一时间点为中心的预定时间区间内最大灰度值曲线与灰度拟合曲线所围成的曲线面积。
然而，在某些造影影像中，造影剂充盈后灰度值比充盈前的灰度值大，上述灰度值变化为上升趋势，即所得到灰度拟合曲线为上升曲线，此时，则需要检测灰度拟合曲线中上升最快的时间点为第一时间点；获取以第一时间点为中心的前后二分之一心动周期内灰度拟合曲线的最小值；以该最小值为纵坐标做出最小灰度值曲线；计算以第一时间点为中心的预定时间区间内最小灰度值曲线与灰度拟合曲线所围成的曲线面积。在该情况下，计算灰度拟合曲线上每点的斜率，获得斜率为正值，所述斜率的值最大的点即为灰度值上升最快的时间点。
实施例2
参见图1，冠脉造影利用人体软组织和造影剂对射线吸收程度的不同，在造影图像上形成了血管和周围组织之间的不同高对比度。造影图像中每一个像素的颜色深度由灰度值表示，灰度值越大像素越亮。灰度直方图是数字图像中最简单且有用的工具，它表示图像中具有某种灰度级的像素的个数，其横坐标是灰度值，取值范围优选为0-255，纵坐标表示图像中该灰度值出现的个数，取值范围优选为0-N，N为图像像素点个数。
如图2所示，我们选择包含狭窄的血管作为感兴趣区域，所述感兴趣区域包括注入造影剂的主支血管及其分支。血管在没有打入造影剂之前(图A)灰度值较高，无法从周围软组织中区分出来。打入造影剂之后(图B)造影剂随血流扩散，由于造影剂对射线的吸收能力更强，感兴趣区域的灰度值随之减小，血管颜色变深。几个心动周期之后造影剂被稀释，感兴趣区域灰度值升高。因此感兴趣区域灰度值的变化速度反映了管腔内的血流速度。
每个心动周期内平均血流速度大小近似，但不同时间段的选取对于计算平均血流速度有很大影响。如图3所示，用Doppler导丝法直接测量的不同心动周期内血流速度曲线。时间间隔相同的不同时间段T1和T2所求解出的平均血流速度相差很大，因此，为了确保计算值精确，优选的选取整数个心动周期进行血流速度均值计算，如整一个心动周期。
如图4所示，提取造影感兴趣区域每一帧的灰度值，并拟合出灰度拟合曲线g(t)。例如，优先选择3个以上(即N＞＝3)心动周期，起始时间为造影剂充盈之前，提取每帧造影中感兴趣区域内灰度直方图，通过所述灰度直方图计算出每帧下感兴趣区域内的灰度值，并根据灰度值拟合出灰度拟合曲线g(t)，所述拟合公式为多项式拟合：
g(t)＝a0+a1t+a2t2+…+antn；其中a0，a1，a2，…an为拟合系数，t为造影剂充盈血管的时间，所述时间是从第一帧影像采集开始计算出的时间。
一般情况，可以根据心电图数据得到心动周期。在缺少心电图数据时，可根据直方图计算得到的原始灰度曲线的峰峰值之间的帧数m，并通过公式T＝m/f计算获得心动周期，其中f代表造影的帧频率。
如图5，求出灰度值下降过程中斜率绝对值最大的点(t0，g(t0))，所述点确定为第一时间点。计算第一时间点前后分别二分之一心动周期[t1，t2]内曲线g(t)与最大灰度值曲线g(t1)所包围的区域(如图中的阴影区域)面积值S。其中，所述最大灰度值曲线g(t1)是以[t1，t2]时间段内曲线g(t)的最大值为纵坐标所做的曲线；该阴影面积值S与一个心动周期血流量Q成正比，即S∝Q。
实施例3
需要指出的是，本发明的实施例中所采用的X射线造影可以是心脏冠脉造影，外周动脉造影如肾动脉造影、颈动脉造影等，或者是有肿瘤治疗前后的造影。基于上述不同的造影方式所获得单位时间血流量或血流速度可用于不同的病症情况分析中的关键参数指标，并获得比现有技术中的参数指标更好的准确性及精准度。例如，基于心脏冠脉造影计算血流速度可以用来评价血管狭窄对血流速度的影响，以及后续计算狭窄血管的压力差或血流储备分数(FFR)值；基于肾动脉造影可以在肾动脉交感神经消融方面用于评估肾动脉单位时间血流量在交感神经消融前后的变化，或者实时评估消融时单位时间血流量和血流速度的变化情况，以提示消融的效果；基于有肿瘤区域的血管造影计算单位时间血流量或血流速度可以评估肿瘤治疗前跟治疗后供血的变化，以提示治疗效果。
在一具体实施例中，本发明提供了一种某段血管的血流储备分数FFR的计算方法，基于本发明中的计算血流速度的方法所获得的平均血流速度或最大平均血流速度，并结合该段血管的其它几何参数，通过相应的计算公式获得血管的压力降或FFR值。所述方法包括：接收该段血管的几何参数，该血管包括近端终点和远端终点，所述几何参数包括第一几何参数，代表该血管段近端横截面的面积(或直径)；第二几何参数，代表该血管段远端横截面的面积(或直径)；第三几何参数，代表该血管段位于近端终点和远端终点之间的第一位置的横截面面积(或直径)；以近端终点为参考点，基于所述几何参数和血管段上的点到参考点的距离，计算出参考管腔直径函数和几何参数差异函数；在多个尺度下对几何参数差异函数求导，得到多个尺度对应的差值导数函数；所述尺度指的是分辨率，即数值计算导数时相邻两点之间的距离；利用本发明实施例中的计算血流速度的方法获得该血管段在常规冠脉造影的平均血流速度计算出其对应的最大平均血流速度；基于多尺度差值导数函数和最大平均血流速度，获得所述血管第一位置处的第二血流压力与近端终点处的第一血流压力之间的比值，即血流储备分数。
本发明的创新点之一就在于，以灰度拟合曲线变化最快位置为中心，计算整一个心动周期时间内的血流量，从而更准确地计算出单位时间血流量与血流速度，有效地避免了不恰当时间区间进行计算所导致的误差。本发明的有益效果在于，该技术方案提供了一种新的单位时间血流量与血流速度计算方法，保证了计算得到的单位时间血流量与血流速度是整数个心动周期血流速度均值，从而有效地避免了不恰当时间段的选择导致计算的单位时间血流量与平均血流速度是非整数心动周期内均值造成的计算错误。而利用感兴趣区域图像灰度值随时间的变化求出血流速度，不仅达到了无创诊断，同时能够有选择性地增加或者排除边支血流，以适应不同的应用。
虽然本发明已以较佳实施例揭示如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的修改和完善，因此本发明的保护范围当以权利要求书所界定的为准。
Claims (20)
1.一种血管单位时间血流量的计算方法，该方法包括：
确定血管的感兴趣区域；
计算并拟合出所述感兴趣区域内的灰度拟合曲线；
确定一预定时间区间内的最大灰度值曲线或最小灰度值曲线；
计算预定时间区间内最大灰度值曲线或最小灰度值曲线与灰度拟合曲线所包围的区域面积值；
基于该区域面积值，获取该面积值所对应的单位时间血流量。
2.一种血管血流速度的计算方法，该方法包括：
确定血管的感兴趣区域；
计算并拟合出所述感兴趣区域内的灰度拟合曲线；
确定一预定时间区间内的最大灰度值曲线或最小灰度值曲线；
计算预定时间区间内最大灰度值曲线或最小灰度值曲线与灰度拟合曲线所包围的区域面积值；
基于该区域面积值，获取该面积值所对应的单位时间血流量；
基于所述单位时间血流量以及血管管腔面积，获得所述血管的血流速度。
3.根据权利要求1或2所述的方法，其特征在于，所述感兴趣区域包括注入造影剂的主支血管及其分支。
4.根据权利要求1或2所述的方法，其特征在于，通过目标图像跟踪检测感兴趣区域在不同心跳时刻下的位置的变化，从而获得最佳的感兴趣区域。
5.根据权利要求1或2所述的方法，其特征在于，所述方法进一步包括：
接受血管的X射线造影影像序列，选择感兴趣区域；选择起始时间为造影剂充盈之前，计算每帧造影中感兴趣区域内总灰度值，并根据灰度值拟合出灰度随时间变化的灰度拟合曲线。
6.根据权利要求1或2所述的方法，其特征在于，所述方法进一步包括：
确定第一时间点，以及以第一时间点为中心的一预定时间区间内灰度拟合曲线的最大值和最小值；
优选的，所述最大灰度值曲线是以预定时间区间内灰度拟合曲线的最大值为纵坐标所做的曲线；所述最小灰度值曲线是以预定时间区间内灰度拟合曲线的最小值为纵坐标所做的曲线。
7.根据权利要求6所述的方法，其特征在于，当所述灰度拟合曲线变化趋势为下降时，所述第一时间点为灰度拟合曲线中灰度值下降最快的时间点；当所述灰度拟合曲线变化趋势为上升时，所述第一时间点为灰度拟合曲线中灰度值上升最快的时间点。
8.根据权利要求7所述的方法，其特征在于，当灰度拟合曲线变化趋势为下降时，计算灰度拟合曲线上每点的斜率，获取斜率为负值，且斜率绝对值最大的点，所述点即为灰度值下降最快的时间点；当灰度拟合曲线变化趋势为上升时，计算灰度拟合曲线上每点的斜率，获取斜率为正值，所述斜率的值最大的点即为灰度值上升最快的时间点。
9.根据权利要求7或8所述的方法，其特征在于，所述区域面积值的计算过程进一步包括：
当灰度拟合曲线变化趋势为下降时，获取所述灰度拟合曲线中的第一时间点，在所述以第一时间点为中心的预定时间区间内对灰度拟合曲线进行积分，计算该预定时间区间内最大灰度值曲线与所述灰度拟合曲线所包围的区域面积值；
当灰度拟合曲线变化趋势为上升时，获取所述灰度拟合曲线中的第一时间点，在所述以第一时间点为中心的预定时间区间内对灰度拟合曲线进行积分，计算该预定时间区间内最小灰度值曲线与所述灰度拟合曲线所包围的区域面积值。
10.根据权利要求1-9中任一所述的方法，其特征在于，所述预定时间区间为整数个心动周期，所述整数大于等于1。
11.根据权利要求10所述的方法，其特征在于，所述预定时间区间为一个心动周期，包括第一时间点前后的各二分之一心动周期；其中，所述第一时间点前后二分之一心动周期时间区间为造影剂开始充盈血管感兴趣区域后，在没有完全充盈感兴趣区域前的这一时间段。
12.根据权利要求1-11中任一所述的方法，其特征在于，通过查对应表，可获得该区域面积值所对应的单位时间血流量，所述对应表为不同面积值与单位时间不同血流量之间的对应表格。
13.根据权利要求1-12中任一所述的方法，其特征在于，所述血管管腔面积可由三维定量测量方法获得。
14.根据权利要求1-13中任一所述的方法，其特征在于，所述灰度拟合曲线的拟合公式为：
g(t)＝a0+a1t+a2t2+…+antn；其中a0，a1，a2，…an为拟合系数，t为时间。
15.根据权利要求1-14中任一所述的方法，其特征在于，根据心电图确定心动周期，或根据拟合前的原始灰度变化曲线计算得到的灰度值峰峰值之间的帧数m，计算心动周期T＝m/f，其中f代表造影的帧频率。
16.根据权利要求15所述的方法，其特征在于，所述原始灰度变化曲线是由每帧造影影像中感兴趣区域内灰度直方图计算得到的灰度值所作的原始数据曲线。
17.根据权利要求2-16中任一所述的方法，其特征在于，通过接受心脏冠脉造影，获得血管感兴趣区域后，所获得的血流速度可用于评价血管狭窄对血流速度的影响，或者后续计算狭窄血管血流储备分数FFR值。
18.根据权利要求1-16中任一所述的方法，其特征在于，通过接受肾动脉造影影像，获得血管感兴趣区域后，所获得单位时间血流量或血流速度在肾动脉交感神经消融方面用于评估肾动脉交感神经消融前后的变化，或者用于实时评估消融时单位时间血流量和血流速度的变化，以提示消融治疗效果。
19.根据权利要求1-16中任一所述的方法，其特征在于，通过接受基于有肿瘤区域的血管造影，获得血管感兴趣区域后，所获得的单位时间血流量或血流速度可用于评估肿瘤治疗前跟治疗后供血的变化，以提示治疗效果。
20.一种某段血管的压力降或血流储备分数(FFR)的计算方法，其基于所述权利要求2-16中任一所述的计算血流速度方法计算该段血管的平均血流速度，并结合该段血管的其它几何参数，通过相应的计算公式获得该段血管的压力降或FFR值。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510916119.8A CN105559810B (zh) | 2015-12-10 | 2015-12-10 | 血管单位时间血流量与血流速度的计算方法 |
| PCT/CN2016/104655 WO2017097074A1 (zh) | 2015-12-10 | 2016-11-04 | 血管单位时间血流量与血流速度的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510916119.8A CN105559810B (zh) | 2015-12-10 | 2015-12-10 | 血管单位时间血流量与血流速度的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105559810A true CN105559810A (zh) | 2016-05-11 |
| CN105559810B CN105559810B (zh) | 2017-08-08 |
Family
ID=55870860
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510916119.8A Active CN105559810B (zh) | 2015-12-10 | 2015-12-10 | 血管单位时间血流量与血流速度的计算方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN105559810B (zh) |
| WO (1) | WO2017097074A1 (zh) |
Cited By (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106537392A (zh) * | 2014-04-22 | 2017-03-22 | 西门子保健有限责任公司 | 用于冠状动脉中的血液动力学计算的方法和系统 |
| CN106805989A (zh) * | 2017-03-13 | 2017-06-09 | 博动医学影像科技（上海）有限公司 | 用于动脉造影的图像处理系统及交感神经状态变化的快速测量系统 |
| WO2017097074A1 (zh) * | 2015-12-10 | 2017-06-15 | 博动医学影像科技（上海）有限公司 | 血管单位时间血流量与血流速度的计算方法 |
| CN108245178A (zh) * | 2018-01-11 | 2018-07-06 | 苏州润迈德医疗科技有限公司 | 一种基于x射线冠脉造影图像的血液流动速度计算方法 |
| CN108550388A (zh) * | 2018-01-12 | 2018-09-18 | 深圳科亚医疗科技有限公司 | 基于血管造影的计算血管血流参数的装置和系统 |
| CN108550189A (zh) * | 2018-05-03 | 2018-09-18 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
| CN109512450A (zh) * | 2018-10-18 | 2019-03-26 | 深圳市孙逸仙心血管医院（深圳市心血管病研究所） | 测量血管血流速度的方法 |
| CN109620199A (zh) * | 2018-11-30 | 2019-04-16 | 博动医学影像科技（上海）有限公司 | 建立血管截面函数、血管压力差和血管应力的方法及装置 |
| CN109717861A (zh) * | 2018-12-28 | 2019-05-07 | 华中科技大学鄂州工业技术研究院 | 一种血流速度的激光散斑检测方法及检测装置 |
| CN110226923A (zh) * | 2018-03-05 | 2019-09-13 | 苏州润迈德医疗科技有限公司 | 一种无需血管扩张剂测量血流储备分数的方法 |
| CN110367965A (zh) * | 2018-09-19 | 2019-10-25 | 苏州润迈德医疗科技有限公司 | 便捷测量冠状动脉血管评定参数的方法、装置及系统 |
| CN110384493A (zh) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | 测量微循环阻力指数的系统以及冠脉分析系统 |
| CN110811593A (zh) * | 2019-11-07 | 2020-02-21 | 天津中医药大学第二附属医院 | 一种高血压中医证型采取肝血流动力学特点的判别方法 |
| CN111227822A (zh) * | 2018-11-28 | 2020-06-05 | 苏州润心医疗器械有限公司 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
| CN111369519A (zh) * | 2020-03-02 | 2020-07-03 | 博动医学影像科技（上海）有限公司 | 冠状动脉的血流速度的计算方法、装置及电子设备 |
| CN111627002A (zh) * | 2020-05-25 | 2020-09-04 | 上海杏脉信息科技有限公司 | 一种冠脉微血管阻力指数计算装置及方法 |
| WO2021017327A1 (zh) * | 2019-07-31 | 2021-02-04 | 苏州润迈德医疗科技有限公司 | 基于造影图像间隔时间修正血流速度的方法及装置 |
| CN113180631A (zh) * | 2021-04-29 | 2021-07-30 | 博动医学影像科技（上海）有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
| RU2756367C1 (ru) * | 2018-02-15 | 2021-09-29 | Университа' Дельи Студи Ди Рома "Ла Сапиенца" | Способ и система для измерения показателя текущего периферического сопротивления |
| CN114376546A (zh) * | 2020-12-28 | 2022-04-22 | 深圳北芯生命科技股份有限公司 | 一种支持双诊断模式的系统 |
| CN115337096A (zh) * | 2021-05-13 | 2022-11-15 | 苏州润迈德医疗科技有限公司 | 用于内脏动脉去交感神经的消融方法、系统及存储介质 |
| CN115363597A (zh) * | 2022-08-30 | 2022-11-22 | 杭州朗博康医疗科技有限公司 | 一种医疗图像显示方法、装置、电子设备及存储介质 |
| CN115998275A (zh) * | 2022-12-26 | 2023-04-25 | 广东省新黄埔中医药联合创新研究院 | 血流速度检测校准方法、装置、设备及可读存储介质 |
| CN116342603A (zh) * | 2023-05-30 | 2023-06-27 | 杭州脉流科技有限公司 | 获得动脉输入函数的方法 |
| WO2023130197A1 (en) * | 2022-01-04 | 2023-07-13 | Shenzhen Xpectvision Technology Co., Ltd. | Flow speed measurements using imaging systems |
| CN118078330A (zh) * | 2024-04-24 | 2024-05-28 | 天津市胸科医院 | 基于深度学习的超声影像血流识别分析方法及系统 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107949311B (zh) | 2015-04-16 | 2021-04-16 | Gentuity有限责任公司 | 用于神经病学的微光探针 |
| US10631718B2 (en) | 2015-08-31 | 2020-04-28 | Gentuity, Llc | Imaging system includes imaging probe and delivery devices |
| EP3700406A4 (en) | 2017-11-28 | 2021-12-29 | Gentuity LLC | Imaging system |
| CN113012108B (zh) * | 2018-07-24 | 2024-02-20 | 上海博动医疗科技股份有限公司 | 血管图像的处理方法、装置及成像设备 |
| CN115984301B (zh) * | 2022-12-13 | 2023-07-11 | 徐州医科大学 | 一种基于布加综合征血管图像三维重建精度对比方法 |
| CN116664564B (zh) * | 2023-07-28 | 2023-10-31 | 杭州脉流科技有限公司 | 基于颅内医学影像获取血流量的方法和装置 |
| CN118314116A (zh) * | 2024-04-30 | 2024-07-09 | 中国海洋大学 | 一种基于图像灰度值的侏儒蛤器官重量预测方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050178732A1 (en) * | 1994-09-16 | 2005-08-18 | Transonic Systems, Inc. | Measurement of a blood flow rate in hemodialysis shunts |
| US20070135716A1 (en) * | 2005-12-13 | 2007-06-14 | Up Management Gmbh & Co Med-Systems Kg | Dilution apparatus, method and computer program |
| WO2015041312A1 (ja) * | 2013-09-20 | 2015-03-26 | 国立大学法人旭川医科大学 | 血管内血流動態の画像処理方法及びシステム |
| CN104688259A (zh) * | 2015-02-16 | 2015-06-10 | 袁小东 | Ct扫描灌注方法及装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001204700A (ja) * | 2000-01-28 | 2001-07-31 | Star Medical Kk | 血流量測定方法と血流量測定装置 |
| DE60212917T2 (de) * | 2001-10-16 | 2007-03-01 | Kabushiki Kaisha Toshiba | Vorrichtung zur Berechnung eines Index von örtlichen Blutflüssen |
| JP3920140B2 (ja) * | 2002-05-13 | 2007-05-30 | 株式会社東芝 | Ｍｒｉ装置及びフロー定量化装置 |
| JP4013128B2 (ja) * | 2002-09-12 | 2007-11-28 | 株式会社日立メディコ | 血流動態解析装置、方法、及び画像診断装置 |
| CN105559810B (zh) * | 2015-12-10 | 2017-08-08 | 博动医学影像科技（上海）有限公司 | 血管单位时间血流量与血流速度的计算方法 |
-
2015
- 2015-12-10 CN CN201510916119.8A patent/CN105559810B/zh active Active
-
2016
- 2016-11-04 WO PCT/CN2016/104655 patent/WO2017097074A1/zh active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050178732A1 (en) * | 1994-09-16 | 2005-08-18 | Transonic Systems, Inc. | Measurement of a blood flow rate in hemodialysis shunts |
| US20070135716A1 (en) * | 2005-12-13 | 2007-06-14 | Up Management Gmbh & Co Med-Systems Kg | Dilution apparatus, method and computer program |
| WO2015041312A1 (ja) * | 2013-09-20 | 2015-03-26 | 国立大学法人旭川医科大学 | 血管内血流動態の画像処理方法及びシステム |
| CN104688259A (zh) * | 2015-02-16 | 2015-06-10 | 袁小东 | Ct扫描灌注方法及装置 |
Cited By (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106537392B (zh) * | 2014-04-22 | 2019-07-26 | 西门子保健有限责任公司 | 用于冠状动脉中的血液动力学计算的方法和系统 |
| US10134129B2 (en) | 2014-04-22 | 2018-11-20 | Siemens Healthcare Gmbh | Method and system for hemodynamic computation in coronary arteries |
| CN106537392A (zh) * | 2014-04-22 | 2017-03-22 | 西门子保健有限责任公司 | 用于冠状动脉中的血液动力学计算的方法和系统 |
| WO2017097074A1 (zh) * | 2015-12-10 | 2017-06-15 | 博动医学影像科技（上海）有限公司 | 血管单位时间血流量与血流速度的计算方法 |
| CN106805989A (zh) * | 2017-03-13 | 2017-06-09 | 博动医学影像科技（上海）有限公司 | 用于动脉造影的图像处理系统及交感神经状态变化的快速测量系统 |
| CN108245178A (zh) * | 2018-01-11 | 2018-07-06 | 苏州润迈德医疗科技有限公司 | 一种基于x射线冠脉造影图像的血液流动速度计算方法 |
| CN108550388A (zh) * | 2018-01-12 | 2018-09-18 | 深圳科亚医疗科技有限公司 | 基于血管造影的计算血管血流参数的装置和系统 |
| RU2756367C1 (ru) * | 2018-02-15 | 2021-09-29 | Университа' Дельи Студи Ди Рома "Ла Сапиенца" | Способ и система для измерения показателя текущего периферического сопротивления |
| CN110226923A (zh) * | 2018-03-05 | 2019-09-13 | 苏州润迈德医疗科技有限公司 | 一种无需血管扩张剂测量血流储备分数的方法 |
| CN108550189A (zh) * | 2018-05-03 | 2018-09-18 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
| WO2019210553A1 (zh) * | 2018-05-03 | 2019-11-07 | 苏州润迈德医疗科技有限公司 | 基于造影图像和流体力学模型的微循环阻力指数计算方法 |
| CN110367965A (zh) * | 2018-09-19 | 2019-10-25 | 苏州润迈德医疗科技有限公司 | 便捷测量冠状动脉血管评定参数的方法、装置及系统 |
| CN110384493A (zh) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | 测量微循环阻力指数的系统以及冠脉分析系统 |
| CN109512450A (zh) * | 2018-10-18 | 2019-03-26 | 深圳市孙逸仙心血管医院（深圳市心血管病研究所） | 测量血管血流速度的方法 |
| CN111227822A (zh) * | 2018-11-28 | 2020-06-05 | 苏州润心医疗器械有限公司 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
| CN111227822B (zh) * | 2018-11-28 | 2022-02-11 | 苏州润迈德医疗科技有限公司 | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 |
| CN109620199A (zh) * | 2018-11-30 | 2019-04-16 | 博动医学影像科技（上海）有限公司 | 建立血管截面函数、血管压力差和血管应力的方法及装置 |
| US11445923B2 (en) | 2018-11-30 | 2022-09-20 | Shanghai Pulse Medical Technology, Inc. | Method and device for establishing blood vessel cross-section function, blood stress vessel pressure difference and blood vessel stress |
| CN109717861A (zh) * | 2018-12-28 | 2019-05-07 | 华中科技大学鄂州工业技术研究院 | 一种血流速度的激光散斑检测方法及检测装置 |
| WO2021017327A1 (zh) * | 2019-07-31 | 2021-02-04 | 苏州润迈德医疗科技有限公司 | 基于造影图像间隔时间修正血流速度的方法及装置 |
| CN110811593A (zh) * | 2019-11-07 | 2020-02-21 | 天津中医药大学第二附属医院 | 一种高血压中医证型采取肝血流动力学特点的判别方法 |
| CN111369519A (zh) * | 2020-03-02 | 2020-07-03 | 博动医学影像科技（上海）有限公司 | 冠状动脉的血流速度的计算方法、装置及电子设备 |
| CN111627002A (zh) * | 2020-05-25 | 2020-09-04 | 上海杏脉信息科技有限公司 | 一种冠脉微血管阻力指数计算装置及方法 |
| CN111627002B (zh) * | 2020-05-25 | 2023-07-18 | 上海杏脉信息科技有限公司 | 一种冠脉微血管阻力指数计算装置及方法 |
| CN114376546A (zh) * | 2020-12-28 | 2022-04-22 | 深圳北芯生命科技股份有限公司 | 一种支持双诊断模式的系统 |
| CN114376546B (zh) * | 2020-12-28 | 2024-03-12 | 深圳北芯生命科技股份有限公司 | 一种支持双诊断模式的系统 |
| CN113180631A (zh) * | 2021-04-29 | 2021-07-30 | 博动医学影像科技（上海）有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
| CN115337096A (zh) * | 2021-05-13 | 2022-11-15 | 苏州润迈德医疗科技有限公司 | 用于内脏动脉去交感神经的消融方法、系统及存储介质 |
| WO2023130197A1 (en) * | 2022-01-04 | 2023-07-13 | Shenzhen Xpectvision Technology Co., Ltd. | Flow speed measurements using imaging systems |
| CN115363597A (zh) * | 2022-08-30 | 2022-11-22 | 杭州朗博康医疗科技有限公司 | 一种医疗图像显示方法、装置、电子设备及存储介质 |
| CN115998275B (zh) * | 2022-12-26 | 2023-10-24 | 广东省新黄埔中医药联合创新研究院 | 血流速度检测校准方法、装置、设备及可读存储介质 |
| CN115998275A (zh) * | 2022-12-26 | 2023-04-25 | 广东省新黄埔中医药联合创新研究院 | 血流速度检测校准方法、装置、设备及可读存储介质 |
| CN116342603A (zh) * | 2023-05-30 | 2023-06-27 | 杭州脉流科技有限公司 | 获得动脉输入函数的方法 |
| CN116342603B (zh) * | 2023-05-30 | 2023-08-29 | 杭州脉流科技有限公司 | 获得动脉输入函数的方法 |
| CN118078330A (zh) * | 2024-04-24 | 2024-05-28 | 天津市胸科医院 | 基于深度学习的超声影像血流识别分析方法及系统 |
| CN118078330B (zh) * | 2024-04-24 | 2024-06-21 | 天津市胸科医院 | 基于深度学习的超声影像血流识别分析方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2017097074A1 (zh) | 2017-06-15 |
| CN105559810B (zh) | 2017-08-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105559810B (zh) | 血管单位时间血流量与血流速度的计算方法 | |
| CN107978371B (zh) | 快速计算微循环阻力的方法及系统 | |
| US10803994B2 (en) | Vascular flow assessment | |
| CN111134651B (zh) | 基于腔内影像计算血流储备分数的方法、装置、系统以及计算机存储介质 | |
| JP7136493B2 (ja) | 圧力センサと造影画像に基づき冠血流予備量比を計算する方法 | |
| JP7118464B2 (ja) | 血管圧差を取得する方法及び装置 | |
| JP2017518844A (ja) | 血流予備量比値を算出するための機器 | |
| JPH03184531A (ja) | 血流量測定装置 | |
| CN108742587B (zh) | 基于病史信息获取血流特征值的方法及装置 | |
| CN112384136A (zh) | 基于造影图像获取血管评定参数的方法、装置及系统 | |
| CN111067494A (zh) | 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法 | |
| US20190076196A1 (en) | Vessel geometry and additional boundary conditions for hemodynamic ffr/ifr simulations from intravascular imaging | |
| US20160029901A1 (en) | System and Method for Determining Arterial Compliance and Stiffness | |
| US20160302672A1 (en) | System and Method for Determining Arterial Compliance and Stiffness | |
| CN109009061B (zh) | 基于血压修正获取血流特征值的计算方法及装置 | |
| CN112673433A (zh) | 基于心肌呈色特性计算用于虚拟FFR和iFR计算的边界条件 | |
| US20200390345A1 (en) | Measurement guidance for coronary flow estimation from bernoulli's principle | |
| WO2022237161A1 (zh) | 脑动脉波强度和波功率的测量方法、终端设备及存储介质 | |
| CN108742570B (zh) | 基于冠脉优势类型获取血管压力差的装置 | |
| WO2016174010A1 (en) | Fractional flow reserve determination | |
| CN111091913A (zh) | 基于血流储备分数和冠状动脉ct造影图像的微循环阻力计算方法 | |
| CN114664455A (zh) | 一种冠状动脉血流储备分数计算方法及装置 | |
| JP7025536B2 (ja) | 血管灌流イメージングによる血流の評価のためのデバイス、システム、及び方法 | |
| CN113180614B (zh) | 无导丝ffr、无导丝imr和无导丝cfr的检测方法 | |
| JP2020534058A (ja) | シミュレートされた血行動態のための血管分岐部への流れの推定 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20160926 Address after: 200240, 41 floor, building 3, 398 Heqing Road, Minhang District, Shanghai, Applicant after: Medical imaging technology (Shanghai) Co., Ltd. Address before: 200240 Dongchuan Road, Shanghai, No. 800, No. Applicant before: Shanghai Jiao Tong University Applicant before: Medical imaging technology (Shanghai) Co., Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 200233 16th floor, building 82, No. 1198, Qinzhou North Road, Xuhui District, Shanghai Patentee after: Shanghai Bodong Medical Technology Co.,Ltd. Address before: 200240 room g3054, floor 3, building 41, No. 398, Heqing Road, Minhang District, Shanghai Patentee before: PULSE MEDICAL IMAGING TECHNOLOGY (SHANGHAI) Co.,Ltd. |