图像分割MATLAB胸部CT肺质提取.rar

训练神经网络需要大量数据，目标检测中需要先进行标记。肺结节的检测非常普遍。现分享一个自己在胸部CT肺质标记的方法。

%  （图像分割）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