基于奇异值分解（SVD）的背景抑制 --Opencv实现

最新推荐文章于 2022-12-13 11:27:50 发布

沧浪之水清兮

最新推荐文章于 2022-12-13 11:27:50 发布

阅读量1.4k

点赞数 1

分类专栏： Opencv

本文链接：https://blog.csdn.net/sinat_36412790/article/details/80104836

版权

Opencv 专栏收录该内容

36 篇文章 4 订阅

订阅专栏

int main()
{
	const char* infile = "F:\\文\\%d.jpg";
	char inPath[100] = { 0 };

	for (int i = 8400; i <= 12000; i++)
	{
		sprintf(inPath, infile, i);
		Mat image = imread(inPath, 0);                 // 读入原图
		double t1 = getTickCount();
		Mat temp(image.size(), CV_32FC1, Scalar(0));
		image.convertTo(image, CV_32FC1);
		Mat U, W, V;
		SVD::compute(image, W, U, V);//opencv得到的V与MATLAB相比已经经过转置了，要想再转置一遍可以用V=V.t();
		image.convertTo(image, CV_8UC1);
		imshow("0", image);
		Mat w(image.rows, image.rows, CV_32FC1, Scalar(0));//opencv进行SVD分解后得到的奇异值不是放入对角矩阵，而是一个列向量中，所以需要自己将其变换为对角矩阵
		for (int i = 0; i < 10; i++)
			w.ptr<float>(i)[i] = W.ptr<float>(i)[0];
		temp = U*w*V;
		temp.convertTo(temp, CV_8UC1);
		temp = temp - image;
		double minv = 0.0, maxv = 0.0;
		double* minp = &minv;
		double* maxp = &maxv;

		minMaxIdx(temp, minp, maxp);
		double T = minv + 0.33*(maxv - minv);
		temp = temp > T;
		double t2 = getTickCount();
		std::cout << i << ": " << (t2 - t1) / getTickFrequency() * 1000 << endl;
		imshow("1", temp);
		waitKey(1);
	}
	return 0;
}