function[ im2]=smallest(imt)
[m,n]=size(imt);
for i=1:m-3
for j=1:n-3
%一号模板
G1=[imt(i,j),imt(i+1,j),imt(i+2,j);imt(i,j+1),imt(i+1,j+1),imt(i+2,j+1);0,imt(i+1,j+2),0];
mean1 =(imt(i,j)+imt(i+1,j)+imt(i+2,j)+imt(i,j+1)+imt(i+1,j+1)+imt(i+2,j+1)+imt(i+1,j+2))/7;
variance1=(1/mean1^2)*((imt(i,j)-mean1)^2+(imt(i+1,j)-mean1)^2+(imt(i+2,j)-mean1)^2+(imt(i,j+1)-mean1)^2+(imt(i+1,j+1)-mean1)^2+(imt(i+2,j+1)-mean1)^2+(imt(i+1,j+2)-mean1)^2);
%二号模板
G2=[0,imt(i+1,j),0;imt(i,j+1),imt(i+1,j+1),imt(i+2,j+1);imt(i,j+2),imt(i+1,j+2),imt(i+2,j+2)];
mean2 =(imt(i+1,j)+imt(i+1,j)+imt(i+2,j+2)+imt(i,j+1)+imt(i+1,j+1)+imt(i+2,j+1)+imt(i+1,j+2))/7;
variance2=(1/mean2^2)*((imt(i+1,j)-mean2)^2+(imt(i+1,j)-mean2)^2+(imt(i+2,j+2)-mean2)^2+(imt(i,j+1)-mean2)^2+(imt(i+1,j+1)-mean2)^2+(imt(i+2,j+1)-mean2)^2+(imt(i+1,j+2)-mean2)^2);
%三号模板
G3=[imt(i,j),imt(i+1,j),0;imt(i,j+1),imt(i+1,j+1),imt(i+2,j+1);imt(i,j+2),imt(i+1,j+2),0];
mean3 =(imt(i,j)+imt(i+1,j)+imt(i,j+2)+imt(i,j+1)+imt(i+1,j+1)+imt(i+2,j+1)+imt(i+1,j+2))/7;
variance3=(1/mean3^2)*((imt(i,j)-mean3)^2+(imt(i+1,j)-mean3)^2+(imt(i,j+2)-mean3)^2+(imt(i,j+1)-mean3)^2+(imt(i+1,j+1)-mean3)^2+(imt(i+2,j+1)-mean3)^2+(imt(i+1,j+2)-mean3)^2);
%四号模板
G4=[0,imt(i+1,j),imt(i+2,j);imt(i,j+1),imt(i+1,j+1),imt(i+2,j+1);0,imt(i+1,j+2),imt(i+2,j+2)];
mean4 =(imt(i+2,j)+imt(i+1,j)+imt(i+2,j+2)+imt(i,j+1)+imt(i+1,j+1)+imt(i+2,j+1)+imt(i+1,j+2))/7;
variance4=(1/mean4^2)*((imt(i+2,j)-mean4)^2+(imt(i+1,j)-mean4)^2+(imt(i+2,j+2)-mean4)^2+(imt(i,j+1)-mean4)^2+(imt(i+1,j+1)-mean4)^2+(imt(i+2,j+1)-mean4)^2+(imt(i+1,j+2)-mean4)^2);
%对应模板取最小方差
A=[variance1,variance2,variance3,variance4];
minvariance=min(A);
%取最小方差模板的均值替换目标位置
if minvariance==variance1
im2(i,j)=mean1;
end
if minvariance==variance2
im2(i,j)=mean2;
end
if minvariance==variance3
im2(i,j)=mean3;
end
if minvariance==variance4
im2(i,j)=mean4;
end
end
end
end