FootprintSet.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_DETECTION_FOOTPRINT_SET_H)
#define LSST_DETECTION_FOOTPRINT_SET_H

#include "lsst/afw/geom.h"
#include "lsst/afw/detection/Threshold.h"
#include "lsst/afw/detection/Footprint.h"
#include "lsst/afw/detection/FootprintCtrl.h"
#include "lsst/afw/image/MaskedImage.h"
#include "lsst/afw/table/Source.h"

namespace lsst {
namespace afw {
namespace detection {

typedef boost::uint64_t FootprintIdPixel;

/************************************************************************************************************/
class FootprintSet : public lsst::daf::base::Citizen {
public:

    typedef std::vector<Footprint::Ptr> FootprintList;

#ifndef SWIG
    template <typename ImagePixelT>
    FootprintSet(image::Image<ImagePixelT> const& img,
                 Threshold const& threshold,
                 int const npixMin=1, bool const setPeaks=true);

    template <typename MaskPixelT>
    FootprintSet(image::Mask<MaskPixelT> const& img,
                 Threshold const& threshold,
                 int const npixMin=1);

    template <typename ImagePixelT, typename MaskPixelT>
    FootprintSet(image::MaskedImage<ImagePixelT, MaskPixelT> const& img,
                 Threshold const& threshold,
                 std::string const& planeName = "",
                 int const npixMin=1, bool const setPeaks=true);

    template <typename ImagePixelT, typename MaskPixelT>
    FootprintSet(image::MaskedImage<ImagePixelT, MaskPixelT> const& img,
                 Threshold const& threshold,
                 int x,
                 int y,
                 std::vector<PTR(Peak)> const* peaks = NULL);

#else // workaround for https://github.com/swig/swig/issues/245
    // if that bug is fixed then you may update footprintset.i by uncommenting two lines
    // and removing this section. However, you must continue to provide SWIG
    // the alternate version of the template <typename MaskPixelT> constructor, because
    // SWIG cannot disambiguate that from the template <typename ImagePixelT> constructor.

    FootprintSet(image::Mask<image::MaskPixel> const& img,
                 Threshold const& threshold,
                 int const npixMin=1);

    FootprintSet(image::Image<boost::uint16_t> const& img,
                 Threshold const& threshold,
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<boost::uint16_t, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 std::string const& planeName = "",
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<boost::uint16_t, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 int x,
                 int y,
                 std::vector<PTR(Peak)> const* peaks = NULL);

    FootprintSet(image::Image<int> const& img,
                 Threshold const& threshold,
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<int, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 std::string const& planeName = "",
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<int, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 int x,
                 int y,
                 std::vector<PTR(Peak)> const* peaks = NULL);

    FootprintSet(image::Image<float> const& img,
                 Threshold const& threshold,
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<float, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 std::string const& planeName = "",
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<float, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 int x,
                 int y,
                 std::vector<PTR(Peak)> const* peaks = NULL);

    FootprintSet(image::Image<double> const& img,
                 Threshold const& threshold,
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<double, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 std::string const& planeName = "",
                 int const npixMin=1, bool const setPeaks=true);
    FootprintSet(image::MaskedImage<double, image::MaskPixel> const& img,
                 Threshold const& threshold,
                 int x,
                 int y,
                 std::vector<PTR(Peak)> const* peaks = NULL);

#endif

    FootprintSet(geom::Box2I region);
    FootprintSet(FootprintSet const&);
    FootprintSet(FootprintSet const& set, int rGrow, FootprintControl const& ctrl);
    FootprintSet(FootprintSet const& set, int rGrow, bool isotropic=true);
    FootprintSet(FootprintSet const& footprints1, 
                 FootprintSet const& footprints2,
                 bool const includePeaks);

    FootprintSet& operator=(FootprintSet const& rhs);

    void swap(FootprintSet& rhs) {
        using std::swap;                    // See Meyers, Effective C++, Item 25        
        swap(*_footprints, *rhs.getFootprints());
        geom::Box2I rhsRegion = rhs.getRegion();
        rhs.setRegion(getRegion());
        setRegion(rhsRegion);
    }
    
    PTR(FootprintList) getFootprints() { return _footprints; } 

    void setFootprints(PTR(FootprintList) footprints) { _footprints = footprints; }

    CONST_PTR(FootprintList) const getFootprints() const { return _footprints; }
    
    void makeSources(afw::table::SourceCatalog & catalog) const;

    void setRegion(geom::Box2I const& region);

    geom::Box2I const getRegion() const { return _region; } 

    PTR(image::Image<FootprintIdPixel>) insertIntoImage(
        const bool relativeIDs
        ) const;

    template <typename MaskPixelT>
    void setMask(
        image::Mask<MaskPixelT> *mask, 
        std::string const& planeName   
    ) {
        setMaskFromFootprintList(
            mask, 
            _footprints,                // calling getFootprints() confuses clang++ 3.0 and leaks memory
            image::Mask<MaskPixelT>::getPlaneBitMask(planeName)
        );        
    }

    template <typename MaskPixelT>
    void setMask(
        PTR(image::Mask<MaskPixelT>) mask, 
        std::string const& planeName   
    ) {
        setMask(mask.get(), planeName);
    }

    void merge(FootprintSet const& rhs, int tGrow=0, int rGrow=0, bool isotropic=true);

    template <typename ImagePixelT, typename MaskPixelT>
    void makeHeavy(image::MaskedImage<ImagePixelT, MaskPixelT> const& mimg,
                   HeavyFootprintCtrl const* ctrl=NULL
                  );
private:
    boost::shared_ptr<FootprintList> _footprints;        
    geom::Box2I _region;                
};

}}}

#endif