GaussianPsf.h

Go to the documentation of this file.

// -*- LSST-C++ -*-
/*
 * LSST Data Management System
 * Copyright 2008-2014 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/>.
 */
#ifndef LSST_AFW_DETECTION_GaussianPsf_h_INCLUDED
#define LSST_AFW_DETECTION_GaussianPsf_h_INCLUDED

#include "lsst/afw/geom/Extent.h"
#include "lsst/afw/geom/Point.h"
#include "lsst/afw/detection/Psf.h"

namespace lsst { namespace afw { namespace detection {

class GaussianPsf : public afw::table::io::PersistableFacade<GaussianPsf>, public Psf {
public:

    GaussianPsf(int width, int height, double sigma);

    GaussianPsf(geom::Extent2I const & dimensions, double sigma);

    virtual PTR(afw::detection::Psf) clone() const;

    geom::Extent2I getDimensions() const { return _dimensions; }

    double getSigma() const { return _sigma; }

    virtual bool isPersistable() const { return true; }

protected:

    virtual std::string getPersistenceName() const;

    virtual std::string getPythonModule() const;

    virtual void write(OutputArchiveHandle & handle) const;

private:

#if 0 // We could reimplement this more efficiently than what's in the base class,
      // but it's tricky to get the position right in all corner cases, and it's
      // not actually performance-critical, so we should just wait for #3116.
    virtual PTR(Image) doComputeImage(
        geom::Point2D const & position, image::Color const& color
    ) const;
#endif

    virtual PTR(Image) doComputeKernelImage(
        geom::Point2D const & position, image::Color const & color
    ) const;

    virtual double doComputeApertureFlux(
        double radius, geom::Point2D const & position, image::Color const & color
    ) const;

    virtual geom::ellipses::Quadrupole doComputeShape(
        geom::Point2D const & position, image::Color const & color
    ) const;

    virtual geom::Box2I doComputeBBox(
        geom::Point2D const & position, image::Color const & color
    ) const;

    geom::Extent2I _dimensions;
    double _sigma;
};

}}} // namespace lsst::afw::detection

#endif // !LSST_AFW_DETECTION_Psf_h_INCLUDED