New in Wolfram Mathematica 7: Integrated Geodesy & GIS  previous | next 
Translate Between Coordinate Systems
Easily convert {lat, long} coordinates to the XYZ coordinate system.
In[1]:=

Click for copyable input
globe = ParametricPlot3D[

   GeoPositionXYZ[GeoPosition[{lat, long}]][[1]], {lat, -90, 

    90}, {long, -180, 180}, MaxRecursion -> 4, Axes -> None];
In[2]:=

Click for copyable input
Show[Graphics3D[{Yellow, Thick, 

   Line /@ Map[GeoPositionXYZ[GeoPosition[#]][[1]] &, 

     Join[CountryData[#, "Coordinates"] & /@ {"UnitedStates", 

        "Canada", "Mexico"}], {-2}]}], globe, Boxed -> False]
Out[2]=