[Proj] invgeod results for small distances
Charles Karney
charles.karney at sri.com
Thu Feb 23 07:57:43 EST 2012
I recommend using the Geod program provided with GeographicLib
http://geographiclib.sf.net
Thus:
Geod -i -b -d -p 6 <<EOF
0 0 0 .000001
0 0 0 .0000001
-30 0 29.9 179.9
EOF
=>
090d00'00.0000000" 270d00'00.0000000" 0.111319
090d00'00.0000000" 270d00'00.0000000" 0.011132
170d59'40.4518759" 188d59'46.8569451" 19992090.302327
The last example is another case where invgeod gets the wrong answer.
(Disclaimer: I'm the author of GeographicLib.)
--Charles
On 02/22/12 23:59, Hermann Peifer wrote:
> # This looks reasonable for a difference in the 6th decimal
> [peifer at whitefish:~]> echo 0 0 0 .000001 | invgeod +ellps=WGS84
> 90d -90d 0.134
>
>
> # This should then be 0.013 for the 7th decimal, but it isn't:-(
> [peifer at whitefish:~]> echo 0 0 0 .0000001 | invgeod +ellps=WGS84
> -2147483648d-2147483648'-nan" -2147483648d-2147483648'-nan" -nan
>
>
> # This isn't true in real life, but at least better than Not-A-Number
> [peifer at whitefish:~]> echo 0 0 0 .0000001 | invgeod +ellps=sphere
> 90d -90d 0.000
More information about the Proj
mailing list