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>
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 getSipBp()
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.
int getNPoints() const
Return the number of points in the catalogue.
Eigen::MatrixXd const getSipB()
std::shared_ptr< CreateWcsWithSip > Ptr
Eigen::MatrixXd const getSipAp()
std::shared_ptr< afw::geom::SkyWcs > getNewWcs()
std::shared_ptr< CreateWcsWithSip const > ConstPtr
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.
Measure the distortions in an image plane and express them a SIP polynomials.
Eigen::MatrixXd const getSipA()
geom::Angle getScatterOnSky() const
Compute the median on-sky separation between items in this object's match list.
int getNGrid() const
Return the number of grid points (on each axis) used in inverse SIP transform.
double getLinearScatterInPixels() const
Compute the median radial separation between items in this object's match list.