| public MakeCircle ( System.Point point, double distance ) : Circle | ||
| point | System.Point | |
| distance | double | |
| Результат | Circle | |
public static void testCircleReset(SpatialContext ctx)
{
Circle c = ctx.MakeCircle(3, 4, 5);
Circle c2 = ctx.MakeCircle(5, 6, 7);
c2.Reset(3, 4, 5); // to c1
Assert.Equal(c, c2);
Assert.Equal(c.GetBoundingBox(), c2.GetBoundingBox());
}
public void TestSimpleCircle(SpatialContext ctx)
{
base.ctx = ctx;
double[] theXs = new double[] { -10, 0, 10 };
foreach (double x in theXs)
{
double[] theYs = new double[] { -20, 0, 20 };
foreach (double y in theYs)
{
TestCircle(x, y, 0);
TestCircle(x, y, 5);
}
}
//INTERSECTION:
//Start with some static tests that have shown to cause failures at some point:
Assert.Equal( /*"getX not getY",*/
SpatialRelation.INTERSECTS,
ctx.MakeCircle(107, -81, 147).Relate(ctx.MakeRect(92, 121, -89, 74), ctx));
TestCircleIntersect();
}
public void TestGeoCircle(SpatialContext ctx)
{
base.ctx = ctx;
Assert.Equal("Circle(Pt(x=10.0,y=20.0), d=30.0° 3335.85km)", ctx.MakeCircle(10, 20, 30).ToString());
double v = 200 * (random.NextDouble() > 0.5 ? -1 : 1);
Assert.Throws<InvalidShapeException>(() => ctx.MakeCircle(v,0,5));
Assert.Throws<InvalidShapeException>(() => ctx.MakeCircle(0, v, 5));
Assert.Throws<InvalidShapeException>(() => ctx.MakeCircle(random.Next(-180, 180), random.Next(-90, 90), v));
//--Start with some static tests that once failed:
//Bug: numeric edge at pole, fails to init
ctx.MakeCircle(110, -12, 90 + 12);
//Bug: horizXAxis not in enclosing rectangle, assertion
ctx.MakeCircle(-44, 16, 106);
ctx.MakeCircle(-36, -76, 14);
ctx.MakeCircle(107, 82, 172);
// TODO need to update this test to be valid
//{
// //Bug in which distance was being confused as being in the same coordinate system as x,y.
// double distDeltaToPole = 0.001;//1m
// double distDeltaToPoleDEG = ctx.getDistCalc().distanceToDegrees(distDeltaToPole);
// double dist = 1;//1km
// double distDEG = ctx.getDistCalc().distanceToDegrees(dist);
// Circle c = ctx.makeCircle(0, 90 - distDeltaToPoleDEG - distDEG, dist);
// Rectangle cBBox = c.getBoundingBox();
// Rectangle r = ctx.makeRect(cBBox.getMaxX() * 0.99, cBBox.getMaxX() + 1, c.getCenter().getY(), c.getCenter().getY());
// assertEquals(INTERSECTS, c.getBoundingBox().relate(r, ctx));
// assertEquals("dist != xy space", INTERSECTS, c.relate(r, ctx));//once failed here
//}
assertEquals("edge rounding issue 2", SpatialRelation.INTERSECTS, ctx.MakeCircle(84, -40, 136).Relate(ctx.MakeRectangle(-150, -80, 34, 84)));
assertEquals("edge rounding issue", SpatialRelation.CONTAINS, ctx.MakeCircle(0, 66, 156).Relate(ctx.MakePoint(0, -90)));
assertEquals("nudge back circle", SpatialRelation.CONTAINS, ctx.MakeCircle(-150, -90, 122).Relate(ctx.MakeRectangle(0, -132, 32, 32)));
assertEquals("wrong estimate", SpatialRelation.DISJOINT, ctx.MakeCircle(-166, 59, kmToDeg(5226.2)).Relate(ctx.MakeRectangle(36, 66, 23, 23)));
assertEquals("bad CONTAINS (dateline)", SpatialRelation.INTERSECTS, ctx.MakeCircle(56, -50, kmToDeg(12231.5)).Relate(ctx.MakeRectangle(108, 26, 39, 48)));
assertEquals("bad CONTAINS (backwrap2)", SpatialRelation.INTERSECTS,
ctx.MakeCircle(112, -3, 91).Relate(ctx.MakeRectangle(-163, 29, -38, 10)));
assertEquals("bad CONTAINS (r x-wrap)", SpatialRelation.INTERSECTS,
ctx.MakeCircle(-139, 47, 80).Relate(ctx.MakeRectangle(-180, 180, -3, 12)));
assertEquals("bad CONTAINS (pwrap)", SpatialRelation.INTERSECTS,
ctx.MakeCircle(-139, 47, 80).Relate(ctx.MakeRectangle(-180, 179, -3, 12)));
assertEquals("no-dist 1", SpatialRelation.WITHIN,
ctx.MakeCircle(135, 21, 0).Relate(ctx.MakeRectangle(-103, -154, -47, 52)));
assertEquals("bbox <= >= -90 bug", SpatialRelation.CONTAINS,
ctx.MakeCircle(-64, -84, 124).Relate(ctx.MakeRectangle(-96, 96, -10, -10)));
//The horizontal axis line of a geo circle doesn't necessarily pass through c's ctr.
assertEquals("c's horiz axis doesn't pass through ctr", SpatialRelation.INTERSECTS,
ctx.MakeCircle(71, -44, 40).Relate(ctx.MakeRectangle(15, 27, -62, -34)));
assertEquals("pole boundary", SpatialRelation.INTERSECTS,
ctx.MakeCircle(-100, -12, 102).Relate(ctx.MakeRectangle(143, 175, 4, 32)));
assertEquals("full circle assert", SpatialRelation.CONTAINS,
ctx.MakeCircle(-64, 32, 180).Relate(ctx.MakeRectangle(47, 47, -14, 90)));
//--Now proceed with systematic testing:
assertEquals(ctx.GetWorldBounds(), ctx.MakeCircle(0, 0, 180).GetBoundingBox());
//assertEquals(ctx.makeCircle(0,0,distToOpposeSide/2 - 500).getBoundingBox());
double[] theXs = new double[] { -180, -45, 90 };
foreach (double x in theXs)
{
double[] theYs = new double[] { -90, -45, 0, 45, 90 };
foreach (double y in theYs)
{
TestCircle(x, y, 0);
TestCircle(x, y, kmToDeg(500));
TestCircle(x, y, 90);
TestCircle(x, y, 180);
}
}
TestCircleIntersect();
}
public void TestGeoCircle(SpatialContext ctx)
{
base.ctx = ctx;
//--Start with some static tests that once failed:
//Bug: numeric edge at pole, fails to init
ctx.MakeCircle(110, -12, ctx.GetDistCalc().DegreesToDistance(90 + 12));
//Bug: horizXAxis not in enclosing rectangle, assertion
ctx.MakeCircle(-44, 16, DegToDist(106));
ctx.MakeCircle(-36, -76, DegToDist(14));
ctx.MakeCircle(107, 82, DegToDist(172));
// TODO need to update this test to be valid
//{
// //Bug in which distance was being confused as being in the same coordinate system as x,y.
// double distDeltaToPole = 0.001;//1m
// double distDeltaToPoleDEG = ctx.getDistCalc().distanceToDegrees(distDeltaToPole);
// double dist = 1;//1km
// double distDEG = ctx.getDistCalc().distanceToDegrees(dist);
// Circle c = ctx.makeCircle(0, 90 - distDeltaToPoleDEG - distDEG, dist);
// Rectangle cBBox = c.getBoundingBox();
// Rectangle r = ctx.makeRect(cBBox.getMaxX() * 0.99, cBBox.getMaxX() + 1, c.getCenter().getY(), c.getCenter().getY());
// assertEquals(INTERSECTS, c.getBoundingBox().relate(r, ctx));
// assertEquals("dist != xy space", INTERSECTS, c.relate(r, ctx));//once failed here
//}
Assert.Equal(/*"nudge back circle", */ SpatialRelation.CONTAINS, ctx.MakeCircle(-150, -90, DegToDist(122)).Relate(ctx.MakeRect(0, -132, 32, 32), ctx));
Assert.Equal(/* "wrong estimate", */ SpatialRelation.DISJOINT, ctx.MakeCircle(-166, 59, 5226.2).Relate(ctx.MakeRect(36, 66, 23, 23), ctx));
Assert.Equal(/*"bad CONTAINS (dateline)",*/ SpatialRelation.INTERSECTS,
ctx.MakeCircle(56, -50, 12231.5).Relate(ctx.MakeRect(108, 26, 39, 48), ctx));
Assert.Equal(/*"bad CONTAINS (backwrap2)",*/ SpatialRelation.INTERSECTS,
ctx.MakeCircle(112, -3, DegToDist(91)).Relate(ctx.MakeRect(-163, 29, -38, 10), ctx));
Assert.Equal(/*"bad CONTAINS (r x-wrap)",*/ SpatialRelation.INTERSECTS,
ctx.MakeCircle(-139, 47, DegToDist(80)).Relate(ctx.MakeRect(-180, 180, -3, 12), ctx));
Assert.Equal(/*"bad CONTAINS (pwrap)",*/ SpatialRelation.INTERSECTS,
ctx.MakeCircle(-139, 47, DegToDist(80)).Relate(ctx.MakeRect(-180, 179, -3, 12), ctx));
Assert.Equal(/*"no-dist 1",*/ SpatialRelation.WITHIN,
ctx.MakeCircle(135, 21, 0).Relate(ctx.MakeRect(-103, -154, -47, 52), ctx));
Assert.Equal(/*"bbox <= >= -90 bug",*/ SpatialRelation.CONTAINS,
ctx.MakeCircle(-64, -84, DegToDist(124)).Relate(ctx.MakeRect(-96, 96, -10, -10), ctx));
//The horizontal axis line of a geo circle doesn't necessarily pass through c's ctr.
Assert.Equal(/*"c's horiz axis doesn't pass through ctr",*/ SpatialRelation.INTERSECTS,
ctx.MakeCircle(71, -44, DegToDist(40)).Relate(ctx.MakeRect(15, 27, -62, -34), ctx));
Assert.Equal(/*"pole boundary",*/ SpatialRelation.INTERSECTS,
ctx.MakeCircle(-100, -12, DegToDist(102)).Relate(ctx.MakeRect(143, 175, 4, 32), ctx));
Assert.Equal(/*"full circle assert",*/ SpatialRelation.CONTAINS,
ctx.MakeCircle(-64, 32, DegToDist(180)).Relate(ctx.MakeRect(47, 47, -14, 90), ctx));
//--Now proceed with systematic testing:
double distToOpposeSide = ctx.GetUnits().EarthRadius() * Math.PI;
Assert.Equal(ctx.GetWorldBounds(), ctx.MakeCircle(0, 0, distToOpposeSide).GetBoundingBox());
//assertEquals(ctx.makeCircle(0,0,distToOpposeSide/2 - 500).getBoundingBox());
double[] theXs = new double[] { -180, -45, 90 };
foreach (double x in theXs)
{
double[] theYs = new double[] { -90, -45, 0, 45, 90 };
foreach (double y in theYs)
{
TestCircle(x, y, 0);
TestCircle(x, y, 500);
TestCircle(x, y, DegToDist(90));
TestCircle(x, y, ctx.GetUnits().EarthRadius() * 6);
}
}
TestCircleIntersect();
}
