ExposurePatch.h

Go to the documentation of this file.

// -*- LSST-C++ -*-

/*
 * LSST Data Management System
 * Copyright 2008, 2009, 2010 LSST Corporation.
 *
 * This product includes software developed by the
 * LSST Project (http://www.lsst.org/).
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the LSST License Statement and
 * the GNU General Public License along with this program.  If not,
 * see <http://www.lsstcorp.org/LegalNotices/>.
 */

#if !defined(LSST_MEAS_ALGORITHMS_EXPOSURE_PATCH_H)
#define LSST_MEAS_ALGORITHMS_EXPOSURE_PATCH_H

#include "lsst/base.h"
#include "lsst/geom/AffineTransform.h"
#include "lsst/geom/Point.h"
#include "lsst/afw/detection/Footprint.h"
#include "lsst/afw/geom/SkyWcs.h"

namespace lsst {
namespace meas {
namespace algorithms {

template <typename ExposureT>
class ExposurePatch {
public:
    typedef unsigned char FlagT;  
    typedef PTR(ExposurePatch) Ptr;
    typedef CONST_PTR(ExposurePatch) ConstPtr;

    ExposurePatch(CONST_PTR(ExposureT) exp,                   
                  CONST_PTR(afw::detection::Footprint) foot,  
                  geom::Point2D const& center                 
                  )
            : _exp(exp), _foot(foot), _center(center), _fromStandard(), _toStandard() {}
    ExposurePatch(CONST_PTR(ExposureT) exp,                       
                  afw::detection::Footprint const& standardFoot,  
                  geom::Point2D const& standardCenter,            
                  afw::geom::SkyWcs const& standardWcs            
                  )
            : _exp(exp) {
        afw::geom::SkyWcs const& expWcs = *exp->getWcs();
        auto const sky = standardWcs.pixelToSky(standardCenter);
        const_cast<CONST_PTR(afw::detection::Footprint)&>(_foot) =
                standardFoot.transform(standardWcs, expWcs, exp->getBBox());
        const_cast<geom::Point2D&>(_center) = expWcs.skyToPixel(sky);
        const_cast<geom::AffineTransform&>(_fromStandard) =
                standardWcs.linearizePixelToSky(sky) * expWcs.linearizeSkyToPixel(sky);
        const_cast<geom::AffineTransform&>(_toStandard) =
                expWcs.linearizePixelToSky(sky) * standardWcs.linearizeSkyToPixel(sky);
    }

    CONST_PTR(ExposureT) const getExposure() const { return _exp; }
    CONST_PTR(afw::detection::Footprint) const getFootprint() const { return _foot; }
    geom::Point2D const& getCenter() const { return _center; }
    geom::AffineTransform const& fromStandard() const { return _fromStandard; }
    geom::AffineTransform const& toStandard() const { return _toStandard; }

    void setCenter(geom::Point2D const& center) { _center = center; }

private:
    CONST_PTR(ExposureT) const _exp;                   
    CONST_PTR(afw::detection::Footprint) const _foot;  
    geom::Point2D _center;                             
    geom::AffineTransform const _fromStandard;         
    geom::AffineTransform const _toStandard;           
};

template <typename ExposureT>
PTR(ExposurePatch<ExposureT>)
makeExposurePatch(CONST_PTR(ExposureT) exp, CONST_PTR(afw::detection::Footprint) foot,
                  geom::Point2D const& center) {
    return std::make_shared<ExposurePatch<ExposureT> >(exp, foot, center);
}
template <typename ExposureT>
PTR(ExposurePatch<ExposureT>)
makeExposurePatch(CONST_PTR(ExposureT) exp, afw::detection::Footprint const& standardFoot,
                  geom::Point2D const& standardCenter, afw::geom::SkyWcs const& standardWcs) {
    return std::make_shared<ExposurePatch<ExposureT> >(exp, standardFoot, standardCenter, standardWcs);
}

}  // namespace algorithms
}  // namespace meas
}  // namespace lsst

#endif