CN108242075A - 一种基于x射线冠脉造影图像的多角度血管重建方法 - Google Patents
一种基于x射线冠脉造影图像的多角度血管重建方法 Download PDFInfo
- Publication number
- CN108242075A CN108242075A CN201810010079.4A CN201810010079A CN108242075A CN 108242075 A CN108242075 A CN 108242075A CN 201810010079 A CN201810010079 A CN 201810010079A CN 108242075 A CN108242075 A CN 108242075A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- segmented
- coronary angiography
- dimensional
- angiography 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于X射线冠脉造影图像的多角度血管重建方法,包括如下步骤:S1:基于两个不同角度的X射线冠脉造影图像,分别获得两造影图像上血管的2D结构数据及特征点;S2:根据血管的特征点将其对应的造影图像上的血管进行分段,并得到每段血管的2D结构数据;S3:将两个不同角度的X射线冠脉造影图像上、且呈映射关系的两分段血管的2D结构数据进行三维重建,并得到该分段血管的3D结构数据;S4:重复步骤S3直到所有分段血管三维重建完成,再将重建后的分段血管合并,即得到完整的三维血管;本发明的优点在于,依据特征点将整段血管分割成小段,然后依据每个小段之间的对应关系分别做三维合成,从而降低了整段血管直接合成造成的误差。
Description
技术领域
本发明属于数字图像处理与医学图像的交叉领域,可应用于临床研究中X射冠脉造影图像分析,尤其是涉及一种基于X射线冠脉造影图像的多角度血管重建方法。
背景技术
目前,国内大多数医院都用的是X射线单臂造影系统对病人做X射线造影,通过旋转造影臂得到一个对应于不同造影角度的造影图像序列。单臂造影可以很方便对病人进行不同角度的造影,但是对重建来说个缺点就是我们无法得到同一时刻的不同视角的造影图像,这给重建带来很大的困难。要重建血管的真实三维空间结构,需要得到血管至少两个不同角度的投影信息。传统方法首先提取出血管的骨架,然后通过不同视角空间约束关系,对不同视角投影图像的血管像素点进行正确匹配并重建,当血管几何形变不大,几何关系明显的时候,才能大体上恢复出血管的三维空间结构。
血管树从整体上可以看作在空间中弯曲延展的管状系统,其骨架是具有树状结构的连续空间曲线,反映了血管树的整体形态特征。对于心血管双视角造影图像的匹配问题,传统的方法是通过外极线约束来找两幅造影图像中匹配的血管点对,但是外极线约束极易发生误匹配。这种方法误差较大,在血管缺失或者重叠的情况下,并不能保证双侧造影图像中过的血管段和取样点是一一对应的关系。常见的匹配方法对图像质量要求高,耗时长,在面对不同的血管的个体差异,鲁棒性差。
鉴于传统血管匹配方法中存在的问题,很多研究者都通过在序列图像中选择处于心脏运动周期中同一时刻的造影图像作为重建的参考图像,比如舒张末期。但这种方法依然解决不了这个问题:在两个视角中,同一根血管不同段的伸缩和扭曲不同。
发明内容
本发明目的是:提供一种基于X射线冠脉造影图像的多角度血管重建方法,通过分析两个角度血管的特征点,并依据特征点将整段血管分割成小段,然后依据每个小段之间的对应关系分别做三维合成,从而降低了整段血管直接合成造成的误差。
本发明的技术方案是:一种基于X射线冠脉造影图像的多角度血管重建方法,包括如下步骤:
步骤S1:基于两个不同角度的X射线冠脉造影图像,分别获得两造影图像上血管的2D结构数据及特征点;
步骤S2:根据血管的特征点将其对应的造影图像上的血管进行分段,并得到每段血管的2D结构数据;
步骤S3:将两个不同角度的X射线冠脉造影图像上、且呈映射关系的两分段血管的2D结构数据进行三维重建,并得到该分段血管的3D结构数据;
步骤S4:重复步骤S3直到所有分段血管三维重建完成,再将重建后的分段血管合并,即得到完整的三维血管。
作为优选的技术方案,步骤S1中血管的2D结构数据包括中心线、半径和角度;
作为优选的技术方案,步骤S1血管的特征点包括起始点、结束点、多个狭窄点和多个分叉点。
作为优选的技术方案,步骤S3中呈映射关系的两分段血管进行三维重建的具体方法如下:
步骤a:通过矩阵变换将两分段血管的2D结构数据分别变换到三维空间的一个正方形面上;
步骤b:将两分段血管的2D结构数据所在的正方形面与空间原点分别构成两个四棱锥体,两个四棱锥体相交的空间即为三维血管的实际限定空间;
步骤c:将两分段血管的中心线与空间原点分别构成两个3D空间曲面,两个3D空间曲面相交形成一条相交线,位于步骤b中实际限定空间内的相交线即为三维血管的中心线。
本发明的优点是:
1.本发明通过分析两个角度血管的特征点,并依据特征点将整段血管分割成小段,然后依据每个小段之间的对应关系分别做三维合成,然后再将分段血管合并,得到完整的三维血管,从而降低了整段血管直接合成造成的误差;
2.本发明解决多角度X射线冠脉造影图像的多角度血管匹配问题,为临床医学心血管疾病辅助检测提供更为有效的手段,提高了血管匹配的可靠性和精确度,从而可以有效增大血管造影三维重建的精度。
附图说明
下面结合附图及实施例对本发明作进一步描述:
图1为本发明的流程图;
图2为X射线冠脉造影图像一示意图;
图3为X射线冠脉造影图像二示意图;
图4为对应造影图像一角度的三维血管合成示意图;
图5为对应造影图像二角度的三维血管合成示意图;
图6为X角度的三维血管合成示意图;
图7为Y角度的三维血管合成示意图;
图8为Z角度的三维血管合成示意图;
图9为血管分叉点对应示意图;
图10为血管狭窄点对应示意图。
具体实施方式
实施例:参照图1所示:本发明基于X射线冠脉造影图像的多角度血管重建方法,包括如下步骤:
步骤S1:基于两个不同角度的X射线冠脉造影图像(参照图2和图3所示),分别获得两造影图像上血管的2D结构数据及特征点,2D结构数据包括中心线、半径和角度;特征点包括起始点、结束点、多个狭窄点和多个分叉点;
步骤S2:根据血管的特征点将其对应的造影图像上的血管进行分段,并得到每段血管的2D结构数据;
步骤S3:将两个不同角度的X射线冠脉造影图像上、且呈映射关系的两分段血管的2D结构数据进行三维重建,并得到该分段血管的3D结构数据;呈映射关系的两分段血管进行三维重建的具体方法如下:
步骤a:通过矩阵变换将两分段血管的2D结构数据分别变换到三维空间的一个正方形面上;
步骤b:将两分段血管的2D结构数据所在的正方形面与空间原点分别构成两个四棱锥体,两个四棱锥体相交的空间即为三维血管的实际限定空间;
步骤c:将两分段血管的中心线与空间原点分别构成两个3D空间曲面,两个3D空间曲面相交形成一条相交线,位于步骤b中实际限定空间内的相交线即为三维血管的中心线;
步骤S4:重复步骤S3直到所有分段血管三维重建完成,再将重建后的分段血管合并,即得到完整的三维血管,参照图4至图8为多角度血管合成示意图。
参照图9,该图为血管分叉点对应示意图,a1b1是平面A上的一条血管中心线,其中c1是该血管上的一个分叉点;a2b2是平面B上的一条血管中心线,其中c2是该血管上的一个分叉点;经过人工确认后,在这两条中心线上,点与点之间是相互对应的,a1对应a2,b1对应b2,c1对应c2,此时,分叉点将血管分成两段,a1c1对应a2c2,c1b1对应c2b2;同上所述在每一小段之间进行线性插值得到N个点,并将平面A和平面B上的分段血管一一对应。
参照图10,该图为血管狭窄点对应示意图,a1b1是平面A上的一条血管中心线,其中d1是该血管上的一个狭窄点;a2b2是平面B上的一条血管中心线,其中d2是该血管上的一个狭窄点;经过人工确认后,在这两条中心线上,点与点之间是相互对应的,a1对应a2,b1对应b2,d1对应d2,此时,狭窄点将血管分成两段,a1d1对应a2d2,d1b1对应d2b2;同上所述在每一小段之间进行线性插值得到N个点,并将平面A和平面B上的分段血管一一对应。
通常在一条血管上会同时出现N个(N>=1)分叉点和狭窄点,所以会将该血管分成N+1段,然后将两X射线冠脉造影图像上的每一段血管一一对应。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (4)
1.一种基于X射线冠脉造影图像的多角度血管重建方法,其特征在于,包括如下步骤:
步骤S1:基于两个不同角度的X射线冠脉造影图像,分别获得两造影图像上血管的2D结构数据及特征点;
步骤S2:根据血管的特征点将其对应的造影图像上的血管进行分段,并得到每段血管的2D结构数据;
步骤S3:将两个不同角度的X射线冠脉造影图像上、且呈映射关系的两分段血管的2D结构数据进行三维重建,并得到该分段血管的3D结构数据;
步骤S4:重复步骤S3直到所有分段血管三维重建完成,再将重建后的分段血管合并,即得到完整的三维血管。
2.根据权利要求1所述的基于X射线冠脉造影图像的多角度血管重建方法,其特征在于,步骤S1中血管的2D结构数据包括中心线、半径和角度。
3.根据权利要求1所述的基于X射线冠脉造影图像的多角度血管重建方法,其特征在于,步骤S1血管的特征点包括起始点、结束点、多个狭窄点和多个分叉点。
4.根据权利要求1所述的基于X射线冠脉造影图像的多角度血管重建方法,其特征在于,步骤S3中呈映射关系的两分段血管进行三维重建的具体方法如下:
步骤a:通过矩阵变换将两分段血管的2D结构数据分别变换到三维空间的一个正方形面上;
步骤b:将两分段血管的2D结构数据所在的正方形面与空间原点分别构成两个四棱锥体,两个四棱锥体相交的空间即为三维血管的实际限定空间;
步骤c:将两分段血管的中心线与空间原点分别构成两个3D空间曲面,两个3D空间曲面相交形成一条相交线,位于步骤b中实际限定空间内的相交线即为三维血管的中心线。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810010079.4A CN108242075A (zh) | 2018-01-05 | 2018-01-05 | 一种基于x射线冠脉造影图像的多角度血管重建方法 |
PCT/CN2018/091165 WO2019134347A1 (zh) | 2018-01-05 | 2018-06-13 | 一种基于x射线冠脉造影图像的多角度血管重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810010079.4A CN108242075A (zh) | 2018-01-05 | 2018-01-05 | 一种基于x射线冠脉造影图像的多角度血管重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108242075A true CN108242075A (zh) | 2018-07-03 |
Family
ID=62699384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810010079.4A Pending CN108242075A (zh) | 2018-01-05 | 2018-01-05 | 一种基于x射线冠脉造影图像的多角度血管重建方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN108242075A (zh) |
WO (1) | WO2019134347A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109872336A (zh) * | 2019-03-13 | 2019-06-11 | 数坤(北京)网络科技有限公司 | 一种血管分割方法、设备及计算机存储介质 |
CN111166316A (zh) * | 2018-11-13 | 2020-05-19 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法 |
CN111166315A (zh) * | 2018-11-13 | 2020-05-19 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算瞬时无波型比率和静息态舒张期压力比率的方法 |
CN112132882A (zh) * | 2019-11-19 | 2020-12-25 | 苏州润迈德医疗科技有限公司 | 从冠状动脉二维造影图像中提取血管中心线的方法和装置 |
CN112419279A (zh) * | 2020-11-25 | 2021-02-26 | 苏州润迈德医疗科技有限公司 | 二维图像选取及三维血管合成的方法和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763642A (zh) * | 2009-12-31 | 2010-06-30 | 华中科技大学 | 一种用于冠脉造影三维重建的匹配方法 |
CN105894445A (zh) * | 2016-03-31 | 2016-08-24 | 北京思创贯宇科技开发有限公司 | 一种冠脉图像处理方法和装置 |
CN105877767A (zh) * | 2016-03-31 | 2016-08-24 | 北京思创贯宇科技开发有限公司 | 一种冠脉图像造影方法和装置 |
CN107392994A (zh) * | 2017-06-30 | 2017-11-24 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5575497B2 (ja) * | 2009-02-04 | 2014-08-20 | 株式会社東芝 | X線診断装置及び画像処理装置 |
CN101799935B (zh) * | 2009-12-31 | 2011-09-14 | 华中科技大学 | 一种单臂x射线造影图像的动态三维重建方法 |
-
2018
- 2018-01-05 CN CN201810010079.4A patent/CN108242075A/zh active Pending
- 2018-06-13 WO PCT/CN2018/091165 patent/WO2019134347A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763642A (zh) * | 2009-12-31 | 2010-06-30 | 华中科技大学 | 一种用于冠脉造影三维重建的匹配方法 |
CN105894445A (zh) * | 2016-03-31 | 2016-08-24 | 北京思创贯宇科技开发有限公司 | 一种冠脉图像处理方法和装置 |
CN105877767A (zh) * | 2016-03-31 | 2016-08-24 | 北京思创贯宇科技开发有限公司 | 一种冠脉图像造影方法和装置 |
CN107392994A (zh) * | 2017-06-30 | 2017-11-24 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111166316A (zh) * | 2018-11-13 | 2020-05-19 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法 |
CN111166315A (zh) * | 2018-11-13 | 2020-05-19 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算瞬时无波型比率和静息态舒张期压力比率的方法 |
WO2020098140A1 (zh) * | 2018-11-13 | 2020-05-22 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法 |
CN111166316B (zh) * | 2018-11-13 | 2023-03-21 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法 |
CN111166315B (zh) * | 2018-11-13 | 2023-03-28 | 苏州润迈德医疗科技有限公司 | 基于造影图像计算瞬时无波型比率和静息态舒张期压力比率的方法 |
CN109872336A (zh) * | 2019-03-13 | 2019-06-11 | 数坤(北京)网络科技有限公司 | 一种血管分割方法、设备及计算机存储介质 |
CN112132882A (zh) * | 2019-11-19 | 2020-12-25 | 苏州润迈德医疗科技有限公司 | 从冠状动脉二维造影图像中提取血管中心线的方法和装置 |
CN112419279A (zh) * | 2020-11-25 | 2021-02-26 | 苏州润迈德医疗科技有限公司 | 二维图像选取及三维血管合成的方法和存储介质 |
WO2022109904A1 (zh) * | 2020-11-25 | 2022-06-02 | 苏州润迈德医疗科技有限公司 | 二维图像选取及三维血管合成的方法和存储介质 |
CN112419279B (zh) * | 2020-11-25 | 2024-02-23 | 苏州润迈德医疗科技有限公司 | 二维图像选取及三维血管合成的方法和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
WO2019134347A1 (zh) | 2019-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108242075A (zh) | 一种基于x射线冠脉造影图像的多角度血管重建方法 | |
Torosdagli et al. | Deep geodesic learning for segmentation and anatomical landmarking | |
Nelson et al. | Visualization of 3D ultrasound data | |
JP3570576B2 (ja) | マルチモダリティに対応した3次元画像合成表示装置 | |
US9980691B2 (en) | Method and apparatus for three dimensional viewing of images | |
US9473766B2 (en) | Method and apparatus for three dimensional viewing of images | |
US4737921A (en) | Three dimensional medical image display system | |
Cong et al. | Quantitative analysis of deformable model-based 3-D reconstruction of coronary artery from multiple angiograms | |
US9349183B1 (en) | Method and apparatus for three dimensional viewing of images | |
Prause et al. | Binary reconstruction of the heart chambers from biplane angiographic image sequences | |
US20050080328A1 (en) | Method and apparatus for medical intervention procedure planning and location and navigation of an intervention tool | |
CN100361626C (zh) | 分布式心脏图像多维重建与交互可视化方法 | |
US20110007138A1 (en) | Global camera path optimization | |
CN103561655B (zh) | 医用图像诊断装置、医用图像处理装置以及方法 | |
CN101006933A (zh) | 用于显示三维对象的方法和设备 | |
CN111009032B (zh) | 基于改进外极线约束匹配的血管三维重建方法 | |
CN111815766A (zh) | 基于2d-dsa图像重建血管三维模型处理方法及系统 | |
CN115619790B (zh) | 一种基于双目定位的混合透视方法、系统及设备 | |
WO2023142781A1 (zh) | 图像三维重建方法、装置、电子设备及存储介质 | |
Wu et al. | Reconstructing 3D lung shape from a single 2D image during the deaeration deformation process using model-based data augmentation | |
Ropinski et al. | Multimodal vessel visualization of mouse aorta PET/CT scans | |
US20110050692A1 (en) | Interpolating and rendering sub-phases of a 4d dataset | |
CN102892015A (zh) | 图像处理装置、图像处理方法以及医用图像诊断装置 | |
US20060103678A1 (en) | Method and system for interactive visualization of locally oriented structures | |
CN113256754B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180703 |
|
RJ01 | Rejection of invention patent application after publication |