/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2012 Excensus LLC.
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Licence as published
 * by the Free Software Foundation. 
 * See the COPYING file for more information.
 *
 **********************************************************************
 *
 * Last port: triangulate/DelaunayTriangulationBuilder.java rev. r524
 *
 **********************************************************************/

#include <geos/triangulate/DelaunayTriangulationBuilder.h>

#include <algorithm>

#include <geos/geom/GeometryFactory.h>
#include <geos/geom/Coordinate.h>
#include <geos/geom/CoordinateSequence.h>
#include <geos/triangulate/IncrementalDelaunayTriangulator.h>
#include <geos/triangulate/quadedge/QuadEdgeSubdivision.h>

namespace geos {
namespace triangulate { //geos.triangulate

using namespace geos::geom;

CoordinateSequence*
DelaunayTriangulationBuilder::extractUniqueCoordinates(
		const Geometry& geom)
{
	geom::CoordinateSequence *coords = geom.getCoordinates();
	unique(*coords);
	return coords;
}

void
DelaunayTriangulationBuilder::unique(CoordinateSequence& coords)
{
	std::vector<Coordinate> coordVector;
	coords.toVector(coordVector);
	std::sort(coordVector.begin(), coordVector.end(), geos::geom::CoordinateLessThen());
	coords.setPoints(coordVector);
	coords.removeRepeatedPoints();
}

IncrementalDelaunayTriangulator::VertexList*
DelaunayTriangulationBuilder::toVertices(
		const CoordinateSequence &coords)
{
	IncrementalDelaunayTriangulator::VertexList* vertexList =
		new IncrementalDelaunayTriangulator::VertexList();

	for(size_t iter=0; iter < coords.size(); ++iter)
	{
		vertexList->push_back(quadedge::Vertex(coords.getAt(iter)));
	}
	return vertexList;
}

DelaunayTriangulationBuilder::DelaunayTriangulationBuilder() :
	siteCoords(NULL), tolerance(0.0), subdiv(NULL)
{
}

DelaunayTriangulationBuilder::~DelaunayTriangulationBuilder() 
{
	if(siteCoords)
		delete siteCoords;
	if(subdiv)
		delete subdiv;
}

void
DelaunayTriangulationBuilder::setSites(const Geometry& geom)
{
	if(siteCoords)
		delete siteCoords;
	// remove any duplicate points (they will cause the triangulation to fail)
	siteCoords = extractUniqueCoordinates(geom);
}

void
DelaunayTriangulationBuilder::setSites(const CoordinateSequence& coords)
{
	if(siteCoords)
		delete siteCoords;
	siteCoords = coords.clone();
	// remove any duplicate points (they will cause the triangulation to fail)
	unique(*siteCoords);
}

void
DelaunayTriangulationBuilder::create()
{
	if(subdiv != NULL || siteCoords == NULL)
		return;

	Envelope siteEnv;
	siteCoords ->expandEnvelope(siteEnv);
	IncrementalDelaunayTriangulator::VertexList* vertices = toVertices(*siteCoords);
	subdiv = new quadedge::QuadEdgeSubdivision(siteEnv, tolerance);
	IncrementalDelaunayTriangulator triangulator = IncrementalDelaunayTriangulator(subdiv);
	triangulator.insertSites(*vertices);
	delete vertices;
}

quadedge::QuadEdgeSubdivision&
DelaunayTriangulationBuilder::getSubdivision()
{
	create();
	return *subdiv;
}

std::auto_ptr<MultiLineString>
DelaunayTriangulationBuilder::getEdges(
    const GeometryFactory& geomFact)
{
	create();
	return subdiv->getEdges(geomFact);
}

std::auto_ptr<geom::GeometryCollection>
DelaunayTriangulationBuilder::getTriangles(
		const geom::GeometryFactory& geomFact)
{
	create();
	return subdiv->getTriangles(geomFact);
}

geom::Envelope
DelaunayTriangulationBuilder::envelope(const geom::CoordinateSequence& coords)
{
	Envelope env;
	std::vector<Coordinate> coord_vector;
	coords.toVector(coord_vector);
	for(std::vector<Coordinate>::iterator it= coord_vector.begin() ; it!=coord_vector.end() ; ++it)
	{   
		const Coordinate& coord = *it;
		env.expandToInclude(coord);
	}   
	return env;
}


} //namespace geos.triangulate
} //namespace goes