运动目标检测_运动目标图像检测_matlab

for i=448:458 j=i+1; a1=imread([num2str(i),'.jpg']); b1=imread([num2str(j),'.jpg']); %a1=imread('448.jpg'); %b1=imread('449.jpg'); c1=a1-b1; c1=a1+1/(1+0.99)*(b1-a1); d1=b1-c1;%最小二乘法 %figure,imshow(d1) d2=im2bw(d1,8/255);%im2bw 将灰度图像转换为二进制图像 %figure,imshow(d2) bo=im2uint8(d2); bo=edge(d1,'sobel',0.004,'both'); k1=sum(bo'); %q=size(k1) n1=1; n2=1; for m1=1:240 if k1(m1)>1 n1=n1+1; if n1>35 g1(n2)=m1; end elseif n1>30 n2=n2+1; n1=1; else n1=1; end end %g1 n1=1; n2=1; for m1=240:-1:1 if k1(m1)>3 n1=n1+1; if n1>20 g2(n2)=m1; end elseif n1>20 n2=n2+1; n1=1; else n1=1; end end %g2 k2=sum(bo); n1=1; n2=1; for m1=1:320 if k2(m1)>4 n1=n1+1; if n1>20 g3(n2)=m1; end elseif n1>20 n2=n2+1; n1=1; else n1=1; end end %g3 n1=1; n2=1; for m1=320:-1:1 if k2(m1)>4 n1=n1+1; if n1>20 g4(n2)=m1; end elseif n1>20 n2=n2+1; n1=1; else n1=1; end end %g4 for m2=g2(2):g1(1) b1(m2,g4(2)-28)=0; b1(m2,g3(1)+8)=0; end %figure,imshow(b1) for m2=g2(1)+2:g1(2)+2 b1(m2,g4(1)-24)=0; b1(m2,g3(2)+13)=0; end %figure,imshow(b1) for m2=g4(2)-28:g3(1)+8 b1(g2(2),m2)=0; b1(g1(1),m2)=0; end for m2=g4(1)-24:g3(2)+13 b1(g2(1)+2,m2)=0; b1(g1(2)+2,m2)=0; end figure,imshow(b1) end