vector_to_fundamental_matrixT_vector_to_fundamental_matrixVectorToFundamentalMatrixVectorToFundamentalMatrix (Operator)

Name

vector_to_fundamental_matrixT_vector_to_fundamental_matrixVectorToFundamentalMatrixVectorToFundamentalMatrix — Compute the fundamental matrix given a set of image point correspondences and reconstruct 3D points.

Signature

Description

For a stereo configuration with unknown camera parameters the geometric relation between the two images is defined by the fundamental matrix. The operator vector_to_fundamental_matrixvector_to_fundamental_matrixVectorToFundamentalMatrixVectorToFundamentalMatrixVectorToFundamentalMatrix determines the fundamental matrix FMatrixFMatrixFMatrixFMatrixFMatrix from given point correspondences (Rows1Rows1Rows1Rows1rows1,Cols1Cols1Cols1Cols1cols1), (Rows2Rows2Rows2Rows2rows2,Cols2Cols2Cols2Cols2cols2), that fulfill the epipolar constraint:

Note the column/row ordering in the point coordinates: since the fundamental matrix encodes the projective relation between two stereo images embedded in 3D space, the x/y notation must be compliant with the camera coordinate system. Therefore, (x,y) coordinates correspond to (column,row) pairs.

For a general relative orientation of the two cameras the minimum number of required point correspondences is eight. Then, MethodMethodMethodMethodmethod is chosen to be 'normalized_dlt'"normalized_dlt""normalized_dlt""normalized_dlt""normalized_dlt" or 'gold_standard'"gold_standard""gold_standard""gold_standard""gold_standard". If left and right camera are identical and the relative orientation between them is a pure translation then choose MethodMethodMethodMethodmethod equal to 'trans_normalized_dlt'"trans_normalized_dlt""trans_normalized_dlt""trans_normalized_dlt""trans_normalized_dlt" or 'trans_gold_standard'"trans_gold_standard""trans_gold_standard""trans_gold_standard""trans_gold_standard". In this special case the minimum number of correspondences is only two. The typical application of the motion beeing a pure translation is that of a single fixed camera looking onto a moving conveyor belt.

The fundamental matrix is determined by minimizing a cost function. To minimize the respective error different algorithms are available, and the user can choose between the direct-linear-transformation ('normalized_dlt') and the gold-standard-algorithm ('gold_standard'). Like the motion case, the algorithm can be selected with the parameter MethodMethodMethodMethodmethod. For MethodMethodMethodMethodmethod = 'normalized_dlt'"normalized_dlt""normalized_dlt""normalized_dlt""normalized_dlt" or 'trans_normalized_dlt'"trans_normalized_dlt""trans_normalized_dlt""trans_normalized_dlt""trans_normalized_dlt", a linear algorithm minimizes an algebraic error based on the above epipolar constraint. This algorithm offers a good compromise between speed and accuracy. For MethodMethodMethodMethodmethod = 'gold_standard'"gold_standard""gold_standard""gold_standard""gold_standard" or 'trans_gold_standard'"trans_gold_standard""trans_gold_standard""trans_gold_standard""trans_gold_standard", a mathematically optimal, but slower optimization is used, which minimizes the geometric backprojection error of reconstructed projective 3D points. In this case, in addition to the fundamental matrix its covariance matrix CovFMatCovFMatCovFMatCovFMatcovFMat is output, along with the projective coordinates (XXXXx,YYYYy,ZZZZz,WWWWw) of the reconstructed points and their covariances CovXYZWCovXYZWCovXYZWCovXYZWcovXYZW. Let n be the number of points. Then the concatenated covariances are stored in a 16xn tuple.

If an optimal gold-standard-algorithm is chosen the covariances of the image points (CovRR1CovRR1CovRR1CovRR1covRR1, CovRC1CovRC1CovRC1CovRC1covRC1, CovCC1CovCC1CovCC1CovCC1covCC1, CovRR2CovRR2CovRR2CovRR2covRR2, CovRC2CovRC2CovRC2CovRC2covRC2, CovCC2CovCC2CovCC2CovCC2covCC2) can be incorporated in the computation. They can be provided for example by the operator points_foerstnerpoints_foerstnerPointsFoerstnerPointsFoerstnerPointsFoerstner. If the point covariances are unknown, which is the default, empty tuples are input. In this case the optimization algorithm internally assumes uniform and equal covariances for all points.

The value ErrorErrorErrorErrorerror indicates the overall quality of the optimization procedure and is the mean euclidian distance in pixels between the points and their corresponding epipolar lines.

If the correspondence between the points are not known, match_fundamental_matrix_ransacmatch_fundamental_matrix_ransacMatchFundamentalMatrixRansacMatchFundamentalMatrixRansacMatchFundamentalMatrixRansac should be used instead.

Execution Information

• Multithreading type: reentrant (runs in parallel with non-exclusive operators).
• Multithreading scope: global (may be called from any thread).
• Processed without parallelization.

Parameters

Possible Predecessors

match_fundamental_matrix_ransacmatch_fundamental_matrix_ransacMatchFundamentalMatrixRansacMatchFundamentalMatrixRansacMatchFundamentalMatrixRansac

Possible Successors

gen_binocular_proj_rectificationgen_binocular_proj_rectificationGenBinocularProjRectificationGenBinocularProjRectificationGenBinocularProjRectification

Alternatives

vector_to_essential_matrixvector_to_essential_matrixVectorToEssentialMatrixVectorToEssentialMatrixVectorToEssentialMatrix, vector_to_rel_posevector_to_rel_poseVectorToRelPoseVectorToRelPoseVectorToRelPose

References

Richard Hartley, Andrew Zisserman: “Multiple View Geometry in Computer Vision”; Cambridge University Press, Cambridge; 2000.
Olivier Faugeras, Quang-Tuan Luong: “The Geometry of Multiple Images: The Laws That Govern the Formation of Multiple Images of a Scene and Some of Their Applications”; MIT Press, Cambridge, MA; 2001.

Module

3D Metrology