# opencv实现图像任意角度旋转的算法解析及代码实现

976人阅读 评论(0)

#include "CXCORE.h"
#include "cv.h"
#include "highgui.h"
using namespace cv;
IplImage*m_img;
double max(double a,double b)
{
return (a > b) ? a : b;
}
void ImgRotate(bool direction)
{
int oldWidth = m_img->width;
int oldHeight = m_img->height;

// 源图四个角的坐标（以图像中心为坐标系原点）
float fSrcX1, fSrcY1, fSrcX2, fSrcY2, fSrcX3, fSrcY3, fSrcX4, fSrcY4;
fSrcX1 = (float)(-(oldWidth) / 2);
fSrcY1 = (float)((oldHeight) / 2);//第二象限
fSrcX2 = (float)((oldWidth ) / 2);
fSrcY2 = (float)((oldHeight) / 2);//第一象限
fSrcX3 = (float)(-(oldWidth) / 2);
fSrcY3 = (float)(-(oldHeight) / 2);//第四象限
fSrcX4 = (float)((oldWidth) / 2);
fSrcY4 = (float)(-(oldHeight) / 2);//第三象限

// 旋转后四个角的坐标（以图像中心为坐标系原点）
float fDstX1, fDstY1, fDstX2, fDstY2, fDstX3, fDstY3, fDstX4, fDstY4;
float theta = 0.5*CV_PI*direction;
fDstX1 = cos(theta) * fSrcX1 + sin(theta) * fSrcY1;
fDstY1 = -sin(theta) * fSrcX1 + cos(theta) * fSrcY1;
fDstX2 = cos(theta) * fSrcX2 + sin(theta) * fSrcY2;
fDstY2 = -sin(theta) * fSrcX2 + cos(theta) * fSrcY2;
fDstX3 = cos(theta) * fSrcX3 + sin(theta) * fSrcY3;
fDstY3 = -sin(theta) * fSrcX3 + cos(theta) * fSrcY3;
fDstX4 = cos(theta) * fSrcX4 + sin(theta) * fSrcY4;
fDstY4 = -sin(theta) * fSrcX4 + cos(theta) * fSrcY4;
//新的宽度和高度
int newWidth =  (max(fabs(fDstX4 - fDstX1), fabs(fDstX3 - fDstX2)) );
int newHeight = (max(fabs(fDstY4 - fDstY1), fabs(fDstY3 - fDstY2)) );

IplImage*dst = cvCreateImage(cvSize(newWidth, newHeight), m_img->depth, m_img->nChannels);
//见旋转原理，行列式中第三行第一二列
float dx = -0.5*newWidth*cos(theta) - 0.5*newHeight*sin(theta) + 0.5*oldWidth;
float dy = 0.5*newWidth*sin(theta) - 0.5*newHeight*cos(theta) + 0.5*oldHeight;
//读取新图片每一个点
for (int height = 0; height<newHeight; height++)
{
uchar*ptrNow = (uchar*)(dst->imageData + height*dst->widthStep);//指向图片的第height行第0列
for (int width = 0; width<newWidth; width++)
{
int x = float(width)*cos(theta) + float(height)*sin(theta) + dx;
int y = float(-width)*sin(theta) + float(height)*cos(theta) + dy;//对应原始图像的点
uchar*ptrOld = (uchar*)(m_img->imageData + y*m_img->widthStep);//指向原始图像的第y行第0列
if ((x<0) || (x >= oldWidth) || (y<0) || (y >= oldHeight))//判断是否是没有图像的区域
{
if (m_img->nChannels == 3)//判断图像是几通道图像
{
ptrNow[3 * width + 1] = 0;
ptrNow[3 * width + 2] = 0;
ptrNow[3 * width + 3] = 0;
}
else if (m_img->nChannels == 1)
{
ptrNow[width] = 0;
}
}
else
{
if (m_img->nChannels == 3)
{
ptrNow[3 * width + 1] = ptrOld[3 * x + 1];
ptrNow[3 * width + 2] = ptrOld[3 * x + 2];
ptrNow[3 * width + 3] = ptrOld[3 * x + 3];
}
else if (m_img->nChannels == 1)
{
ptrNow[width] = ptrOld[x];
}
}
}
}
cvZero(m_img);
m_img = cvCreateImage(cvSize(newWidth, newHeight), dst->depth, dst->nChannels);
cvCopy(dst, m_img);
cvReleaseImage(&dst);
}

int main(int argc, char*argv[])
{
ImgRotate(true);
cvNamedWindow("show_image", 0);
cvShowImage("show_image", m_img);
cvWaitKey(0);
return 0;
}

0
0

* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
个人资料
• 访问：2786次
• 积分：58
• 等级：
• 排名：千里之外
• 原创：3篇
• 转载：0篇
• 译文：0篇
• 评论：0条
文章分类
文章存档