From 111ce96fcadd9bd700c7000d6f22fea60211d8c7 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Sun, 15 Apr 2018 18:35:38 +0200 Subject: Update +proj=aeqd +guam test case with data from EPSG Note 7-2 --- test/gie/builtins.gie | 34 ++++++++++++---------------------- 1 file changed, 12 insertions(+), 22 deletions(-) diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index fb8a53c3..63e56967 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -124,30 +124,20 @@ expect 3_860_398.3783 5_430_089.0490 roundtrip 100 ------------------------------------------------------------------------------- -Test equatorial aspect of the ellipsoidal azimuthal equidistant. Test data from -Snyder pp. 196-197, table 30. +# Test the azimuthal equidistant modified for Guam. Test data from the EPSG +Guidance Note 7 part 2, September 2016, p. 85 ------------------------------------------------------------------------------- -operation +proj=aeqd +ellps=GRS80 +guam +operation +proj=aeqd +guam +ellps=clrk66 +x_0=50000.00 +y_0=50000.00 + +lon_0=144.74875069444445 +lat_0=13.47246633333333 ------------------------------------------------------------------------------- -tolerance 1 m -accept 0 0 -expect 0.0000 0.0000 -roundtrip 100 -accept 90 90 -except 0.0000 10_001_965.7292 -roundtrip 100 -accept 0 90 -expect 0.0000 10_001_965.7292 -roundtrip 100 -accept 90 90 -expect 0.0000 10_001_965.7292 -roundtrip 100 -accept 45 45 -expect 3548107.5793 5970183.542 -#roundtrip 100 -accept -45 -45 -expect -3548107.5793 -5970183.542 -#roundtrip 100 +tolerance 1 cm +accept 144.635331291666660 13.33903846111111 +expect 37712.48 35242.00 +roundtrip 100 + +direction inverse +accept 37712.48 35242.00 +expect 144.635331291666660 13.33903846111111 ------------------------------------------------------------------------------- Test northern polar aspect of the ellipsoidal azimuthal equidistant. Test data -- cgit v1.2.3 From b0fa1a508d4827032f81ae03d29a54e3091ae47b Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Sun, 15 Apr 2018 18:41:48 +0200 Subject: Update aeqd docs with note of area of use for +guam option --- docs/source/operations/projections/aeqd.rst | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/docs/source/operations/projections/aeqd.rst b/docs/source/operations/projections/aeqd.rst index 32165e49..d4baa8b6 100644 --- a/docs/source/operations/projections/aeqd.rst +++ b/docs/source/operations/projections/aeqd.rst @@ -12,7 +12,10 @@ Azimuthal Equidistant +---------------------+----------------------------------------------------------+ | **Options** | +---------------------+----------------------------------------------------------+ -| `+guam` | Use Guam elliptical formulas. Defaults to false. | +| `+guam` | Use Guam elliptical formulas. Only accurate near the | +| | Island of Guam (:math:`\lambda\approx 144.5^{\circ}`, | +| | :math:`\phi\approx 13.5^{\circ}`) | +| | Defaults to false. | +---------------------+----------------------------------------------------------+ -- cgit v1.2.3