Cell.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2018-2025 ISciences, LLC
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Public Licence as published
 * by the Free Software Foundation.
 * See the COPYING file for more information.
 *
 **********************************************************************/
 
#pragma once
 
#include <geos/geom/Envelope.h>
#include <geos/geom/Geometry.h>
 
#include <geos/operation/grid/Side.h>
#include <geos/operation/grid/Traversal.h>
 
namespace geos::operation::grid {
 
class Cell
{
 
  public:
    Cell(double xmin, double ymin, double xmax, double ymax)
      : m_box{ xmin, ymin, xmax, ymax }
    {
    }
 
    explicit Cell(const geom::Envelope& b)
      : m_box{ b }
    {
    }
 
    const geom::Envelope& box() const { return m_box; }
 
    double getWidth() const;
 
    double getHeight() const;
 
    double getArea() const;
 
    void forceExit();
 
    bool isDetermined() const;
 
    double getTraversalLength() const;
 
    double getCoveredFraction() const;
 
    std::unique_ptr<geom::Geometry> getCoveredPolygons(const geom::GeometryFactory&) const;
 
    Traversal& getLastTraversal();
 
    bool take(const geom::CoordinateXY& c, const geom::CoordinateXY* prev_original = nullptr);
 
  private:
    std::vector<const std::vector<geom::CoordinateXY>*> getCoordLists() const;
 
    enum class Location
    {
        INSIDE,
        OUTSIDE,
        BOUNDARY
    };
 
    geom::Envelope m_box;
 
    std::vector<Traversal> m_traversals;
 
    Side side(const geom::CoordinateXY& c) const;
 
    Location location(const geom::CoordinateXY& c) const;
 
    Traversal& traversal_in_progress();
};
 
}