透视单应性的问题

时间:2013-01-17 20:50:52

标签: python image-processing opencv perspectivecamera

我试图在OpenCV和Python中使用透视扭曲来纠正图像。我知道相机的方向(X,Y,Z角度)参考需要翘曲的平面。我知道最简单的方法是根据已知点计算单应矩阵,但是当这些信息不可用时,我试图做同样的事情。我正在使用的代码创建一个旋转矩阵,然后组合平移矩阵和内在矩阵。目前,该代码仅适用于对z轴的操作。对x和y轴的任何操纵都会导致图像的奇怪失真。我的代码基于以下帖子底部的答案:Perspective Warping in OpenCV based on know camera orientation

附上原始图像和标准单应法中的扭曲图像。

Original Image Perspective Warp

from numpy import *
import cv

x = float(0)
y = float(5)
z = float(0)
f = 1

im = cv.LoadImage("Homography_test.jpg")
cv.NamedWindow("Distorted")
cv.NamedWindow("undistorted")

h, w = cv.GetSize(im)

x1 = x * (pi / 180)
y1 = y * (pi / 180)
z1 = z * (pi / 180)

# Create a rotation matrix
R_array = array([[x1], [y1], [z1]])
R_Vec = cv.fromarray(R_array)
R = cv.CreateMat(3, 3, cv.CV_64FC1)

cv.Rodrigues2(R_Vec, R)

#Create and combine with translation matrix
Trans_Mat = array([[[1], [0], [-w/2]],
                    [[0], [1], [-h/2]],
                    [[0], [0], [1]]])

Trans_Mat2 = cv.fromarray(Trans_Mat)
R_T_Mat = dot(R, Trans_Mat2)

#Create and combine with camera matrix
Intrinsic_Mat = array([[[f], [0], [w/2]],
                       [[0], [f], [h/2]],
                       [[0], [0], [1]]])

Int_Mat = cv.fromarray(Intrinsic_Mat)
H = dot(Int_Mat, R_T_Mat)
H2 = cv.fromarray(H)

persp = cv.CreateImage(cv.GetSize(im), cv.IPL_DEPTH_8U, 3)
cv.WarpPerspective(im, persp, H2)

cv.ShowImage("Distorted", im)
cv.ShowImage("undistorted", persp)

cv.WaitKey(0)
cv.DestroyWindow("Distorted")
cv.DestroyWindow("undistorted")

1 个答案:

答案 0 :(得分:1)

你似乎错过了我的解决方案中的第4步,即“将图像向下移动z轴”,这就是我以前所做的那条线。

//4
trans(2,2) += image.rows;

我随意选择了image.rows,它对我在示例中使用的旋转给出了很好的比例感。因为你没有这样做你的z坐标被固定为1.我想这是当你围绕x和y旋转时失真的原因,主要是因为当它围绕x或y旋转时图像非常靠近相机透视失真是巨大的。沿z轴向下平移越大,失真应​​该越小。如果您不希望图像在沿z轴向下移动时收缩,只需增加焦距即可。您可能会注意到,在我的示例中,我将焦距设置为image.rows。