clear;
close all;
I=imread('D:\瓷砖.bmp');
[M,N]=size(I);
for m=2:M-1
for n=2:N-1
D1(m,n)=abs(I(m-1,n)-I(m+1,n)); % 0度方向上的差分矩阵
D2(m,n)=abs(I(m,n-1)-I(m,n+1)); % 90度方向上的差分矩阵
D3(m,n)=abs(I(m+1,n-1)-I(m-1,n+1)); % 45度方向上的差分矩阵
D4(m,n)=abs(I(m-1,n-1)-I(m+1,n+1)); % 135度方向上的差分矩阵
end
end