CascadedPolygonUnion.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2011 Sandro Santilli <strk@kbt.io>
 * Copyright (C) 2006 Refractions Research Inc.
 *
 * 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.
 *
 **********************************************************************
 *
 * Last port: operation/union/CascadedPolygonUnion.java r487 (JTS-1.12+)
 * Includes custom code to deal with https://trac.osgeo.org/geos/ticket/837
 *
 **********************************************************************/
 
#pragma once
 
#include <vector>
 
#include <geos/export.h>
 
#include <geos/operation/union/UnionStrategy.h>
#include <geos/util/Progress.h>
 
// Forward declarations
namespace geos {
namespace geom {
class GeometryFactory;
class Geometry;
class Polygon;
class MultiPolygon;
class Envelope;
}
}
 
namespace geos {
namespace operation { // geos::operation
namespace geounion {  // geos::operation::geounion
 

class GEOS_DLL ClassicUnionStrategy : public UnionStrategy {
 
public:
 
    ClassicUnionStrategy() {};

    std::unique_ptr<geom::Geometry> Union(const geom::Geometry*, const geom::Geometry*) override;

    bool isFloatingPrecision() const override;
 
private:

    std::unique_ptr<geom::Geometry> unionPolygonsByBuffer(const geom::Geometry* g0, const geom::Geometry* g1);
 
};
 
 

class GEOS_DLL CascadedPolygonUnion {
private:
    std::vector<geom::Polygon*>* inputPolys;
    geom::GeometryFactory const* geomFactory;

    static int const STRTREE_NODE_CAPACITY = 4;

    static std::unique_ptr<geom::Geometry> restrictToPolygons(std::unique_ptr<geom::Geometry> g);
 
public:
    CascadedPolygonUnion();

    static std::unique_ptr<geom::Geometry> Union(std::vector<geom::Polygon*>* polys);
    static std::unique_ptr<geom::Geometry> Union(std::vector<geom::Polygon*>* polys, UnionStrategy* unionFun, geos::util::ProgressFunction* progressFunction);

    template <class T>
    static std::unique_ptr<geom::Geometry>
    Union(T start, T end, UnionStrategy *unionStrategy, geos::util::ProgressFunction* progressFunction)
    {
        std::vector<geom::Polygon*> polys;
        for(T i = start; i != end; ++i) {
            const geom::Polygon* p = dynamic_cast<const geom::Polygon*>(*i);
            polys.push_back(const_cast<geom::Polygon*>(p));
        }
        return Union(&polys, unionStrategy, progressFunction);
    }

    static std::unique_ptr<geom::Geometry> Union(const geom::MultiPolygon* polys, geos::util::ProgressFunction* progressFunction);

    CascadedPolygonUnion(std::vector<geom::Polygon*>* polys)
        : inputPolys(polys)
        , geomFactory(nullptr)
        , unionFunction(&defaultUnionFunction)
    {}
 
    CascadedPolygonUnion(std::vector<geom::Polygon*>* polys, UnionStrategy* unionFun)
        : inputPolys(polys)
        , geomFactory(nullptr)
        , unionFunction(unionFun)
    {}

    std::unique_ptr<geom::Geometry> Union(geos::util::ProgressFunction* progressFunction);
 
private:
 
    UnionStrategy* unionFunction;
    ClassicUnionStrategy defaultUnionFunction;

    std::unique_ptr<geom::Geometry> binaryUnion(
        const std::vector<const geom::Geometry*> & geoms,
        std::size_t start,
        std::size_t end,
        std::function<void()>* unitProgress);

    std::unique_ptr<geom::Geometry> unionSafe(const geom::Geometry* g0, const geom::Geometry* g1) const;
 
    std::unique_ptr<geom::Geometry> unionSafe(std::unique_ptr<geom::Geometry> &&, std::unique_ptr<geom::Geometry> &&);

    std::unique_ptr<geom::Geometry> unionActual(const geom::Geometry* g0, const geom::Geometry* g1) const;
 
    std::unique_ptr<geom::Geometry> unionActual(std::unique_ptr<geom::Geometry> &&, std::unique_ptr<geom::Geometry> &&) const;
};
 
 
 
 
 
} // namespace geos::operation::union
} // namespace geos::operation
} // namespace geos