LSST Applications g2079a07aa2+86d27d4dc4,g2305ad1205+a659bff248,g2bbee38e9b+3c60f8fe34,g337abbeb29+3c60f8fe34,g33d1c0ed96+3c60f8fe34,g3502564af9+d77d6d1350,g3a166c0a6a+3c60f8fe34,g487adcacf7+25d9892218,g4be5004598+d77d6d1350,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+4d81263f9a,g5cd07815a0+980d2b1c3b,g607f77f49a+d77d6d1350,g858d7b2824+d77d6d1350,g88963caddf+83e433e629,g99cad8db69+a4d3c48eeb,g9ddcbc5298+9a081db1e4,ga1e77700b3+bcf1af89ad,ga57fefb910+9a39d7b2d7,gae0086650b+585e252eca,gb065fddaf9+4f9fd82a2c,gb0e22166c9+60f28cb32d,gb363559e06+d84b1d3d07,gb3b7280ab2+4563d032e1,gb4b16eec92+babe958938,gba4ed39666+c2a2e4ac27,gbb8dafda3b+ed6854b564,gc120e1dc64+b72d212f87,gc28159a63d+3c60f8fe34,gc3e9b769f7+921dbcd359,gcf0d15dbbd+9a39d7b2d7,gdaeeff99f8+f9a426f77a,gddc38dedce+585e252eca,ge79ae78c31+3c60f8fe34,w.2024.21
LSST Data Management Base Package
Loading...
Searching...
No Matches
Public Member Functions | Protected Attributes | List of all members
lsst::jointcal::SimplePolyMapping Class Reference

Mapping implementation for a polynomial transformation. More...

#include <SimpleAstrometryMapping.h>

Inheritance diagram for lsst::jointcal::SimplePolyMapping:
lsst::jointcal::SimpleAstrometryMapping lsst::jointcal::AstrometryMapping

Public Member Functions

 ~SimplePolyMapping ()=default
 
 SimplePolyMapping (AstrometryTransformLinear CenterAndScale, AstrometryTransformPolynomial const &transform)
 The transformation will be initialized to transform, so that the effective transformation reads transform*CenterAndScale.
 
 SimplePolyMapping (SimplePolyMapping const &)=delete
 No copy or move: there is only ever one instance of a given mapping (i.e.. per ccd+visit)
 
 SimplePolyMapping (SimplePolyMapping &&)=delete
 
SimplePolyMappingoperator= (SimplePolyMapping const &)=delete
 
SimplePolyMappingoperator= (SimplePolyMapping &&)=delete
 
void positionDerivative (Point const &where, Eigen::Matrix2d &derivative, double epsilon) const override
 The derivative w.r.t. position.
 
void computeTransformAndDerivatives (FatPoint const &where, FatPoint &outPoint, Eigen::MatrixX2d &H) const override
 Calls the transforms and implements the centering and scaling of coordinates.
 
void transformPosAndErrors (FatPoint const &where, FatPoint &outPoint) const override
 The same as above but without the parameter derivatives (used to evaluate chi^2)
 
AstrometryTransform const & getTransform () const override
 Access to the (fitted) transform.
 
virtual void freezeErrorTransform ()
 
std::size_t getNpar () const override
 Number of parameters in total.
 
void getMappingIndices (IndexVector &indices) const override
 Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has enough space reserved.
 
void offsetParams (Eigen::VectorXd const &delta) override
 Remember the error scale and freeze it.
 
Eigen::Index getIndex () const
 position of the parameters within the grand fitting scheme
 
void setIndex (Eigen::Index i)
 Set the index of this mapping in the grand fit.
 
bool getToBeFit () const
 Get whether this mapping is fit as part of a Model.
 
void setToBeFit (bool value)
 Set whether this Mapping is to be fit as part of a Model.
 
void print (std::ostream &out) const override
 Print a string representation of the contents of this mapping, for debugging.
 

Protected Attributes

bool toBeFit
 
Eigen::Index index {}
 
std::shared_ptr< AstrometryTransformtransform
 
std::shared_ptr< AstrometryTransformerrorProp
 
std::unique_ptr< AstrometryTransformLinearlin
 

Detailed Description

Mapping implementation for a polynomial transformation.

Definition at line 118 of file SimpleAstrometryMapping.h.

Constructor & Destructor Documentation

◆ ~SimplePolyMapping()

lsst::jointcal::SimplePolyMapping::~SimplePolyMapping ( )
default

◆ SimplePolyMapping() [1/3]

lsst::jointcal::SimplePolyMapping::SimplePolyMapping ( AstrometryTransformLinear CenterAndScale,
AstrometryTransformPolynomial const & transform )

The transformation will be initialized to transform, so that the effective transformation reads transform*CenterAndScale.

Definition at line 75 of file SimpleAstrometryMapping.cc.

77 : SimpleAstrometryMapping(transform), _centerAndScale(std::move(CenterAndScale)) {
78 // We assume that the initialization was done properly, for example that
79 // transform = pixToTangentPlane*CenterAndScale.inverted(), so we do not touch transform.
80 /* store the (spatial) derivative of _centerAndScale. For the extra
81 diagonal terms, just copied the ones in positionDerivatives */
82 preDer(0, 0) = _centerAndScale.getCoefficient(1, 0, 0);
83 preDer(1, 0) = _centerAndScale.getCoefficient(0, 1, 0);
84 preDer(0, 1) = _centerAndScale.getCoefficient(1, 0, 1);
85 preDer(1, 1) = _centerAndScale.getCoefficient(0, 1, 1);
86
87 // check of matrix indexing (once for all)
88 MatrixX2d H(3, 2);
89 assert((&H(1, 0) - &H(0, 0)) == 1);
90}
Eigen::Matrix< double, Eigen::Dynamic, 2 > MatrixX2d
Definition Eigenstuff.h:33
double getCoefficient(std::size_t powX, std::size_t powY, std::size_t whichCoord) const
Get the coefficient of a given power in x and y, for either the x or y coordinate.
SimpleAstrometryMapping(AstrometryTransform const &astrometryTransform, bool toBeFit=true)
T move(T... args)

◆ SimplePolyMapping() [2/3]

lsst::jointcal::SimplePolyMapping::SimplePolyMapping ( SimplePolyMapping const & )
delete

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

◆ SimplePolyMapping() [3/3]

lsst::jointcal::SimplePolyMapping::SimplePolyMapping ( SimplePolyMapping && )
delete

Member Function Documentation

◆ computeTransformAndDerivatives()

void lsst::jointcal::SimplePolyMapping::computeTransformAndDerivatives ( FatPoint const & where,
FatPoint & outPoint,
Eigen::MatrixX2d & H ) const
overridevirtual

Calls the transforms and implements the centering and scaling of coordinates.

Reimplemented from lsst::jointcal::SimpleAstrometryMapping.

Definition at line 109 of file SimpleAstrometryMapping.cc.

110 {
111 FatPoint mid;
112 _centerAndScale.transformPosAndErrors(where, mid);
113 transform->transformPosAndErrors(mid, outPoint);
114 FatPoint tmp;
115 errorProp->transformPosAndErrors(mid, tmp);
116 outPoint.vx = tmp.vx;
117 outPoint.vy = tmp.vy;
118 outPoint.vxy = tmp.vxy;
119 transform->paramDerivatives(mid, &H(0, 0), &H(0, 1));
120}
virtual void transformPosAndErrors(const FatPoint &in, FatPoint &out) const override
a mix of apply and Derivative
std::shared_ptr< AstrometryTransform > errorProp

◆ freezeErrorTransform()

virtual void lsst::jointcal::SimpleAstrometryMapping::freezeErrorTransform ( )
inlinevirtualinherited

Definition at line 56 of file SimpleAstrometryMapping.h.

56 {
57 // from there on, updating the transform does not change the errors.
58 errorProp = transform->clone();
59 }

◆ getIndex()

Eigen::Index lsst::jointcal::SimpleAstrometryMapping::getIndex ( ) const
inlineinherited

position of the parameters within the grand fitting scheme

Definition at line 84 of file SimpleAstrometryMapping.h.

◆ getMappingIndices()

void lsst::jointcal::SimpleAstrometryMapping::getMappingIndices ( IndexVector & indices) const
overridevirtualinherited

Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has enough space reserved.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 34 of file SimpleAstrometryMapping.cc.

34 {
35 if (indices.size() < getNpar()) {
36 indices.resize(getNpar());
37 }
38 for (std::size_t k = 0; k < getNpar(); ++k) {
39 indices[k] = index + k;
40 }
41}
std::size_t getNpar() const override
Number of parameters in total.
T resize(T... args)
T size(T... args)

◆ getNpar()

std::size_t lsst::jointcal::SimpleAstrometryMapping::getNpar ( ) const
inlineoverridevirtualinherited

Number of parameters in total.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 62 of file SimpleAstrometryMapping.h.

62 {
63 if (toBeFit)
64 return transform->getNpar();
65 else
66 return 0;
67 }

◆ getToBeFit()

bool lsst::jointcal::SimpleAstrometryMapping::getToBeFit ( ) const
inlineinherited

Get whether this mapping is fit as part of a Model.

Definition at line 97 of file SimpleAstrometryMapping.h.

97{ return toBeFit; }

◆ getTransform()

AstrometryTransform const & lsst::jointcal::SimplePolyMapping::getTransform ( ) const
overridevirtual

Access to the (fitted) transform.

Reimplemented from lsst::jointcal::SimpleAstrometryMapping.

Definition at line 133 of file SimpleAstrometryMapping.cc.

133 {
134 // Cannot fail given the contructor:
135 const auto *fittedPoly =
136 dynamic_cast<const AstrometryTransformPolynomial *>(&(*transform));
137 actualResult = (*fittedPoly) * _centerAndScale;
138 return actualResult;
139}

◆ offsetParams()

void lsst::jointcal::SimpleAstrometryMapping::offsetParams ( Eigen::VectorXd const & delta)
inlineoverridevirtualinherited

Remember the error scale and freeze it.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 79 of file SimpleAstrometryMapping.h.

79 {
80 if (toBeFit) transform->offsetParams(delta);
81 }

◆ operator=() [1/2]

SimplePolyMapping & lsst::jointcal::SimplePolyMapping::operator= ( SimplePolyMapping && )
delete

◆ operator=() [2/2]

SimplePolyMapping & lsst::jointcal::SimplePolyMapping::operator= ( SimplePolyMapping const & )
delete

◆ positionDerivative()

void lsst::jointcal::SimplePolyMapping::positionDerivative ( Point const & where,
Eigen::Matrix2d & derivative,
double epsilon ) const
overridevirtual

The derivative w.r.t. position.

Reimplemented from lsst::jointcal::SimpleAstrometryMapping.

Definition at line 92 of file SimpleAstrometryMapping.cc.

93 {
94 Point tmp = _centerAndScale.apply(where);
95 errorProp->computeDerivative(tmp, *lin, epsilon);
96 derivative(0, 0) = lin->getCoefficient(1, 0, 0);
97 //
98 /* This does not work : it was proved by rotating the frame
99 see the compilation switch ROTATE_T2 in constrainedAstrometryModel.cc
100 derivative(1,0) = lin->getCoefficient(1,0,1);
101 derivative(0,1) = lin->getCoefficient(0,1,0);
102 */
103 derivative(1, 0) = lin->getCoefficient(0, 1, 0);
104 derivative(0, 1) = lin->getCoefficient(1, 0, 1);
105 derivative(1, 1) = lin->getCoefficient(0, 1, 1);
106 derivative = preDer * derivative;
107}
void apply(double xIn, double yIn, double &xOut, double &yOut) const override
std::unique_ptr< AstrometryTransformLinear > lin

◆ print()

void lsst::jointcal::SimpleAstrometryMapping::print ( std::ostream & out) const
overridevirtualinherited

Print a string representation of the contents of this mapping, for debugging.

This string representation can be very verbose, as it contains all of the parameters of all of the transforms in this mapping.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 73 of file SimpleAstrometryMapping.cc.

73{ out << *transform; }
std::shared_ptr< AstrometryTransform > transform

◆ setIndex()

void lsst::jointcal::SimpleAstrometryMapping::setIndex ( Eigen::Index i)
inlineinherited

Set the index of this mapping in the grand fit.

Definition at line 87 of file SimpleAstrometryMapping.h.

87{ index = i; }

◆ setToBeFit()

void lsst::jointcal::SimpleAstrometryMapping::setToBeFit ( bool value)
inlineinherited

Set whether this Mapping is to be fit as part of a Model.

Definition at line 99 of file SimpleAstrometryMapping.h.

◆ transformPosAndErrors()

void lsst::jointcal::SimplePolyMapping::transformPosAndErrors ( FatPoint const & where,
FatPoint & outPoint ) const
overridevirtual

The same as above but without the parameter derivatives (used to evaluate chi^2)

Reimplemented from lsst::jointcal::SimpleAstrometryMapping.

Definition at line 122 of file SimpleAstrometryMapping.cc.

122 {
123 FatPoint mid;
124 _centerAndScale.transformPosAndErrors(where, mid);
125 transform->transformPosAndErrors(mid, outPoint);
126 FatPoint tmp;
127 errorProp->transformPosAndErrors(mid, tmp);
128 outPoint.vx = tmp.vx;
129 outPoint.vy = tmp.vy;
130 outPoint.vxy = tmp.vxy;
131}

Member Data Documentation

◆ errorProp

std::shared_ptr<AstrometryTransform> lsst::jointcal::SimpleAstrometryMapping::errorProp
protectedinherited

Definition at line 111 of file SimpleAstrometryMapping.h.

◆ index

Eigen::Index lsst::jointcal::SimpleAstrometryMapping::index {}
protectedinherited

Definition at line 106 of file SimpleAstrometryMapping.h.

106{};

◆ lin

std::unique_ptr<AstrometryTransformLinear> lsst::jointcal::SimpleAstrometryMapping::lin
protectedinherited

Definition at line 114 of file SimpleAstrometryMapping.h.

◆ toBeFit

bool lsst::jointcal::SimpleAstrometryMapping::toBeFit
protectedinherited

Definition at line 105 of file SimpleAstrometryMapping.h.

◆ transform

std::shared_ptr<AstrometryTransform> lsst::jointcal::SimpleAstrometryMapping::transform
protectedinherited

Definition at line 109 of file SimpleAstrometryMapping.h.


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