| [4949] | 1 | """ |
|---|
| 2 | Convert Easting Northing to lon/lat for ultimately get kml file |
|---|
| 3 | """ |
|---|
| 4 | |
|---|
| 5 | def alter_file(filename,filenameout): |
|---|
| 6 | |
|---|
| 7 | from anuga.coordinate_transforms.lat_long_UTM_conversion import UTMtoLL |
|---|
| 8 | |
|---|
| 9 | fid = open(filename) |
|---|
| 10 | lines = fid.readlines() |
|---|
| 11 | fid.close() |
|---|
| 12 | |
|---|
| 13 | fid_out = open(filenameout, 'w') |
|---|
| 14 | s = '%s,%s,%s\n' %('lon','lat','stage') |
|---|
| 15 | fid_out.write(s) |
|---|
| 16 | for line in lines[1:]: |
|---|
| 17 | fields = line.split(',') |
|---|
| 18 | easting = float(fields[0]) |
|---|
| 19 | northing = float(fields[1]) |
|---|
| 20 | stage = float(fields[3]) |
|---|
| 21 | lat, lon = UTMtoLL(northing,easting, 50) |
|---|
| 22 | s = '%f,%f,%f\n' %(lon,lat,stage) |
|---|
| 23 | fid_out.write(s) |
|---|
| 24 | |
|---|
| 25 | fid_out.close() |
|---|
| 26 | |
|---|
| 27 | return |
|---|
| 28 | |
|---|
| 29 | import project |
|---|
| 30 | from os import sep |
|---|
| 31 | |
|---|
| 32 | # Dampier |
|---|
| 33 | #timedir = '20070618_063824_run_final_2.4_onslow_nbartzis' |
|---|
| 34 | timedir = '20070620_044413_run_final_2.4_exmouth_nbartzis' |
|---|
| 35 | #timedir = '20070419_065050_run_final_2.4_dampier_nbartzi |
|---|
| 36 | |
|---|
| 37 | # Karratha |
|---|
| 38 | |
|---|
| 39 | filename = project.outputdir + timedir + sep + 'stage0_interp.csv' |
|---|
| 40 | #filename = project.outputdir + timedir + sep + 'stage5_interp.csv' |
|---|
| 41 | #filename = project.outputdir+ timedir + sep + 'stage20_interp.csv' |
|---|
| 42 | #filename = project.outputdir+ timedir + sep + 'stage80_interp.csv' |
|---|
| 43 | filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv' |
|---|
| 44 | #filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv' |
|---|
| 45 | #filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv' |
|---|
| 46 | #filename_out = project.outputdir + timedir + sep + 'stage80_convert.csv' |
|---|
| 47 | alter_file(filename,filename_out) |
|---|
