LSST Applications g063fba187b+fee0456c91,g0f08755f38+ea96e5a5a3,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g1a2382251a+90257ff92a,g20f6ffc8e0+ea96e5a5a3,g217e2c1bcf+937a289c59,g28da252d5a+daa7da44eb,g2bbee38e9b+253935c60e,g2bc492864f+253935c60e,g3156d2b45e+6e55a43351,g32e5bea42b+31359a2a7a,g347aa1857d+253935c60e,g35bb328faa+a8ce1bb630,g3a166c0a6a+253935c60e,g3b1af351f3+a8ce1bb630,g3e281a1b8c+c5dd892a6c,g414038480c+416496e02f,g41af890bb2+afe91b1188,g599934f4f4+0db33f7991,g7af13505b9+e36de7bce6,g80478fca09+da231ba887,g82479be7b0+a4516e59e3,g858d7b2824+ea96e5a5a3,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,ga5288a1d22+bc6ab8dfbd,gb58c049af0+d64f4d3760,gc28159a63d+253935c60e,gcab2d0539d+3f2b72788c,gcf0d15dbbd+4ea9c45075,gda6a2b7d83+4ea9c45075,gdaeeff99f8+1711a396fd,ge79ae78c31+253935c60e,gef2f8181fd+3031e3cf99,gf0baf85859+c1f95f4921,gfa517265be+ea96e5a5a3,gfa999e8aa5+17cd334064,w.2024.50
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)