详解Python计算机视觉 图像扭曲(仿射扭曲)


Posted in Python onMarch 27, 2019

对图像块应用仿射变换,我们将其称为图像扭曲(或者仿射扭曲)。该操作不仅经常应用在计算机图形学中,而且经常出现在计算机视觉算法中。

一、仿射变换原理

仿射变换能够保持图像的“平直性”,包括旋转,缩放,平移,错切操作。对于三个点,仿射变换可以将一副图像进行扭曲,使得三对对应点对可以完美地匹配上。仿射变换具有6个自由度,有三个对应点对可以给出6个约束条件(对于这三个对应点对,x和y坐标必须都要匹配)

仿射变换是在几何上定义为两个向量空间之间的一个仿射变换或者仿射映射。由一个非奇异的线性变换(运用一次函数进行的变换)接上一个平移变换组成。在有限维的情况,每个仿射变换可以由一个矩阵A和一个向量b给出,它可以写作A和一个附加的列b。一个仿射变换对应于一个矩阵和一个向量的乘法,而仿射变换的复合对应于普通的矩阵乘法,只要加入一个额外的行到矩阵的底下,这一行全部是0除了最右边是一个1,而列向量的底下要加上一个1。

详解Python计算机视觉 图像扭曲(仿射扭曲)

二、图像中的图像

仿射扭曲简单的一个例子是,将图像或者图像的一部分放置在另一幅图像中,使得它们能够和指定的区域或者标记物对齐。
以下是代码:

from PCV.geometry import warp, homography
from PIL import Image
from pylab import *
from scipy import ndimage

# 仿射扭曲im1到im2的例子

im1 = array(Image.open('jida.jpg').convert('L'))
im2 = array(Image.open('beijing.jpg').convert('L'))
# 选定一些目标点
tp = array([[420,830,830,420],[400,350,1060,1000],[1,1,1,1]])
#标记物的坐标tp是用齐次坐标意义下的坐标表示的
#array为 第一张图片贴到第二张图片的四个角点的坐标
#tp = array([[675,826,826,677],[55,52,281,277],[1,1,1,1]])
im3 = warp.image_in_image(im1,im2,tp)
#image_in_image()函数的输入参数为两幅图像和一个坐标
figure()
gray()
subplot(141)
axis('off')
imshow(im1)
subplot(142)
axis('off')
imshow(im2)
subplot(143)
axis('off')
imshow(im3)

# 选定im1角上的一些点
m,n = im1.shape[:2]
fp = array([[0,m,m,0],[0,0,n,n],[1,1,1,1]])
# 第一个三角形
tp2 = tp[:,:3]
fp2 = fp[:,:3]
# 计算H
H = homography.Haffine_from_points(tp2,fp2)
im1_t = ndimage.affine_transform(im1,H[:2,:2],
(H[0,2],H[1,2]),im2.shape[:2])
# 三角形的alpha
alpha = warp.alpha_for_triangle(tp2,im2.shape[0],im2.shape[1])
im3 = (1-alpha)*im2 + alpha*im1_t
# 第二个三角形
tp2 = tp[:,[0,2,3]]
fp2 = fp[:,[0,2,3]]
# 计算H
H = homography.Haffine_from_points(tp2,fp2)
#Haffine_from_points()返回给定对应点对的最优仿射变换
im1_t = ndimage.affine_transform(im1,H[:2,:2],
(H[0,2],H[1,2]),im2.shape[:2])
# 三角形的alpha
alpha = warp.alpha_for_triangle(tp2,im2.shape[0],im2.shape[1])
im4 = (1-alpha)*im3 + alpha*im1_t
subplot(144)
imshow(im4)
axis('off')
show()

详解Python计算机视觉 图像扭曲(仿射扭曲)

在实验过程中,也出现了一些问题。

问题一:

详解Python计算机视觉 图像扭曲(仿射扭曲)

为解决问题,先去PCV\PCV-master\PCV\geometry中找到wary.py和homegraphy.py文件,把print后面的语句都加上括号。但是在检查了好几遍括号以后,问题都没有解决,所以重装了一次pcv以后,异常就解决了。(重装pcv的方法在前面的博客有记录)

详解Python计算机视觉 图像扭曲(仿射扭曲)

问题二:

详解Python计算机视觉 图像扭曲(仿射扭曲)

因为是matplotlib发生了异常,所以我把matplotlib卸载了(在终端输入pip uninstall matplotlib),然后重装matplotlib,我在终端输入pip install matplotlib以后,安装无法成功,显示Could not install packages due to an EnvironmentError: [Errno 13] Permission denied。

这时候解决方法如下,输入:pip install matplotlib --user便可成功安装(注意:是有两个-)。

详解Python计算机视觉 图像扭曲(仿射扭曲)

详解Python计算机视觉 图像扭曲(仿射扭曲)

解决以上两个问题后,代码可正常运行。

以上所述是小编给大家介绍的Python计算机视觉 图像扭曲(仿射扭曲)详解整合,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对三水点靠木网站的支持!

Python 相关文章推荐
Python数组条件过滤filter函数使用示例
Jul 22 Python
在Python中使用pngquant压缩png图片的教程
Apr 09 Python
JPype实现在python中调用JAVA的实例
Jul 19 Python
python网络爬虫之如何伪装逃过反爬虫程序的方法
Nov 23 Python
使用apidocJs快速生成在线文档的实例讲解
Feb 07 Python
python2.6.6如何升级到python2.7.14
Apr 08 Python
对pycharm代码整体左移和右移缩进快捷键的介绍
Jul 16 Python
Python Numpy数组扩展repeat和tile使用实例解析
Dec 09 Python
Django 拼接两个queryset 或是两个不可以相加的对象实例
Mar 28 Python
Python实现电视里的5毛特效实例代码详解
May 15 Python
Pymysql实现往表中插入数据过程解析
Jun 02 Python
详解Python函数print用法
Jun 18 Python
python实现向微信用户发送每日一句 python实现微信聊天机器人
Mar 27 #Python
Pandas读写CSV文件的方法示例
Mar 27 #Python
使用Python的SymPy库解决数学运算问题的方法
Mar 27 #Python
超简单使用Python换脸实例
Mar 27 #Python
python爬虫爬取微博评论案例详解
Mar 27 #Python
Python实现查找字符串数组最长公共前缀示例
Mar 27 #Python
详解python中@的用法
Mar 27 #Python
You might like
php smarty模版引擎中的缓存应用
2009/12/02 PHP
php中在PDO中使用事务(Transaction)
2011/05/14 PHP
php中全局变量global的使用演示代码
2011/05/18 PHP
php自动识别文件编码并转换为UTF-8的方法
2014/06/12 PHP
PHP中使用匿名函数操作数据库的例子
2014/11/17 PHP
php启用sphinx全文搜索的实现方法
2014/12/24 PHP
PHP自定义序列化接口Serializable用法分析
2017/12/29 PHP
js获取url参数值的两种方式
2013/09/10 Javascript
html5 canvas js(数字时钟)实例代码
2013/12/23 Javascript
什么是cookie?js手动创建和存储cookie
2014/05/27 Javascript
JavaScript实现生成GUID(全局统一标识符)
2014/09/05 Javascript
JS拖拽组件学习使用
2016/01/19 Javascript
如何用JS/HTML将时间戳转换为“xx天前”的形式
2017/02/06 Javascript
JS日程管理插件FullCalendar中文说明文档
2017/02/06 Javascript
Java设计中的Builder模式的介绍
2018/03/22 Javascript
30分钟快速入门掌握ES6/ES2015的核心内容(下)
2018/04/18 Javascript
解决vue-cli脚手架打包后vendor文件过大的问题
2018/09/27 Javascript
nodejs基础之多进程实例详解
2018/12/27 NodeJs
JS使用new操作符创建对象的方法分析
2019/05/30 Javascript
深入解析koa之中间件流程控制
2019/06/17 Javascript
Vue.js实现大转盘抽奖总结及实现思路
2019/10/09 Javascript
[02:17]《辉夜杯》TRG战队巡礼
2015/10/26 DOTA
Python使用Srapy框架爬虫模拟登陆并抓取知乎内容
2016/07/02 Python
Python闭包执行时值的传递方式实例分析
2018/06/04 Python
使用Python向DataFrame中指定位置添加一列或多列的方法
2019/01/29 Python
Django 对IP访问频率进行限制的例子
2019/08/30 Python
Python控制鼠标键盘代码实例
2020/12/08 Python
德国网上超市:myTime.de
2019/08/26 全球购物
毕业生的自我评价分享
2013/12/18 职场文书
2014年民政局关于保密工作整改措施
2014/09/19 职场文书
授权委托书范本(单位)
2014/09/28 职场文书
导游词之西安大清真寺
2019/12/17 职场文书
python+opencv实现视频抽帧示例代码
2021/06/11 Python
vite+vue3.0+ts+element-plus快速搭建项目的实现
2021/06/24 Vue.js
Python使用openpyxl模块处理Excel文件
2022/06/05 Python
java实现面板之间切换功能
2022/06/10 Java/Android