#include?“cv.h“
#include?“highgui.h“
#include?
#include?
#include?
#include?
#include?

//?部分跟蹤參數（秒）
const?double?MHI_DURATION?=?0.5;
const?double?MAX_TIME_DELTA?=?0.5;
const?double?MIN_TIME_DELTA?=?0.05;
const?int?N?=?3;

//
const?int?CONTOUR_MAX_AERA?=?400;

//?ring?image?buffer
IplImage?**buf?=?0;
int?last?=?0;

//?當前視頻幀
IplImage?*mhi?=?0;?//?MHI:?motion?history?image

CvFilter?filter?=?CV_GAUSSIAN_5x5;
CvConnectedComp?*cur_comp?min_comp;
CvConnectedComp?comp;
CvMemStorage?*storage;
CvPoint?pt[4];

//??參數：
//??img?–?輸入視頻幀
//??dst?–?檢測結果
void??update_mhi（?IplImage*?img?IplImage*?dst?int?diff_threshold?）
{
????double?timestamp?=?clock（）/100.;?//?get?current?time?in?seconds
????CvSize?size?=?cvSize（img->widthimg->height）;?//?get?current?frame?size
????int?i?j?idx1?idx2;
????IplImage*?silh;
????uchar?val;
????float?temp;
????IplImage*?pyr?=?cvCreateImage（?cvSize（（size.width?&?-2）/2?（size.height?&?-2）/2）?8?1?）;
????CvMemStorage?*stor;
????CvSeq?*cont?*result?*squares;
????CvSeqReader?reader;

????if（?!mhi?||?mhi->width?!=?size.width?||?mhi->height?!=?size.height?）?
????{
????????if（?buf?==?0?）?
????????{
????????????buf?=?（IplImage**）malloc（N*sizeof（buf[0]））;
????????????memset（?buf?0?N*sizeof（buf[0]））;
????????}
????????
????????for（?i?=?0;?i?????????{
????????????cvReleaseImage（?&buf[i]?）;
????????????buf[i]?=?cvCreateImage（?size?IPL_DEPTH_8U?1?）;
????????????cvZero（?buf[i]?）;
????????}
????????cvReleaseImage（?&mhi?）;
????????mhi?=?cvCreateImage（?size?IPL_DEPTH_32F?1?）;
????????cvZero（?mhi?）;?//?clear?MHI?at?the?beginning
????}?//?end?of?if（mhi）

????cvCvtColor（?img?buf[last]?CV_BGR2GRAY?）;?//?convert?frame?to?grayscale

????idx1?=?last;
????idx2?=?（last?+?1）?%?N;?//?index?of?（last?-?（N-1））th?frame?
????last?=?idx2;

????//?做幀差
????silh?=?buf[idx2];
????cvAbsDiff（?buf[idx1]?buf[idx2]?silh?）;?//?get?difference?between?frames
????
????//?對差圖像做二值化
????cvThreshold（?silh?silh?30?255?CV_THRESH_BINARY?）;?//?and?threshold?it
????
????cvUpdateMotionHistory（?silh?mhi?timestamp?MHI_DURATION?）;?//?update?MHI
????cvCvtScale（?mhi?dst?255./MHI_DURATION?
??????（MHI_DURATION?-?timestamp）*255./MHI_DURATION?）;????
????cvCvtScale（?mhi?dst?255./MHI_DURATION?0?）;????
????
????//?中值濾波，消除小的噪聲
????cvSmooth（?dst?dst?CV_MEDIAN?3?0?0?0?）;
????
????//?向下采樣，去掉噪聲
????cvPyrDown（?dst?pyr?CV_GAUSSIAN_5x5?）;
????cvDilate（?pyr?pyr?0?1?）;??//?做膨脹操作，消除目標的不連續空洞
????cvPyrUp（?pyr?dst?CV_GAUSSIAN_5x5?）;
????//
????//?下面的程序段用來找到輪廓
????//
????//?建立點陣序列保存輪廓點陣.
????stor?=?cvCreateMemStorage（0）;
????cont?=?cvCreateSeq（CV_SEQ_ELTYPE_POINT?sizeof（CvSeq）?sizeof（CvPoint）??stor）;
????
????//?找到所有輪廓
????cvFindContours（?dst?stor?&cont?sizeof（CvContour）?
????????????????????CV_RETR_LIST?CV_CHAIN_APPROX_SIMPLE?cvPoint（00））;

????//?直接使用CONTOUR中的矩形來畫輪廓
????for（;cont;cont?=?cont->h_next）
????{
??????????????CvRect?r?=?（（CvContour*）cont）->rect;
??????????????if（r.height?*?r.width?>?CONTOUR_M

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????188502??2012-03-09?14:39??多運動物體跟蹤\Debug\motiondetect.exe

?????文件??????20352??2012-03-09?14:39??多運動物體跟蹤\Debug\motiondetect.obj

?????文件?????476160??2012-03-09?14:39??多運動物體跟蹤\Debug\motiondetect.pdb

?????文件??????77824??2012-03-09?14:39??多運動物體跟蹤\Debug\vc60.pdb

?????文件????????140??2007-02-24?14:01??多運動物體跟蹤\Makefile

?????文件???????4686??2012-03-09?14:39??多運動物體跟蹤\motiondetect.c

?????文件???????4220??2007-03-07?20:41??多運動物體跟蹤\motiondetect.dsp

?????文件????????549??2007-03-07?20:41??多運動物體跟蹤\motiondetect.dsw

?????文件??????50176??2012-03-09?14:39??多運動物體跟蹤\motiondetect.ncb

?????文件??????48640??2012-03-09?14:39??多運動物體跟蹤\motiondetect.opt

?????文件???????1968??2012-03-09?14:39??多運動物體跟蹤\motiondetect.plg

?????文件????1746620??2008-06-04?11:57??多運動物體跟蹤\test.avi

?????文件???????4687??2012-03-08?15:43??多運動物體跟蹤\代碼.txt

?????目錄??????????0??2012-04-19?23:19??多運動物體跟蹤\Debug

?????目錄??????????0??2012-03-09?14:39??多運動物體跟蹤

-----------?---------??----------?-----??----

??????????????2624524????????????????????15