Changeset 7251


Ignore:
Timestamp:
Jun 23, 2009, 2:39:22 PM (15 years ago)
Author:
jgriffin
Message:

script cleaned up for DVD

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/new_south_wales/gosford/project.py

    r7225 r7251  
    6767# Used in build_elevation.py
    6868# 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??
     69ascii_grid_filenames = []   
    7370                   
    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 list
    83 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 
    9071# 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
     72point_filenames = []
     73       
     74
     75# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    11476# Used in build_elevation.py
    11577# Format for points easting,northing (no header)
     
    230192# check various directories/files that must exist
    231193anuga_folder = join(home, state, scenario_folder, 'anuga')
    232 topographies_folder = join(anuga_folder, 'topographies/original')
     194topographies_folder = join(anuga_folder, 'topographies')
    233195polygons_folder = join(anuga_folder, 'polygons')
    234196boundaries_folder = join(anuga_folder, 'boundaries')
Note: See TracChangeset for help on using the changeset viewer.