Ignore:
Timestamp:
Feb 16, 2010, 9:52:30 AM (14 years ago)
Author:
fountain
Message:

updates to mandurah storm surge scenario

Location:
anuga_work/production/mandurah_storm_surge_2009
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/mandurah_storm_surge_2009/build_elevation.py

    r7614 r7627  
    8282print 'project.elevation_clip_box', project.elevation_clip_box
    8383
    84 G1 = geospatial_data['topog'].clip(project.elevation_clip_box)
    85 G2 = geospatial_data['man10m_ss'].clip_outside(project.elevation_clip_box)
     84# G1 = geospatial_data['Man_25m'].clip(project.elevation_clip_box)
     85# G2 = geospatial_data['ph10m_ss'].clip_outside(project.elevation_clip_box)
    8686
    87 G = G1 + G2
     87G1 = geospatial_data['m_peel_aoi'].clip(project.elevation_clip_box)
     88G2 = geospatial_data['ph10m_ss'].clip_outside(project.elevation_clip_box)
     89G3 = geospatial_data['MA-46893-SNDS_AHD.csv']
     90G4 = geospatial_data['MS0205HY_AHD.csv']
     91G5 = geospatial_data['MS0404_AHD.csv']
     92G6 = geospatial_data['YU0403HY_AHD.csv']
     93G7 = geospatial_data['original_data_ss.csv']
     94
     95G = G1 + G2 + G3 + G4 + G5 + G6 + G7
    8896       
    8997print 'Export combined DEM file'
  • anuga_work/production/mandurah_storm_surge_2009/project.py

    r7614 r7627  
    7373# Used in build_elevation.py
    7474# Format for ascii grids, as produced in ArcGIS + a projection file
    75 ##ascii_grid_filenames = ['m_peel.asc',
    76 ##                        'm_harvey.asc',
    77 ##                        'm_peel_10m.asc'
    78 ##                        'm_harvey_10m.asc']
    79 ascii_grid_filenames = ['Man_25m',              # this is the latest 25m DEM from GEMS
    80                                                                 'man10m_ss']    # this is to fill in areas not covered by the
     75
     76# This set of ascii files is utilising the original input data, not the GEMS elevation model
     77ascii_grid_filenames = ['m_peel_aoi',           # Original 1m lidar for PEEL, clipped to AOI
     78                        'ph10m_ss']             # Mosiaced PEEL-HARVEY Lidar, resampled to 10m,
     79                                                # clipped to bounding polygon
     80 
     81# This set of ascii files utilises the GEMS elevation model - for comparison with GEMS inundation
     82# ascii_grid_filenames = ['Man_25m',            # this is the latest 25m DEM from GEMS
     83                                                                # 'ph10m_ss']   # this is to fill in areas not covered by the
    8184                                                                # GEMS grid and is derived from the 10m resampled LiDAR from DoT
    8285
    8386# Format for point is x,y,elevation (with header)
    84 point_filenames = []
    85 ##                      ['MA-46893-SNDS_AHD.csv',
    86 ##                   'MS0205HY_AHD.csv',
    87 ##                   'MS0404_AHD.csv',
    88 ##                   'YU0403HY_AHD.csv']
     87# Don't use these for comparison with GEMS inundation
     88point_filenames = ['MA-46893-SNDS_AHD.csv',     # These files contain inf-fill bathymetry for
     89                  'MS0205HY_AHD.csv',           # the Peel & Harvey estuaries, as well as around
     90                  'MS0404_AHD.csv',             # the canals
     91                  'YU0403HY_AHD.csv',
     92                  'original_data_ss.csv']       # This data is from what we had available for the 2007
     93                                                # tsunami simulations for SW WA, clipped to bounding polygon
    8994
    9095# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
     
    128133# Used in get_timeseries.py. 
    129134# Format easting,northing,name,elevation (with header)
    130 gauges_filename = 'gauges.csv'
     135gauges_filename = 'gauges_mandurah_edited.csv'
    131136
    132137# BUILDINGS EXPOSURE - for identifying inundated houses
     
    155160# Elevation clip box - used when there are multiple overlapping datasets and you want to use
    156161# one over another in an area of overlap
    157 elevation_clip_box_filename = 'elevation_clip_box.csv'
     162#elevation_clip_box_filename = 'gems_elevation_clip_box.csv'
     163elevation_clip_box_filename = 'area_of_interest.csv'
    158164
    159165################################################################################
Note: See TracChangeset for help on using the changeset viewer.