[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