OpenCV C++(三)---- 几何变换

几何变换:放大缩小、旋转等操作,改变了原图各区域的空间关系。

  • 首先, 需要一个算法实现空间坐标变换, 用它描述每个像素如何从初始位置移动到终止位置;
  • 其次, 还需要一个插值算法完成输出图像的每个像素的灰度值。

3.1 仿射变换

image.png
image.png
image.png

通常称A为仿射变换矩阵。

3.1.1、平移

image.png

其中,若tx> 0,则表示沿x轴正方向移动;若tx< 0,则表示沿x轴负方向移动。ty与之类似。

3.1.2、放缩

image.png

以上介绍的是以原点(0, 0)为中心的缩放变换。
其余情况可以将该变换过程理解为先将原点平移到中心点, 再以原点为中心进行缩放, 然后移回坐标原点, 用矩阵形式可以表示为:

image.png

注:等式右边的计算是从右向左进行的。

3.1.3、旋转

以原点为中心,顺时针旋转α

image.png

以原点为中心,逆时针旋转α

image.png

如果(x,y)绕任意一点(xo,yo)逆时针旋转α,则首先将原点移到旋转中心,然后绕原点旋转,最后移回坐标原点, 即:

image.png

3.1.4、计算仿射矩阵

1、方程法

仿射变换矩阵有六个未知数, 所以只需要三组对应位置坐标, 构造出由六个方程组成的方程组即可解六个未知数。
(1)坐标保存在Point2f数组中

Point2f src[]={Point2f(0,0),Point2f(200,0),Point2f(0,200)};
Point2f dst[]={Point2f(0,0),Point2f(200,0),Point2f(0,200)};
Mat A=getAffineTransform(src,dst);

需要注意的是,数据类型是, CV_64F而不是CV_32F
(2)将原位置坐标和对应的变换后的坐标保存在Mat中, 每一行代表一个坐 标, 数据类型必须是CV_32F, 否则会报错

Mat src=(Mat_<float<(3,2)<<0,0,200,0,0,200);
Mat dst=(Mat_<float<(3,2)<<0,0,200,0,0,200);
Mat A=getAffineTransform(src,dst);

2、矩阵法
对于使用矩阵相乘法计算仿射矩阵,前提是需要知道基本仿射变换步骤,即如果 (x, y)先缩放再平移, 则变换后的矩阵形式为:

image.png

使用矩阵基本操作可以计算出上面的矩阵,注意等式右边的运算是从右向左进行的。
如果是等比例缩放的, 即sx=sy, 则两个仿射变换矩阵是相等的。对于这种等比例缩放的仿射变换, OpenCV提供了函数:

Mat A=getRotationMatrix2D(Point2f(40,50),30,0.5);

数据类型是CV_64F

3.1.5、插值算法

如何利用已知的整数坐标处的函数值估算非整数坐标处的函数值。

1、最近邻插值
最近邻插值就是从(x, y) 的四个相邻整数坐标中找到离它最近的一个 。
使用最近邻插值方法完成图像几何变换, 输出图像会出现锯齿状外观,对图像放大处理的效果会更明显。 为了得到更好的效果,应使用更多的信息, 而不仅仅使用最近像素的灰度值 。

2、双线性插值

  • 第一步: |x-[x]|是点(x, y) 和([x], [y]) 的水平距离, |[x]+1-x|是点(x, y) 和 ([x]+1, [y]) 的水平距离, 显然0< |x-[x]|< 1, 0< |[x]+1-x|<1且|x-[x]|+| [x]+1-x|=1。 为了表示方便, 记a=|x-[x]|, 通过以下线性关系估计fI 在(x, [y ]) 处的值:
image.png
  • 第二步: 同第一步类似,|x-[x]|是点(x, y)和([x], [y ]+1)的水平距离, |[x]+1- x|是点(x, y)和([x]+1, [y]+1)的水平距离,通过以下线性关系估计fI 在(x, [y]+1)处的值
image.png
  • 第三步: 通过第一步和第二步分别得到了fI 在(x, [y]+1)和(x, [y])处的函数 值,(x, y) 和(x, [y ]) 的垂直距离为|y-[y ]|,(x, y) 和(x, [y]+1)处的垂直距 离为|[y ]+1-y |, 通过以下线性关系估计fI 在(x, y) 处的函数值:
image.png

令b=|y-[y]|, 简化上式得到:

image.png

这样对于非整数坐标处的函数值, 就可以利用它的邻域的四个整数坐标处的函数值 进行插值计算而得到。 当然, 如果先进行两次垂直方向上的插值, 然后再进行水平方向上的插值, 得到的结果是一样的。

image.png
image.png

利用opencv实现仿射变换一般会涉及到warpAffinegetRotationMatrix2D两个函数,其中warpAffine可以实现一些简单的重映射,而getRotationMatrix2D可以获得旋转矩阵。

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

Mat cv::getRotationMatrix2D     (   Point2f     center,
        double      angle,
        double      scale 
    )   

使用函数warpAffine对图像进行缩放, 需要先创建缩放仿射矩阵。为了使用更方便, 对于图像的缩放,OpenCV还提供了另一个函数:

void resize(InputArray src,OutputArray dst,Size dsize,double fx=0,double fy=0,int interpolation=INTER_LINEAR);

//第一种方法:利用warpAffine进行缩放
Mat s=(Mat_<float>(2,3)<<0.5,0,0,0,0.5,0);
Mat dst1;
warpAffine(img,dst1,s,Size(img.cols/2,img.rows/2));
//第二种方法:利用resize缩放
Mat dst2;
resize(img,dst2,Size(img.cols/2,img.rows/2),0.5,0.5);

在OpenCV 3.X中完成图像矩阵顺时针旋转90°、 180°、 270° ,如下:

void rotate(InputArray src,OutputArray dst,int rotateCode)

3.2、投影变换

在对仿射变换的讨论中, 校正物体都是在二维空间中完成的, 如果物体在三维空间中发生了旋转, 那么这种变换通常被称为投影变换。 由于可能出现阴影或者遮挡, 所以此投影变换是很难修正的。 但是如果物体是平面的,那么就能通过二维投影变换对此物体三维变换进行模型化, 这就是专用的二维投影变换,可由如下公式描述:

image.png

与用方程法计算仿射变换矩阵的函数getAffineTransform类似(使用方法也类似,有两种),OpenCV提供了函数:

getPerspectiveTransform(src,dst);

3.3、极坐标变换

通常利用极坐标变换来校正图像中的圆形物体或被包含在圆环中的物体。

3.3.1、将笛卡尔坐标转换为极坐标

Mat x=(Mat_<float>(3,3)<<0,1,2,0,1,2,0,1,2)-1;
Mat y=(Mat_<float>(3,3)<<0,0,0,1,1,1,2,2,2)-1;
Mat r,theta;
cartToPolar(x,y,r,theta,true);

3.3.2、将极坐标转换为笛卡尔坐标

Mat angle=(Mat_<float>(2,2)<<30,31,30,31);
Mat r=(Mat_<float>(2,2)<<10,10,11,11);
Mat x,y;
polatToCart(r,angle,x,y true);
x+=-12;
y+=15;

3.3.3、利用极坐标变换对图像进行变换

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 219,366评论 6 508
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,521评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 165,689评论 0 356
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,925评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,942评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,727评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,447评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,349评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,820评论 1 317
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,990评论 3 337
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,127评论 1 351
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,812评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,471评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,017评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,142评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,388评论 3 373
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,066评论 2 355

推荐阅读更多精彩内容