LSSTApplications  19.0.0+10,19.0.0+80,19.0.0-1-g20d9b18+31,19.0.0-1-g49a97f9+4,19.0.0-1-g8c57eb9+31,19.0.0-1-g9a028c0+13,19.0.0-1-ga72da6b+3,19.0.0-1-gb77924a+15,19.0.0-1-gbfe0924+66,19.0.0-1-gc3516c3,19.0.0-1-gefe1d0d+49,19.0.0-1-gf8cb8b4+3,19.0.0-14-g7511ce4+6,19.0.0-16-g3dc1a33c+6,19.0.0-17-g59f0e8a+4,19.0.0-17-g9c22e3c+9,19.0.0-18-g35bb99870+2,19.0.0-19-g2772d4a+9,19.0.0-2-g260436e+53,19.0.0-2-g31cdcee,19.0.0-2-g9675b69+3,19.0.0-2-gaa2795f,19.0.0-2-gbcc4de1,19.0.0-2-gd6f004e+6,19.0.0-2-gde8e5e3+5,19.0.0-2-gff6972b+15,19.0.0-21-ge306cd8,19.0.0-21-gf122e69+2,19.0.0-3-g2d43a51+2,19.0.0-3-gf3b1435+6,19.0.0-4-g56feb96,19.0.0-4-gac56cce+17,19.0.0-49-gce872c1+1,19.0.0-5-g66946eb+6,19.0.0-5-gd8897ba+6,19.0.0-51-gfc4a647e,19.0.0-7-g686a884+5,19.0.0-7-g886f805+5,19.0.0-8-g305ff64,w.2020.17
LSSTDataManagementBasePackage
ellipse.cc
Go to the documentation of this file.
1 /*
2  * LSST Data Management System
3  * See COPYRIGHT file at the top of the source tree.
4  *
5  * This product includes software developed by the
6  * LSST Project (http://www.lsst.org/).
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 
24 #include <memory>
25 
26 #include "lsst/sphgeom/Box.h"
27 #include "lsst/sphgeom/Box3d.h"
28 #include "lsst/sphgeom/Circle.h"
30 #include "lsst/sphgeom/Ellipse.h"
31 #include "lsst/sphgeom/Matrix3d.h"
32 #include "lsst/sphgeom/Region.h"
34 
37 
38 namespace py = pybind11;
39 using namespace pybind11::literals;
40 
41 namespace lsst {
42 namespace sphgeom {
43 namespace {
44 
46  uint8_t const *buffer = reinterpret_cast<uint8_t const *>(
47  PYBIND11_BYTES_AS_STRING(bytes.ptr()));
48  size_t n = static_cast<size_t>(PYBIND11_BYTES_SIZE(bytes.ptr()));
49  return Ellipse::decode(buffer, n);
50 }
51 
53  py::module::import("lsst.sphgeom.region");
54 
55  py::class_<Ellipse, std::unique_ptr<Ellipse>, Region> cls(mod, "Ellipse");
56 
57  cls.attr("TYPE_CODE") = py::int_(Ellipse::TYPE_CODE);
58 
59  cls.def_static("empty", &Ellipse::empty);
60  cls.def_static("full", &Ellipse::full);
61 
62  cls.def(py::init<>());
63  cls.def(py::init<Circle const &>(), "circle"_a);
64  cls.def(py::init<UnitVector3d const &, Angle>(), "center"_a,
65  "angle"_a = Angle(0.0));
66  cls.def(py::init<UnitVector3d const &, UnitVector3d const &, Angle>(),
67  "focus1"_a, "focus2"_a, "alpha"_a);
68  cls.def(py::init<UnitVector3d const &, Angle, Angle, Angle>(), "center"_a,
69  "alpha"_a, "beta"_a, "orientation"_a);
70  cls.def(py::init<Ellipse const &>(), "ellipse"_a);
71 
72  cls.def("__eq__", &Ellipse::operator==, py::is_operator());
73  cls.def("__ne__", &Ellipse::operator!=, py::is_operator());
74 
75  cls.def("isEmpty", &Ellipse::isEmpty);
76  cls.def("isFull", &Ellipse::isFull);
77  cls.def("isGreatCircle", &Ellipse::isGreatCircle);
78  cls.def("isCircle", &Ellipse::isCircle);
79  cls.def("getTransformMatrix", &Ellipse::getTransformMatrix);
80  cls.def("getCenter", &Ellipse::getCenter);
81  cls.def("getF1", &Ellipse::getF1);
82  cls.def("getF2", &Ellipse::getF2);
83  cls.def("getAlpha", &Ellipse::getAlpha);
84  cls.def("getBeta", &Ellipse::getBeta);
85  cls.def("getGamma", &Ellipse::getGamma);
86  cls.def("complement", &Ellipse::complement);
87  cls.def("complemented", &Ellipse::complemented);
88 
89  // Note that the Region interface has already been wrapped.
90 
91  // The lambda is necessary for now; returning the unique pointer
92  // directly leads to incorrect results and crashes.
93  cls.def_static("decode",
94  [](py::bytes bytes) { return decode(bytes).release(); },
95  "bytes"_a);
96 
97  cls.def("__str__", [](Ellipse const &self) {
98  return py::str("Ellipse({!s}, {!s}, {!s})")
99  .format(self.getF1(), self.getF2(), self.getAlpha());
100  });
101  cls.def("__repr__", [](Ellipse const &self) {
102  return py::str("Ellipse({!r}, {!r}, {!r})")
103  .format(self.getF1(), self.getF2(), self.getAlpha());
104  });
105  cls.def(py::pickle(
106  [](const Ellipse &self) { return python::encode(self); },
107  [](py::bytes bytes) { return decode(bytes).release(); }));
108 }
109 
110 } // <anonymous>
111 } // sphgeom
112 } // lsst
This file declares a class for representing circular regions on the unit sphere.
This file declares a class for representing longitude/latitude angle boxes on the unit sphere...
PYBIND11_MODULE(ellipse, mod)
Definition: ellipse.cc:42
This file declares a class for representing axis-aligned bounding boxes in ℝ³.
This file declares a class for representing unit vectors in ℝ³.
A base class for image defects.
This file defines an interface for spherical regions.
STL class.
table::Key< table::Array< std::uint8_t > > bytes
Definition: Transform.cc:199
This file declares a class for representing elliptical regions on the unit sphere.
lsst::geom::Angle Angle
Definition: misc.h:33
This file declares a class for representing convex polygons with great circle edges on the unit spher...
This file contains a class representing 3x3 real matrices.