CN112472112B - 调节血管狭窄区间的方法、系统及存储介质 - Google Patents
调节血管狭窄区间的方法、系统及存储介质 Download PDFInfo
- Publication number
- CN112472112B CN112472112B CN202011337599.XA CN202011337599A CN112472112B CN 112472112 B CN112472112 B CN 112472112B CN 202011337599 A CN202011337599 A CN 202011337599A CN 112472112 B CN112472112 B CN 112472112B
- Authority
- CN
- China
- Prior art keywords
- stenosis
- pipe diameter
- point
- curve
- interval
- 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 49
- 206010057469 Vascular stenosis Diseases 0.000 title claims abstract description 48
- 230000001105 regulatory effect Effects 0.000 title claims abstract description 6
- 208000031481 Pathologic Constriction Diseases 0.000 claims abstract description 117
- 230000036262 stenosis Effects 0.000 claims abstract description 117
- 208000037804 stenosis Diseases 0.000 claims abstract description 117
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 48
- 230000003902 lesion Effects 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 9
- 230000036285 pathological change Effects 0.000 claims description 9
- 231100000915 pathological change Toxicity 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000002586 coronary angiography Methods 0.000 claims description 6
- 230000002966 stenotic effect Effects 0.000 claims description 5
- 230000001575 pathological effect Effects 0.000 claims description 4
- 230000003993 interaction Effects 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 11
- 210000004351 coronary vessel Anatomy 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 238000002591 computed tomography 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
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000002792 vascular 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
- 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
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 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
- 238000003384 imaging method Methods 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
- 230000004089 microcirculation 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
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Computing Systems (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Pulmonology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请提供了一种调节血管狭窄区间的方法、系统及存储介质，包括：获取血管的修正管径曲线；获取血管狭窄起始点、结束点以及血管狭窄区间；在所述起始点、所述结束点上分别设置标记；在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。本申请能有效避免对血管狭窄区间的直接定义失真，通过滑动杆将人因工程的概念融入调节血管狭窄区间的操作和设计，人工交互与自动化判断的结合在一起，有效提高了准确性，降低单纯使用自动化判断算法的风险，提高整个系统的鲁棒性。
Description
技术领域
本发明涉及冠状动脉医学技术领域，特别是涉及一种调节血管狭窄区间的方法、系统及存储介质。
背景技术
人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块，继而导致血管狭窄；特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足，诱发冠心病、心绞痛等病症，对人类的健康造成严重威胁。据统计，我国现有冠心病患者约1100万人，心血管介入手术治疗患者数量每年增长大于10％。
冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度，但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性，1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve，FFR)，经过长期的基础与临床研究，FFR已成为冠脉狭窄功能性评价的金标准。
FFR是冠状动脉血管评定参数的一种，微循环阻力指数IMR等属于冠状动脉血管评定参数。
现有技术总是试图通过算法自动判断血管的狭窄区间，但由于成像设备带来的误差，以及患者血管血管情况的多样性等多种原有，单纯通过像素信息无法每次都得到正确的狭窄区间，造成使用者在通过算法自动判断血管的狭窄区间的过程中，必须对结果结合自己的经验进行再次确认。
发明内容
本发明提供了一种调节血管狭窄区间的方法、系统及存储介质，以解决单纯通过算法无法自动精准判断血管的狭窄区间，给使用者提供一款合适的人工交互工具，对结果进行调整并确认。
为实现上述目的，第一方面，本申请提供了一种调节血管狭窄区间的方法，包括：
获取血管的修正管径曲线；
获取血管狭窄起始点、结束点以及血管狭窄区间；
在所述起始点、所述结束点上分别设置标记；
在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；
如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。
可选地，上述的调节血管狭窄区间的方法，所述获取血管的修正管径曲线的方法包括：
从所述冠脉造影图像上获取真实管径曲线；
拟合正常血管管径，获取拟合管径曲线；
根据所述拟合管径曲线和真实管径曲线，获取狭窄病变血管段；
从所述狭窄病变血管段获取的管径曲线为修正管径曲线。
可选地，上述的调节血管狭窄区间的方法，所述根据所述拟合管径曲线和真实管径曲线，获取狭窄病变血管段的方法包括：
根据所述拟合管径曲线和真实管径曲线，获取第一狭窄病变区间；
从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间；
根据所述第二狭窄病变区间重新拟合管径区间曲线，获取狭窄点；
根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间、狭窄起始点和结束点；
所述第三狭窄病变区间内获取的管径曲线为修正管径曲线。
可选地，上述的调节血管狭窄区间的方法，所述拟合正常血管管径，获取拟合管径曲线的方法包括：
根据拟合代价函数，获取拟合管径，具体公式为
其中，i表示第i个管径的曲线采样点；n表示管径曲线采样数量的总和；xi表示第i个管径的曲线采样点的长度；yi表示在xi处的管径；
将每个所述拟合管径对应到坐标系中获得相应的管径点，将所述管径点依次平滑连接，获取所述拟合管径曲线。
可选地，上述的调节血管狭窄区间的方法，所述根据所述拟合管径曲线和真实管径曲线，获取第一狭窄区间的方法包括：
获取所述血管的真实管径；
将所述真实管径对应到拟合管径曲线的坐标系中；
获取所述真实管径曲线，以及所述真实管径曲线与所述拟合管径曲线的交叉点；
如果交叉点前一点的所述真实管径大于所述拟合管径，则所述交叉点为所述狭窄区域的第一入口点，反之，所述交叉点为所述狭窄区域的第一出口点；
所述第一入口点和所述第一出口点之间的曲线为初步判定的狭窄位置，即所述第一狭窄病变区间。
可选地，上述的调节血管狭窄区间的方法，所述从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间的方法包括；
计算狭窄度；
计算所述第一狭窄病变区间的血管中心线长度L；
根据所述狭窄度和/或血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间。
可选地，上述的调节血管狭窄区间的方法，所述计算狭窄度的方法包括：
其中，A表示血管的狭窄度，Dmin表示第一入口点与第一出口点之间的血管的最小管径，D入、D出分别表示第一入口点的血管管径和第一出口点的血管管径。
可选地，上述的调节血管狭窄区间的方法，所述根据所述狭窄度和/或血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间的方法包括：
如果A＜0.2，则判断为误判区域，采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线；
如果L＜5mm，则判断为误判区域，采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线；
去除误判区域后重新获得的狭窄区域，即为所述第二狭窄病变区间。
可选地，上述的调节血管狭窄区间的方法，所述根据所述第二狭窄病变区间重新拟合管径区间曲线，获取狭窄点的方法包括：
取所述第一入口点前1～3cm，所述第一出口点后1～3cm之间的区域，根据所述拟合代价函数重新拟合管径区间曲线；
获取所述重新拟合管径区间曲线中管径最小的点为狭窄点。
可选地，上述的调节血管狭窄区间的方法，所述根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间的方法包括：
在所述狭窄点两侧，所述拟合管径区间曲线与所述真实管径曲线交叉的两个点为第二入口点和第二出口点，所述第二入口点和所述第二出口点之间的区间为第三狭窄病变区间，所述第二入口点为狭窄起始点，所述第二出口点为狭窄出口点。
第二方面，本申请提供了一种调节血管狭窄区间的系统，包括：管径曲线获取装置、血管狭窄区间获取装置、标记装置和滑动装置；
所述管径曲线获取装置，用于获取获取血管的修正管径曲线；
所述血管狭窄区间获取装置，与所述管径曲线获取装置连接，用于获取血管狭窄起始点、结束点以及血管狭窄区间；
所述标记装置，与所述血管狭窄区间获取装置来接，用于在所述起始点、所述结束点上分别设置标记；
所述滑动装置，与所述管径曲线获取装置、所述血管狭窄区间获取装置连接、所述标记装置连接，用于在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。
可选地，上述的调节血管狭窄区间的系统，所述血管狭窄区间获取装置包括第一狭窄病变单元、第二狭窄病变单元、狭窄点单元和第三狭窄病变单元；
所述第一狭窄病变单元，与所述管径曲线获取装置连接，用于根据所述拟合管径曲线和真实管径曲线，获取第一狭窄病变区间；
所述第二狭窄病变单元，与所述第一狭窄病变单元连接，用于从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间；
所述狭窄点单元，与所述管径曲线获取装置连接，用于根据拟合管径区间曲线获取狭窄点；
所述第三狭窄病变单元，与所述狭窄点单元、所述管径曲线获取装置连接，用于根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间。
第三方面，本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的调节血管狭窄区间的方法。
本申请实施例提供的方案带来的有益效果至少包括：
本申请提供了一种调节血管狭窄区间的方法，通过滑动杆实现人工交互与自动化算法相配合的方法，将人因工程的概念融入调节血管狭窄区间的操作和设计，给使用者能够快速确认血管狭窄区间，降低操作难度，提高工作效率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本申请的调节血管狭窄区间的方法的流程图；
图2为本申请的S100的流程图；
图3为本申请的S120的流程图；
图4为本申请的S140的流程图；
图5为本申请的S141的流程图；
图6为本申请的S142的流程图；
图7为本申请的S143的流程图；
图8为本申请的调节血管狭窄区间的系统的结构框图；
图9为本申请的血管狭窄区间获取装置200的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下将以图式揭露本发明的多个实施方式，为明确说明起见，许多实务上的细节将在以下叙述中一并说明。然而，应了解到，这些实务上的细节不应用以限制本发明。也就是说，在本发明的部分实施方式中，这些实务上的细节是非必要的。此外，为简化图式起见，一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。
在冠状动脉造影图像中，需要结合流体力学分析计算冠状动脉血管评定参数，而现有技术中没有用于流体力学分析的血管三维网格模型。
实施例1：
如图1所示，本申请提供了一种调节血管狭窄区间的方法，包括：
S100，如图2所示，获取血管的修正管径曲线，包括：
S110，从所述冠脉造影图像上获取真实管径曲线；
S120，如图3所示，拟合正常血管管径，获取拟合管径曲线，包括：
S121，根据拟合代价函数，获取拟合管径，具体公式为：
其中，i表示第i个管径的曲线采样点；n表示管径曲线采样数量的总和；xi表示第i个管径的曲线采样点的长度；yi表示在xi处的管径；
S122，将每个所述拟合管径对应到坐标系中获得相应的管径点，将所述管径点依次平滑连接，获取所述拟合管径曲线。
S130，根据所述拟合管径曲线和真实管径曲线，获取狭窄病变血管段；
S140，如图4所示，从所述狭窄病变血管段获取的管径曲线为修正管径曲线，包括：
S141，如图5所示，根据所述拟合管径曲线和真实管径曲线，获取第一狭窄病变区间，包括：
S1411，获取所述血管的真实管径；
S1412，将所述真实管径对应到拟合管径曲线的坐标系中；
S1413，获取所述真实管径曲线，以及所述真实管径曲线与所述拟合管径曲线的交叉点；
S1414，如果交叉点前一点的所述真实管径大于所述拟合管径，则所述交叉点为所述狭窄区域的第一入口点，反之，所述交叉点为所述狭窄区域的第一出口点；
S1415，所述第一入口点和所述第一出口点之间的曲线为初步判定的狭窄位置，即所述第一狭窄病变区间。
S142，如图6所示，从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间，包括：
S1421，计算狭窄度，具体公式：其中，A表示血管的狭窄度，Dmin表示第一入口点与第一出口点之间的血管的最小管径，D入、D出分别表示第一入口点的血管管径和第一出口点的血管管径。
S1422，计算所述第一狭窄病变区间的血管中心线长度L；
S1423，根据所述狭窄度和/或血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间，包括：
I)如果A＜0.2，则判断为误判区域，采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线；
II)如果L＜5mm，则判断为误判区域，采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线；
III)去除误判区域后重新获得的狭窄区域，即为所述第二狭窄病变区间。
S143，如图7所示，根据所述第二狭窄病变区间重新拟合管径区间曲线，获取狭窄点，包括：
S1431，取所述第一入口点前1～3cm，所述第一出口点后1～3cm之间的区域，根据所述拟合代价函数重新拟合管径区间曲线；
S1432，获取所述重新拟合管径区间曲线中管径最小的点为狭窄点。
S144，根据所述狭窄点获取第三狭窄病变区间即为精确狭窄病变区间、狭窄起始点和结束点，包括：在所述狭窄点两侧，所述拟合管径区间曲线与所述真实管径曲线交叉的两个点为第二入口点和第二出口点，所述第二入口点和所述第二出口点之间的区间为第三狭窄病变区间，所述第二入口点为狭窄起始点，所述第二出口点为狭窄出口点。
S145，所述第三狭窄病变区间内获取的管径曲线为修正管径曲线。
S200，获取血管狭窄起始点、结束点以及血管狭窄区间；
S300，在所述起始点、所述结束点上分别设置标记；
S400，在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；
S500，如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。
如图8所示，本申请提供了一种调节血管狭窄区间的系统，包括：管径曲线获取装置100、血管狭窄区间获取装置200、标记装置300和滑动装置400；所述管径曲线获取装置100，用于获取获取血管的修正管径曲线；所述血管狭窄区间获取装置200，与所述管径曲线获取装置100连接，用于获取血管狭窄起始点、结束点以及血管狭窄区间；所述标记装置300，与所述血管狭窄区间获取装置200来接，用于在所述起始点、所述结束点上分别设置标记；所述滑动装置400，与所述管径曲线获取装置100、所述血管狭窄区间获取装置200连接、所述标记装置300连接，用于在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。
如图9所示，本申请的一个实施例中，血管狭窄区间获取装置200包括第一狭窄病变单元210、第二狭窄病变单元220、狭窄点单元230和第三狭窄病变单元240；所述第一狭窄病变单元210，与所述管径曲线获取装置100连接，用于根据所述拟合管径曲线和真实管径曲线，获取第一狭窄病变区间；所述第二狭窄病变单元220，与所述第一狭窄病变单元210连接，用于从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间；所述狭窄点单元230，与所述管径曲线获取装置连100接，用于根据拟合管径区间曲线获取狭窄点；所述第三狭窄病变单元240，与所述狭窄点单元230、所述管径曲线获取装置100连接，用于根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间。
本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的调节血管狭窄区间的方法。
所属技术领域的技术人员知道，本发明的各个方面可以实现为系统、方法或计算机程序产品。因此，本发明的各个方面可以具体实现为以下形式，即：完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等)，或硬件和软件方面结合的实施方式，这里可以统称为“电路”、“模块”或“系统”。此外，在一些实施例中，本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式，该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或系统的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。
例如，可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件，可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作系统执行的多个软件指令。在本发明的示例性实施例中，由数据处理器来执行如本文的根据方法和/或系统的示例性实施例的一个或多个任务，诸如用于执行多个指令的计算平台。可选地，该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器，例如，磁硬盘和/或可移动介质。可选地，也提供了一种网络连接。可选地也提供显示器和/或用户输入设备，诸如键盘或鼠标。
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项：
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括(但不限于)无线、有线、光缆、RF等等，或者上述的任意合适的组合。
例如，可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码，包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言，诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。
也可以把这些计算机程序指令存储在计算机可读介质中，这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作，从而，存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。
还可将计算机程序指令加载到计算机(例如，冠状动脉分析系统)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程，使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。
本发明的以上的具体实例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (11)
1.一种调节血管狭窄区间的方法，其特征在于，包括：
从冠脉造影图像上获取真实管径曲线；
拟合正常血管管径，获取拟合管径曲线；
根据所述拟合管径曲线和真实管径曲线，获取第一狭窄病变区间；
从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间；
根据所述第二狭窄病变区间重新拟合管径区间曲线，获取狭窄点；
根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间、狭窄起始点和结束点；
所述第三狭窄病变区间内获取的管径曲线为修正管径曲线；
获取血管的修正管径曲线；
获取血管狭窄起始点、结束点以及血管狭窄区间；
在所述起始点、所述结束点上分别设置标记；
在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；
如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。
2.根据权利要求1所述的调节血管狭窄区间的方法，其特征在于，所述拟合正常血管管径，获取拟合管径曲线的方法包括：
根据拟合代价函数，获取拟合管径，具体公式为
其中，i表示第i个管径的曲线采样点；n表示管径曲线采样数量的总和；xi表示第i个管径的曲线采样点的长度；yi表示在xi处的管径；
将每个所述拟合管径对应到坐标系中获得相应的管径点，将所述管径点依次平滑连接，获取所述拟合管径曲线。
3.根据权利要求2所述的调节血管狭窄区间的方法，其特征在于，所述根据所述拟合管径曲线和真实管径曲线，获取第一狭窄区间的方法包括：
获取所述血管的真实管径；
将所述真实管径对应到拟合管径曲线的坐标系中；
获取所述真实管径曲线，以及所述真实管径曲线与所述拟合管径曲线的交叉点；
如果交叉点前一点的所述真实管径大于所述拟合管径，则所述交叉点为所述狭窄区域的第一入口点，反之，所述交叉点为所述狭窄区域的第一出口点；
所述第一入口点和所述第一出口点之间的曲线为初步判定的狭窄位置，即所述第一狭窄病变区间。
4.根据权利要求3所述的调节血管狭窄区间的方法，其特征在于，所述从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间的方法包括；
计算狭窄度；
计算所述第一狭窄病变区间的血管中心线长度L；
根据所述狭窄度和/或血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间。
5.根据权利要求4所述的调节血管狭窄区间的方法，其特征在于，所述计算狭窄度的方法包括：
其中，A表示血管的狭窄度，Dmin表示第一入口点与第一出口点之间的血管的最小管径，D入、D出分别表示第一入口点的血管管径和第一出口点的血管管径。
6.根据权利要求5所述的调节血管狭窄区间的方法，其特征在于，所述根据所述狭窄度和/或血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间的方法包括：
如果A＜0.2，则判断为误判区域，采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线；
如果L＜5mm，则判断为误判区域，采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线；
去除误判区域后重新获得的狭窄区域，即为所述第二狭窄病变区间。
7.根据权利要求6所述的调节血管狭窄区间的方法，其特征在于，所述根据所述第二狭窄病变区间重新拟合管径区间曲线，获取狭窄点的方法包括：
取所述第一入口点前1～3cm，所述第一出口点后1～3cm之间的区域，根据所述拟合代价函数重新拟合管径区间曲线；
获取所述重新拟合管径区间曲线中管径最小的点为狭窄点。
8.根据权利要求7所述的调节血管狭窄区间的方法，其特征在于，所述根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间的方法包括：
在所述狭窄点两侧，所述拟合管径区间曲线与所述真实管径曲线交叉的两个点为第二入口点和第二出口点，所述第二入口点和所述第二出口点之间的区间为第三狭窄病变区间，所述第二入口点为狭窄起始点，所述第二出口点为狭窄出口点。
9.一种调节血管狭窄区间的系统，其特征在于，包括：管径曲线获取装置、血管狭窄区间获取装置、标记装置和滑动装置；
所述管径曲线获取装置，用于获取获取血管的修正管径曲线；
所述血管狭窄区间获取装置，与所述管径曲线获取装置连接，用于获取血管狭窄起始点、结束点以及血管狭窄区间；
所述标记装置，与所述血管狭窄区间获取装置来接，用于在所述起始点、所述结束点上分别设置标记；
所述滑动装置，与所述管径曲线获取装置、所述血管狭窄区间获取装置连接、所述标记装置连接，用于在所述修正管径曲线上设置滑动杆，所述滑动杆一端对应所述起始点，所述滑动杆一端对应所述结束点；如果滑动滑动杆，则所述起始点和所述结束点分别对应的标记在所述修正管径曲线上均发生相应距离的滑动，即调节了所述血管狭窄区间。
10.根据权利要求9所述的调节血管狭窄区间的系统，其特征在于，所述血管狭窄区间获取装置包括第一狭窄病变单元、第二狭窄病变单元、狭窄点单元和第三狭窄病变单元；
所述第一狭窄病变单元，与所述管径曲线获取装置连接，用于根据拟合管径曲线和真实管径曲线，获取第一狭窄病变区间；
所述第二狭窄病变单元，与所述第一狭窄病变单元连接，用于从所述第一狭窄病变区间去除误判的狭窄区域，获得第二狭窄病变区间；
所述狭窄点单元，与所述管径曲线获取装置连接，用于根据拟合管径区间曲线获取狭窄点；
所述第三狭窄病变单元，与所述狭窄点单元、所述管径曲线获取装置连接，用于根据所述狭窄点获取第三狭窄病变区间，即为精确狭窄病变区间。
11.一种计算机存储介质，其特征在于，计算机程序被处理器执行时实现权利要求1～8任一项所述的调节血管狭窄区间的方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011337599.XA CN112472112B (zh) | 2020-11-25 | 2020-11-25 | 调节血管狭窄区间的方法、系统及存储介质 |
| PCT/CN2020/131704 WO2022109908A1 (zh) | 2020-11-25 | 2020-11-26 | 调节血管狭窄区间的方法、系统及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011337599.XA CN112472112B (zh) | 2020-11-25 | 2020-11-25 | 调节血管狭窄区间的方法、系统及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112472112A CN112472112A (zh) | 2021-03-12 |
| CN112472112B true CN112472112B (zh) | 2024-02-27 |
Family
ID=74934529
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011337599.XA Active CN112472112B (zh) | 2020-11-25 | 2020-11-25 | 调节血管狭窄区间的方法、系统及存储介质 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112472112B (zh) |
| WO (1) | WO2022109908A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114359128A (zh) * | 2021-09-10 | 2022-04-15 | 数坤(北京)网络科技股份有限公司 | 一种血管狭窄的检测方法、装置及计算机可读介质 |
| CN117179894B (zh) * | 2023-11-07 | 2024-03-15 | 北京唯迈医疗设备有限公司 | 介入手术支架辅助选型和定位的装置、系统和存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105380637A (zh) * | 2014-09-02 | 2016-03-09 | 株式会社东芝 | 医用图像处理装置以及医用图像处理方法 |
| CN106539622A (zh) * | 2017-01-28 | 2017-03-29 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入方法和系统 |
| WO2020190976A2 (en) * | 2019-03-17 | 2020-09-24 | Lightlab Imaging, Inc. | Arterial imaging and assessment systems and methods and related user interface based-workflows |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101991420B (zh) * | 2009-08-27 | 2013-01-16 | 上海西门子医疗器械有限公司 | 一种疑似血管狭窄区域的定位方法和系统 |
| WO2018060529A1 (en) * | 2016-09-30 | 2018-04-05 | Koninklijke Philips N.V. | Apparatus for determining a functional index for stenosis assessment |
| CN108932714B (zh) * | 2018-07-23 | 2021-11-23 | 苏州润迈德医疗科技有限公司 | 冠脉ct图像的斑块分类方法 |
-
2020
- 2020-11-25 CN CN202011337599.XA patent/CN112472112B/zh active Active
- 2020-11-26 WO PCT/CN2020/131704 patent/WO2022109908A1/zh active Application Filing
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105380637A (zh) * | 2014-09-02 | 2016-03-09 | 株式会社东芝 | 医用图像处理装置以及医用图像处理方法 |
| CN106539622A (zh) * | 2017-01-28 | 2017-03-29 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入方法和系统 |
| WO2020190976A2 (en) * | 2019-03-17 | 2020-09-24 | Lightlab Imaging, Inc. | Arterial imaging and assessment systems and methods and related user interface based-workflows |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022109908A1 (zh) | 2022-06-02 |
| CN112472112A (zh) | 2021-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US12082912B2 (en) | Lumen morphology and vascular resistance measurements data collection systems apparatus and methods | |
| CN110889896B (zh) | 获取血管狭窄病变区间及三维合成方法、装置和系统 | |
| JP7203077B2 (ja) | 血管内画像化システムの作動方法 | |
| CN112164020B (zh) | 精确提取血管中心线的方法、装置、分析系统和存储介质 | |
| US10431328B1 (en) | Method and system for anatomical tree structure analysis | |
| CN112472112B (zh) | 调节血管狭窄区间的方法、系统及存储介质 | |
| EP3125764A1 (en) | Processing apparatus and method for processing cardiac data of a living being | |
| CN112132949B (zh) | 具有狭窄病变区间的血管三维建模方法、装置和系统 | |
| CN112419276B (zh) | 调节血管轮廓及中心线的方法及存储介质 | |
| US20200188027A1 (en) | Virtual stent placement apparatus, virtual stent placement method, and virtual stent placement program | |
| CN113040812B (zh) | 医学图像分析方法、装置、计算机设备和存储介质 | |
| JP2020510451A (ja) | 生物物理的シミュレーション用の血管樹標準化及び／又はプルーニングされた部分用の拡張シミュレーション | |
| CN114677492A (zh) | 冠状动脉微循环血流动力学仿真分析方法及装置 | |
| CN112419280B (zh) | 精确获取狭窄病变区间的方法及存储介质 | |
| US20230172575A1 (en) | Method, system and computer program product for determining ischemia region of the organ | |
| CN111815584B (zh) | 基于ct序列图像获取心脏重心的方法和系统 | |
| US20200170713A1 (en) | Virtual stent placement apparatus, virtual stent placement method, and virtual stent placement program | |
| CN110929604B (zh) | 基于造影图像的流速的筛选方法、装置、系统和存储介质 | |
| CN112487342B (zh) | 精确获取狭窄病变区间的方法、系统及存储介质 | |
| CN113902693A (zh) | 图像分割方法及计算机可读存储介质 | |
| KR20180097037A (ko) | 관상동맥의 시작점을 자동으로 추출하기 위한 방법 및 그 장치 |
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 |