|
Last change
on this file since 5401 was
4297,
checked in by jakeman, 17 years ago
|
|
Commmit new boxingday_2007 folder
|
|
File size:
843 bytes
|
| Line | |
|---|
| 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.
