Changeset 2442

Timestamp:

Feb 22, 2006, 6:36:40 PM (19 years ago)

Author:

ole

Message:

Test that point duplication is OK

File:

• inundation/utilities/test_polygon.py (modified) (1 diff)

inundation/utilities/test_polygon.py

@@ -241 +241 @@
         assert not inside_polygon( (0.5, 0.5), polygon )
 
+
+
+    def test_duplicate_points_being_ok(self):
+
+
+        #Simplest case: Polygon is the unit square
+        polygon = [[0,0], [1,0], [1,0], [1,0], [1,1], [0,1], [0,0]]
+
+        assert inside_polygon( (0.5, 0.5), polygon )
+        assert not inside_polygon( (0.5, 1.5), polygon )
+        assert not inside_polygon( (0.5, -0.5), polygon )
+        assert not inside_polygon( (-0.5, 0.5), polygon )
+        assert not inside_polygon( (1.5, 0.5), polygon )
+
+        #Try point on borders
+        assert inside_polygon( (1., 0.5), polygon, closed=True)
+        assert inside_polygon( (0.5, 1), polygon, closed=True)
+        assert inside_polygon( (0., 0.5), polygon, closed=True)
+        assert inside_polygon( (0.5, 0.), polygon, closed=True)
+
+        assert not inside_polygon( (0.5, 1), polygon, closed=False)
+        assert not inside_polygon( (0., 0.5), polygon, closed=False)
+        assert not inside_polygon( (0.5, 0.), polygon, closed=False)
+        assert not inside_polygon( (1., 0.5), polygon, closed=False)
+
+        #From real example
+        polygon = [[20,20], [40,20], [40,40], [20,40]]
+        points = [ [40, 50] ]
+        res = inside_polygon(points, polygon)
+        assert len(res) == 0
+
+        
+
     def test_inside_polygon_vector_version(self):
         #Now try the vector formulation returning indices