from anuga.geospatial_data.geospatial_data import * from anuga.coordinate_transforms.redfearn import * from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts,xya2pts #xya2pts('thaicoas_9', verbose=True) print'create G1' G1 = Geospatial_data(file_name = 'thaicoas_9' + '.xya') G = G1 G.export_points_file('boxingday' + '.pts') # Define extent and print into csv file points = [(6.0,97.5),(10.0,97.5),(10.0,100.5),(6.0,100.5)] points, zone = convert_from_latlon_to_utm(points) #print 'points', points #print 'zone', zone file_name = 'extent' +'.csv' f = file(file_name, 'w') for i in range(len(points)): #print '%f, %f' %(points[i][0],points[i][1]) f.write('%f, %f\n' %(points[i][0],points[i][1])) f.close # Define Interior regions and print into corresponding csv files # Coast # Island 1 # Island 2