clear all;
clc;
>> I = imread('trees.tif');
>> figure,imshow(I);title('原始图像');
>> J = im2bw(I,graythresh(I));
>> figure,imshow(J),title('二值图像');
>> J = ~J;
>> figure,imshow(J);title('取反二值图像');
>> %去除边界白色非目标
>> J = imclearborder(J,4);%imclearborder:support light structures connected to image border,抑制图像和边界相连的亮对象
>> figure,imshow(J);title('目标图像');
>> %腐蚀膨胀操作,去除小的白色点
>> se = strel('disk',1);
>> J = imerode(J,se); %图像腐蚀
>> J = imdilate(J,se); %图像膨胀
>> figure,imshow(J);title('最终图像');
注:graythresh():最大类间方差找到合适的阈值
imclearborder():抑制图像和边界相连的亮对象
strel('disk',1):创建形态学结构元素,使用一个1×1的圆形结构元素对象
imerode(J,se); %图像腐蚀
imdilate(J,se); %图像膨胀