SimplePolyModel.h

Go to the documentation of this file.

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

#include "memory"

#include "lsst/jointcal/Eigenstuff.h"

#include "lsst/jointcal/AstrometryModel.h"
#include "lsst/jointcal/Gtransfo.h"
#include "lsst/jointcal/SimplePolyMapping.h"
#include "lsst/jointcal/ProjectionHandler.h"
#include <map>

namespace lsst {
namespace jointcal {

class CcdImage;

/* We deal here with coordinate transforms which are fitted
   and/or necessary to AstrometryFit. The classes SimplePolyModel
and SimplePolyMapping implement a model where  there is one
separate transfrom per CcdImage. One could chose other setups.

*/

/* This modeling of distortions can even accommodate images set mixing instruments */
class SimplePolyModel : public AstrometryModel {
public:
    SimplePolyModel(CcdImageList const &ccdImageList, ProjectionHandler const *projectionHandler,
                    bool initFromWCS, unsigned nNotFit = 0, unsigned degree = 3);

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

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

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

    // dispaches the offsets after a fit step into the actual locations of parameters
    void offsetParams(Eigen::VectorXd const &delta);

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

    virtual void freezeErrorTransform();

    Gtransfo const &getTransfo(CcdImage const &ccdImage) const;

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

    ~SimplePolyModel(){};

private:
    typedef std::map<const CcdImage *, std::unique_ptr<SimpleGtransfoMapping>> mapType;
    mapType _myMap;
    const ProjectionHandler *_sky2TP;
};
}  // namespace jointcal
}  // namespace lsst

#endif  // LSST_JOINTCAL_SIMPLE_POLY_MODEL_H