///// <summary>
///// Assumes input is valid
///// (e.g. <paramref name="start" /> minor or equals to <paramref name="end" />).
///// </summary>
///// <param name="start"></param>
///// <param name="end"></param>
///// <returns></returns>
//private ILineString ComputeLine(LinearLocation start, LinearLocation end)
//{
// var coordinates = _line.Coordinates;
// var newCoordinates = new CoordinateList();
// var startSegmentIndex = start.SegmentIndex;
// if (start.SegmentFraction > 0.0)
// startSegmentIndex += 1;
// var lastSegmentIndex = end.SegmentIndex;
// if (end.SegmentFraction == 1.0)
// lastSegmentIndex += 1;
// if (lastSegmentIndex >= coordinates.Length)
// lastSegmentIndex = coordinates.Length - 1;
// // not needed - LinearLocation values should always be correct
// // Assert.IsTrue(end.SegmentFraction <= 1.0, "invalid segment fraction value");
// if (!start.IsVertex)
// newCoordinates.Add(start.GetCoordinate(_line));
// for (var i = startSegmentIndex; i <= lastSegmentIndex; i++)
// newCoordinates.Add(coordinates[i]);
// if (!end.IsVertex)
// newCoordinates.Add(end.GetCoordinate(_line));
// // ensure there is at least one coordinate in the result
// if (newCoordinates.Count <= 0)
// newCoordinates.Add(start.GetCoordinate(_line));
// var newCoordinateArray = newCoordinates.ToCoordinateArray();
// /*
// * Ensure there is enough coordinates to build a valid line.
// * Make a 2-point line with duplicate coordinates, if necessary.
// * There will always be at least one coordinate in the coordList.
// */
// if (newCoordinateArray.Length <= 1)
// newCoordinateArray = new[] { newCoordinateArray[0], newCoordinateArray[0] };
// return _line.Factory.CreateLineString(newCoordinateArray);
//}
/// <summary>
/// Assumes input is valid
/// (e.g. <paramref name="start" /> minor or equals to <paramref name="end" />).
/// </summary>
/// <param name="start"></param>
/// <param name="end"></param>
/// <returns></returns>
private IGeometry ComputeLinear(LinearLocation start, LinearLocation end)
{
var builder = new LinearGeometryBuilder(_line.Factory);
builder.FixInvalidLines = true;
if (!start.IsVertex)
{
builder.Add(start.GetCoordinate(_line));
}
for (var it = new LinearIterator(_line, start); it.HasNext(); it.Next())
{
if (end.CompareLocationValues(it.ComponentIndex, it.VertexIndex, 0.0) < 0)
{
break;
}
var pt = it.SegmentStart;
builder.Add(pt);
if (it.IsEndOfLine)
{
builder.EndLine();
}
}
if (!end.IsVertex)
{
builder.Add(end.GetCoordinate(_line));
}
return(builder.GetGeometry());
}
/// <summary>
/// Computes the <see cref="Coordinate" /> for the point
/// on the line at the given index.
/// If the index is out of range the first or last point on the
/// line will be returned.
/// </summary>
/// <remarks>
/// The Z-ordinate of the computed point will be interpolated from
/// the Z-ordinates of the line segment containing it, if they exist.
/// </remarks>
/// <param name="index">The index of the desired point.</param>
/// <returns>The <see cref="Coordinate" /> at the given index.</returns>
public Coordinate ExtractPoint(double index)
{
LinearLocation loc = LengthLocationMap.GetLocation(_linearGeom, index);
return(loc.GetCoordinate(_linearGeom));
}
/// <summary>
/// Computes the <see cref="Coordinate" />for the point on the line at the given index.
/// If the <paramref name="index" /> is out of range,
/// the first or last point on the line will be returned.
/// </summary>
/// <remarks>
/// The Z-ordinate of the computed point will be interpolated from
/// the Z-ordinates of the line segment containing it, if they exist.
/// </remarks>
/// <param name="index">The index of the desired point.</param>
/// <returns>The <see cref="Coordinate" /> at the given index.</returns>
public Coordinate ExtractPoint(LinearLocation index)
{
return(index.GetCoordinate(_linearGeom));
}
