294 lines
8.1 KiB
JavaScript
294 lines
8.1 KiB
JavaScript
// See http://www.movable-type.co.uk/scripts/latlong.html
|
|
// FIXME add tests for crossTrackDistance
|
|
// FIXME add tests for maximumLatitude
|
|
// FIXME add tests for offset
|
|
|
|
goog.provide('ol.test.Sphere');
|
|
|
|
|
|
describe('ol.Sphere', function() {
|
|
|
|
var sphere = new ol.Sphere(6371);
|
|
var expected = [
|
|
{
|
|
c1: [0, 0],
|
|
c2: [0, 0],
|
|
cosineDistance: 0,
|
|
equirectangularDistance: 0,
|
|
finalBearing: 180,
|
|
haversineDistance: 0,
|
|
initialBearing: 0,
|
|
midpoint: [0, 0]
|
|
},
|
|
{
|
|
c1: [0, 0],
|
|
c2: [45, 45],
|
|
cosineDistance: 6671.695598673524,
|
|
equirectangularDistance: 6812.398372654371,
|
|
finalBearing: 54.735610317245346,
|
|
haversineDistance: 6671.695598673525,
|
|
initialBearing: 35.264389682754654,
|
|
midpoint: [18.434948822922006, 24.0948425521107]
|
|
},
|
|
{
|
|
c1: [0, 0],
|
|
c2: [-45, 45],
|
|
cosineDistance: 6671.695598673524,
|
|
equirectangularDistance: 6812.398372654371,
|
|
finalBearing: 305.26438968275465,
|
|
haversineDistance: 6671.695598673525,
|
|
initialBearing: -35.264389682754654,
|
|
midpoint: [-18.434948822922006, 24.0948425521107]
|
|
},
|
|
{
|
|
c1: [0, 0],
|
|
c2: [-45, -45],
|
|
cosineDistance: 6671.695598673524,
|
|
equirectangularDistance: 6812.398372654371,
|
|
finalBearing: 234.73561031724535,
|
|
haversineDistance: 6671.695598673525,
|
|
initialBearing: -144.73561031724535,
|
|
midpoint: [-18.434948822922006, -24.0948425521107]
|
|
},
|
|
{
|
|
c1: [0, 0],
|
|
c2: [45, -45],
|
|
cosineDistance: 6671.695598673524,
|
|
equirectangularDistance: 6812.398372654371,
|
|
finalBearing: 125.26438968275465,
|
|
haversineDistance: 6671.695598673525,
|
|
initialBearing: 144.73561031724535,
|
|
midpoint: [18.434948822922006, -24.0948425521107]
|
|
},
|
|
{
|
|
c1: [45, 45],
|
|
c2: [45, 45],
|
|
cosineDistance: 0,
|
|
equirectangularDistance: 0,
|
|
finalBearing: 180,
|
|
haversineDistance: 0,
|
|
initialBearing: 0,
|
|
midpoint: [45.00000000000005, 45]
|
|
},
|
|
{
|
|
c1: [45, 45],
|
|
c2: [-45, 45],
|
|
cosineDistance: 6671.695598673525,
|
|
equirectangularDistance: 7076.401799751738,
|
|
finalBearing: 234.73561031724535,
|
|
haversineDistance: 6671.695598673525,
|
|
initialBearing: -54.73561031724535,
|
|
midpoint: [0, 54.735610317245346]
|
|
},
|
|
{
|
|
c1: [45, 45],
|
|
c2: [-45, -45],
|
|
cosineDistance: 13343.391197347048,
|
|
equirectangularDistance: 14152.803599503475,
|
|
finalBearing: 234.73561031724535,
|
|
haversineDistance: 13343.391197347048,
|
|
initialBearing: -125.26438968275465,
|
|
midpoint: [0, 0]
|
|
},
|
|
{
|
|
c1: [45, 45],
|
|
c2: [45, -45],
|
|
cosineDistance: 10007.543398010284,
|
|
equirectangularDistance: 10007.543398010286,
|
|
finalBearing: 180,
|
|
haversineDistance: 10007.543398010286,
|
|
initialBearing: 180,
|
|
midpoint: [45.00000000000005, 0]
|
|
},
|
|
{
|
|
c1: [-45, 45],
|
|
c2: [-45, 45],
|
|
cosineDistance: 0,
|
|
equirectangularDistance: 0,
|
|
finalBearing: 180,
|
|
haversineDistance: 0,
|
|
initialBearing: 0,
|
|
midpoint: [-45.00000000000005, 45]
|
|
},
|
|
{
|
|
c1: [-45, 45],
|
|
c2: [-45, -45],
|
|
cosineDistance: 10007.543398010284,
|
|
equirectangularDistance: 10007.543398010286,
|
|
finalBearing: 180,
|
|
haversineDistance: 10007.543398010286,
|
|
initialBearing: 180,
|
|
midpoint: [-45.00000000000005, 0]
|
|
},
|
|
{
|
|
c1: [-45, 45],
|
|
c2: [45, -45],
|
|
cosineDistance: 13343.391197347048,
|
|
equirectangularDistance: 14152.803599503475,
|
|
finalBearing: 125.26438968275465,
|
|
haversineDistance: 13343.391197347048,
|
|
initialBearing: 125.26438968275465,
|
|
midpoint: [0, 0]
|
|
},
|
|
{
|
|
c1: [-45, -45],
|
|
c2: [-45, -45],
|
|
cosineDistance: 0,
|
|
equirectangularDistance: 0,
|
|
finalBearing: 180,
|
|
haversineDistance: 0,
|
|
initialBearing: 0,
|
|
midpoint: [-45.00000000000005, -45]
|
|
},
|
|
{
|
|
c1: [-45, -45],
|
|
c2: [45, -45],
|
|
cosineDistance: 6671.695598673525,
|
|
equirectangularDistance: 7076.401799751738,
|
|
finalBearing: 54.735610317245346,
|
|
haversineDistance: 6671.695598673525,
|
|
initialBearing: 125.26438968275465,
|
|
midpoint: [0, -54.735610317245346]
|
|
},
|
|
{
|
|
c1: [45, -45],
|
|
c2: [45, -45],
|
|
cosineDistance: 0,
|
|
equirectangularDistance: 0,
|
|
finalBearing: 180,
|
|
haversineDistance: 0,
|
|
initialBearing: 0,
|
|
midpoint: [45.00000000000005, -45]
|
|
}
|
|
];
|
|
|
|
describe('cosineDistance', function() {
|
|
|
|
it('results match Chris Veness\'s reference implementation', function() {
|
|
var e, i;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
expect(sphere.cosineDistance(e.c1, e.c2)).to.roughlyEqual(
|
|
e.cosineDistance, 1e-3);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('equirectangularDistance', function() {
|
|
|
|
it('results match Chris Veness\'s reference implementation', function() {
|
|
var e, i;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
expect(sphere.equirectangularDistance(e.c1, e.c2)).to.roughlyEqual(
|
|
e.equirectangularDistance, 1e-9);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('finalBearing', function() {
|
|
|
|
it('results match Chris Veness\'s reference implementation', function() {
|
|
var e, i;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
expect(sphere.finalBearing(e.c1, e.c2)).to.roughlyEqual(
|
|
e.finalBearing, 1e-9);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('haversineDistance', function() {
|
|
|
|
it('results match Chris Veness\'s reference implementation', function() {
|
|
var e, i;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
expect(sphere.haversineDistance(e.c1, e.c2)).to.roughlyEqual(
|
|
e.haversineDistance, 1e-9);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('initialBearing', function() {
|
|
|
|
it('results match Chris Veness\'s reference implementation', function() {
|
|
var e, i;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
expect(sphere.initialBearing(e.c1, e.c2)).to.roughlyEqual(
|
|
e.initialBearing, 1e-9);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('interpolate', function() {
|
|
|
|
it('results match at the start, midpoint, and end', function() {
|
|
var e, i;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
var c1 = sphere.interpolate(e.c1, e.c2, 0);
|
|
expect(c1[0]).to.roughlyEqual(e.c1[0], 1e-9);
|
|
expect(c1[1]).to.roughlyEqual(e.c1[1], 1e-9);
|
|
var midpoint = sphere.interpolate(e.c1, e.c2, 0.5);
|
|
expect(midpoint[0]).to.roughlyEqual(e.midpoint[0], 1e-9);
|
|
expect(midpoint[1]).to.roughlyEqual(e.midpoint[1], 1e-5);
|
|
var c2 = sphere.interpolate(e.c1, e.c2, 1);
|
|
expect(c2[0]).to.roughlyEqual(e.c2[0], 1e-9);
|
|
expect(c2[1]).to.roughlyEqual(e.c2[1], 1e-5);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('midpoint', function() {
|
|
|
|
it('results match Chris Veness\'s reference implementation', function() {
|
|
var e, i, midpoint;
|
|
for (i = 0; i < expected.length; ++i) {
|
|
e = expected[i];
|
|
midpoint = sphere.midpoint(e.c1, e.c2);
|
|
// Test modulo 360 to avoid unnecessary expensive modulo operations
|
|
// in our implementation.
|
|
expect(goog.math.modulo(midpoint[0], 360)).to.roughlyEqual(
|
|
goog.math.modulo(e.midpoint[0], 360), 1e-9);
|
|
expect(midpoint[1]).to.roughlyEqual(e.midpoint[1], 1e-9);
|
|
}
|
|
});
|
|
|
|
});
|
|
|
|
describe('Vincenty area', function() {
|
|
var geometry;
|
|
before(function(done) {
|
|
afterLoadText('spec/ol/format/wkt/illinois.wkt', function(wkt) {
|
|
try {
|
|
var format = new ol.format.WKT();
|
|
geometry = format.readGeometry(wkt);
|
|
} catch (e) {
|
|
done(e);
|
|
}
|
|
done();
|
|
});
|
|
});
|
|
|
|
it('results match the expected area of Ilinois', function() {
|
|
var coords = geometry.getPolygon(0).getLinearRing(0).getCoordinates();
|
|
expect(ol.sphere.WGS84.geodesicArea(coords)).to.equal(145978332359.37125);
|
|
});
|
|
});
|
|
|
|
});
|
|
|
|
|
|
goog.require('goog.math');
|
|
goog.require('ol.Sphere');
|
|
goog.require('ol.sphere.WGS84');
|
|
goog.require('ol.format.WKT');
|