边缘检测后计算区域周长

周长length

clear all;
clc;
I=imread('F:\\bi\\original2.jpg');
P=rgb2gray(I);
itedge = edge(P,'canny');
[row,col]=size(P);
M=zeros(row,col);
length=0;
for ii=1:row
    for jj=1:col
     if  ( itedge(ii,jj) == 1 )
       if ( itedge(ii,jj+1))
          length=length+1;
            if ( itedge(ii+1,jj) == 1)
                length=length+1;
                if ( itedge(ii+1,jj-1) == 1)
                    length=length+sqrt(2);
                end
            end 
       end
     end
    end
end