Changeset 3183

Timestamp:

Jun 20, 2006, 9:16:41 AM (19 years ago)

Author:

duncan

Message:

bug fixing

File:

• inundation/pmesh/mesh_interface.py (modified) (3 diffs)

inundation/pmesh/mesh_interface.py

@@ -73 +73 @@
     # geo_ref
     if interior_regions is not None:        
-     
+        # Test that all the interior polygons are inside the bounding_poly
         for interior_polygon, res in interior_regions:
             indices = inside_polygon(interior_polygon, bounding_polygon,
@@ -138 +138 @@
     if interior_regions is not None:    
         for polygon, res in interior_regions:
-            #polygon = convert_points_from_latlon_to_utm(P, refzone)
             m.add_region_from_polygon(polygon,
                                       geo_reference=poly_geo_reference)
@@ -147 +146 @@
             else:
                 polygon_absolute = polygon
-            [inner_point] = populate_polygon(polygon_absolute, 1)
+            # This isn't working it seems    
+            #[inner_point] = populate_polygon(polygon_absolute, 1)
+            inner_point = point_in_polygon(polygon_absolute)
             region = m.add_region(inner_point[0], inner_point[1])
             region.setMaxArea(res)