Last change
on this file since 5175 was
4297,
checked in by jakeman, 18 years ago
|
Commmit new boxingday_2007 folder
|
File size:
843 bytes
|
Rev | Line | |
---|
[4297] | 1 | from anuga.geospatial_data.geospatial_data import * |
---|
| 2 | from anuga.coordinate_transforms.redfearn import * |
---|
| 3 | from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts,xya2pts |
---|
| 4 | |
---|
| 5 | #xya2pts('thaicoas_9', verbose=True) |
---|
| 6 | print'create G1' |
---|
| 7 | G1 = Geospatial_data(file_name = 'thaicoas_9' + '.xya') |
---|
| 8 | G = G1 |
---|
| 9 | G.export_points_file('boxingday' + '.pts') |
---|
| 10 | |
---|
| 11 | |
---|
| 12 | # Define extent and print into csv file |
---|
| 13 | points = [(6.0,97.5),(10.0,97.5),(10.0,100.5),(6.0,100.5)] |
---|
| 14 | points, zone = convert_from_latlon_to_utm(points) |
---|
| 15 | #print 'points', points |
---|
| 16 | #print 'zone', zone |
---|
| 17 | |
---|
| 18 | file_name = 'extent' +'.csv' |
---|
| 19 | f = file(file_name, 'w') |
---|
| 20 | for i in range(len(points)): |
---|
| 21 | #print '%f, %f' %(points[i][0],points[i][1]) |
---|
| 22 | f.write('%f, %f\n' %(points[i][0],points[i][1])) |
---|
| 23 | f.close |
---|
| 24 | |
---|
| 25 | # Define Interior regions and print into corresponding csv files |
---|
| 26 | # Coast |
---|
| 27 | |
---|
| 28 | # Island 1 |
---|
| 29 | |
---|
| 30 | # Island 2 |
---|
Note: See
TracBrowser
for help on using the repository browser.