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
Loading...
Searching...
No Matches
SipTransform.h
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2016 LSST/AURA
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24#ifndef LSST_MEAS_ASTROM_SipTransform_INCLUDED
25#define LSST_MEAS_ASTROM_SipTransform_INCLUDED
26
27#include "lsst/geom/Extent.h"
28#include "lsst/geom/Point.h"
34
35namespace lsst {
36namespace meas {
37namespace astrom {
38
47public:
51 geom::Point2D const& getPixelOrigin() const { return _pixelOrigin; }
52
56 geom::LinearTransform const& getCdMatrix() const { return _cdMatrix; }
57
61 PolynomialTransform const& getPoly() const { return _poly; }
62
63protected:
76 SipTransformBase(geom::Point2D const& pixelOrigin, geom::LinearTransform const& cdMatrix,
78 : _pixelOrigin(pixelOrigin), _cdMatrix(cdMatrix), _poly(poly) {}
79
80 SipTransformBase(SipTransformBase const& other) = default;
84
85 void swap(SipTransformBase& other) {
88 _poly.swap(other._poly);
89 }
90
92
96};
97
137public:
145 static SipForwardTransform convert(PolynomialTransform const& poly, geom::Point2D const& pixelOrigin,
146 geom::LinearTransform const& cdMatrix);
147
156 geom::Point2D const& pixelOrigin,
157 geom::LinearTransform const& cdMatrix);
158
167
175 SipForwardTransform(geom::Point2D const& pixelOrigin, geom::LinearTransform const& cdMatrix,
176 PolynomialTransform const& forwardSipPoly)
177 : SipTransformBase(pixelOrigin, cdMatrix, forwardSipPoly) {}
178
180
182
184
186
188
193
197 geom::Point2D operator()(geom::Point2D const& uv) const;
198
204};
205
247public:
255 static SipReverseTransform convert(PolynomialTransform const& poly, geom::Point2D const& pixelOrigin,
256 geom::LinearTransform const& cdMatrix);
257
266 geom::Point2D const& pixelOrigin,
267 geom::LinearTransform const& cdMatrix);
268
277
285 SipReverseTransform(geom::Point2D const& pixelOrigin, geom::LinearTransform const& cdMatrix,
286 PolynomialTransform const& reverseSipPoly)
287 : SipTransformBase(pixelOrigin, cdMatrix, reverseSipPoly), _cdInverse(cdMatrix.inverted()) {}
288
290
292
294
296
299 std::swap(_cdInverse, other._cdInverse);
300 }
301
306
310 geom::Point2D operator()(geom::Point2D const& xy) const;
311
317
318private:
321 geom::LinearTransform _cdInverse;
322};
323
339 SipReverseTransform const& sipReverse,
340 geom::SpherePoint const& skyOrigin);
341
359 geom::AffineTransform const& s);
360
370 geom::Extent2I const& dimensions);
371
372} // namespace astrom
373} // namespace meas
374} // namespace lsst
375
376#endif // !LSST_MEAS_ASTROM_SipTransform_INCLUDED
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition SkyWcs.h:117
An affine coordinate transformation consisting of a linear transformation and an offset.
A 2D linear coordinate transformation.
Point in an unspecified spherical coordinate system.
Definition SpherePoint.h:57
A 2-d coordinate transform represented by a pair of standard polynomials (one for each coordinate).
void swap(PolynomialTransform &other)
Lightweight swap.
A 2-d coordinate transform represented by a lazy composition of an AffineTransform,...
A transform that maps pixel coordinates to intermediate world coordinates according to the SIP conven...
SipForwardTransform(SipForwardTransform &&other)=default
SipForwardTransform & operator=(SipForwardTransform &&other)=default
SipForwardTransform(SipForwardTransform const &other)=default
SipForwardTransform transformPixels(geom::AffineTransform const &s) const
Return a new forward SIP transform that includes a transformation of the pixel coordinate system by t...
geom::AffineTransform linearize(geom::Point2D const &in) const
Return an approximate affine transform at the given point.
void swap(SipForwardTransform &other)
SipForwardTransform & operator=(SipForwardTransform const &other)=default
geom::Point2D operator()(geom::Point2D const &uv) const
Apply the transform to a point.
static SipForwardTransform convert(PolynomialTransform const &poly, geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix)
Convert a PolynomialTransform to an equivalent SipForwardTransform.
SipForwardTransform(geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix, PolynomialTransform const &forwardSipPoly)
Construct a SipForwardTransform from its components.
A transform that maps intermediate world coordinates to pixel coordinates according to the SIP conven...
SipReverseTransform & operator=(SipReverseTransform const &other)=default
SipReverseTransform(geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix, PolynomialTransform const &reverseSipPoly)
Construct a SipReverseTransform from its components.
SipReverseTransform transformPixels(geom::AffineTransform const &s) const
Return a new reverse SIP transform that includes a transformation of the pixel coordinate system by t...
geom::AffineTransform linearize(geom::Point2D const &in) const
Return an approximate affine transform at the given point.
SipReverseTransform(SipReverseTransform const &other)=default
void swap(SipReverseTransform &other)
SipReverseTransform(SipReverseTransform &&other)=default
static SipReverseTransform convert(PolynomialTransform const &poly, geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix)
Convert a PolynomialTransform to an equivalent SipReverseTransform.
geom::Point2D operator()(geom::Point2D const &xy) const
Apply the transform to a point.
SipReverseTransform & operator=(SipReverseTransform &&other)=default
Base class for SIP transform objects.
geom::Point2D const & getPixelOrigin() const
Return the pixel origin (CRPIX, but zero-indexed) of the transform.
PolynomialTransform const & getPoly() const
Return the polynomial component of the transform (A,B) or (AP,BP).
SipTransformBase(SipTransformBase &&other)=default
SipTransformBase & operator=(SipTransformBase const &other)=default
void transformPixelsInPlace(geom::AffineTransform const &s)
geom::LinearTransform const & getCdMatrix() const
Return the CD matrix of the transform.
void swap(SipTransformBase &other)
SipTransformBase(SipTransformBase const &other)=default
SipTransformBase & operator=(SipTransformBase &&other)=default
SipTransformBase(geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix, PolynomialTransform const &poly)
Construct a SipTransformBase from its components.
Low-level polynomials (including special polynomials) in C++.
std::shared_ptr< afw::geom::SkyWcs > transformWcsPixels(afw::geom::SkyWcs const &wcs, geom::AffineTransform const &s)
Create a new SkyWcs whose pixel coordinate system has been transformed via an affine transform.
std::shared_ptr< afw::geom::SkyWcs > rotateWcsPixelsBy90(afw::geom::SkyWcs const &wcs, int nQuarter, geom::Extent2I const &dimensions)
Return a new SkyWcs that represents a rotation of the image it corresponds to about the image's cente...
std::shared_ptr< afw::geom::SkyWcs > makeWcs(SipForwardTransform const &sipForward, SipReverseTransform const &sipReverse, geom::SpherePoint const &skyOrigin)
Create a new TAN SIP Wcs from a pair of SIP transforms and the sky origin.
T swap(T... args)