fenxingweishu.rar_fractal image_分形 图像_分形图像_分形维数 图像

I=imread(' E:\学习\大论文\测试图片\前车灯提取示例\09img_contour.jpg'); subplot(121); imshow(I); % I=rgb2gray(I); IM=edge(I); subplot(122); imshow(IM); [W,H]=size(I); Nn1=0; Nn2=0; Nn3=0; Nn4=0; Nn5=0; Nn6=0; Nn7=0; Nn8=0; Y=0; Tn=1; for i=1:Tn:W-Tn for j=1:Tn:H-Tn Y=find( IM(i:i+Tn-1,j:j+Tn-1)); if Y~=0 Nn1=Nn1+1; end end end Tn2=Tn*2; for i=1:Tn2:W-Tn2 for j=1:Tn2:H-Tn2 Y=find(IM(i:i+Tn2-1,j:j+Tn2-1)); if Y~=0 Nn2=Nn2+1; end end end Tn3=Tn2*2; for i=1:Tn3:W-Tn3 for j=1:Tn3:H-Tn3 Y=find( IM(i:i+Tn3-1,j:j+Tn3-1)); if Y~=0 Nn3=Nn3+1; end end end Tn4=Tn3*2; for i=1:Tn4:W-Tn4 for j=1:Tn4:H-Tn4 Y=find( IM(i:i+Tn4-1,j:j+Tn4-1)); if Y~=0 Nn4=Nn4+1; end end end Tn5=Tn4*2; for i=1:Tn5:W-Tn5 for j=1:Tn5:H-Tn5 Y=find( IM(i:i+Tn5-1,j:j+Tn5-1)); if Y~=0 Nn5=Nn5+1; end end end Tn6=Tn5*2; for i=1:Tn6:W-Tn6 for j=1:Tn6:H-Tn6 Y=find( IM(i:i+Tn6-1,j:j+Tn6-1)); if Y~=0 Nn6=Nn6+1; end end end Tn7=Tn6*2; for i=1:Tn7:W-Tn7 for j=1:Tn7:H-Tn7 Y=find( IM(i:i+Tn7-1,j:j+Tn7-1)); if Y~=0 Nn7=Nn7+1; end end end Tn8=Tn7*2; for i=1:Tn8:W-Tn8 for j=1:Tn8:H-Tn8 Y=find( IM(i:i+Tn8-1,j:j+Tn8-1)); if Y~=0 Nn8=Nn8+1; end end end Z=[Nn1,Nn2,Nn3,Nn4,Nn5,Nn6,Nn7,Nn8]; X=[-log(Tn),-log(Tn2),-log(Tn3),-log(Tn4),-log(Tn5),-log(Tn6),-log(Tn7),-log(Tn8)]; Y=[log(Nn1),log(Nn2),log(Nn3),log(Nn4),log(Nn5),log(Nn6),log(Nn7),log(Nn8)]; D=polyfit (X,Y,1); fxws = D(1)