clear all;
clc;
M=6;
N=9; %solarcell的行数与列数
I=imread('.\testpic\milkdrop.gif');
%I=imread('.\testpic\peppers.gif');
%I=imread('.\testpic\left.png');
%I = imread('.\testpic\bacteria.BMP');
%I = imread('.\use\IMAGE15.bmp');
%I = imread('.\IRPIC\IR_all_s_g.BMP');
[m,n]=size(I);
imshow(I);
%中值滤波
R=zeros(m,n);
for i=2:m-1
for j=2:n-1
aa=[I(i-1,j-1) I(i-1,j) I(i-1,j+1) I(i,j-1) I(i,j) I(i,j+1) I(i+1,j-1) I(i+1,j) I(i+1,j+1)];
ii=0;
for jj=1:9-ii
for ii=1:8
if(aa(ii)>aa(ii+1))
cc=aa(ii);
aa(ii)=aa(ii+1);
aa(ii+1)=cc;
end
end
end
R(i,j)=aa(5);
end
end
S=uint8(R);
figure;
imshow(S);
%otsu 自动阈值计算
Y=zeros(1,256);
thresholdValue=1; %阈值初值
for i=2:m-1
for j=2:n-1
k=S(i,j);
Y(1,k+1)=Y(1,k+1)+1;
end
end
sum=0;
csum=0;
n0=0;
for k=1:256
sum=sum+k*double(Y(1,k));
n0=n0+Y(k);
end
fmax=-1;
n1=0;
for k=1:256
n1=n1+Y(k);
if(n1==0)
continue;
end
n2=n0-n1;
if(n2==0)
break;
end
csum=csum+k*double(Y(k));
m1=csum/n1;
m2=(sum-csum)/n2;
sb=n1*n2*(m1-m2)^2;
if(sb>fmax)
fmax=sb;
thresholdValue=k;
end
end
%图像二值化
U=zeros(m,n);
for i=1:m
for j=1:n
if(S(i,j)>=thresholdValue)
U(i,j)=0;
else
U(i,j)=1;
end
end
end
figure;
imshow(U);
%采用腐蚀的方法进行边缘检测
B=[1 1 1;1 1 1;1 1 1];
qsum=0;
QQ=zeros(m,n);
for i=2:m-1
for j=2:n-1
for ii=1:3
for jj=1:3
qsum=qsum+~xor(B(ii,jj),(1-U(i+ii-2,j+jj-2)));
end
end
if(qsum==9)
QQ(i,j)=0;
else
QQ(i,j)=1;
end
qsum=0;
end
end
SS=abs(QQ-U);
figure;
imshow(SS);
title('腐蚀处理后得到的边界图像');
SS1=SS;
for i=1:m
for j=1:n
if(i<=5 | j<=5 | i>=m-5 | j>=n-5)
SS1(i,j)=0;
end
end
end
figure;
imshow(SS1);
title('纯背景边界图像');
C=[0 1 0;1 1 1;0 1 0];
X=zeros(m,n);
flag_x=0;
for i=1:m
if(flag_x==1)
break;
end
for j=1:n
if(SS1(i,j)==1)
X(i,j)=1;
flag_x=1;
break;
end
end
end
figure;
imshow(X);
qsum=0;
k=1;
QQ=zeros(m,n);
X2=zeros(m,n);
ysum=1;
while(ysum~=0)
if(k>1)
X=X2;
end
k=k+1;
for i=2:m-1
for j=2:n-1
for ii=1:3
for jj=1:3
qsum=qsum+C(ii,jj)*X(i+ii-2,j+jj-2);
end
end
if(qsum>=1)
QQ(i,j)=1;
else
QQ(i,j)=0;
end
qsum=0;
end
end
X2=QQ&SS1;
Y2=(X2~=X);
ysum=0;
for iii=1:m
for jjj=1:n
ysum=ysum+Y2(iii,jjj);
end
end
end
figure;
imshow(X2);
figure;
imshow(X)
title('膨胀处理后的图像');
