LocationIndexedLine.h

/**********************************************************************
 *
 * 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: linearref/LocationIndexedLine.java r466
 *
 **********************************************************************/
 
#ifndef GEOS_LINEARREF_LOCATIONINDEXEDLINE_H
#define GEOS_LINEARREF_LOCATIONINDEXEDLINE_H
 
#include <geos/export.h>
#include <geos/geom/Coordinate.h>
#include <geos/geom/Geometry.h>
#include <geos/geom/Lineal.h>
#include <geos/linearref/LinearLocation.h>
#include <geos/linearref/LocationIndexOfPoint.h>
#include <geos/linearref/LocationIndexOfLine.h>
#include <geos/util/IllegalArgumentException.h>
 
namespace geos {
namespace linearref { // geos::linearref

class GEOS_DLL LocationIndexedLine
{
private:
  const geom::Geometry *linearGeom;
 
  void checkGeometryType()
  {
    if ( ! dynamic_cast<const geom::Lineal*>(linearGeom) )
      throw util::IllegalArgumentException("Input geometry must be linear");
  }
 
public:

  LocationIndexedLine(const geom::Geometry *linearGeom)
    : linearGeom(linearGeom)
  {
    checkGeometryType();
  }

  geom::Coordinate extractPoint(const LinearLocation& index) const
  {
    return index.getCoordinate(linearGeom);
  }
 

  geom::Coordinate extractPoint(const LinearLocation& index,
                                               double offsetDistance) const
  {
    geom::Coordinate ret;
    index.getSegment(linearGeom)->pointAlongOffset(
                index.getSegmentFraction(), offsetDistance, ret
           );
    return ret;
  }

  geom::Geometry *extractLine(const LinearLocation& startIndex,
                              const LinearLocation& endIndex) const
  {
    return ExtractLineByLocation::extract(linearGeom, startIndex, endIndex);
  }
 

  LinearLocation indexOf(const geom::Coordinate& pt) const
  {
    return LocationIndexOfPoint::indexOf(linearGeom, pt);
  }

  LinearLocation indexOfAfter(const geom::Coordinate& pt,
                                const LinearLocation& minIndex) const
  {
    return LocationIndexOfPoint::indexOfAfter(linearGeom, pt, &minIndex);
  }

  LinearLocation* indicesOf(const geom::Geometry *subLine) const
  {
    return LocationIndexOfLine::indicesOf(linearGeom, subLine);
  }
 

  LinearLocation project(const geom::Coordinate& pt) const
  {
    return LocationIndexOfPoint::indexOf(linearGeom, pt);
  }

  LinearLocation getStartIndex() const
  {
    return LinearLocation();
  }

  LinearLocation getEndIndex() const
  {
    return LinearLocation::getEndLocation(linearGeom);
  }

  bool isValidIndex(const LinearLocation& index) const
  {
    return index.isValid(linearGeom);
  }
 

  LinearLocation clampIndex(const LinearLocation& index) const
  {
    LinearLocation loc = index;
    loc.clamp(linearGeom);
    return loc;
  }
};
 
} // geos::linearref
} // geos
#endif