CN112419279B - 二维图像选取及三维血管合成的方法和存储介质 - Google Patents
二维图像选取及三维血管合成的方法和存储介质 Download PDFInfo
- Publication number
- CN112419279B CN112419279B CN202011342206.4A CN202011342206A CN112419279B CN 112419279 B CN112419279 B CN 112419279B CN 202011342206 A CN202011342206 A CN 202011342206A CN 112419279 B CN112419279 B CN 112419279B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- blood vessel
- vessel
- image
- line
- 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 40
- 210000004204 blood vessel Anatomy 0.000 title claims description 120
- 230000002194 synthesizing effect Effects 0.000 title claims description 10
- 230000002792 vascular Effects 0.000 claims abstract description 41
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 30
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 14
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 14
- 238000002586 coronary angiography Methods 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 11
- 238000010586 diagram Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000002285 radioactive effect Effects 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 3
- 239000008280 blood Substances 0.000 claims description 2
- 210000004369 blood Anatomy 0.000 claims description 2
- 238000001308 synthesis method Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 7
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 230000036262 stenosis 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
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002107 myocardial effect Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 238000004364 calculation method Methods 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
- 238000003708 edge detection Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000000605 extraction Methods 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
- 210000004165 myocardium Anatomy 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Quality & Reliability (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请提供了一种二维图像选取及三维血管合成的方法和存储介质，包括：加载多组冠状动脉二维造影图像序列组；选择一组所述冠状动脉二维造影图像序列组；每组所述冠状动脉二维造影图像序列组内显示一系列的缩略图；从至少两组所述冠状动脉二维造影图像序列组内选择感兴趣的所述缩略图，所述缩略图会显示出放大图，即获得用于三维血管合成的冠状动脉二维造影图像。本申请实现了一种二维图像选取及三维血管合成的方法和存储介质，这种方法能有效降低血管信息获取的难度和误差，提高血管信息获得的精度，具有更好的实用性。
Description
技术领域
本发明涉及冠状动脉医学技术领域，特别是涉及一种二维图像选取及三维血管合成的方法和存储介质。
背景技术
人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块，继而导致血管狭窄；特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足，诱发冠心病、心绞痛等病症，对人类的健康造成严重威胁。据统计，我国现有冠心病患者约1100万人，心血管介入手术治疗患者数量每年增长大于10％。
冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度，但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性，1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve，FFR)，经过长期的基础与临床研究，FFR已成为冠脉狭窄功能性评价的金标准。
血流储备分数(FFR)通常是指心肌血流储备分数，定义为病变冠脉能为心肌提供的最大血流与该冠脉完全正常时最大供血流量之比，研究表明，在冠脉最大充血状态下，血流量的比值可以用压力值来代替。即FFR值的测量可在冠脉最大充血状态下，通过压力传感器对冠脉远端狭窄处的压力和冠脉狭窄近端压力进行测定继而计算得出。
现有技术中通常都是将图像处理算法应用在二维造影序列图像中，先采取全自动、半自动或者手动的方法计算出血管的路径，再用边缘检测方法得到血管的轮廓，从而进一步算出血管的直径。但该方法存在的问题在于血管周围各种器官的重叠和交叉，这些情况都会给算法造成不确定因素。
发明内容
本发明提供了一种二维图像选取及三维血管合成的方法和存储介质，以避免在二维造影序列图像中血管周围各种器官的重叠和交叉给算法造成不确定因素。
为实现上述目的，
第一方面，本申请提供了一种三维血管的合成方法，包括：
加载多组冠状动脉二维造影图像序列组；
选择一组所述冠状动脉二维造影图像序列组；
每组所述冠状动脉二维造影图像序列组内显示一系列的缩略图；
从至少两组所述冠状动脉二维造影图像序列组内选择感兴趣的所述缩略图，所述缩略图会显示出放大图，即获得用于三维血管合成的冠状动脉二维造影图像；
其中，每组所述冠状动脉二维造影图像序列组在页面显示时，默认选取中间帧作为缩略图显示；
获取至少两幅拍摄角度不同、感兴趣的冠状动脉二维造影图像的图像信息，包括拍摄角度和拍摄距离；
根据所述冠状动脉二维造影图像的图像信息获取三维血管中心线和三维血管半径；
依据所述三维血管中心线和所述三维血管半径合成三维血管；
根据所述冠状动脉二维造影图像的图像信息获取三维血管半径包括：
从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线；
根据所述二维血管中心线获取二维血管轮廓线；
根据所述二维血管轮廓线获取每幅所述感兴趣的冠状动脉二维造影图像中的二维血管半径；
根据所述二维血管半径获取所述三维血管半径，具体公式为：
其中，R表示三维血管半径，r1、r2、rn分别表示第一幅、第二幅、第n幅感兴趣的二维造影图像的二维血管半径；
所述根据所述二维血管中心线获取二维血管轮廓线包括：
根据冠状动脉二维造影图像提取二维血管中心线；
根据所述二维血管中心线，获得拉直血管图像；
在所述拉直血管图像上，设定血管直径阈值D阈；
根据所述D阈，在血管中心直线两侧生成血管预设轮廓线；
将所述血管预设轮廓线向所述血管中心直线逐级靠拢，获取拉直后血管的轮廓线；
将所述拉直后血管的轮廓线投射回提取所述二维血管中心线的图像上，获得二维血管轮廓线。
优选的，根据所述冠状动脉二维造影图像的图像信息获取三维血管中心线的方法，包括：
所述从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线；
将放射源投影到三维空间内形成放射点；
所述二维血管中心线投影到三维空间内；
所有三维空间内的点均与所述放射点连线，会产生一系列的交叉点；
将所述交叉点依次连接，得到所述三维血管中心线。
优选的，所述从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线的方法，包括：
读取冠状动脉二维造影图像；
获取感兴趣的血管段；
拾取所述感兴趣的血管段的起始点、种子点和结束点；
分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割，得到至少两个局部血管区域图；
从每个所述局部血管区域图中提取至少一条血管局部路径线；
将每个所述局部血管区域图上相对应的血管局部路径线连接，获得至少一条血管路径线；
选取一条所述血管路径线作为所述二维血管中心线。
优选的，所述从每个所述局部血管区域图上提取至少一条血管局部路径线的方法包括：
对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图；
对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线。
优选的，所述对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线的方法包括：
对所述粗略血管图进行网格划分；
沿着所述起始点至所述结束点的血管延伸方向，搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点，搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点，所述第三个点重复上述步骤，直至最短时间路径到达结束点，其中，n为大于等于1的正整数；
按照搜索顺序，从所述起始点至所述结束点的血管延伸方向连线，获得至少一条血管局部路径线。
优选的，所述选取一条所述血管路径线作为所述二维血管中心线的方法包括：
如果血管路径线为两条或两条以上，则对每条血管路径线从所述起始点至所述结束点所用的时间求和；
取用时最少的所述血管路径线作为所述二维血管中心线。
优选的，所述依据所述三维血管中心线和所述三维血管半径合成三维血管的方法包括：
每个所述三维血管中心线上的点，均沿着对应的三维血管半径在所述三维空间内画图，得到多个边缘点，将所述边缘点依次连接，得到一个近似圆的多边形；
将相邻两个所述多边形上的点按照直角三角形的形式依次连接，得到三维血管。
第二方面，本申请提供了一种计算机存储介质，该计算机存储介质中包含计算机程序，计算机程序被处理器执行时用于实现权利要求1～7任一所述的一种三维血管的合成方法。
本申请实施例提供的方案带来的有益效果至少包括：
本申请提供了一种三维血管的合成方法，在由冠状动脉二维图像计算得到血管信息的过程中，通过降低其他部分对血管的影响，能够更加容易的获得血管半径，且提高了血管半径的精确度，提高算法的鲁棒性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
下面对附图标记进行说明：
图1为本申请的用于三维血管合成的冠状动脉二维图像选取的方法的一个实施例的流程图；
图2为本申请的用于三维血管合成的冠状动脉二维图像选取的方法的另一实施例的流程图；
图3为本申请的S600的流程图；
图4为本申请的S620的流程图；
图5为本申请的S630的流程图；
图6为本申请的S700的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下将以图式揭露本发明的多个实施方式，为明确说明起见，许多实务上的细节将在以下叙述中一并说明。然而，应了解到，这些实务上的细节不应用以限制本发明。也就是说，在本发明的部分实施方式中，这些实务上的细节是非必要的。此外，为简化图式起见，一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。
现有技术中通过血管三维模型计算血管评价参数时常常需要提取血管轮廓线，由于血管存在卷曲、且边缘不清晰的问题，导致血管轮廓提取特别困难，且运算数据庞大、繁冗，因此如何快速提取血管轮廓线，以及提取的准确度一直是技术人员需要解决的问题。
实施例1：
如图1所示，本申请为了解决上述问题，本申请提供了一种用于三维血管合成的冠状动脉二维图像选取的方法，包括：
S100，加载多组冠状动脉二维造影图像序列组；
S200，选择一组所述冠状动脉二维造影图像序列组；
S300，每组所述冠状动脉二维造影图像序列组内显示一系列的缩略图；
S400，从至少两组所述冠状动脉二维造影图像序列组内选择感性趣的所述缩略图，所述缩略图会显示出放大图，即获得用于三维血管合成的冠状动脉二维造影图像。
本申请的一个实施例中，每组所述冠状动脉二维造影图像序列组在页面显示时，默认选取中间帧作为缩略图显示。
实施例2：
如图2所示，本申请提供了一种三维血管的合成方法，包括：
S500，根据实施例1的方法，获取至少两幅拍摄角度不同、感兴趣的冠状动脉二维造影图像的图像信息，包括拍摄角度和拍摄距离；
S600，如图3所示，根据所述冠状动脉二维造影图像的图像信息获取三维血管中心线和三维血管半径，包括：
S610，从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线，包括：
S611，读取冠状动脉二维造影图像；
S612，获取感兴趣的血管段；
S613，拾取所述感兴趣的血管段的起始点、种子点和结束点；
S614，分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割，得到至少两个局部血管区域图；
S615，从每个所述局部血管区域图中提取至少一条血管局部路径线，包括：
对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图，包括：在每幅所述局部血管区域图中，以所述感兴趣的血管段作为前景，其他区域作为背景，强化所述前景，弱化所述背景，得到对比强烈的所述粗略血管图。
对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线，包括：
对所述粗略血管图进行网格划分；
沿着所述起始点至所述结束点的血管延伸方向，搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点，搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点，所述第三个点重复上述步骤，直至最短时间路径到达结束点，其中，n为大于等于1的正整数；
按照搜索顺序，从所述起始点至所述结束点的血管延伸方向连线，获得至少一条血管局部路径线。
S616，选取一条所述血管路径线作为所述二维血管中心线，包括：
如果血管路径线为两条或两条以上，则对每条血管路径线从所述起始点至所述结束点所用的时间求和；
取用时最少的所述血管路径线作为所述二维血管中心线。
S620，如图4所示，根据每幅所述冠状动脉二维造影图像的拍摄角度，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线，包括：
S621，将放射源投影到所述三维空间内形成放射点；
S622，所述二维血管中心线投影到三维空间内；
S623，所有三维空间内的点均与所述放射点连线，会产生一系列的交叉点；
S624，将所述交叉点依次连接，得到所述三维血管中心线；
S630，如图5所示，根据所述二维血管中心线获取二维血管轮廓线，包括：
S631，根据冠状动脉二维造影图像提取二维血管中心线；
S632，根据所述二维血管中心线，获得拉直血管图像，包括：
将所述二维血管中心线拉直，获得血管中心直线；
沿着所述起始点至所述结束点的血管延伸方向，将所述局部血管区域图分为x个单元，其中x为正整数；
将每个所述单元的二维血管中心线沿着所述血管中心直线对应设置；
对应设置后的图像为所述拉直血管图像。
S633，在所述拉直血管图像上，设定血管直径阈值D阈；
S634，根据所述D阈，在血管中心直线两侧生成血管预设轮廓线；
S635，将所述血管预设轮廓线向所述血管中心直线逐级靠拢，获取拉直后血管的轮廓线，包括：
将所述血管预设轮廓线分成y个单元，其中y为正整数；
获取每个所述单元的位于每条所述血管预设轮廓线上的z个点；
沿着垂直于所述血管中心直线方向，将z个点分别向所述血管中心直线分级靠拢，产生z个靠拢点，其中z为正整数；
设定RGB差值阈值为ΔRGB阈，沿着垂直于所述血管中心直线方向，每次靠拢均将所述靠拢点的RGB值与所述血管中心直线上的点的RGB值作比较，当差值小于等于ΔRGB阈时，则所述靠拢点停止向所述血管中心直线靠拢；
获取所述靠拢点作为轮廓点；
依次连接所述轮廓点形成的平滑曲线即为所述拉直后血管的轮廓线。
S636，将所述拉直后血管的轮廓线投射回提取所述二维血管中心线的图像上，获得二维血管轮廓线。
S640，根据所述二维血管轮廓线获取每幅所述感兴趣的冠状动脉二维造影图像中的二维血管半径；
S650，根据所述二维血管半径获取所述三维血管半径，具体公式为：
其中，R表示三维血管半径，r1、r2、rn分别表示第一幅、第二幅、第n幅感兴趣的二维造影图像的二维血管半径。
S700，如图6所示，依据所述三维血管中心线和所述三维血管半径合成三维血管，包括：
S710，每个所述三维血管中心线上的点，均沿着对应的三维血管半径在所述三维空间内画图，得到多个边缘点，将所述边缘点依次连接，得到一个近似圆的多边形；
S720，将相邻两个所述多边形上的点按照直角三角形的形式依次连接，得到三维血管。
本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的用于三维血管合成的冠状动脉二维图像选取的方法。
所属技术领域的技术人员知道，本发明的各个方面可以实现为***、方法或计算机程序产品。因此，本发明的各个方面可以具体实现为以下形式，即：完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等)，或硬件和软件方面结合的实施方式，这里可以统称为“电路”、“模块”或“***”。此外，在一些实施例中，本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式，该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或***的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。
例如，可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件，可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作***执行的多个软件指令。在本发明的示例性实施例中，由数据处理器来执行如本文的根据方法和/或***的示例性实施例的一个或多个任务，诸如用于执行多个指令的计算平台。可选地，该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器，例如，磁硬盘和/或可移动介质。可选地，也提供了一种网络连接。可选地也提供显示器和/或用户输入设备，诸如键盘或鼠标。
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的***、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项：
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行***、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行***、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括(但不限于)无线、有线、光缆、RF等等，或者上述的任意合适的组合。
例如，可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码，包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言，诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。
也可以把这些计算机程序指令存储在计算机可读介质中，这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作，从而，存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。
还可将计算机程序指令加载到计算机(例如，冠状动脉分析***)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程，使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。
本发明的以上的具体实例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种三维血管的合成方法，其特征在于，包括：
加载多组冠状动脉二维造影图像序列组；
选择一组所述冠状动脉二维造影图像序列组；
每组所述冠状动脉二维造影图像序列组内显示一系列的缩略图；
从至少两组所述冠状动脉二维造影图像序列组内选择感兴趣的所述缩略图，所述缩略图会显示出放大图，即获得用于三维血管合成的冠状动脉二维造影图像；
其中，每组所述冠状动脉二维造影图像序列组在页面显示时，默认选取中间帧作为缩略图显示；
获取至少两幅拍摄角度不同、感兴趣的冠状动脉二维造影图像的图像信息，包括拍摄角度和拍摄距离；
根据所述冠状动脉二维造影图像的图像信息获取三维血管中心线和三维血管半径；
依据所述三维血管中心线和所述三维血管半径合成三维血管；
根据所述冠状动脉二维造影图像的图像信息获取三维血管半径包括：
从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线；
根据所述二维血管中心线获取二维血管轮廓线；
根据所述二维血管轮廓线获取每幅所述感兴趣的冠状动脉二维造影图像中的二维血管半径；
根据所述二维血管半径获取所述三维血管半径，具体公式为：
其中，R表示三维血管半径，r1、r2、rn分别表示第一幅、第二幅、第n幅感兴趣的二维造影图像的二维血管半径；
所述根据所述二维血管中心线获取二维血管轮廓线包括：
根据冠状动脉二维造影图像提取二维血管中心线；
根据所述二维血管中心线，获得拉直血管图像；
在所述拉直血管图像上，设定血管直径阈值D阈；
根据所述D阈，在血管中心直线两侧生成血管预设轮廓线；
将所述血管预设轮廓线向所述血管中心直线逐级靠拢，获取拉直后血管的轮廓线；
将所述拉直后血管的轮廓线投射回提取所述二维血管中心线的图像上，获得二维血管轮廓线。
2.根据权利要求1所述的一种三维血管的合成方法，其特征在于，根据所述冠状动脉二维造影图像的图像信息获取三维血管中心线的方法，包括：
从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线；
将放射源投影到三维空间内形成放射点；
所述二维血管中心线投影到三维空间内；
所有三维空间内的点均与所述放射点连线，会产生一系列的交叉点；
将所述交叉点依次连接，得到所述三维血管中心线。
3.根据权利要求2所述的一种三维血管的合成方法，其特征在于，所述从每幅所述感兴趣的冠状动脉二维造影图像中提取一条二维血管中心线的方法，包括：
读取冠状动脉二维造影图像；
获取感兴趣的血管段；
拾取所述感兴趣的血管段的起始点、种子点和结束点；
分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割，得到至少两个局部血管区域图；
从每个所述局部血管区域图中提取至少一条血管局部路径线；
将每个所述局部血管区域图上相对应的血管局部路径线连接，获得至少一条血管路径线；
选取一条所述血管路径线作为所述二维血管中心线。
4.根据权利要求3所述的一种三维血管的合成方法，其特征在于，所述从每个所述局部血管区域图上提取至少一条血管局部路径线的方法包括：
对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图；
对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线。
5.根据权利要求4所述的一种三维血管的合成方法，其特征在于，所述对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线的方法包括：
对所述粗略血管图进行网格划分；
沿着所述起始点至所述结束点的血管延伸方向，搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点，搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点，所述第三个点重复上述步骤，直至最短时间路径到达结束点，其中，n为大于等于1的正整数；
按照搜索顺序，从所述起始点至所述结束点的血管延伸方向连线，获得至少一条血管局部路径线。
6.根据权利要求5所述的一种三维血管的合成方法，其特征在于，所述选取一条所述血管路径线作为所述二维血管中心线的方法包括：
如果血管路径线为两条或两条以上，则对每条血管路径线从所述起始点至所述结束点所用的时间求和；
取用时最少的所述血管路径线作为所述二维血管中心线。
7.根据权利要求6所述的一种三维血管的合成方法，其特征在于，所述依据所述三维血管中心线和所述三维血管半径合成三维血管的方法包括：
每个所述三维血管中心线上的点，均沿着对应的三维血管半径在所述三维空间内画图，得到多个边缘点，将所述边缘点依次连接，得到一个近似圆的多边形；
将相邻两个所述多边形上的点按照直角三角形的形式依次连接，得到三维血管。
8.一种计算机存储介质，其特征在于，该计算机存储介质中包含计算机程序，计算机程序被处理器执行时用于实现权利要求1～7任一所述的一种三维血管的合成方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011342206.4A CN112419279B (zh) | 2020-11-25 | 2020-11-25 | 二维图像选取及三维血管合成的方法和存储介质 |
| PCT/CN2020/131700 WO2022109904A1 (zh) | 2020-11-25 | 2020-11-26 | 二维图像选取及三维血管合成的方法和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011342206.4A CN112419279B (zh) | 2020-11-25 | 2020-11-25 | 二维图像选取及三维血管合成的方法和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112419279A CN112419279A (zh) | 2021-02-26 |
| CN112419279B true CN112419279B (zh) | 2024-02-23 |
Family
ID=74842431
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011342206.4A Active CN112419279B (zh) | 2020-11-25 | 2020-11-25 | 二维图像选取及三维血管合成的方法和存储介质 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112419279B (zh) |
| WO (1) | WO2022109904A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117058328B (zh) * | 2023-10-11 | 2024-01-09 | 杭州脉流科技有限公司 | 冠状动脉血管树分级方法、设备、存储介质和程序产品 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6047080A (en) * | 1996-06-19 | 2000-04-04 | Arch Development Corporation | Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images |
| CN108242075A (zh) * | 2018-01-05 | 2018-07-03 | 苏州润迈德医疗科技有限公司 | 一种基于x射线冠脉造影图像的多角度血管重建方法 |
| CN110287956A (zh) * | 2019-06-13 | 2019-09-27 | 北京理工大学 | 血管中心线自动匹配方法及装置 |
| CN110889896A (zh) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | 获取血管狭窄病变区间及三维合成方法、装置和*** |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9934566B2 (en) * | 2015-07-14 | 2018-04-03 | Siemens Healthcare Gmbh | 3-D vessel tree surface reconstruction method |
| CN108133478B (zh) * | 2018-01-11 | 2021-11-23 | 苏州润迈德医疗科技有限公司 | 一种提取冠状动脉血管中心线的方法 |
| CN109461138B (zh) * | 2018-09-29 | 2020-10-27 | 深圳睿心智能医疗科技有限公司 | 参数计算方法、***、可读存储介质及计算机设备 |
| CN111161342B (zh) * | 2019-12-09 | 2023-08-29 | 杭州脉流科技有限公司 | 基于冠脉造影图像获取血流储备分数的方法、装置、设备、***及可读存储介质 |
| CN111754506B (zh) * | 2020-07-01 | 2024-02-06 | 杭州脉流科技有限公司 | 基于腔内影像的冠脉狭窄率计算方法、装置、***和计算机存储介质 |
-
2020
- 2020-11-25 CN CN202011342206.4A patent/CN112419279B/zh active Active
- 2020-11-26 WO PCT/CN2020/131700 patent/WO2022109904A1/zh active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6047080A (en) * | 1996-06-19 | 2000-04-04 | Arch Development Corporation | Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images |
| CN108242075A (zh) * | 2018-01-05 | 2018-07-03 | 苏州润迈德医疗科技有限公司 | 一种基于x射线冠脉造影图像的多角度血管重建方法 |
| CN110287956A (zh) * | 2019-06-13 | 2019-09-27 | 北京理工大学 | 血管中心线自动匹配方法及装置 |
| CN110889896A (zh) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | 获取血管狭窄病变区间及三维合成方法、装置和*** |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022109904A1 (zh) | 2022-06-02 |
| CN112419279A (zh) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112116615B (zh) | 根据血管中心线获取血管轮廓线的方法和装置 | |
| CN112164020B (zh) | 精确提取血管中心线的方法、装置、分析***和存储介质 | |
| US9375191B2 (en) | Method and apparatus for determining three-dimensional reconstruction of an object | |
| CN112132949B (zh) | 具有狭窄病变区间的血管三维建模方法、装置和*** | |
| CN112419484B (zh) | 三维血管合成方法、***及冠状动脉分析***和存储介质 | |
| US11076770B2 (en) | System for vascular assessment | |
| JP5129480B2 (ja) | 管状臓器の３次元再構成を行うシステム及び血管撮像装置の作動方法 | |
| CN112419462B (zh) | 三维血管的渲染合成方法、***及存储介质 | |
| JP2004243117A (ja) | 身体構造の物理パラメータを求める方法 | |
| WO2022109905A1 (zh) | 调节血管轮廓及中心线的方法及存储介质 | |
| EP3606433B1 (en) | Standardized coronary artery disease metric | |
| WO2021097821A1 (zh) | 从冠状动脉二维造影图像中提取血管中心线的方法和装置 | |
| Van Hamersvelt et al. | Diagnostic performance of on-site coronary CT angiography–derived fractional flow reserve based on patient-specific lumped parameter models | |
| CN112419279B (zh) | 二维图像选取及三维血管合成的方法和存储介质 | |
| CN112419277B (zh) | 三维血管中心线合成方法、***及存储介质 | |
| EP3886702B1 (en) | Most relevant x-ray image selection for hemodynamic simulation | |
| CN111815584B (zh) | 基于ct序列图像获取心脏重心的方法和*** | |
| EP3607527B1 (en) | Quantitative evaluation of time-varying data | |
| JP7408381B2 (ja) | 画像処理装置、プログラム及び方法 |
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 |