opencv 透视变换

来看透视变换是想用在 OpenCV 中的，但是由于太懒了，在 Matlab 中实现之后就没有再看，现在终于又看了看 OpenCV 的 API ，感觉比 Matlab 容易一些，在此写出与大家分享，样例图片跟在 Matlab 中的一样，这里就不再贴出来了，如果需要看的话可以看博主的前一篇文章在Matlab中实现透视变换的方法，另外本文为博主原创，转载请注明，谢谢~

样例点的设定还是参见博主前一篇文章，用的是一样的点对，这里给出 OpenCV 中的源代码和结果，这里用的是 C++ 版本的代码，有需要的朋友可以自行查找相关的 C 函数。

// warpPerspective.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include "cv.h"
#include <cstdio>
using namespace std;
using namespace cv;


int _tmain(int argc, _TCHAR* argv[])
{
	Point2f src[4]; //源图像中的点
	Point2f dst[4]; //目标图像中的点

	src[0].x = 1;
	src[0].y = 1;
	src[1].x = 1;
	src[1].y = 480;
	src[2].x = 640;
	src[2].y = 97;
	src[3].x = 640;
	src[3].y = 384;

	dst[0].x = 1;
	dst[0].y = 1;
	dst[1].x = 1;
	dst[1].y = 480;
	dst[2].x = 640;
	dst[2].y = 1;
	dst[3].x = 640;
	dst[3].y = 480;

	Mat t;  //变换矩阵
	t = getPerspectiveTransform(src,dst);

       // warpPerspective(src_gray, src_gray_dst, warp_mat, src_gray.size());

	printf("transform matrix\n");
	for(int i =0;i<3;i++)
	{
		printf("% .4f ",t.at<double>(0,i));
		printf("% .4f ",t.at<double>(1,i));
		printf("% .4f \n",t.at<double>(2,i));
	}

	Mat sample = (Mat_<double>(3,1)<<400,240,1); //样例点 (400,240)
	Mat r = t*sample;  //变换矩阵乘以样例点，得到结果点的其次坐标

	double s = r.at<double>(2,0);  //归一化尺度

	printf("\nresult...\n");
	printf("x = % .0f \n",ceil(r.at<double>(0,0)/s));  //别忘了除以 s 哦
	printf("y = % .0f ",ceil(r.at<double>(1,0)/s));

	return 0;
}

</pre><pre name="code" class="cpp"><pre name="code" class="cpp">#include<cv.h>
#include<highgui.h>

#pragma comment(lib, "cv.lib")
#pragma comment(lib, "cxcore.lib")
#pragma comment(lib, "highgui.lib")

int main()
{
	CvPoint2D32f srcTri[4], dstTri[4];
	CvMat*       warp_mat = cvCreateMat (3, 3, CV_32FC1);
	IplImage*    src = NULL;
	IplImage*    dst = NULL;

	src = cvLoadImage ("test.png", 1);
	dst = cvCloneImage (src);
	dst->origin = src->origin;
	cvZero (dst);

	srcTri[0].x = 0;
	srcTri[0].y = 0;
	srcTri[1].x = src->width - 1;
	srcTri[1].y = 0;
	srcTri[2].x = 0;
	srcTri[2].y = src->height - 1;
	srcTri[3].x = src->width - 1;
	srcTri[3].y = src->height - 1;

	dstTri[0].x = src->width * 0.05;
	dstTri[0].y = src->height * 0.33;
	dstTri[1].x = src->width * 0.9;
	dstTri[1].y = src->height * 0.25;
	dstTri[2].x = src->width * 0.2;
	dstTri[2].y = src->height * 0.7;
	dstTri[3].x = src->width * 0.8;
	dstTri[3].y = src->height * 0.9;

	cvGetPerspectiveTransform (srcTri, dstTri, warp_mat);
	cvWarpPerspective (src, dst, warp_mat);

	cvNamedWindow("src", 1);
	cvShowImage("src", src);
	cvNamedWindow ("Affine_Transform", 1);
	cvShowImage ("Affine_Transform", dst);

	cvWaitKey (0);

	cvReleaseImage (&src);
	cvReleaseImage (&dst);
	cvReleaseMat (&warp_mat);

	return 0;
}