LSST Applications  21.0.0-147-g0e635eb1+1acddb5be5,22.0.0+052faf71bd,22.0.0+1ea9a8b2b2,22.0.0+6312710a6c,22.0.0+729191ecac,22.0.0+7589c3a021,22.0.0+9f079a9461,22.0.1-1-g7d6de66+b8044ec9de,22.0.1-1-g87000a6+536b1ee016,22.0.1-1-g8e32f31+6312710a6c,22.0.1-10-gd060f87+016f7cdc03,22.0.1-12-g9c3108e+df145f6f68,22.0.1-16-g314fa6d+c825727ab8,22.0.1-19-g93a5c75+d23f2fb6d8,22.0.1-19-gb93eaa13+aab3ef7709,22.0.1-2-g8ef0a89+b8044ec9de,22.0.1-2-g92698f7+9f079a9461,22.0.1-2-ga9b0f51+052faf71bd,22.0.1-2-gac51dbf+052faf71bd,22.0.1-2-gb66926d+6312710a6c,22.0.1-2-gcb770ba+09e3807989,22.0.1-20-g32debb5+b8044ec9de,22.0.1-23-gc2439a9a+fb0756638e,22.0.1-3-g496fd5d+09117f784f,22.0.1-3-g59f966b+1e6ba2c031,22.0.1-3-g849a1b8+f8b568069f,22.0.1-3-gaaec9c0+c5c846a8b1,22.0.1-32-g5ddfab5d3+60ce4897b0,22.0.1-4-g037fbe1+64e601228d,22.0.1-4-g8623105+b8044ec9de,22.0.1-5-g096abc9+d18c45d440,22.0.1-5-g15c806e+57f5c03693,22.0.1-7-gba73697+57f5c03693,master-g6e05de7fdc+c1283a92b8,master-g72cdda8301+729191ecac,w.2021.39
LSST Data Management Base Package
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"
29 #include "lsst/geom/SpherePoint.h"
32 #include "lsst/afw/geom/SkyWcs.h"
34 
35 namespace lsst {
36 namespace meas {
37 namespace astrom {
38 
47 public:
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 
63 protected:
76  SipTransformBase(geom::Point2D const& pixelOrigin, geom::LinearTransform const& cdMatrix,
78  : _pixelOrigin(pixelOrigin), _cdMatrix(cdMatrix), _poly(poly) {}
79 
80  SipTransformBase(SipTransformBase const& other) = default;
81  SipTransformBase(SipTransformBase&& other) = default;
82  SipTransformBase& operator=(SipTransformBase const& other) = default;
84 
85  void swap(SipTransformBase& other) {
88  _poly.swap(other._poly);
89  }
90 
92 
96 };
97 
137 public:
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 
179  SipForwardTransform(SipForwardTransform const& other) = default;
180 
182 
184 
186 
188 
193 
197  geom::Point2D operator()(geom::Point2D const& uv) const;
198 
204 };
205 
247 public:
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 
289  SipReverseTransform(SipReverseTransform const& other) = default;
290 
292 
294 
296 
297  void swap(SipReverseTransform& other) {
298  SipTransformBase::swap(other);
299  std::swap(_cdInverse, other._cdInverse);
300  }
301 
306 
310  geom::Point2D operator()(geom::Point2D const& xy) const;
311 
317 
318 private:
319  friend class PolynomialTransform;
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
afw::table::PointKey< int > dimensions
Definition: GaussianPsf.cc:48
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:66
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...
Definition: SipTransform.h:136
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...
Definition: SipTransform.cc:95
SipForwardTransform & operator=(SipForwardTransform const &other)=default
geom::AffineTransform linearize(geom::Point2D const &in) const
Return an approximate affine transform at the given point.
Definition: SipTransform.cc:85
void swap(SipForwardTransform &other)
Definition: SipTransform.h:187
geom::Point2D operator()(geom::Point2D const &uv) const
Apply the transform to a point.
Definition: SipTransform.cc:90
static SipForwardTransform convert(PolynomialTransform const &poly, geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix)
Convert a PolynomialTransform to an equivalent SipForwardTransform.
Definition: SipTransform.cc:52
SipForwardTransform(geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix, PolynomialTransform const &forwardSipPoly)
Construct a SipForwardTransform from its components.
Definition: SipTransform.h:175
A transform that maps intermediate world coordinates to pixel coordinates according to the SIP conven...
Definition: SipTransform.h:246
SipReverseTransform(geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix, PolynomialTransform const &reverseSipPoly)
Construct a SipReverseTransform from its components.
Definition: SipTransform.h:285
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)
Definition: SipTransform.h:297
SipReverseTransform(SipReverseTransform &&other)=default
SipReverseTransform & operator=(SipReverseTransform &&other)=default
static SipReverseTransform convert(PolynomialTransform const &poly, geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix)
Convert a PolynomialTransform to an equivalent SipReverseTransform.
SipReverseTransform & operator=(SipReverseTransform const &other)=default
geom::Point2D operator()(geom::Point2D const &xy) const
Apply the transform to a point.
Base class for SIP transform objects.
Definition: SipTransform.h:46
SipTransformBase(SipTransformBase &&other)=default
void transformPixelsInPlace(geom::AffineTransform const &s)
Definition: SipTransform.cc:40
geom::LinearTransform const & getCdMatrix() const
Return the CD matrix of the transform.
Definition: SipTransform.h:56
void swap(SipTransformBase &other)
Definition: SipTransform.h:85
SipTransformBase(SipTransformBase const &other)=default
geom::Point2D const & getPixelOrigin() const
Return the pixel origin (CRPIX, but zero-indexed) of the transform.
Definition: SipTransform.h:51
SipTransformBase & operator=(SipTransformBase const &other)=default
geom::LinearTransform _cdMatrix
Definition: SipTransform.h:94
SipTransformBase & operator=(SipTransformBase &&other)=default
SipTransformBase(geom::Point2D const &pixelOrigin, geom::LinearTransform const &cdMatrix, PolynomialTransform const &poly)
Construct a SipTransformBase from its components.
Definition: SipTransform.h:76
PolynomialTransform const & getPoly() const
Return the polynomial component of the transform (A,B) or (AP,BP).
Definition: SipTransform.h:61
Low-level polynomials (including special polynomials) in C++.
Definition: Basis1d.h:26
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.
A base class for image defects.
T swap(T... args)