透视变换(单应性?)能提供更大的灵活性,但是一个透视投影并不是线性变换,因此所采用的映射矩阵是3*3,且控点变为4个,其他方面与仿射变换完全类似,下面的例程是针对密集变换,稀疏图像变换则采用cvPerspectiveTransform函数来处理。
------------------------------------------------------------------------------------------------
WarpPerspective
对图像进行透视变换
void cvWarpPerspective( const CvArr* src, CvArr* dst,const CvMat* map_matrix,
int flags=CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS,
CvScalar fillval=cvScalarAll(0) );
src
输入图像.
dst
输出图像.
map_matrix
3×3 变换矩阵
flags
插值方法和以下开关选项的组合:
· CV_WARP_FILL_OUTLIERS- 填充所有缩小图像的象素。如果部分象素落在输入图像的边界外,那么它们的值设定为 fillval.
· CV_WARP_INVERSE_MAP- 指定 matrix 是输出图像到输入图像的反变换,因此可以直接用来做象素插值。否则, 函数从 map_matrix 得到反变换。
fillval
用来填充边界外面的值
函数 cvWarpPerspective 利用下面指定矩阵变换输入图像:
要变换稀疏矩阵,使用 cxcore 中的函数 cvTransform 。
-----------------------------
GetPerspectiveTransform
由四对点计算透射变换
CvMat* cvGetPerspectiveTransform( const CvPoint2D32f*src, const CvPoint2D32f* dst,
CvMat*map_matrix );
#define cvWarpPerspectiveQMatrixcvGetPerspectiveTransform
src
输入图像的四边形顶点坐标。
dst
输出图像的相应的四边形顶点坐标。
map_matrix
指向3×3输出矩阵的指针。
函数cvGetPerspectiveTransform计算满足以下关系的透射变换矩阵:
这里,dst(i)= (x'i,y'i),src(i)= (xi,yi),i = 0..3.
------------------------------------------------------------------------------------------------
/*code*/
#include <highgui.h> #include <cv.h> int main(int argc, char** argv) { CvPoint2D32f srcTri[4], dstTri[4]; //二维坐标下的点,类型为浮点 //CvMat* rot_mat = cvCreateMat( 2, 3, CV_32FC1 ); //多通道矩阵 CvMat* warp_mat = cvCreateMat( 3, 3, CV_32FC1 ); IplImage *src, *dst; if( argc == 2 && ( ( src = cvLoadImage( argv[1], 1 ) ) != 0 ) ) { dst = cvCloneImage( src ); //制作图像的完整拷贝 dst ->origin = src ->origin; /* int origin; /* 0 - 顶—左结构, 1 - 底—左结构 (Windows bitmaps 风格) */ 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; //bot right 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( "Perspective Warp", 1 ); cvShowImage( "Perspective Warp", dst ); //最终是输出dst cvWaitKey(); } cvReleaseImage( &dst ); cvReleaseMat( &warp_mat ); return 0; }
source image
Perspective Transform