[发明专利]基于超体素和图割算法的三维肝脏CT图像自动分割方法有效

专利信息
申请号: 201510173307.6 申请日: 2015-04-13
公开(公告)号: CN104809723B 公开(公告)日: 2018-01-19
发明(设计)人: 吴薇薇;周著黄;吴水才;白燕萍 申请(专利权)人: 北京工业大学
主分类号: G06T7/00 分类号: G06T7/00;G06T7/11;G06T7/136;G06T7/155
代理公司: 北京思海天达知识产权代理有限公司11203 代理人: 沈波
地址: 100124 *** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 基于 超体素 算法 三维 肝脏 ct 图像 自动 分割 方法
【权利要求书】:

1.基于超体素和图割算法的三维肝脏CT图像自动分割方法,其特征在于:该方法包括以下步骤,

1.1.对输入的一例腹部CT图像I进行直方图分析,从而自适应增强图像对比度,得到增强对比度后的CT图像I′;

1.2.采用自适应阈值及形态学方法对图像I逐层进行肝脏初始轮廓分割,选取最大肝脏切片lmax,计算并对增强后图像I′提取肝脏感兴趣区域I′ROI

1.3.在最大肝脏切片lmax上根据肝脏初始轮廓选取肝脏种子点SF和背景种子点SB,采用基于高斯混合模型的方法计算肝脏和背景区域的灰度概率分布,得到高斯混合模型模型Pfg和Pbkg

1.4.对I′ROI利用SLIC算法进行过分割生成超体素;

1.5.以超体素为顶点构建无向加权图G,利用图割算法对图G做切割,得到肝脏区域分割结果二值图像Imask

1.6.利用形态学开运算、中值滤波对Imask做后处理,得到平滑的肝脏分割结果。

2.根据权利要求1所述的基于超体素和图割算法的三维肝脏CT图像自动分割方法,其特征在于:所述的步骤1.1包括,

1.1.1.分析图像直方图中峰值个数及峰值CT值,若存在2个较明显的峰值则为高对比度图像Ihigh,若仅存在1个峰值则为低对比度图像Ilow

1.1.2.计算自适应对比度拉伸的阈值范围[Imin,Imax],对于高对比度图像Ihigh有Imin=(peak1+peak2)/2-60,Imax=peak2+250;对于低对比图图像Ilow有Imin=peak1-60,Imax=peak1+250,其中peak1、peak2表示峰值CT值;

1.1.3.做对比度增强,I′=255(I-Imin)/(Imax-Imin)。

3.根据权利要求1所述的基于超体素和图割算法的三维肝脏CT图像自动分割方法,其特征在于:所述的步骤1.2包括,

1.2.1.对于每一层切片l利用自适应阈值[peak1,400Hu]方法提取肝脏可能区域lpossible

1.2.2.对lpossible依次采用形态学腐蚀、保留最大子区域、孔洞填充、形态学膨胀方法,得到肝脏区域的初始轮廓lliver

1.2.3.选取lliver中肝脏区域面积最大的作为最大肝脏切片lmax

1.2.4.根据每层lliver的二维肝脏区域边界框计算得到能包围整个三维肝脏感兴趣区域的最小边界框,从增强后图像I′中提取肝脏感兴趣区域I′ROI

4.根据权利要求1所述的基于超体素和图割算法的三维肝脏CT图像自动分割方法,其 特征在于:所述的步骤1.3包括,

1.3.1.对于最大肝脏切片lmax,在肝脏初始轮廓lliver内部按一定间隔选取肝脏种子点SF,lliver边界框以外区域按一定间隔选取背景种子点SB

1.3.2.根据种子点SF,SB计算高斯混合模型Pfg,Pbkg

其中,Ip为体素对应的灰度值,k为高斯混合模型的分量数,ω,μ和σ2分别为高斯分量的权重、均值和方差。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京工业大学,未经北京工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201510173307.6/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top