Box.h

Go to the documentation of this file.

// Box.h: Class representing a 2-dim box
//
// Copyright (C) 2020 ASTRON (Netherlands Institute for Radio Astronomy)
// SPDX-License-Identifier: GPL-3.0-or-later
 
 
#ifndef LOFAR_PARMDB_BOX_H
#define LOFAR_PARMDB_BOX_H
 
#include <casacore/casa/BasicMath/Math.h>
 
#include <cassert>
#include <utility>
#include <vector>
 
namespace dp3 {
namespace parmdb {
 
typedef std::pair<double, double> Point;
 
class Box;
Box unite(const Box& lhs, const Box& rhs);
Box intersect(const Box& lhs, const Box& rhs);
 
 
class Box {
 public:
  Box() : itsStart(0, 0), itsEnd(0, 0) {}
 
  Box(const Point& start, const Point& end) : itsStart(start), itsEnd(end) {
    assert(start.first <= end.first && start.second <= end.second);
  }
 
  Box(const std::vector<double>&);
 
  Box(double x1, double x2, double y1, double y2, bool asStartEnd = false);
 
  bool operator==(const Box& that) const {
    return itsStart == that.itsStart && itsEnd == that.itsEnd;
  }
  bool operator!=(const Box& that) const { return !operator==(that); }
 
  const Point& lower() const { return itsStart; }
  const Point& upper() const { return itsEnd; }
  double lowerX() const { return itsStart.first; }
  double lowerY() const { return itsStart.second; }
  double upperX() const { return itsEnd.first; }
  double upperY() const { return itsEnd.second; }
 
  double widthX() const { return itsEnd.first - itsStart.first; }
  double widthY() const { return itsEnd.second - itsStart.second; }
 
  bool intersects(const Box& other) const {
    return (other.itsStart.first < itsEnd.first &&
            !casacore::near(other.itsStart.first, itsEnd.first) &&
            other.itsEnd.first > itsStart.first &&
            !casacore::near(other.itsEnd.first, itsStart.first) &&
            other.itsStart.second < itsEnd.second &&
            !casacore::near(other.itsStart.second, itsEnd.second) &&
            other.itsEnd.second > itsStart.second &&
            !casacore::near(other.itsEnd.second, itsStart.second));
  }
 
  bool contains(const Box& other) const {
    return ((other.itsStart.first >= itsStart.first ||
             casacore::near(other.itsStart.first, itsStart.first)) &&
            (other.itsEnd.first <= itsEnd.first ||
             casacore::near(other.itsEnd.first, itsEnd.first)) &&
            (other.itsStart.second >= itsStart.second ||
             casacore::near(other.itsStart.second, itsStart.second)) &&
            (other.itsEnd.second <= itsEnd.second ||
             casacore::near(other.itsEnd.second, itsEnd.second)));
  }
 
  bool empty() const {
    return (casacore::near(itsStart.first, itsEnd.first) ||
            casacore::near(itsStart.second, itsEnd.second));
  }
 
  Box operator&(const Box& that) const { return intersect(*this, that); }
 
  Box operator|(const Box& that) const { return unite(*this, that); }
 
  bool operator<(const Box& that) const {
    return itsStart.second < that.itsStart.second ||
           (itsStart.second == that.itsStart.second &&
            itsStart.first < that.itsStart.first);
  }
  bool operator>(const Box& that) const {
    return itsStart.second > that.itsStart.second ||
           (itsStart.second == that.itsStart.second &&
            itsStart.first > that.itsStart.first);
  }
 
  void print() const;
 
 private:
  Point itsStart;
  Point itsEnd;
};
 
 
}  // namespace parmdb
}  // namespace dp3
 
#endif