MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: GeoDistance and transformation from WGS84 to WGS84

  • To: mathgroup at
  • Subject: [mg101559] Re: GeoDistance and transformation from WGS84 to WGS84
  • From: ragfield <ragfield at>
  • Date: Fri, 10 Jul 2009 06:44:17 -0400 (EDT)
  • References: <h340nr$gjs$>

On Jul 9, 12:55 am, "Sjoerd C. de Vries" <sjoerd.c.devr... at>
> Anyone an idea why the following doesn't work? It calculates the
> distance between two points on the earth with given latlong
> coordinates. WGS84 is the datum used for GPS, if I'm not mistaken.
> In[86]:= GeoDistance[GeoPosition[{1, 0}, "WGS84"],  GeoPosition[{0,
> 0}, "WGS84"]]
> During evaluation of In[86]:= GeoDistance::datum: No direct
> transformation between datums "WGS84" and "WGS84" is available. >>
> Out[86]= GeoDistance[GeoPosition[{1, 0}, "WGS84"],
>  GeoPosition[{0, 0}, "WGS84"]]
> The "More Information" box of GeoPosition specifically states that
> WGS84 is one of the allowed datums. However, GeodesyData["Datum"] and
> GeodesyData["ReferenceEllipsoid"] make clear that WGS84 is not a
> datum, but a reference ellipsoid (WGS72 seems to be both).

I don't know *why* it doesn't work, but if you change WGS84 to
WGS84Original it does work:

In[255]:= GeoDistance[GeoPosition[{1, 0}, "WGS84Original"],
 GeoPosition[{0, 0}, "WGS84Original"]]

Out[255]= 110574.


  • Prev by Date: Re: Selecting Rows Based on Column Values
  • Next by Date: Re: Problem with fitting
  • Previous by thread: GeoDistance and transformation from WGS84 to WGS84
  • Next by thread: Re: RecurrenceTable can't finish