LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
|
Measure the distortions in an image plane and express them a SIP polynomials. More...
#include <CreateWcsWithSip.h>
Public Types | |
typedef std::shared_ptr< CreateWcsWithSip > | Ptr |
typedef std::shared_ptr< CreateWcsWithSip const > | ConstPtr |
Public Member Functions | |
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. | |
std::shared_ptr< afw::geom::SkyWcs > | getNewWcs () |
double | getScatterInPixels () const |
Compute the median separation, in pixels, between items in this object's match list. | |
geom::Angle | getScatterOnSky () const |
Compute the median on-sky separation between items in this object's match list. | |
double | getLinearScatterInPixels () const |
Compute the median radial 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,. | |
int | getOrder () const |
Return the number of terms in the SIP matrix. | |
int | getNPoints () const |
Return the number of points in the catalogue. | |
int | getNGrid () const |
Return the number of grid points (on each axis) used in inverse SIP transform. | |
Eigen::MatrixXd const | getSipA () |
Eigen::MatrixXd const | getSipB () |
Eigen::MatrixXd const | getSipAp () |
Eigen::MatrixXd const | getSipBp () |
Measure the distortions in an image plane and express them a SIP polynomials.
Given a list of matching sources between a catalogue and an image, and a linear Wcs that describes the mapping from pixel space in the image and ra/dec space in the catalogue, calculate discrepancies between the two and compute SIP distortion polynomials to describe the discrepancy
SIP polynomials are defined in Shupe at al. (2005) ASPC 347 491.
Note that the SIP standard insists (although it is only mentioned obliquly between Eqns 3 and 4) that the lowest three terms in the distortion polynomials be zero (A00, A10, A01, B00, etc.). To achieve this, we need to adjust the values of CD and CRPIX from the input wcs. This may not be the behaviour you expect.
A Wcs may be created in a variety of ways (e.g. lsst::meas::astrom::net::GlobalAstrometrySolution ), and the list of matched sources (matches) can be generated with the matchRaDec function.
Note that the matches must be one-to-one; this is ensured by passing closest=true to matchRaDec.
Definition at line 76 of file CreateWcsWithSip.h.
typedef std::shared_ptr<CreateWcsWithSip const> lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::ConstPtr |
Definition at line 79 of file CreateWcsWithSip.h.
typedef std::shared_ptr<CreateWcsWithSip> lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::Ptr |
Definition at line 78 of file CreateWcsWithSip.h.
lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::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.
Constructor.
[in] | matches | list of matches |
[in] | linearWcs | initial WCS, typically pure TAN but need not be |
[in] | order | SIP order for fit WCS |
[in] | bbox | bounding box over which to compute the reverse SIP transform. If empty then a bounding box is computed based on the matches, extended a bit to allow for the fact that the sources will not necessarily reach to each edge of the image. Specifially the box is grown by dimensions/sqrt(number of matches). |
[in] | ngrid | number of points along x or y for the grid of points on which the reverse SIP transform is computed |
Definition at line 102 of file CreateWcsWithSip.cc.
double lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::getLinearScatterInPixels | ( | ) | const |
Compute the median radial separation between items in this object's match list.
For each match, project the reference object coord to pixels using the initial "linearWcs" WCS, and measure the radial separation to the source centroid.
Definition at line 358 of file CreateWcsWithSip.cc.
geom::Angle lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::getLinearScatterOnSky | ( | ) | const |
Compute the median on-sky separation between items in this object's match list,.
For each match, project the source centroid to RA,Dec using the initial "linearWcs" WCS, and measure the on-sky angular separation to the reference source coord.
Definition at line 370 of file CreateWcsWithSip.cc.
|
inline |
Definition at line 98 of file CreateWcsWithSip.h.
|
inline |
Return the number of grid points (on each axis) used in inverse SIP transform.
Definition at line 137 of file CreateWcsWithSip.h.
|
inline |
|
inline |
Return the number of terms in the SIP matrix.
Definition at line 133 of file CreateWcsWithSip.h.
double lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::getScatterInPixels | ( | ) | const |
Compute the median separation, in pixels, between items in this object's match list.
For each match, project the reference object coord to pixels using the fit TAN-SIP WCS, and measure the radial separation to the source centroid
Definition at line 352 of file CreateWcsWithSip.cc.
geom::Angle lsst::meas::astrom::sip::CreateWcsWithSip< MatchT >::getScatterOnSky | ( | ) | const |
Compute the median on-sky separation between items in this object's match list.
For each match, project the source centroid to RA,Dec using the fit TAN-SIP WCS, and measure the on-sky angular separation to the reference source coord.
Definition at line 364 of file CreateWcsWithSip.cc.
|
inline |
Definition at line 140 of file CreateWcsWithSip.h.
|
inline |
Definition at line 144 of file CreateWcsWithSip.h.
|
inline |
Definition at line 142 of file CreateWcsWithSip.h.
|
inline |
Definition at line 146 of file CreateWcsWithSip.h.