基于視頻的交通事故汽車車速計(jì)算研究汽車服務(wù)畢業(yè)論-資料下載頁(yè)
2025-06-24 17:39本頁(yè)面
【正文】 u=(e(i)+f(i)*p(1)p(2)*p(1))/(p(1)^2+1) v=(p(1)*e(i)+f(i)*p(1)^2p(2)*p(1)^2)/(p(1)^2+1)+p(2) c=[c u] d=[d v] plot(c(i),d(i),39。g.39。)。hold onendglobal c。global d。global len。 global im。 ；global bdd ；global x1；global x2；global x3；global y1；global y2；global y3global DE；global EF；global AB；global BCAB=str2num(get(,39。string39。))/2。L1=length(AB)。if L11 errordlg(39。未輸入車長(zhǎng)或軸距39。)。 return。endBC=AB 。bdd=[]。axes()A=rgb2gray(im)。imshow(A)。button = 1。n = 0。a = []。b = []。hold on。while button == 1 % 輸入數(shù)據(jù)點(diǎn)，以右鍵作結(jié) [xi,yi,button] = ginput(1)。 n = n + 1。 if n 1 plot([xi,a(n1,1)], [yi,b(n1,1)], 39。rs39。,39。LineWidth39。,2,... 39。MarkerEdgeColor39。,39。k39。,... 39。MarkerFaceColor39。,39。g39。,... 39。MarkerSize39。,10)。 else plot(xi, yi, 39。s39。,39。LineWidth39。,2,... 39。MarkerEdgeColor39。,39。k39。,... 39。MarkerFaceColor39。,39。g39。,... 39。MarkerSize39。,10)。 end a(n,1) = xi。 b(n,1) = yi。endplot([xi,a(1,1)], [yi,b(1,1)], 39。rs39。,39。LineWidth39。,2,... 39。MarkerEdgeColor39。,39。k39。,... 39。MarkerFaceColor39。,39。g39。,... 39。MarkerSize39。,10)。hold off。if length(a) = 3 va = round(a)。 vb = round(b)。 x = 1:size(A,1)。 y = 1:size(A,2)。 [Y,X] = meshgrid(y,x)。 in = inpolygon(Y(:),X(:),va,vb)。 AA = reshape(A, [size(A,1)*size(A,2),size(A,3)])。 BB = repmat(uint8(255),size(AA))。 BB(in,:) = AA(in,:)。 B = reshape(BB, [size(A,1),size(A,2),size(A,3)])。 bw=im2bw(B,)。endbw1 = ~bw。%黑白互換bw1=bwareaopen(bw1,500)。[r, h]=find(bw1==1)。[rectx, recty, area, perimeter] = minboundrect(h, r, 39。a39。)。%作最小矩形l=((rectx(1)rectx(2))^2+(recty(2)recty(1))^2)^。%矩形長(zhǎng)度da=str2num(get(,39。string39。))。%修正系數(shù)l=l+da。%修正后的長(zhǎng)度imshow(bw1)。 hold online(rectx, recty, 39。color39。, 39。r39。)。for i=1:len plot(c(i),d(i),39。.39。)。end [p s]= polyfit(c,d,1)。x=0::1000。y=p(1)*x+p(2)。plot(x,y)axis oh = imline(gca,[10 100], [100 100])。 %做用于標(biāo)定的線setColor(h,[0 1 0])。id = addNewPositionCallback(h,@(pos) title(mat2str(pos,3)))。api = iptgetapi(h)。pausepos = () u1=pos(1,1)。v1=pos(1,2)。%[x1,y1] = ginput(1)。 x1=(u1+v1*p(1)p(2)*p(1))/(p(1)^2+1) y1=(p(1)*u1+v1*p(1)^2p(2)*p(1)^2)/(p(1)^2+1)+p(2)plot(x1,y1,39。r.39。)%D點(diǎn)bdd=[bdd [x1。y1]]。u3=pos(2,1)。v3=pos(2,2)%[x3,y3]=ginput(1)。x3=(u3+v3*p(1)p(2)*p(1))/(p(1)^2+1)y3=(p(1)*u3+v3*p(1)^2p(2)*p(1)^2)/(p(1)^2+1)+p(2)plot(x3,y3,39。r.39。)%F點(diǎn)bdd=[bdd [x3。y3]]。z1=str2num(get(,39。string39。))。z=z1/(1z1)。x2=(x1+z*x3)/(1+z)。y2=(y1+z*y3)/(1+z)。plot(x2,y2,39。y.39。)plot(x2,y2,39。r.39。)bdd=[bdd [x2。y2]]。DE=[(x2x1)^2+(y2y1)^2]^。EF=[(x2x3)^2+(y2y3)^2]^。d1=[(x3x1)^2+(y3y1)^2]^。wc=100*(ld1)/lset(,39。string39。,num2str(wc))。global y3。global DE。global EF。global AB。global BC。global len。global c。global d 。global x1。global x2。global x3。global y1。global y2。global y3。 sk=[]。sd=[]for j=2:len1 t=j1。s=j+1。 %t=1。s=3 MD=[(x1c(t))^2+(y1d(t))^2]^if c(t)x1NA1=2*AB*EF*MD/(EF*DE+DE*DE+MD*DEMD*EF)。%N在最下方，實(shí)現(xiàn)endif c(t)x1 NA1=2*AB*EF*MD/(EF*DE+DE*DEMD*DE+MD*EF)。 %M在近間斷 endMD=[(x1c(s))^2+(y1d(s))^2]^。if c(s)x1 NA2=2*AB*EF*MD/(EF*DE+DE*DE+MD*DEMD*EF)。%N在最下方，實(shí)現(xiàn)endif c(s)x3 NA2=2*AB*EF*MD/(EF*MD+DE*EF+DE*DEMD*DE)。%M在F最遠(yuǎn)處 %NA=AB*EF*MD*(AB+BC)/(EF*MD*ABBC*DE*MD+BC*DE*DE+BC*EF*DE)。endif c(s)x1 c NA2=2*AB*EF*MD/(EF*DE+DE*DEMD*DE+MD*EF)。 %M在近間斷End Dmn=abs(NA2NA1)tt=str2num(get(,39。string39。))。sj=*tt*abs(ts)p=(t+s)/2。D=Dmn/1000。v=*D/(sj) sk=[sk p]；sd=[sd v]；%速度和時(shí)刻endxmin=min(sk)。xmax=max(sk)+2。ymin=min(sd)10。ymax=max(sd)+10。figure,axis([xmin,xmax,ymin,ymax])for jj=1:len2 plot(sk(jj),sd(jj),39。.39。),hold onendmeanv=mean(sd)set(,39。string39。,num2str(meanv))。grid on 3
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1