[Proj] What is the expected Proj4 result when an input is out of range?

Mikael Rittri Mikael.Rittri at carmenta.com
Fri Feb 23 05:27:52 PST 2018


Hello, list.

I am confused about the result I get from the Proj4 inverse projection of Lambert Azimuthal Equal Area,
when the original projected input point is outside the valid circle.

I thought that in cases like these, the result point should have coordinates that are HUGE_VAL ,  that is, 1e+300.
There seems to be code in pj_inv.c that checks error conditions and sets coordinates to HUGE_VAL, but these
checks aren't triggered in my example.
Instead, my debugger says that I get a point with coordinates  -nan(ind).
I suppose that is an Indefinite NaN.  I am not too familiar with the various flavours of NaN, but it's certainly different from 1e+300.

Here is a concrete example that I have run in the console application, using Proj 4.9.3 and EPSG:3035 which is ETRS89 / LAEA Europe:

C:\> proj -I +init=epsg:3035
4321000 3210000
10dE    52dN
16350000 -1090000
-2147483648d-2147483648'-1.#IO"E        -2147483648d-2147483648'-1.#IO"N

I think I remember that invalid coordinates are usually written as asterisks by the console application.

Avoiding the default DMS formatting of the output, I get the following instead:

C:\>proj -I +init=epsg:3035 -f "%.5f"
4321000 3210000
10.00000        52.00000
16350000 -1090000
-1.#IND0        -1.#IND0

So my question is, is this a kind of bug in the implementation of inverse LAEA, or is it to be expected that invalid coordinates
can be represented by either HUGE_VAL or NaN depending on the projection?

Best regards,

Mikael Rittri
Carmenta
Sweden
http://www.carmenta.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20180223/62aa5364/attachment.html>


More information about the Proj mailing list