[发明专利]一种CT去病床方法有效
申请号: | 202110446942.2 | 申请日: | 2021-04-25 |
公开(公告)号: | CN113177963B | 公开(公告)日: | 2022-05-20 |
发明(设计)人: | 应关荣;叶宏伟;张辉 | 申请(专利权)人: | 明峰医疗系统股份有限公司 |
主分类号: | G06T7/136 | 分类号: | G06T7/136;G06T5/00;G06T7/62;G06T7/13;G06T5/30 |
代理公司: | 浙江千克知识产权代理有限公司 33246 | 代理人: | 周希良 |
地址: | 311215 浙江省杭州市杭*** | 国省代码: | 浙江;33 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 ct 病床 方法 | ||
1.一种CT去病床方法,其特征在于,包括以下步骤:
S1、对待去病床的原始图像进行阈值分割,生成mask图像;
S2、对mask图像逐层使用孔洞填充算法进行填充,生成填充图像;
S3、对填充图像采用三维连通区域标记算法标记连通区域;
S4、丢弃像素个数小于目标个数的连通区域,并对剩余的连通区域逐个分析,转至步骤S5;
S5、计算每个连通区域的像素点个数大小,并将像素点个数最大的连通区域的像素点个数标记为Max;
S6、判断连通区域的面积是否大于Max/a,a为预设系数;若是,则转至步骤S7;若否,则转至步骤S9;
S7、提取连通区域的边界,采用连通区域的边界从原始图像中获取对应像素点的CT值,并基于CT值进行直方图统计,抓取直方图的极大值点对应的尖峰;
S8、判断尖峰数量为1或2或大于2;若为1,则转至步骤S9;若为2,则转至步骤S11;若大于2,则保留连通区域;
S9、通过连通区域,统计原始图像对应的CT值,计算其中大于骨CT值的像素点个数,得到骨像素个数;
S10、判断骨像素个数是否大于N/b,N为连通区域的像素点个数,b为预设参数;若是,则转至步骤S15;若否,则连通区域为干扰区域,丢弃;
S11、判断两个尖峰落在病床区或皮肤区;若两个尖峰均落在病床区,则连通区域为病床,丢弃;若两个尖峰均落在皮肤区,则转至步骤S15;若一个尖峰落在病床区,另一个尖峰落在皮肤区,则转至步骤S12;
S12、提取两个尖峰之间的峰谷,并将峰谷对应的CT值作为分割阈值,以分割阈值分割连通区域边界,得到病床区种子和人体区种子;
S13、利用病床区种子在连通区域区域生长出病床区域,并从连通区域中减去病床区域,获得目标区域;
S14、利用人体区种子在目标区域中区域生长出人体区域,转至步骤S15;
S15、相应的区域即为人体mask,保留;
S16、待所有连通区域分析完成之后,得到最终的CT去病床图像。
2.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S1中,阈值分割的CT值范围为-300~10000Hu。
3.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S4中,目标个数的取值为300~800。
4.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S6中,a取值为15。
5.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S7中,在直方图统计之后、抓取尖峰之前,还对直方图进行中值滤波平滑。
6.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S10中,b取值为20。
7.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S13中,生长出病床区域之后,并过滤像素点个数小于阈值的病床区域。
8.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S14中,对人体区域还进行形态学滤波、腐蚀以及膨胀处理。
9.根据权利要求1所述的一种CT去病床方法,其特征在于,所述步骤S14中,生长出人体区域之后,并过滤像素点个数小于阈值的人体区域。
10.根据权利要求1所述的一种CT去病床方法,其特征在于,所述标记连通区域过程以及区域生长过程中均采用GPU加速。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于明峰医疗系统股份有限公司,未经明峰医疗系统股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110446942.2/1.html,转载请声明来源钻瓜专利网。