/*
* Copyright (c) 2016 Martin Davis.
*
* All rights reserved. This program and the accompanying materials
* are made available under the terms of the Eclipse Public License v1.0
* and Eclipse Distribution License v. 1.0 which accompanies this distribution.
* The Eclipse Public License is available at http://www.eclipse.org/legal/epl-v10.html
* and the Eclipse Distribution License is available at
*
* http://www.eclipse.org/org/documents/edl-v10.php.
*/
package org.locationtech.jts.operation.distance3d;
import org.locationtech.jts.geom.Geometry;
import org.locationtech.jts.geom.GeometryFactory;
import org.locationtech.jts.io.ParseException;
import org.locationtech.jts.io.WKTReader;
import junit.framework.TestCase;
import junit.textui.TestRunner;
public class WithinDistance3DTest extends TestCase
{
static GeometryFactory geomFact = new GeometryFactory();
static WKTReader rdr = new WKTReader();
public static void main(String args[]) {
TestRunner.run(WithinDistance3DTest.class);
}
public WithinDistance3DTest(String name) {
super(name);
}
public void testEmpty()
{
checkWithinDistance( "POINT EMPTY", "POINT EMPTY", 0);
checkWithinDistance( "LINESTRING EMPTY", "POINT (0 0 0)", 1, true);
}
public void testPointPoint() {
checkWithinDistance( "POINT (0 0 0 )",
"POINT (0 0 1 )",
1);
checkWithinDistance( "POINT (0 0 0 )",
"POINT (0 0 1 )",
0.5, false);
checkWithinDistance( "POINT (10 10 1 )",
"POINT (11 11 2 )",
1.733);
checkWithinDistance( "POINT (10 10 0 )",
"POINT (10 20 10 )",
14.143);
}
public void testPointSeg() {
checkWithinDistance( "LINESTRING (0 0 0, 10 10 10 )",
"POINT (5 5 5 )",
0);
checkWithinDistance( "LINESTRING (10 10 10, 20 20 20 )",
"POINT (11 11 10 )",
0.8, false);
}
public void testCrossSegmentsFlat() {
checkWithinDistance( "LINESTRING (0 0 0, 10 10 0 )",
"LINESTRING (10 0 0, 0 10 0 )",
0);
checkWithinDistance( "LINESTRING (0 0 10, 30 10 10 )",
"LINESTRING (10 0 10, 0 10 10 )",
0);
}
public void testCrossSegments() {
checkWithinDistance( "LINESTRING (0 0 0, 10 10 0 )",
"LINESTRING (10 0 1, 0 10 1 )",
1);
checkWithinDistance( "LINESTRING (0 0 0, 20 20 0 )",
"LINESTRING (10 0 1, 0 10 1 )",
1);
checkWithinDistance( "LINESTRING (20 10 20, 10 20 10 )",
"LINESTRING (10 10 20, 20 20 10 )",
0);
}
public void testTSegmentsFlat() {
checkWithinDistance( "LINESTRING (10 10 0, 10 20 0 )",
"LINESTRING (20 15 0, 25 15 0 )",
10);
}
public void testParallelSegmentsFlat() {
checkWithinDistance( "LINESTRING (10 10 0, 20 20 0 )",
"LINESTRING (10 20 0, 20 30 0 )",
7.0710678118654755);
}
public void testParallelSegments() {
checkWithinDistance( "LINESTRING (0 0 0, 1 0 0 )",
"LINESTRING (0 0 1, 1 0 1 )",
1);
checkWithinDistance( "LINESTRING (10 10 0, 20 10 0 )",
"LINESTRING (10 20 10, 20 20 10 )",
14.142135623730951);
checkWithinDistance( "LINESTRING (10 10 0, 20 20 0 )",
"LINESTRING (10 20 10, 20 30 10 )",
12.24744871391589 );
// = distance from LINESTRING (10 10 0, 20 20 0 ) to POINT(10 20 10)
// = hypotenuse(7.0710678118654755, 10)
}
public void testLineLine()
{
checkWithinDistance( "LINESTRING (0 1 2, 1 1 1, 1 0 2 )",
"LINESTRING (0 0 0.1, .5 .5 0, 1 1 0, 1.5 1.5 0, 2 2 0 )",
1);
checkWithinDistance( "LINESTRING (10 10 20, 20 20 30, 20 20 1, 30 30 5 )",
"LINESTRING (1 80 10, 0 39 5, 39 0 5, 80 1 20)",
0.7071067811865476);
}
public void testPointPolygon()
{
// point above poly
checkWithinDistance( "POINT (150 150 10)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
10);
// point below poly
checkWithinDistance( "POINT (150 150 -10)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
10);
// point right of poly in YZ plane
checkWithinDistance( "POINT (10 150 150)",
"POLYGON ((0 100 200, 0 200 200, 0 200 100, 0 100 100, 0 100 200))",
10);
}
public void testPointPolygonFlat()
{
// inside
checkWithinDistance( "POINT (150 150 0)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
0);
// outside
checkWithinDistance( "POINT (250 250 0)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
70.71067811865476);
// on
checkWithinDistance( "POINT (200 200 0)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
0);
}
public void testLinePolygonFlat()
{
// line inside
checkWithinDistance( "LINESTRING (150 150 0, 160 160 0)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
0);
// line outside
checkWithinDistance( "LINESTRING (200 250 0, 260 260 0)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
50);
// line touching
checkWithinDistance( "LINESTRING (200 200 0, 260 260 0)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
0);
}
public void testLinePolygonSimple()
{
// line crossing inside
checkWithinDistance( "LINESTRING (150 150 10, 150 150 -10)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
0);
// vertical line above
checkWithinDistance( "LINESTRING (200 200 10, 260 260 100)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
10);
// vertical line touching
checkWithinDistance( "LINESTRING (200 200 0, 260 260 100)",
"POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0))",
0);
}
String polyHoleFlat = "POLYGON ((100 200 0, 200 200 0, 200 100 0, 100 100 0, 100 200 0), (120 180 0, 180 180 0, 180 120 0, 120 120 0, 120 180 0))";
public void testLinePolygonHoleFlat()
{
// line crossing hole
checkWithinDistance( "LINESTRING (150 150 10, 150 150 -10)", polyHoleFlat, 20, false);
// line crossing interior
checkWithinDistance( "LINESTRING (110 110 10, 110 110 -10)", polyHoleFlat, 0);
}
public void testPointPolygonHoleFlat()
{
// point above poly hole
checkWithinDistance( "POINT (130 130 10)", polyHoleFlat, 14.143);
// point below poly hole
checkWithinDistance( "POINT (130 130 -10)", polyHoleFlat, 14.143);
// point above poly
checkWithinDistance( "POINT (110 110 100)", polyHoleFlat, 100);
}
public void testMultiPoint()
{
checkWithinDistance(
"MULTIPOINT ((0 0 0), (0 0 100), (100 100 100))",
"MULTIPOINT ((100 100 99), (50 50 50), (25 100 33))",
1
);
}
public void testMultiLineString()
{
checkWithinDistance(
"MULTILINESTRING ((0 0 0, 10 10 10), (0 0 100, 25 25 25, 40 40 50), (100 100 100, 100 101 102))",
"MULTILINESTRING ((100 100 99, 100 100 99), (100 102 102, 200 200 20), (25 100 33, 25 100 35))",
1
);
}
public void testMultiPolygon()
{
checkWithinDistance(
// Polygons parallel to XZ plane
"MULTIPOLYGON ( ((120 120 -10, 120 120 100, 180 120 100, 180 120 -10, 120 120 -10)), ((120 200 -10, 120 200 190, 180 200 190, 180 200 -10, 120 200 -10)) )",
// Polygons parallel to XY plane
"MULTIPOLYGON ( ((100 200 200, 200 200 200, 200 100 200, 100 100 200, 100 200 200)), ((100 200 210, 200 200 210, 200 100 210, 100 100 210, 100 200 210)) )",
10
);
}
//==========================================================
// Convenience methods
//==========================================================
private void checkWithinDistance(String wkt1, String wkt2, double distance)
{
checkWithinDistance(wkt1, wkt2, distance, true);
}
private void checkWithinDistance(String wkt1, String wkt2, double distance, boolean expectedResult)
{
Geometry g1;
Geometry g2;
try {
g1 = rdr.read(wkt1);
} catch (ParseException e) {
throw new RuntimeException(e.toString());
}
try {
g2 = rdr.read(wkt2);
} catch (ParseException e) {
throw new RuntimeException(e.toString());
}
// check both orders for arguments
checkWithinDistance(g1, g2, distance, expectedResult);
checkWithinDistance(g2, g1, distance, expectedResult);
}
private void checkWithinDistance(Geometry g1, Geometry g2, double distance, boolean expectedResult)
{
boolean isWithinDist = Distance3DOp.isWithinDistance(g1, g2, distance);
assertEquals(expectedResult, isWithinDist);
}
}