[Proj] Natural Earth projection

OvV_HN ovv at hetnet.nl
Fri Mar 9 06:23:20 EST 2012

[Proj] Natural Earth projection
Bernhard Jenny jennyb at geo.oregonstate.edu
Thu Mar 8 18:46:20 EST 2012

Hi Berhard,

Regarding a possible error in the inverse projection function Natural Earth:

Thanks for looking into this problem!
Luckily it was my fault, excuses for that.

The problem was not quite standard (JavaScript) code in the original, which 
I cleaned up. In this process a small typographic error was made, with large 
Your code will undoubtedly pass many JavaScript interpreters, but not all.
Especially the following two constructs are likely to give somewhere a 
    for(; ; )
    yc -= tol = f / fder;

My corrected code now goes all the way up to the poles, for instance:

R = 6371008.7714;
lat = 90.0d, lon = 120.0d
gives x = 6394608.652918302, y = 9062062.394736718
back to lat, lon: 90.00000000000014, 120.00000000000055 d
in 7 Newton iteration steps.

If someone could check the provided PROJ C-code, the Natural Earth 
projection code could be added to the PROJ package, if so desired!

Oscar van Vlijmen

More information about the Proj mailing list