单应性矩阵
来源:互联网 发布:苹果电脑mac地址 编辑:程序博客网 时间:2024/05/01 14:12
单应性:
在计算机视觉中:对单应性最感兴趣的部分只是其他意义的一个子集。平面的单应性被定义为从一个平面到另一个平面的投影映射。比如,一个二维平面上的点映射到摄像机成像仪上的映射就是平面单应性的例子。
在一个三维空间中的平面,投影到多个影像之间存在一特定关系,此关系称为 planar homography
H=⎡⎣⎢⎢h11h21h31h12h22h32h13h23h33⎤⎦⎥⎥
矩阵H 会将一幅图像上的一个点的坐标a=(x,y,1) 映射成另一幅图像上的点的坐标b=(x1,y1,1) ,也就是说,我们已知a 和b ,它们是在同一平面上。 则有下面的公式:
b=HaT(1)
即:
⎧⎩⎨⎪⎪x1=h11x+h12y+h13y1=h21x+h22y+h231=h31x+h32y+h33(2)
由上面这个公式中的1=h31x+h32y+h33 可得到下面两个等式
⎧⎩⎨⎪⎪x1=x11=h11x+h12y+h13h31x+h32y+h33y1=y11=h21x+h22y+h23h31x+h32y+h33(3)
⇒
{h11x+h12y+h13=h31xx1+h32yx1+h33x1h21x+h22y+h23=h31xy1+h32yy1+h33y1(4)
⇒
{0=h31xx1+h32yx1+h33x1−(h11x+h12y+h13)0=h31xy1+h32yy1+h33y1−(h21x+h22y+h23)(5)
对于方程(???) ,可写成一个矩阵与一个向量相乘,即:
[−x0−y0−100−x0−y0−1xx1xy1yx1yy1x1y1]h=0(6)
其中,h=[h11,h12,h13,h21,h22,h23,h31,h32,h33]T ,是一个9维的列向量。若令:
A=[−x0−y0−100−x0−y0−1xx1xy1yx1yy1x1y1](7)
则(6) 可以记为
Ah=0(8)
这里的A∈R2×9 。这只是1对点所得到的矩阵A ,若有4对点,则得到的矩阵A∈R8×9 。如何求解向量h 呢?方法很简单,真接对A 进行SVD分解,即
U∗Σ∗VT(9)
然后取V 的最后一列出来作为求解h 。因为矩阵A 是行满秩,即只有一个自由度。
具体实现时,先要得到两幅图,然后在两幅图之间找到4对点的坐标,由此得到矩阵A
在一个三维空间中的平面,投影到多个影像之间存在一特定关系,此关系称为 planar homography
单应性(homography)变换的推导
矩阵的一个重要作用是将空间中的点变换到另一个空间中。这个作用在国内的《线性代数》教学中基本没有介绍。要能形像地理解这一作用,比较直观的方法就是图像变换,图像变换的方法很多,单应性变换是其中一种方法,单应性变换会涉及到单应性矩阵。单应性变换的目标是通过给定的几个点(通常是4对点)来得到单应性矩阵。下面单应性矩阵的推导过程。
矩阵
即:
由上面这个公式中的
对于方程
其中,
则
这里的
然后取
具体实现时,先要得到两幅图,然后在两幅图之间找到4对点的坐标,由此得到矩阵
,然后在matlab中可以这样实现:
[U,S,V]=svd(A);
h=V(:,9);
H= reshape(h,3,3);
由单应性矩阵可以得到仿射变换,还可以在单应性矩阵上做图像拼接。
0 0
- 单应性矩阵、本质矩阵
- 单应性矩阵,本征矩阵,基础矩阵
- 单应性矩阵
- 单应性矩阵
- 四点求解单应性矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 矩阵
- 1101. Quick Sort (25)
- vim tab设置为4个空格
- GsonFromat使用快速生成实体对象
- Python进程和线程操作
- Python学习day5作业
- 单应性矩阵
- JS页面时钟
- Hello! CSDN
- Light OJ 1125 Divisible Group Sums 背包DP 2017/1/23
- tomcat8配置https ARP方式
- 【c++】继承
- 链表List的代码,但是缺少保护
- 在VS2015中使用SDL2.0
- 已知二叉树的后序遍历和中序遍历求前序遍历(二叉树)