CN112419277B - 三维血管中心线合成方法、系统及存储介质 - Google Patents
三维血管中心线合成方法、系统及存储介质 Download PDFInfo
- Publication number
- CN112419277B CN112419277B CN202011337612.1A CN202011337612A CN112419277B CN 112419277 B CN112419277 B CN 112419277B CN 202011337612 A CN202011337612 A CN 202011337612A CN 112419277 B CN112419277 B CN 112419277B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- point
- blood vessel
- vessel
- image
- 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
- 230000002792 vascular Effects 0.000 title claims abstract description 46
- 238000001308 synthesis method Methods 0.000 title claims abstract description 13
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 35
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 25
- 238000002601 radiography Methods 0.000 claims abstract description 10
- 210000004204 blood vessel Anatomy 0.000 claims description 94
- 238000000034 method Methods 0.000 claims description 46
- 238000002586 coronary angiography Methods 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 238000010586 diagram Methods 0.000 claims description 12
- 238000003786 synthesis reaction Methods 0.000 claims description 12
- 230000002285 radioactive effect Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 9
- 230000005855 radiation Effects 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 7
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000017531 blood circulation Effects 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
- 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
- 230000002107 myocardial effect Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 238000002583 angiography Methods 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
- 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
- 210000004165 myocardium Anatomy 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
- G06T7/11—Region-based segmentation
-
- 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)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请提供了三维血管中心线合成方法、系统及存储介质,包括:获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息;从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线;根据每幅所述冠状动脉二维造影图像的图像信息,将每条所述二维血管中心线投影到三维空间内,合成所述三维血管中心线。本申请的三维血管中心线合成方法摆脱了对多角度冠状动脉二维造影图像成像时间相同、空间相同的要求,在保证精度的前提下,适用范围更广。
Description
技术领域
本发明涉及冠状动脉医学技术领域,特别是涉及三维血管中心线合成方法、系统及存储介质。
背景技术
人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块,继而导致血管狭窄;特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足,诱发冠心病、心绞痛等病症,对人类的健康造成严重威胁。据统计,我国现有冠心病患者约1100万人,心血管介入手术治疗患者数量每年增长大于10%。
冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度,但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性,1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve,FFR),经过长期的基础与临床研究,FFR已成为冠脉狭窄功能性评价的金标准。
血流储备分数(FFR)通常是指心肌血流储备分数,定义为病变冠脉能为心肌提供的最大血流与该冠脉完全正常时最大供血流量之比,研究表明,在冠脉最大充血状态下,血流量的比值可以用压力值来代替。即FFR值的测量可在冠脉最大充血状态下,通过压力传感器对冠脉远端狭窄处的压力和冠脉狭窄近端压力进行测定继而计算得出。
现有存在的问题在于多角度冠状动脉二维造影图像成像的过程中,大部分的DSA成像设备,在同一时间只能拍同一角度的冠状动脉二维造影图像成像,也就是说,多个不同角度的冠状动脉二维造影图像成像必须由不同时间得到。不符合直接三维血管中心线合成方法对图像的要求,因此合成结果大多失真。
发明内容
本发明提供了一种三维血管中心线合成方法、系统及存储介质,以解决三维血管中心线合成要求多角度冠状动脉二维造影图像成像时间和空间必须都相同的问题。
为实现上述目的,第一方面,本申请提供了一种三维血管中心线的合成方法,包括:
获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息;
从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线;
根据每幅所述冠状动脉二维造影图像的图像信息,将每条所述二维血管中心线投影到三维空间内,合成所述三维血管中心线。
可选地,上述的三维血管中心线的合成方法,所述获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息的方法包括:
获取至少两组拍摄角度不同的冠状动脉二维造影图像组;
读取每组所述冠状动脉二维造影图像组的图像信息,包括拍摄角度和探测距离;
根据所述探测距离,分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。
可选地,上述的三维血管中心线的合成方法,所述根据每幅所述冠状动脉二维造影图像的图像信息,将每条所述二维血管中心线投影到三维空间内,合成所述三维血管中心线的方法包括:
以心脏为坐标原点,建立三维坐标系;
获取每幅所述感兴趣的二维造影图像的左右角度α、前后角度β,以及人体与平板探测器之间的距离S,所述二维血管中心线上各点的坐标为(x,y);
将各点(x,y)投影到三维空间内,获得一系列三维坐标点P,坐标为(x’、y’、z’);
将放射源投影到所述三维空间内形成放射点R;
将每个所述三维坐标点P均与所述放射点R进行连线,从PR连线上获得所述三维血管中心线上的点,将所述三维血管中心线上的点依次连接获得所述三维血管中心线。
可选地,上述的三维血管中心线的合成方法,所述将各点(x,y)投影到三维空间内,获得一系列三维坐标点P,坐标为(x’、y’、z’)的方法,包括:
将各点(x,y)绕y轴旋转获得(x’、y’、z’)系列点,具体公式为:
将(x’、y’、z’)系列点绕x轴旋转获得一系列三维坐标点P,坐标为(x’、y’、z’);
可选地,上述的三维血管中心线的合成方法,所述将将放射源投影到所述三维空间内形成放射点R的方法,包括:
获得每幅所述感兴趣的二维造影图像中人体与放射源之间的距离S’;
根据公式
其中,R在三维空间内的坐标为(a,b,c)。
可选地,上述的三维血管中心线的合成方法,所述将每个所述三维坐标点P均与所述放射点R进行连线,从PR连线上获得所述三维血管中心线上的点,将所述三维血管中心线上的点依次连接获得所述三维血管中心线的方法,包括:
将从同一幅所述感兴趣的二维造影图像中获得的一系列三维坐标点P和放射点R对应连接,获得多条PR直线;
获取所述血管的同一位置的两条PR直线之间的最小距离的点,分别为A点和B点;
将所述A点与所述B点连线,获取所述AB线段的中点作为所述三维血管的中心线上的点;
将获得的一系列的所述三维血管的中心线上的点依次连接,获得所述三维血管中心线。
可选地,上述的三维血管中心线的合成方法,所述分别从每幅所述冠状动脉二维造影图像中提取一条二维血管中心线的方法包括:
读取冠状动脉二维造影图像;
获取感兴趣的血管段;
拾取所述感兴趣的血管段的起始点、种子点和结束点;
分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割,得到至少两个局部血管区域图;
从每个所述局部血管区域图中提取至少一条血管局部路径线;
将每个所述局部血管区域图上相对应的血管局部路径线连接,获得至少一条所述血管路径线;
选取一条所述血管路径线作为所述二维血管中心线。
可选地,上述的三维血管中心线的合成方法,所述从每个所述局部血管区域图上提取至少一条血管局部路径线的方法包括:
对所述局部血管区域图做图像增强处理,得到对比强烈的粗略血管图;
对所述粗略血管图做网格划分,沿着所述起始点至所述结束点方向,提取至少一条血管局部路径线。
可选地,上述的三维血管中心线的合成方法,所述对所述局部血管区域图做图像增强处理,得到对比强烈的粗略血管图的方法,包括:
在每幅所述局部血管区域图中,以所述感兴趣的血管段作为前景,其他区域作为背景,强化所述前景,弱化所述背景,得到对比强烈的所述粗略血管图。
可选地,上述的三维血管中心线的合成方法,所述对所述粗略血管图做网格划分,沿着所述起始点至所述结束点方向,提取至少一条血管局部路径线的方法包括:
对所述粗略血管图进行网格划分;
沿着所述起始点至所述结束点的血管延伸方向,搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点,搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点,所述第三个点重复上述步骤,直至最短时间路径到达结束点,其中,n为大于等于1的正整数;
按照搜索顺序,从所述起始点至所述结束点的血管延伸方向连线,获得至少一条血管局部路径线。
可选地,上述的三维血管中心线的合成方法,所述选取一条所述血管路径线作为所述二维血管中心线的方法包括:
如果血管路径线为两条或两条以上,则对每条血管路径线从所述起始点至所述结束点所用的时间求和;
取用时最少的所述血管路径线作为所述二维血管中心线。
第二方面,本申请提供了一种三维血管中心线合成系统,用于上述的三维血管中心线的合成方法,包括:图像读取装置,与所述图像读取装置连接的二维血管中心线提取装置和所述三维血管中心线获取装置,所述二维血管中心线提取装置与所述三维血管中心线获取装置连接;
所述图像读取装置,用于获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息;
所述二维血管中心线提取装置,用于接收所述图像读取装置发送的冠状动脉二维造影图像,从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线;
所述三维血管中心线获取装置,用于接收所述二维血管中心线提取装置发送的二维血管中心线,以及接收所述图像读取装置发送的冠状动脉二维造影图像和图像信息,根据每幅所述冠状动脉二维造影图像的图像信息,将每条所述二维血管中心线投影到三维空间内,合成所述三维血管中心线。
可选地,上述的三维血管合成系统,所述图像读取装置包括:依次连接的二维图像获取结构、图像信息获取结构、筛选感兴趣图像结构,所述筛选感兴趣图像结构与所述二维图像获取结构连接;
所述二维图像获取结构,用于获取至少两组拍摄角度不同的冠状动脉二维造影图像组;
所述图像信息获取结构,用于读取所述二维图像获取结构传递的每组所述冠状动脉二维造影图像组的图像信息,包括拍摄角度和探测距离;
所述筛选感兴趣图像结构,用于根据所述探测距离,分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。
第三方面,本申请提供了一种冠状动脉分析系统,包括:上述的三维血管合成系统。
第四方面,本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述的三维血管的合成方法。
本申请实施例提供的方案带来的有益效果至少包括:
本申请提供了一种三维血管的合成方法,摆脱了对多角度冠状动脉二维造影图像成像时间相同、空间相同的要求,在保证精度的前提下,适用范围更广。该方法适应于现有部署在实际场景中的成像设备。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本申请的三维血管中心线的合成方法的流程图;
图2为本申请的S100的流程图;
图3为本申请的S200的流程图;
图4为本申请的S260的流程图;
图5为本申请的S262的流程图;
图6为本申请的S270的流程图;
图7为本申请的S300的流程图;
图8为本申请的S350的流程图;
图9为本申请的三维血管中心线合成系统的结构框图;
图10为本申请的图像读取装置100的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下将以图式揭露本发明的多个实施方式,为明确说明起见,许多实务上的细节将在以下叙述中一并说明。然而,应了解到,这些实务上的细节不应用以限制本发明。也就是说,在本发明的部分实施方式中,这些实务上的细节是非必要的。此外,为简化图式起见,一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。
实施例1:
如图1所示,本申请为了解决上述问题,本申请提供了一种三维血管中心线的合成方法,包括:
S100,获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息;
S200,从每幅感兴趣的二维造影图像中提取一条二维血管中心线;
S300,根据每幅冠状动脉二维造影图像的图像信息,将每条二维血管中心线投影到三维空间内,合成三维血管中心线。
实施例2:
本申请为了解决上述问题,本申请提供了一种三维血管中心线的合成方法,包括:
S100,如图2所示,获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息,包括:
S110,获取至少两组拍摄角度不同的冠状动脉二维造影图像组;
S120,读取每组冠状动脉二维造影图像组的图像信息,包括拍摄角度和探测距离;
S130,根据探测距离,分别从每组冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。
S200,如图3所示,从每幅感兴趣的二维造影图像中提取一条二维血管中心线,包括:
S210,读取冠状动脉二维造影图像;
S220,获取感兴趣的血管段;
S230,拾取感兴趣的血管段的起始点、种子点和结束点;
S240,分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割,得到至少两个局部血管区域图;
S250,从每个局部血管区域图中提取至少一条血管局部路径线;
S260,如图4所示,将每个局部血管区域图上相对应的血管局部路径线连接,获得至少一条血管路径线,包括:
S261,对局部血管区域图做图像增强处理,得到对比强烈的粗略血管图,包括:在每幅局部血管区域图中,以感兴趣的血管段作为前景,其他区域作为背景,强化前景,弱化背景,得到对比强烈的粗略血管图。
S262,如图5所示,对粗略血管图做网格划分,沿着起始点至结束点方向,提取至少一条血管局部路径线,包括:
S2621,对粗略血管图进行网格划分;
S2622,沿着起始点至结束点的血管延伸方向,搜索起始点与周边n个网格上的交叉点的最短时间路径作为第二个点,搜索第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点,第三个点重复上述步骤,直至最短时间路径到达结束点,其中,n为大于等于1的正整数;
S2623,按照搜索顺序,从起始点至结束点的血管延伸方向连线,获得至少一条血管局部路径线。
S270,如图6所示,选取一条血管路径线作为二维血管中心线,包括:
S271,如果血管路径线为两条或两条以上,则对每条血管路径线从起始点至结束点所用的时间求和;
S272,取用时最少的血管路径线作为二维血管中心线。
S300,如图7所示,根据每幅冠状动脉二维造影图像的图像信息,将每条二维血管中心线投影到三维空间内,合成三维血管中心线,包括:
S310,以心脏为坐标原点,建立三维坐标系;
S320,获取每幅感兴趣的二维造影图像的左右角度α、前后角度β,以及人体与平板探测器之间的距离S,二维血管中心线上各点的坐标为(x,y);
S330,将各点(x,y)投影到三维空间内,获得一系列三维坐标点P,坐标为(x’、y’、z’),包括:
将各点(x,y)绕y轴旋转获得(x’、y’、z’)系列点,具体公式为:
将(x’、y’、z’)系列点绕x轴旋转获得一系列三维坐标点P,坐标为(x’、y’、z’);
S340,将放射源投影到三维空间内形成放射点R,包括:
获得每幅感兴趣的二维造影图像中人体与放射源之间的距离S’;
根据公式
其中,R在三维空间内的坐标为(a,b,c)。
S350,如图8所示,将每个三维坐标点P均与放射点R进行连线,从PR连线上获得三维血管中心线上的点,将三维血管中心线上的点依次连接获得三维血管中心线,包括:
S351,将从同一幅感兴趣的二维造影图像中获得的一系列三维坐标点P和放射点R对应连接,获得多条PR直线;
S352,获取血管的同一位置的两条PR直线之间的最小距离的点,分别为A点和B点;
S353,将A点与B点连线,获取AB线段的中点作为三维血管的中心线上的点;
S354,将获得的一系列的三维血管的中心线上的点依次连接,获得三维血管中心线。
如图9所示,本申请提供了一种三维血管中心线合成系统,用于上述的三维血管中心线的合成方法,包括:图像读取装置100,与图像读取装置100连接的二维血管中心线提取装置200和三维血管中心线获取装置300,二维血管中心线提取装置200与三维血管中心线获取装置300连接;图像读取装置100,用于获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息;二维血管中心线提取装置200用于接收图像读取装置发送的冠状动脉二维造影图像,从每幅感兴趣的二维造影图像中提取一条二维血管中心线;三维血管中心线获取装置300用于接收二维血管中心线提取装置200发送的二维血管中心线,以及接收图像读取装置100发送的冠状动脉二维造影图像和图像信息,根据每幅冠状动脉二维造影图像的图像信息,将每条二维血管中心线投影到三维空间内,合成三维血管中心线。
如图10所示,本申请的一个实施例中,图像读取装置100包括:依次连接的二维图像获取结构110、图像信息获取结构120、筛选感兴趣图像结构130,筛选感兴趣图像结构130与二维图像获取结构110连接;二维图像获取结构110用于获取至少两组拍摄角度不同的冠状动脉二维造影图像组;图像信息获取结构120用于读取二维图像获取结构传递的每组冠状动脉二维造影图像组的图像信息,包括拍摄角度和探测距离;筛选感兴趣图像结构130用于根据探测距离,分别从每组冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。
本申请提供了一种冠状动脉分析系统,包括:上述的三维血管合成系统。
本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述的三维血管的合成方法。
所属技术领域的技术人员知道,本发明的各个方面可以实现为系统、方法或计算机程序产品。因此,本发明的各个方面可以具体实现为以下形式,即:完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等),或硬件和软件方面结合的实施方式,这里可以统称为“电路”、“模块”或“系统”。此外,在一些实施例中,本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式,该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或系统的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。
例如,可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件,可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作系统执行的多个软件指令。在本发明的示例性实施例中,由数据处理器来执行如本文的根据方法和/或系统的示例性实施例的一个或多个任务,诸如用于执行多个指令的计算平台。可选地,该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器,例如,磁硬盘和/或可移动介质。可选地,也提供了一种网络连接。可选地也提供显示器和/或用户输入设备,诸如键盘或鼠标。
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项:
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括(但不限于)无线、有线、光缆、RF等等,或者上述的任意合适的组合。
例如,可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码,包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言,诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。
也可以把这些计算机程序指令存储在计算机可读介质中,这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作,从而,存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。
还可将计算机程序指令加载到计算机(例如,冠状动脉分析系统)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程,使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。
本发明的以上的具体实例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种三维血管中心线的合成方法,其特征在于,包括:
获取至少两组拍摄角度不同的冠状动脉二维造影图像组;
读取每组所述冠状动脉二维造影图像组的图像信息,包括拍摄角度和探测距离;
根据所述探测距离,分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像;
从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线;
根据每幅所述冠状动脉二维造影图像的图像信息,将每条所述二维血管中心线投影到三维空间内,合成所述三维血管中心线的方法包括:
以心脏为坐标原点,建立三维坐标系;
获取每幅所述感兴趣的二维造影图像的左右角度α、前后角度β,以及人体与平板探测器之间的距离S,所述二维血管中心线上各点的坐标为(x,y);
将各点(x,y)投影到三维空间内,获得一系列三维坐标点P,坐标为(x’、y’、z’);
将放射源投影到所述三维空间内形成放射点R;
将每个所述三维坐标点P均与所述放射点R进行连线,从PR连线上获得所述三维血管中心线上的点,将所述三维血管中心线上的点依次连接获得所述三维血管中心线;
所述将各点(x,y)投影到三维空间内,获得一系列三维坐标点P,坐标为(x’、y’、z’)的方法,包括:
将各点(x,y)绕y轴旋转获得(x’、y’、z’)系列点,具体公式为:
将(x’、y’、z’)系列点绕x轴旋转获得一系列三维坐标点P,坐标为(x’、y’、z’),具体公式为:
2.根据权利要求1所述的一种三维血管中心线的合成方法,其特征在于,所述将放射源投影到所述三维空间内形成放射点R的方法,包括:
获得每幅所述感兴趣的二维造影图像中人体与放射源之间的距离S’;
根据公式
其中,R在三维空间内的坐标为(a,b,c)。
3.根据权利要求2所述的一种三维血管中心线的合成方法,其特征在于,所述将每个所述三维坐标点P均与所述放射点R进行连线,从PR连线上获得所述三维血管中心线上的点,将所述三维血管中心线上的点依次连接获得所述三维血管中心线的方法,包括:
将从同一幅所述感兴趣的二维造影图像中获得的一系列三维坐标点P和放射点R对应连接,获得多条PR直线;
获取所述血管的同一位置的两条PR直线之间的最小距离的点,分别为A点和B点;
将所述A点与所述B点连线,获取AB线段的中点作为所述三维血管的中心线上的点;
将获得的一系列的所述三维血管的中心线上的点依次连接,获得所述三维血管中心线。
4.根据权利要求1所述的一种三维血管中心线的合成方法,其特征在于,所述分别从每幅所述冠状动脉二维造影图像中提取一条二维血管中心线的方法包括:
读取冠状动脉二维造影图像;
获取感兴趣的血管段;
拾取所述感兴趣的血管段的起始点、种子点和结束点;
分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割,得到至少两个局部血管区域图;
从每个所述局部血管区域图中提取至少一条血管局部路径线;
将每个所述局部血管区域图上相对应的血管局部路径线连接,获得至少一条所述血管路径线;
选取一条所述血管路径线作为所述二维血管中心线。
5.根据权利要求4所述的一种三维血管中心线的合成方法,其特征在于,所述从每个所述局部血管区域图上提取至少一条血管局部路径线的方法包括:
对所述局部血管区域图做图像增强处理,得到对比强烈的粗略血管图;
对所述粗略血管图做网格划分,沿着所述起始点至所述结束点方向,提取至少一条血管局部路径线。
6.根据权利要求5所述的一种三维血管中心线的合成方法,其特征在于,所述对所述局部血管区域图做图像增强处理,得到对比强烈的粗略血管图的方法,包括:
在每幅所述局部血管区域图中,以所述感兴趣的血管段作为前景,其他区域作为背景,强化所述前景,弱化所述背景,得到对比强烈的所述粗略血管图。
7.根据权利要求6所述的一种三维血管中心线的合成方法,其特征在于,所述对所述粗略血管图做网格划分,沿着所述起始点至所述结束点方向,提取至少一条血管局部路径线的方法包括:
对所述粗略血管图进行网格划分;
沿着所述起始点至所述结束点的血管延伸方向,搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点,搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点,所述第三个点重复上述步骤,直至最短时间路径到达结束点,其中,n为大于等于1的正整数;
按照搜索顺序,从所述起始点至所述结束点的血管延伸方向连线,获得至少一条血管局部路径线。
8.根据权利要求7所述的一种三维血管中心线的合成方法,其特征在于,所述选取一条所述血管路径线作为所述二维血管中心线的方法包括:
如果血管路径线为两条或两条以上,则对每条血管路径线从所述起始点至所述结束点所用的时间求和;
取用时最少的所述血管路径线作为所述二维血管中心线。
9.一种三维血管中心线合成系统,用于权利要求1~8任一项所述的三维血管中心线的合成方法,其特征在于,包括:图像读取装置,与所述图像读取装置连接的二维血管中心线提取装置和所述三维血管中心线获取装置,所述二维血管中心线提取装置与所述三维血管中心线获取装置连接;
所述图像读取装置,用于获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息;
所述二维血管中心线提取装置,用于接收所述图像读取装置发送的冠状动脉二维造影图像,从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线;
所述三维血管中心线获取装置,用于接收所述二维血管中心线提取装置发送的二维血管中心线,以及接收所述图像读取装置发送的冠状动脉二维造影图像和图像信息,根据每幅所述冠状动脉二维造影图像的图像信息,将每条所述二维血管中心线投影到三维空间内,合成所述三维血管中心线。
10.根据权利要求9所述的一种三维血管中心线合成系统,其特征在于,所述图像读取装置包括:依次连接的二维图像获取结构、图像信息获取结构、筛选感兴趣图像结构,所述筛选感兴趣图像结构与所述二维图像获取结构连接;
所述二维图像获取结构,用于获取至少两组拍摄角度不同的冠状动脉二维造影图像组;
所述图像信息获取结构,用于读取所述二维图像获取结构传递的每组所述冠状动脉二维造影图像组的图像信息,包括拍摄角度和探测距离;
所述筛选感兴趣图像结构,用于根据所述探测距离,分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。
11.一种冠状动脉分析系统,其特征在于,包括:权利要求9~10任一项所述的一种三维血管中心线合成系统。
12.一种计算机存储介质,其特征在于,计算机程序被处理器执行时实现权利要求1~8任一项所述的一种三维血管中心线的合成方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011337612.1A CN112419277B (zh) | 2020-11-25 | 2020-11-25 | 三维血管中心线合成方法、系统及存储介质 |
PCT/CN2020/131697 WO2022109902A1 (zh) | 2020-11-25 | 2020-11-26 | 三维血管中心线合成方法、系统及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011337612.1A CN112419277B (zh) | 2020-11-25 | 2020-11-25 | 三维血管中心线合成方法、系统及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112419277A CN112419277A (zh) | 2021-02-26 |
CN112419277B true CN112419277B (zh) | 2024-02-23 |
Family
ID=74843327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011337612.1A Active CN112419277B (zh) | 2020-11-25 | 2020-11-25 | 三维血管中心线合成方法、系统及存储介质 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112419277B (zh) |
WO (1) | WO2022109902A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115965750B (zh) * | 2023-03-16 | 2023-06-16 | 深圳微创踪影医疗装备有限公司 | 血管重建方法、装置、计算机设备和可读存储介质 |
CN116863146B (zh) * | 2023-06-09 | 2024-03-08 | 强联智创(北京)科技有限公司 | 用于对血管瘤特征进行提取的方法、设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH105203A (ja) * | 1996-06-21 | 1998-01-13 | Toshiba Corp | 医用診断システム,医用診断情報生成方法及び3次元画像再構成方法 |
US6047080A (en) * | 1996-06-19 | 2000-04-04 | Arch Development Corporation | Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images |
CN110287956A (zh) * | 2019-06-13 | 2019-09-27 | 北京理工大学 | 血管中心线自动匹配方法及装置 |
CN110889896A (zh) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | 获取血管狭窄病变区间及三维合成方法、装置和系统 |
CN111161342A (zh) * | 2019-12-09 | 2020-05-15 | 杭州脉流科技有限公司 | 基于冠脉造影图像获取血流储备分数的方法、装置、设备、系统及可读存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108133478B (zh) * | 2018-01-11 | 2021-11-23 | 苏州润迈德医疗科技有限公司 | 一种提取冠状动脉血管中心线的方法 |
CN108186038B (zh) * | 2018-02-11 | 2020-11-17 | 杭州脉流科技有限公司 | 基于动脉造影影像计算冠脉血流储备分数的系统 |
-
2020
- 2020-11-25 CN CN202011337612.1A patent/CN112419277B/zh active Active
- 2020-11-26 WO PCT/CN2020/131697 patent/WO2022109902A1/zh active Application Filing
Patent Citations (5)
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 |
JPH105203A (ja) * | 1996-06-21 | 1998-01-13 | Toshiba Corp | 医用診断システム,医用診断情報生成方法及び3次元画像再構成方法 |
CN110287956A (zh) * | 2019-06-13 | 2019-09-27 | 北京理工大学 | 血管中心线自动匹配方法及装置 |
CN110889896A (zh) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | 获取血管狭窄病变区间及三维合成方法、装置和系统 |
CN111161342A (zh) * | 2019-12-09 | 2020-05-15 | 杭州脉流科技有限公司 | 基于冠脉造影图像获取血流储备分数的方法、装置、设备、系统及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112419277A (zh) | 2021-02-26 |
WO2022109902A1 (zh) | 2022-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112164020B (zh) | 精确提取血管中心线的方法、装置、分析系统和存储介质 | |
CN110889896B (zh) | 获取血管狭窄病变区间及三维合成方法、装置和系统 | |
CN112419484B (zh) | 三维血管合成方法、系统及冠状动脉分析系统和存储介质 | |
JP5129480B2 (ja) | 管状臓器の3次元再構成を行うシステム及び血管撮像装置の作動方法 | |
US20220277447A1 (en) | Method and apparatus for acquiring contour line of blood vessel according to centerline of blood vessel | |
US20130064343A1 (en) | Method and Apparatus for Determining Three-Dimensional Reconstruction of an Object | |
CN112132949B (zh) | 具有狭窄病变区间的血管三维建模方法、装置和系统 | |
CN112419462B (zh) | 三维血管的渲染合成方法、系统及存储介质 | |
CN112419277B (zh) | 三维血管中心线合成方法、系统及存储介质 | |
CN112446866B (zh) | 血流参数的计算方法、装置、设备及存储介质 | |
CN1875886B (zh) | 医用图像显示装置和系统、以及x射线计算机断层摄影装置 | |
CN112419276B (zh) | 调节血管轮廓及中心线的方法及存储介质 | |
WO2021087961A1 (zh) | 测量舒张期血流速度的方法、装置、系统及存储介质 | |
CN112446867A (zh) | 血流参数的确定方法、装置、设备及存储介质 | |
JP2014174654A (ja) | 核医学骨画像解析技術 | |
US20050185831A1 (en) | Analysis of a multi-dimensional structure | |
CN108805815B (zh) | 基于x射线血管造影图像的血管拉直重建方法 | |
CN112419279B (zh) | 二维图像选取及三维血管合成的方法和存储介质 | |
CN112419280B (zh) | 精确获取狭窄病变区间的方法及存储介质 | |
US7116808B2 (en) | Method for producing an image sequence from volume datasets | |
WO2022000730A1 (zh) | 基于ct序列图像获取心脏重心的方法和系统 | |
US10977792B2 (en) | Quantitative evaluation of time-varying data | |
JP6386629B2 (ja) | 核医学骨画像解析技術 | |
JP2022059493A (ja) | モデル生成方法、モデル生成装置、画像処理方法及び画像処理装置 | |
CN112487342B (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 |