BuildSingleKernelVisitor.h

Go to the documentation of this file.

// -*- lsst-c++ -*-
#ifndef LSST_IP_DIFFIM_BUILDSINGLEKERNELVISITOR_H
#define LSST_IP_DIFFIM_BUILDSINGLEKERNELVISITOR_H
 
#include <memory>
 
#include "lsst/afw/image.h"
#include "lsst/afw/math.h"
 
#include "lsst/daf/base/PropertySet.h"
 
#include "lsst/ip/diffim/ImageStatistics.h"
 
namespace lsst {
namespace ip {
namespace diffim {
namespace detail {
 
    template<typename PixelT>
    class BuildSingleKernelVisitor : public lsst::afw::math::CandidateVisitor {
        typedef lsst::afw::image::MaskedImage<PixelT> MaskedImageT;
    public:
        typedef std::shared_ptr<BuildSingleKernelVisitor<PixelT> > Ptr;
 
        BuildSingleKernelVisitor(
            lsst::afw::math::KernelList const& basisList,
            lsst::daf::base::PropertySet const& ps
            );
        BuildSingleKernelVisitor(
            lsst::afw::math::KernelList const& basisList,
            lsst::daf::base::PropertySet const& ps,
            Eigen::MatrixXd const& hMat
            );
        virtual ~BuildSingleKernelVisitor() {};
 
        /*
           Don't reprocess candidate if its already been build.  The use
           case for this functionality is : when iterating over all Cells
           and rejecting bad Kernels, we need to re-visit *all* Cells to
           build the next candidate in the list.  Without this flag we would
           unncessarily re-build all the good Kernels.
        */
        void setSkipBuilt(bool skip)      {_skipBuilt = skip;}
 
        int getNRejected()    {return _nRejected;}
        int getNProcessed()   {return _nProcessed;}
        void reset()          {_nRejected = 0; _nProcessed = 0;}
 
        void processCandidate(lsst::afw::math::SpatialCellCandidate *candidate);
 
    private:
        lsst::afw::math::KernelList const _basisList; 
        lsst::daf::base::PropertySet::Ptr _ps; 
        Eigen::MatrixXd const _hMat;          
        ImageStatistics<PixelT> _imstats;     
        bool _skipBuilt;                      
        int _nRejected;                       
        int _nProcessed;                      
        bool _useRegularization;              
 
        bool _useCoreStats;                   
        int _coreRadius;                      
    };
 
    template<typename PixelT>
    std::shared_ptr<BuildSingleKernelVisitor<PixelT> >
    makeBuildSingleKernelVisitor(
        lsst::afw::math::KernelList const& basisList,
        lsst::daf::base::PropertySet const& ps
        ) {
 
        return std::shared_ptr<BuildSingleKernelVisitor<PixelT>>(
            new BuildSingleKernelVisitor<PixelT>(basisList, ps)
            );
    }
 
    template<typename PixelT>
    std::shared_ptr<BuildSingleKernelVisitor<PixelT> >
    makeBuildSingleKernelVisitor(
        lsst::afw::math::KernelList const& basisList,
        lsst::daf::base::PropertySet const& ps,
        Eigen::MatrixXd const & hMat
        ) {
 
        return std::shared_ptr<BuildSingleKernelVisitor<PixelT>>(
            new BuildSingleKernelVisitor<PixelT>(basisList, ps, hMat)
            );
    }
 
}}}} // end of namespace lsst::ip::diffim::detail
 
#endif