LSST Applications g0aad566f14+98bdd3b178,g0dd44d6229+e0c2448475,g2079a07aa2+86d27d4dc4,g2305ad1205+28192cebb5,g2bbee38e9b+047b288a59,g337abbeb29+047b288a59,g33d1c0ed96+047b288a59,g3a166c0a6a+047b288a59,g3d1719c13e+98bdd3b178,g487adcacf7+08349b74a1,g4be5004598+98bdd3b178,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+742623b69a,g63cd9335cc+585e252eca,g858d7b2824+98bdd3b178,g88963caddf+0cb8e002cc,g88c60e1ba4+10d1129650,g99cad8db69+26979ce5a5,g9ddcbc5298+9a081db1e4,ga1e77700b3+a912195c07,gae0086650b+585e252eca,gb0e22166c9+60f28cb32d,gb1fda7adca+90d62c6f66,gb3a676b8dc+b4feba26a1,gb4b16eec92+950562e714,gba4ed39666+c2a2e4ac27,gbb8dafda3b+2fe9a6a885,gc120e1dc64+0553034b2f,gc28159a63d+047b288a59,gc3e9b769f7+f8a3f20ab6,gcf0d15dbbd+e0c2448475,gdaeeff99f8+f9a426f77a,gdb0af172c8+f61b8374e1,ge79ae78c31+047b288a59,w.2024.20
LSST Data Management Base Package
Loading...
Searching...
No Matches
Frame.cc
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2/*
3 * This file is part of jointcal.
4 *
5 * Developed for the LSST Data Management System.
6 * This product includes software developed by the LSST Project
7 * (https://www.lsst.org).
8 * See the COPYRIGHT file at the top-level directory of this distribution
9 * for details of code ownership.
10 *
11 * This program is free software: you can redistribute it and/or modify
12 * it under the terms of the GNU General Public License as published by
13 * the Free Software Foundation, either version 3 of the License, or
14 * (at your option) any later version.
15 *
16 * This program is distributed in the hope that it will be useful,
17 * but WITHOUT ANY WARRANTY; without even the implied warranty of
18 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19 * GNU General Public License for more details.
20 *
21 * You should have received a copy of the GNU General Public License
22 * along with this program. If not, see <https://www.gnu.org/licenses/>.
23 */
24
25#include <iostream>
26
27#include "lsst/jointcal/Frame.h"
28
29namespace lsst {
30namespace jointcal {
31
32using namespace std;
33
34/****************** Frame class methods ***********************/
35Frame::Frame(const Point &lowerLeft, const Point &upperRight) {
36 *this = Frame(lowerLeft.x, lowerLeft.y, upperRight.x, upperRight.y);
37}
38
39Frame::Frame(double xmin, double ymin, double xmax, double ymax) {
40 xMin = min(xmin, xmax);
41 xMax = max(xmin, xmax);
42 yMin = min(ymin, ymax);
43 yMax = max(ymin, ymax);
44}
45
46Frame::Frame() { xMin = xMax = yMin = yMax = 0; }
47
48/* positive if inside, negative if outside */
49double Frame::minDistToEdges(const Point &point) const {
50 return min(min(point.x - xMin, xMax - point.x) /* minx */,
51 min(point.y - yMin, yMax - point.y) /* miny */);
52}
53
55 Frame result = *this;
56 result *= right;
57 return result;
58}
59
61 Frame rightCopy = right;
62 // make sure that coordinates are properly ordered
63 order();
64 rightCopy.order();
65 xMin = max(xMin, rightCopy.xMin);
66 xMax = min(xMax, right.xMax);
67 yMin = max(yMin, right.yMin);
68 yMax = min(yMax, right.yMax);
69 // check for an actual overlap. Why was this check added?
70 if (xMin > xMax || yMin > yMax) *this = Frame();
71 return *this;
72}
73
75 Frame result = *this;
76 result += right;
77 return result;
78}
79
81 Frame rightCopy = right;
82 // make sure that coordinates are properly ordered
83 order();
84 rightCopy.order();
85 xMin = min(xMin, rightCopy.xMin);
86 xMax = max(xMax, right.xMax);
87 yMin = min(yMin, right.yMin);
88 yMax = max(yMax, right.yMax);
89 return *this;
90}
91
92void Frame::order() {
93 if (xMin > xMax) swap(xMin, xMax);
94 if (yMin > yMax) swap(yMin, yMax);
95}
96
97bool Frame::operator==(const Frame &right) const {
98 return ((xMin == right.xMin) && (xMax == right.xMax) && (yMin == right.yMin) && (yMax == right.yMax));
99}
100
101void Frame::cutMargin(const double marginX, const double marginY) {
102 xMin += marginX;
103 yMin += marginY;
104 xMax -= marginX;
105 yMax -= marginY;
106}
107
108void Frame::cutMargin(const double marginSize) { cutMargin(marginSize, marginSize); }
109
110Frame Frame::rescale(const double factor) const {
111 double hxsize = fabs(factor * 0.5 * (xMax - xMin));
112 double xcenter = 0.5 * (xMax + xMin);
113 double hysize = fabs(factor * 0.5 * (yMax - yMin));
114 double ycenter = 0.5 * (yMax + yMin);
115 return Frame(xcenter - hxsize, ycenter - hysize, xcenter + hxsize, ycenter + hysize);
116}
117
118double Frame::getArea() const { return fabs((xMax - xMin) * (yMax - yMin)); }
119
120bool Frame::inFrame(double x, double y) const {
121 return ((x <= xMax) && (y <= yMax) && (x >= xMin) && (y >= yMin));
122}
123
124void Frame::print(ostream &stream) const {
125 stream << "(xMin,yMin), (xMax,yMax): (" << std::setprecision(12) << xMin << ", " << yMin << "), (" << xMax
126 << ", " << yMax << ")";
127}
128} // namespace jointcal
129} // namespace lsst
py::object result
Definition _schema.cc:429
int min
int max
int y
Definition SpanSet.cc:48
int xmax
Definition SpanSet.cc:48
int xmin
Definition SpanSet.cc:48
rectangle with sides parallel to axes.
Definition Frame.h:38
void print(std::ostream &out) const
Definition Frame.cc:124
void cutMargin(double marginSize)
shrinks the frame (if marginSize>0), enlarges it (if marginSize<0).
Definition Frame.cc:108
bool inFrame(double x, double y) const
inside?
Definition Frame.cc:120
Frame operator*(const Frame &right) const
intersection of Frame's.
Definition Frame.cc:54
double minDistToEdges(const Point &point) const
distance to closest boundary.
Definition Frame.cc:49
double xMin
coordinate of boundary.
Definition Frame.h:41
Frame()
Default constructor.
Definition Frame.cc:46
bool operator==(const Frame &right) const
necessary for comparisons (!= is defined from this one implicitely)
Definition Frame.cc:97
Frame rescale(double factor) const
rescale it. The center does not move.
Definition Frame.cc:110
double getArea() const
Definition Frame.cc:118
Frame & operator*=(const Frame &right)
intersection of Frame's
Definition Frame.cc:60
Frame & operator+=(const Frame &right)
union of Frames
Definition Frame.cc:80
Frame operator+(const Frame &right) const
union of Frames
Definition Frame.cc:74
A point in a plane.
Definition Point.h:37
double x
coordinate
Definition Point.h:42
T fabs(T... args)
T right(T... args)
STL namespace.
T setprecision(T... args)
T swap(T... args)