%*********************************************
%***************交通視頻車輛跟蹤***************
%******************幀間差分法******************
%*********************************************
clcclearclose?all;
%讀入視頻，獲取視頻文件信息
videofile=VideoReader（‘C:\Users\陳璇\Desktop\v4.avi‘）;
get（videofile）;
nframes=get（videofile‘NumberOfframes‘）;
I=read（videofile1）;%讀取視頻中的第1幀
taggedCars=zeros（[size（I1）?size（I2）?3?nframes]class（I））;

????singleframe=read（videofile13）;
????I1=rgb2gray（singleframe）;%將圖像轉化為灰度圖
????figure;imshow（I1）title（‘第13幀灰度圖‘）;print（‘-dpng‘‘第13幀灰度圖‘）;
????I2=read（videofile15）;
????I2=rgb2gray（I2）;%將圖像轉化為灰度圖
????figure;imshow（I2）title（‘第15幀灰度圖‘）;print（‘-dpng‘‘第15幀灰度圖‘）;
????
????
????
????car1=medfilt2（I1）;%中值濾波去噪
????figure;imshow（car1）title（‘第13幀中值濾波去噪‘）;print（‘-dpng‘‘第13幀中值濾波去噪‘）;
????car2=medfilt2（I2）;%中值濾波去噪
????figure;imshow（car2）title（‘第15幀中值濾波去噪‘）;print（‘-dpng‘‘第15幀中值濾波去噪‘）;
????car1=im2double（car1）;
????car2=im2double（car2）;
????
????cars=car2-car1;%差分?
????figure;imshow（cars）title（‘13、15幀差分圖‘）;print（‘-dpng‘‘13、15幀差分圖‘）;
????cars=histeq（cars）;%直方圖均衡化
????figure;imshow（cars）title（‘直方圖均衡化‘）;print（‘-d