LSST Applications g05325017a3+be04ab5bc8,g2079a07aa2+3e9fd84d81,g20cdd03214+600de29eb2,g2305ad1205+5e1b0c67c7,g2bbee38e9b+69e041a28d,g337abbeb29+69e041a28d,g3a166c0a6a+69e041a28d,g4322eb9e3a+9670871960,g4e26c3aaff+9125fe5495,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+608c162c81,g5c3423f6d4+600de29eb2,g607f77f49a+600de29eb2,g69c4c69b53+69e041a28d,g858d7b2824+600de29eb2,g99cad8db69+a7f8c65a2d,g9ddcbc5298+9a081db1e4,ga1e77700b3+2cbb763275,ga8d6e2f9c9+7a336fa28a,gae0086650b+585e252eca,gb0e22166c9+0e73c8378f,gb7c702c636+aec2f3b930,gbb8dafda3b+7d1409958d,gc120e1dc64+acd995d0d8,gc28159a63d+69e041a28d,gcdd4ae20e8+5ef4c697da,gcf0d15dbbd+5ef4c697da,gd1535ee943+1ef0950dc9,gdaeeff99f8+f9a426f77a,gddc38dedce+585e252eca,ge79ae78c31+69e041a28d,gf905fbdb4b+471d9f2c03,gf907d663e3+01e3c6626b,gfbcc870c63+b310236976,w.2024.25
LSST Data Management Base Package
Loading...
Searching...
No Matches
_box.cc
Go to the documentation of this file.
1/*
2 * This file is part of sphgeom.
3 *
4 * Developed for the LSST Data Management System.
5 * This product includes software developed by the LSST Project
6 * (http://www.lsst.org).
7 * See the COPYRIGHT file at the top-level directory of this distribution
8 * for details of code ownership.
9 *
10 * This software is dual licensed under the GNU General Public License and also
11 * under a 3-clause BSD license. Recipients may choose which of these licenses
12 * to use; please see the files gpl-3.0.txt and/or bsd_license.txt,
13 * respectively. If you choose the GPL option then the following text applies
14 * (but note that there is still no warranty even if you opt for BSD instead):
15 *
16 * This program is free software: you can redistribute it and/or modify
17 * it under the terms of the GNU General Public License as published by
18 * the Free Software Foundation, either version 3 of the License, or
19 * (at your option) any later version.
20 *
21 * This program is distributed in the hope that it will be useful,
22 * but WITHOUT ANY WARRANTY; without even the implied warranty of
23 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
24 * GNU General Public License for more details.
25 *
26 * You should have received a copy of the GNU General Public License
27 * along with this program. If not, see <http://www.gnu.org/licenses/>.
28 */
29#include "pybind11/pybind11.h"
30#include "pybind11/numpy.h"
31
32#include "lsst/sphgeom/python.h"
33
35#include "lsst/sphgeom/Box.h"
36#include "lsst/sphgeom/Box3d.h"
37#include "lsst/sphgeom/Circle.h"
40#include "lsst/sphgeom/LonLat.h"
42#include "lsst/sphgeom/Region.h"
44
47
48namespace py = pybind11;
49using namespace pybind11::literals;
50
51namespace lsst {
52namespace sphgeom {
53
54template <>
55void defineClass(py::class_<Box, std::unique_ptr<Box>, Region> &cls) {
56 cls.attr("TYPE_CODE") = py::int_(Box::TYPE_CODE);
57
58 cls.def_static("fromDegrees", &Box::fromDegrees, "lon1"_a, "lat1"_a,
59 "lon2"_a, "lat2"_a);
60 cls.def_static("fromRadians", &Box::fromRadians, "lon1"_a, "lat1"_a,
61 "lon2"_a, "lat2"_a);
62 cls.def_static("empty", &Box::empty);
63 cls.def_static("full", &Box::full);
64 cls.def_static("halfWidthForCircle", &Box::halfWidthForCircle, "radius"_a,
65 "lat"_a);
66 cls.def_static("allLongitudes", &Box::allLongitudes);
67 cls.def_static("allLatitudes", &Box::allLatitudes);
68
69 cls.def(py::init<>());
70 cls.def(py::init<LonLat const &>(), "point"_a);
71 cls.def(py::init<LonLat const &, LonLat const &>(), "point1"_a, "point2"_a);
72 cls.def(py::init<LonLat const &, Angle, Angle>(), "center"_a, "width"_a,
73 "height"_a);
74 cls.def(py::init<NormalizedAngleInterval const &, AngleInterval const &>(),
75 "lon"_a, "lat"_a);
76 cls.def(py::init<Box const &>(), "box"_a);
77
78 cls.def("__eq__", (bool (Box::*)(Box const &) const) & Box::operator==,
79 py::is_operator());
80 cls.def("__eq__", (bool (Box::*)(LonLat const &) const) & Box::operator==,
81 py::is_operator());
82 cls.def("__ne__", (bool (Box::*)(Box const &) const) & Box::operator!=,
83 py::is_operator());
84 cls.def("__ne__", (bool (Box::*)(LonLat const &) const) & Box::operator!=,
85 py::is_operator());
86 cls.def("__contains__",
87 (bool (Box::*)(LonLat const &) const) & Box::contains,
88 py::is_operator());
89 cls.def("__contains__", (bool (Box::*)(Box const &) const) & Box::contains,
90 py::is_operator());
91 // Rewrap this base class method since there are overloads in this subclass
92 cls.def("__contains__",
93 (bool (Box::*)(UnitVector3d const &) const) & Box::contains,
94 py::is_operator());
95
96 cls.def("getLon", &Box::getLon);
97 cls.def("getLat", &Box::getLat);
98 cls.def("isEmpty", &Box::isEmpty);
99 cls.def("isFull", &Box::isFull);
100 cls.def("getCenter", &Box::getCenter);
101 cls.def("getWidth", &Box::getWidth);
102 cls.def("getHeight", &Box::getHeight);
103 cls.def("contains", (bool (Box::*)(LonLat const &) const) & Box::contains);
104 cls.def("contains", (bool (Box::*)(Box const &) const) & Box::contains);
105 // Rewrap these base class methods since there are overloads in this subclass
106 cls.def("contains",
107 (bool (Box::*)(UnitVector3d const &) const) & Box::contains);
108 cls.def("contains", py::vectorize((bool (Box::*)(double, double, double) const)&Box::contains),
109 "x"_a, "y"_a, "z"_a);
110 cls.def("contains", py::vectorize((bool (Box::*)(double, double) const)&Box::contains),
111 "lon"_a, "lat"_a);
112 cls.def("isDisjointFrom",
113 (bool (Box::*)(LonLat const &) const) & Box::isDisjointFrom);
114 cls.def("isDisjointFrom",
115 (bool (Box::*)(Box const &) const) & Box::isDisjointFrom);
116 cls.def("intersects",
117 (bool (Box::*)(LonLat const &) const) & Box::intersects);
118 cls.def("intersects", (bool (Box::*)(Box const &) const) & Box::intersects);
119 cls.def("isWithin", (bool (Box::*)(LonLat const &) const) & Box::isWithin);
120 cls.def("isWithin", (bool (Box::*)(Box const &) const) & Box::isWithin);
121 cls.def("clipTo", (Box & (Box::*)(LonLat const &)) & Box::clipTo);
122 cls.def("clipTo", (Box & (Box::*)(Box const &)) & Box::clipTo);
123 cls.def("clippedTo", (Box(Box::*)(LonLat const &) const) & Box::clippedTo);
124 cls.def("clippedTo", (Box(Box::*)(Box const &) const) & Box::clippedTo);
125 cls.def("expandTo", (Box & (Box::*)(LonLat const &)) & Box::expandTo);
126 cls.def("expandTo", (Box & (Box::*)(Box const &)) & Box::expandTo);
127 cls.def("expandedTo",
128 (Box(Box::*)(LonLat const &) const) & Box::expandedTo);
129 cls.def("expandedTo", (Box(Box::*)(Box const &) const) & Box::expandedTo);
130 cls.def("dilateBy", (Box & (Box::*)(Angle)) & Box::dilateBy, "angle"_a);
131 cls.def("dilateBy", (Box & (Box::*)(Angle, Angle)) & Box::dilateBy,
132 "width"_a, "height"_a);
133 cls.def("dilatedBy", (Box(Box::*)(Angle) const) & Box::dilatedBy,
134 "angle"_a);
135 cls.def("dilatedBy", (Box(Box::*)(Angle, Angle) const) & Box::dilatedBy,
136 "width"_a, "height"_a);
137 cls.def("erodeBy", (Box & (Box::*)(Angle)) & Box::erodeBy, "angle"_a);
138 cls.def("erodeBy", (Box & (Box::*)(Angle, Angle)) & Box::erodeBy, "width"_a,
139 "height"_a);
140 cls.def("erodedBy", (Box(Box::*)(Angle) const) & Box::erodedBy, "angle"_a);
141 cls.def("erodedBy", (Box(Box::*)(Angle, Angle) const) & Box::erodedBy,
142 "width"_a, "height"_a);
143 cls.def("getArea", &Box::getArea);
144 cls.def("relate",
145 (Relationship(Box::*)(LonLat const &) const) & Box::relate,
146 "point"_a);
147 // Rewrap this base class method since there are overloads in this subclass
148 cls.def("relate",
149 (Relationship(Box::*)(Region const &) const) & Box::relate,
150 "region"_a);
151
152 // Note that the Region interface has already been wrapped.
153
154 cls.def("__str__", [](Box const &self) {
155 return py::str("Box({!s}, {!s})").format(self.getLon(), self.getLat());
156 });
157 cls.def("__repr__", [](Box const &self) {
158 return py::str("Box({!r}, {!r})").format(self.getLon(), self.getLat());
159 });
160 cls.def(py::pickle(&python::encode, &python::decode<Box>));
161}
162
163} // sphgeom
164} // lsst
This file defines a class for representing angle intervals.
This file declares a class for representing axis-aligned bounding boxes in ℝ³.
This file declares a class for representing circular regions on the unit sphere.
This file declares a class for representing convex polygons with great circle edges on the unit spher...
This file contains a class representing spherical coordinates.
This file declares a class representing closed intervals of normalized angles, i.e.
This file defines an interface for spherical regions.
This file declares a class for representing unit vectors in ℝ³.
Angle represents an angle in radians.
Definition Angle.h:50
Box represents a rectangle in spherical coordinate space that contains its boundary.
Definition Box.h:61
Box clippedTo(LonLat const &x) const
clippedTo returns the intersection of this box and x.
Definition Box.h:242
AngleInterval const & getLat() const
getLat returns the latitude interval of this box.
Definition Box.h:155
Box dilatedBy(Angle r) const
Definition Box.h:280
static NormalizedAngle halfWidthForCircle(Angle r, Angle lat)
halfWidthForCircle computes the half-width of bounding boxes for circles with radius r and centers at...
Definition Box.cc:50
Box & erodeBy(Angle r)
Definition Box.h:299
NormalizedAngle getWidth() const
getWidth returns the width in longitude angle of this box.
Definition Box.h:172
double getArea() const
getArea returns the area of this box in steradians.
Definition Box.cc:122
static Box fromRadians(double lon1, double lat1, double lon2, double lat2)
Definition Box.h:71
Box expandedTo(LonLat const &x) const
Definition Box.h:270
NormalizedAngleInterval const & getLon() const
getLon returns the longitude interval of this box.
Definition Box.h:152
Relationship relate(LonLat const &p) const
Definition Box.h:304
bool isDisjointFrom(LonLat const &x) const
Definition Box.h:193
static Box empty()
Definition Box.h:76
bool isFull() const
isFull returns true if this box contains all points on the unit sphere.
Definition Box.h:162
static constexpr uint8_t TYPE_CODE
Definition Box.h:63
Box & dilateBy(Angle r)
dilateBy minimally expands this Box to include all points within angular separation r of its boundary...
Definition Box.cc:85
Angle getHeight() const
getHeight returns the height in latitude angle of this box.
Definition Box.h:176
LonLat getCenter() const
getCenter returns the center of this box.
Definition Box.h:166
static NormalizedAngleInterval allLongitudes()
allLongitudes returns a normalized angle interval containing all valid longitude angles.
Definition Box.h:88
bool isWithin(LonLat const &x) const
Definition Box.h:213
Box erodedBy(Angle r) const
Definition Box.h:301
Box & expandTo(LonLat const &x)
Definition Box.h:253
static AngleInterval allLatitudes()
allLatitudes returns an angle interval containing all valid latitude angles.
Definition Box.h:94
bool contains(LonLat const &x) const
Definition Box.h:181
static Box fromDegrees(double lon1, double lat1, double lon2, double lat2)
Definition Box.h:66
bool isEmpty() const
isEmpty returns true if this box does not contain any points.
Definition Box.h:158
Box & clipTo(LonLat const &x)
clipTo shrinks this box until it contains only x.
Definition Box.h:224
static Box full()
Definition Box.h:78
bool intersects(LonLat const &x) const
Definition Box.h:201
LonLat represents a spherical coordinate (longitude/latitude angle) pair.
Definition LonLat.h:55
Region is a minimal interface for 2-dimensional regions on the unit sphere.
Definition Region.h:86
UnitVector3d is a unit vector in ℝ³ with components stored in double precision.
pybind11::bytes encode(Region const &self)
Encode a Region as a pybind11 bytes object.
Definition utils.h:60
void defineClass(Pybind11Class &cls)
std::bitset< 3 > Relationship
Relationship describes how two sets are related.
This file declares a class for representing longitude/latitude angle boxes on the unit sphere.
This file declares a class for representing elliptical regions on the unit sphere.