- Timestamp:
- Jun 9, 2010, 12:28:24 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/source/anuga/file_conversion/test_sww2dem.py
r7776 r7814 132 132 133 133 #Export to ascii/prj files 134 sww2dem(self.domain.get_name(), 134 sww2dem(self.domain.get_name()+'.sww', 135 self.domain.get_name()+'_elevation.asc', 135 136 quantity = 'elevation', 136 137 cellsize = cellsize, 137 138 number_of_decimal_places = 9, 138 verbose = self.verbose, 139 format = 'asc') 139 verbose = self.verbose) 140 140 141 141 #Check prj (meta data) … … 220 220 os.remove(ascfile) 221 221 222 ascfile = self.domain.get_name() + '_depth.asc' 223 prjfile = self.domain.get_name() + '_depth.prj' 224 222 225 #Export to ascii/prj files 223 sww2dem(self.domain.get_name(), 226 sww2dem(self.domain.get_name()+'.sww', 227 ascfile, 224 228 quantity = 'depth', 225 229 cellsize = cellsize, 226 230 number_of_decimal_places = 9, 227 verbose = self.verbose, 228 format = 'asc') 231 verbose = self.verbose) 229 232 230 233 #Check asc file 231 ascfile = self.domain.get_name() + '_depth.asc'232 prjfile = self.domain.get_name() + '_depth.prj'233 234 ascid = open(ascfile) 234 235 lines = ascid.readlines() … … 353 354 354 355 #Export to ascii/prj files 355 sww2dem(domain.get_name(), 356 sww2dem(domain.get_name() + '.sww', 357 domain.get_name() + '_elevation.asc', 356 358 quantity = 'elevation', 357 359 cellsize = cellsize, 358 360 number_of_decimal_places = 9, 359 361 verbose = self.verbose, 360 format = 'asc',361 362 block_size=2) 362 363 … … 525 526 526 527 #Export to ascii/prj files 527 sww2dem(domain.get_name(), 528 sww2dem(domain.get_name()+'.sww', 529 domain.get_name() + '_elevation.asc', 528 530 quantity = 'elevation', 529 531 cellsize = cellsize, 530 532 number_of_decimal_places = 9, 531 533 verbose = self.verbose, 532 format = 'asc',533 534 block_size=2) 534 535 … … 720 721 721 722 # Export to ascii/prj files 722 sww2dem(domain.get_name(), 723 sww2dem(domain.get_name() + '.sww', 724 domain.get_name() + '_elevation.asc', 723 725 quantity = 'elevation', 724 726 cellsize = cellsize, … … 728 730 northing_min = 6189050, 729 731 northing_max = 6189100, 730 verbose = self.verbose, 731 format = 'asc') 732 verbose = self.verbose) 732 733 733 734 fid.close() … … 867 868 868 869 #Export to ascii/prj files 869 sww2dem(self.domain.get_name(), 870 sww2dem(self.domain.get_name() + '.sww', 871 self.domain.get_name() + '_stage.asc', 870 872 quantity = 'stage', 871 873 cellsize = cellsize, 872 874 number_of_decimal_places = 9, 873 reduction = min, 874 format = 'asc', 875 verbose=self.verbose) 875 reduction = min) 876 876 877 877 … … 974 974 975 975 #Export to ascii/prj files 976 sww2dem(self.domain.get_name(), 976 sww2dem(self.domain.get_name() + '.sww', 977 self.domain.get_name() + '_stage.asc', 977 978 quantity = 'stage', 978 979 cellsize = cellsize, 979 980 number_of_decimal_places = 9, 980 981 reduction = 1, 981 format = 'asc',982 982 verbose=self.verbose) 983 983 … … 1082 1082 1083 1083 #Export to ascii/prj files 1084 sww2dem(self.domain.get_name() ,1085 basename_out = 'datatest_depth',1084 sww2dem(self.domain.get_name()+'.sww', 1085 name_out = 'datatest_depth.asc', 1086 1086 quantity = 'stage - elevation', 1087 1087 cellsize = cellsize, 1088 1088 number_of_decimal_places = 9, 1089 1089 reduction = min, 1090 format = 'asc',1091 1090 verbose = self.verbose) 1092 1091 … … 1234 1233 1235 1234 #Export to ascii/prj files 1236 sww2dem(domain.get_name(), 1235 sww2dem(domain.get_name()+'.sww', 1236 domain.get_name()+'_elevation.asc', 1237 1237 quantity = 'elevation', 1238 1238 cellsize = cellsize, 1239 1239 number_of_decimal_places = 9, 1240 verbose = self.verbose, 1241 format = 'asc') 1240 verbose = self.verbose) 1242 1241 1243 1242 … … 1343 1342 1344 1343 #Export to ers files 1345 sww2dem(self.domain.get_name(), 1344 outname = self.domain.get_name() + '_elevation.ers' 1345 sww2dem(self.domain.get_name() + '.sww', 1346 outname, 1346 1347 quantity = 'elevation', 1347 1348 cellsize = cellsize, 1348 1349 number_of_decimal_places = 9, 1349 1350 NODATA_value = NODATA_value, 1350 verbose = self.verbose, 1351 format = 'ers') 1351 verbose = self.verbose) 1352 1352 1353 1353 #Check header data 1354 1354 from ermapper_grids import read_ermapper_header, read_ermapper_data 1355 1355 1356 header = read_ermapper_header( self.domain.get_name() + '_elevation.ers')1357 #print header 1356 header = read_ermapper_header(outname) 1357 1358 1358 assert header['projection'].lower() == '"utm-56"' 1359 1359 assert header['datum'].lower() == '"wgs84"'
Note: See TracChangeset
for help on using the changeset viewer.