Python - 从旋转角度对OpenCV进行透视变换

use*_*754 9 python opencv stereo-3d perspectivecamera

我正在努力depth map with OpenCV.我可以获得它,但它是从左相机原点重建的,后者有一点倾斜,如图所示,深度"移位"(深度应该接近,没有水平渐变):

在此输入图像描述

我想以零角度表达它,我尝试使用warp透视函数,如下所示,但我获得了一个空字段...

P = np.dot(cam,np.dot(Transl,np.dot(Rot,A1)))

dst = cv2.warpPerspective(depth, P, (2048, 2048))
Run Code Online (Sandbox Code Playgroud)

用:

#Projection 2D -> 3D matrix
A1 = np.zeros((4,3))
A1[0,0] = 1 
A1[0,2] = -1024
A1[1,1] = 1
A1[1,2] = -1024
A1[3,2] = 1

#Rotation matrice around the Y axis
theta = np.deg2rad(5) 
Rot = np.zeros((4,4))
Rot[0,0] = np.cos(theta)
Rot[0,2] = -np.sin(theta)
Rot[1,1] = 1
Rot[2,0] = np.sin(theta)
Rot[2,2] = np.cos(theta)
Rot[3,3] = 1

#Translation matrix on the X axis 
dist = 0
Transl = np.zeros((4,4))
Transl[0,0] = 1 
Transl[0,2] = dist
Transl[1,1] = 1
Transl[2,2] = 1
Transl[3,3] = 1

#Camera Intrisecs matrix 3D -> 2D
cam = np.concatenate((C1,np.zeros((3,1))),axis=1)
cam[2,2] = 1
P = np.dot(cam,np.dot(Transl,np.dot(Rot,A1)))

dst = cv2.warpPerspective(Z0_0, P, (2048*3, 2048*3))
Run Code Online (Sandbox Code Playgroud)

稍后编辑:

您可以在此处下载32MB字段数据集:https://filex.ec-lille.fr/get? k = cCBoyoV4tbmkzSV5bi6.然后,加载并查看图像:

from matplotlib import pyplot as plt
import numpy as np

img = np.load('testZ0.npy')
plt.imshow(img)
plt.show()
Run Code Online (Sandbox Code Playgroud)

Jer*_*uke 6

我已经有了一个粗略的解决方案。您可以稍后修改它。

我使用 OpenCV 中可用的鼠标处理操作来裁剪给定热图中的感兴趣区域。

(我刚才是不是说我用鼠标裁剪了该区域?)是的,我做到了。要了解有关 OpenCV 中鼠标功能的更多信息,请参阅此。此外,还有许多其他问题可以在这方面为您提供帮助。:)

使用这些函数我能够获得以下内容:

在此输入图像描述

现在谈谈消除倾斜的问题。我通过获取上面图像的角点并将其用于确定尺寸的“白色”图像来使用单应性原理。我cv2.findHomography()为此使用了该功能。

现在使用cv2.warpPerspective()OpenCV 中的函数,我能够获得以下结果:

在此输入图像描述

现在您可以根据需要缩放此图像。

代码:

我还附上了一些代码片段供您阅读:

#First I created an image of white color of a definite size
back = np.ones((435, 379, 3)) # size
back[:] = (255, 255, 255)     # white color
Run Code Online (Sandbox Code Playgroud)

pts_src接下来我获得了下面倾斜图像上的角点:

pts_src = np.array([[25.0, 2.0],[403.0,22.0],[375.0,436.0],[6.0,433.0]])
Run Code Online (Sandbox Code Playgroud)

我希望将上面的点映射到下面给出的点“pts_dst”:

pts_dst = np.array([[2.0, 2.0], [379.0, 2.0], [379.0, 435.0],[2.0, 435.0]])
Run Code Online (Sandbox Code Playgroud)

现在我使用单应性原理:

h, status = cv2.findHomography(pts_src, pts_dst)
Run Code Online (Sandbox Code Playgroud)

最后,我使用透视变换将原始图像映射到白色图像。

fin = cv2.warpPerspective(img, h, (back.shape[1],back.shape[0]))
# img -> original tilted image.
# back -> image of white color.
Run Code Online (Sandbox Code Playgroud)

希望这可以帮助!我也从这个问题中学到了很多东西。

注意:输入到“cv2.findHomography()”的点必须位于float. 有关同应性的更多信息,请访问此页面