[Proj] Any access to geodetic[sic] functions ??
Clifford J Mugnier
cjmce at lsu.edu
Wed Nov 12 12:33:47 EST 2008
In regard to terminology,
The shortest distance on the surface of a solid is generally termed a geodesic, be it an ellipsoid of revolution, aposphere, etc. On a sphere, the geodesic is termed a Great Circle.
HOWEVER, when computing the distance between two points using a projected coordinate system, that is a conformal projection such as Transverse Mercator, Oblique Mercator, Normal Mercator, Stereographic, or Lambert Conformal Conic - that then is a GRID distance which can be converted to an equivalent GEODETIC distance using the function for "Scale Factor at a Point." The conversion is then termed "Grid Distance to Geodetic Distance," even though it will not be as exactly correct as a true ellipsoidal geodesic. Closer to the truth with a TM than with a Lambert or other conformal projection, but still not exactly "on."
So, it can be termed "geodetic distance" or a "geodesic distance," depending on just how you got there ...
LOUISIANA STATE UNIVERSITY
From: proj-bounces at lists.maptools.org on behalf of Gerald I. Evenden
Sent: Wed 12-Nov-08 10:44
To: PROJ.4 and general Projections Discussions
Subject: Re: [Proj] Any access to geodetic[sic] functions ??
I've tried to stay out of this thread, but ... .
Unfortunately it keeps coming up.
BTW: when referring to distance calculations it is "geodesic" not geodetic.
On Wednesday 12 November 2008 10:50:37 am Frank Warmerdam wrote:
> Benoît Andrieu wrote:
> > Thank you so much for your answer, Frank.
> > So, this leads to some questions (that are somehow related...) :
> > - how are computed distance calculation between two points by softwares
> > like Postgis for example ?
> I'm not sure.
> > - is the best way to compute distance calculation to use UTM systems ? I
> > have to compute distances that could be very long (hundreds of
> > kilometers), strong accuracy is not required but I don't want to have
> > more than 1% accuracy error.
> Other have addressed this, but generally this is not a very good approach
> except for fairly local areas.
Very true and even if short distances are computed this way then one needs to
refer to the local scale factor and appropriately correct the pythagorean
computation. If you are using lproj to determine the UTM (or tmerc)
coordinates, turn the -V switch on to get the scale factor. Taking the mean
value for the two points should be adequate.
If you want to be lazy *and* do it the right way then for a very good
procedure to calculate inter point distances see:
and algorithm details published in:
Note: the first reference will point you to an online procedure to compute
We need to put the above in a standard location so that we can paste it into a
response any query to this group related to the geodesic.
> > - in what way are those functions distinct from the proj4 ? Am I wrong
> > to use them ? Are there others functions in proj4 that would allow me to
> > do the same things that I can do with those functions (distance
> > calculation between two points and point positionning given initial
> > point, distance and azimuth) ?
> The functions are not a core part of PROJ.4 as I understand it. They
> aren't normally in the library - just linked with the geod command line
At one time a program called 'geod' was distributed with my versions of proj4
but it was based upon a poorer algorithm has has been dropped. I believe it
my still be floating around in some distributions (and is probably good
enough for most applications).
I made a C version of the Vincente procedures but have not distributed it.
> I'm not opposed to treating them as a core capability but I do think some
> consideration should be given to the API before this would occur. Ideally
> this would be something someone else than me might take on (possibly with
> some advice from me).
> Best regards,
The whole religious complexion of the modern world is due
to the absence from Jerusalem of a lunatic asylum.
-- Havelock Ellis (1859-1939) British psychologist
Proj mailing list
Proj at lists.maptools.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Proj