一种CT影像体表提取方法及系统技术方案

本发明专利技术提出了一种CT影像体表提取方法及系统，采用粗‑精结合的体表分割方法，粗分割采用基于孤立连接的区域增长区分体表内和体表外，从而大致获得体表结构，细分割过程针对粗分割三维体表的轴位轮廓进行修复和细节保持平滑处理，从而实现了三维体表的快速、精确和鲁棒性分割提取。

【技术实现步骤摘要】

本专利技术涉及医学影像处理
，特别涉及一种CT影像体表提取方法及系统。
技术介绍
体表分割属于图像分割的范畴，常用的图像分割方法主要包括以下几类：基于区域的方法，基于边缘的方法，基于分类器、聚类、统计学的方法。依据待分割对象的不同特性，所采用或设计的分割方法也千差外别。针对体表的分割，在理想情况下使用边界跟踪或合理的阈值分割等方法就可以得到体表区域的边缘。但是临床的医学影像中，由于人体内外结构的异质性、CT成像过程的畸变或噪声的影响及其它不确定干扰因素，使得实际处理的医学影像体表往往边缘模糊、信号强度不均匀甚至噪声较大。另外，人体体表的结构比较复杂而且个体外部形态和生理差异较大，这也给医学影像的体表分割带来了很大的困难。Manning Wang和Zhijian Song在Pattern Recognition Letters提出了采用全自动或手动确定阈值的分割方法提取体表结构。Giovanni Fattori等在IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING采用了基于Marching Cube分割的体表提取方法，主要选择体表的信号强度值本文档来自技高网...

【技术保护点】
一种CT影像体表提取方法，其特征在于，包括下述步骤：获取三维CT影像；对所述CT影像体表进行体表粗分割；对所述CT影像体表进行体表细分割；输出经细分割后的CT影像体表影像；其中，体表粗分割包括下述步骤：分别选取所述三维CT影像中体表的边缘和外部的两个种子点，分别记为第一种子点和第二种子点；基于所述第一种子点和第二种子点采用孤立连接区域增长的方法区分体表内和体表外，获取体表的三维二值图像；其中，体表细分割包括下述步骤：采用边缘跟踪算法对所述三维二值图像的轴位影像分别逐层追踪最外边缘的轮廓，获得每一层切片的边缘轮廓坐标；遍历所述三维二值图像所有层面的轴位影像进行逐层轮廓跟踪提取，获得所有层面的轮廓...

【技术特征摘要】
1.一种CT影像体表提取方法，其特征在于，包括下述步骤：获取三维CT影像；对所述CT影像体表进行体表粗分割；对所述CT影像体表进行体表细分割；输出经细分割后的CT影像体表影像；其中，体表粗分割包括下述步骤：分别选取所述三维CT影像中体表的边缘和外部的两个种子点，分别记为第一种子点和第二种子点；基于所述第一种子点和第二种子点采用孤立连接区域增长的方法区分体表内和体表外，获取体表的三维二值图像；其中，体表细分割包括下述步骤：采用边缘跟踪算法对所述三维二值图像的轴位影像分别逐层追踪最外边缘的轮廓，获得每一层切片的边缘轮廓坐标；遍历所述三维二值图像所有层面的轴位影像进行逐层轮廓跟踪提取，获得所有层面的轮廓坐标集合。2.如权利要求1所述的CT影像体表提取方法，其特征在于，所述三维CT影像为胸腹部三维CT影像。3.如权利要求2所述的CT影像体表提取方法，其特征在于，其中，体表细分割中，边缘跟踪算法为8连通区域的轮廓跟踪算法。4.如权利要求3所述的CT影像体表提取方法，其特征在于，所述8连通区域的轮廓跟踪算法，包括下述步骤：步骤S110：从上到下，从左到右的顺序扫描所述三维二值图像的轴位影像，寻找最外围的第一个边界起始点，记为A0，其中，A0是具有最小行或列值的边缘点；步骤S120：定义一个扫描方向变量dir,该变量用于记录上一步中沿着前一个边界点到当前边界点的移动方向，8连通区域时其初始化时取值为dir＝7；步骤S130：按逆时针方向搜索当前像素的3×3邻域，其起始搜索方向设定如下：对8连通区域，若dir为奇数取(dir+7)取余8；若dir为偶数为(dir+6)取余8；步骤S140：在3×3邻域中搜索到的第一个与当前像素值相同的像素...

【专利技术属性】
技术研发人员：周武，张丽娟，梁长虹，黄美萍，庄建，
申请(专利权)人：中国科学院深圳先进技术研究院，广东省心血管病研究所，
类型：发明
国别省市：广东;44