1、将原图扩增几圈
w=3; %结构元素取3x3方阵
Imgg=zeros(M+w-1,N+w-1); %将原图扩增几圈
for i=1+(w-1)/2:M+(w-1)/2
for j=1+(w-1)/2:N+(w-1)/2
Imgg(i,j)=Img(i-(w-1)/2,j-(w-1)/2);
end
end
for i=1+(w-1)/2:M+(w-1)/2
for j=1+(w-1)/2:N+(w-1)/2
temp=zeros(1,w*w);
k=1;
for x=-(w-1)/2:(w-1)/2
for y=-(w-1)/2:(w-1)/2
ii=(i-x);
jj=(j-y);
temp(k)=Imgg(ii,jj);
k=k+1;
end
end
e=max(temp)-min(temp);
g=(1+e)^(-2);