Loading [MathJax]/extensions/tex2jax.js
LSST Applications g0fba68d861+05816baf74,g1ec0fe41b4+f536777771,g1fd858c14a+a9301854fb,g35bb328faa+fcb1d3bbc8,g4af146b050+a5c07d5b1d,g4d2262a081+6e5fcc2a4e,g53246c7159+fcb1d3bbc8,g56a49b3a55+9c12191793,g5a012ec0e7+3632fc3ff3,g60b5630c4e+ded28b650d,g67b6fd64d1+ed4b5058f4,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g8352419a5c+fcb1d3bbc8,g87b7deb4dc+7b42cf88bf,g8852436030+e5453db6e6,g89139ef638+ed4b5058f4,g8e3bb8577d+d38d73bdbd,g9125e01d80+fcb1d3bbc8,g94187f82dc+ded28b650d,g989de1cb63+ed4b5058f4,g9d31334357+ded28b650d,g9f33ca652e+50a8019d8c,gabe3b4be73+1e0a283bba,gabf8522325+fa80ff7197,gb1101e3267+d9fb1f8026,gb58c049af0+f03b321e39,gb665e3612d+2a0c9e9e84,gb89ab40317+ed4b5058f4,gcf25f946ba+e5453db6e6,gd6cbbdb0b4+bb83cc51f8,gdd1046aedd+ded28b650d,gde0f65d7ad+941d412827,ge278dab8ac+d65b3c2b70,ge410e46f29+ed4b5058f4,gf23fb2af72+b7cae620c0,gf5e32f922b+fcb1d3bbc8,gf67bdafdda+ed4b5058f4,w.2025.16
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
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,
77 PolynomialTransform const& poly)
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).
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
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.
Point< double, 2 > Point2D
Definition Point.h:324
Extent< int, 2 > Extent2I
Definition Extent.h:397
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)