D8算法Matlab实现

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