bubianju.txt

这是一个不变距算法的程序

function m=bubianju2(nargin)
% M is the seven bubianju


iw=imread(nargin);
ir=IM2BW(iw,0.3);

% imshow(ir);
[i,j]=find(ir==1);
longer=nnz(i);
m01=0;m10=0;m00=0;m20=0;m11=0;
for ai=min(i):max(i);
for aj=min(j):max(j);
m01=m01+ai^0*aj^1*ir(ai,aj);
m10=m10+ai^1*aj^0*ir(ai,aj);
m00=m00+ai^0*aj^0*ir(ai,aj);
m20=m20+ai^2*aj^0*ir(ai,aj);
m11=m11+ai^1*aj^1*ir(ai,aj);
end
end 
n02=0;n20=0;n12=0;n21=0;n03=0;n30=0;n11=0;n00=0;n01=0;n10=0;
for bi=min(i):max(i);
for bj=min(j):max(j);
n00=n00+((bi-m10/m00)^0)*((bj-m01/m00)^0)*ir(bi,bj);
n01=n01+((bi-m10/m00)^0)*((bj-m01/m00)^1)*ir(bi,bj);
n10=n10+((bi-m10/m00)^1)*((bj-m01/m00)^0)*ir(bi,bj);
n02=n02+((bi-m10/m00)^0)*((bj-m01/m00)^2)*ir(bi,bj);
n20=n20+((bi-m10/m00)^2)*((bj-m01/m00)^0)*ir(bi,bj);
n12=n12+((bi-m10/m00)^1)*((bj-m01/m00)^2)*ir(bi,bj);
n21=n21+((bi-m10/m00)^2)*((bj-m01/m00)^1)*ir(bi,bj);
n03=n03+((bi-m10/m00)^0)*((bj-m01/m00)^3)*ir(bi,bj);
n30=n30+((bi-m10/m00)^3)*((bj-m01/m00)^0)*ir(bi,bj);
n11=n11+((bi-m10/m00)^1)*((bj-m01/m00)^1)*ir(bi,bj);
end
end 
p02=n02/n00;
p20=n20/n00;
p12=n12/n00^1.25;
p21=n21/n00^1.25;
p03=n03/n00^1.25;
p30=n30/n00^1.25;
p11=n11/n00;
q1=p20+p02
q2=(p20-p02)^2+4*p11
q3=(p30-3*p21)^2+(3*p21-p03)^2
q4=(p30+p12)^2+(p21+p03)
q5=(p30-3*p21)*(p21+p30)*((p30+p21)^2-3*(p21+p03)^2)+(3*p21-p03)*(p21+p03)*(3*(p30+p12)^2-(p21+p03)^2)
q6=(p20-p02)*((p30+p12)^2-(p21+p03)^2)+4*p11*(p30+p12)*(p21+p03)
q7=(3*p21-p03)*(p30+p12)*((p30+p12)^2-3*(p21+p03)^2)+(p30-3*p12)*(p21+p03)*(3*(p30+p12)^2-(p21+p03)^2)
L=[q1,q2,q3,q4,q5,q6,q7];