数字图像处理-旋转

#include<opencv2/opencv.hpp>
#include<opencv2/core.hpp>
cv::Mat myRotationbYCenter(const cv::Mat& src, double angle)
{
	angle = angle * CV_PI / 180.0; //角度转弧度
	int ow = src.cols;//图片宽
	int oh = src.rows;//图片高
	int dst_w = round(fabs(ow * cos(angle)) + fabs(oh * sin(angle)));
	int dst_h = round(fabs(ow * sin(angle)) + fabs(oh * cos(angle)));
	cv::Mat dst(dst_w, dst_h, CV_8UC3);
	dst.setTo(255);
	cv::Mat T1 = (cv::Mat_<double>(3, 3) << 1.0, 0.0, 0.0, 0.0, -1.0, 0.0, -0.5 * src.cols, 0.5 * src.rows, 1.0); 
	cv::Mat T2 = (cv::Mat_<double>(3, 3) << cos(angle), -sin(angle), 0.0, sin(angle), cos(angle), 0.0, 0.0, 0.0, 1.0); 
	cv::Mat T3 = (cv::Mat_<double>(3, 3) << 1.0,0.0,0.0,0.0,-1.0,0.0,0.5*src.cols,0.5*src.rows,1.0);
	cv::Mat T = T1 * T2 * T3;  //目标像素相对于原像素的转移矩阵
	cv::Mat TT = T.inv();
	if (src.type() == CV_8UC3) //三通道
	{
		for (int x = 0; x < dst.cols; x++)
		{
			for (int y = 0; y < dst.rows; y++)
			{
				cv::Mat temp = (cv::Mat_<double>(1, 3) << x, y, 1);
				temp = temp*TT;  //计算原始坐标旋转后的坐标
				int newW = round(temp.at<double>(0));
				int newH = round(temp.at<double>(1));
				if (newW >= src.cols || newH >= src.rows || newW<0|| newH<0)
				{
					continue;
				}
				dst.at<cv::Vec3b>(y, x)[0] = src.at<cv::Vec3b>(newH, newW)[0];
				dst.at<cv::Vec3b>(y, x)[1] = src.at<cv::Vec3b>(newH, newW)[1];
				dst.at<cv::Vec3b>(y, x)[2] = src.at<cv::Vec3b>(newH, newW)[2];
			}
		}
	}
	return dst;
}
int main()
{
	cv::Mat img = cv::imread("C:/Users/Administrator/Desktop/test.jpg");
	cv::Mat dst = myRotationbYCenter(img,-20);
	cv::imshow("dst",dst);
	cv::waitKey(0);
}