基于CTA图像的动态心脏各腔室分割方法技术

本发明专利技术公开了一种基于CTA图像的动态心脏各腔室分割方法，包括以下步骤：S1，获取心脏CTA数据和CVH数据集；S2，针对心脏CTA图像，获得CT值的灰度直方图，利用阈值法获取整个心脏血液分布的分割数据，实现心脏血液分布的最大连通区域；S3，截断主动脉根部与左心室的连通关系，实现冠状动脉的分割；S4，利用CVH数据集获得的心脏纤维骨架三维模型实现对心脏各个腔室的空间划分；S5，将心脏纤维骨架与心脏CTA数据的主动脉根部进行配准，获得CVH空间到心脏CTA图像的变换矩阵；S6，在CVH数据集上定义好瓣环的空间坐标，然后将变换矩阵应用到CTA数据集上，实现心脏各个腔室的体数据分割；S7，对心脏CTA数据至少两个时相断层影像重复S1到S6，实现动态心脏的影像分割。

【技术实现步骤摘要】
基于CTA图像的动态心脏各腔室分割方法
本专利技术涉及医学图像处理
，具体涉及一种基于CTA图像的动态心脏各腔室分割方法。
技术介绍
心脏是人体心血管系统的关键器官，内部组织结构较复杂，活体处于不断的舒缩中，而目前心脏的影像资料多为正交的二维平面影像像，每个平面都是在某个时刻的心脏运动过程中的某个采样，临床医师较难通过二维图像去想象病变部位的三维形态和空间毗邻关系，这就决定了心脏疾病在影像学上定位和定性诊断较困难。中国可视化人(CVH)项目，于2002年10月，由第三军医大学发起，该项目是从尸体获得包含丰富解剖学信息的高分辨率、连续断层切片图像的数据集。此数据集的优点是，它不仅可以清楚地观察到的心脏纤维骨架(该数据集中连续断层图像中各个组织呈自然的纹理颜色，容易区分不同组织的边界)，而且保持了解剖结构的原始空间关系。CVH数据集层间距小于临床连续断层影像。水平分辨率越大，层间距越小，能够呈现的解剖细节就越丰富，特别是二尖瓣、三尖瓣等结构需要在心脏轴面来观察的结构。通过数字化人体数据集建立高精度的心脏可视化解剖模型可以反应心脏毗邻解剖结构原位的空间关系[1]，但是如何将这些信息用于活体的断层解剖图像，以指导临床的诊断和资料，这是心脏可视化面临的一个挑战。借助医学影像技术，定性和定量分析心脏的解剖及运动规律，进而研究其与心脏疾病之间的关系是现代医学研究的热点。利用计算机模型对形态进行仿真，模拟真实心脏的运动过程及病理状态，进一步认识心血管系统的运动规律和本质，深入了解心血管系统疾病和两次表现之间的关系。关键问题是形变模型要适应心脏组织的周期性变化和不同的心脏个体，即以模型为引导的分割，用模型去适应图像特征的变化。心脏CTA(CardiacComputedTomographyArteriography，计算机断层造影术)以成像速度快、分辨率高、信息量大，作为心脏检查的重要手段。基于CTA图像的量化评估对心功能的评估上具有重要的作用，特别是在动态定量评估方面能通过较高的时间分辨率捕捉心动周期内心脏内部形态变化特征，进而评估心脏生理功能。因此，基于CTA的心脏量化评估引起国内外学者的高度关注。定量评估的核心是对心脏区域的分割。心脏分割的目的主要是从复杂的心脏图像中提取出心脏感兴趣区域的轮廓曲线，然而由于医学图像的复杂性与多样性，成像质量常受到噪声、运动伪影、局部体效应的干扰，在很多场合下还是会选择人工手动分割的方式进行处理[2]。目前国内外学者提出了许多针对心脏中某些解剖结构的分割算法，其主要关注点在于心室的分割[3]。由于心脏形态结构的复杂性，心脏分割的准确性和实时性一直是研究的热点。在CTA图像中,单一心腔的分割相对比较容易特别是左心室，但是对于四个腔室的同时分割比较困难。其原因是心腔在解剖结构上具有特殊性，使医学图像呈现出一些不容易分割的状态：1)单纯依靠图像特征难以完全实现解剖结构的分割。心腔的CTA图像中的灰度变化与腔体的分隔线并非完全对应，某些在灰度上连通的区域在心脏解剖结构上是分离的，而心肌与周围组织的灰度近似，仅依靠灰度难以区分；2)心脏一直处于运动状态，成像伪影及噪声会影响成像效果；3)由于心脏结构的特殊性，心脏图像中会出现曲率很大的凹陷区域，这部分图像难以分割[4]。因此，如何借助解剖知识实现心脏的分割是进行心脏各腔室量化评估的关键问题。心脏纤维骨架在心脏运动中起支点和稳定的作用，心脏纤维骨架既是几何形态的中心，又是位于各方向运动的中心，其“双中心”的地位对理解心脏运动规律有极其重要的意义。超声能实时捕捉到瓣膜的运动但是对周围解剖结构的显示却不理想，CT、MRI能显示更详细的解剖结构，但是重建出瓣膜的三维形态非常困难，关键在于其形态走行复杂从无法准确判断二维的影像上起止点，特别是在心肌附着处难以区分其边界。因此，建立人体心脏原位的心脏纤维骨架的三维模型是非常困难的。人们尝试着使用动物和使用不同的技术，包括解剖或影像学检查人的研究，已有学者研究了二尖瓣环和三尖瓣环的形态学特征及动态特征[5]，然而已有的研究中少有对同一心脏四个瓣膜的情况进行分析的。对当前的成像后处理方式而言，同时定位和识别4个瓣膜仍然面临困难。数字化人体是高精度的彩色纹理连续断层图像，对心腔内部精细结构的鉴别和区分具有较大优势，在此基础上建立的数字化模型维持了心脏内部毗邻解剖结构的原位空间关系，组织间没有相对的位移和形变，便于和临床影像有机结合。可以将已建立的详细的解剖结构模型结合CTA、MRI等多种模式下的可视化方式，从而全面了解心脏在生理、病理下的状态。1.郭燕丽,人体心脏薄层断面解剖和三维可视化与多平面TEE的对照研究.2003,第三军医大学.2.付增良,陈晓军,叶铭,李峰,王成焘,心脏CT图像分割方法.计算机工程,2009.35(12):p.189-191.3.Zhou,Y.,W.R.Shi,W.Chen,Y.L.Chen,Y.Li,etal.,Activecontoursdrivenbylocalizingregionandedge-basedintensityfittingenergywithapplicationtosegmentationoftheleftventricleincardiacCTimages.Neurocomputing,2015.156(C):p.199-210.4.贺娟,基于CT图像的全心脏分割算法研究.2015,浙江大学.5.Maffessanti,F.,P.Gripari,G.Pontone,D.Andreini,E.Bertella,etal.,Three-dimensionaldynamicassessmentoftricuspidandmitralannuliusingcardiovascularmagneticresonance.EuropeanHeartJournalCardiovascularImaging,2013.14(10):p.986-995.
技术实现思路
鉴于以上技术问题，本专利技术提出了一种基于CTA图像的动态心脏各腔室分割方法，该专利技术通过结合解剖学与生理学知识，用CVH高分辨率切片重建的心脏纤维骨架三维模型确定其在心脏三维空间中的位置，实现连续影像数据心脏平面的快速构建，以解决对心脏四个腔室同时进行分割及不同个体CTA图像在可视化量化分析的一致性问题。本专利技术是通过如下技术方案来实现的：一种基于CTA图像的动态心脏各腔室分割方法，包括以下步骤：S1，获取心脏CTA数据和CVH数据集；S2，针对获得的心脏CTA图像，获得CT值的灰度直方图，针对该直方图设定一个阈值，利用阈值法获取整个心脏血液分布的分割数据，实现心脏血液分布的最大连通；S3，截断主动脉根部与左心室的连通关系，实现冠状动脉的分割；S4，利用CVH数据集获得的心脏纤维骨架三维模型实现对心脏各个腔室的空间进行划分；S5，将心脏纤维骨架与心脏CTA数据的主动脉根部之间进行位置、大小、角度等的变换，实现两个图像的配准，获得CVH空间到心脏CTA图像的几何变换矩阵；S6，在CVH数据集上定义好四个瓣环的空间坐标，然后将上述步骤S5获得的几何变换矩阵应用到CTA数据集上，获得定义好的各个腔室的空间坐标，并本文档来自技高网...

【技术保护点】
一种基于CTA图像的动态心脏各腔室分割方法，其特征在于，包括以下步骤：S1，获取心脏CTA数据和CVH数据集；S2，针对获得的心脏CTA图像，获得CT值的灰度直方图，针对该直方图设定一个阈值，利用阈值法获取整个心脏血液分布的分割数据，实现心脏血液分布的最大连通；S3，截断主动脉根部与左心室的连通关系，实现冠状动脉的分割；S4，利用CVH数据集获得的心脏纤维骨架三维模型实现对心脏各个腔室的空间进行划分；S5，将心脏纤维骨架与心脏CTA数据的主动脉根部之间进行位置、大小、角度等的变换，实现两个图像的配准，获得CVH空间到心脏CTA图像的几何变换矩阵；S6，在CVH数据集上定义好四个瓣环的空间坐标，然后将上述步骤S5获得的几何变换矩阵应用到CTA数据集上，获得定义好的各个腔室的空间坐标，并通过三维数据集上的重采样重建出具有解剖细节的CTA数据集的四个瓣环平面图像，实现心脏各个腔室的体数据分割；S7，对心脏CTA数据至少两个时相断层影像的每个时相重复步骤S1到步骤S6，实现动态心脏的影像分割。

【技术特征摘要】
1.一种基于CTA图像的动态心脏各腔室分割方法，其特征在于，包括以下步骤：S1，获取心脏CTA数据和CVH数据集；S2，针对获得的心脏CTA图像，获得CT值的灰度直方图，针对该直方图设定一个阈值，利用阈值法获取整个心脏血液分布的分割数据，实现心脏血液分布的最大连通；S3，截断主动脉根部与左心室的连通关系，实现冠状动脉的分割；S4，利用CVH数据集获得的心脏纤维骨架三维模型实现对心脏各个腔室的空间进行划分；S5，将心脏纤维骨架与心脏CTA数据的主动脉根部之间进行位置、大小、角度等的变换，实现两个图像的配准，获得CVH空间到心脏CTA图像的几何变换矩阵；S6，在CVH数据集上定义好四个瓣环的空间坐标，然后将上述步骤S5获得的几何变换矩阵应用到CTA数据集上，获得定义好的各个腔室的空间坐标，并通过三维数据集上的重采样重建出具有解剖细节的CTA数据集的四个瓣环平面图像，实现心脏各个腔室的体数据分割；S7，对心脏CTA数据至少两个时相断层影像的每个时相重复步骤S1到步骤S6，实现动态心脏的影像分割。2.根据权利要求1所述的基于CTA图像的动态心脏各腔室分割方法，其特征在于：所述步骤S3具体为，截断主动脉根部和左心室的连接部分，对主动脉根部使用连通域的方法，分割出主动脉根部图像，获得冠状动脉的分割图像。3.根据权利要求1所述的基于CTA图像的动态心脏各腔室分割方法，其特征在于：所述步骤S4具体为，针对CVH数据集，人工进行心脏纤维骨架解剖结构的分割，获得心脏纤维骨架的三维模型，实现对心脏各个腔室空间的划分。4.根据权利要求1所述的基于CTA图像的动态心脏各腔室分割方法，其特征在于：所述步骤S5具体包括以下步骤，S51,由于心脏CTA成像的结果与心脏纤维骨架数据重建模型呈三维的镜像关系，所以针对心脏CTA数据，首先需要通过镜像变换调整主动脉根部图像的空间位置，使之与心脏纤维骨架的位置相对应；S52,以主动脉瓣环三皇冠三个顶点为准拟合得到的圆环(环1)与心脏纤维骨架环主动脉瓣环三皇冠三个底点拟合圆环(环2)之间的半径比作为图像尺度变换的参数；S53,以环1所在平面与环2所在平面之间的夹角作为旋转变换的参数，以此来消除由于个体差异和心脏运动状态的不同导致的角度变化；S54,通过平移变换实现对整个冠状动脉和纤维骨架空间位置的对应...

【专利技术属性】
技术研发人员：李颖，张绍祥，谭立文，席阳，张小勤，
申请(专利权)人：中国人民解放军第三军医大学，
类型：发明
国别省市：重庆,50