Changeset 2903

Timestamp:

May 18, 2006, 10:19:51 AM (19 years ago)

Author:

sexton

Message:

check interior_regions fall in bounding_polygon - tells you which is the offending polygon (not the point)

File:

• production/pt_hedland_2006/run_pt_hedland.py (modified) (1 diff)

production/pt_hedland_2006/run_pt_hedland.py

@@ -125 +125 @@
 
 print 'number of interior regions', len(interior_regions)
+
+from utilities.polygon import inside_polygon
+
+bounding_polygon = project.polyAll
+count = 0
+for i in range(len(interior_regions)):
+    region = interior_regions[i]
+    interior_polygon = region[0]
+    if len(inside_polygon(interior_polygon, bounding_polygon,
+                   closed = True, verbose = False)) <> len(interior_polygon):
+        print 'WARNING: interior polygon %d is outside bounding polygon' %(i)
+        count += 1
+
+if count == 0:
+    print 'interior regions OK'
+else:
+    print 'check out your interior polygons'
+    import sys; sys.exit()
+    
 
 from caching import cache