25 #ifndef CREATE_WCS_WITH_SIP
26 #define CREATE_WCS_WITH_SIP
75 template <
class MatchT>
140 Eigen::MatrixXd
const getSipA() {
return _sipA; }
142 Eigen::MatrixXd
const getSipB() {
return _sipB; }
144 Eigen::MatrixXd
const getSipAp() {
return _sipAp; }
146 Eigen::MatrixXd
const getSipBp() {
return _sipBp; }
155 int const _sipOrder, _reverseSipOrder;
157 Eigen::MatrixXd _sipA, _sipB;
158 Eigen::MatrixXd _sipAp, _sipBp;
162 void _calculateForwardMatrices();
163 void _calculateReverseMatrices();
167 template <
class MatchT>
Measure the distortions in an image plane and express them a SIP polynomials.
Eigen::MatrixXd const getSipA()
Eigen::MatrixXd const getSipBp()
std::shared_ptr< afw::geom::SkyWcs > getNewWcs()
double getLinearScatterInPixels() const
Compute the median radial separation between items in this object's match list.
std::shared_ptr< CreateWcsWithSip > Ptr
std::shared_ptr< CreateWcsWithSip const > ConstPtr
int getNGrid() const
Return the number of grid points (on each axis) used in inverse SIP transform.
Eigen::MatrixXd const getSipB()
int getNPoints() const
Return the number of points in the catalogue.
int getOrder() const
Return the number of terms in the SIP matrix.
double getScatterInPixels() const
Compute the median separation, in pixels, between items in this object's match list.
Eigen::MatrixXd const getSipAp()
geom::Angle getScatterOnSky() const
Compute the median on-sky separation between items in this object's match list.
geom::Angle getLinearScatterOnSky() const
Compute the median on-sky separation between items in this object's match list,.
CreateWcsWithSip(std::vector< MatchT > const &matches, afw::geom::SkyWcs const &linearWcs, int const order, geom::Box2I const &bbox=geom::Box2I(), int const ngrid=0)
Construct a CreateWcsWithSip.
CreateWcsWithSip< MatchT > makeCreateWcsWithSip(std::vector< MatchT > const &matches, afw::geom::SkyWcs const &linearWcs, int const order, geom::Box2I const &bbox=geom::Box2I(), int const ngrid=0)
Factory function for CreateWcsWithSip.