function?r=funx（path）
tic
P=imread（path）;
[rowcol]=size（P）;%row?為行數，col為列數
result=zeros（rowcol）;?%初始化result為0
r=result;
DB=zeros（18）;%第一列為E、二列為SE、三列S、四列SW、五列W、六列NW、七列N、八列NE方向
for?i=1:row
????for?j=1:col??
?????if?i~=1
?????????DB（7）=P（ij）-P（i-1j）;?%北方時-------------N
?????else
?????????DB（7）=-1;
?????end
?????if?j~=col
?????????DB（1）=P（ij）-P（ij+1）;%東方時-------------E
?????else
?????????DB（1）=-1;
?????end
??????if?i~=row?&j~=col
?????????DB（2）=（P（ij）-P（i+1j+1））/sqrt（2）;%東南方時-----SE
?????else
?????????DB（2）=-1;
??????end
???????if?i~=row
?????????DB（3）=P（ij）-P（i+1j）;?%南方時，與南方那點的距離-------S
??????else
?????????DB（3）=-1;
???????end
????????if?i~=row?&?j~=1
?????????DB（4）=（P（ij）-P（i+1j-1））/sqrt（2）;%西南方時----SW
?????else
?????????DB（4）=-1;
?????end
????