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/afw/detection/Footprint.h"
#include "lsst/afw/image/Wcs.h"
#include "lsst/afw/geom/AffineTransform.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, 
                  afw::geom::Point2D const& center           
        ): _exp(exp), _foot(foot), _center(center), _fromStandard(), _toStandard() {}
    ExposurePatch(CONST_PTR(ExposureT) exp, 
                  afw::detection::Footprint const& standardFoot, 
                  afw::geom::Point2D const& standardCenter,  
                  afw::image::Wcs const& standardWcs         
        ) : _exp(exp) {
        afw::image::Wcs const& expWcs = *exp->getWcs();
        afw::coord::Coord::ConstPtr sky = standardWcs.pixelToSky(standardCenter);
        const_cast<CONST_PTR(afw::detection::Footprint)&>(_foot) = standardFoot.transform(standardWcs, expWcs,
                                                                                          exp->getBBox());
        const_cast<afw::geom::Point2D&>(_center) = expWcs.skyToPixel(*sky);
        const_cast<afw::geom::AffineTransform&>(_fromStandard) = standardWcs.linearizePixelToSky(*sky) *
            expWcs.linearizeSkyToPixel(*sky);
        const_cast<afw::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; }
    afw::geom::Point2D const& getCenter() const { return _center; }
    afw::geom::AffineTransform const& fromStandard() const { return _fromStandard; }
    afw::geom::AffineTransform const& toStandard() const { return _toStandard; }

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

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

template<typename ExposureT>
PTR(ExposurePatch<ExposureT>) makeExposurePatch(
    CONST_PTR(ExposureT) exp,
    CONST_PTR(afw::detection::Footprint) foot,
    afw::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,
    afw::geom::Point2D const& standardCenter,
    afw::image::Wcs const& standardWcs
    ) {
    return std::make_shared<ExposurePatch<ExposureT> >(exp, standardFoot, standardCenter, standardWcs);
}

}}} // namespace

#endif