差分法求解运动物体轨迹

差分法求解运动物体轨迹

数据集中的物体横向运动，下面使用差分法检测物体的运动轨迹，计算运动物体的重心坐标。

实验环境

VS2019+OpenCV
自制数据集中一共80张32*32的彩色图像，二值化阈值设定为15（可以根据需求自行调整）

#include<iostream>
#include<fstream>
#include "opencv2/opencv.hpp"
#include <opencv2\highgui\highgui.hpp>
#include<string>
using namespace std;
using namespace cv;
int main() {
	char fileName[500];
	IplImage* img1;
	IplImage* ImgGray1;
	IplImage* img2;
	IplImage* ImgGray2;
	ofstream ou("result.csv");
	ou <<"id"<<","<< "area" << "," << "x" << "," << "y" << endl;
	img1 = cvLoadImage(".//picture//0.jpg", 1);//读取彩色图像
	ImgGray1 = cvCreateImage(cvGetSize(img1), IPL_DEPTH_8U, 1);//创建灰度图像
	cvCvtColor(img1, ImgGray1, CV_BGR2GRAY);//彩色图像转换为灰度图像
	for (int i = 1; i < 80; i++) {
		sprintf_s(fileName, ".//picture//%d.jpg", i);//提取出文价夹下面的每张图片路径，放到fileName
			img2 = cvLoadImage(fileName, 1);//读取彩色图像
			ImgGray2 = cvCreateImage(cvGetSize(img2), IPL_DEPTH_8U, 1);//创建灰度图像
			cvCvtColor(img2, ImgGray2, CV_BGR2GRAY);//彩色图像转换为灰度图像
			IplImage* img = cvCreateImage(cvGetSize(ImgGray2), ImgGray2->depth, ImgGray2->nChannels);
			cvSub(ImgGray2,ImgGray1, img);//ImgGray2减去ImgGray1给img
			cvThreshold(img, img, 15, 255, CV_THRESH_BINARY);//转化为二值图像，分割阈值为15
			//计算重心
			CvMoments* moments = (CvMoments*)malloc(sizeof(CvMoments));
			cvMoments(img, moments, 1);
			double area = cvGetSpatialMoment(moments, 0, 0);
			if (area) {
				double moment10 = cvGetSpatialMoment(moments, 1, 0);
				double moment01 = cvGetSpatialMoment(moments, 0, 1);
				double posX = moment10 / area;
				double posY = moment01 / area;
				//显示重心
				CvPoint center; center.x = int(posX); center.y = int(posY);
				cout << "area=" << area << "\t" << "x= " << center.x << "\t" << "y= " << center.y << endl;
				ou <<i<<","<<area<<","<< center.x << "," << center.y << endl;
			}
			else
			{
				cout <<"area=" << 0 <<"\t"<< "x= " << 0 << "\t" << "y= " << 0 << endl;
				ou << i << "," << 0 << "," << 0 << "," << 0 << endl;
			}
			img1 = cvCloneImage(img2);//复制图像
			ImgGray1 = cvCloneImage(ImgGray2);//复制图像
	}
	ou.close();
	return(0);
}

运动轨迹图

上图显示运动物体重心坐标变化情况，可以看出物体在横向移动。