LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
LSSTDataManagementBasePackage
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"
26 #include "lsst/geom/LinearTransform.h"
27 #include "lsst/afw/geom/SkyWcs.h"
28 #include "lsst/meas/astrom/SipTransform.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
SkyWcs.h
SipTransform.h
lsst::meas::astrom::PYBIND11_MODULE
PYBIND11_MODULE(sipTransform, mod)
Definition: sipTransform.cc:104
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
lsst::meas::astrom::makeWcs
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.
Definition: SipTransform.cc:148
lsst::meas::astrom::transformWcsPixels
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...
Definition: SipTransform.cc:173
lsst::meas::astrom::rotateWcsPixelsBy90
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&#39;s cente...
Definition: SipTransform.cc:179
lsst::geom::LinearTransform
A 2D linear coordinate transformation.
Definition: LinearTransform.h:69
Generated on Sat Nov 16 2019 09:14:48 for LSSTApplications by   doxygen 1.8.13