Changeset 2376 for inundation/pmesh
- Timestamp:
- Feb 10, 2006, 4:37:46 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/pmesh/create_mesh.py
r2351 r2376 13 13 14 14 if refzone is None: 15 zone = refzone15 refzone = zone 16 16 else: 17 17 assert zone == refzone … … 68 68 resolution, 69 69 filename = None, 70 interior_regions = None): 70 interior_regions = None, 71 UTM = False, 72 refzone = None): 71 73 """Create mesh from bounding polygon, tags for all segments and resolution. 72 74 … … 99 101 from coordinate_transforms.redfearn import redfearn 100 102 from utilities.polygon import populate_polygon 101 102 m = Mesh()103 from utilities.numerical_tools import ensure_numeric 104 from coordinate_transforms.geo_reference import Geo_reference 103 105 104 106 107 bounding_polygon = ensure_numeric(bounding_polygon) 108 #print 'refzone', refzone 109 105 110 #Convert to UTM 106 bounding_polygon, refzone = convert_points_from_latlon_to_utm(bounding_polygon) 111 if not UTM: 112 bounding_polygon, refzone = convert_points_from_latlon_to_utm(bounding_polygon) 113 else: 114 assert refzone is not None, 'Refzone must be specified, got %s' %refzone 115 116 117 #FIXME: Must give georeference here. 118 #Why can't Mesh work out the ll corners? 119 120 xllcorner = min(bounding_polygon[:,0]) 121 yllcorner = min(bounding_polygon[:,1]) 122 123 geo = Geo_reference(xllcorner = xllcorner, 124 yllcorner = yllcorner, 125 zone = refzone) 126 127 print 'geo reference derived from bounding polygon', geo 128 m = Mesh(geo_reference=geo) 129 107 130 108 131 … … 116 139 excluded_polygons = [] 117 140 for P, res in interior_regions: 118 polygon, _ = convert_points_from_latlon_to_utm(P, refzone) 141 if not UTM: 142 polygon, _ = convert_points_from_latlon_to_utm(P, refzone) 143 else: 144 polygon = P 119 145 excluded_polygons.append( polygon ) 120 146 else: … … 133 159 if interior_regions is not None: 134 160 for P, res in interior_regions: 135 polygon, _ = convert_points_from_latlon_to_utm(P, refzone) 161 if not UTM: 162 polygon, _ = convert_points_from_latlon_to_utm(P, refzone) 163 else: 164 polygon = P 165 136 166 region_dict = create_region(polygon, None, refzone) 137 167
Note: See TracChangeset
for help on using the changeset viewer.