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)
