Move height_func to RoamPoint and add lat-lon cache

[grits] / src / gis-world.c

diff --git a/src/gis-world.c b/src/gis-world.c
index 6e9beaadbe88987a8cbd7f1cffe08b4fc0e09b0c..c4d74b836e302cfdf32cd42773b034eb1fbe6240 100644 (file)
--- a/src/gis-world.c
+++ b/src/gis-world.c
@@ -151,6 +151,13 @@ void xyz2lle(gdouble x, gdouble y, gdouble z,
        *elev = rad2elev(rad);
 }
 
+void xyz2ll(gdouble x, gdouble y, gdouble z,
+               gdouble *lat, gdouble *lon)
+{
+       gdouble rad = sqrt(x*x + y*y + z*z);
+       *lat = incl2lat(acos(y / rad));
+       *lon = azim2lon(atan2(x,z));
+}
 
 gdouble ll2m(gdouble lon_dist, gdouble lat)
 {