www.pudn.com > rice.rar > rice.m

I=imread('rice.jpg'); 
I_en=wiener2(I); 
I_en=medfilt2(I_en,[5,5]); 
f=[1,2,1;2,4,2;1,2,1]/16; 
I_en=conv2(I_en,f); 
I_en=uint8(I_en); 
subplot(2,2,1),imshow(I_en); 
title({'维纳-中值-高斯','增强的图像'}); 
 
se = strel('disk', 20); 
Itop = imtophat (I_en, se); 
Ibot = imbothat (I_en, se); 
I_en = imsubtract(imadd (Itop, I_en), Ibot); 
I_en = imcomplement(I_en); 
subplot(2,2,2), imshow (I_en), title ('增强对比度后反相'); 
 
Iemin = imextendedmin (I_en, 80); 
Iimpose = imimposemin (I_en, Iemin); 
subplot(2,2,3), imshow(Iimpose), title('将谷底置零'); 
 
wat = watershed (Iimpose); 
stats = regionprops(wat, 'Area', 'Orientation') 
area = [ stats(:).Area ]; 
x1=length(find(area(:)<2000)); 
subplot(2,2,4), imshow(wat), title({'分水岭算法结果',['米粒个数为',num2str(x1)]});