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 Types | Public Member Functions | Public Attributes | List of all members
lsst::afw::geom::ellipses::BaseCore::Transformer Class Referencefinal

A temporary-only expression object for ellipse core transformations. More...

#include <Transformer.h>

Public Types

using DerivativeMatrix = Eigen::Matrix3d
 Matrix type for derivative with respect to input ellipse parameters. More...
 
using TransformDerivativeMatrix = Eigen::Matrix< double, 3, 4 >
 Matrix type for derivative with respect to transform parameters. More...
 

Public Member Functions

 Transformer (BaseCore &input_, lsst::geom::LinearTransform const &transform_)
 Standard constructor. More...
 
std::shared_ptr< BaseCore > copy () const
 Return a new transformed ellipse core. More...
 
void inPlace ()
 Transform the ellipse core in-place. More...
 
void apply (BaseCore &result) const
 
DerivativeMatrix d () const
 Return the derivative of transformed core with respect to input core. More...
 
TransformDerivativeMatrix dTransform () const
 Return the derivative of transformed core with respect to transform parameters. More...
 

Public Attributes

BaseCore & input
 input core to be transformed More...
 
lsst::geom::LinearTransform const & transform
 transform object More...
 

Detailed Description

A temporary-only expression object for ellipse core transformations.

Transformer simply provides a clean syntax for transform-related operations, including in-place and new-object transformations, derivatives of the transformations, and implicit conversion to a shared_ptr to a new transformed core.

Definition at line 49 of file Transformer.h.

Member Typedef Documentation

◆ DerivativeMatrix

using lsst::afw::geom::ellipses::BaseCore::Transformer::DerivativeMatrix = Eigen::Matrix3d

Matrix type for derivative with respect to input ellipse parameters.

Definition at line 52 of file Transformer.h.

◆ TransformDerivativeMatrix

using lsst::afw::geom::ellipses::BaseCore::Transformer::TransformDerivativeMatrix = Eigen::Matrix<double, 3, 4>

Matrix type for derivative with respect to transform parameters.

Definition at line 55 of file Transformer.h.

Constructor & Destructor Documentation

◆ Transformer()

lsst::afw::geom::ellipses::BaseCore::Transformer::Transformer ( BaseCore &  input_,
lsst::geom::LinearTransform const &  transform_ 
)
inline

Standard constructor.

Definition at line 58 of file Transformer.h.

59  : input(input_), transform(transform_) {}
lsst::afw::geom::ellipses::BaseCore::Transformer::transform
lsst::geom::LinearTransform const & transform
transform object
Definition: Transformer.h:76
lsst::afw::geom::ellipses::BaseCore::Transformer::input
BaseCore & input
input core to be transformed
Definition: Transformer.h:75

Member Function Documentation

◆ apply()

void lsst::afw::geom::ellipses::BaseCore::Transformer::apply ( BaseCore &  result) const

Definition at line 41 of file Transformer.cc.

41  {
42  Eigen::Matrix2d m;
43  input._assignToQuadrupole(m(0, 0), m(1, 1), m(0, 1));
44  m(1, 0) = m(0, 1);
45  m = transform.getMatrix() * m * transform.getMatrix().transpose();
46  result._assignFromQuadrupole(m(0, 0), m(1, 1), m(0, 1));
47 }
result
py::object result
Definition: _schema.cc:429
m
int m
Definition: SpanSet.cc:48
lsst::afw::geom::ellipses::BaseCore::_assignToQuadrupole
virtual void _assignToQuadrupole(double &ixx, double &iyy, double &ixy) const =0
Return the size of the bounding box for the ellipse core.
lsst::geom::LinearTransform::getMatrix
Matrix const & getMatrix() const noexcept
Definition: LinearTransform.h:151

◆ copy()

std::shared_ptr< BaseCore > lsst::afw::geom::ellipses::BaseCore::Transformer::copy ( ) const

Return a new transformed ellipse core.

Definition at line 33 of file Transformer.cc.

33  {
34  std::shared_ptr<BaseCore> r(input.clone());
35  apply(*r);
36  return r;
37 }
lsst::afw::geom::ellipses::BaseCore::clone
std::shared_ptr< BaseCore > clone() const
Deep-copy the Core.
Definition: BaseCore.h:82

◆ d()

BaseCore::Transformer::DerivativeMatrix lsst::afw::geom::ellipses::BaseCore::Transformer::d ( ) const

Return the derivative of transformed core with respect to input core.

Definition at line 49 of file Transformer.cc.

49  {
50  std::shared_ptr<BaseCore> output(input.clone());
51  Eigen::Matrix2d m;
52  Jacobian rhs = input._dAssignToQuadrupole(m(0, 0), m(1, 1), m(0, 1));
53  m(1, 0) = m(0, 1);
54  m = transform.getMatrix() * m * transform.getMatrix().transpose();
55  Jacobian lhs = output->_dAssignFromQuadrupole(m(0, 0), m(1, 1), m(0, 1));
56  Jacobian mid = Jacobian::Zero();
57  mid(0, 0) = transform[lsst::geom::LinearTransform::XX] * transform[lsst::geom::LinearTransform::XX];
58  mid(0, 1) = transform[lsst::geom::LinearTransform::XY] * transform[lsst::geom::LinearTransform::XY];
59  mid(0, 2) = 2 * transform[lsst::geom::LinearTransform::XY] * transform[lsst::geom::LinearTransform::XX];
60  mid(1, 0) = transform[lsst::geom::LinearTransform::YX] * transform[lsst::geom::LinearTransform::YX];
61  mid(1, 1) = transform[lsst::geom::LinearTransform::YY] * transform[lsst::geom::LinearTransform::YY];
62  mid(1, 2) = 2 * transform[lsst::geom::LinearTransform::YY] * transform[lsst::geom::LinearTransform::YX];
63  mid(2, 0) = transform[lsst::geom::LinearTransform::YX] * transform[lsst::geom::LinearTransform::XX];
64  mid(2, 1) = transform[lsst::geom::LinearTransform::YY] * transform[lsst::geom::LinearTransform::XY];
65  mid(2, 2) = transform[lsst::geom::LinearTransform::XX] * transform[lsst::geom::LinearTransform::YY] +
66  transform[lsst::geom::LinearTransform::XY] * transform[lsst::geom::LinearTransform::YX];
67  return lhs * mid * rhs;
68 }
lsst::afw::geom::ellipses::BaseCore::_dAssignToQuadrupole
virtual Jacobian _dAssignToQuadrupole(double &ixx, double &iyy, double &ixy) const =0
Return the size of the bounding box for the ellipse core.
lsst::afw::geom::ellipses::BaseCore::Jacobian
Eigen::Matrix3d Jacobian
Parameter Jacobian matrix type.
Definition: BaseCore.h:64

◆ dTransform()

BaseCore::Transformer::TransformDerivativeMatrix lsst::afw::geom::ellipses::BaseCore::Transformer::dTransform ( ) const

Return the derivative of transformed core with respect to transform parameters.

Definition at line 70 of file Transformer.cc.

70  {
71  std::shared_ptr<BaseCore> output(input.clone());
72  Eigen::Matrix2d m;
73  input._assignToQuadrupole(m(0, 0), m(1, 1), m(0, 1));
74  Eigen::Matrix<double, 3, 4> mid = Eigen::Matrix<double, 3, 4>::Zero();
75  m(1, 0) = m(0, 1);
76  mid(0, lsst::geom::LinearTransform::XX) = 2.0 * (transform[lsst::geom::LinearTransform::XX] * m(0, 0) +
77  transform[lsst::geom::LinearTransform::XY] * m(0, 1));
78  mid(0, lsst::geom::LinearTransform::XY) = 2.0 * (transform[lsst::geom::LinearTransform::XX] * m(0, 1) +
79  transform[lsst::geom::LinearTransform::XY] * m(1, 1));
80  mid(1, lsst::geom::LinearTransform::YX) = 2.0 * (transform[lsst::geom::LinearTransform::YX] * m(0, 0) +
81  transform[lsst::geom::LinearTransform::YY] * m(0, 1));
82  mid(1, lsst::geom::LinearTransform::YY) = 2.0 * (transform[lsst::geom::LinearTransform::YX] * m(0, 1) +
83  transform[lsst::geom::LinearTransform::YY] * m(1, 1));
84  mid(2, lsst::geom::LinearTransform::XX) = transform[lsst::geom::LinearTransform::YX] * m(0, 0) +
85  transform[lsst::geom::LinearTransform::YY] * m(0, 1);
86  mid(2, lsst::geom::LinearTransform::XY) = transform[lsst::geom::LinearTransform::YX] * m(0, 1) +
87  transform[lsst::geom::LinearTransform::YY] * m(1, 1);
88  mid(2, lsst::geom::LinearTransform::YX) = transform[lsst::geom::LinearTransform::XX] * m(0, 0) +
89  transform[lsst::geom::LinearTransform::XY] * m(0, 1);
90  mid(2, lsst::geom::LinearTransform::YY) = transform[lsst::geom::LinearTransform::XX] * m(0, 1) +
91  transform[lsst::geom::LinearTransform::XY] * m(1, 1);
92  m = transform.getMatrix() * m * transform.getMatrix().transpose();
93  Jacobian lhs = output->_dAssignFromQuadrupole(m(0, 0), m(1, 1), m(0, 1));
94  return lhs * mid;
95 }

◆ inPlace()

void lsst::afw::geom::ellipses::BaseCore::Transformer::inPlace ( )

Transform the ellipse core in-place.

Definition at line 39 of file Transformer.cc.

39 { apply(input); }

Member Data Documentation

◆ input

BaseCore& lsst::afw::geom::ellipses::BaseCore::Transformer::input

input core to be transformed

Definition at line 75 of file Transformer.h.

◆ transform

lsst::geom::LinearTransform const& lsst::afw::geom::ellipses::BaseCore::Transformer::transform

transform object

Definition at line 76 of file Transformer.h.

The documentation for this class was generated from the following files:
Generated on Fri Nov 26 2021 10:16:11 for LSST Applications by doxygen 1.9.1