火焰圆形度计算

火焰圆形度，定义是

clear all;
clc;
I=imread('F:\\bi\\original2.jpg');
P=rgb2gray(I);
K = imresize (P ,0.2);
[row,col]=size(K);
S=0;
//累加1值点，即火焰总面积
for ii=1:row
    for jj=1:col
    if  ( K(ii,jj) == 1 ) 
        S=S+1;
    end
    end
end
itedge = edge(K,'canny');
//求周长
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

C=4*pi*S/length;