Changeset 3056 for inundation/geospatial_data/geospatial_data.py

Timestamp:

Jun 2, 2006, 3:10:48 PM (19 years ago)

Author:

duncan

Message:

Fixing mesh_interface tests. geospatial_data.py now has ensure_geospatial.

File:

• inundation/geospatial_data/geospatial_data.py (modified) (1 diff)

inundation/geospatial_data/geospatial_data.py

@@ -704 +704 @@
     return points
      
-        
+
+def ensure_geospatial(points, geo_reference = None):
+    """
+    This function inputs several formats and
+    outputs one format. - a geospatial_data instance.
+
+    Inputed formats are;
+    points: List or numeric array of coordinate pairs [xi, eta] of
+              points or geospatial object
+
+    mesh_origin: A geo_reference object or 3-tuples consisting of
+                 UTM zone, easting and northing.
+                 If specified vertex coordinates are assumed to be
+                 relative to their respective origins.
+    """
+    if isinstance(points,Geospatial_data):
+        #points = points.get_data_points( \
+        #        absolute = True)
+        msg = "Use a Geospatial_data object or a mesh origin. Not both."
+        assert geo_reference == None, msg
+            
+    else:
+        points = ensure_numeric(points, Float)
+    if geo_reference is None:
+        geo = None #Geo_reference()
+    else:
+        if isinstance(geo_reference, Geo_reference):
+            geo = geo_reference
+        else:
+            geo = Geo_reference(geo_reference[0],
+                                geo_reference[1],
+                                geo_reference[2])
+        points = Geospatial_data(data_points=points, geo_reference=geo)
+        #points = geo.get_absolute(points)
+    return points
+