求帧差

本文转自http://blog.csdn.net/abcjennifer/article/details/6621986

1. //**********************************************************************          
2. //@2011/7/21 10:05:21   
3. //@Author: 钱建 & 张睿卿   
4. //**********************************************************************  
5.   
6. #include"iostream"   
7. #include "highgui.h"   
8. #include <conio.h>   
9. #include <stdlib.h>   
10. #include "math.h"   
11. #include "cv.h"   
12. #include "highgui.h"   
13. #include "cxcore.h"   
14. using namespace std;  
15.   
16. int Width;  
17. int Height;  
18. int T;  
19. IplImage* IplImageFrameSub;  
20. unsigned char *templast,*tempcur,*temp2y;  
21. CvScalar s,s1;  
22.   
23. void FD()  
24. {  
25.     int i,j;  
26.   
27.     IplImage *file_input_src1=cvLoadImage("C:\\Users\\Zrq\\Desktop\\c1.bmp",1);  
28.     IplImage *file_input_src2=cvLoadImage("C:\\Users\\Zrq\\Desktop\\c2.bmp",1);  
29.   
30.     Width=min(file_input_src1->width,file_input_src2->width);  
31.     Height=min(file_input_src1->height,file_input_src2->height);  
32.   
33.     templast=(unsigned char*)malloc(Height *Width*3/2*sizeof(unsigned char*));  
34.     tempcur=(unsigned char*)malloc(Height *Width*3/2*sizeof(unsigned char*));  
35.     temp2y=(unsigned char*)malloc(Height *Width*3/2*sizeof(unsigned char*));  
36. IplImageFrameSub=cvCloneImage(file_input_src1);  
37.   
38.     for (i=0; i<Height; i++)  
39.         for (j=0; j<Width; j++)  
40.         {  
41.             s=cvGet2D(file_input_src1,i,j);  
42.             templast[i*Width+j]=unsigned char(s.val[0]);  
43.             s=cvGet2D(file_input_src2,i,j);  
44.             tempcur[i*Width+j]=unsigned char(s.val[0]);  
45.         }  
46.         cout<<"Load Succeed!"<<endl;  
47.   
48.         for (j=0; j<Height; j++)  
49.             for (i=0; i<Width; i++)  
50.             {  
51.                 temp2y[j*Width+i]=((abs(templast[j*Width+i]-tempcur[j*Width+i]))<=T ? 0: 255);  
52.             }  
53.   
54.             for (j = 0;j < Height;j++)  
55.                 for(i = 0;i < Width;i++)  
56.                     templast[j*Width+i] = tempcur[j * Width + i];  
57.   
58.             cout<<"Begin Save BMP"<<endl;  
59.             for (j=0; j<Height; j++)  
60.                 for (i=0; i<Width; i++)  
61.                 {  
62.                     s1.val[0]=temp2y[j*Width+i];  
63.                     cvSet2D(IplImageFrameSub,j,i,s1);  
64.                 }  
65.                 cout<<"OK!"<<endl;  
66.                 cvSaveImage("D:\\数字图象处理\\MVD-break61(帧差1).bmp",IplImageFrameSub);  
67. }  
68.   
69. void main()  
70. {  
71.     cout<<"input T"<<endl;  
72.     cin>>T;  
73.     FD();  
74. }