LSST Applications g04c3c9f7ca+2075667efa,g1e125bf412+5f448d5fcf,g2079a07aa2+3e9fd84d81,g2305ad1205+b635cf1488,g2bbee38e9b+6c6beb4891,g337abbeb29+6c6beb4891,g33d1c0ed96+6c6beb4891,g3a166c0a6a+6c6beb4891,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+42f171e1e6,g5c3423f6d4+d536b04327,g607f77f49a+d536b04327,g6f43f06aed+ca1339dc19,g858d7b2824+d536b04327,g8ee334c5b4+d7f9608c2f,g9963eaa53e+b3dc1655d3,g998f4353bf+d536b04327,g99cad8db69+8ef2408349,g9ddcbc5298+9a081db1e4,ga1e77700b3+2cbb763275,gadfd92a7e4+aec2f3b930,gae0086650b+585e252eca,gb0e22166c9+0e73c8378f,gb3b7280ab2+cb5fdb229e,gbb8dafda3b+a327199e22,gc120e1dc64+88074880ea,gc28159a63d+6c6beb4891,gcdd4ae20e8+bd241b2308,gcde1bda545+903e937d91,gcf0d15dbbd+bd241b2308,gdaeeff99f8+f9a426f77a,gddc38dedce+585e252eca,ge79ae78c31+6c6beb4891,gfbcc870c63+b310236976,w.2024.23
LSST Data Management Base Package
Loading...
Searching...
No Matches
Public Member Functions | Protected Attributes | List of all members
lsst::jointcal::TanSipPixelToRaDec Class Reference

Implements the (forward) SIP distorsion scheme. More...

#include <AstrometryTransform.h>

Inheritance diagram for lsst::jointcal::TanSipPixelToRaDec:
lsst::jointcal::BaseTanWcs lsst::jointcal::AstrometryTransform

Public Member Functions

 TanSipPixelToRaDec (AstrometryTransformLinear const &pixToTan, Point const &tangentPoint, const AstrometryTransformPolynomial *corrections=nullptr)
 pixToTan describes the transform from pix to tangent plane (degrees).
 
AstrometryTransformPolynomial getPixelToTangentPlane () const
 the transformation from pixels to tangent plane (degrees)
 
virtual void pixToTangentPlane (double xPixel, double yPixel, double &xTangentPlane, double &yTangentPlane) const
 transforms from pixel space to tangent plane (degrees)
 
 TanSipPixelToRaDec ()
 
std::unique_ptr< AstrometryTransforminverseTransform (double precision, const Frame &region) const
 Inverse transform: returns a TanRaDecToPixel if there are no corrections, or the iterative solver if there are.
 
std::unique_ptr< AstrometryTransformclone () const
 returns a copy (allocated by new) of the transformation.
 
void print (std::ostream &out) const
 prints the transform coefficients to stream.
 
double fit (StarMatchList const &starMatchList)
 Not implemented yet, because we do it otherwise.
 
void apply (double xIn, double yIn, double &xOut, double &yOut) const
 Transform pixels to ICRS RA, Dec in degrees.
 
virtual void apply (double xIn, double yIn, double &xOut, double &yOut) const=0
 
void apply (Point const &in, Point &out) const
 applies the tranfo to in and writes into out. Is indeed virtual.
 
Point apply (Point const &in) const
 All these apply(..) shadow the virtual one in derived classes, unless one writes "using AstrometryTransform::apply".
 
Frame apply (Frame const &inputframe, bool inscribed) const
 Transform a bounding box, taking either the inscribed or circumscribed box.
 
Point getTangentPoint () const
 Get the sky origin (CRVAL in FITS WCS terminology) in degrees.
 
AstrometryTransformLinear getLinPart () const
 The Linear part (corresponding to CD's and CRPIX's)
 
const AstrometryTransformPolynomialgetCorr () const
 Get a non-owning pointer to the correction transform polynomial.
 
void setCorrections (std::unique_ptr< AstrometryTransformPolynomial > corrections)
 Assign the correction polynomial (what it means is left to derived classes)
 
Point getCrPix () const
 Get the pixel origin of the WCS (CRPIX in FITS WCS terminology, but zero-based)
 
std::string __str__ () const
 
void transformStar (FatPoint &in) const
 
virtual double getJacobian (Point const &point) const
 returns the local jacobian.
 
virtual double getJacobian (double x, double y) const
 returns the local jacobian.
 
virtual std::unique_ptr< AstrometryTransformcomposeAndReduce (AstrometryTransform const &right) const
 Return a reduced composition of newTransform = this(right()), or nullptr if it cannot be reduced.
 
virtual void computeDerivative (Point const &where, AstrometryTransformLinear &derivative, double step=0.01) const
 Computes the local Derivative of a transform, w.r.t.
 
virtual AstrometryTransformLinear linearApproximation (Point const &where, double step=0.01) const
 linear (local) approximation.
 
virtual void transformPosAndErrors (const FatPoint &in, FatPoint &out) const
 
virtual void transformErrors (Point const &where, const double *vIn, double *vOut) const
 transform errors (represented as double[3] in order V(xx),V(yy),Cov(xy))
 
void getParams (double *params) const
 params should be at least Npar() long
 
void offsetParams (Eigen::VectorXd const &delta)
 
virtual double paramRef (Eigen::Index i) const
 
virtual double & paramRef (Eigen::Index i)
 
virtual void paramDerivatives (Point const &where, double *dx, double *dy) const
 Derivative w.r.t parameters.
 
virtual std::unique_ptr< AstrometryTransformroughInverse (const Frame &region) const
 Rough inverse.
 
virtual std::size_t getNpar () const
 returns the number of parameters (to compute chi2's)
 
virtual std::shared_ptr< ast::MappingtoAstMap (jointcal::Frame const &domain) const
 Create an equivalent AST mapping for this transformation, including an analytic inverse if possible.
 
void write (const std::string &fileName) const
 
virtual void write (std::ostream &stream) const
 

Protected Attributes

AstrometryTransformLinear linPixelToTan
 
std::unique_ptr< AstrometryTransformPolynomialcorr
 
double ra0 {}
 
double dec0 {}
 
double cos0 {}
 
double sin0 {}
 

Detailed Description

Implements the (forward) SIP distorsion scheme.

Definition at line 672 of file AstrometryTransform.h.

Constructor & Destructor Documentation

◆ TanSipPixelToRaDec() [1/2]

lsst::jointcal::TanSipPixelToRaDec::TanSipPixelToRaDec ( AstrometryTransformLinear const & pixToTan,
Point const & tangentPoint,
const AstrometryTransformPolynomial * corrections = nullptr )

pixToTan describes the transform from pix to tangent plane (degrees).

TangentPoint in degrees. Corrections are applied before Lin.

Definition at line 1589 of file AstrometryTransform.cc.

1591 : BaseTanWcs(pixToTan, tangentPoint, corrections) {}
BaseTanWcs(AstrometryTransformLinear const &pixToTan, Point const &tangentPoint, const AstrometryTransformPolynomial *corrections=nullptr)

◆ TanSipPixelToRaDec() [2/2]

lsst::jointcal::TanSipPixelToRaDec::TanSipPixelToRaDec ( )

Definition at line 1594 of file AstrometryTransform.cc.

1594: BaseTanWcs(AstrometryTransformLinear(), Point(0, 0), nullptr) {}

Member Function Documentation

◆ __str__()

std::string lsst::jointcal::AstrometryTransform::__str__ ( ) const
inlineinherited

Definition at line 94 of file AstrometryTransform.h.

94 {
96 print(s);
97 return s.str();
98 }
virtual void print(std::ostream &out) const =0
prints the transform coefficients to stream.

◆ apply() [1/5]

void lsst::jointcal::BaseTanWcs::apply ( double xIn,
double yIn,
double & xOut,
double & yOut ) const
virtualinherited

Transform pixels to ICRS RA, Dec in degrees.

Implements lsst::jointcal::AstrometryTransform.

Reimplemented in lsst::jointcal::TanPixelToRaDec.

Definition at line 1426 of file AstrometryTransform.cc.

1426 {
1427 double l, m; // radians in the tangent plane
1428 pixToTangentPlane(xIn, yIn, l, m); // l, m in degrees.
1429 l = deg2rad(l);
1430 m = deg2rad(m); // now in radians
1431 // Code inspired from worldpos.c in wcssubs (ancestor of the wcslib)
1432 /* At variance with wcslib, it collapses the projection to a plane
1433 and expression of sidereal cooordinates into a single set of
1434 operations. */
1435 double dect = cos0 - m * sin0;
1436 if (dect == 0) {
1437 LOGL_WARN(_log, "No sidereal coordinates at pole!");
1438 xOut = 0;
1439 yOut = 0;
1440 return;
1441 }
1442 double rat = ra0 + atan2(l, dect);
1443 dect = atan(std::cos(rat - ra0) * (m * cos0 + sin0) / dect);
1444 if (rat - ra0 > M_PI) rat -= (2. * M_PI);
1445 if (rat - ra0 < -M_PI) rat += (2. * M_PI);
1446 if (rat < 0.0) rat += (2. * M_PI);
1447 // convert to degree
1448 xOut = rad2deg(rat);
1449 yOut = rad2deg(dect);
1450}
#define M_PI
Definition ListMatch.cc:31
#define LOGL_WARN(logger, message...)
Log a warn-level message using a varargs/printf style interface.
Definition Log.h:547
int m
Definition SpanSet.cc:48
T atan2(T... args)
T atan(T... args)
virtual void pixToTangentPlane(double xPixel, double yPixel, double &xTangentPlane, double &yTangentPlane) const =0
Transform from pixels to tangent plane (degrees)
T cos(T... args)

◆ apply() [2/5]

virtual void lsst::jointcal::AstrometryTransform::apply ( double xIn,
double yIn,
double & xOut,
double & yOut ) const
virtualinherited

◆ apply() [3/5]

Frame lsst::jointcal::AstrometryTransform::apply ( Frame const & inputframe,
bool inscribed ) const
inherited

Transform a bounding box, taking either the inscribed or circumscribed box.

Parameters
[in]inputframeThe frame to be transformed.
[in]inscribedReturn the inscribed (true) or circumscribed (false) box.
Returns
The transformed frame.

Definition at line 89 of file AstrometryTransform.cc.

77 {
78 // 2 opposite corners
79 double xtmin1, xtmax1, ytmin1, ytmax1;
80 apply(inputframe.xMin, inputframe.yMin, xtmin1, ytmin1);
81 apply(inputframe.xMax, inputframe.yMax, xtmax1, ytmax1);
82 Frame fr1(std::min(xtmin1, xtmax1), std::min(ytmin1, ytmax1), std::max(xtmin1, xtmax1),
83 std::max(ytmin1, ytmax1));
84 // 2 other corners
85 double xtmin2, xtmax2, ytmin2, ytmax2;
86 apply(inputframe.xMin, inputframe.yMax, xtmin2, ytmax2);
87 apply(inputframe.xMax, inputframe.yMin, xtmax2, ytmin2);
88 Frame fr2(std::min(xtmin2, xtmax2), std::min(ytmin2, ytmax2), std::max(xtmin2, xtmax2),
89 std::max(ytmin2, ytmax2));
90
91 if (inscribed) return fr1 * fr2;
92 return fr1 + fr2;
93}
void apply(double xIn, double yIn, double &xOut, double &yOut) const
Transform pixels to ICRS RA, Dec in degrees.
T max(T... args)
T min(T... args)

◆ apply() [4/5]

Point lsst::jointcal::AstrometryTransform::apply ( Point const & in) const
inlineinherited

All these apply(..) shadow the virtual one in derived classes, unless one writes "using AstrometryTransform::apply".

Definition at line 75 of file AstrometryTransform.h.

75 {
76 double xout, yout;
77 apply(in.x, in.y, xout, yout);
78 return Point(xout, yout);
79 }

◆ apply() [5/5]

void lsst::jointcal::AstrometryTransform::apply ( Point const & in,
Point & out ) const
inlineinherited

applies the tranfo to in and writes into out. Is indeed virtual.

Definition at line 71 of file AstrometryTransform.h.

71{ apply(in.x, in.y, out.x, out.y); }

◆ clone()

std::unique_ptr< AstrometryTransform > lsst::jointcal::TanSipPixelToRaDec::clone ( ) const
virtual

returns a copy (allocated by new) of the transformation.

Implements lsst::jointcal::AstrometryTransform.

Definition at line 1630 of file AstrometryTransform.cc.

1630 {
1633}
std::unique_ptr< AstrometryTransformPolynomial > corr
AstrometryTransformLinear getLinPart() const
The Linear part (corresponding to CD's and CRPIX's)
Point getTangentPoint() const
Get the sky origin (CRVAL in FITS WCS terminology) in degrees.

◆ composeAndReduce()

std::unique_ptr< AstrometryTransform > lsst::jointcal::AstrometryTransform::composeAndReduce ( AstrometryTransform const & right) const
virtualinherited

Return a reduced composition of newTransform = this(right()), or nullptr if it cannot be reduced.

"Reduced" in this context means that they are capable of being merged into a single transform, for example, for two polynomials:

\[ f(x) = 1 + x^2, g(x) = -1 + 3x \]

we would have h = f.composeAndReduce(g) == 2 - 6x + 9x^2.

To be overloaded by derived classes if they can properly reduce the composition.

Parameters
rightThe transform to apply first.
Returns
The new reduced and composed AstrometryTransform, or nullptr if no such reduction is possible.

Reimplemented in lsst::jointcal::AstrometryTransformPolynomial, lsst::jointcal::TanPixelToRaDec, and lsst::jointcal::AstrometryTransformIdentity.

Definition at line 95 of file AstrometryTransform.cc.

96 { // by default no way to compose
98}

◆ computeDerivative()

void lsst::jointcal::AstrometryTransform::computeDerivative ( Point const & where,
AstrometryTransformLinear & derivative,
double step = 0.01 ) const
virtualinherited

Computes the local Derivative of a transform, w.r.t.

the Derivative is represented by a AstrometryTransformLinear, in which (hopefully), the offset terms are zero.

position.

Step is used for numerical derivation.

Derivative should transform a vector of offsets into a vector of offsets.

Reimplemented in lsst::jointcal::AstrometryTransformIdentity, lsst::jointcal::AstrometryTransformPolynomial, and lsst::jointcal::AstrometryTransformLinear.

Definition at line 119 of file AstrometryTransform.cc.

120 {
121 double x = where.x;
122 double y = where.y;
123 double xp0, yp0;
124 apply(x, y, xp0, yp0);
125
126 double xp, yp;
127 apply(x + step, y, xp, yp);
128 derivative.a11() = (xp - xp0) / step;
129 derivative.a21() = (yp - yp0) / step;
130 apply(x, y + step, xp, yp);
131 derivative.a12() = (xp - xp0) / step;
132 derivative.a22() = (yp - yp0) / step;
133 derivative.dx() = 0;
134 derivative.dy() = 0;
135}
int const step
int y
Definition SpanSet.cc:48
virtual void apply(double xIn, double yIn, double &xOut, double &yOut) const =0

◆ fit()

double lsst::jointcal::TanSipPixelToRaDec::fit ( StarMatchList const & starMatchList)
virtual

Not implemented yet, because we do it otherwise.

Implements lsst::jointcal::AstrometryTransform.

Definition at line 1644 of file AstrometryTransform.cc.

1644 {
1645 /* OK we could implement this routine, but it is
1646 probably useless since to do the match, we have to
1647 project from sky to tangent plane. When a match is
1648 found, it is easier to carry out the fit in the
1649 tangent plane, rather than going back to the celestial
1650 sphere (and reproject to fit...). Anyway if this
1651 message shows up, we'll think about it.
1652 */
1653 throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
1654 "TanSipPixelToRaDec::fit is NOT implemented (although it is doable)) ");
1655 return -1;
1656}
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48

◆ getCorr()

const AstrometryTransformPolynomial * lsst::jointcal::BaseTanWcs::getCorr ( ) const
inlineinherited

Get a non-owning pointer to the correction transform polynomial.

Definition at line 595 of file AstrometryTransform.h.

595{ return corr.get(); }

◆ getCrPix()

Point lsst::jointcal::BaseTanWcs::getCrPix ( ) const
inherited

Get the pixel origin of the WCS (CRPIX in FITS WCS terminology, but zero-based)

Definition at line 1460 of file AstrometryTransform.cc.

1460 {
1461 /* CRPIX's are defined by:
1462 ( CD1_1 CD1_2 ) (x - crpix1)
1463 transformed = ( ) * ( )
1464 ( CD2_1 CD2_2 ) (y - crpix2)
1465
1466 so that CrPix is the point which transforms to (0,0)
1467 */
1468 const AstrometryTransformLinear inverse = linPixelToTan.inverted();
1469 return Point(inverse.Dx(), inverse.Dy());
1470}
AstrometryTransformLinear inverted() const
returns the inverse: T1 = T2.inverted();
AstrometryTransformLinear linPixelToTan

◆ getJacobian() [1/2]

double lsst::jointcal::AstrometryTransform::getJacobian ( double x,
double y ) const
virtualinherited

returns the local jacobian.

Definition at line 100 of file AstrometryTransform.cc.

100 {
101 double x2, y2;
102 double eps = x * 0.01;
103 if (eps == 0) eps = 0.01;
104 apply(x, y, x2, y2);
105 double dxdx, dydx;
106 apply(x + eps, y, dxdx, dydx);
107 dxdx -= x2;
108 dydx -= y2;
109 double dxdy, dydy;
110 apply(x, y + eps, dxdy, dydy);
111 dxdy -= x2;
112 dydy -= y2;
113 return ((dxdx * dydy - dxdy * dydx) / (eps * eps));
114}

◆ getJacobian() [2/2]

virtual double lsst::jointcal::AstrometryTransform::getJacobian ( Point const & point) const
inlinevirtualinherited

returns the local jacobian.

Definition at line 110 of file AstrometryTransform.h.

110{ return getJacobian(point.x, point.y); }
virtual double getJacobian(Point const &point) const
returns the local jacobian.

◆ getLinPart()

AstrometryTransformLinear lsst::jointcal::BaseTanWcs::getLinPart ( ) const
inherited

The Linear part (corresponding to CD's and CRPIX's)

Definition at line 1454 of file AstrometryTransform.cc.

1454{ return linPixelToTan; }

◆ getNpar()

virtual std::size_t lsst::jointcal::AstrometryTransform::getNpar ( ) const
inlinevirtualinherited

◆ getParams()

void lsst::jointcal::AstrometryTransform::getParams ( double * params) const
inherited

params should be at least Npar() long

Definition at line 217 of file AstrometryTransform.cc.

217 {
218 std::size_t npar = getNpar();
219 for (std::size_t i = 0; i < npar; ++i) params[i] = paramRef(i);
220}
virtual std::size_t getNpar() const
returns the number of parameters (to compute chi2's)
virtual double paramRef(Eigen::Index i) const

◆ getPixelToTangentPlane()

AstrometryTransformPolynomial lsst::jointcal::TanSipPixelToRaDec::getPixelToTangentPlane ( ) const
virtual

the transformation from pixels to tangent plane (degrees)

Implements lsst::jointcal::BaseTanWcs.

Definition at line 1612 of file AstrometryTransform.cc.

1612 {
1613 if (corr)
1614 return AstrometryTransformPolynomial(linPixelToTan) * (*corr);
1615 else
1616 return linPixelToTan;
1617}

◆ getTangentPoint()

Point lsst::jointcal::BaseTanWcs::getTangentPoint ( ) const
inherited

Get the sky origin (CRVAL in FITS WCS terminology) in degrees.

Definition at line 1452 of file AstrometryTransform.cc.

1452{ return Point(rad2deg(ra0), rad2deg(dec0)); }

◆ inverseTransform()

std::unique_ptr< AstrometryTransform > lsst::jointcal::TanSipPixelToRaDec::inverseTransform ( double precision,
const Frame & region ) const
virtual

Inverse transform: returns a TanRaDecToPixel if there are no corrections, or the iterative solver if there are.

Reimplemented from lsst::jointcal::AstrometryTransform.

Definition at line 1606 of file AstrometryTransform.cc.

1608 { /* We have not implemented (yet) the reverse corrections available in SIP */
1609 return std::unique_ptr<AstrometryTransform>(new AstrometryTransformInverse(this, precision, region));
1610}

◆ linearApproximation()

AstrometryTransformLinear lsst::jointcal::AstrometryTransform::linearApproximation ( Point const & where,
double step = 0.01 ) const
virtualinherited

linear (local) approximation.

Reimplemented in lsst::jointcal::AstrometryTransformIdentity, and lsst::jointcal::AstrometryTransformLinear.

Definition at line 137 of file AstrometryTransform.cc.

138 {
139 Point outwhere = apply(where);
140 AstrometryTransformLinear der;
141 computeDerivative(where, der, step);
142 return AstrometryTransformLinearShift(outwhere.x, outwhere.y) * der *
143 AstrometryTransformLinearShift(-where.x, -where.y);
144}
virtual void computeDerivative(Point const &where, AstrometryTransformLinear &derivative, double step=0.01) const
Computes the local Derivative of a transform, w.r.t.

◆ offsetParams()

void lsst::jointcal::AstrometryTransform::offsetParams ( Eigen::VectorXd const & delta)
inherited

Definition at line 222 of file AstrometryTransform.cc.

222 {
223 std::size_t npar = getNpar();
224 for (std::size_t i = 0; i < npar; ++i) paramRef(i) += delta[i];
225}

◆ paramDerivatives()

void lsst::jointcal::AstrometryTransform::paramDerivatives ( Point const & where,
double * dx,
double * dy ) const
virtualinherited

Derivative w.r.t parameters.

Derivatives should be al least 2*NPar long. first Npar, for x, last Npar for y.

Reimplemented in lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 237 of file AstrometryTransform.cc.

237 {
238 throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
239 "AstrometryTransform::paramDerivatives() should never be called ");
240}

◆ paramRef() [1/2]

double & lsst::jointcal::AstrometryTransform::paramRef ( Eigen::Index i)
virtualinherited

Reimplemented in lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 232 of file AstrometryTransform.cc.

232 {
233 throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
234 "AstrometryTransform::paramRef should never be called ");
235}

◆ paramRef() [2/2]

double lsst::jointcal::AstrometryTransform::paramRef ( Eigen::Index i) const
virtualinherited

Reimplemented in lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 227 of file AstrometryTransform.cc.

227 {
228 throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
229 std::string("AstrometryTransform::paramRef should never be called "));
230}

◆ pixToTangentPlane()

void lsst::jointcal::TanSipPixelToRaDec::pixToTangentPlane ( double xPixel,
double yPixel,
double & xTangentPlane,
double & yTangentPlane ) const
virtual

transforms from pixel space to tangent plane (degrees)

Implements lsst::jointcal::BaseTanWcs.

Definition at line 1619 of file AstrometryTransform.cc.

1620 {
1621 // xTangentPlane, yTangentPlane returned in degrees
1622 if (corr) {
1623 double xtmp, ytmp;
1624 corr->apply(xPixel, yPixel, xtmp, ytmp);
1625 linPixelToTan.apply(xtmp, ytmp, xTangentPlane, yTangentPlane);
1626 } else
1627 linPixelToTan.apply(xPixel, yPixel, xTangentPlane, yTangentPlane);
1628}
void apply(double xIn, double yIn, double &xOut, double &yOut) const override

◆ print()

void lsst::jointcal::TanSipPixelToRaDec::print ( std::ostream & out) const
virtual

prints the transform coefficients to stream.

Implements lsst::jointcal::AstrometryTransform.

Definition at line 1635 of file AstrometryTransform.cc.

1635 {
1636 stream << " TanSipPixelToRaDec, lin part :" << endl << linPixelToTan;
1637 Point tp = getTangentPoint();
1638 stream << " tangent point " << tp.x << ' ' << tp.y << endl;
1639 Point crpix = getCrPix();
1640 stream << " crpix " << crpix.x << ' ' << crpix.y << endl;
1641 if (corr) stream << "PV correction: " << endl << *corr;
1642}
table::PointKey< double > crpix
Definition OldWcs.cc:129
Point getCrPix() const
Get the pixel origin of the WCS (CRPIX in FITS WCS terminology, but zero-based)
T endl(T... args)

◆ roughInverse()

std::unique_ptr< AstrometryTransform > lsst::jointcal::AstrometryTransform::roughInverse ( const Frame & region) const
virtualinherited

Rough inverse.

Stored by the numerical inverter to guess starting point for the trials. Just here to enable overloading.

Reimplemented in lsst::jointcal::AstrometryTransformInverse, lsst::jointcal::TanPixelToRaDec, and lsst::jointcal::TanRaDecToPixel.

Definition at line 196 of file AstrometryTransform.cc.

196 {
197 // "in" and "out" refer to the inverse direction.
198 Point centerOut = region.getCenter();
199 Point centerIn = apply(centerOut);
200 AstrometryTransformLinear der;
201 computeDerivative(centerOut, der, std::sqrt(region.getArea()) / 5.);
202 der = der.inverted();
203 der = AstrometryTransformLinearShift(centerOut.x, centerOut.y) * der *
204 AstrometryTransformLinearShift(-centerIn.x, -centerIn.y);
205 return std::unique_ptr<AstrometryTransform>(new AstrometryTransformLinear(der));
206}
T sqrt(T... args)

◆ setCorrections()

void lsst::jointcal::BaseTanWcs::setCorrections ( std::unique_ptr< AstrometryTransformPolynomial > corrections)
inherited

Assign the correction polynomial (what it means is left to derived classes)

Definition at line 1456 of file AstrometryTransform.cc.

1456 {
1457 corr = std::move(corrections);
1458}
T move(T... args)

◆ toAstMap()

virtual std::shared_ptr< ast::Mapping > lsst::jointcal::AstrometryTransform::toAstMap ( jointcal::Frame const & domain) const
inlinevirtualinherited

Create an equivalent AST mapping for this transformation, including an analytic inverse if possible.

Parameters
domainThe domain of the transform, to help find an inverse.
Returns
An AST Mapping that represents this transformation.

Reimplemented in lsst::jointcal::AstrometryTransformIdentity, and lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 189 of file AstrometryTransform.h.

189 {
190 throw std::logic_error("toAstMap is not implemented for this class.");
191 }

◆ transformErrors()

void lsst::jointcal::AstrometryTransform::transformErrors ( Point const & where,
const double * vIn,
double * vOut ) const
virtualinherited

transform errors (represented as double[3] in order V(xx),V(yy),Cov(xy))

Definition at line 163 of file AstrometryTransform.cc.

163 {
164 AstrometryTransformLinear der;
165 computeDerivative(where, der, 0.01);
166 double a11 = der.A11();
167 double a22 = der.A22();
168 double a21 = der.A21();
169 double a12 = der.A12();
170
171 /* (a11 a12) (vxx vxy)
172 M = ( ) and V = ( )
173 (a21 a22) (xvy vyy)
174
175 Vxx = Vin[0], vyy = Vin[1], Vxy = Vin[2];
176 we want to compute M*V*tp(M)
177 A lin alg light package would be perfect...
178 */
179 int xx = 0;
180 int yy = 1;
181 int xy = 2;
182 // M*V :
183
184 double b11 = a11 * vIn[xx] + a12 * vIn[xy];
185 double b22 = a21 * vIn[xy] + a22 * vIn[yy];
186 double b12 = a11 * vIn[xy] + a12 * vIn[yy];
187 double b21 = a21 * vIn[xx] + a22 * vIn[xy];
188
189 // (M*V) * tp(M)
190
191 vOut[xx] = b11 * a11 + b12 * a12;
192 vOut[xy] = b11 * a21 + b12 * a22;
193 vOut[yy] = b21 * a21 + b22 * a22;
194}

◆ transformPosAndErrors()

void lsst::jointcal::AstrometryTransform::transformPosAndErrors ( const FatPoint & in,
FatPoint & out ) const
virtualinherited

Reimplemented in lsst::jointcal::TanRaDecToPixel, and lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 146 of file AstrometryTransform.cc.

146 {
147 FatPoint res; // in case in and out are the same address...
148 res = apply(in);
149 AstrometryTransformLinear der;
150 // could save a call here, since Derivative needs the transform of where that we already have
151 // 0.01 may not be a very good idea in all cases. May be we should provide a way of altering that.
152 computeDerivative(in, der, 0.01);
153 double a11 = der.A11();
154 double a22 = der.A22();
155 double a21 = der.A21();
156 double a12 = der.A12();
157 res.vx = a11 * (a11 * in.vx + 2 * a12 * in.vxy) + a12 * a12 * in.vy;
158 res.vy = a21 * a21 * in.vx + a22 * a22 * in.vy + 2. * a21 * a22 * in.vxy;
159 res.vxy = a21 * a11 * in.vx + a22 * a12 * in.vy + (a21 * a12 + a11 * a22) * in.vxy;
160 out = res;
161}

◆ transformStar()

void lsst::jointcal::AstrometryTransform::transformStar ( FatPoint & in) const
inlineinherited

Definition at line 107 of file AstrometryTransform.h.

107{ transformPosAndErrors(in, in); }
virtual void transformPosAndErrors(const FatPoint &in, FatPoint &out) const

◆ write() [1/2]

void lsst::jointcal::AstrometryTransform::write ( const std::string & fileName) const
inherited

Definition at line 247 of file AstrometryTransform.cc.

247 {
248 ofstream s(fileName.c_str());
249 write(s);
250 bool ok = !s.fail();
251 s.close();
252 if (!ok)
253 throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
254 "AstrometryTransform::write, something went wrong for file " + fileName);
255}
T c_str(T... args)
void write(const std::string &fileName) const

◆ write() [2/2]

void lsst::jointcal::AstrometryTransform::write ( std::ostream & stream) const
virtualinherited

Reimplemented in lsst::jointcal::AstrometryTransformIdentity, and lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 257 of file AstrometryTransform.cc.

257 {
258 throw LSST_EXCEPT(
259 pex::exceptions::InvalidParameterError,
260 "AstrometryTransform::write(ostream), should never be called. MEans that it is missing in some "
261 "derived class ");
262}

Member Data Documentation

◆ corr

std::unique_ptr<AstrometryTransformPolynomial> lsst::jointcal::BaseTanWcs::corr
protectedinherited

Definition at line 616 of file AstrometryTransform.h.

◆ cos0

double lsst::jointcal::BaseTanWcs::cos0 {}
protectedinherited

Definition at line 618 of file AstrometryTransform.h.

618{}, sin0{}; // cos(dec0), sin(dec0)

◆ dec0

double lsst::jointcal::BaseTanWcs::dec0 {}
protectedinherited

Definition at line 617 of file AstrometryTransform.h.

617{}, dec0{}; // sky origin (radians)

◆ linPixelToTan

AstrometryTransformLinear lsst::jointcal::BaseTanWcs::linPixelToTan
protectedinherited

Definition at line 614 of file AstrometryTransform.h.

◆ ra0

double lsst::jointcal::BaseTanWcs::ra0 {}
protectedinherited

Definition at line 617 of file AstrometryTransform.h.

617{}, dec0{}; // sky origin (radians)

◆ sin0

double lsst::jointcal::BaseTanWcs::sin0 {}
protectedinherited

Definition at line 618 of file AstrometryTransform.h.

618{}, sin0{}; // cos(dec0), sin(dec0)

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