/*
* The JTS Topology Suite is a collection of Java classes that
* implement the fundamental operations required to validate a given
* geo-spatial data set to a known topological specification.
*
* Copyright (C) 2001 Vivid Solutions
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this library; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
* For more information, contact:
*
* Vivid Solutions
* Suite #1A
* 2328 Government Street
* Victoria BC V8T 5G5
* Canada
*
* (250)385-6040
* www.vividsolutions.com
*/
package com.revolsys.geometry.linearref;
import com.revolsys.geometry.model.Geometry;
import com.revolsys.geometry.model.Point;
import com.revolsys.geometry.model.coordinates.LineSegmentUtil;
import com.revolsys.geometry.util.Assert;
/**
* Computes the {@link LinearLocation} of the point
* on a linear {@link Geometry} nearest a given {@link Coordinates}.
* The nearest point is not necessarily unique; this class
* always computes the nearest point closest to
* the start of the geometry.
*/
class LocationIndexOfPoint {
public static LinearLocation indexOf(final Geometry linearGeom, final Point inputPt) {
final LocationIndexOfPoint locater = new LocationIndexOfPoint(linearGeom);
return locater.indexOf(inputPt);
}
public static LinearLocation indexOfAfter(final Geometry linearGeom, final Point inputPt,
final LinearLocation minIndex) {
final LocationIndexOfPoint locater = new LocationIndexOfPoint(linearGeom);
return locater.indexOfAfter(inputPt, minIndex);
}
private final Geometry linearGeom;
public LocationIndexOfPoint(final Geometry linearGeom) {
this.linearGeom = linearGeom;
}
/**
* Find the nearest location along a linear {@link Geometry} to a given point.
*
* @param inputPt the coordinate to locate
* @return the location of the nearest point
*/
public LinearLocation indexOf(final Point inputPt) {
return indexOfFromStart(inputPt, null);
}
/**
* Find the nearest {@link LinearLocation} along the linear {@link Geometry}
* to a given {@link Coordinates}
* after the specified minimum {@link LinearLocation}.
* If possible the location returned will be strictly greater than the
* <code>minLocation</code>.
* If this is not possible, the
* value returned will equal <code>minLocation</code>.
* (An example where this is not possible is when
* minLocation = [end of line] ).
*
* @param inputPt the coordinate to locate
* @param minIndex the minimum location for the point location
* @return the location of the nearest point
*/
public LinearLocation indexOfAfter(final Point inputPt, final LinearLocation minIndex) {
if (minIndex == null) {
return indexOf(inputPt);
}
// sanity check for minLocation at or past end of line
final LinearLocation endLoc = LinearLocation.getEndLocation(this.linearGeom);
if (endLoc.compareTo(minIndex) <= 0) {
return endLoc;
}
final LinearLocation closestAfter = indexOfFromStart(inputPt, minIndex);
/**
* Return the minDistanceLocation found.
* This will not be null, since it was initialized to minLocation
*/
Assert.isTrue(closestAfter.compareTo(minIndex) >= 0,
"computed location is before specified minimum location");
return closestAfter;
}
private LinearLocation indexOfFromStart(final Point inputPt, final LinearLocation minIndex) {
double minDistance = Double.MAX_VALUE;
int minComponentIndex = 0;
int minSegmentIndex = 0;
double minFrac = -1.0;
for (final LinearIterator it = new LinearIterator(this.linearGeom); it.hasNext(); it.next()) {
if (!it.isEndOfLine()) {
final Point p0 = it.getSegmentStart();
final Point p1 = it.getSegmentEnd();
final double segDistance = LineSegmentUtil.distanceLinePoint(p0, p1, inputPt);
final double segFrac = LineSegmentUtil.segmentFraction(p0, p1, inputPt);
final int candidateComponentIndex = it.getComponentIndex();
final int candidateSegmentIndex = it.getVertexIndex();
if (segDistance < minDistance) {
// ensure after minLocation, if any
if (minIndex == null || minIndex.compareLocationValues(candidateComponentIndex,
candidateSegmentIndex, segFrac) < 0) {
// otherwise, save this as new minimum
minComponentIndex = candidateComponentIndex;
minSegmentIndex = candidateSegmentIndex;
minFrac = segFrac;
minDistance = segDistance;
}
}
}
}
if (minDistance == Double.MAX_VALUE) {
// no minimum was found past minLocation, so return it
return new LinearLocation(minIndex);
}
// otherwise, return computed location
final LinearLocation loc = new LinearLocation(minComponentIndex, minSegmentIndex, minFrac);
return loc;
}
/**
* Computes the fraction of distance (in <tt>[0.0, 1.0]</tt>)
* that a point occurs along a line segment.
* If the point is beyond either ends of the line segment,
* the closest fractional value (<tt>0.0</tt> or <tt>1.0</tt>) is returned.
*
* @param seg the line segment to use
* @param inputPt the point
* @return the fraction along the line segment the point occurs
*/
/*
* // MD - no longer needed private static double segmentFraction( LineSegment seg, Point inputPt)
* { double segFrac = seg.projectionFactor(inputPt); if (segFrac < 0.0) segFrac = 0.0; else if
* (segFrac > 1.0) segFrac = 1.0; return segFrac; }
*/
}