public void CanDetermineInteriorLineIntersectsXY()
{
var line1 = new Line3D(new Pnt3D(0, 0, 0), new Pnt3D(10, 10, 0));
var line2 = new Line3D(new Pnt3D(0, 10, 0), new Pnt3D(10, 0, 0));
Assert.True(GeomRelationUtils.SourceInteriorIntersectsXY(line1, line2, 0));
Assert.True(GeomRelationUtils.SourceInteriorIntersectsXY(line2, line1, 0));
var touchingLine = new Line3D(line1.EndPoint, line2.StartPoint);
double tolerance = double.Epsilon;
Assert.False(
GeomRelationUtils.SourceInteriorIntersectsXY(line1, touchingLine, tolerance));
var disjointLine = new Line3D(new Pnt3D(5, 5.05, 17), new Pnt3D(6, 18, 2));
Assert.False(GeomRelationUtils.SourceInteriorIntersectsXY(line1, disjointLine, 0));
Assert.True(GeomRelationUtils.SourceInteriorIntersectsXY(line1, disjointLine, 0.1));
var otherLineInteriorIntersected = new Line3D(new Pnt3D(0, 20, 0),
new Pnt3D(20, 0, 0));
// however, this line is not interior-intersected:
Assert.False(
GeomRelationUtils.SourceInteriorIntersectsXY(
line1, otherLineInteriorIntersected, 0.01));
}
public void CanDetermineEqualMultipointXyChangedOrderWithDuplicates()
{
Multipoint <IPnt> mp1 = new Multipoint <IPnt>(
new[]
{
new Pnt3D(1234.5678, 9876.54321, 345.6),
new Pnt3D(234.5678, 987.54321, 355.6),
new Pnt3D(34.5678, 98.54321, 365.6),
new Pnt3D(4.5678, 9.54321, 375.6)
});
Multipoint <IPnt> mp2 = new Multipoint <IPnt>(mp1.GetPoints(0, null, true).Reverse());
mp2.AddPoint(new Pnt3D(34.5678, 98.54321, 365.6));
Assert.IsTrue(GeomRelationUtils.AreMultipointsEqualXY(mp1, mp2, 0.00001));
mp2 = new Multipoint <IPnt>(mp2.GetPoints(0, null, true).Reverse());
Assert.IsTrue(GeomRelationUtils.AreMultipointsEqualXY(mp1, mp2, 0.00001));
mp2.AddPoint(new Pnt3D(34.4678, 98.64321, 365.6));
Assert.IsFalse(GeomRelationUtils.AreMultipointsEqualXY(mp1, mp2, 0.00001));
Assert.IsTrue(GeomRelationUtils.AreBoundsEqual(mp1, mp2, 0.000001));
}
public void CanDetermineTouchesXY()
{
var pnts1 = new List <Pnt3D>();
var pnts2 = new List <Pnt3D>();
var pnts3 = new List <Pnt3D>();
// ring1: horizontal:
pnts1.Add(new Pnt3D(0, 0, 9));
pnts1.Add(new Pnt3D(0, 100, 9));
pnts1.Add(new Pnt3D(100, 100, 9));
pnts1.Add(new Pnt3D(100, 0, 9));
pnts1.Add(new Pnt3D(0, 0, 9));
// ring 2: also horizontal, adjacent, touching along line
pnts2.Add(new Pnt3D(100, 10, 18));
pnts2.Add(new Pnt3D(100, 90, 18));
pnts2.Add(new Pnt3D(200, 90, 18));
pnts2.Add(new Pnt3D(200, 10, 18));
pnts2.Add(new Pnt3D(100, 10, 18));
// ring 3: containing both rings, touching ring1
pnts3.Add(new Pnt3D(0, -10, 0));
pnts3.Add(new Pnt3D(0, 110, 0));
pnts3.Add(new Pnt3D(250, 110, 0));
pnts3.Add(new Pnt3D(250, -10, 0));
pnts3.Add(new Pnt3D(0, -10, 0));
var ring1 = new Linestring(pnts1);
var ring2 = new Linestring(pnts2);
var ring3 = new Linestring(pnts3);
bool disjoint;
Assert.IsTrue(GeomRelationUtils.TouchesXY(ring1, ring2, 0.00001, out disjoint));
Assert.False(disjoint);
// ring3 is exterior ring that contains ring1
Assert.False(GeomRelationUtils.TouchesXY(ring1, ring3, 0.00001, out disjoint, true));
Assert.False(disjoint);
// use orientation information (ring3
Assert.False(GeomRelationUtils.TouchesXY(ring1, ring3, 0.00001, out disjoint));
ring3.ReverseOrientation();
// Now ring3 is interior ring that contains ring1 (i.e. outside -> touches)
Assert.True(GeomRelationUtils.TouchesXY(ring1, ring3, 0.00001, out disjoint));
Assert.False(disjoint);
const bool disregardRingOrientation = true;
Assert.False(GeomRelationUtils.TouchesXY(ring1, ring3, 0.00001, out disjoint,
disregardRingOrientation));
Assert.False(disjoint);
// now ring2 becomes an inner ring and ring1 should not be touching it any more:
pnts2.Reverse();
var innerRing2 = new Linestring(pnts2);
Assert.IsFalse(GeomRelationUtils.TouchesXY(ring1, innerRing2, 0.00001, out disjoint));
}
public void CanReadWriteMultipointXy()
{
var points = new WKSPointZ[4];
points[0] = new WKSPointZ {
X = 2600000, Y = 1200000, Z = double.NaN
};
points[1] = new WKSPointZ {
X = 2600030, Y = 1200020, Z = double.NaN
};
points[2] = new WKSPointZ {
X = 2600020, Y = 1200030, Z = double.NaN
};
points[3] = new WKSPointZ {
X = 2600040, Y = 1200040, Z = double.NaN
};
ISpatialReference sr =
SpatialReferenceUtils.CreateSpatialReference(WellKnownHorizontalCS.LV95);
IMultipoint multipoint = GeometryFactory.CreateMultipoint(points, sr);
GeometryUtils.MakeNonZAware(multipoint);
WkbGeometryWriter writer = new WkbGeometryWriter();
byte[] wkb = writer.WriteMultipoint(multipoint);
// ArcObjects
byte[] arcObjectsWkb = GeometryUtils.ToWkb(multipoint);
Assert.AreEqual(wkb, arcObjectsWkb);
// Wkx
byte[] wkx = ToChristianSchwarzWkb(ToWkxMultipoint(points, Ordinates.Xy));
Assert.AreEqual(wkx, wkb);
// Bonus test: Geom
WkbGeomWriter geomWriter = new WkbGeomWriter();
Multipoint <IPnt> multipnt = GeometryConversionUtils.CreateMultipoint(multipoint);
byte[] wkbGeom = geomWriter.WriteMultipoint(multipnt, Ordinates.Xy);
Assert.AreEqual(wkb, wkbGeom);
WkbGeometryReader reader = new WkbGeometryReader();
IMultipoint restored = reader.ReadMultipoint(new MemoryStream(wkb));
Assert.IsTrue(GeometryUtils.AreEqual(multipoint, restored));
// Geom
WkbGeomReader geomReader = new WkbGeomReader();
Multipoint <IPnt> deserializedPnts =
geomReader.ReadMultiPoint(new MemoryStream(wkbGeom));
Assert.IsTrue(
GeomRelationUtils.AreMultipointsEqualXY(multipnt, deserializedPnts,
double.Epsilon));
}
public void CanDetermineTouchesXYVerticalRing()
{
var verticalRing = new Linestring(new List <Pnt3D>
{
new Pnt3D(0, 0, 0),
new Pnt3D(30, 100, 0),
new Pnt3D(30, 100, 55),
new Pnt3D(0, 0, 55),
new Pnt3D(0, 0, 0)
});
var adjacentRing = new Linestring(new List <Pnt3D>
{
new Pnt3D(0, 0, 0),
new Pnt3D(30, 100, 0),
new Pnt3D(30, 0, 0),
new Pnt3D(0, 0, 0)
});
bool disjoint;
Assert.True(GeomRelationUtils.TouchesXY(adjacentRing, verticalRing, 0.0001,
out disjoint, true,
true));
Assert.False(disjoint);
adjacentRing.ReverseOrientation();
Assert.True(GeomRelationUtils.TouchesXY(adjacentRing, verticalRing, 0.0001,
out disjoint, true,
true));
Assert.False(disjoint);
var containedVerticalRing = new Linestring(new List <Pnt3D>
{
new Pnt3D(0, 0, 0),
new Pnt3D(30, 80, 0),
new Pnt3D(30, 80, 55),
new Pnt3D(0, 0, 55),
new Pnt3D(0, 0, 0)
});
Assert.False(GeomRelationUtils.TouchesXY(adjacentRing, containedVerticalRing, 0.0001,
out disjoint, true, true));
}
public void CanDetermineContainsXY()
{
var ring1 = new List <Pnt3D>();
ring1.Add(new Pnt3D(0, 0, 10));
ring1.Add(new Pnt3D(0, 100, 20));
ring1.Add(new Pnt3D(100, 100, 10));
ring1.Add(new Pnt3D(50, 50, 10));
ring1.Add(new Pnt3D(25, 75, 10));
ring1.Add(new Pnt3D(0, 0, 10));
Linestring l = new Linestring(ring1);
Assert.IsFalse(GeomRelationUtils.PolycurveContainsXY(l, new Pnt3D(70, 50, 2), 0.001));
Assert.IsFalse(GeomRelationUtils.PolycurveContainsXY(l, new Pnt3D(40, 50, 2), 0.001));
Assert.IsTrue(GeomRelationUtils.PolycurveContainsXY(l, new Pnt3D(50, 50, 2), 0.001));
Assert.IsTrue(GeomRelationUtils.PolycurveContainsXY(l, new Pnt3D(100, 100, 2), 0.001));
Assert.IsTrue(GeomRelationUtils.PolycurveContainsXY(l, new Pnt3D(0, 50, 2), 0.001));
}
public void CanDetermineEqualMultipointXyVertical()
{
Multipoint <IPnt> mp1 = new Multipoint <IPnt>(
new[]
{
new Pnt3D(1234.5678, 9876.54321, 345.6),
new Pnt3D(1234.5678, 9876.54321, 355.6),
new Pnt3D(1234.5678, 9876.54321, 365.6),
new Pnt3D(1234.5678, 9876.54321, 375.6)
});
Multipoint <IPnt> mp2 = new Multipoint <IPnt>(mp1.GetPoints(0, null, true));
Assert.IsTrue(GeomRelationUtils.AreMultipointsEqualXY(mp1, mp2, 0.00001));
mp2 = new Multipoint <IPnt>(mp1.GetPoints(0, null, true).Reverse());
Assert.IsTrue(GeomRelationUtils.AreMultipointsEqualXY(mp1, mp2, 0.00001));
}
private static void AssignResultsToFootprintParts(
IEnumerable <RingGroup> resultRingGroups,
IDictionary <RingGroup, List <RingGroup> > splitPartsByFootprintPart,
double tolerance)
{
foreach (RingGroup resultPoly in resultRingGroups)
{
var assignmentCount = 0;
foreach (RingGroup cutFootprintPart in splitPartsByFootprintPart.Keys)
{
bool interiorIntersects;
if (resultPoly.IsVertical(tolerance))
{
// NOTE: Even though the second ring is vertical, do not use
// ring2CanHaveLinearSelfIntersections == true because a vertical
// ring running along the other ring's boundary should be interior
// intersecting and assigned.
interiorIntersects = GeomRelationUtils.InteriorIntersectXY(
cutFootprintPart.ExteriorRing, resultPoly.ExteriorRing,
tolerance, true);
}
else
{
interiorIntersects = GeomRelationUtils.InteriorIntersectXY(
cutFootprintPart, resultPoly, tolerance);
}
if (interiorIntersects)
{
assignmentCount++;
splitPartsByFootprintPart[cutFootprintPart].Add(resultPoly);
}
}
Assert.AreEqual(1, assignmentCount,
"Unexpected number of assignments to footprint parts.");
}
}
public void CanDetermineTouchesIdenticalRings()
{
var ring1 = new List <Pnt3D>();
// ring1: horizontal:
ring1.Add(new Pnt3D(0, 0, 9));
ring1.Add(new Pnt3D(0, 100, 9));
ring1.Add(new Pnt3D(100, 50, 9));
ring1.Add(new Pnt3D(100, 20, 9));
var ring2 = new List <Pnt3D>(ring1);
ring2.Add(ring2[0]);
for (var i = 0; i < 4; i++)
{
Pnt3D[] array1 = ring1.ToArray();
CollectionUtils.Rotate(array1, i);
var rotatedRing = new List <Pnt3D>(array1);
rotatedRing.Add((Pnt3D)rotatedRing[0].Clone());
var linestring1 = new Linestring(rotatedRing);
var linestring2 = new Linestring(ring2);
bool disjoint;
Assert.False(
GeomRelationUtils.TouchesXY(linestring1, linestring2, 0.0001, out disjoint));
Assert.False(disjoint);
linestring2.ReverseOrientation();
Assert.False(GeomRelationUtils.TouchesXY(linestring1, linestring2, 0.0001,
out disjoint,
true));
Assert.False(disjoint);
}
}
public void CanGetContainsXY()
{
var ring1 = new List <Pnt3D>();
ring1.Add(new Pnt3D(0, 0, 10));
ring1.Add(new Pnt3D(0, 100, 20));
ring1.Add(new Pnt3D(100, 100, 10));
ring1.Add(new Pnt3D(50, 50, 10));
ring1.Add(new Pnt3D(25, 75, 10));
ring1.Add(new Pnt3D(0, 0, 10));
Linestring l = new Linestring(ring1);
double tolerance = 0.001;
var pointOutsideRightOfCape = new Pnt3D(70, 50, 2);
Assert.IsFalse(
GeomRelationUtils.AreBoundsDisjoint(l, pointOutsideRightOfCape, tolerance));
Assert.IsFalse(GeomRelationUtils.LinesContainXY(l, pointOutsideRightOfCape, tolerance));
Assert.IsFalse(
GeomRelationUtils.PolycurveContainsXY(l, pointOutsideRightOfCape, tolerance));
Assert.AreEqual(
false, GeomRelationUtils.AreaContainsXY(l, pointOutsideRightOfCape, tolerance));
Assert.AreEqual(
false,
GeomRelationUtils.AreaContainsXY(l, pointOutsideRightOfCape, tolerance, true));
var pointOutsideLeftOfCape = new Pnt3D(40, 50, 2);
Assert.IsFalse(
GeomRelationUtils.AreBoundsDisjoint(l, pointOutsideLeftOfCape, tolerance));
Assert.IsFalse(GeomRelationUtils.LinesContainXY(l, pointOutsideLeftOfCape, tolerance));
Assert.IsFalse(
GeomRelationUtils.PolycurveContainsXY(l, pointOutsideLeftOfCape, tolerance));
Assert.AreEqual(
false, GeomRelationUtils.AreaContainsXY(l, pointOutsideLeftOfCape, tolerance));
Assert.AreEqual(
false,
GeomRelationUtils.AreaContainsXY(l, pointOutsideLeftOfCape, tolerance, true));
var pointInside = new Pnt3D(10, 90, 123);
Assert.IsFalse(GeomRelationUtils.AreBoundsDisjoint(l, pointInside, tolerance));
Assert.IsFalse(GeomRelationUtils.LinesContainXY(l, pointInside, tolerance));
Assert.IsTrue(GeomRelationUtils.PolycurveContainsXY(l, pointInside, tolerance));
Assert.AreEqual(true, GeomRelationUtils.AreaContainsXY(l, pointInside, tolerance));
Assert.AreEqual(
true, GeomRelationUtils.AreaContainsXY(l, pointInside, tolerance, true));
var pointOnTopRightCorner = new Pnt3D(100, 100, 2);
Assert.IsFalse(
GeomRelationUtils.AreBoundsDisjoint(l, pointOnTopRightCorner, tolerance));
Assert.IsTrue(GeomRelationUtils.LinesContainXY(l, pointOnTopRightCorner, tolerance));
Assert.IsTrue(
GeomRelationUtils.PolycurveContainsXY(l, pointOnTopRightCorner, tolerance));
Assert.AreEqual(
null, GeomRelationUtils.AreaContainsXY(l, pointOnTopRightCorner, tolerance));
Assert.AreEqual(
null, GeomRelationUtils.AreaContainsXY(l, pointOnTopRightCorner, tolerance, true));
var pointOnCape = new Pnt3D(50, 50, 321);
Assert.IsFalse(GeomRelationUtils.AreBoundsDisjoint(l, pointOnCape, tolerance));
Assert.IsTrue(GeomRelationUtils.LinesContainXY(l, pointOnCape, tolerance));
Assert.IsTrue(GeomRelationUtils.PolycurveContainsXY(l, pointOnCape, tolerance));
Assert.AreEqual(null, GeomRelationUtils.AreaContainsXY(l, pointOnCape, tolerance));
Assert.AreEqual(
null, GeomRelationUtils.AreaContainsXY(l, pointOnCape, tolerance, true));
var pointOnLeftEdge = new Pnt3D(0, 50, 2);
Assert.IsFalse(GeomRelationUtils.AreBoundsDisjoint(l, pointOnLeftEdge, tolerance));
Assert.IsTrue(GeomRelationUtils.LinesContainXY(l, pointOnLeftEdge, tolerance));
Assert.IsTrue(GeomRelationUtils.PolycurveContainsXY(l, pointOnLeftEdge, tolerance));
Assert.AreEqual(null, GeomRelationUtils.AreaContainsXY(l, pointOnLeftEdge, tolerance));
Assert.AreEqual(
null, GeomRelationUtils.AreaContainsXY(l, pointOnLeftEdge, tolerance, true));
}
