ConstrainedPolyModel.h

Go to the documentation of this file.

// -*- LSST-C++ -*-
#ifndef LSST_JOINTCAL_CONSTRAINED_POLY_MODEL_H
#define LSST_JOINTCAL_CONSTRAINED_POLY_MODEL_H

#include "memory"  // for std::*_ptr

#include "lsst/jointcal/Eigenstuff.h"

class CcdImage;

#include "lsst/jointcal/AstrometryModel.h"
#include "lsst/jointcal/Gtransfo.h"
#include "lsst/jointcal/Frame.h"
#include "lsst/jointcal/SimplePolyMapping.h"
#include "lsst/jointcal/ProjectionHandler.h"
#include "lsst/jointcal/TwoTransfoMapping.h"
#include "lsst/jointcal/CcdImage.h"

#include <map>

namespace lsst {
namespace jointcal {

class ConstrainedPolyModel : public AstrometryModel {
public:
    ConstrainedPolyModel(CcdImageList const &ccdImageList, ProjectionHandler const *projectionHandler,
                         bool initFromWCS, unsigned nNotFit = 0, int chipDegree = 3, int visitDegree = 2);

    ConstrainedPolyModel(ConstrainedPolyModel const &) = delete;
    ConstrainedPolyModel(ConstrainedPolyModel &&) = delete;
    ConstrainedPolyModel &operator=(ConstrainedPolyModel const &) = delete;
    ConstrainedPolyModel &operator=(ConstrainedPolyModel &&) = delete;

    // The following routines are the interface to AstrometryFit
    Mapping const *getMapping(CcdImage const &) const;

    unsigned assignIndices(unsigned firstIndex, std::string const &whatToFit);

    void offsetParams(Eigen::VectorXd const &Delta);

    void freezeErrorTransform();

    Gtransfo const &getChipTransfo(CcdIdType const chip) const;

    Gtransfo const &getVisitTransfo(VisitIdType const &visit) const;

    std::vector<VisitIdType> getVisits() const;

    const Gtransfo *getSky2TP(CcdImage const &ccdImage) const { return _sky2TP->getSky2TP(ccdImage); }

    std::shared_ptr<TanSipPix2RaDec> produceSipWcs(CcdImage const &ccdImage) const;

private:
    typedef std::map<const CcdImage *, std::unique_ptr<TwoTransfoMapping>> mappingMapType;
    mappingMapType _mappings;
    typedef std::map<CcdIdType, std::unique_ptr<SimpleGtransfoMapping>> chipMapType;
    chipMapType _chipMap;
    typedef std::map<VisitIdType, std::unique_ptr<SimpleGtransfoMapping>> visitMapType;
    visitMapType _visitMap;
    const ProjectionHandler *_sky2TP;
    bool _fittingChips, _fittingVisits;
};
}  // namespace jointcal
}  // namespace lsst

#endif  // LSST_JOINTCAL_CONSTRAINED_POLY_MODEL_H