Changeset 7251
- Timestamp:
- Jun 23, 2009, 2:39:22 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/new_south_wales/gosford/project.py
r7225 r7251 67 67 # Used in build_elevation.py 68 68 # Format for ascii grids, as produced in ArcGIS + a projection file 69 ascii_grid_filenames = ['e151_s34_clip2', 70 'aus197_topo2raster_clip2', 71 'hydro_topo2raster_clip2'] 72 #'aus197_topo2raster_terrigal'] #GA data?? 69 ascii_grid_filenames = [] 73 70 74 75 76 #Get long list of Lidar filenames.77 ENV_INUNDATIONHOME = 'INUNDATIONHOME'78 home = join(os.getenv(ENV_INUNDATIONHOME), 'data')79 anuga_folder = join(home, state, scenario_folder, 'anuga')80 topographies_folder = join(anuga_folder, 'topographies/original')81 82 #This reads a file containing lidar filenames and appends to a list83 file_list = join(topographies_folder, 'file_list.csv')84 fid = csv.reader(open(file_list))85 lidar_filenames = []86 for row in fid:87 lidar_filenames.append(','.join(row))88 89 90 71 # Format for point is x,y,elevation (with header) 91 point_filenames = ['estuaries.txt', 92 ## 'XYHI341_HSDB_SD100035103', 93 ## 'XYAHD', 94 ## 'XYSCHOOL_11_99_HSDB_SD100021002_MGA', 95 'ENT5C1S08_03_AHD_prepared.txt', 96 'ENT5C1S09_03_AHD_prepared.txt', 97 'ENT5C1S10_03_AHD_prepared.txt', 98 ## 'XYAHD', 99 ## 'XYHI341_HSDB_SD100035103_MGA', 100 ## 'XYHI339_F_HSDB_SD100031669_MGA', 101 ## 'AUS197_MGA_AHD_1972.txt', 102 'hydro_clip_neg.txt' 103 ] 104 point_filenames.extend(lidar_filenames) 105 #print point_filenames 106 107 108 ### Add csv header list to all files in point_filenames 109 ##headerlist = ['x', 'y', 'elevation'] 110 ##for f in point_filenames: 111 ## add_csv_header(join(topographies_folder, f), headerlist) 112 113 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh 72 point_filenames = [] 73 74 75 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh 114 76 # Used in build_elevation.py 115 77 # Format for points easting,northing (no header) … … 230 192 # check various directories/files that must exist 231 193 anuga_folder = join(home, state, scenario_folder, 'anuga') 232 topographies_folder = join(anuga_folder, 'topographies /original')194 topographies_folder = join(anuga_folder, 'topographies') 233 195 polygons_folder = join(anuga_folder, 'polygons') 234 196 boundaries_folder = join(anuga_folder, 'boundaries')
Note: See TracChangeset
for help on using the changeset viewer.