lsst.meas.astrom gcf9c84bd28+876f2a9e86
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,
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-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
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
SipForwardTransform & operator=(SipForwardTransform const &other)=default
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 & operator=(SipReverseTransform const &other)=default
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
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.
Definition: SipTransform.h:46
geom::Point2D const & getPixelOrigin() const
Return the pixel origin (CRPIX, but zero-indexed) of the transform.
Definition: SipTransform.h:51
PolynomialTransform const & getPoly() const
Return the polynomial component of the transform (A,B) or (AP,BP).
Definition: SipTransform.h:61
SipTransformBase(SipTransformBase &&other)=default
SipTransformBase & operator=(SipTransformBase const &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::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
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)