74 lines
1.9 KiB
C++
74 lines
1.9 KiB
C++
/**
|
|
* @file expressions.h
|
|
* @brief Common expressions for solving geometry/slam/sfm problems
|
|
* @date Oct 1, 2014
|
|
* @author Frank Dellaert
|
|
*/
|
|
|
|
#pragma once
|
|
|
|
#include <gtsam_unstable/nonlinear/expressions.h>
|
|
#include <gtsam/geometry/Cal3_S2.h>
|
|
#include <gtsam/geometry/PinholeCamera.h>
|
|
|
|
namespace gtsam {
|
|
|
|
// 2D Geometry
|
|
|
|
typedef Expression<Point2> Point2_;
|
|
typedef Expression<Rot2> Rot2_;
|
|
typedef Expression<Pose2> Pose2_;
|
|
|
|
Point2_ transform_to(const Pose2_& x, const Point2_& p) {
|
|
Point2(Pose2::*transform)(const Point2& p,
|
|
boost::optional<Matrix23&> H1,
|
|
boost::optional<Matrix2&> H2) const = &Pose2::transform_to;
|
|
|
|
return Point2_(x, transform, p);
|
|
}
|
|
|
|
// 3D Geometry
|
|
|
|
typedef Expression<Point3> Point3_;
|
|
typedef Expression<Rot3> Rot3_;
|
|
typedef Expression<Pose3> Pose3_;
|
|
|
|
Point3_ transform_to(const Pose3_& x, const Point3_& p) {
|
|
|
|
Point3(Pose3::*transform)(const Point3& p,
|
|
boost::optional<Matrix36&> Dpose,
|
|
boost::optional<Matrix3&> Dpoint) const = &Pose3::transform_to;
|
|
|
|
return Point3_(x, transform, p);
|
|
}
|
|
|
|
// Projection
|
|
|
|
typedef Expression<Cal3_S2> Cal3_S2_;
|
|
|
|
Point2_ project(const Point3_& p_cam) {
|
|
return Point2_(PinholeCamera<Cal3_S2>::project_to_camera, p_cam);
|
|
}
|
|
|
|
Point2 project6(const Pose3& x, const Point3& p, const Cal3_S2& K,
|
|
boost::optional<Matrix26&> Dpose, boost::optional<Matrix23&> Dpoint,
|
|
boost::optional<Matrix25&> Dcal) {
|
|
return PinholeCamera<Cal3_S2>(x, K).project(p, Dpose, Dpoint, Dcal);
|
|
}
|
|
|
|
Point2_ project3(const Pose3_& x, const Point3_& p, const Cal3_S2_& K) {
|
|
return Point2_(project6, x, p, K);
|
|
}
|
|
|
|
template<class CAL>
|
|
Point2_ uncalibrate(const Expression<CAL>& K, const Point2_& xy_hat) {
|
|
Point2(CAL::*uncal)(const Point2& p,
|
|
boost::optional<Matrix25&> Dpose,
|
|
boost::optional<Matrix2&> Dpoint) const = &CAL::uncalibrate;
|
|
|
|
return Point2_(K, uncal, xy_hat);
|
|
}
|
|
|
|
} // \namespace gtsam
|
|
|