LSSTApplications  10.0-2-g4f67435,11.0.rc2+1,11.0.rc2+12,11.0.rc2+3,11.0.rc2+4,11.0.rc2+5,11.0.rc2+6,11.0.rc2+7,11.0.rc2+8
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 
32  TanWcs const &tanWcs,
33  geom::XYTransform const &pixelsToTanPixels
34 ) :
35  TanWcs(tanWcs),
36  _pixelsToTanPixelsPtr(pixelsToTanPixels.clone())
37 {
38  if (tanWcs.hasDistortion()) {
39  throw LSST_EXCEPT(pex::exceptions::InvalidParameterError, "tanWcs has distortion terms");
40  }
41 }
42 
43 PTR(Wcs) DistortedTanWcs::clone() const {
44  return PTR(Wcs)(new DistortedTanWcs(*this));
45 }
46 
47 bool DistortedTanWcs::operator==(Wcs const & rhs) const {
48  throw LSST_EXCEPT(pex::exceptions::LogicError, "== is not implemented");
49 }
50 
51 void DistortedTanWcs::flipImage(int flipLR, int flipTB, geom::Extent2I dimensions) const {
52  throw LSST_EXCEPT(pex::exceptions::LogicError, "flipImage is not implemented");
53 }
54 
56  throw LSST_EXCEPT(pex::exceptions::LogicError, "rotateImageBy90 is not implemented");
57 }
58 
59 void DistortedTanWcs::shiftReferencePixel(double dx, double dy) {
60  throw LSST_EXCEPT(pex::exceptions::LogicError, "shiftReferencePixel is not implemented");
61 }
62 
64  auto const tanPos = TanWcs::skyToPixelImpl(sky1, sky2);
65  return _pixelsToTanPixelsPtr->reverseTransform(tanPos);
66 }
67 
68 void DistortedTanWcs::pixelToSkyImpl(double pixel1, double pixel2, geom::Angle sky[2]) const {
69  auto const pos = geom::Point2D(pixel1, pixel2);
70  auto const tanPos = _pixelsToTanPixelsPtr->forwardTransform(pos);
71  TanWcs::pixelToSkyImpl(tanPos[0], tanPos[1], sky);
72 }
73 
74 }}} // namespace lsst::afw::image
virtual void flipImage(int flipLR, int flipTB, lsst::afw::geom::Extent2I dimensions) const
An include file to include the header files for lsst::afw::geom.
#define PTR(...)
Definition: base.h:41
virtual geom::Point2D skyToPixelImpl(geom::Angle sky1, geom::Angle sky2) const
bool operator==(Wcs const &other) const
Implementation of the WCS standard for a any projection.
Definition: Wcs.h:107
virtual geom::Point2D skyToPixelImpl(geom::Angle sky1, geom::Angle sky2) const
table::Key< table::Array< Kernel::Pixel > > image
Definition: FixedKernel.cc:117
DistortedTanWcs(TanWcs const &tanWcs, geom::XYTransform const &pixelsToTanPixels)
Construct a DistortedTanWcs.
afw::table::PointKey< int > dimensions
Implementation of the WCS standard for the special case of the Gnomonic (tangent plane) projection...
Definition: TanWcs.h:68
bool hasDistortion() const
Definition: TanWcs.h:140
virtual void rotateImageBy90(int nQuarter, lsst::afw::geom::Extent2I dimensions) const
virtual void shiftReferencePixel(double dx, double dy)
#define LSST_EXCEPT(type,...)
Definition: Exception.h:46
virtual void pixelToSkyImpl(double pixel1, double pixel2, geom::Angle skyTmp[2]) const
Combination of a TAN WCS and a distortion model.
Virtual base class for 2D transforms.
Definition: XYTransform.h:48
boost::shared_ptr< geom::XYTransform > _pixelsToTanPixelsPtr
Point< double, 2 > Point2D
Definition: Point.h:286
Include files required for standard LSST Exception handling.
virtual void pixelToSkyImpl(double pixel1, double pixel2, geom::Angle skyTmp[2]) const