LSST Applications  21.0.0-172-gfb10e10a+18fedfabac,22.0.0+297cba6710,22.0.0+80564b0ff1,22.0.0+8d77f4f51a,22.0.0+a28f4c53b1,22.0.0+dcf3732eb2,22.0.1-1-g7d6de66+2a20fdde0d,22.0.1-1-g8e32f31+297cba6710,22.0.1-1-geca5380+7fa3b7d9b6,22.0.1-12-g44dc1dc+2a20fdde0d,22.0.1-15-g6a90155+515f58c32b,22.0.1-16-g9282f48+790f5f2caa,22.0.1-2-g92698f7+dcf3732eb2,22.0.1-2-ga9b0f51+7fa3b7d9b6,22.0.1-2-gd1925c9+bf4f0e694f,22.0.1-24-g1ad7a390+a9625a72a8,22.0.1-25-g5bf6245+3ad8ecd50b,22.0.1-25-gb120d7b+8b5510f75f,22.0.1-27-g97737f7+2a20fdde0d,22.0.1-32-gf62ce7b1+aa4237961e,22.0.1-4-g0b3f228+2a20fdde0d,22.0.1-4-g243d05b+871c1b8305,22.0.1-4-g3a563be+32dcf1063f,22.0.1-4-g44f2e3d+9e4ab0f4fa,22.0.1-42-gca6935d93+ba5e5ca3eb,22.0.1-5-g15c806e+85460ae5f3,22.0.1-5-g58711c4+611d128589,22.0.1-5-g75bb458+99c117b92f,22.0.1-6-g1c63a23+7fa3b7d9b6,22.0.1-6-g50866e6+84ff5a128b,22.0.1-6-g8d3140d+720564cf76,22.0.1-6-gd805d02+cc5644f571,22.0.1-8-ge5750ce+85460ae5f3,master-g6e05de7fdc+babf819c66,master-g99da0e417a+8d77f4f51a,w.2021.48
LSST Data Management Base Package
Public Member Functions | Protected Attributes | List of all members
lsst::jointcal::TanPixelToRaDec Class Reference

The transformation that handles pixels to sideral transformations (Gnomonic, possibly with polynomial distortions). More...

#include <AstrometryTransform.h>

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

Public Member Functions

 TanPixelToRaDec (AstrometryTransformLinear const &pixToTan, Point const &tangentPoint, const AstrometryTransformPolynomial *corrections=nullptr)
 pixToTan describes the transform from pix to tangent plane (degrees). More...
 
AstrometryTransformPolynomial getPixelToTangentPlane () const
 the transformation from pixels to tangent plane (degrees) More...
 
virtual void pixToTangentPlane (double xPixel, double yPixel, double &xTangentPlane, double &yTangentPlane) const
 transforms from pixel space to tangent plane (degrees) More...
 
 TanPixelToRaDec ()
 
TanPixelToRaDec operator* (AstrometryTransformLinear const &right) const
 composition with AstrometryTransformLinear More...
 
std::unique_ptr< AstrometryTransformcomposeAndReduce (AstrometryTransformLinear const &right) const
 Return a reduced composition of newTransform = this(right()), or nullptr if it cannot be reduced. More...
 
TanRaDecToPixel inverted () const
 approximate inverse : it ignores corrections; More...
 
std::unique_ptr< AstrometryTransformroughInverse (const Frame &region) const
 Overload the "generic routine" (available for all AstrometryTransform types. More...
 
std::unique_ptr< AstrometryTransforminverseTransform (const double precision, const Frame &region) const
 Inverse transform: returns a TanRaDecToPixel if there are no corrections, or the iterative solver if there are. More...
 
std::unique_ptr< AstrometryTransformclone () const
 returns a copy (allocated by new) of the transformation. More...
 
void print (std::ostream &out) const
 prints the transform coefficients to stream. More...
 
double fit (StarMatchList const &starMatchList)
 Not implemented yet, because we do it otherwise. More...
 
virtual void apply (const double xIn, const 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. More...
 
Point apply (Point const &in) const
 All these apply(..) shadow the virtual one in derived classes, unless one writes "using AstrometryTransform::apply". More...
 
Frame apply (Frame const &inputframe, bool inscribed) const
 Transform a bounding box, taking either the inscribed or circumscribed box. More...
 
virtual std::unique_ptr< AstrometryTransformcomposeAndReduce (AstrometryTransform const &right) const
 Return a reduced composition of newTransform = this(right()), or nullptr if it cannot be reduced. More...
 
void apply (const double xIn, const double yIn, double &xOut, double &yOut) const
 Transform pixels to ICRS RA, Dec in degrees. More...
 
void apply (Point const &in, Point &out) const
 applies the tranfo to in and writes into out. Is indeed virtual. More...
 
Point apply (Point const &in) const
 All these apply(..) shadow the virtual one in derived classes, unless one writes "using AstrometryTransform::apply". More...
 
Frame apply (Frame const &inputframe, bool inscribed) const
 Transform a bounding box, taking either the inscribed or circumscribed box. More...
 
Point getTangentPoint () const
 Get the sky origin (CRVAL in FITS WCS terminology) in degrees. More...
 
AstrometryTransformLinear getLinPart () const
 The Linear part (corresponding to CD's and CRPIX's) More...
 
const AstrometryTransformPolynomialgetCorr () const
 Get a non-owning pointer to the correction transform polynomial. More...
 
void setCorrections (std::unique_ptr< AstrometryTransformPolynomial > corrections)
 Assign the correction polynomial (what it means is left to derived classes) More...
 
Point getCrPix () const
 Get the pixel origin of the WCS (CRPIX in FITS WCS terminology, but zero-based) More...
 
std::string __str__ ()
 
void transformStar (FatPoint &in) const
 
virtual double getJacobian (Point const &point) const
 returns the local jacobian. More...
 
virtual double getJacobian (const double x, const double y) const
 returns the local jacobian. More...
 
virtual std::unique_ptr< AstrometryTransformcomposeAndReduce (AstrometryTransform const &right) const
 Return a reduced composition of newTransform = this(right()), or nullptr if it cannot be reduced. More...
 
virtual void computeDerivative (Point const &where, AstrometryTransformLinear &derivative, const double step=0.01) const
 Computes the local Derivative of a transform, w.r.t. More...
 
virtual AstrometryTransformLinear linearApproximation (Point const &where, const double step=0.01) const
 linear (local) approximation. More...
 
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)) More...
 
void getParams (double *params) const
 params should be at least Npar() long More...
 
void offsetParams (Eigen::VectorXd const &delta)
 
virtual double paramRef (Eigen::Index const i) const
 
virtual double & paramRef (Eigen::Index const i)
 
virtual void paramDerivatives (Point const &where, double *dx, double *dy) const
 Derivative w.r.t parameters. More...
 
virtual std::size_t getNpar () const
 returns the number of parameters (to compute chi2's) More...
 
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. More...
 
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

The transformation that handles pixels to sideral transformations (Gnomonic, possibly with polynomial distortions).

Definition at line 627 of file AstrometryTransform.h.

Constructor & Destructor Documentation

◆ TanPixelToRaDec() [1/2]

lsst::jointcal::TanPixelToRaDec::TanPixelToRaDec ( 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 between Lin and deprojection parts (as in Swarp).

Definition at line 1497 of file AstrometryTransform.cc.

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

◆ TanPixelToRaDec() [2/2]

lsst::jointcal::TanPixelToRaDec::TanPixelToRaDec ( )

Definition at line 1502 of file AstrometryTransform.cc.

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

Member Function Documentation

◆ __str__()

std::string lsst::jointcal::AstrometryTransform::__str__ ( )
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.
T str(T... args)

◆ apply() [1/8]

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

Transform pixels to ICRS RA, Dec in degrees.

Implements lsst::jointcal::AstrometryTransform.

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/8]

virtual void lsst::jointcal::AstrometryTransform::apply

◆ apply() [3/8]

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 76 of file AstrometryTransform.cc.

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

◆ apply() [4/8]

Frame lsst::jointcal::AstrometryTransform::apply

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.

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

◆ apply() [5/8]

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() [6/8]

Point lsst::jointcal::AstrometryTransform::apply
inline

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() [7/8]

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); }

◆ apply() [8/8]

void lsst::jointcal::AstrometryTransform::apply
inline

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::TanPixelToRaDec::clone ( ) const
virtual

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

Implements lsst::jointcal::AstrometryTransform.

Definition at line 1559 of file AstrometryTransform.cc.

1559  {
1561  new TanPixelToRaDec(getLinPart(), getTangentPoint(), corr.get()));
1562 }
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() [1/3]

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::AstrometryTransformIdentity.

Definition at line 94 of file AstrometryTransform.cc.

95  { // by default no way to compose
97 }

◆ composeAndReduce() [2/3]

std::unique_ptr< AstrometryTransform > lsst::jointcal::AstrometryTransform::composeAndReduce

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.

Definition at line 131 of file AstrometryTransform.cc.

95  { // by default no way to compose
97 }

◆ composeAndReduce() [3/3]

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

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.

Definition at line 1504 of file AstrometryTransform.cc.

1505  {
1506  if (right.getOrder() == 1) {
1507  return std::make_unique<TanPixelToRaDec>((*this) * (right));
1508  } else {
1509  return std::unique_ptr<AstrometryTransform>(nullptr);
1510  }
1511 }
T right(T... args)

◆ computeDerivative()

void lsst::jointcal::AstrometryTransform::computeDerivative ( Point const &  where,
AstrometryTransformLinear derivative,
const 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::AstrometryTransformLinear, lsst::jointcal::AstrometryTransformPolynomial, and lsst::jointcal::AstrometryTransformIdentity.

Definition at line 118 of file AstrometryTransform.cc.

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

◆ fit()

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

Not implemented yet, because we do it otherwise.

Implements lsst::jointcal::AstrometryTransform.

Definition at line 1573 of file AstrometryTransform.cc.

1573  {
1574  /* OK we could implement this routine, but it is
1575  probably useless since to do the match, we have to
1576  project from sky to tangent plane. When a match is
1577  found, it is easier to carry out the fit in the
1578  tangent plane, rather than going back to the celestial
1579  sphere (and reproject to fit...). Anyway if this
1580  message shows up, we'll think about it.
1581  */
1582  throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
1583  "TanPixelToRaDec::fit is NOT implemented (although it is doable)) ");
1584  return -1;
1585 }
#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 ( const double  x,
const double  y 
) const
virtualinherited

returns the local jacobian.

Definition at line 99 of file AstrometryTransform.cc.

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

◆ 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 216 of file AstrometryTransform.cc.

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

◆ getPixelToTangentPlane()

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

the transformation from pixels to tangent plane (degrees)

Implements lsst::jointcal::BaseTanWcs.

Definition at line 1541 of file AstrometryTransform.cc.

1541  {
1542  if (corr)
1543  return (*corr) * linPixelToTan;
1544  else
1545  return linPixelToTan;
1546 }

◆ 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::TanPixelToRaDec::inverseTransform ( const 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 1532 of file AstrometryTransform.cc.

1533  {
1534  if (!corr)
1536  new TanRaDecToPixel(getLinPart().inverted(), getTangentPoint()));
1537  else
1538  return std::unique_ptr<AstrometryTransform>(new AstrometryTransformInverse(this, precision, region));
1539 }
TanRaDecToPixel inverted() const
approximate inverse : it ignores corrections;

◆ inverted()

TanRaDecToPixel lsst::jointcal::TanPixelToRaDec::inverted ( ) const

approximate inverse : it ignores corrections;

Definition at line 1519 of file AstrometryTransform.cc.

1519  {
1520  if (corr != nullptr) {
1521  LOGL_WARN(_log, "You are inverting a TanPixelToRaDec with corrections.");
1522  LOGL_WARN(_log, "The inverse you get ignores the corrections!");
1523  }
1524  return TanRaDecToPixel(getLinPart().inverted(), getTangentPoint());
1525 }

◆ linearApproximation()

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

linear (local) approximation.

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

Definition at line 136 of file AstrometryTransform.cc.

137  {
138  Point outwhere = apply(where);
139  AstrometryTransformLinear der;
140  computeDerivative(where, der, step);
141  return AstrometryTransformLinearShift(outwhere.x, outwhere.y) * der *
142  AstrometryTransformLinearShift(-where.x, -where.y);
143 }
virtual void computeDerivative(Point const &where, AstrometryTransformLinear &derivative, const 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 221 of file AstrometryTransform.cc.

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

◆ operator*()

TanPixelToRaDec lsst::jointcal::TanPixelToRaDec::operator* ( AstrometryTransformLinear const &  right) const

composition with AstrometryTransformLinear

Definition at line 1513 of file AstrometryTransform.cc.

1513  {
1514  TanPixelToRaDec result(*this);
1515  result.linPixelToTan = result.linPixelToTan * right;
1516  return result;
1517 }
py::object result
Definition: _schema.cc:429

◆ 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 236 of file AstrometryTransform.cc.

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

◆ paramRef() [1/2]

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

Reimplemented in lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 231 of file AstrometryTransform.cc.

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

◆ paramRef() [2/2]

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

Reimplemented in lsst::jointcal::AstrometryTransformPolynomial.

Definition at line 226 of file AstrometryTransform.cc.

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

◆ pixToTangentPlane()

void lsst::jointcal::TanPixelToRaDec::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 1548 of file AstrometryTransform.cc.

1549  {
1550  // xTangentPlane, yTangentPlane in degrees.
1551  linPixelToTan.apply(xPixel, yPixel, xTangentPlane, yTangentPlane);
1552  if (corr) {
1553  double xtmp = xTangentPlane;
1554  double ytmp = yTangentPlane;
1555  corr->apply(xtmp, ytmp, xTangentPlane, yTangentPlane); // still in degrees.
1556  }
1557 }
void apply(const double xIn, const double yIn, double &xOut, double &yOut) const override

◆ print()

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

prints the transform coefficients to stream.

Implements lsst::jointcal::AstrometryTransform.

Definition at line 1564 of file AstrometryTransform.cc.

1564  {
1565  stream << " TanPixelToRaDec, lin part :" << endl << linPixelToTan;
1566  Point tp = getTangentPoint();
1567  stream << " tangent point " << tp.x << ' ' << tp.y << endl;
1568  Point crpix = getCrPix();
1569  stream << " crpix " << crpix.x << ' ' << crpix.y << endl;
1570  if (corr) stream << "PV correction: " << endl << *corr;
1571 }
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::TanPixelToRaDec::roughInverse ( const Frame region) const
virtual

Overload the "generic routine" (available for all AstrometryTransform types.

Reimplemented from lsst::jointcal::AstrometryTransform.

Definition at line 1527 of file AstrometryTransform.cc.

1527  {
1529  new TanRaDecToPixel(getLinPart().inverted(), getTangentPoint()));
1530 }

◆ 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::AstrometryTransformPolynomial, and lsst::jointcal::AstrometryTransformIdentity.

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 162 of file AstrometryTransform.cc.

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

◆ transformPosAndErrors()

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

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

Definition at line 145 of file AstrometryTransform.cc.

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

◆ 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 246 of file AstrometryTransform.cc.

246  {
247  ofstream s(fileName.c_str());
248  write(s);
249  bool ok = !s.fail();
250  s.close();
251  if (!ok)
252  throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
253  "AstrometryTransform::write, something went wrong for file " + fileName);
254 }
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::AstrometryTransformPolynomial, and lsst::jointcal::AstrometryTransformIdentity.

Definition at line 256 of file AstrometryTransform.cc.

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

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.

◆ dec0

double lsst::jointcal::BaseTanWcs::dec0
protectedinherited

Definition at line 617 of file AstrometryTransform.h.

◆ 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.

◆ sin0

double lsst::jointcal::BaseTanWcs::sin0
protectedinherited

Definition at line 618 of file AstrometryTransform.h.


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