颗粒分析

%颗粒分析
f=imread('plugs.tif');
subplot(2,2,1),imshow(f),title('原图像')

sumpixels=zeros(1,36);
for k=0:35
    se=strel('disk',k);
    fo=imopen(f,se);
    sumpixels(k+1)=sum(fo(:));
end

subplot(2,2,2),plot(0:35,sumpixels),xlabel('k'),ylabel('表面积')

subplot(2,2,3),plot(-diff(sumpixels))
xlabel('k')
ylabel('表面区域减少')

fasf=f
sumpixels=zeros(1,36);
for k=0:35
    se=strel('disk',k);
    fo=imclose(imopen(fasf,se),se);
    sumpixels(k+1)=sum(fo(:));
end

subplot(2,2,4),plot(-diff(sumpixels))
xlabel('k')
ylabel('表面区域减少')