w1=5; k=floor(w1/2); g=double(origpicgray); for r=3:m-2 for c=3:n-2 V1=0;V2=0;V3=0;V4=0; for i=-k:k-1 V1=V1+(g(r,c+i)-g(r,c+i+1))^2; %东西方向 V2=V2+(g(r+i,c+i)-g(r+i+1,c+i+1))^2; %东南方向 V3=V3+(g(r+i,c)-g(r+i+1,c))^2; %南北方向 V4=V4+(g(r-i,c+i)-g(r-i-1,c+i+1))^2; %东北方向 end IV(r,c)=min([V1 V2 V3 V4]);