Я написал простой скрипт для проецирования 3D-точек на изображение на основе внутренних и внешних характеристик камеры. Но когда у меня есть камера в начале координат, указывающая вниз по оси Z, а 3D указывает дальше по оси Z, она кажется позади камеры, а не перед ней. Вот мой скрипт, я проверял его столько раз.
import numpy as np
def project(point, P):
Hp = P.dot(point)
if Hp[-1] < 0:
print 'Point is behind camera'
Hp = Hp / Hp[-1]
print Hp[0][0], 'x', Hp[1][0]
return Hp[0][0], Hp[1][0]
if __name__ == '__main__':
# Rc and C are the camera orientation and location in world coordinates
# Camera posed at origin pointed down the negative z-axis
Rc = np.eye(3)
C = np.array([0, 0, 0])
# Camera extrinsics
R = Rc.T
t = -R.dot(C).reshape(3, 1)
# The camera projection matrix is then:
# P = K [ R | t] which projects 3D world points
# to 2D homogenous image coordinates.
# Example intrinsics dont really matter ...
K = np.array([
[2000, 0, 2000],
[0, 2000, 1500],
[0, 0, 1],
])
# Sample point in front of camera
# i.e. further down the negative x-axis
# should project into center of image
point = np.array([[0, 0, -10, 1]]).T
# Project point into the camera
P = K.dot(np.hstack((R, t)))
# But when projecting it appears to be behind the camera?
project(point,P)
Единственное, о чем я могу думать, это то, что матрица вращения идентификации не соответствует камере, указывающей вниз по отрицательной оси z с вектором вверх в направлении положительной оси y. Но я не понимаю, как это может быть не так, например, я построил Rc из функции, такой как gluLookAt, и дал ему камеру в начале координат, направленную вниз по отрицательной оси z, я получил бы матрицу идентичности.