LSSTApplications  20.0.0
LSSTDataManagementBasePackage
Transformer.cc
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 
3 /*
4  * LSST Data Management System
5  * Copyright 2008, 2009, 2010 LSST Corporation.
6  *
7  * This product includes software developed by the
8  * LSST Project (http://www.lsst.org/).
9  *
10  * This program is free software: you can redistribute it and/or modify
11  * it under the terms of the GNU General Public License as published by
12  * the Free Software Foundation, either version 3 of the License, or
13  * (at your option) any later version.
14  *
15  * This program is distributed in the hope that it will be useful,
16  * but WITHOUT ANY WARRANTY; without even the implied warranty of
17  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18  * GNU General Public License for more details.
19  *
20  * You should have received a copy of the LSST License Statement and
21  * the GNU General Public License along with this program. If not,
22  * see <http://www.lsstcorp.org/LegalNotices/>.
23  */
25 
26 #include "Eigen/LU"
27 
28 namespace lsst {
29 namespace afw {
30 namespace geom {
31 namespace ellipses {
32 
35  apply(*r);
36  return r;
37 }
38 
39 void BaseCore::Transformer::inPlace() { apply(input); }
40 
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 }
48 
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();
67  return lhs * mid * rhs;
68 }
69 
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);
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 }
96 
98  std::shared_ptr<Ellipse> r = std::make_shared<Ellipse>(
99  input.getCore().transform(transform.getLinear()).copy(), transform(input.getCenter()));
100  return r;
101 }
102 
104  input.setCenter(transform(input.getCenter()));
105  input.getCore().transform(transform.getLinear()).inPlace();
106 }
107 
109  DerivativeMatrix r = DerivativeMatrix::Zero();
110  r.block<2, 2>(3, 3) = transform.getLinear().getMatrix();
111  r.block<3, 3>(0, 0) = input.getCore().transform(transform.getLinear()).d();
112  return r;
113 }
114 
116  TransformDerivativeMatrix r = TransformDerivativeMatrix::Zero();
117  r.block<2, 6>(3, 0) = transform.dTransform(input.getCenter());
118  r.block<3, 4>(0, 0) = input.getCore().transform(transform.getLinear()).dTransform();
119  return r;
120 }
121 } // namespace ellipses
122 } // namespace geom
123 } // namespace afw
124 } // namespace lsst
lsst::geom::LinearTransform::XX
@ XX
Definition: LinearTransform.h:71
lsst::afw::geom::ellipses::BaseCore::Transformer::transform
lsst::geom::LinearTransform const & transform
transform object
Definition: Transformer.h:76
std::shared_ptr
STL class.
ellipses
lsst::afw::geom::ellipses::Ellipse::Transformer::d
DerivativeMatrix d() const
Return the derivative of transform output ellipse with respect to input ellipse.
Definition: Transformer.cc:108
lsst::afw::geom::ellipses::Ellipse::Transformer::inPlace
void inPlace()
Transform the ellipse in-place.
Definition: Transformer.cc:103
lsst::afw
Definition: imageAlgorithm.dox:1
lsst::afw::geom::ellipses::Ellipse::Transformer::dTransform
TransformDerivativeMatrix dTransform() const
Return the derivative of transform output ellipse with respect to transform parameters.
Definition: Transformer.cc:115
lsst::afw::geom::ellipses::BaseCore::Transformer::input
BaseCore & input
input core to be transformed
Definition: Transformer.h:75
lsst::geom::LinearTransform::XY
@ XY
Definition: LinearTransform.h:71
lsst::afw::geom::ellipses::BaseCore::Transformer::apply
void apply(BaseCore &result) const
Definition: Transformer.cc:41
lsst::geom::LinearTransform::YX
@ YX
Definition: LinearTransform.h:71
lsst::afw::geom::ellipses::BaseCore
A base class for parametrizations of the "core" of an ellipse - the ellipticity and size.
Definition: BaseCore.h:55
lsst::geom::LinearTransform::YY
@ YY
Definition: LinearTransform.h:71
lsst::afw::geom::ellipses::BaseCore::Transformer::inPlace
void inPlace()
Transform the ellipse core in-place.
Definition: Transformer.cc:39
Transformer.h
lsst::afw::geom::ellipses::BaseCore::Transformer::TransformDerivativeMatrix
Eigen::Matrix< double, 3, 4 > TransformDerivativeMatrix
Matrix type for derivative with respect to transform parameters.
Definition: Transformer.h:55
lsst::afw::geom::ellipses::BaseCore::Transformer::copy
std::shared_ptr< BaseCore > copy() const
Return a new transformed ellipse core.
Definition: Transformer.cc:33
result
py::object result
Definition: _schema.cc:429
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
lsst::afw::geom::ellipses::Ellipse::Transformer::TransformDerivativeMatrix
Eigen::Matrix< double, 5, 6 > TransformDerivativeMatrix
Matrix type for derivative with respect to transform parameters.
Definition: Transformer.h:92
lsst::geom
Definition: geomOperators.dox:4
lsst::afw::geom::ellipses::BaseCore::Jacobian
Eigen::Matrix3d Jacobian
Parameter Jacobian matrix type.
Definition: BaseCore.h:64
lsst::afw::geom::ellipses::BaseCore::transform
Transformer transform(lsst::geom::LinearTransform const &transform)
Return the transform that maps the ellipse to the unit circle.
Definition: Transformer.h:116
lsst::afw::geom::ellipses::BaseCore::Transformer::dTransform
TransformDerivativeMatrix dTransform() const
Return the derivative of transformed core with respect to transform parameters.
Definition: Transformer.cc:70
lsst::afw::geom::ellipses::BaseCore::clone
std::shared_ptr< BaseCore > clone() const
Deep-copy the Core.
Definition: BaseCore.h:82
lsst::afw::geom::ellipses::Ellipse::Transformer::copy
std::shared_ptr< Ellipse > copy() const
Return a new transformed ellipse.
Definition: Transformer.cc:97
lsst::afw::geom::ellipses::BaseCore::Transformer::DerivativeMatrix
Eigen::Matrix3d DerivativeMatrix
Matrix type for derivative with respect to input ellipse parameters.
Definition: Transformer.h:52
lsst::afw::geom::ellipses::BaseCore::Transformer::d
DerivativeMatrix d() const
Return the derivative of transformed core with respect to input core.
Definition: Transformer.cc:49
lsst::afw::geom::ellipses::Ellipse::Transformer::DerivativeMatrix
Eigen::Matrix< double, 5, 5 > DerivativeMatrix
Matrix type for derivative with respect to input ellipse parameters.
Definition: Transformer.h:89
m
int m
Definition: SpanSet.cc:49