火焰圆形度计算

火焰圆形度，定义是

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    endenditedge = 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    endendC=4*pi*S/length;