Cells.h

Go to the documentation of this file.

//=============================================================================
//    Copyright (C) 2021-2024 Wageningen University - All Rights Reserved
//                     Author: Gonzalo Mier
//                        BSD-3 License
//=============================================================================
 
#pragma once
#ifndef FIELDS2COVER_TYPES_CELLS_H_
#define FIELDS2COVER_TYPES_CELLS_H_
 
#include <gdal/ogr_geometry.h>
#include <vector>
#include "fields2cover/types/Geometry.h"
#include "fields2cover/types/Geometries.h"
#include "fields2cover/types/Point.h"
#include "fields2cover/types/LineString.h"
#include "fields2cover/types/Cell.h"
 
namespace f2c::types {
 
struct Cells : public Geometries<Cells, OGRMultiPolygon, wkbMultiPolygon,
    Cell> {
 public:
  using Geometries<Cells, OGRMultiPolygon, wkbMultiPolygon, Cell>::Geometries;
  Cells();
  explicit Cells(const OGRGeometry* geom);
  explicit Cells(const Cell& c);
  ~Cells() = default;
 
  void getGeometry(size_t i, Cell& cell);
  void getGeometry(size_t i, Cell& cell) const;
  Cell getGeometry(size_t i);
  const Cell getGeometry(size_t i) const;
  void setGeometry(size_t i, const Cell& cell);
 
  const Cell getCell(size_t i) const;
  const LinearRing getCellBorder(size_t i) const;
  const LinearRing getInteriorRing(size_t i_cell, size_t i_ring) const;
 
 
  void operator*=(double b);
 
  void addGeometry(const Cell& c);
  void addRing(size_t i, const LinearRing& ring);
 
  size_t size() const;
 
  bool isConvex() const;
 
  Cell convexHull() const;
 
  static Cells intersection(const Cell& cell, const Cell& c);
  Cells intersection(const Cell& c) const;
  Cells intersection(const Cells& c) const;
 
  Cells difference(const Cell& c) const;
  Cells difference(const Cells& c) const;
 
  Cells unionOp(const Cell& c) const;
  Cells unionOp(const Cells& c) const;
 
  Cells unionCascaded() const;
 
  Cells splitByLine(const LineString& line) const;
 
  Cells splitByLine(const MultiLineString& lines) const;
 
  LineString createSemiLongLine(const Point& point, double angle) const;
 
  LineString createStraightLongLine(const Point& point, double angle) const;
 
  MultiLineString getLinesInside(const LineString& line) const;
 
  MultiLineString getLinesInside(const MultiLineString& lines) const;
 
  Cells getCellsInside(const Cells& cell) const;
 
  bool isPointInBorder(const Point& p) const;
 
  bool isPointIn(const Point& p) const;
 
  const Cell getCellWherePoint(const Point& p) const;
 
  LineString createLineUntilBorder(
      const f2c::types::Point& p, double ang) const;
 
  template <class T, OGRwkbGeometryType R>
  static Cells buffer(const Geometry<T, R>& geom, double width,
      int side = 0);
 
  Cells buffer(double width) const;
 
  Point closestPointOnBorderTo(const Point& p) const;
};
 
 
 
template <class T, OGRwkbGeometryType R>
Cells Cells::buffer(const Geometry<T, R>& geom, double width, int side) {
  OGRGeometry* buffer = geom.OGRBuffer(width, side);
  Cells cells {buffer};
  OGRGeometryFactory::destroyGeometry(buffer);
  return cells;
}
 
 
}  // namespace f2c::types
 
#endif  // FIELDS2COVER_TYPES_CELLS_H_