[Proj] Little geod problem

Mikhail Tchernychev misha at mail.geometrics.com
Thu Jan 20 11:21:46 PST 2011


  Hello list,

I came across a small inconsistency in the geod program. What I use is:

geod +ellps=WGS84  -I -f "%.3lf" data_file

With data file:

-12.4989444 130.88649687 -12.4989444 130.8864967

What I got is:

-1.#IO  -1.#IO  -1.#IO

These two points are very close to each other however are not identical. 
I would
not be surprised if geod would return 0 0 0. If fact this is what it 
does if points
are simply identical.  I would suggest to round result to 0, even if it 
is not but say
smaller than some reasonable number (0.01 m?)

  I hope this would help to improve great proj software even more.

Thank you,

Mikhail



More information about the Proj mailing list