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