lsst::jointcal::ConstrainedPolyModel Class Reference

This is the model used to fit mappings as the combination of a transformation depending on the chip number (instrument model) and a transformation per visit (anamorphism). More...

#include <ConstrainedPolyModel.h>

Inheritance diagram for lsst::jointcal::ConstrainedPolyModel:

Public Member Functions
	ConstrainedPolyModel (CcdImageList const &ccdImageList, ProjectionHandler const *projectionHandler, bool initFromWCS, unsigned nNotFit=0, int chipDegree=3, int visitDegree=2)
	ConstrainedPolyModel (ConstrainedPolyModel const &)=delete
	No copy or move: there is only ever one instance of a given model (i.e. per ccd+visit) More...
	ConstrainedPolyModel (ConstrainedPolyModel &&)=delete
ConstrainedPolyModel &	operator= (ConstrainedPolyModel const &)=delete
ConstrainedPolyModel &	operator= (ConstrainedPolyModel &&)=delete
Mapping const *	getMapping (CcdImage const &) const
	Mapping associated to a given CcdImage. More...
unsigned	assignIndices (unsigned firstIndex, std::string const &whatToFit)
	Positions the various parameter sets into the parameter vector, starting at firstIndex. More...
void	offsetParams (Eigen::VectorXd const &Delta)
	Dispaches the offsets after a fit step into the actual locations of parameters. More...
void	freezeErrorTransform ()
	From there on, measurement errors are propagated using the current transfos (and no longer evolve). More...
Gtransfo const &	getChipTransfo (CcdIdType const chip) const
	Access to mappings. More...
Gtransfo const &	getVisitTransfo (VisitIdType const &visit) const
	Access to mappings. More...
std::vector< VisitIdType >	getVisits () const
	Access to array of visits involved in the solution. More...
const Gtransfo *	getSky2TP (CcdImage const &ccdImage) const
	The mapping of sky coordinates (i.e. More...
std::shared_ptr< TanSipPix2RaDec >	produceSipWcs (CcdImage const &ccdImage) const
	Cook up a SIP WCS. More...

Detailed Description

This is the model used to fit mappings as the combination of a transformation depending on the chip number (instrument model) and a transformation per visit (anamorphism).

The two-transformation Mapping required for this model is TwoTransfoMapping. This modeling of distortions is meant for a set of images from a single mosaic imager.

Definition at line 31 of file ConstrainedPolyModel.h.

Constructor & Destructor Documentation

ConstrainedPolyModel() [1/3]

lsst::jointcal::ConstrainedPolyModel::ConstrainedPolyModel	(	CcdImageList const &	ccdImageList,
		ProjectionHandler const *	projectionHandler,
		bool	initFromWCS,
		unsigned	nNotFit = 0,
		int	chipDegree = 3,
		int	visitDegree = 2
	)

Definition at line 30 of file ConstrainedPolyModel.cc.

        : _sky2TP(projectionHandler)

{
    // first loop to initialize all visit  and chip transfos.
    for (auto &ccdImage : ccdImageList) {
        const CcdImage &im = *ccdImage;
        auto visit = im.getVisit();
        auto chip = im.getCcdId();
        auto visitp = _visitMap.find(visit);
        if (visitp == _visitMap.end()) {
            if (_visitMap.size() == 0) {
                _visitMap[visit] =
                        std::unique_ptr<SimpleGtransfoMapping>(new SimpleGtransfoMapping(GtransfoIdentity()));
            } else {
                _visitMap[visit] = std::unique_ptr<SimpleGtransfoMapping>(
                        new SimplePolyMapping(GtransfoLin(), GtransfoPoly(visitDegree)));
            }
        }
        auto chipp = _chipMap.find(chip);
        if (chipp == _chipMap.end()) {
            const Frame &frame = im.getImageFrame();

            GtransfoPoly pol(chipDegree);

            if (initFromWCS) {
                pol = GtransfoPoly(im.getPix2TangentPlane(), frame, chipDegree);
            }
            GtransfoLin shiftAndNormalize = normalizeCoordinatesTransfo(frame);

            _chipMap[chip] = std::unique_ptr<SimplePolyMapping>(
                    new SimplePolyMapping(shiftAndNormalize, pol * shiftAndNormalize.invert()));
        }
    }
    // now, second loop to set the mappings of the CCdImages
    for (auto &ccdImage : ccdImageList) {
        const CcdImage &im = *ccdImage;
        auto visit = im.getVisit();
        auto chip = im.getCcdId();
        // check that the chip_indexed part was indeed assigned
        // (i.e. the reference visit was complete)
        if (_chipMap.find(chip) == _chipMap.end()) {
            LOGLS_WARN(_log, "Chip " << chip << " is missing in the reference exposure, expect troubles.");
            GtransfoLin norm = normalizeCoordinatesTransfo(im.getImageFrame());
            _chipMap[chip] =
                    std::unique_ptr<SimplePolyMapping>(new SimplePolyMapping(norm, GtransfoPoly(chipDegree)));
        }
        _mappings[&im] = std::unique_ptr<TwoTransfoMapping>(
                new TwoTransfoMapping(_chipMap[chip].get(), _visitMap[visit].get()));
    }
    LOGLS_INFO(_log, "Constructor got " << _chipMap.size() << " chip mappings and " << _visitMap.size()
                                        << " visit mappings.");
    // DEBUG
    for (auto i = _visitMap.begin(); i != _visitMap.end(); ++i) LOGLS_DEBUG(_log, i->first);
}

ConstrainedPolyModel() [2/3]

lsst::jointcal::ConstrainedPolyModel::ConstrainedPolyModel ( ConstrainedPolyModel const &  )

delete

No copy or move: there is only ever one instance of a given model (i.e. per ccd+visit)

ConstrainedPolyModel() [3/3]

lsst::jointcal::ConstrainedPolyModel::ConstrainedPolyModel ( ConstrainedPolyModel &&  )

delete

Member Function Documentation

assignIndices()

unsigned lsst::jointcal::ConstrainedPolyModel::assignIndices ( unsigned  firstIndex, std::string const &  whatToFit  )

virtual

Positions the various parameter sets into the parameter vector, starting at firstIndex.

This routine decodes "DistortionsChip" and "DistortionsVisit" in whatToFit.

If whatToFit contains "Distortions" and not Distortions<Something>, it is understood as both chips and visits.

Implements lsst::jointcal::AstrometryModel.

Definition at line 98 of file ConstrainedPolyModel.cc.

                                                                                            {
    unsigned index = firstIndex;
    if (whatToFit.find("Distortions") == std::string::npos) {
        LOGLS_ERROR(_log, "assignIndices was called and Distortions is *not* in whatToFit");
        return 0;
    }
    // if we get here "Distortions" is in whatToFit
    _fittingChips = (whatToFit.find("DistortionsChip") != std::string::npos);
    _fittingVisits = (whatToFit.find("DistortionsVisit") != std::string::npos);
    // If nothing more than "Distortions" is specified, it means all:
    if ((!_fittingChips) && (!_fittingVisits)) {
        _fittingChips = _fittingVisits = true;
    }
    if (_fittingChips)
        for (auto &i : _chipMap) {
            i.second->setIndex(index);
            index += i.second->getNpar();
        }
    if (_fittingVisits)
        for (auto &i : _visitMap) {
            i.second->setIndex(index);
            index += i.second->getNpar();
        }
    // Tell the mappings which derivatives they will have to fill:
    for (auto &i : _mappings) {
        i.second->setWhatToFit(_fittingChips, _fittingVisits);
    }
    return index;
}

freezeErrorTransform()

void lsst::jointcal::ConstrainedPolyModel::freezeErrorTransform ( )

virtual

From there on, measurement errors are propagated using the current transfos (and no longer evolve).

Implements lsst::jointcal::AstrometryModel.

Definition at line 141 of file ConstrainedPolyModel.cc.

                                                {
    for (auto i = _visitMap.begin(); i != _visitMap.end(); ++i) i->second->freezeErrorTransform();
    for (auto i = _chipMap.begin(); i != _chipMap.end(); ++i) i->second->freezeErrorTransform();
}

getChipTransfo()

const Gtransfo & lsst::jointcal::ConstrainedPolyModel::getChipTransfo

(

CcdIdType const

chip

)

const

Access to mappings.

Definition at line 146 of file ConstrainedPolyModel.cc.

                                                                               {
    auto chipp = _chipMap.find(chip);
    if (chipp == _chipMap.end()) {
        std::stringstream errMsg;
        errMsg << "No such chipId: '" << chip << "' found in chipMap of:  " << this;
        throw pexExcept::InvalidParameterError(errMsg.str());
    }
    return chipp->second->getTransfo();
}

getMapping()

const Mapping * lsst::jointcal::ConstrainedPolyModel::getMapping ( CcdImage const &  ) const

virtual

Mapping associated to a given CcdImage.

Implements lsst::jointcal::AstrometryModel.

Definition at line 88 of file ConstrainedPolyModel.cc.

                                                                              {
    mappingMapType::const_iterator i = _mappings.find(&ccdImage);
    if (i == _mappings.end()) return nullptr;
    return (i->second.get());
}

getSky2TP()

const Gtransfo* lsst::jointcal::ConstrainedPolyModel::getSky2TP ( CcdImage const &  ccdImage ) const

inlinevirtual

The mapping of sky coordinates (i.e.

the coordinate system in which fitted stars are reported) onto the Tangent plane (into which the pixel coordinates are transformed).

Implements lsst::jointcal::AstrometryModel.

Definition at line 78 of file ConstrainedPolyModel.h.

{ return _sky2TP->getSky2TP(ccdImage); }

getVisits()

std::vector< VisitIdType > lsst::jointcal::ConstrainedPolyModel::getVisits

(

)

const

Access to array of visits involved in the solution.

Definition at line 157 of file ConstrainedPolyModel.cc.

                                                             {
    std::vector<VisitIdType> res;
    res.reserve(_visitMap.size());
    for (auto i = _visitMap.begin(); i != _visitMap.end(); ++i) res.push_back(i->first);
    return res;
}

getVisitTransfo()

const Gtransfo & lsst::jointcal::ConstrainedPolyModel::getVisitTransfo

(

VisitIdType const &

visit

)

const

Access to mappings.

Definition at line 164 of file ConstrainedPolyModel.cc.

                                                                                    {
    auto visitp = _visitMap.find(visit);
    if (visitp == _visitMap.end()) {
        std::stringstream errMsg;
        errMsg << "No such visitId: '" << visit << "' found in visitMap of: " << this;
        throw pexExcept::InvalidParameterError(errMsg.str());
    }
    return visitp->second->getTransfo();
}

offsetParams()

void lsst::jointcal::ConstrainedPolyModel::offsetParams ( Eigen::VectorXd const &  Delta )

virtual

Dispaches the offsets after a fit step into the actual locations of parameters.

Implements lsst::jointcal::AstrometryModel.

Definition at line 128 of file ConstrainedPolyModel.cc.

                                                                  {
    if (_fittingChips)
        for (auto &i : _chipMap) {
            auto mapping = i.second.get();
            mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
        }
    if (_fittingVisits)
        for (auto &i : _visitMap) {
            auto mapping = i.second.get();
            mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
        }
}

operator=() [1/2]

ConstrainedPolyModel& lsst::jointcal::ConstrainedPolyModel::operator= ( ConstrainedPolyModel const &  )

delete

operator=() [2/2]

ConstrainedPolyModel& lsst::jointcal::ConstrainedPolyModel::operator= ( ConstrainedPolyModel &&  )

delete

produceSipWcs()

std::shared_ptr< TanSipPix2RaDec > lsst::jointcal::ConstrainedPolyModel::produceSipWcs ( CcdImage const &  ccdImage ) const

virtual

Cook up a SIP WCS.

Implements lsst::jointcal::AstrometryModel.

Definition at line 174 of file ConstrainedPolyModel.cc.

                                                                                                 {
    const TwoTransfoMapping *mapping;
    try {
        mapping = _mappings.at(&ccdImage).get();
    } catch (std::out_of_range &) {
        LOGLS_ERROR(_log, "CcdImage with ccd/visit " << ccdImage.getCcdId() << "/" << ccdImage.getVisit()
                                                     << " not found in constrainedPolyModel mapping list.");
        std::ostringstream os;
        for (auto const &i : _mappings) os << i.first << ",";
        LOGLS_ERROR(_log, "Available CcdImages: " << os.str());
        return nullptr;
    }

    GtransfoPoly pix2Tp;
    const GtransfoPoly &t1 = dynamic_cast<const GtransfoPoly &>(mapping->getTransfo1());
    // TODO: This line produces a warning on clang (t1 is always valid: a failed dynamic_cast of a reference
    // raises bad_cast instead of returning nullptr like a failed pointer cast), but I'll deal with it as
    // part of DM-10524 (hopefully removing the necessity of the casts).
    if (!(&t1)) {
        LOGLS_ERROR(_log, "Problem with transform 1 of ccd/visit " << ccdImage.getCcdId() << "/"
                                                                   << ccdImage.getVisit() << ": T1 "
                                                                   << mapping->getTransfo1());
        return nullptr;
    }
    // NOTE: we currently expect T2 to be an identity for the first visit, so we have to treat it separately.
    // TODO: We are aware that this is a hack, but it will be fixed as part of DM-10524.
    try {
        const GtransfoIdentity &t2 = dynamic_cast<const GtransfoIdentity &>(mapping->getTransfo2());
        pix2Tp = t1;
    } catch (std::bad_cast &) {
        try {
            const GtransfoPoly &t2_poly = dynamic_cast<const GtransfoPoly &>(mapping->getTransfo2());
            pix2Tp = t2_poly * t1;
        } catch (std::bad_cast &) {
            LOGLS_ERROR(_log, "Problem with transform 2 of ccd/visit " << ccdImage.getCcdId() << "/"
                                                                       << ccdImage.getVisit() << ": T2 "
                                                                       << mapping->getTransfo2());
            return nullptr;
        }
    }
    const TanRaDec2Pix *proj = dynamic_cast<const TanRaDec2Pix *>(getSky2TP(ccdImage));
    if (!proj) {
        LOGLS_ERROR(_log, "Problem with projection of ccd/visit " << ccdImage.getCcdId() << "/"
                                                                  << ccdImage.getVisit() << ": projection "
                                                                  << getSky2TP(ccdImage));
        return nullptr;
    }

    // should be the identity, but who knows? So, let us incorporate it into the pix2TP part.
    const GtransfoLin &projLinPart = proj->getLinPart();
    GtransfoPoly wcsPix2Tp = GtransfoPoly(projLinPart.invert()) * pix2Tp;

    // compute a decent approximation, if higher order corrections get ignored
    GtransfoLin cdStuff = wcsPix2Tp.linearApproximation(ccdImage.getImageFrame().getCenter());

    // wcsPix2TP = cdStuff*sip , so
    GtransfoPoly sip = GtransfoPoly(cdStuff.invert()) * wcsPix2Tp;
    Point tangentPoint(proj->getTangentPoint());
    return std::make_shared<TanSipPix2RaDec>(cdStuff, tangentPoint, &sip);
}

---

The documentation for this class was generated from the following files:

• /home/jenkins-slave/snowflake/release/lsstsw/stack/Linux64/jointcal/15.0-6-ge2d9597+10/include/lsst/jointcal/ConstrainedPolyModel.h
• /home/jenkins-slave/snowflake/release/lsstsw/stack/Linux64/jointcal/15.0-6-ge2d9597+10/src/ConstrainedPolyModel.cc