/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2011 Sandro Santilli <strk@keybit.net
 *
 * 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/PointGeometryUnion.java r320 (JTS-1.12)
 *
 **********************************************************************/

#include <memory> // for auto_ptr
#include <cassert> // for assert
#include <algorithm> // for copy
#include <geos/operation/union/PointGeometryUnion.h> 
#include <geos/geom/Coordinate.h> 
#include <geos/geom/Puntal.h> 
#include <geos/geom/Point.h>
#include <geos/geom/MultiPoint.h>
#include <geos/geom/Geometry.h>
#include <geos/geom/Location.h>
#include <geos/geom/GeometryFactory.h>
#include <geos/geom/util/GeometryCombiner.h>
#include <geos/algorithm/PointLocator.h>

namespace geos {
namespace operation { // geos::operation
namespace geounion {  // geos::operation::geounion

/* public */
std::auto_ptr<geom::Geometry>
PointGeometryUnion::Union() const
{
  using namespace geom;
  using algorithm::PointLocator;
  using geom::util::GeometryCombiner;

  PointLocator locater;
  // use a set to eliminate duplicates, as required for union
  std::set<Coordinate> exteriorCoords;

  for (std::size_t i=0, n=pointGeom.getNumGeometries(); i<n; ++i) {
      const Point* point = dynamic_cast<const Point*>(pointGeom.getGeometryN(i));
      assert(point);
      const Coordinate* coord = point->getCoordinate();
      int loc = locater.locate(*coord, &otherGeom);
      if (loc == Location::EXTERIOR)
              exteriorCoords.insert(*coord);
  }

  // if no points are in exterior, return the other geom
  if (exteriorCoords.empty())
          return std::auto_ptr<Geometry>(otherGeom.clone());

  // make a puntal geometry of appropriate size
  std::auto_ptr<Geometry> ptComp;

  if (exteriorCoords.size() == 1) {
    ptComp.reset( geomFact->createPoint(*(exteriorCoords.begin())) );
  }
  else
  {
    std::vector<Coordinate> coords(exteriorCoords.size());
    std::copy(exteriorCoords.begin(), exteriorCoords.end(), coords.begin());
    ptComp.reset( geomFact->createMultiPoint(coords) );
  }

  // add point component to the other geometry
  return std::auto_ptr<Geometry> (
    GeometryCombiner::combine(ptComp.get(), &otherGeom)
  );
}

/* public  static */
std::auto_ptr<geom::Geometry>
PointGeometryUnion::Union(const geom::Puntal& pointGeom,
      const geom::Geometry& otherGeom)
{
  PointGeometryUnion unioner(pointGeom, otherGeom);
  return unioner.Union();
}

/* public */
PointGeometryUnion::PointGeometryUnion(const geom::Puntal& pointGeom_,
                     const geom::Geometry& otherGeom_)
  :
  pointGeom ( pointGeom_ ),
  otherGeom ( otherGeom_ )
{
  geomFact = otherGeom.getFactory();
}

} // namespace geos::operation::union
} // namespace geos::operation
} // namespace geos