matlab实现肺实质初始分割及采用滚球法

实现肺实质初始分割及采用滚球法（形态学开运算）进行修补肺实质轮廓凹陷。

%%滚球法修补肺实质
% 两个参数：滚球半径、去除气管参数，可自行选择
%椭圆半径15，高度10，计算时间1.102
%  没有气管
%
%


clear;
clc;
tic
A=imread（‘IMG-0018-00001.jpg‘）;%读取原图像
B=rgb2gray（A）;%将原图像转换为灰度图像
t=graythresh（B）;%计算阈值t
C=im2bw（Bt）;%根据阈值二值化图像
D=imfill（C8‘holes‘）;%对二值化后的图像填充肺实质
E=D-C;%得到肺实质的图像E
F=imfill（E8‘holes‘）;%填充肺实质空洞
FMask=bwareaopen（F1000）;%去除面积小于T的部分（气管），取值可变
figure（）imshow（FMask[]）;%得到掩膜
r_ball=30;%可变的，取值为
se_ball=strel（‘ball‘r_ball10）;%椭圆体半径40，高度10
r_disk=ceil（r_ball/6）;%圆整r_ball/6得到大于或等于它的最接近整数，为5
    if r_disk==0;
        r_disk=1;%最小为1
    end
se_erode=strel（‘disk‘r_disk）; %圆形半径
%mask=imopen（FMaskr_disk）;
%figure（）imshow（mask[]）;
L=bwlabel（FMask）;   %二值图像标记，将图像分成多个区域
%stat = regionprops（FMask）;%计算图像区域特征，区域连通，object为二值图像，
  [rowcol]=size（B）;
mask_leftlung=im2bw（zeros（rowco

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-12-02 22:24  gunqiufa\
     文件      145267  2015-11-28 15:40  gunqiufa\000084.jpg
     目录           0  2015-12-07 10:53  gunqiufa\1\
     文件      145267  2015-11-28 15:40  gunqiufa\1\000084.jpg
     文件       37839  2015-12-07 10:53  gunqiufa\1\gunqiufa-15-1.102-000084-比较.jpg
     文件       14203  2015-11-28 17:14  gunqiufa\1\gunqiufa-15-1.102-000084.jpg
     文件        2739  2015-11-28 17:17  gunqiufa\1\gunqiufa.m
     目录           0  2015-12-07 10:34  gunqiufa\2\
     文件      150458  2015-11-28 15:40  gunqiufa\2\IMG-0001-00001.jpg
     文件        7793  2015-11-28 17:04  gunqiufa\2\gunqiufa-30-2.243-.IMG-0001-00001.jpg
     文件        7912  2015-11-28 16:56  gunqiufa\2\gunqiufa-40-2.844-IMG-0001-00001.jpg
     文件       18558  2015-12-07 11:01  gunqiufa\2\gunqiufa-50-4.091IMG-0001-00001-比较.jpg
     文件        7709  2015-11-28 17:06  gunqiufa\2\gunqiufa-50-4.091IMG-0001-00001.jpg
     文件        2753  2015-11-28 17:18  gunqiufa\2\gunqiufa.m
     目录           0  2015-11-28 17:26  gunqiufa\3\
     文件      143019  2015-11-28 15:40  gunqiufa\3\IMG-0001-00002.jpg
     文件        8517  2015-11-28 17:20  gunqiufa\3\gunqiufa-15-1.093-IMG-0001-00002.jpg
     文件        8620  2015-11-28 17:21  gunqiufa\3\gunqiufa-20-1.471-IMG-0001-00002.jpg
     文件        8566  2015-11-28 17:24  gunqiufa\3\gunqiufa-30-1.940-IMG-0001-00002.jpg
     文件       19409  2015-12-07 11:19  gunqiufa\3\gunqiufa-40-3.061-IMG-0001-00002.jpg
     文件        2864  2015-11-28 17:25  gunqiufa\3\gunqiufa.m
     目录           0  2015-11-28 18:53  gunqiufa\4\
     文件      142265  2015-11-28 15:40  gunqiufa\4\IMG-0018-00001.jpg
     文件       11036  2015-11-28 17:28  gunqiufa\4\gunqiufa-20-1.516-IMG-0018-00001.jpg
     文件       10792  2015-11-28 17:30  gunqiufa\4\gunqiufa-25-1.760-IMG-0018-00001.jpg
     文件       10993  2015-11-28 17:31  gunqiufa\4\gunqiufa-30-2.004-IMG-0018-00001.jpg
     文件       10747  2015-11-28 17:34  gunqiufa\4\gunqiufa-35-2.637-IMG-0018-00001.jpg
     文件        2863  2015-11-28 18:53  gunqiufa\4\gunqiufa.m
     文件      150458  2015-11-28 15:40  gunqiufa\IMG-0001-00001.jpg
     文件      143019  2015-11-28 15:40  gunqiufa\IMG-0001-00002.jpg
     文件      142265  2015-11-28 15:40  gunqiufa\IMG-0018-00001.jpg
............此处省略8个文件信息