LSSTApplications  11.0-22-g33de520,13.0+153,14.0+52,14.0+57,14.0-1-g013352c+36,14.0-1-g13ef843+9,14.0-1-g4b114ac+14,14.0-1-g7257b6a+12,14.0-1-g8b7e855+51,14.0-13-g7a60b79+2,14.0-14-g87d16e8+10,14.0-14-gbf7a6f8a,14.0-17-g4f4ea82+5,14.0-2-g319577b+11,14.0-2-ga5af9b6+10,14.0-22-gc48c03f+3,14.0-3-g20413be+3,14.0-46-g76222d5f+3,14.0-47-g0a51fac97,14.0-5-g744ff5f+2,14.0-5-g86eb1bd+31,14.0-6-gd5b81a9+6,14.0-6-ge2c9487+42,14.0-8-g7f6dd6b+6,14.0-8-gb81b6e9+4,14.0-9-g11010eb,14.0-9-g330837b+5
LSSTDataManagementBasePackage
DistortedTanWcs.cc
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008, 2009, 2010 LSST Corporation.
5  *
6  * This product includes software developed by the
7  * LSST Project (http://www.lsst.org/).
8  *
9  * This program is free software: you can redistribute it and/or modify
10  * it under the terms of the GNU General Public License as published by
11  * the Free Software Foundation, either version 3 of the License, or
12  * (at your option) any later version.
13  *
14  * This program is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17  * GNU General Public License for more details.
18  *
19  * You should have received a copy of the LSST License Statement and
20  * the GNU General Public License along with this program. If not,
21  * see <http://www.lsstcorp.org/LegalNotices/>.
22  */
23 #include "lsst/pex/exceptions.h"
24 #include "lsst/afw/geom.h"
26 
27 namespace lsst {
28 namespace afw {
29 namespace image {
30 
31 DistortedTanWcs::DistortedTanWcs(TanWcs const &tanWcs, Transform const &pixelsToTanPixels)
32  : TanWcs(tanWcs), _pixelsToTanPixelsPtr(std::make_shared<Transform>(pixelsToTanPixels)) {
33  if (tanWcs.hasDistortion()) {
34  throw LSST_EXCEPT(pex::exceptions::InvalidParameterError, "tanWcs has distortion terms");
35  }
36 }
37 
40 
42  return std::shared_ptr<Wcs>(new DistortedTanWcs(*this));
43 }
44 
45 bool DistortedTanWcs::operator==(Wcs const &rhs) const {
46  throw LSST_EXCEPT(pex::exceptions::LogicError, "== is not implemented");
47 }
48 
49 void DistortedTanWcs::flipImage(int flipLR, int flipTB, geom::Extent2I dimensions) const {
50  throw LSST_EXCEPT(pex::exceptions::LogicError, "flipImage is not implemented");
51 }
52 
54  throw LSST_EXCEPT(pex::exceptions::LogicError, "rotateImageBy90 is not implemented");
55 }
56 
57 void DistortedTanWcs::shiftReferencePixel(double dx, double dy) {
58  throw LSST_EXCEPT(pex::exceptions::LogicError, "shiftReferencePixel is not implemented");
59 }
60 
62  auto const tanPos = TanWcs::skyToPixelImpl(sky1, sky2);
63  return _pixelsToTanPixelsPtr->applyInverse(tanPos);
64 }
65 
66 void DistortedTanWcs::pixelToSkyImpl(double pixel1, double pixel2, geom::Angle sky[2]) const {
67  auto const pos = geom::Point2D(pixel1, pixel2);
68  auto const tanPos = _pixelsToTanPixelsPtr->applyForward(pos);
69  TanWcs::pixelToSkyImpl(tanPos[0], tanPos[1], sky);
70 }
71 }
72 }
73 } // namespace lsst::afw::image
virtual void pixelToSkyImpl(double pixel1, double pixel2, geom::Angle skyTmp[2]) const
Worker routine for skyToPixel.
virtual std::shared_ptr< Wcs > clone() const
Polymorphic deep-copy.
virtual void shiftReferencePixel(double dx, double dy)
void pixelToSkyImpl(double pixel1, double pixel2, geom::Angle skyTmp[2]) const override
Worker routine for pixelToSky.
STL namespace.
Implementation of the WCS standard for a any projection.
Definition: Wcs.h:104
geom::Point2D skyToPixelImpl(geom::Angle sky1, geom::Angle sky2) const override
Worker routine for skyToPixel.
bool operator==(Wcs const &other) const
virtual geom::Point2D skyToPixelImpl(geom::Angle sky1, geom::Angle sky2) const
Worker routine for pixelToSky.
bool hasDistortion() const override
Definition: TanWcs.h:125
afw::table::PointKey< int > dimensions
Definition: GaussianPsf.cc:44
virtual void flipImage(int flipLR, int flipTB, lsst::afw::geom::Extent2I dimensions) const
A class representing an angle.
Definition: Angle.h:102
A base class for image defects.
Definition: cameraGeom.dox:3
Implementation of the WCS standard for the special case of the Gnomonic (tangent plane) projection...
Definition: TanWcs.h:68
Point< double, 2 > Point2D
Definition: Point.h:300
Include files required for standard LSST Exception handling.
virtual void rotateImageBy90(int nQuarter, lsst::afw::geom::Extent2I dimensions) const
#define LSST_EXCEPT(type,...)
Create an exception with a given type and message and optionally other arguments (dependent on the ty...
Definition: Exception.h:46
Combination of a TAN WCS and a distortion model.
std::shared_ptr< Transform > _pixelsToTanPixelsPtr
Transform LSST spatial data, such as Point2D and IcrsCoord, using an AST transform.
Definition: Transform.h:69
DistortedTanWcs(TanWcs const &tanWcs, Transform const &pixelsToTanPixels)
Construct a DistortedTanWcs.