Changeset 7627
- Timestamp:
- Feb 16, 2010, 9:52:30 AM (15 years ago)
- 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 82 82 print 'project.elevation_clip_box', project.elevation_clip_box 83 83 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) 86 86 87 G = G1 + G2 87 G1 = geospatial_data['m_peel_aoi'].clip(project.elevation_clip_box) 88 G2 = geospatial_data['ph10m_ss'].clip_outside(project.elevation_clip_box) 89 G3 = geospatial_data['MA-46893-SNDS_AHD.csv'] 90 G4 = geospatial_data['MS0205HY_AHD.csv'] 91 G5 = geospatial_data['MS0404_AHD.csv'] 92 G6 = geospatial_data['YU0403HY_AHD.csv'] 93 G7 = geospatial_data['original_data_ss.csv'] 94 95 G = G1 + G2 + G3 + G4 + G5 + G6 + G7 88 96 89 97 print 'Export combined DEM file' -
anuga_work/production/mandurah_storm_surge_2009/project.py
r7614 r7627 73 73 # Used in build_elevation.py 74 74 # 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 77 ascii_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 81 84 # GEMS grid and is derived from the 10m resampled LiDAR from DoT 82 85 83 86 # 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 88 point_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 89 94 90 95 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh … … 128 133 # Used in get_timeseries.py. 129 134 # Format easting,northing,name,elevation (with header) 130 gauges_filename = 'gauges .csv'135 gauges_filename = 'gauges_mandurah_edited.csv' 131 136 132 137 # BUILDINGS EXPOSURE - for identifying inundated houses … … 155 160 # Elevation clip box - used when there are multiple overlapping datasets and you want to use 156 161 # 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' 163 elevation_clip_box_filename = 'area_of_interest.csv' 158 164 159 165 ################################################################################
Note: See TracChangeset
for help on using the changeset viewer.