by HPC_ZY
具体步骤
- 阈值分割
- 提取人体部分
- 提取疑似肺质
- 去除非肺质
- 其他优化操作
代码实现
- 阈值分割
这里分别提供全局阈值分割、基于边缘像素分割、基于有效像素分割。但这一步我们只是要区别人体和背景,所以用最简单的全局分割就可以了。
% 全局分割
imbi0 = imbinarize(im);
% 基于轮廓像素分割
Eidx = edge(im,'log'); % 拉普拉斯获取边缘
E1 = im(Eidx); % 获取边缘像素
E1 = E1(E1>0);
imbi1 = imbinarize(im,graythresh(E1));
% 基于有效像素分割
E2 = im(im>0.02); % 获取非0像素
imbi2 = imbinarize(im