public static CalculateReshapeLinesResponse CalculateReshapeLines(
[NotNull] CalculateReshapeLinesRequest request,
[CanBeNull] ITrackCancel trackCancel)
{
GetFeatures(request.SourceFeatures, request.TargetFeatures,
request.ClassDefinitions,
out IList <IFeature> sourceFeatures, out IList <IFeature> targetFeatures);
IList <CutSubcurve> resultLines = CalculateReshapeLines(
sourceFeatures, targetFeatures, request, trackCancel,
out ReshapeAlongCurveUsability usability);
CalculateReshapeLinesResponse result =
PackCalculateReshapeLinesResponse(usability, resultLines);
return(result);
}
public override async Task <CalculateReshapeLinesResponse> CalculateReshapeLines(
[NotNull] CalculateReshapeLinesRequest request,
[NotNull] ServerCallContext context)
{
Stopwatch watch = _msg.DebugStartTiming();
Func <ITrackCancel, CalculateReshapeLinesResponse> func =
trackCancel =>
ChangeAlongServiceUtils.CalculateReshapeLines(request, trackCancel);
CalculateReshapeLinesResponse response =
await GrpcServerUtils.ExecuteServiceCall(func, context, _staTaskScheduler, true) ??
new CalculateReshapeLinesResponse();
_msg.DebugStopTiming(
watch,
"Calculated reshape lines for peer {0} ({1} source features, {2}, target features)",
context.Peer, request.SourceFeatures.Count, request.TargetFeatures.Count);
return(response);
}
private static CalculateReshapeLinesResponse PackCalculateReshapeLinesResponse(
ReshapeAlongCurveUsability usability,
[NotNull] IList <CutSubcurve> resultLines)
{
Stopwatch watch = Stopwatch.StartNew();
var result = new CalculateReshapeLinesResponse
{
ReshapeLinesUsability = (int)usability
};
foreach (CutSubcurve resultCurve in resultLines)
{
result.ReshapeLines.Add(ToReshapeLineMsg(resultCurve));
}
_msg.DebugStopTiming(
watch, "Packed {0} reshape along calculation results into response",
resultLines.Count);
return(result);
}
public void CanReshapeAlongFilteredLinesByNoTargetOverlap()
{
GetOverlappingPolygons(out GdbFeature sourceFeature, out GdbFeature targetFeature);
CalculateReshapeLinesRequest calculationRequest =
CreateCalculateReshapeLinesRequest(sourceFeature, targetFeature);
//
// Filter by 'no target overlaps':
//
calculationRequest.FilterOptions = new ReshapeLineFilterOptionsMsg();
calculationRequest.FilterOptions.ExcludeResultingInOverlaps = true;
CalculateReshapeLinesResponse calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual(ReshapeAlongCurveUsability.CanReshape,
(ReshapeAlongCurveUsability)calculateResponse.ReshapeLinesUsability);
// 1 inside-source line and 1 (reshapable, but filtered) outside-source line.
AssertReshapeLineCount(calculateResponse.ReshapeLines, 2, 2, 1);
int insideLineIndex =
GetInsideReshapeLineIndex(calculateResponse.ReshapeLines, sourceFeature.Shape);
IPolyline insideLine =
(IPolyline)ProtobufGeometryUtils.FromShapeMsg(
calculateResponse.ReshapeLines[insideLineIndex].Path);
Assert.NotNull(insideLine);
Assert.AreEqual(1000, insideLine.Length);
//
// Reshape the default-side:
ApplyReshapeLinesRequest applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = false;
applyRequest.UseNonDefaultReshapeSide = false;
ApplyReshapeLinesResponse applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
GdbObjectMsg updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
IGeometry updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.AreEqual(1000 * 1000 * 3 / 4, ((IArea)updatedGeometry).Area);
// Check the new reshape line - there should be 1 reshapable, but filtered line
AssertReshapeLineCount(applyResponse.NewReshapeLines, 1, 1, 1);
// Technically the curve could be used to reshape (filtered out is not evaluated for usability)
Assert.AreEqual(ReshapeAlongCurveUsability.CanReshape,
(ReshapeAlongCurveUsability)applyResponse.ReshapeLinesUsability);
// Reshape the non-default side. This should probably not be possible but as it is currently
// just labelled a reshape line filter, it can be justified to reshape anyway...
applyRequest.UseNonDefaultReshapeSide = true;
applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.AreEqual(1000 * 1000 * 1 / 4, ((IArea)updatedGeometry).Area);
// Check the new reshape line - there should be 1 reshapable, but filtered line
AssertReshapeLineCount(applyResponse.NewReshapeLines, 1, 1, 1);
// Technically the curve could be used to reshape (filtered out is not evaluated for usability)
Assert.AreEqual(ReshapeAlongCurveUsability.CanReshape,
(ReshapeAlongCurveUsability)applyResponse.ReshapeLinesUsability);
}
public void CanReshapeAlongFilteredLinesByExtent()
{
GetOverlappingPolygons(out GdbFeature sourceFeature, out GdbFeature targetFeature);
CalculateReshapeLinesRequest calculationRequest =
CreateCalculateReshapeLinesRequest(sourceFeature, targetFeature);
IEnvelope visibleExtent = sourceFeature.Shape.Envelope;
visibleExtent.Expand(10, 10, false);
calculationRequest.FilterOptions = new ReshapeLineFilterOptionsMsg();
calculationRequest.FilterOptions.ClipLinesOnVisibleExtent = true;
calculationRequest.FilterOptions.VisibleExtents.Add(
ProtobufGeometryUtils.ToEnvelopeMsg(visibleExtent));
CalculateReshapeLinesResponse calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
calculateResponse.ReshapeLinesUsability);
// 1 inside-line and 2 extent-intersecting dangles to the outside within the extent.
AssertReshapeLineCount(calculateResponse.ReshapeLines, 3, 1);
int insideLineIndex =
GetInsideReshapeLineIndex(calculateResponse.ReshapeLines, sourceFeature.Shape);
IPolyline insideLine =
(IPolyline)ProtobufGeometryUtils.FromShapeMsg(
calculateResponse.ReshapeLines[insideLineIndex].Path);
Assert.NotNull(insideLine);
Assert.AreEqual(1000, insideLine.Length);
//
// Reshape the non-default side (should be possible):
//
ApplyReshapeLinesRequest applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = false;
applyRequest.UseNonDefaultReshapeSide = true;
ApplyReshapeLinesResponse applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
GdbObjectMsg updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
GdbObjectReference resultFeatureObjRef = new GdbObjectReference(
updatedFeatureMsg.ClassHandle, updatedFeatureMsg.ObjectId);
Assert.AreEqual(new GdbObjectReference(sourceFeature), resultFeatureObjRef);
IGeometry updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.AreEqual(1000 * 1000 * 1 / 4, ((IArea)updatedGeometry).Area);
// Check the new reshape line - there should be 2 non-reshapable dangles
AssertReshapeLineCount(applyResponse.NewReshapeLines, 2, 0);
Assert.AreEqual(ReshapeAlongCurveUsability.InsufficientOrAmbiguousReshapeCurves,
(ReshapeAlongCurveUsability)applyResponse.ReshapeLinesUsability);
}
public void CanReshapeAlongNonDefaultSide()
{
GetOverlappingPolygons(out GdbFeature sourceFeature, out GdbFeature targetFeature);
CalculateReshapeLinesRequest calculationRequest =
CreateCalculateReshapeLinesRequest(sourceFeature, targetFeature);
CalculateReshapeLinesResponse calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
calculateResponse.ReshapeLinesUsability);
AssertReshapeLineCount(calculateResponse.ReshapeLines, 2, 2);
List <ShapeMsg> reshapePathMsgs =
calculateResponse.ReshapeLines.Select(l => l.Path).ToList();
List <IPolyline> resultLines =
ProtobufGeometryUtils.FromShapeMsgList <IPolyline>(reshapePathMsgs);
int insideLineIndex = resultLines
.Select((reshapePath, index) => (reshapePath, index))
.First(rl => GeometryUtils.InteriorIntersects(
rl.reshapePath, sourceFeature.Shape)).index;
var insideLines =
resultLines.Where(rl => GeometryUtils.InteriorIntersects(rl, sourceFeature.Shape))
.ToList();
Assert.AreEqual(1, insideLines.Count);
Assert.AreEqual(1000, (insideLines[0]).Length);
//
// Reshape the default side:
//
ApplyReshapeLinesRequest applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = false;
applyRequest.UseNonDefaultReshapeSide = false;
ApplyReshapeLinesResponse applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
GdbObjectMsg updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
GdbObjectReference resultFeatureObjRef = new GdbObjectReference(
updatedFeatureMsg.ClassHandle, updatedFeatureMsg.ObjectId);
Assert.AreEqual(new GdbObjectReference(sourceFeature), resultFeatureObjRef);
IGeometry updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.AreEqual(1000 * 1000 * 3 / 4, ((IArea)updatedGeometry).Area);
// Check the new reshape line:
AssertReshapeLineCount(applyResponse.NewReshapeLines, 1, 1);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
applyResponse.ReshapeLinesUsability);
//
// Reshape the non-default side:
//
applyRequest.UseNonDefaultReshapeSide = true;
applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
resultFeatureObjRef = new GdbObjectReference(
updatedFeatureMsg.ClassHandle, updatedFeatureMsg.ObjectId);
Assert.AreEqual(new GdbObjectReference(sourceFeature), resultFeatureObjRef);
updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.AreEqual((double)1000 * 1000 * 1 / 4, ((IArea)updatedGeometry).Area);
}
public void CanReshapeAlongBufferedTarget()
{
GetOverlappingPolygons(out GdbFeature sourceFeature, out GdbFeature targetFeature);
CalculateReshapeLinesRequest calculationRequest =
CreateCalculateReshapeLinesRequest(sourceFeature, targetFeature);
double minSegmentLength = 2;
calculationRequest.TargetBufferOptions = new TargetBufferOptionsMsg();
calculationRequest.TargetBufferOptions.BufferDistance = 10;
calculationRequest.TargetBufferOptions.BufferMinimumSegmentLength = minSegmentLength;
CalculateReshapeLinesResponse calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
calculateResponse.ReshapeLinesUsability);
AssertReshapeLineCount(calculateResponse.ReshapeLines, 4, 4);
List <ShapeMsg> reshapePathMsgs =
calculateResponse.ReshapeLines.Select(l => l.Path).ToList();
List <IPolyline> resultLines =
ProtobufGeometryUtils.FromShapeMsgList <IPolyline>(reshapePathMsgs);
int insideLineIndex = resultLines
.Select((reshapePath, index) => (reshapePath, index))
.OrderBy(rl => rl.reshapePath.Length)
.First(rl => GeometryUtils.InteriorIntersects(
rl.reshapePath, sourceFeature.Shape)).index;
var insideLines =
resultLines.Where(rl => GeometryUtils.InteriorIntersects(rl, sourceFeature.Shape))
.ToList();
Assert.AreEqual(2, insideLines.Count);
foreach (IPolyline resultLine in resultLines)
{
Assert.AreNotEqual(1000, resultLine.Length);
Assert.AreEqual(
0, GeometryUtils.GetShortSegments(resultLine, minSegmentLength).Count);
}
//
// Reshape the default side:
//
ApplyReshapeLinesRequest applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = false;
applyRequest.UseNonDefaultReshapeSide = false;
ApplyReshapeLinesResponse applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
GdbObjectMsg updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
GdbObjectReference resultFeatureObjRef = new GdbObjectReference(
updatedFeatureMsg.ClassHandle, updatedFeatureMsg.ObjectId);
Assert.AreEqual(new GdbObjectReference(sourceFeature), resultFeatureObjRef);
IGeometry updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
// The shortest reshape line results in a greater remaining area:
Assert.IsNotNull(updatedGeometry);
Assert.Greater(((IArea)updatedGeometry).Area, 1000 * 1000 * 3 / (double)4);
// Check the new reshape line:
AssertReshapeLineCount(applyResponse.NewReshapeLines, 3, 3);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
applyResponse.ReshapeLinesUsability);
}
public void CanReshapeAlongMinimalTolerance()
{
GetOverlappingPolygons(out GdbFeature sourceFeature, out GdbFeature targetFeature);
// Insert an extra point:
IPoint targetAlmostIntersectPoint = GeometryFactory.CreatePoint(
2601000.001, 1200500.0, sourceFeature.Shape.SpatialReference);
bool pointInserted = ReshapeUtils.EnsurePointsExistInTarget(
targetFeature.Shape,
new[] { targetAlmostIntersectPoint }, 0.001);
Assert.IsTrue(pointInserted);
Assert.AreEqual(6, GeometryUtils.GetPointCount(targetFeature.Shape));
CalculateReshapeLinesRequest calculationRequest =
CreateCalculateReshapeLinesRequest(sourceFeature, targetFeature);
CalculateReshapeLinesResponse calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
calculateResponse.ReshapeLinesUsability);
AssertReshapeLineCount(calculateResponse.ReshapeLines, 2, 2);
int insideLineIndex =
GetInsideReshapeLineIndex(calculateResponse.ReshapeLines, sourceFeature.Shape);
IPolyline insideLine =
(IPolyline)ProtobufGeometryUtils.FromShapeMsg(
calculateResponse.ReshapeLines[insideLineIndex].Path);
Assert.NotNull(insideLine);
Assert.AreNotEqual(1000.0, (insideLine).Length);
// NOTE: If IntersectionUtils.UseCustomIntersect = true the from point is exactly on the inserted point
// Otherwise it's somewhere between - so we cannot do the exact comparison
Assert.AreNotEqual(insideLine.FromPoint.X, 2601000.0);
Assert.AreNotEqual(insideLine.ToPoint.X, 2601000.0);
//
// Reshape the default side:
ApplyReshapeLinesRequest applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = false;
applyRequest.UseNonDefaultReshapeSide = false;
ApplyReshapeLinesResponse applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
GdbObjectMsg updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
IGeometry updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.Greater(((IArea)updatedGeometry).Area, 1000.0 * 1000.0 * 3 / 4);
// Check the new reshape line:
AssertReshapeLineCount(applyResponse.NewReshapeLines, 1, 1);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
applyResponse.ReshapeLinesUsability);
//
//
// The same using the minimum tolerance:
//
calculationRequest.Tolerance = 0;
calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
calculateResponse.ReshapeLinesUsability);
AssertReshapeLineCount(calculateResponse.ReshapeLines, 2, 2);
insideLineIndex =
GetInsideReshapeLineIndex(calculateResponse.ReshapeLines, sourceFeature.Shape);
insideLine =
(IPolyline)ProtobufGeometryUtils.FromShapeMsg(
calculateResponse.ReshapeLines[insideLineIndex].Path);
Assert.NotNull(insideLine);
Assert.AreEqual(1000.0, insideLine.Length, 0.000000001);
// NOTE: If IntersectionUtils.UseCustomIntersect = true the from point is exactly on the inserted point
// Otherwise it's somewhere between - so we cannot do the exact comparison
Assert.AreEqual(2601000.0, insideLine.FromPoint.X);
Assert.AreEqual(1200500.0, insideLine.FromPoint.Y);
//
// Apply with minimum tolerance:
applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = false;
applyRequest.UseNonDefaultReshapeSide = false;
applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(1, applyResponse.ResultFeatures.Count);
updatedFeatureMsg = applyResponse.ResultFeatures[0].Update;
updatedGeometry = ProtobufGeometryUtils.FromShapeMsg(updatedFeatureMsg.Shape);
Assert.IsNotNull(updatedGeometry);
Assert.AreEqual(1000 * 1000 * 3 / 4, ((IArea)updatedGeometry).Area);
// Check the new remaining reshape line (outside)
AssertReshapeLineCount(applyResponse.NewReshapeLines, 1, 1);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
applyResponse.ReshapeLinesUsability);
IPolyline outsideLine = (IPolyline)ProtobufGeometryUtils.FromShapeMsg(
applyResponse.NewReshapeLines[0].Path);
Assert.NotNull(outsideLine);
Assert.AreEqual(3000.0, outsideLine.Length, 0.000000001);
// NOTE: If IntersectionUtils.UseCustomIntersect = true the from point is exactly on the inserted point
// Otherwise it's somewhere between - so we cannot do the exact comparison
Assert.AreEqual(2601000.0, outsideLine.ToPoint.X);
Assert.AreEqual(1200500.0, outsideLine.ToPoint.Y);
}
public void CanReshapeAlongInsertTargetVertices()
{
GetOverlappingPolygons(out GdbFeature sourceFeature, out GdbFeature targetFeature);
CalculateReshapeLinesRequest calculationRequest =
CreateCalculateReshapeLinesRequest(sourceFeature, targetFeature);
CalculateReshapeLinesResponse calculateResponse =
ChangeAlongServiceUtils.CalculateReshapeLines(calculationRequest, null);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
calculateResponse.ReshapeLinesUsability);
AssertReshapeLineCount(calculateResponse.ReshapeLines, 2, 2);
int insideLineIndex =
GetInsideReshapeLineIndex(calculateResponse.ReshapeLines, sourceFeature.Shape);
IPolyline insideLine =
(IPolyline)ProtobufGeometryUtils.FromShapeMsg(
calculateResponse.ReshapeLines[insideLineIndex].Path);
Assert.NotNull(insideLine);
Assert.AreEqual(1000.0, insideLine.Length);
//
// Reshape the default side:
//
ApplyReshapeLinesRequest applyRequest = new ApplyReshapeLinesRequest();
applyRequest.ReshapeLines.Add(calculateResponse.ReshapeLines[insideLineIndex]);
applyRequest.CalculationRequest = calculationRequest;
applyRequest.InsertVerticesInTarget = true;
applyRequest.UseNonDefaultReshapeSide = false;
ApplyReshapeLinesResponse applyResponse =
ChangeAlongServiceUtils.ApplyReshapeLines(applyRequest, null);
Assert.AreEqual(2, applyResponse.ResultFeatures.Count);
ResultObjectMsg sourceResultMsg =
applyResponse.ResultFeatures.First(
f => f.Update.ObjectId == sourceFeature.OID);
IGeometry updatedSourceGeometry =
ProtobufGeometryUtils.FromShapeMsg(sourceResultMsg.Update.Shape);
Assert.IsNotNull(updatedSourceGeometry);
Assert.AreEqual(1000 * 1000 * 3 / 4, ((IArea)updatedSourceGeometry).Area);
ResultObjectMsg targetResultMsg =
applyResponse.ResultFeatures.First(
f => f.Update.ObjectId == targetFeature.OID);
IGeometry updatedTargetGeometry =
ProtobufGeometryUtils.FromShapeMsg(targetResultMsg.Update.Shape);
Assert.AreEqual(GeometryUtils.GetPointCount(targetFeature.Shape) + 2,
GeometryUtils.GetPointCount(updatedTargetGeometry));
// Check the new reshape line:
AssertReshapeLineCount(applyResponse.NewReshapeLines, 1, 1);
Assert.AreEqual((int)ReshapeAlongCurveUsability.CanReshape,
applyResponse.ReshapeLinesUsability);
}
