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
sipTransform.cc
Go to the documentation of this file.
1 /*
2  * LSST Data Management System
3  *
4  * This product includes software developed by the
5  * LSST Project (http://www.lsst.org/).
6  * See the COPYRIGHT file
7  *
8  * This program is free software: you can redistribute it and/or modify
9  * it under the terms of the GNU General Public License as published by
10  * the Free Software Foundation, either version 3 of the License, or
11  * (at your option) any later version.
12  *
13  * This program is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16  * GNU General Public License for more details.
17  *
18  * You should have received a copy of the LSST License Statement and
19  * the GNU General Public License along with this program. If not,
20  * see <https://www.lsstcorp.org/LegalNotices/>.
21  */
22 #include "pybind11/pybind11.h"
23 #include "pybind11/stl.h"
24 
25 #include "lsst/geom/Point.h"
27 #include "lsst/afw/geom/SkyWcs.h"
29 
30 namespace py = pybind11;
31 using namespace pybind11::literals;
32 
33 namespace lsst {
34 namespace meas {
35 namespace astrom {
36 namespace {
37 
38 static void declareSipTransformBase(py::module &mod) {
39  py::class_<SipTransformBase, std::shared_ptr<SipTransformBase>> cls(mod, "_SipTransformBase");
40 
41  cls.def("getPixelOrigin", &SipTransformBase::getPixelOrigin, py::return_value_policy::copy);
42  cls.def("getCdMatrix", &SipTransformBase::getCdMatrix, py::return_value_policy::copy);
43  cls.def("getPoly", &SipTransformBase::getPoly, py::return_value_policy::copy);
44 }
45 
46 static void declareSipForwardTransform(py::module &mod) {
47  py::class_<SipForwardTransform, std::shared_ptr<SipForwardTransform>, SipTransformBase> cls(
48  mod, "SipForwardTransform");
49 
50  cls.def(py::init<geom::Point2D const &, geom::LinearTransform const &, PolynomialTransform const &>(),
51  "pixelOrigin"_a, "cdMatrix"_a, "forwardSipPoly"_a);
52  cls.def(py::init<SipForwardTransform const &>(), "other"_a);
53 
54  cls.def_static("convert",
55  (SipForwardTransform(*)(PolynomialTransform const &, geom::Point2D const &,
56  geom::LinearTransform const &)) &
57  SipForwardTransform::convert,
58  "poly"_a, "pixelOrigin"_a, "cdMatrix"_a);
59  cls.def_static("convert",
60  (SipForwardTransform(*)(ScaledPolynomialTransform const &, geom::Point2D const &,
61  geom::LinearTransform const &)) &
62  SipForwardTransform::convert,
63  "scaled"_a, "pixelOrigin"_a, "cdMatrix"_a);
64  cls.def_static("convert",
65  (SipForwardTransform(*)(ScaledPolynomialTransform const &)) & SipForwardTransform::convert,
66  "scaled"_a);
67 
68  cls.def("__call__", &SipForwardTransform::operator(), "in"_a);
69  cls.def("transformPixels", &SipForwardTransform::transformPixels, "s"_a);
70 
71  cls.def("linearize", &SipForwardTransform::linearize);
72 }
73 
74 static void declareSipReverseTransform(py::module &mod) {
75  py::class_<SipReverseTransform, std::shared_ptr<SipReverseTransform>, SipTransformBase> cls(
76  mod, "SipReverseTransform");
77 
78  cls.def(py::init<geom::Point2D const &, geom::LinearTransform const &, PolynomialTransform const &>(),
79  "pixelOrigin"_a, "cdMatrix"_a, "reverseSipPoly"_a);
80  cls.def(py::init<SipReverseTransform const &>(), "other"_a);
81 
82  cls.def_static("convert",
83  (SipReverseTransform(*)(PolynomialTransform const &, geom::Point2D const &,
84  geom::LinearTransform const &)) &
85  SipReverseTransform::convert,
86  "poly"_a, "pixelOrigin"_a, "cdMatrix"_a);
87  cls.def_static("convert",
88  (SipReverseTransform(*)(ScaledPolynomialTransform const &, geom::Point2D const &,
89  geom::LinearTransform const &)) &
90  SipReverseTransform::convert,
91  "scaled"_a, "pixelOrigin"_a, "cdMatrix"_a);
92  cls.def_static("convert",
93  (SipReverseTransform(*)(ScaledPolynomialTransform const &)) & SipReverseTransform::convert,
94  "scaled"_a);
95 
96  cls.def("__call__", &SipReverseTransform::operator(), "in"_a);
97  cls.def("transformPixels", &SipReverseTransform::transformPixels, "s"_a);
98 
99  cls.def("linearize", &SipReverseTransform::linearize);
100 }
101 
102 } // namespace
103 
104 PYBIND11_MODULE(sipTransform, mod) {
105  declareSipTransformBase(mod);
106  declareSipForwardTransform(mod);
107  declareSipReverseTransform(mod);
108 
109  mod.def("makeWcs", makeWcs, "sipForward"_a, "sipReverse"_a, "skyOrigin"_a);
110  mod.def("transformWcsPixels", transformWcsPixels, "wcs"_a, "s"_a);
111  mod.def("rotateWcsPixelsBy90", rotateWcsPixelsBy90, "wcs"_a, "nQuarter"_a, "dimensions"_a);
112 }
113 
114 } // namespace astrom
115 } // namespace meas
116 } // namespace lsst
A 2D linear coordinate transformation.
std::shared_ptr< afw::geom::SkyWcs > transformWcsPixels(afw::geom::SkyWcs const &wcs, geom::AffineTransform const &s)
Create a new SkyWcs whose pixel coordinate system has been transformed via an affine transform.
PYBIND11_MODULE(sipTransform, mod)
std::shared_ptr< afw::geom::SkyWcs > rotateWcsPixelsBy90(afw::geom::SkyWcs const &wcs, int nQuarter, geom::Extent2I const &dimensions)
Return a new SkyWcs that represents a rotation of the image it corresponds to about the image's cente...
std::shared_ptr< afw::geom::SkyWcs > makeWcs(SipForwardTransform const &sipForward, SipReverseTransform const &sipReverse, geom::SpherePoint const &skyOrigin)
Create a new TAN SIP Wcs from a pair of SIP transforms and the sky origin.
A base class for image defects.