25 #include "wcslib/wcs.h"
27 namespace lsst {
namespace meas {
namespace algorithms {
30 PTR(afw::image::Wcs) parent,
35 afw::image::Wcs(*parent),
36 _parent(parent), _xBin(xBin), _yBin(yBin), _xy0(xy0),
37 _binnedToOriginal(afw::geom::AffineTransform::makeTranslation(afw::geom::Extent2D(_xy0))*
38 afw::geom::AffineTransform::makeScaling(_xBin, _yBin)),
39 _originalToBinned(_binnedToOriginal.invert())
42 afw::geom::Point2D
const crpix = _originalToBinned(parent->getPixelOrigin());
43 _wcsInfo->crpix[0] = crpix.getX() + 1;
44 _wcsInfo->crpix[1] = crpix.getY() + 1;
50 PTR(afw::coord::Coord) coord = _parent->pixelToSky(_binnedToOriginal(afw::geom::Point2D(x, y)));
51 skyTmp[0] = coord->getLongitude();
52 skyTmp[1] = coord->getLatitude();
57 return _originalToBinned(_parent->skyToPixel(*makeCorrectCoord(sky1, sky2)));
virtual void pixelToSkyImpl(double pixel1, double pixel2, afw::geom::Angle skyTmp[2]) const
BinnedWcs(boost::shared_ptr< afw::image::Wcs > parent, unsigned int xBin, unsigned int yBin, afw::geom::Point2I xy0)
virtual afw::geom::Point2D skyToPixelImpl(afw::geom::Angle sky1, afw::geom::Angle sky2) const