283 lines
9.9 KiB
C++
283 lines
9.9 KiB
C++
/* ----------------------------------------------------------------------------
|
|
|
|
* GTSAM Copyright 2010, Georgia Tech Research Corporation,
|
|
* Atlanta, Georgia 30332-0415
|
|
* All Rights Reserved
|
|
* Authors: Frank Dellaert, et al. (see THANKS for the full author list)
|
|
|
|
* See LICENSE for the license information
|
|
|
|
* -------------------------------------------------------------------------- */
|
|
|
|
/**
|
|
* @file triangulation.h
|
|
* @brief Functions for triangulation
|
|
* @date July 31, 2013
|
|
* @author Chris Beall
|
|
*/
|
|
|
|
#pragma once
|
|
|
|
#include <gtsam/geometry/Pose3.h>
|
|
#include <gtsam/geometry/Point2.h>
|
|
#include <gtsam/geometry/PinholeCamera.h>
|
|
#include <gtsam/inference/Symbol.h>
|
|
#include <gtsam/nonlinear/NonlinearFactorGraph.h>
|
|
#include <gtsam/nonlinear/LevenbergMarquardtOptimizer.h>
|
|
#include <gtsam/slam/ProjectionFactor.h>
|
|
#include <gtsam/slam/PriorFactor.h>
|
|
#include <gtsam_unstable/base/dllexport.h>
|
|
|
|
#include <boost/foreach.hpp>
|
|
#include <boost/assign.hpp>
|
|
#include <boost/assign/std/vector.hpp>
|
|
|
|
#include <vector>
|
|
|
|
namespace gtsam {
|
|
|
|
/// Exception thrown by triangulateDLT when SVD returns rank < 3
|
|
class TriangulationUnderconstrainedException: public std::runtime_error {
|
|
public:
|
|
TriangulationUnderconstrainedException() :
|
|
std::runtime_error("Triangulation Underconstrained Exception.") {
|
|
}
|
|
};
|
|
|
|
/// Exception thrown by triangulateDLT when landmark is behind one or more of the cameras
|
|
class TriangulationCheiralityException: public std::runtime_error {
|
|
public:
|
|
TriangulationCheiralityException() :
|
|
std::runtime_error(
|
|
"Triangulation Cheirality Exception: The resulting landmark is behind one or more cameras.") {
|
|
}
|
|
};
|
|
|
|
/**
|
|
* DLT triangulation: See Hartley and Zisserman, 2nd Ed., page 312
|
|
* @param projection_matrices Projection matrices (K*P^-1)
|
|
* @param measurements 2D measurements
|
|
* @param rank_tol SVD rank tolerance
|
|
* @return Triangulated Point3
|
|
*/
|
|
GTSAM_UNSTABLE_EXPORT Point3 triangulateDLT(const std::vector<Matrix>& projection_matrices, const std::vector<Point2>& measurements, double rank_tol);
|
|
|
|
// Frank says: putting priors on poses and then optimizing is a terrible idea: we turn a 3dof problem into a much more difficult problem
|
|
// We should have a projectionfactor that knows pose is fixed
|
|
|
|
///
|
|
/**
|
|
* Create a factor graph with projection factors from poses and one calibration
|
|
* @param poses Camera poses
|
|
* @param sharedCal shared pointer to single calibration object
|
|
* @param measurements 2D measurements
|
|
* @param landmarkKey to refer to landmark
|
|
* @param initialEstimate
|
|
* @return graph and initial values
|
|
*/
|
|
template<class CALIBRATION>
|
|
std::pair<NonlinearFactorGraph, Values> triangulationGraph(
|
|
const std::vector<Pose3>& poses, boost::shared_ptr<CALIBRATION> sharedCal,
|
|
const std::vector<Point2>& measurements, Key landmarkKey,
|
|
const Point3& initialEstimate) {
|
|
Values values;
|
|
values.insert(landmarkKey, initialEstimate); // Initial landmark value
|
|
NonlinearFactorGraph graph;
|
|
static SharedNoiseModel unit2(noiseModel::Unit::Create(2));
|
|
static SharedNoiseModel prior_model(noiseModel::Isotropic::Sigma(6, 1e-6));
|
|
for (size_t i = 0; i < measurements.size(); i++) {
|
|
const Pose3& pose_i = poses[i];
|
|
graph.push_back(GenericProjectionFactor<Pose3, Point3, CALIBRATION> //
|
|
(measurements[i], unit2, i, landmarkKey, sharedCal));
|
|
graph.push_back(PriorFactor<Pose3>(i, pose_i, prior_model));
|
|
values.insert(i, pose_i);
|
|
}
|
|
return std::make_pair(graph, values);
|
|
}
|
|
|
|
/**
|
|
* Create a factor graph with projection factors from pinhole cameras
|
|
* (each camera has a pose and calibration)
|
|
* @param cameras pinhole cameras
|
|
* @param measurements 2D measurements
|
|
* @param landmarkKey to refer to landmark
|
|
* @param initialEstimate
|
|
* @return graph and initial values
|
|
*/
|
|
template<class CALIBRATION>
|
|
std::pair<NonlinearFactorGraph, Values> triangulationGraph(
|
|
const std::vector<PinholeCamera<CALIBRATION> >& cameras,
|
|
const std::vector<Point2>& measurements, Key landmarkKey,
|
|
const Point3& initialEstimate) {
|
|
Values values;
|
|
values.insert(landmarkKey, initialEstimate); // Initial landmark value
|
|
NonlinearFactorGraph graph;
|
|
static SharedNoiseModel unit2(noiseModel::Unit::Create(2));
|
|
static SharedNoiseModel prior_model(noiseModel::Isotropic::Sigma(6, 1e-6));
|
|
for (size_t i = 0; i < measurements.size(); i++) {
|
|
const PinholeCamera<CALIBRATION>& camera_i = cameras[i];
|
|
boost::shared_ptr<CALIBRATION> // Seems wasteful to create new object
|
|
sharedCal(new CALIBRATION(camera_i.calibration()));
|
|
graph.push_back(GenericProjectionFactor<Pose3, Point3, CALIBRATION> //
|
|
(measurements[i], unit2, i, landmarkKey, sharedCal));
|
|
const Pose3& pose_i = camera_i.pose();
|
|
graph.push_back(PriorFactor<Pose3>(i, pose_i, prior_model));
|
|
values.insert(i, pose_i);
|
|
}
|
|
return std::make_pair(graph, values);
|
|
}
|
|
|
|
///
|
|
/**
|
|
* Optimize for triangulation
|
|
* @param graph nonlinear factors for projection
|
|
* @param values initial values
|
|
* @param landmarkKey to refer to landmark
|
|
* @return refined Point3
|
|
*/
|
|
GTSAM_UNSTABLE_EXPORT Point3 optimize(const NonlinearFactorGraph& graph, const Values& values, Key landmarkKey);
|
|
|
|
/**
|
|
* Given an initial estimate , refine a point using measurements in several cameras
|
|
* @param poses Camera poses
|
|
* @param sharedCal shared pointer to single calibration object
|
|
* @param measurements 2D measurements
|
|
* @param initialEstimate
|
|
* @return refined Point3
|
|
*/
|
|
template<class CALIBRATION>
|
|
Point3 triangulateNonlinear(const std::vector<Pose3>& poses,
|
|
boost::shared_ptr<CALIBRATION> sharedCal,
|
|
const std::vector<Point2>& measurements, const Point3& initialEstimate) {
|
|
|
|
// Create a factor graph and initial values
|
|
Values values;
|
|
NonlinearFactorGraph graph;
|
|
boost::tie(graph, values) = triangulationGraph(poses, sharedCal, measurements,
|
|
Symbol('p', 0), initialEstimate);
|
|
|
|
return optimize(graph, values, Symbol('p', 0));
|
|
}
|
|
|
|
/**
|
|
* Given an initial estimate , refine a point using measurements in several cameras
|
|
* @param cameras pinhole cameras
|
|
* @param measurements 2D measurements
|
|
* @param initialEstimate
|
|
* @return refined Point3
|
|
*/
|
|
template<class CALIBRATION>
|
|
Point3 triangulateNonlinear(
|
|
const std::vector<PinholeCamera<CALIBRATION> >& cameras,
|
|
const std::vector<Point2>& measurements, const Point3& initialEstimate) {
|
|
|
|
// Create a factor graph and initial values
|
|
Values values;
|
|
NonlinearFactorGraph graph;
|
|
boost::tie(graph, values) = triangulationGraph(cameras, measurements,
|
|
Symbol('p', 0), initialEstimate);
|
|
|
|
return optimize(graph, values, Symbol('p', 0));
|
|
}
|
|
|
|
/**
|
|
* Function to triangulate 3D landmark point from an arbitrary number
|
|
* of poses (at least 2) using the DLT. The function checks that the
|
|
* resulting point lies in front of all cameras, but has no other checks
|
|
* to verify the quality of the triangulation.
|
|
* @param poses A vector of camera poses
|
|
* @param sharedCal shared pointer to single calibration object
|
|
* @param measurements A vector of camera measurements
|
|
* @param rank tolerance, default 1e-9
|
|
* @param optimize Flag to turn on nonlinear refinement of triangulation
|
|
* @return Returns a Point3
|
|
*/
|
|
template<class CALIBRATION>
|
|
Point3 triangulatePoint3(const std::vector<Pose3>& poses,
|
|
boost::shared_ptr<CALIBRATION> sharedCal,
|
|
const std::vector<Point2>& measurements, double rank_tol = 1e-9,
|
|
bool optimize = false) {
|
|
|
|
assert(poses.size() == measurements.size());
|
|
if (poses.size() < 2)
|
|
throw(TriangulationUnderconstrainedException());
|
|
|
|
// construct projection matrices from poses & calibration
|
|
std::vector<Matrix> projection_matrices;
|
|
BOOST_FOREACH(const Pose3& pose, poses) {
|
|
projection_matrices.push_back(
|
|
sharedCal->K() * sub(pose.inverse().matrix(), 0, 3, 0, 4));
|
|
}
|
|
|
|
// Triangulate linearly
|
|
Point3 point = triangulateDLT(projection_matrices, measurements, rank_tol);
|
|
|
|
// The n refine using non-linear optimization
|
|
if (optimize)
|
|
point = triangulateNonlinear(poses, sharedCal, measurements, point);
|
|
|
|
#ifdef GTSAM_THROW_CHEIRALITY_EXCEPTION
|
|
// verify that the triangulated point lies infront of all cameras
|
|
BOOST_FOREACH(const Pose3& pose, poses) {
|
|
const Point3& p_local = pose.transform_to(point);
|
|
if (p_local.z() <= 0)
|
|
throw(TriangulationCheiralityException());
|
|
}
|
|
#endif
|
|
|
|
return point;
|
|
}
|
|
|
|
/**
|
|
* Function to triangulate 3D landmark point from an arbitrary number
|
|
* of poses (at least 2) using the DLT. This function is similar to the one
|
|
* above, except that each camera has its own calibration. The function
|
|
* checks that the resulting point lies in front of all cameras, but has
|
|
* no other checks to verify the quality of the triangulation.
|
|
* @param cameras pinhole cameras
|
|
* @param measurements A vector of camera measurements
|
|
* @param rank tolerance, default 1e-9
|
|
* @param optimize Flag to turn on nonlinear refinement of triangulation
|
|
* @return Returns a Point3
|
|
*/
|
|
template<class CALIBRATION>
|
|
Point3 triangulatePoint3(
|
|
const std::vector<PinholeCamera<CALIBRATION> >& cameras,
|
|
const std::vector<Point2>& measurements, double rank_tol = 1e-9,
|
|
bool optimize = false) {
|
|
|
|
size_t m = cameras.size();
|
|
assert(measurements.size()==m);
|
|
|
|
if (m < 2)
|
|
throw(TriangulationUnderconstrainedException());
|
|
|
|
// construct projection matrices from poses & calibration
|
|
typedef PinholeCamera<CALIBRATION> Camera;
|
|
std::vector<Matrix> projection_matrices;
|
|
BOOST_FOREACH(const Camera& camera, cameras)
|
|
projection_matrices.push_back(
|
|
camera.calibration().K()
|
|
* sub(camera.pose().inverse().matrix(), 0, 3, 0, 4));
|
|
|
|
Point3 point = triangulateDLT(projection_matrices, measurements, rank_tol);
|
|
|
|
// The n refine using non-linear optimization
|
|
if (optimize)
|
|
point = triangulateNonlinear(cameras, measurements, point);
|
|
|
|
#ifdef GTSAM_THROW_CHEIRALITY_EXCEPTION
|
|
// verify that the triangulated point lies infront of all cameras
|
|
BOOST_FOREACH(const Camera& camera, cameras) {
|
|
const Point3& p_local = camera.pose().transform_to(point);
|
|
if (p_local.z() <= 0)
|
|
throw(TriangulationCheiralityException());
|
|
}
|
|
#endif
|
|
|
|
return point;
|
|
}
|
|
|
|
} // \namespace gtsam
|
|
|