Changeset 1144
- Timestamp:
- Mar 24, 2005, 4:22:35 PM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/pyvolution/data_manager.py
r1143 r1144 1102 1102 volumes = fid.variables['volumes'][:] 1103 1103 1104 vertex_points = concatenate ((x[:, NewAxis] ,y[:, NewAxis]), axis = 1)1105 assert len(vertex_points.shape) == 21106 1107 1104 ymin = min(y); ymax = max(y) 1108 1105 xmin = min(x); xmax = max(x) … … 1151 1148 prjid.close() 1152 1149 1153 #Create grid and update xll/yll corner 1150 #Create grid and update xll/yll corner and x,y 1154 1151 if verbose: print 'Creating grid' 1155 1152 ncols = int((xmax-xmin)/cellsize)+1 1156 1153 nrows = int((ymax-ymin)/cellsize)+1 1157 1154 1158 xllcorner = xmin+xllcorner 1159 yllcorner = ymin+yllcorner 1160 1155 newxllcorner = xmin+xllcorner 1156 newyllcorner = ymin+yllcorner 1157 1158 x = x+xllcorner-newxllcorner 1159 y = y+yllcorner-newyllcorner 1160 1161 vertex_points = concatenate ((x[:, NewAxis] ,y[:, NewAxis]), axis = 1) 1162 assert len(vertex_points.shape) == 2 1163 1161 1164 1162 1165 from Numeric import zeros, Float … … 1198 1201 ascid.write('ncols %d\n' %ncols) 1199 1202 ascid.write('nrows %d\n' %nrows) 1200 ascid.write('xllcorner %d\n' % xllcorner)1201 ascid.write('yllcorner %d\n' % yllcorner)1203 ascid.write('xllcorner %d\n' %newxllcorner) 1204 ascid.write('yllcorner %d\n' %newyllcorner) 1202 1205 ascid.write('cellsize %f\n' %cellsize) 1203 1206 ascid.write('NODATA_value %d\n' %NODATA_value) … … 1209 1212 #print grid_points 1210 1213 #print grid_values 1214 1215 inside_indices = inside_polygon(grid_points, P) 1211 1216 1212 1217 for i in range(nrows): 1213 if verbose and i%(( i+10)/10)==0:1218 if verbose and i%((nrows+1)/10)==0: 1214 1219 print 'Doing row %d of %d' %(i, nrows) 1215 1220 … … 1217 1222 index = (nrows-i-1)*ncols+j 1218 1223 1219 if inside_polygon(grid_points[index], P): 1224 #if inside_polygon(grid_points[index], P): 1225 if index in inside_indices: 1220 1226 ascid.write('%f ' %grid_values[index]) 1221 1227 else:
Note: See TracChangeset
for help on using the changeset viewer.