%??（圖像分割）MATLAB胸部CT肺質提取
%??demo
%??by?HPC_ZY?20190718
clear;?close?all;?clc

load?imset?%?數據集
showflag?=?0;?%?是否顯示中間圖像
imidx?=?1;?%?測試圖像序號

im?=?im2double（imset{imidx}）;

%%?1.閾值分割
%?全局分割
imbi0?=?imbinarize（im）;

%?基于輪廓像素分割
Eidx?=?edge（im‘log‘）;?%?拉普拉斯獲取邊緣
E1?=?im（Eidx）;?%?獲取邊緣像素
E1?=?E1（E1>0）;
imbi1?=?imbinarize（imgraythresh（E1））;

%?基于有效像素分割
E2?=?im（im>0.02）;?%?獲取非0像素
imbi2?=?imbinarize（imgraythresh（E2））;

if?showflag
????figure
????subplot（221）imshow（im）title（‘\fontsize{16}原圖‘）
????subplot（222）imshow（imbi0）title（‘\fontsize{16}全局分割‘）
????subplot（223）imshow（imbi1）title（‘\fontsize{16}基于輪廓像素分割‘）
????subplot（224）imshow（imbi2）title（‘\fontsize{16}基于有效像素分割‘）
end


%%?2.提取人體部分
%?計算連通分量
[labelnum]?=?bwlabel（imbi0）;
%?計算最大連通分量
MAX?=?0;
for?k?=?1:num
????maxtmp?=?sum（find（label==k））;
????if?maxtmp>MAX
????????IDX?=?k;
????????MAX?=?maxtmp;
????end
end
imbi?=?label==IDX;

if?showflag
????figure
????subplot（121）imshow（imbi0）title（‘\fontsize{16}二值圖‘）
????subplot（122）imshow（imbi）title（‘\fontsize{16}胸腔‘）
end


%%?3.提取疑似肺質
imbiFull?=?imfill（imbi‘hole‘）;?%?填充
objtmp?=?imbiFull-imbi;

if?showflag
????figure
????subplot（131）imshow（imbi）title（‘\fontsize{16}胸腔‘）
????subplot（132）imshow（imbiFull）title（‘\fontsize{16}填充‘）
????subplot（133）imshow（objtmp）title（‘\fontsize{16}疑似肺質‘）
end


%%?4.去除非肺質
P?=?2000;
MASK?=?bwareaopen（objtmpP4）;??%?刪除面積小于P的連通分量
if?showflag
????figure
????subplot（131）imshow（objtmp）title（‘\fontsize{16}疑似肺質‘）
????subplot（132）imshow（MASK）title（‘\fontsize{16}肺質MASK‘）
????subplot（133）imshow（im）title（‘\fontsize{16}原圖‘）
end


%%?5.其他優化操作

%?略

figure
subplot（221）imshow（im）title（‘\fontsize{16}原圖‘）
subplot（222）imshow（imbi）title（‘\fontsize{16}胸腔‘）
subplot（223）imshow（objtmp）title（‘\fontsize{16}疑似肺質‘）
subplot（224）imshow（MASK）title（‘\fontsize{16}肺質MASK‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????320835??2019-07-18?11:11??imset.mat

?????文件???????2072??2019-07-18?11:18??demo.m

-----------?---------??----------?-----??----

???????????????322907????????????????????2