仿射变换详解 warpAffine

今天遇到一个问题是关于仿射变换的，但是由于没有将仿射变换的具体原理型明白，看别人的代码看的很费解，最后终于在师兄的帮助下将原理弄明白了，我觉得最重要的是理解仿射变换可以看成是几种简单变换的复合实现，

具体实现形式即将几种简单变换的变换矩阵M相乘，这样就很容易理解啦

定义：仿射变换的功能是从二维坐标到二维坐标之间的线性变换，且保持二维图形的“平直性”和“平行性”。仿射变换可以通过一系列的原子变换的复合来实现，包括平移，缩放，翻转，旋转和剪切。

这类变换可以用一个3*3的矩阵M来表示，其最后一行为（0，0，1）。该变换矩阵将原坐标为（x,y)变换为新坐标（x',y')，

即

$\begin{bmatrix} x'\\ y'\\ 1 \end{bmatrix}= \begin{bmatrix} m_{00} &m_{01} & m_{02}\\ m_{10}& m_{11} & m_{12} \\ 0 & 0 &1 \end{bmatrix} \begin{bmatrix} x\\ y\\ 1 \end{bmatrix}$

opencv中相应的函数是：

void warpAffine(InputArray src, OutputArray dst, InputArray M, Size dsize, int flags=INTER_LINEAR, int borderMode=BORDER_CONSTANT, const Scalar& borderValue=Scalar())¶

Parameters:

src – input image.
dst – output image that has the size dsize and the same type as src .
M – $2\times 3$ transformation matrix，最重要的东东了，本文中着重讲M的构造
dsize – size of the output image.
flags – combination of interpolation methods (see resize() ) and the optional flag WARP_INVERSE_MAP that means that M is the inverse transformation ( $\texttt{dst}\rightarrow\texttt{src}$ ).
borderMode – pixel extrapolation method (see borderInterpolate()); when borderMode=BORDER_TRANSPARENT , it means that the pixels in the destination image corresponding to the “outliers” in the source image are not modified by the function.
borderValue – value used in case of a constant border; by default, it is 0.

下面介绍一些典型的仿射变换：

（1）平移，将每一点移到到（x+t , y+t)，变换矩阵为

$\begin{bmatrix} 1 &0 &t_{x} \\ 0& 1 &t_{y} \\ 0& 0& 1 \end{bmatrix}$

(2)缩放变换 将每一点的横坐标放大或缩小s_x倍，纵坐标放大（缩小）到s_y倍，变换矩阵为

$\begin{bmatrix} 1 &0 &s_{x} \\ 0& 1 &s_{y} \\ 0& 0& 1 \end{bmatrix}$

（3）旋转变换原点：目标图形围绕原点顺时针旋转Θ 弧度，变换矩阵为

$\begin{bmatrix} cos\Theta &-sin\Theta &0 \\ sin\Theta & cos\Theta &0 \\ 0& 0& 1 \end{bmatrix}$

(4) 旋转变换：目标图形以（x , y ）为轴心顺时针旋转θ弧度，变换矩阵为

$\begin{bmatrix} cos\Theta &-sin\Theta &x-x\cdot cos\Theta+y\cdot sin\Theta \\ sin\Theta & cos\Theta &y-x\cdot sin\Theta-y\cdot cos\Theta \\ 0& 0& 1 \end{bmatrix}$

相当于两次平移与一次原点旋转变换的复合，即先将轴心（x,y)移到到原点，然后做旋转变换，最后将图片的左上角置为图片的原点,即

$M=\begin{bmatrix} 1 & 0 &-x \\ 0& - & -y\\ 0 &0 & 1 \end{bmatrix} \begin{bmatrix} cos\Theta &-sin\Theta & 0\\ sin\Theta &cos\Theta & 0\\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} 1 &0 &x \\ 0 & 1&y \\ 0&0 &1 \end{bmatrix}$

有的人可能会说为什么这么复杂呢，那是因为在opencv的图像处理中，所有对图像的处理都是从原点进行的，而图像的原点默认为图像的左上角，而我们对图像作旋转处理时一般以图像的中点为轴心，因此就需要做如下处理

如果你觉得这样很麻烦，可以使用opencv中自带的Mat getRotationMatrix2D(Point2f center, double angle, double scale)函数获得变换矩阵M，

　　center:旋转中心

　　angle：旋转弧度，一定要将角度转换成弧度

　　scale:缩放尺度

它得到的矩阵是：

$\begin{bmatrix} \alpha &-\beta &center.x\cdot {(\alpha - 1)}-center.y\cdot \beta \\ \beta &\alpha &center.x\cdot \beta + {(\beta-1)\cdot center.y} \\ 0&0 & 1 \end{bmatrix}$

其中α = scale * cos( angle ) , β = scale * sing( angle ) , ( center.x , center.y ) 表示旋转轴心

但是不得不说opencv的文档以及相关书籍中都把这个矩阵写错了，如下：

$\begin{bmatrix} \alpha &\beta &center.x\cdot {(1-\alpha )}-center.y\cdot \beta \\ -\beta &\alpha &center.x\cdot \beta + {(1-\alpha)\cdot center.y} \\ 0&0 & 1 \end{bmatrix}$

建议大家自己通过下式验证一下，即首先将轴心(x,y)移到原点，然后做旋转平绽放变换，最后再将图像的左上角转换为原点

$M= \begin{bmatrix} 1 & 0 &-x \\ 0&1 & -y\\ 0& 0 &1 \end{bmatrix} \begin{bmatrix} cos\Theta &-sin\Theta &0 \\ sin\Theta & cos\Theta & 0\\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} s & 0 & 0\\ 0&s & 0\\ 0 &0 & 1 \end{bmatrix} \begin{bmatrix} 1 & 0 & x\\ 0 & 1 &y \\ 0& 0 & 1 \end{bmatrix}$

没有去研究该函数的源码，不晓得源码中到底怎么写的，但是在别人的博客中看到这个函数貌似需要修正

opencv中还有一个函数：Mat getAffineTransform(InputArray src, InputArray dst)¶

　它通过三组点对就可以获得它们之间的仿射变换，如果我们在一组图像变换中知道变换后的三组点，那么我们就可以利用该函数求得变换矩阵，然后对整张图片进行仿射变换

还有一种与仿射变换经常混淆的变换为透视变换，透视变换需要四组点对才能确定变换矩阵，由于仿射变换保持“平直性”与“平行性”，因此只需要三组点对，而透视变换没有这种约束，故需要四组点对

posted @ 2014-10-27 20:08 dupuleng 阅读(35352) 评论(3) 编辑收藏举报

刷新页面返回顶部