OpenCV仿射变换——旋转

一、原理与数学推导

1.1图示

1、旋转前和旋转后位置关系图，θ是旋转的角度

2、首先需要平移一次，将图像平移到所要得到的旋转后矩阵的中心，这里需要通过计算旋转后的cv::Mat对象大小

3、通过中心点为旋转中心，旋转图像

1.2公式推导：

假设旋转角度为θ，旋转中心为图片对象的中心。

cv::Mat src = cv::imread(path);//path为图片路径
int col=src.cols;
int row=src.rows;

通过公式证明，可以得到新的矩阵大小：

1.2.1 平移
平移可以参照之前的文章OpenCV仿射变换——平移.
1.2.3 旋转
Opencv中的旋转是以Y轴向下为正，逆时针转动：

OA转到OB，角度为θ
OB转到OC，角度为β
知道B的坐标B(x,y),求C的坐标C(X’,y’)
由图可知;

由于：

所以可以得到：

同理：

写成矩阵格式为：

函数调用形式：
cv::Mat getRotationMatrix2D(cv::Point2f center, double angle, double scale)

参数：
Point2f center：表示旋转的中心点
double angle：表示旋转的角度
double scale：图像缩放因子

二、实现代码

#include<iostream>
#include<opencv2/core.hpp>
#include<opencv2/imgproc.hpp>
#include<opencv2/highgui.hpp>

int main() {
	cv::Mat src = cv::imread("D:\\Windows\\WorkSpace\\imageTest\\OIP.jpg",cv::IMREAD_COLOR);
	if (src.empty()) {
		std::cout << "failed to read image." << std::endl;
		return EXIT_FAILURE;
	}

	int col = src.cols;
	int row = src.rows;
	//旋转角度
	double angle = 30;

	double angle_radian = angle * CV_PI / 180;

	//计算新生成的图像矩阵的大小
	int _new_col_ = col * std::sin(angle_radian) + row * std::cos(angle_radian);
	int _new_row_ = col * std::cos(angle_radian) + row * std::sin(angle_radian);
	
	cv::Point2f pt1[3], pt2[3];
	float xoffset = _new_col_ / 2 - col / 2;
	float yoffset = _new_row_/2 - row / 2;

	/*平移前的位置*/
	pt1[0] = cv::Point2f(0, 0);//            0  0
	pt1[1] = cv::Point2f(xoffset, 0);//      x  0
	pt1[2] = cv::Point2f(0, yoffset);//      0  y
	/*平移后的位置*/
	pt2[0] = cv::Point2f(xoffset, yoffset);//         x  y
	pt2[1] = cv::Point2f(xoffset * 2, yoffset);//    2x  y
	pt2[2] = cv::Point2f(xoffset, yoffset * 2);//     x  2y

	//生成只有平移的防射矩阵
	cv::Mat T;
	T = cv::getAffineTransform(pt1, pt2);
	cv::Mat res = cv::Mat::zeros(_new_col_, _new_row_, CV_32FC1);
	//平移
	cv::warpAffine(src, res, T, res.size());


	//定义旋转中心
	cv::Point2f center(_new_col_ / 2, _new_row_ / 2);

	//确定旋转后的缩放大小
	double scale = 1;//缩放比例不变

	//定义用于旋转的仿射矩阵
	cv::Mat M;
	M = cv::getRotationMatrix2D(center, angle, scale);
	std::cout << "M= " << M << std::endl;
	
	//生成新的显示结果的矩阵（图）
	cv::Mat dst = cv::Mat::zeros(_new_col_, _new_row_, CV_32FC1);
	
	//旋转
	cv::warpAffine(res, dst, M, dst.size());
	
	cv::namedWindow("旋转前", cv::WINDOW_AUTOSIZE);
	cv::imshow("旋转前", src);
	cv::namedWindow("平移后", cv::WINDOW_AUTOSIZE);
	cv::imshow("平移后", res);
	cv::namedWindow("旋转后", cv::WINDOW_AUTOSIZE);
	cv::imshow("旋转后", dst);
	cv::waitKey(0);
}

三、运行结果

3.1图片效果：

3.2仿射矩阵M