CN112494020B - 获取感兴趣的主动脉压力曲线的方法及存储介质 - Google Patents
获取感兴趣的主动脉压力曲线的方法及存储介质 Download PDFInfo
- Publication number
- CN112494020B CN112494020B CN202011348579.2A CN202011348579A CN112494020B CN 112494020 B CN112494020 B CN 112494020B CN 202011348579 A CN202011348579 A CN 202011348579A CN 112494020 B CN112494020 B CN 112494020B
- Authority
- CN
- China
- Prior art keywords
- pressure
- aortic
- acquiring
- list
- curve
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000035488 systolic blood pressure Effects 0.000 claims description 21
- 230000035487 diastolic blood pressure Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000004321 preservation Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 239000003086 colorant Substances 0.000 claims description 2
- 210000000709 aorta Anatomy 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 210000004351 coronary vessel Anatomy 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000036262 stenosis Effects 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000002586 coronary angiography Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000000544 hyperemic effect Effects 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 150000002632 lipids Chemical group 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Vascular Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本申请提供了一种获取感兴趣的主动脉压力曲线的方法及存储介质，包括：获取主动脉压力Pa数据；根据所述主动脉压力Pa提取平稳压力波形；根据所述平稳压力波形生成主动脉压力波形图像；获取感兴趣的时刻；根据所述感兴趣的时刻从所述主动脉压力波形图像上选取相对应的主动脉压力曲线。本申请能从连续的主动脉压力曲线中得到感兴趣的主动脉压力曲线部分，通过这种方式，本发明能够避免客观环境的影响，同时克服了临床上分析主动脉压力曲线只能依靠观察和经验的弊端。
Description
技术领域
本发明涉及冠状动脉医学技术领域，特别是涉及一种获取感兴趣的主动脉压力曲线的方法及存储介质。
背景技术
人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块，继而导致血管狭窄；特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足，诱发冠心病、心绞痛等病症，对人类的健康造成严重威胁。据统计，我国现有冠心病患者约1100万人，心血管介入手术治疗患者数量每年增长大于10％。
冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度，但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性，1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve，FFR)，经过长期的基础与临床研究，FFR已成为冠脉狭窄功能性评价的金标准。
血流储备分数(FFR)通常是指心肌血流储备分数，设置为病变冠脉能为心肌提供的最大血流与该冠脉完全正常时最大供血流量之比，研究表明，在冠脉最大充血状态下，血流量的比值可以用压力值来代替。即FFR值的测量可在冠脉最大充血状态下，通过压力传感器对冠脉远端狭窄处的压力和冠脉狭窄近端压力进行测定继而计算得出。
在当前的设备中，主动脉压力曲线是通过压力传感器获取并数模转换后得到的。实际使用场景下，由于客观环境的影响，通过传感器直接获取平滑、稳定的主动脉压力曲线是不可能的，总会有一些不稳定的因素存在，导致压力曲线波动、失真或者间断。想要得到能够进行参数运算的压力波形，必须通过一系列的手段对波形进行处理，才能获得到感兴趣的主动脉压力曲线。
发明内容
本发明提供了一种获取感兴趣的主动脉压力曲线的方法及存储介质，以解决在客观环境下无法直接通过传感器直接获取平滑、稳定的主动脉压力曲线的问题。
为实现上述目的，第一方面，本申请提供了一种获取感兴趣的主动脉压力曲线的方法，包括：
获取主动脉压力Pa数据；
根据所述主动脉压力Pa提取平稳压力波形；
根据所述平稳压力波形生成主动脉压力波形图像；
获取感兴趣的时刻；
根据所述感兴趣的时刻从所述主动脉压力波形图像上选取相对应的主动脉压力曲线。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，根据所述感兴趣的时刻从所述主动脉压力波形图像上选取相对应的主动脉压力曲线的方法，包括：
从所述感兴趣的时刻开始，在所述主动脉压力波形图像上向前获取若干心跳周期内的压力波形段；
获取所述若干心跳周期内的压力波形段上的主动脉压力的平均值；
根据所述平均值去除偏差较大的压力值点，根据所述贝塞尔曲线重新绘制平滑的主动脉压力曲线。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述若干心跳周期内的压力波形段的为6～10个心跳周期内的压力波形段。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述平均值去除偏差较大的压力值点，根据所述贝塞尔曲线重新绘制平滑的主动脉压力曲线的方法包括：
根据所述平均值去除偏差较大的压力值点，从所述重新绘制平滑的主动脉压力曲线上保留压力最平稳的一个心跳周期或者两个心跳周期的压力曲线，作为感兴趣的主动脉压力曲线。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述主动脉压力Pa提取平稳压力波形的方法包括：
设置多个索引列表，包括：暂存列表、保存列表、均值点索引列表、极大值点索引列表、极小值点索引列表；
设置生理参数阈值，包括：心跳周期数量C、舒张压至收缩压范围、最小压力幅度阈值、收缩压或舒张压的波动差、心率范围、压力个数N阈、文件数量上限、截止压力阈值P截、截止心率阈值HR截、压力传输速率阈值V阈；
根据所述索引列表、所述生理参数阈值，所述主动脉压力数据计算平均压力；
根据所述平均压力获取压力均值点；
根据所述压力均值点分割数据周期，获取每个心跳周期内的压力参数；
根据所述压力参数，获取平稳压力波形。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述索引列表、所述生理参数阈值，所述主动脉压力数据计算平均压力的方法包括：
根据所述心率阈值获得一个心跳周期时间阈值；
根据一个所述心跳周期数量C和压力传输速率阈值获取一个所述心跳周期内的数据点个数范围；
根据所述数据点个数范围划定一个时间周期，每结束一个时间周期即对所有数据点个数进行统计，且将所有获取的主动脉压力与获取的时间一一对应的存储到所述暂存列表内；
如果所述暂存列表内的数据个数是N阈的整数倍，且大于数据点个数范围的最小阈值，且小于数据点个数范围的最大阈值，则计算所有暂存列表内的主动脉压力的平均值，即为平均压力。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述索引列表、所述生理参数阈值，所述主动脉压力数据计算平均压力的方法还包括：
如果所述暂存列表内连续n个主动脉压力Pa均小于截止压力阈值P截，则从所述暂存列表内去除所述n个主动脉压力Pa。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述平均压力获取压力均值点的方法，包括：
设置均值点索引列表；
向所述均值点索引列表内加入压力均值点；
根据压力均值点生成压力均值点曲线；
过滤所述压力均值点曲线上的压力均值点；
判断所述均值点索引列表内的所述压力均值点的数量是否满足心跳周期要求，如果不满足则调整n的大小，直至满足心跳周期要求。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述向所述均值点索引列表内加入压力均值点的方法包括：
将所述暂存列表内的主动脉压力Pa分别与所述平均压力进行比较；
如果且/>则将PK加入所述均值点索引列表中，其中PK、PK+1分别表示所述暂存列表内的第K个、第K+1个主动脉压力；
如果且/>则将PK加入所述均值点索引列表中。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述过滤所述压力均值点曲线上的压力均值点的方法包括：
根据Dn＝Pn+1-Pn，其中Pn、Pn+1分别表示所述均值点索引列表内第n个、第n+1个主动脉压力，所述Dn表示Pn+1、Pn的压力差；
根据HRn＝30/(Dn/n)，获取所述第n个半心跳周期内的心率频率值HRn；
如果HRn＜HR截，则从所述均值点索引列表中删除所述HRn对应的均值压力点；
如果且在Pn后连续m个主动脉压力均大于/>且/>则从所述均值点索引列表中删除Pn后连续m个压力均值点。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述判断所述均值点索引列表内的所述压力均值点的数量是否满足心跳周期要求，如果不满足则重复所述权利要求6，调整n的大小，直至满足心跳周期要求的方法包括：
如果均值点列表内的数据个数N均＜2C+1，则不满足心跳周期要求，调整n的大小，直至N均≥2C+1。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述压力均值点分割数据周期，获取每个心跳周期内的压力参数的方法，包括：
如果从所述暂存列表内提取Pn所在心跳周期内的最大压力Pmax和所述最大压力Pmax所在的位置Smax，将所述Smax、Pmax一一对应的均加入舒张压列表中；
获取所述舒张压列表中的Pmax的平均值，获得舒张压P舒；
如果从所述暂存列表内提取Pn所在心跳周期内的最小压力Pmin和所述最小压力Pmin所在的位置Smin，将所述Smin、Pmin一一对应的均加入收缩压列表中；
获取所述收缩压列表中的Pmin的平均值，获得收缩压P收。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述压力参数，获取平稳压力波形的方法，包括：
设置偏移量为V阈/a，截取起始点O、结束点E；
将所述均值点索引列表、舒张压列表、收缩压列表全部按照上述偏移量、起始点O、结束点E执行，输出所述平稳压力波形、舒张压波形和收缩压波形。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述生成主动脉压力波形图像的方法包括：
设置所述平稳压力波形的显示区间；
设置所述平稳压力波形的数据范围；
根据所述显示区间、数据范围以及平稳压力波形生成主动脉压力波形图像。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述设置所述平稳压力波形的显示区间的方法包括：
设置所述生成的主动脉压力波形图像的长度W1和宽度H1；
设置所述主动脉压力波形图像的外边框长度W2和宽度H2；
根据W＝W1-2W2，H＝H1-2H2，获得所述平稳压力波形的显示区间的长度W和宽度H。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述设置所述平稳压力波形的数据范围的方法包括：
设置所述平稳压力波形的X轴占空比阈值，以及平稳压力波形与所述X轴的2个交叉点之间的插值数；设置采样间隔时间；通过调整所述插值数和所述采样间隔时间，获得所述平稳压力波形的X轴占空比大于等于所述X轴占空比阈值的X轴的数据范围；和/或
设置Y轴占空比阈值为r，设置所述数据范围占所述显示空间的比例为a；则所述Y轴数据长度为L＝(P舒-P收)+(P舒-P收)/r；所述Y轴数据最小值ymin＝P收-(P舒-P收)/ar；所述Y轴数据最大值ymax＝P舒+(P舒-P收)/ar；所述平稳压力的Y轴数据范围为[ymin，ymax]。
可选地，上述的一种获取感兴趣的主动脉压力曲线的方法，所述根据所述显示区间、数据范围以及平稳压力波形生成主动脉压力波形图像的方法包括：
新建长度为W和宽度为H的图片空间；
设置所述图片空间的背景颜色或背景图片、添加边框；
设置(x，y)像素点的颜色；
X轴显示像素范围为[W2,W1-1-W2]，Y轴显示像素范围为[H2,H1-1-H2]；
根据所述(x，y)像素点生成主动脉压力波形图像。
第二方面，本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的获取感兴趣的主动脉压力曲线的方法。
本申请实施例提供的方案带来的有益效果至少包括：
本申请通过对原始主动脉压力波形的计算和判断，最终得到了用户需要的主动脉压力波形图像。本发明可以获取感兴趣的主动脉压力曲线，克服了临床上分析主动脉压力曲线只能依靠观察和经验的弊端。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本申请的获取感兴趣的主动脉压力曲线的方法的流程图；
图2为本申请S200的流程图；
图3为本申请S220的流程图；
图4为本申请S230的流程图；
图5为本申请S232的流程图；
图6为本申请S234的流程图；
图7为本申请S240的流程图；
图8为本申请S250的流程图；
图9为本申请S300的流程图；
图10为本申请S310的流程图；
图11为本申请S320的流程图；
图12为本申请S330的流程图；
图13为本申请S500的流程图；
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下将以图式揭露本发明的多个实施方式，为明确说明起见，许多实务上的细节将在以下叙述中一并说明。然而，应了解到，这些实务上的细节不应用以限制本发明。也就是说，在本发明的部分实施方式中，这些实务上的细节是非必要的。此外，为简化图式起见，一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。
如图1所示，本申请为了解决上述问题，提供了一种本申请提供了一种获取感兴趣的主动脉压力曲线的方法，包括：
本申请提供了一种获取感兴趣的主动脉压力曲线的方法，包括：
S100，获取主动脉压力Pa数据，包括：
通过一次性血压传感器获取主动脉入口压力。
S200，如图2所示，根据主动脉压力Pa提取平稳压力波形，包括：
S210，设置多个索引列表和生理参数阈值；
多个索引列表包括：暂存列表、保存列表、均值点索引列表、极大值点索引列表、极小值点索引列表。
生理参数阈值包括：心跳周期数量C、舒张压至收缩压范围、最小压力幅度阈值、收缩压或舒张压的波动差、心率范围、压力个数N阈、文件数量上限、截止压力阈值P截、截止心率阈值HR截、压力传输速率阈值V阈。
S220，如图3所示，根据索引列表、生理参数阈值，主动脉压力数据计算平均压力，包括：
S221，根据心率阈值获得一个心跳周期时间阈值；
S222，根据一个心跳周期数量C和压力传输速率阈值获取一个心跳周期内的数据点个数范围；
S223，根据数据点个数范围划定一个时间周期，每结束一个时间周期即对所有数据点个数进行统计，且将所有获取的主动脉压力与获取的时间一一对应的存储到暂存列表内；
如果暂存列表内连续n个主动脉压力Pa均小于截止压力阈值P截，则从暂存列表内去除n个主动脉压力Pa，降低计算量，提高运行速度。
S224，如果暂存列表内的数据个数是N阈的整数倍，且大于数据点个数范围的最小阈值，且小于数据点个数范围的最大阈值，则计算所有暂存列表内的主动脉压力的平均值，即为平均压力。
S230，如图4所示，根据平均压力获取压力均值点，包括：
S231，设置均值点索引列表；
S232，如图5所示，向均值点索引列表内加入压力均值点，包括：
S2321，将暂存列表内的主动脉压力Pa分别与平均压力进行比较；
S2322，如果且/>则将PK加入均值点索引列表中，其中PK、PK+1分别表示暂存列表内的第K个、第K+1个主动脉压力；
S2323，如果且/>则将PK加入均值点索引列表中。
S233，根据压力均值点生成压力均值点曲线；
S234，如图6所示，过滤压力均值点曲线上的压力均值点，包括：
S2341，根据Dn＝Pn+1-Pn，其中Pn、Pn+1分别表示均值点索引列表内第n个、第n+1个主动脉压力，Dn表示Pn+1、Pn的压力差；
S2342，根据HRn＝30/(Dn/n)，获取第n个半心跳周期内的心率频率值HRn；
S2343，如果HRn＜HR截，则从均值点索引列表中删除HRn对应的均值压力点；
S2344，如果且在Pn后连续m个主动脉压力均大于/>且/> 则从均值点索引列表中删除Pn后连续m个压力均值点。
S235，判断均值点索引列表内的压力均值点的数量是否满足心跳周期要求，如果不满足则重复上述方法，调整n的大小，直至满足心跳周期要求，包括：
如果均值点列表内的数据个数N均＜2C+1，则不满足心跳周期要求，重复上述方法；
调整n的大小，直至N均≥2C+1。
S240，如图7所示，根据压力均值点分割数据周期，获取每个心跳周期内的压力参数，包括：
S241，如果从暂存列表内提取Pn所在心跳周期内的最大压力Pmax和最大压力Pmax所在的位置Smax，将Smax、Pmax一一对应的均加入舒张压列表中；
S242，获取舒张压列表中的Pmax的平均值，获得舒张压P舒；
如果P舒不在舒张压至收缩压范围内，则重新获取暂存列表，包括：将暂存列表的前半部分截取至保存列表中，如果数据个数大于保存列表数量上限，则将采用保存列表的暂存列表的前半部分数据直接替换到暂存列表内，并清空保存列表。
S243，如果从暂存列表内提取Pn所在心跳周期内的最小压力Pmin和最小压力Pmin所在的位置Smin，将Smin、Pmin一一对应的均加入收缩压列表中；
S244，获取收缩压列表中的Pmin的平均值，获得收缩压P收。
S250，如图8所示，根据压力参数，获取平稳压力波形，包括：
S251，设置偏移量为V阈/a，截取起始点O、结束点E；
S252，将均值点索引列表、舒张压列表、收缩压列表全部按照上述偏移量、起始点O、结束点E执行，输出平稳压力波形、舒张压波形和收缩压波形。
S300，如图9所示，根据平稳压力波形生成主动脉压力波形图像，包括：
S310，如图10所示，设置平稳压力波形的显示区间，包括：
S311，设置生成的主动脉压力波形图像的长度W1和宽度H1；
S312，设置主动脉压力波形图像的外边框长度W2和宽度H2；
S313，根据W＝W1-2W2，H＝H1-2H2，获得平稳压力波形的显示区间的长度W和宽度H。
S320，如图11所示，设置平稳压力波形的数据范围，包括：
S321，设置平稳压力波形的X轴占空比阈值，以及平稳压力波形与X轴的2个交叉点之间的插值数；设置采样间隔时间；通过调整插值数和采样间隔时间，获得平稳压力波形的X轴占空比大于等于X轴占空比阈值的X轴的数据范围；
S322，设置Y轴占空比阈值为r，设置数据范围占显示空间的比例为a；则Y轴数据长度为L＝(P舒-P收)+(P舒-P收)/r；Y轴数据最小值ymin＝P收-(P舒-P收)/ar；Y轴数据最大值ymax＝P舒+(P舒-P收)/ar；平稳压力的Y轴数据范围为[ymin，ymax]。
S330，如图12所示，根据显示区间、数据范围以及平稳压力波形生成主动脉压力波形图像，包括：
S331，新建长度为W和宽度为H的图片空间；
S332，设置图片空间的背景颜色或背景图片、添加边框；
S333，设置(x，y)像素点的颜色；
S334，X轴显示像素范围为[W2,W1-1-W2]，Y轴显示像素范围为[H2,H1-1-H2]；
S335，根据(x，y)像素点生成主动脉压力波形图像。
S400，获取感兴趣的时刻；
S500，如图13所示，根据感兴趣的时刻从主动脉压力波形图像上选取相对应的主动脉压力曲线，包括：
S510，从感兴趣的时刻开始，在主动脉压力波形图像上向前获取若干心跳周期内的压力波形段；
S520，获取若干心跳周期内的压力波形段上的主动脉压力的平均值，为了运算全面和且降低运算量，本申请通常选择6～10个心跳周期内的压力波形段；
S530，根据平均值去除偏差较大的压力值点，根据贝塞尔曲线重新绘制平滑的主动脉压力曲线，包括：
根据平均值去除偏差较大的压力值点，从重新绘制平滑的主动脉压力曲线上保留压力最平稳的一个心跳周期或者两个心跳周期的压力曲线，作为感兴趣的主动脉压力曲线。
本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的获取感兴趣的主动脉压力曲线的方法。
所属技术领域的技术人员知道，本发明的各个方面可以实现为***、方法或计算机程序产品。因此，本发明的各个方面可以具体实现为以下形式，即：完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等)，或硬件和软件方面结合的实施方式，这里可以统称为“电路”、“模块”或“***”。此外，在一些实施例中，本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式，该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或***的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。
例如，可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件，可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作***执行的多个软件指令。在本发明的示例性实施例中，由数据处理器来执行如本文的根据方法和/或***的示例性实施例的一个或多个任务，诸如用于执行多个指令的计算平台。可选地，该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器，例如，磁硬盘和/或可移动介质。可选地，也提供了一种网络连接。可选地也提供显示器和/或用户输入设备，诸如键盘或鼠标。
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的***、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项：
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行***、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行***、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括(但不限于)无线、有线、光缆、RF等等，或者上述的任意合适的组合。
例如，可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码，包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言，诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。
也可以把这些计算机程序指令存储在计算机可读介质中，这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作，从而，存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。
还可将计算机程序指令加载到计算机(例如，冠状动脉分析***)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程，使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。
本发明的以上的具体实例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (16)
1.一种获取感兴趣的主动脉压力曲线的方法，其特征在于，包括：
获取主动脉压力Pa数据；
设置多个索引列表，包括：暂存列表、保存列表、均值点索引列表、极大值点索引列表、极小值点索引列表；
设置生理参数阈值，包括：
心跳周期数量C、舒张压至收缩压范围、最小压力幅度阈值、心率范围，用于收集心率数据；
压力个数N阈，为通过心率和心跳周期C获得的时间内采集到的相应压力个数，其中，采集的心跳周期C大于等于4；
文件数量上限，用于限制收集的压力个数总量；
其中，截止压力阈值P截、截止心率阈值HR截、压力传输速率阈值V阈，用于收集心率压力波形，小于截止压力阈值P截、截止心率阈值HR截、压力传输速率阈值V阈的压力波形为无效压力波形；
根据所述索引列表、所述生理参数阈值，所述主动脉压力Pa数据计算平均压力的方法包括：
根据一个所述心跳周期数量C和压力传输速率阈值V阈获取一个所述心跳周期内的数据点个数范围；
根据所述数据点个数范围划定一个时间周期，每结束一个时间周期即对所有数据点个数进行统计，且将所有获取的主动脉压力与获取的时间一一对应的存储到所述暂存列表内；
如果所述暂存列表内的数据个数是N阈的整数倍，且大于数据点个数范围的最小阈值，且小于数据点个数范围的最大阈值，则计算所有暂存列表内的主动脉压力的平均值，即为平均压力；
根据所述平均压力获取压力均值点；
根据所述压力均值点分割数据周期，获取每个心跳周期内的压力参数；
根据所述压力参数，获取平稳压力波形图像；
获取感兴趣的时刻；
根据所述感兴趣的时刻从所述主动脉压力波形图像上选取相对应的主动脉压力曲线。
2.根据权利要求1所述的一种获取感兴趣的主动脉压力曲线的方法，其特征在于，根据所述感兴趣的时刻从所述主动脉压力波形图像上选取相对应的主动脉压力曲线的方法，包括：
从所述感兴趣的时刻开始，在所述主动脉压力波形图像上向前获取若干心跳周期内的压力波形段；
获取所述若干心跳周期内的压力波形段上的主动脉压力的平均值；
根据所述平均值去除心跳周期偏差最大的压力值点，根据贝塞尔曲线重新绘制平滑的主动脉压力曲线。
3.根据权利要求2所述的一种获取感兴趣的主动脉压力曲线的方法，其特征在于，所述若干心跳周期内的压力波形段的为6～10个心跳周期内的压力波形段。
4.根据权利要求2所述的一种获取感兴趣的主动脉压力曲线的方法，其特征在于，所述根据所述平均值去除心跳周期偏差最大的压力值点，根据所述贝塞尔曲线重新绘制平滑的主动脉压力曲线的方法包括：
根据所述平均值去除偏差最大的压力值点，从所述重新绘制平滑的主动脉压力曲线上保留压力最平稳的一个心跳周期或者两个心跳周期的压力曲线，作为感兴趣的主动脉压力曲线。
5.根据权利要求1所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述根据所述索引列表、所述生理参数阈值，所述主动脉压力数据计算平均压力的方法还包括：
如果所述暂存列表内连续n个主动脉压力Pa均小于截止压力阈值P截，则从所述暂存列表内去除所述n个主动脉压力Pa。
6.根据权利要求5所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述根据所述平均压力获取压力均值点的方法，包括：
设置均值点索引列表；
向所述均值点索引列表内加入压力均值点；
根据压力均值点生成压力均值点曲线；
过滤所述压力均值点曲线上的压力均值点；
判断所述均值点索引列表内的所述压力均值点的数量是否满足心跳周期要求，如果不满足则调整n的大小，直至满足心跳周期要求。
7.根据权利要求6所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述向所述均值点索引列表内加入压力均值点的方法包括：
将所述暂存列表内的主动脉压力Pa分别与所述平均压力P进行比较；
如果PK＞P且PK+1＜P，则将PK加入所述均值点索引列表中，其中PK、PK+1分别表示所述暂存列表内的第K个、第K+1个主动脉压力；
如果PK＜P且PK+1＞P，则将PK加入所述均值点索引列表中。
8.根据权利要求7所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述过滤所述压力均值点曲线上的压力均值点的方法包括：
根据Dn＝Pn+1-Pn，其中Pn、Pn+1分别表示所述均值点索引列表内第n个、第n+1个主动脉压力，所述Dn表示Pn+1、Pn的压力差；
根据HRn＝30/(Dn/n)，获取所述第n个半心跳周期内的心率频率值HRn；
如果HRn＜HR截，则从所述均值点索引列表中删除所述HRn对应的均值压力点；
如果Pn＞P，且在Pn后连续m个主动脉压力均大于P，且Pn+m+1＜P，则从所述均值点索引列表中删除Pn后连续m个压力均值点。
9.根据权利要求8所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述判断所述均值点索引列表内的所述压力均值点的数量是否满足心跳周期要求，如果不满足则重复所述权利要求1，调整n的大小，直至满足心跳周期要求的方法包括：
如果均值点列表内的数据个数N均＜2C+1，则不满足心跳周期要求，调整n的大小，直至N均≥2C+1。
10.根据权利要求9所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述根据所述压力均值点分割数据周期，获取每个心跳周期内的压力参数的方法，包括：
如果Pn＞P，从所述暂存列表内提取Pn所在心跳周期内的最大压力Pmax和所述最大压力Pmax所在的位置Smax，将所述Smax、Pmax一一对应的均加入舒张压列表中；
获取所述舒张压列表中的Pmax的平均值，获得舒张压P舒；
如果Pn＜P，从所述暂存列表内提取Pn所在心跳周期内的最小压力Pmin和所述最小压力Pmin所在的位置Smin，将所述Smin、Pmin一一对应的均加入收缩压列表中；
获取所述收缩压列表中的Pmin的平均值，获得收缩压P收。
11.根据权利要求10所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述根据所述压力参数，获取平稳压力波形的方法，包括：
设置偏移量为V阈/a，其中a＝4或者a＝5，在V阈波形范围内的起点和终点偏移，截取起始点O、结束点E；
将所述均值点索引列表、舒张压列表、收缩压列表全部按照上述偏移量、起始点O、结束点E执行，输出所述平稳压力波形、舒张压波形和收缩压波形。
12.根据权利要求11所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述生成主动脉压力波形图像的方法包括：
设置所述平稳压力波形的显示区间；
设置所述平稳压力波形的数据范围；
根据所述显示区间、数据范围以及平稳压力波形生成主动脉压力波形图像。
13.根据权利要求12所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述设置所述平稳压力波形的显示区间的方法包括：
设置所述生成的主动脉压力波形图像的长度W1和宽度H1；
设置所述主动脉压力波形图像的外边框长度W2和宽度H2；
根据W＝W1-2W2，H＝H1-2H2，获得所述平稳压力波形的显示区间的长度W和宽度H。
14.根据权利要求13所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述设置所述平稳压力波形的数据范围的方法包括：
设置所述平稳压力波形的X轴占空比阈值，以及平稳压力波形与所述X轴的2个交叉点之间的插值数；设置采样间隔时间；通过调整所述插值数和所述采样间隔时间，获得所述平稳压力波形的X轴占空比大于等于所述X轴占空比阈值的X轴的数据范围；和/或
设置Y轴占空比阈值为r，设置所述数据范围占显示空间的比例为a；则所述Y轴数据长度为L＝(P舒-P收)+(P舒-P收)/r；所述Y轴数据最小值ymin＝P收-(P舒-P收)/ar；所述Y轴数据最大值ymax＝P舒+(P舒-P收)/ar；所述平稳压力的Y轴数据范围为[ymin，ymax]。
15.根据权利要求14所述的获取感兴趣的主动脉压力曲线的方法，其特征在于，所述根据所述显示区间、数据范围以及平稳压力波形生成主动脉压力波形图像的方法包括：
新建长度为W和宽度为H的图片空间；
设置所述图片空间的背景颜色或背景图片、添加边框；
设置(x，y)像素点的颜色；
X轴显示像素范围为[W2,W1-1-W2]，Y轴显示像素范围为[H2,H1-1-H2]；
根据所述(x，y)像素点生成主动脉压力波形图像。
16.一种计算机存储介质，其特征在于，计算机程序被处理器执行时实现权利要求1～15任一项所述的获取感兴趣的主动脉压力曲线的方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011348579.2A CN112494020B (zh) | 2020-11-26 | 2020-11-26 | 获取感兴趣的主动脉压力曲线的方法及存储介质 |
| PCT/CN2020/132275 WO2022110020A1 (zh) | 2020-11-26 | 2020-11-27 | 获取感兴趣的主动脉压力曲线的方法及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011348579.2A CN112494020B (zh) | 2020-11-26 | 2020-11-26 | 获取感兴趣的主动脉压力曲线的方法及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112494020A CN112494020A (zh) | 2021-03-16 |
| CN112494020B true CN112494020B (zh) | 2024-01-23 |
Family
ID=74966271
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011348579.2A Active CN112494020B (zh) | 2020-11-26 | 2020-11-26 | 获取感兴趣的主动脉压力曲线的方法及存储介质 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112494020B (zh) |
| WO (1) | WO2022110020A1 (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112494019A (zh) * | 2020-11-26 | 2021-03-16 | 苏州润迈德医疗科技有限公司 | 基于主动脉压力提取平稳压力波形的方法及存储介质 |
| CN112494021A (zh) * | 2020-11-26 | 2021-03-16 | 苏州润迈德医疗科技有限公司 | 主动脉压力波形图像的生成方法及存储介质 |
| CN115271256B (zh) * | 2022-09-20 | 2022-12-16 | 华东交通大学 | 一种多维分类下的智能化订货方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109475309A (zh) * | 2016-07-13 | 2019-03-15 | 悦享趋势科技（北京）有限责任公司 | 检测生理状态的方法和装置 |
| CN109770888A (zh) * | 2019-03-19 | 2019-05-21 | 苏州润迈德医疗科技有限公司 | 基于压力传感器和造影图像计算瞬时无波形比率的方法 |
| CN109805949A (zh) * | 2019-03-19 | 2019-05-28 | 苏州润迈德医疗科技有限公司 | 基于压力传感器和造影图像计算血流储备分数的方法 |
| CN110226923A (zh) * | 2018-03-05 | 2019-09-13 | 苏州润迈德医疗科技有限公司 | 一种无需血管扩张剂测量血流储备分数的方法 |
| CN110353645A (zh) * | 2019-08-19 | 2019-10-22 | 苏州润迈德医疗科技有限公司 | 测量无波形期压力、比率方法、装置、***及存储介质 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9314170B2 (en) * | 2010-05-07 | 2016-04-19 | Atcor Medical Pty Ltd | Brachial cuff |
| CN103284703B (zh) * | 2013-06-18 | 2015-05-13 | 中国科学院合肥物质科学研究院 | 基于上肢动脉信息的主动脉脉搏波传递时间测量方法 |
| CN104983412B (zh) * | 2015-05-28 | 2017-08-11 | 中国科学院合肥物质科学研究院 | 一种中心动脉收缩期平均归一化血流波形模型及基于其获取主动脉脉搏波传递时间的方法 |
| CN107992452B (zh) * | 2017-12-12 | 2020-12-18 | 北京动亮健康科技有限公司 | 计算中心血流动力学指标的方法、装置、存储介质及设备 |
| CN110786840B (zh) * | 2019-11-04 | 2021-06-08 | 苏州润迈德医疗科技有限公司 | 基于生理参数获取血管评定参数的方法、装置及存储介质 |
-
2020
- 2020-11-26 CN CN202011348579.2A patent/CN112494020B/zh active Active
- 2020-11-27 WO PCT/CN2020/132275 patent/WO2022110020A1/zh active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109475309A (zh) * | 2016-07-13 | 2019-03-15 | 悦享趋势科技（北京）有限责任公司 | 检测生理状态的方法和装置 |
| CN110226923A (zh) * | 2018-03-05 | 2019-09-13 | 苏州润迈德医疗科技有限公司 | 一种无需血管扩张剂测量血流储备分数的方法 |
| CN109770888A (zh) * | 2019-03-19 | 2019-05-21 | 苏州润迈德医疗科技有限公司 | 基于压力传感器和造影图像计算瞬时无波形比率的方法 |
| CN109805949A (zh) * | 2019-03-19 | 2019-05-28 | 苏州润迈德医疗科技有限公司 | 基于压力传感器和造影图像计算血流储备分数的方法 |
| CN110353645A (zh) * | 2019-08-19 | 2019-10-22 | 苏州润迈德医疗科技有限公司 | 测量无波形期压力、比率方法、装置、***及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022110020A1 (zh) | 2022-06-02 |
| CN112494020A (zh) | 2021-03-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112494020B (zh) | 获取感兴趣的主动脉压力曲线的方法及存储介质 | |
| CN108186038B (zh) | 基于动脉造影影像计算冠脉血流储备分数的*** | |
| CN108511075B (zh) | 一种非侵入式获取血流储备分数的方法和*** | |
| CN108992057B (zh) | 一种确定冠状动脉血流储备分数ffr的方法和装置 | |
| JP6685319B2 (ja) | 定量的フロー分析のための方法および装置 | |
| US20210236000A1 (en) | Method, device and system for acquiring blood vessel evaluation parameters based on angiographic image | |
| CN112164020A (zh) | 精确提取血管中心线的方法、装置、分析***和存储介质 | |
| CN110786841B (zh) | 基于微循环阻力指数调节最大充血状态流速的方法及装置 | |
| CN110786842B (zh) | 测量舒张期血流速度的方法、装置、***及存储介质 | |
| CN110786840B (zh) | 基于生理参数获取血管评定参数的方法、装置及存储介质 | |
| CN112419276B (zh) | 调节血管轮廓及中心线的方法及存储介质 | |
| EP3606433B1 (en) | Standardized coronary artery disease metric | |
| WO2022110095A1 (zh) | 获取冠状动脉血管评价参数的方法及存储介质 | |
| CN112132882A (zh) | 从冠状动脉二维造影图像中提取血管中心线的方法和装置 | |
| CN112472112B (zh) | 调节血管狭窄区间的方法、***及存储介质 | |
| WO2022110019A1 (zh) | 主动脉压力波形图像的生成方法及存储介质 | |
| WO2020083390A1 (zh) | 获取心表大动脉的血流量的方法、装置、***及存储介质 | |
| CN110929604B (zh) | 基于造影图像的流速的筛选方法、装置、***和存储介质 | |
| CN111815584B (zh) | 基于ct序列图像获取心脏重心的方法和*** | |
| WO2022110021A1 (zh) | 基于主动脉压力提取平稳压力波形的方法及存储介质 | |
| CN113017667A (zh) | 一种血管狭窄的量化方法、装置、设备和可读存储介质 | |
| WO2022109909A1 (zh) | 精确获取狭窄病变区间的方法、***及存储介质 | |
| JP2015217113A (ja) | 血管解析装置、医用画像診断装置、血管解析方法及び血管解析プログラム | |
| EP4005472B1 (en) | Method and apparatus for correcting blood flow velocity on the basis of interval time between angiographic images | |
| CN117809852A (zh) | 一种冠脉血流储备分数的确定方法、装置、设备及介质 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |