Changeset 3054 for inundation/utilities/test_polygon.py

Timestamp:

Jun 2, 2006, 3:02:36 PM (19 years ago)

Author:

duncan

Message:

removing geo-ref from inside_polygon

File:

• inundation/utilities/test_polygon.py (modified) (3 diffs)

inundation/utilities/test_polygon.py

@@ -5 +5 @@
 from Numeric import zeros, array, allclose
 from math import sqrt, pi
+from utilities.numerical_tools import ensure_numeric
 
 from polygon import *
 from coordinate_transforms.geo_reference import Geo_reference
+from geospatial_data.geospatial_data import Geospatial_data
 
 def test_function(x, y):
@@ -631 +633 @@
 
         
-    def test_inside_polygon_geo_ref(self):
+    def test_inside_polygon_geospatial(self):
+
+
+        polygon_absolute = [[0,0], [1,0], [1,1], [0,1]]
+        poly_geo_ref = Geo_reference(57,100,100)
+        
+
 
 
@@ -638 +646 @@
         poly_geo_ref = Geo_reference(57,100,100)
         polygon = poly_geo_ref.change_points_geo_ref(polygon_absolute)
-
+        poly_spatial = Geospatial_data(polygon,
+                                       geo_reference=poly_geo_ref)
+        
         points_absolute = (0.5, 0.5)
         points_geo_ref = Geo_reference(57,78,-56)
         points = points_geo_ref.change_points_geo_ref(points_absolute)
+        points_spatial = Geospatial_data(points,
+                                         geo_reference=points_geo_ref) 
+        
+        assert is_inside_polygon(points_absolute, polygon_absolute)
+        assert is_inside_polygon(ensure_numeric(points_absolute),
+                                 ensure_numeric(polygon_absolute))
+        assert is_inside_polygon(points_absolute, poly_spatial)
+        assert is_inside_polygon(points_spatial, poly_spatial)
+        assert is_inside_polygon(points_spatial, polygon_absolute)
 
         assert is_inside_polygon(points_absolute, polygon_absolute)
         
-        points_inside = inside_polygon(points, polygon,
-                               polygon_geo_ref=poly_geo_ref,
-                               points_geo_ref=points_geo_ref)
-        #print "points_inside",points_inside 
-        #assert 1 == len(inside_polygon( points, polygon,
-        #                       polygon_geo_ref=poly_geo_ref,
-        #                       points_geo_ref=points_geo_ref))
-        assert not is_inside_polygon((0.5, 1.5), polygon,
-                                     polygon_geo_ref=poly_geo_ref )
-        assert not is_inside_polygon((0.5, -0.5), polygon,
-                                     polygon_geo_ref=poly_geo_ref )
-        assert not is_inside_polygon((-0.5, 0.5), polygon,
-                                     polygon_geo_ref=poly_geo_ref )
-        assert not is_inside_polygon((1.5, 0.5), polygon,
-                                     polygon_geo_ref=poly_geo_ref )
+        
 #-------------------------------------------------------------
 if __name__ == "__main__":