![]() |
University of Murcia, Spain ![]() |
Projective GeometryFunctions related to Projective Geometry.
More... |
Functions | |
QVMatrix | ComputeProjectiveHomography (const QList< QPair< QPointF, QPointF > > &matchings) |
Obtains a planar homography from two lists of corresponding points. | |
QVMatrix | ComputeEuclideanHomography (const QList< QPair< QPointF, QPointF > > &matchings) |
QPointF | ApplyHomography (const QVMatrix &homography, const QPointF &point) |
Obtains the homogenized coordinates of a point mapped using an homography transformation. | |
QList< QPointF > | ApplyHomography (const QVMatrix &homography, const QList< QPointF > &sourcePoints) |
Obtains the homogenized coordinates of a set of points mapped using an homography transformation. | |
double | HomographyTestError (const QVMatrix &homography) |
Function to test if a 3x3 matrix corresponds to an homography. | |
void | GetExtrinsicCameraMatrixFromHomography (const QVMatrix &K, const QVMatrix &H, QVMatrix &M4x4) |
void | GetDirectIntrinsicCameraMatrixFromHomography (const QVMatrix &H, QVMatrix &K) |
Diagonal intrinsic camera matrix calibration. | |
void | CalibrateCameraFromPlanarHomography (const QVMatrix &H, QVMatrix &K, QVMatrix &Rt) |
Diagonal intrinsic camera matrix calibration. | |
void | GetPinholeCameraIntrinsicsFromPlanarHomography (const QVMatrix &H, QVMatrix &K, const int iterations=100, const double maxGradientNorm=1e-3, const double step=0.01, const double tol=1e-4) |
Obtains the intrinsic camera matrix from a planar homography. |
Projective Geometry models the principles of perspective, and provides important results for Computer Vision.
QList<QPointF> ApplyHomography | ( | const QVMatrix & | homography, | |
const QList< QPointF > & | sourcePoints | |||
) |
Obtains the homogenized coordinates of a set of points mapped using an homography transformation.
This function applies a transformation homography over a set of points. For further info about this mapping see function ApplyHomography(const QVMatrix &, const QPointF &).
homography | The homography transformation matrix | |
sourcePoints | Points to apply the homography transformation |
Definition at line 297 of file qvprojective.cpp.
QPointF ApplyHomography | ( | const QVMatrix & | homography, | |
const QPointF & | point | |||
) |
Obtains the homogenized coordinates of a point mapped using an homography transformation.
This function obtains the point which holds the following equation:
Where and
are the input transformation homography and point respectively. Matrix
must be of size
.
homography | The homography transformation matrix | |
point | Point to apply the homography transformation |
Definition at line 286 of file qvprojective.cpp.
Diagonal intrinsic camera matrix calibration.
This function obtains a diagonal intrinsic camera matrix, consisting on the focal distance only. This matrix is such that
Where is a rotation matrix. This matrix is obtained with a minimization process, so its result is more precise than that obtained with GetDirectIntrinsicCameraMatrixFromHomography function.
Definition at line 343 of file qvprojective.cpp.
QVMatrix ComputeEuclideanHomography | ( | const QList< QPair< QPointF, QPointF > > & | matchings | ) |
This function obtains an euclidean mapping between the source an destination locations of a set of point matchings. The mapping is returned as a matrix which can be multiplied to the source location of each mapping point in homogeneous coordinates, obtaining a location close to the destination location from the matching.
Definition at line 229 of file qvprojective.cpp.
QVMatrix ComputeProjectiveHomography | ( | const QList< QPair< QPointF, QPointF > > & | matchings | ) |
Obtains a planar homography from two lists of corresponding points.
This function obtains the planar homography which maps the location of the points from a source plane to a destination plane. One of the planes can be the image plane as well.
The function returns the matrix corresponding to the planar homography, from a list of four or more point correspondences between the location of those points at the source plane and their location in the destination plane.
This function can be used to compute a rectification homography, which can be used to obtain a frontal view of any planar figure appearing in an image. The following code and images shows this effect:
Image on the left is the original picture obtained with a pin-hole camera. Image on the right was obtained aplying a rectification homography to the pixels in the original image. This rectification homography is obtained using the function ComputeProjectiveHomography with the four point correspondences between the points in the images featured with yellow dots:
QList< QPair<QPointF, QPointF> > matchings; matchings.append(QPair<QPointF, QPointF>(QPointF(-171,109), QPointF(-100,+100))); matchings.append(QPair<QPointF, QPointF>(QPointF(-120,31), QPointF(-100,-100))); matchings.append(QPair<QPointF, QPointF>(QPointF(117,53), QPointF(+100,-100))); matchings.append(QPair<QPointF, QPointF>(QPointF(11,115), QPointF(+100,+100))); const QVMatrix H = ComputeProjectiveHomography(matchings);
The obtained matrix $H$ holds the following equation:
Also for each point matching contained in the matchings list, the following equation holds.
or
In C++ source code:
QPair<QPointF, QPointF> matching; foreach(matching, matchings) { const QVVector p_0 = QVVector::homogeneousCoordinates(matching.first), p_1 = QVVector::homogeneousCoordinates(matching.second); std::cout << "p1 x H p0 =" << (p_1 ^ H*p_0) << std::endl; }
The printed values must be close to 0.
The location at the rectified version of an image corresponding to pixel
at the original version of the image can be obtained given the rectification homography
using the following C++ code:
QPointF p_1 = H * QVVector::homogeneousCoordinates(p_0);
Or the ApplyHomography functions.
matchings | list of point matchings from the original location, to the destination location. |
Definition at line 90 of file qvprojective.cpp.
Diagonal intrinsic camera matrix calibration.
This function obtains a diagonal intrinsic camera matrix, consisting on the focal distance only. This matrix is such that
Where is a rotation matrix. This function returns a direct approximation for the
matrix. For a more precise approach, use function GetIntrinsicCameraMatrixFromHomography.
Definition at line 327 of file qvprojective.cpp.
void GetPinholeCameraIntrinsicsFromPlanarHomography | ( | const QVMatrix & | H, | |
QVMatrix & | K, | |||
const int | iterations = 100 , |
|||
const double | maxGradientNorm = 1e-3 , |
|||
const double | step = 0.01 , |
|||
const double | tol = 1e-4 | |||
) |
Obtains the intrinsic camera matrix from a planar homography.
This functions obtains the intrinsic calibration matrix corresponding to a simple pinhole camera model. The intrinsic camera matrix has only one free parameter, related to the focal distance of the camera:
This function should receive a planar homography corresponding to the mapping of a set of know template points, to an image frame captured with the camera containing a view of that template.
The following is an example of a full intrinsic camera calibration, knowing a set of point matchings between the template image and the input image read from the camera:
[...] QList< QPair<QPointF, QPointF> > matchings; matchings.append(QPair<QPointF, QPointF>(QPointF(-171,109), QPointF(-100,+100))); matchings.append(QPair<QPointF, QPointF>(QPointF(-120,31), QPointF(-100,-100))); matchings.append(QPair<QPointF, QPointF>(QPointF(117,53), QPointF(+100,-100))); matchings.append(QPair<QPointF, QPointF>(QPointF(11,115), QPointF(+100,+100))); QVMatrix H, K; H = ComputeProjectiveHomography(matchings); GetPinholeCameraIntrinsicsFromPlanarHomography(H, K); [...]
For further understanding of the planar homography calibration process, see ComputeProjectiveHomography function documentation.
H | Input planar homography. | |
K | Matrix to store the intrinsic camera matrix | |
iterations | Cumber of iterations to perform camera calibration | |
maxGradientNorm | Minimal value of the gradient size (norm 2) to stop the minimization when reached. | |
step | Corresponds to parameter step for the gsl_multimin_fdfminimizer_set function. | |
tol | Corresponds to parameter tol for the gsl_multimin_fdfminimizer_set function. |
Definition at line 400 of file qvprojective.cpp.
double HomographyTestError | ( | const QVMatrix & | homography | ) |
Function to test if a 3x3 matrix corresponds to an homography.
Every matrix which holds a perspective deformation from one plane to another should validate some constrains. The most important is that its two first columns should be perpendicular, and with a similar size, because for the matrix to correspond to an homography, they should be contained in a base of a rotated coordinate system. Given the following homography matrix:
The error value returned by this function for it will be:
That error is a measure of the difference of sizes between the norm of the two column vectors of the homography, corresponding to the two first columns of the rotation matrix, and their dot product. When both values are close to zero, the matrix corresponds to an homography, and it won't otherwise.
A good method to stablish that a matrix corresponds to an homography using this function, can be done testing the return value with a threshold of aproximately 0.3. If the return value of this function for a matrix is lower than this threshold, that matrix is likely to correspond to an homography, and is not likely to correspond otherwise.
homography | a possible homography matrix. |
Definition at line 305 of file qvprojective.cpp.