OpenCV中图像旋转(warpAffine)算法的实现过程
在OpenCV中，目前并没有现成的函数直接用来实现图像旋转，它是用仿射变换函数cv::warpAffine来实现的，此函数目前支持4种插值算法，最近邻、双线性、双三次、兰索斯插值，如果传进去的参数为基于像素区域关系插值算法(INTER_AREA)，则按双线性插值。
通常使用2*3矩阵来表示仿射变换：
其中，T相当于变换前的原始图像，x,y为变换后的图像坐标。
对于cv::getRotationMatrix2D函数的实现公式为：
其中scale为缩放因子(x、y方向保持一致)，angle为旋转角度(弧长)，centerx，centery为旋转中心。
以lena.jpg图像旋转45度为例：
采用最近邻插值算法的实现代码为：
cv::Mat matSrc = cv::imread("lena.jpg", 2 | 4);if (matSrc.empty()) return;const double degree = 45;double angle = degree * CV_PI / 180.;double alpha = cos(angle);double beta = sin(angle);int iWidth = matSrc.cols;int iHeight = matSrc.rows;int iNewWidth = cvRound(iWidth * fabs(alpha) + iHeight * fabs(beta));int iNewHeight = cvRound(iHeight * fabs(alpha) + iWidth * fabs(beta));double m[6];m[0] = alpha;m[1] = beta;m[2] = (1 - alpha) * iWidth / 2. - beta * iHeight / 2.;m[3] = -m[1];m[4] = m[0];m[5] = beta * iWidth / 2. + (1 - alpha) * iHeight / 2.;cv::Mat M = cv::Mat(2, 3, CV_64F, m);cv::Mat matDst1 = cv::Mat(cv::Size(iNewWidth, iNewHeight), matSrc.type(), cv::Scalar::all(0));double D = m[0]*m[4] - m[1]*m[3];D = D != 0 ? 1./D : 0;double A11 = m[4]*D, A22 = m[0]*D;m[0] = A11; m[1] *= -D;m[3] *= -D; m[4] = A22;double b1 = -m[0]*m[2] - m[1]*m[5];double b2 = -m[3]*m[2] - m[4]*m[5];m[2] = b1; m[5] = b2;int round_delta = 512;//nearestfor (int y=0; y<iNewHeight; ++y){for (int x=0; x<iNewWidth; ++x){//int tmpx = cvFloor(m[0] * x + m[1] * y + m[2]);//int tmpy = cvFloor(m[3] * x + m[4] * y + m[5]);int adelta = cv::saturate_cast<int>(m[0] * x * 1024);int bdelta = cv::saturate_cast<int>(m[3] * x * 1024);int X0 = cv::saturate_cast<int>((m[1] * y + m[2]) * 1024) + round_delta;int Y0 = cv::saturate_cast<int>((m[4] * y + m[5]) * 1024) + round_delta;int X = (X0 + adelta) >> 10;int Y = (Y0 + bdelta) >> 10;if ((unsigned)X < iWidth && (unsigned)Y < iHeight){matDst1.at<cv::Vec3b>(y, x) = matSrc.at<cv::Vec3b>(Y, X);}}}cv::imwrite("rotate_nearest_1.jpg", matDst1);M = cv::getRotationMatrix2D(cv::Point2f(iWidth / 2., iHeight / 2.), degree, 1);cv::Mat matDst2;cv::warpAffine(matSrc, matDst2, M, cv::Size(iNewWidth, iNewHeight), 0, 0, 0);cv::imwrite("rotate_nearest_2.jpg", matDst2);采用双线性插值算法的实现代码为：
cv::Mat matSrc = cv::imread("lena.jpg", 2 | 4);if (matSrc.empty()) return;const double degree = 45;double angle = degree * CV_PI / 180.;double alpha = cos(angle);double beta = sin(angle);int iWidth = matSrc.cols;int iHeight = matSrc.rows;int iNewWidth = cvRound(iWidth * fabs(alpha) + iHeight * fabs(beta));int iNewHeight = cvRound(iHeight * fabs(alpha) + iWidth * fabs(beta));double m[6];m[0] = alpha;m[1] = beta;m[2] = (1 - alpha) * iWidth / 2. - beta * iHeight / 2.;m[3] = -m[1];m[4] = m[0];m[5] = beta * iWidth / 2. + (1 - alpha) * iHeight / 2.;cv::Mat M = cv::Mat(2, 3, CV_64F, m);cv::Mat matDst1 = cv::Mat(cv::Size(iNewWidth, iNewHeight), matSrc.type(), cv::Scalar::all(0));double D = m[0]*m[4] - m[1]*m[3];D = D != 0 ? 1./D : 0;double A11 = m[4]*D, A22 = m[0]*D;m[0] = A11; m[1] *= -D;m[3] *= -D; m[4] = A22;double b1 = -m[0]*m[2] - m[1]*m[5];double b2 = -m[3]*m[2] - m[4]*m[5];m[2] = b1; m[5] = b2;for (int y=0; y<iNewHeight; ++y){for (int x=0; x<iNewWidth; ++x){//int tmpx = cvFloor(m[0] * x + m[1] * y + m[2]);//int tmpy = cvFloor(m[3] * x + m[4] * y + m[5]);float fx = m[0] * x + m[1] * y + m[2];float fy = m[3] * x + m[4] * y + m[5];int sy = cvFloor(fy);fy -= sy;//sy = std::min(sy, iHeight-2);//sy = std::max(0, sy);if (sy < 0 || sy >= iHeight) continue;short cbufy[2];cbufy[0] = cv::saturate_cast<short>((1.f - fy) * 2048);cbufy[1] = 2048 - cbufy[0];int sx = cvFloor(fx);fx -= sx;//if (sx < 0) {// fx = 0, sx = 0;//}//if (sx >= iWidth - 1) {// fx = 0, sx = iWidth - 2;//}if (sx < 0 || sx >= iWidth) continue;short cbufx[2];cbufx[0] = cv::saturate_cast<short>((1.f - fx) * 2048);cbufx[1] = 2048 - cbufx[0];for (int k=0; k<matSrc.channels(); ++k){if (sy == iHeight - 1 || sx == iWidth - 1) {continue;} else {matDst1.at<cv::Vec3b>(y, x)[k] = (matSrc.at<cv::Vec3b>(sy, sx)[k] * cbufx[0] * cbufy[0] +matSrc.at<cv::Vec3b>(sy+1, sx)[k] * cbufx[0] * cbufy[1] +matSrc.at<cv::Vec3b>(sy, sx+1)[k] * cbufx[1] * cbufy[0] +matSrc.at<cv::Vec3b>(sy+1, sx+1)[k] * cbufx[1] * cbufy[1]) >> 22;}}}}cv::imwrite("rotate_bilinear_1.jpg", matDst1);M = cv::getRotationMatrix2D(cv::Point2f(iWidth / 2., iHeight / 2.), degree, 1);cv::Mat matDst2;cv::warpAffine(matSrc, matDst2, M, cv::Size(iNewWidth, iNewHeight), 1, 0, 0);cv::imwrite("rotate_bilinear_2.jpg", matDst2);
其它插值算法的实现代码与双线性类似，可参考 http://blog.csdn.net/fengbingchun/article/details/17335477
相关文章：
10亿美元续命！OpenAI获微软投资，意在通用人工智能？
来源 | OpenAI官博译者 | 孙薇编辑 | 一一出品 | AI科技大本营（ID:rgznai100） 7 月 22 日，微软宣布将对非营利人工智能研究组织 OpenAI 投资 10 亿美元，用于通用人工智能（AGI）的开发。双方将以微软原有的公…
TrayIcon 类 添加系统托盘不显示托盘图标
为什么80%的码农都做不了架构师？>>> 好久不碰 java swing最近写了一个swing 程序 添加托盘时，怎么也不显示图标，就一空白 ，在网上搜了老半天，大部无效。 边看帖子边看 java api ,结合理解,有一属性设置了一…
Crystal Report 加载模板报错 无法在c++ 堆栈中打开由jrc 引擎处理的文档
2019独角兽企业重金招聘Python工程师标准>>> 纠结了很久, 尝试过录入一个错误的路径,文件读取也是包相同的错误,也就是表示找不到路径文件而已,并不是开发环境的问题 于是设置一个最简单的路径,放置下去rpt模板,代码后续没有报错; 再次之前也修复了一个关于引用的dl…
Ubuntu下makefile及gcc生成静态库动态库的简单使用举例
环境：Ubuntu-13.10 32位(虚拟机)、gcc4.8.1 首先创建一个test_makefile_gcc文件夹，此test_makefile_gcc文件夹下包括：src文件夹用于存放源文件； include文件夹用于存放头文件；bin文件夹用于存放生成的动态库.so文件&…
Exchange Server 2013 安装完成后配置外部URL
Exchange Server 2013 安装完成后配置外部URL 比如 mail.contoso.com 1、转到 EAC → “服务器”，然后单击“配置外部访问域”。2、在“选择要与外部 URL 一起使用的客户端访问服务器”下面，单击“添加”3、选择您要配置的客户端访问服务器，…
Ubuntu下CodeBlocks的安装、配置及静态库动态库的简单使用举例
1、 从Ubuntu Software Center中搜索Code::Blocks并安装； 2、 在第一次启动时选择GNU GCC Compiler作为默认的编译器； 3、 生成静态库并调用操作步骤，代码同 http://blog.csdn.net/fengbingchun/article/details/17994489 3.1、New fil…
React 打怪笔记
介绍 本文为学习react中的记录。 Tips: 当组件的props或state有变化，执行render函数。无论是使用函数或是类来声明一个组件，它决不能修改它自己的propsReact 可以将多个setState() 调用合并成一个调用来提高性能。无状态函数式组件 (stateless functiona…
新闻智能分类练习赛开始报名啦！最先达到80分就可以领GPU，技术书籍！
现代信息爆炸般地产生，信息如海如潮。信息分类，不仅有利于加快信息检索速度，且有利于提高查准率。Internet是信息的重要载体，深入地研究与探讨网上信息自动分类的方法、技术和理论，已成为时代的迫切需求和新的研究热点…
项目经理应该具备的技能
作为一个优秀的项目经理应该具备五个方面的技能: 项目管理知识体系 应用领域的相关知识、标准和规则项目环境知识一般管理知识软技能/人际关系技能项目管理知识体系 就是要掌握常说的9大知识领域：范围、时间、成本、质量、人力资源、风险、沟通、采购再加上集成…
SSE2 Intrinsics各函数介绍
SIMD相关头文件包括： //#include <ivec.h>//MMX//#include <fvec.h>//SSE(also include ivec.h)//#include <dvec.h>//SSE2(also include fvec.h)#include <mmintrin.h> //MMX#include <xmmintrin.h> //SSE(include mmintrin.h)#incl…
中国城市道路名图鉴
作者| AlfredWu来源 | Alfred数据室（ID：Alfred_Lab）不知道大家出差或旅游的时候有没有发现，有些城市特别喜欢使用其它城市或者省份的名字作为道路名，特别是青岛市，这会儿还在徐州路，走过一个街区…
中国电子信息产业发展研究院主办的2018中国软件大会上大快搜索“又双叒叕”获奖了...
大快搜索自荣获“2018中国大数据企业50强”殊荣，12月20日在由工信部指导，中国电子信息产业化发展研究院主办的2018中国软件大会上，大快搜索获评“2018中国大数据基础软件领域领军企业”称号，入选中国数字化转型TOP100服务商&#…
SystemCenter2012SP1实践(15)共享库服务器和ISO
用过HyperV的同学都知道，HyperV调用ISO作为启动光盘的时候，必须保存在本地才行。网络共享下的一概不认。在SCVMM下，我们可以通过一些设置，让SCVMM下创建的虚拟机，支持调用不在同一台主机上的ISO文件。我是分隔线首先要…
XLNet：公平PK，BERT你已经被超过！
作者 | XLNet Team译者 | 孙薇责编 | Jane出品 | AI科技大本营（ID: rgznai100）【导语】几周前，XLNet 团队发布了新型预训练语言模型 XLNet，这个新模型在各项基准测试中都优于谷歌之前发布的BERT模型，其中模型 XLNet-La…
C/C++中switch用法的一种替换方式
在C/C中，switch语句是经常被用到的，当switch内的case语句较多时程序有时显得比较繁乱，此种情况下可以用另外一种实现方式替代switch。详细用法见例子： #include "stdafx.h"float AddFunc(float a, float b) {return (a …
My excellent 2018
又到一年年底时，除了感叹一句时间过得真快之外，也非常庆幸自己能够渡过了这精彩的一年。 工作 首先是工作上面的变化。自16年毕业以后就在招联消费金融有限公司上班，一直到今年的五月份。在招联的期间，有幸得到老大的赏识…
POJ 2955 Brackets （区间DP）
题目链接：http://poj.org/problem?id2955 BracketsTime Limit: 1000MS Memory Limit: 65536KTotal Submissions: 1977 Accepted: 1012Description We give the following inductive definition of a “regular brackets” sequence: the empty sequence is a regul…
从芯片到AI智能芯片，一文了解它的前世今生
作者 | 元宵大师，Python高级工程师，致力于推动人工智能、大数据分析在金融量化交易领域中的应用。欢迎大家关注我的个人公众号《元宵大师带你用Python量化交易》。责编 | 胡巍巍来源 | CSDN（ID：CSDNnews）经过长期的发展…
Windows7下OpenGL简单使用举例
1、 从http://www.opengl.org/resources/libraries/glut/glut_downloads.php下载glut相关头文件和库glutdlls37beta.zip ，(默认的windows机子上并没有glut头文件及相应的库，它主要用来打开窗口、开发和管理菜单，以及管理事件等)，…
Snagit9-12注册码
SnagIt 9 注册码： AM5SC-8LWML-MVMWU-DTLGE-ERMBE SnagIt 10 注册码： 5HCAK-DEGMZ-EYABA-M4LCC-ACBE2 DFKDA-JZ5FC-TGLAA-CM5DM-MFEBD CMCFH-93DCD-SFZYC-K5KCM-C7CA7 SnagIt 11 注册码： 7CTCC-5WQCS-98AY8-V8F2M-76258 NCTCC-5WFCK-98A28-V8…
Strut2访问
访问HelloWorld应用的路径的设置 在struts2中，访问struts2中action的URL路径由两部份组成： 包的命名空间action的名称 例如: 访问本例子HelloWorldAction的URL路径为： /l6n/helloWorldAction (注意：完整路径为：http:/…
单v100 GPU，4小时搜索到一个鲁棒的网络结构
作者 | Slumbers，毕业于中山大学，深度学习工程师，主要方向是目标检测，语义分割，GAN责编 | JaneNAS最近也很火，正好看到了这篇论文，解读一下，这篇论文是基于DAG（directed …
关于pyecharts 地图显示添加数据的问题
echarts ： 香港地区显示（人口密集的人口数目） http://echarts.baidu.com/examples/editor.html?cmap-HK series: [ { name: 香港18区人口密度, type: map, mapType: HK …
MMX Intrinsics各函数介绍
SIMD相关头文件包括： //#include <ivec.h>//MMX //#include <fvec.h>//SSE(also include ivec.h) //#include <dvec.h>//SSE2(also include fvec.h)#include <mmintrin.h> //MMX #include <xmmintrin.h> //SSE(include mmintrin.h) #…
大数据中台向AI中台演进是大势所趋？
来源 | 宜信技术学院（ID:CE_TECH）导读：自从阿里巴巴提出“中台”的概念之后，这个词汇就成为各领域企业关注的焦点，很多人在考虑建设自己的中台。然而，构建中台是否真有必要？是否所有的企业都要建…
WordPress标签
1、分类目录调用函数： <?php wp_list_cats();?> 2、调用页面函数： <?php wp_nav_menu( array( theme_location > ast-menu-primary, container > false ) ); ?> 3、转载于:https://blog.51cto.com/okowo/…
SSE3和SSSE3 Intrinsics各函数介绍
SIMD相关头文件包括： //#include <ivec.h>//MMX //#include <fvec.h>//SSE(also include ivec.h) //#include <dvec.h>//SSE2(also include fvec.h)#include <mmintrin.h> //MMX #include <xmmintrin.h> //SSE(include mmintrin.h) #…
kubernetes学习笔记之十三：基于calico的网络策略入门
一、.安装calico [rootk8s-master01 ~]# kubectl apply -f https://docs.projectcalico.org/v3.3/getting-started/kubernetes/installation/hosted/canal/rbac.yaml clusterrole.rbac.authorization.k8s.io "calico" created clusterrole.rbac.authorization.k8s.i…
设计模式之抽象工厂模式(Abstract Factory)摘录
面向对象系统的分析和设计实际上追求的就是两点：高内聚(Cohesion)和低耦合(Coupling). 23种GOF设计模式一般分为三大类：创建型模式、结构型模式、行为模式。 创建型模式包括：1、FactoryMethod(工厂方法模式)；2、Abstract Factor…
grub设置密码的方法
grub设置密码的方法：一、grub设置明文口令的方法：修改/etc/grub.conf配置文件就可以了.[root RedHat ~] # vi /etc/grub.conf #注：此为链接文件，指向 /boot/grub/grub.conf#boot/dev/hdbdefault0timeout5splashimage(hd0,0)/gr…