package org.apache.lucene.spatial;
import junit.framework.TestCase;
import org.apache.lucene.spatial.tier.InvalidGeoException;
/**
* Licensed to the Apache Software Foundation (ASF) under one or more
* contributor license agreements. See the NOTICE file distributed with
* this work for additional information regarding copyright ownership.
* The ASF licenses this file to You under the Apache License, Version 2.0
* (the "License"); you may not use this file except in compliance with
* the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
/**
*
*
**/
public class DistanceUtilsTest extends TestCase {
public void testBoxCorner() throws Exception {
double[] zero = new double[]{0, 0};
double[] zeroOne = new double[]{0, 1};
double[] oneOne = new double[]{1, 1};
double[] pt1 = new double[]{1.5, 110.3};
double[] result = DistanceUtils.vectorBoxCorner(zero, null, Math.sqrt(2), true);
assertEquals(1.0, result[0], 0);
assertEquals(1.0, result[1], 0);
result = DistanceUtils.vectorBoxCorner(zero, null, Math.sqrt(2), false);
assertEquals(-1.0, result[0], 0);
assertEquals(-1.0, result[1], 0);
result = DistanceUtils.vectorBoxCorner(oneOne, null, Math.sqrt(2), true);
assertEquals(2.0, result[0], 0);
assertEquals(2.0, result[1], 0);
result = DistanceUtils.vectorBoxCorner(zeroOne, null, Math.sqrt(2), true);
assertEquals(1.0, result[0], 0);
assertEquals(2.0, result[1], 0);
result = DistanceUtils.vectorBoxCorner(pt1, null, Math.sqrt(2), true);
assertEquals(2.5, result[0], 0.1);
assertEquals(111.3, result[1], 0.1);
result = DistanceUtils.vectorBoxCorner(pt1, null, Math.sqrt(2), false);
assertEquals(0.5, result[0], 0.1);
assertEquals(109.3, result[1], 0.1);
}
public void testNormLatLon() throws Exception {
}
public void testLatLonCorner() throws Exception {
double[] zero = new double[]{0, 0};
double[] zero45 = new double[]{0, DistanceUtils.DEG_45_AS_RADS};
double[] result;
// 00°38′09″N, 000°38′09″E
//Verify at http://www.movable-type.co.uk/scripts/latlong.html
result = DistanceUtils.latLonCorner(zero[0], zero[1], 100, null, true, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(0.63583 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(0.63583 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
result = DistanceUtils.latLonCornerDegs(zero[0], zero[1], 100, null, true, DistanceUtils.EARTH_MEAN_RADIUS_KM);
// 00°38′09″N, 000°38′09″E
assertEquals(0.63583, result[0], 0.001);
assertEquals(0.63583, result[1], 0.001);
result = DistanceUtils.latLonCornerDegs(zero[0], zero[1], 100, null, false, DistanceUtils.EARTH_MEAN_RADIUS_KM);
// 00°38′09″N, 000°38′09″E
assertEquals(-0.63583, result[0], 0.001);
assertEquals(-0.63583, result[1], 0.001);
//test some edge cases
//89°16′02″N, 060°12′35″E
result = DistanceUtils.latLonCornerDegs(89.0, 0, 100, null, true, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(89.26722, result[0], 0.001);
assertEquals(60.20972, result[1], 0.001);
result = DistanceUtils.latLonCornerDegs(0, -179.0, 100, null, true, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(0.63583, result[0], 0.001);
assertEquals(-178.36417, result[1], 0.001);
}
public void testPointBearing() throws Exception {
double[] zero = new double[]{0, 0};
double[] zero45 = new double[]{40 * DistanceUtils.DEGREES_TO_RADIANS, DistanceUtils.DEG_45_AS_RADS};
double[] result;
// 00°38′09″N, 000°38′09″E
//Verify at http://www.movable-type.co.uk/scripts/latlong.html
result = DistanceUtils.pointOnBearing(zero[0], zero[1], 100, DistanceUtils.DEG_45_AS_RADS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(0.63583 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(0.63583 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
//should be above the current point at 0.8994°,0.0000°
result = DistanceUtils.pointOnBearing(zero[0], zero[1], 100, 0, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(0.8994 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(0, result[1], 0.001);
//directly below
result = DistanceUtils.pointOnBearing(zero[0], zero[1], 100, DistanceUtils.DEG_180_AS_RADS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(-0.8994 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(0, result[1], 0.001);
//0.7183°,0.5414° -- 37 deg bearing
result = DistanceUtils.pointOnBearing(zero[0], zero[1], 100, 37 * DistanceUtils.DEGREES_TO_RADIANS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(0.7183 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(0.5414 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
result = DistanceUtils.pointOnBearing(zero45[0], zero45[1], 100, DistanceUtils.DEG_45_AS_RADS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
//40.6328°,45.8381°
assertEquals(40.6328 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(45.8381 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
result = DistanceUtils.pointOnBearing(1 * DistanceUtils.DEGREES_TO_RADIANS, 1 * DistanceUtils.DEGREES_TO_RADIANS, 100, DistanceUtils.DEG_90_AS_RADS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
//0.9997°,1.8994°
assertEquals(0.9997 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(1.8994 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
result = DistanceUtils.pointOnBearing(-10 * DistanceUtils.DEGREES_TO_RADIANS, -150 * DistanceUtils.DEGREES_TO_RADIANS, 15, 205*DistanceUtils.DEGREES_TO_RADIANS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
//-10.1222°,-150.0578°
assertEquals(-10.1222 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(-150.0578 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
result = DistanceUtils.pointOnBearing(-10 * DistanceUtils.DEGREES_TO_RADIANS, -150 * DistanceUtils.DEGREES_TO_RADIANS, 200, 63*DistanceUtils.DEGREES_TO_RADIANS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
//-9.1797°,-148.3767°
assertEquals(-9.1797 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(-148.3767 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
result = DistanceUtils.pointOnBearing(-10 * DistanceUtils.DEGREES_TO_RADIANS, -150 * DistanceUtils.DEGREES_TO_RADIANS, 3000, 63*DistanceUtils.DEGREES_TO_RADIANS, null, DistanceUtils.EARTH_MEAN_RADIUS_KM);
//2.7561°,-126.1281°
assertEquals(2.7561 * DistanceUtils.DEGREES_TO_RADIANS, result[0], 0.001);
assertEquals(-126.1281 * DistanceUtils.DEGREES_TO_RADIANS, result[1], 0.001);
}
public void testVectorDistance() throws Exception {
double[] zero = new double[]{0, 0};
double[] zeroOne = new double[]{0, 1};
double[] oneZero = new double[]{1, 0};
double[] oneOne = new double[]{1, 1};
double distance;
distance = DistanceUtils.vectorDistance(zero, zeroOne, 2);
assertEquals(1.0, distance, 0);
distance = DistanceUtils.vectorDistance(zero, oneZero, 2);
assertEquals(1.0, distance, 0);
distance = DistanceUtils.vectorDistance(zero, oneOne, 2);
assertEquals(Math.sqrt(2), distance, 0.001);
distance = DistanceUtils.squaredEuclideanDistance(zero, oneOne);
assertEquals(2, distance, 0.001);
}
public void testHaversine() throws Exception {
double distance;
//compare to http://www.movable-type.co.uk/scripts/latlong.html
distance = DistanceUtils.haversine(0, 0, Math.PI / 4.0, Math.PI / 4.0, DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(6672.0, distance, 0.5);
distance = DistanceUtils.haversine(0, 0, Math.toRadians(20), Math.toRadians(20), DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(3112, distance, 0.5);
distance = DistanceUtils.haversine(0, 0, Math.toRadians(1), Math.toRadians(1), DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(157.2, distance, 0.5);
//Try some around stuff
distance = DistanceUtils.haversine(Math.toRadians(1), Math.toRadians(-1),
Math.toRadians(1), Math.toRadians(1), DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(222.4, distance, 0.5);
distance = DistanceUtils.haversine(Math.toRadians(89), Math.toRadians(-1),
Math.toRadians(89), Math.toRadians(179), DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(222.4, distance, 0.5);
distance = DistanceUtils.haversine(Math.toRadians(89), Math.toRadians(-1),
Math.toRadians(49), Math.toRadians(179), DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(4670, distance, 0.5);
distance = DistanceUtils.haversine(Math.toRadians(0), Math.toRadians(-179),
Math.toRadians(0), Math.toRadians(179), DistanceUtils.EARTH_MEAN_RADIUS_KM);
assertEquals(222.4, distance, 0.5);
}
public void testParse() throws Exception {
String[] parse;
parse = DistanceUtils.parsePoint(null, "89.0,73.2", 2);
assertEquals(2, parse.length);
assertEquals("89.0", parse[0]);
assertEquals("73.2", parse[1]);
parse = DistanceUtils.parsePoint(null, "89.0,73.2,-92.3", 3);
assertEquals(3, parse.length);
assertEquals("89.0", parse[0]);
assertEquals("73.2", parse[1]);
assertEquals("-92.3", parse[2]);
parse = DistanceUtils.parsePoint(null, " 89.0 , 73.2 , -92.3 ", 3);
assertEquals(3, parse.length);
assertEquals("89.0", parse[0]);
assertEquals("73.2", parse[1]);
assertEquals("-92.3", parse[2]);
String[] foo = DistanceUtils.parsePoint(parse, "89.0 , 73.2 , -92.3", 3);
//should be same piece of memory
assertTrue(foo == parse);
assertEquals(3, parse.length);
assertEquals("89.0", parse[0]);
assertEquals("73.2", parse[1]);
assertEquals("-92.3", parse[2]);
//array should get automatically resized
parse = DistanceUtils.parsePoint(new String[1], "89.0 , 73.2 , -92.3", 3);
assertEquals(3, parse.length);
assertEquals("89.0", parse[0]);
assertEquals("73.2", parse[1]);
assertEquals("-92.3", parse[2]);
try {
parse = DistanceUtils.parsePoint(null, "89.0 , ", 3);
assertTrue(false);
} catch (InvalidGeoException e) {
}
try {
parse = DistanceUtils.parsePoint(null, " , 89.0 ", 3);
assertTrue(false);
} catch (InvalidGeoException e) {
}
try {
parse = DistanceUtils.parsePoint(null, "", 3);
assertTrue(false);
} catch (InvalidGeoException e) {
}
double[] dbls = DistanceUtils.parsePointDouble(null, "89.0 , 73.2 , -92.3", 3);
assertEquals(3, dbls.length);
assertEquals(89.0, dbls[0], 0);
assertEquals(73.2, dbls[1], 0.1);
assertEquals(-92.3, dbls[2], 0.1);
try {
dbls = DistanceUtils.parsePointDouble(null, "89.0 , foo , -92.3", 3);
assertTrue(false);
} catch (NumberFormatException e) {
}
dbls = DistanceUtils.parseLatitudeLongitude(null, "89.0 , 73.2 ");
assertEquals(2, dbls.length);
assertEquals(89.0, dbls[0], 0.1);
assertEquals(73.2, dbls[1], 0.1);
//test some bad lat/long pairs
try {
dbls = DistanceUtils.parseLatitudeLongitude(null, "189.0 , 73.2 ");
assertTrue(false);
} catch (InvalidGeoException e) {
}
try {
dbls = DistanceUtils.parseLatitudeLongitude(null, "89.0 , 273.2 ");
assertTrue(false);
} catch (InvalidGeoException e) {
}
}
}