/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2009-2010  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/overlay/snap/GeometrySnapper.java r320 (JTS-1.12)
 *
 **********************************************************************/

#pragma once

#include <geos/geom/Coordinate.h>

#include <memory>
#include <vector>

// Forward declarations
namespace geos {
namespace geom {
//class PrecisionModel;
class Geometry;
class CoordinateSequence;
struct GeomPtrPair;
}
}

namespace geos {
namespace operation { // geos::operation
namespace overlay { // geos::operation::overlay
namespace snap { // geos::operation::overlay::snap

/** \brief
 * Snaps the vertices and segments of a {@link geom::Geometry}
 * to another Geometry's vertices.
 *
 * Where possible, this operation tries to avoid creating invalid geometries;
 * however, it does not guarantee that output geometries will be valid.  It is
 * the responsibility of the caller to check for and handle invalid geometries.
 *
 * Because too much snapping can result in invalid
 * topology being created, heuristics are used to determine the number and
 * location of snapped vertices that are likely safe to snap. These heuristics
 * may omit some potential snaps that are otherwise within the tolerance.
 *
 */
class GEOS_DLL GeometrySnapper {

public:

    typedef std::unique_ptr<geom::Geometry> GeomPtr;

    /**
     * Snaps two geometries together with a given tolerance.
     *
     * @param g0 a geometry to snap
     * @param g1 a geometry to snap
     * @param snapTolerance the tolerance to use
     * @param ret the snapped geometries as a pair of smart pointers
     *            (output parameter)
     */
    static void snap(const geom::Geometry& g0,
                     const geom::Geometry& g1,
                     double snapTolerance, geom::GeomPtrPair& ret);

    static GeomPtr snapToSelf(const geom::Geometry& g0,
                              double snapTolerance, bool cleanResult);

    /**
     * Creates a new snapper acting on the given geometry
     *
     * @param g the geometry to snap
     */
    GeometrySnapper(const geom::Geometry& g)
        :
        srcGeom(g)
    {
    }

    /** \brief
     * Snaps the vertices in the component {@link geom::LineString}s
     * of the source geometry to the vertices of the given snap geometry
     * with a given snap tolerance
     *
     * @param g a geometry to snap the source to
     * @param snapTolerance
     * @return a new snapped Geometry
     */
    std::unique_ptr<geom::Geometry> snapTo(const geom::Geometry& g,
                                           double snapTolerance);

    /** \brief
     * Snaps the vertices in the component {@link geom::LineString}s
     * of the source geometry to the vertices of itself
     * with a given snap tolerance and optionally cleaning the result.
     *
     * @param snapTolerance
     * @param cleanResult clean the result
     * @return a new snapped Geometry
     */
    std::unique_ptr<geom::Geometry> snapToSelf(double snapTolerance,
            bool cleanResult);

    /** \brief
     * Estimates the snap tolerance for a Geometry, taking into account
     * its precision model.
     *
     * @param g a Geometry
     * @return the estimated snap tolerance
     */
    static double computeOverlaySnapTolerance(const geom::Geometry& g);

    static double computeSizeBasedSnapTolerance(const geom::Geometry& g);

    /** \brief
     * Computes the snap tolerance based on input geometries;
     */
    static double computeOverlaySnapTolerance(const geom::Geometry& g1,
            const geom::Geometry& g2);


private:

    // eventually this will be determined from the geometry topology
    //static const double snapTol; //  = 0.000001;

    static const double snapPrecisionFactor; //  = 10e-10

    const geom::Geometry& srcGeom;

    /// Extract target (unique) coordinates
    std::unique_ptr<geom::Coordinate::ConstVect> extractTargetCoordinates(
        const geom::Geometry& g);

    // Declare type as noncopyable
    GeometrySnapper(const GeometrySnapper& other) = delete;
    GeometrySnapper& operator=(const GeometrySnapper& rhs) = delete;
};


} // namespace geos::operation::overlay::snap
} // namespace geos::operation::overlay
} // namespace geos::operation
} // namespace geos