Changeset 4058 for anuga_work/production
- Timestamp:
- Dec 5, 2006, 5:30:41 PM (18 years ago)
- Location:
- anuga_work/production
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/newcastle_2006/project_slide.py
r4046 r4058 27 27 # onshore data provided by NSW LPI, Krishna merged with existing 28 28 # PMD data to make 100m DEM 29 on_offshore_name = '' 29 nsw100 = 'dem_utm' 30 30 31 # AHO 32 offshore_name1 = 'XY100011610' 31 # survey data 32 offshore_name2 = 'surveyAreaB' 33 offshore_name3 = 'surveyAreaC' 34 35 # AHO data 36 offshore_name4 = '1000003611export' 37 offshore_name5 = '1000003613export' 38 offshore_name6 = '1000003614export' 39 offshore_name7 = '1000003627export' 40 offshore_name8 = '1000003628export' 41 offshore_name9 = 'AHDexport' 33 42 34 43 #swollen/ all data output … … 54 63 # Necessary if using point datasets, rather than grid 55 64 on_offshore_dem_name = datadir + on_offshore_name 56 offshore_dem_name1 = datadir + offshore_name1 65 offshore_dem_name2 = datadir + offshore_name2 66 offshore_dem_name3 = datadir + offshore_name3 67 offshore_dem_name4 = datadir + offshore_name4 68 offshore_dem_name5 = datadir + offshore_name5 69 offshore_dem_name6 = datadir + offshore_name6 70 offshore_dem_name7 = datadir + offshore_name7 71 offshore_dem_name8 = datadir + offshore_name8 72 offshore_dem_name9 = datadir + offshore_name9 57 73 combined_dem_name = datadir + 'newcastle_combined_elevation' 58 74 … … 82 98 # Clipping regions for export to asc and regions for clipping data 83 99 ################################################################### 100 101 # clipping for nsw grid for conversion 102 eastingmin_nsw = 546000 103 eastingmax_nsw = 350390 104 northingmin_nsw = 6253970 105 northingmax_nsw = 6424600 84 106 85 107 # exporting asc grid -
anuga_work/production/newcastle_2006/run_newcastle_slide.py
r4036 r4058 52 52 53 53 # filenames 54 on_offshore_dem_name = project.on_offshore_dem_name55 meshname = project .meshname+'.msh'54 nsw_dem_name = project_slide.nsw_dem_name 55 meshname = project_slide.meshname+'.msh' 56 56 57 57 # creates DEM from asc data 58 convert_dem_from_ascii2netcdf( on_offshore_dem_name, use_cache=True, verbose=True)58 convert_dem_from_ascii2netcdf(nsw_dem_name, use_cache=True, verbose=True) 59 59 60 60 #creates pts file for onshore DEM 61 dem2pts(on_offshore_dem_name, use_cache=True, verbose=True) 61 dem2pts(nsw_dem_name, 62 easting_min=project_slide.eastingmin_nsw, 63 easting_max=project_slide.eastingmax_nsw, 64 northing_min=project_slide.northingmin_nsw, 65 northing_max= project_slide.northingmax_nsw, 66 use_cache=True, verbose=True) 62 67 63 68 print 'create offshore' 64 G1 = Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')+\ 65 Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')+\ 66 Geospatial_data(file_name = project.offshore_dem_name3 + '.xya')+\ 67 Geospatial_data(file_name = project.offshore_dem_name4 + '.xya')+\ 68 Geospatial_data(file_name = project.offshore_dem_name5 + '.xya')+\ 69 Geospatial_data(file_name = project.offshore_dem_name6 + '.xya')+\ 70 Geospatial_data(file_name = project.offshore_dem_name7 + '.xya')+\ 71 Geospatial_data(file_name = project.offshore_dem_name8 + '.xya')+\ 72 Geospatial_data(file_name = project.offshore_dem_name9 + '.xya')+\ 73 Geospatial_data(file_name = project.offshore_dem_name10 + '.xya')+\ 74 Geospatial_data(file_name = project.offshore_dem_name11 + '.xya')+\ 75 Geospatial_data(file_name = project.offshore_dem_name12 + '.xya')+\ 76 Geospatial_data(file_name = project.offshore_dem_name13 + '.xya')+\ 77 Geospatial_data(file_name = project.offshore_dem_name14 + '.xya')+\ 78 Geospatial_data(file_name = project.offshore_dem_name15 + '.xya')+\ 79 Geospatial_data(file_name = project.offshore_dem_name16 + '.xya')+\ 80 Geospatial_data(file_name = project.offshore_dem_name17 + '.xya')+\ 81 Geospatial_data(file_name = project.offshore_dem_name18 + '.xya')+\ 82 Geospatial_data(file_name = project.offshore_dem_name19 + '.xya')+\ 83 Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')+\ 84 Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\ 85 Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\ 86 Geospatial_data(file_name = project.offshore_interp_dem_name + '.pts') 69 G11 = Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\ 70 Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya') 71 G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\ 72 Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\ 73 Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\ 74 Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\ 75 Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\ 76 Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya') 87 77 print 'create onshore' 88 G 2 = Geospatial_data(file_name = project.on_offshore_dem_name + '.pts')78 G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts') 89 79 print 'add' 90 G = G1 + G2 + G3 80 G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\ 81 G12.clip(Geospatial_data(project_slide.polyAll)) 82 (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)) 91 83 print 'export points' 92 G.export_points_file(project.combined_dem_name + '.pts') 93 G.export_points_file(project.combined_dem_name + '.xya') 84 G.export_points_file(project_slide.combined_dem_name + '.pts') 85 #G.export_points_file(project_slide.combined_dem_name + '.xya') 86 94 87 95 88 #---------------------------------------------------------------------------- -
anuga_work/production/sydney_2006/project_slide.py
r4056 r4058 31 31 nsw100 = 'dem_utm' 32 32 33 # AHO33 # survey data 34 34 offshore_name1 = 'surveyAreaA' 35 35 offshore_name2 = 'surveyAreaB' 36 36 offshore_name3 = 'surveyAreaC' 37 38 # AHO data 39 offshore_name4 = '1000003611export' 40 offshore_name5 = '1000003613export' 41 offshore_name6 = '1000003614export' 42 offshore_name7 = '1000003627export' 43 offshore_name8 = '1000003628export' 44 offshore_name9 = 'AHDexport' 37 45 38 46 #swollen/ all data output … … 63 71 offshore_dem_name2 = datadir + offshore_name2 64 72 offshore_dem_name3 = datadir + offshore_name3 73 offshore_dem_name4 = datadir + offshore_name4 74 offshore_dem_name5 = datadir + offshore_name5 75 offshore_dem_name6 = datadir + offshore_name6 76 offshore_dem_name7 = datadir + offshore_name7 77 offshore_dem_name8 = datadir + offshore_name8 78 offshore_dem_name9 = datadir + offshore_name9 65 79 combined_dem_name = datadir + 'sydney_combined_elevation' 66 80 … … 134 148 bulli_density = 1.46 135 149 bulli_slope = 4.0 136 bulli_alpha = 126.0 +90.0150 bulli_alpha = 126.0 - 90.0 137 151 138 152 shovel_depth = 968.0 … … 142 156 shovel_density = 1.49 143 157 shove_slope = 4.0 144 shovel_alpha = 118.0 +90.0158 shovel_alpha = 118.0 - 90.0 145 159 146 160 yacaaba_depth = 1119.0 … … 150 164 yacaaba_density = 1.48 151 165 yacaaba_slope = 3.0 152 yacaaba_alpha = 133.0 +90.0166 yacaaba_alpha = 133.0 - 90.0 -
anuga_work/production/sydney_2006/run_sydney_slide.py
r4056 r4058 83 83 84 84 print 'create offshore' 85 G1 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\ 86 Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\ 87 Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya')+\ 85 G11 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\ 86 Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\ 87 Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya') 88 G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\ 89 Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\ 90 Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\ 91 Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\ 92 Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\ 93 Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya') 88 94 print 'create onshore' 89 95 G2 = Geospatial_data(file_name = project_slide.on_offshore25_dem_name + '.pts') … … 91 97 G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts') 92 98 print 'add' 93 G = G1.clip(Geospatial(project_slide.poly_surveyclip)) +\ 99 G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\ 100 G12.clip(Geospatial_data(project_slide.polyAll)) 94 101 G2.clip(Geospatial_data(project_slide.poly_25mclip)) +\ 95 102 G3.clip(Geospatial_data(project_slide.poly_origsyd)) +\ 96 (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)) 103 (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)).clip_outside(Geospatial_data(project_slide.poly_origsyd)) 97 104 print 'export points' 98 105 G.export_points_file(project_slide.combined_dem_name + '.pts') -
anuga_work/production/wollongong_2006/project_slide.py
r4046 r4058 27 27 # onshore data provided by NSW LPI, Krishna merged with existing 28 28 # PMD data to make 100m DEM 29 on_offshore_name = '' 29 on_offshore10_name = 'wollongong_10' 30 nsw100 = 'dem_utm' 30 31 31 # AHO 32 offshore_name1 = 'XY100011610' 32 # survey data 33 offshore_name1 = 'surveyAreaA' 34 35 # AHO data 36 offshore_name4 = '1000003611export' 37 offshore_name5 = '1000003613export' 38 offshore_name6 = '1000003614export' 39 offshore_name7 = '1000003627export' 40 offshore_name8 = '1000003628export' 41 offshore_name9 = 'AHDexport' 33 42 34 43 #swollen/ all data output … … 55 64 on_offshore_dem_name = datadir + on_offshore_name 56 65 offshore_dem_name1 = datadir + offshore_name1 66 offshore_dem_name4 = datadir + offshore_name4 67 offshore_dem_name5 = datadir + offshore_name5 68 offshore_dem_name6 = datadir + offshore_name6 69 offshore_dem_name7 = datadir + offshore_name7 70 offshore_dem_name8 = datadir + offshore_name8 71 offshore_dem_name9 = datadir + offshore_name9 57 72 combined_dem_name = datadir + 'gong_combined_elevation' 58 73 … … 62 77 63 78 # bounding polygon for study area 64 polyAll = read_polygon(polygondir+'extentA .csv')79 polyAll = read_polygon(polygondir+'extentA_up.csv') 65 80 66 81 print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0 … … 70 85 ############################### 71 86 72 # broome digitized polygons 73 poly_gong1 = read_polygon(polygondir+'.csv') 74 poly_gong2 = read_polygon(polygondir+'.csv') 75 poly_gong3 = read_polygon(polygondir+'.csv') 87 # areaA digitized polygons 88 poly_local = read_polygon(polygondir+'local.csv') 89 poly_gong = read_polygon(polygondir+'gong.csv') 76 90 77 print 'Area of local polygon', polygon_area(poly_gong1)/1000000.0 78 print 'Area of close polygon', polygon_area(poly_gong2)/1000000.0 79 print 'Area of coastal polygon', polygon_area(poly_gong3)/1000000.0 91 print 'Area of local polygon', polygon_area(poly_local)/1000000.0 80 92 81 93 ################################################################### … … 83 95 ################################################################### 84 96 97 poly_10mclip = read_polygon(polygondir+'10mclip.csv') 98 poly_surveyclip = read_polygon(polygondir+'surveyclip.csv') 99 100 # clipping for nsw grid for conversion 101 eastingmin_nsw = 272900 102 eastingmax_nsw = 463550 103 northingmin_nsw = 6112950 104 northingmax_nsw = 6269750 105 85 106 # exporting asc grid 86 eastingmin =87 eastingmax =88 northingmin =89 northingmax =107 #eastingmin = 108 #eastingmax = 109 #northingmin = 110 #northingmax = 90 111 91 112 ################################################################### … … 103 124 bulli_density = 1.46 104 125 bulli_slope = 4.0 105 bulli_alpha = 126.0 +90.0126 bulli_alpha = 126.0 - 90.0 106 127 107 128 shovel_depth = 968.0 … … 111 132 shovel_density = 1.49 112 133 shove_slope = 4.0 113 shovel_alpha = 118.0 +90.0134 shovel_alpha = 118.0 - 90.0 114 135 115 136 yacaaba_depth = 1119.0 … … 119 140 yacaaba_density = 1.48 120 141 yacaaba_slope = 3.0 121 yacaaba_alpha = 133.0 +90.0142 yacaaba_alpha = 133.0 - 90.0 122 143 -
anuga_work/production/wollongong_2006/run_gong_slide.py
r4036 r4058 52 52 53 53 # filenames 54 on_offshore_dem_name = project.on_offshore_dem_name 55 meshname = project.meshname+'.msh' 54 on_offshore10_dem_name = project_slide.on_offshore10_dem_name 55 nsw_dem_name = project_slide.nsw_dem_name 56 meshname = project_slide.meshname+'.msh' 56 57 57 58 # creates DEM from asc data 58 convert_dem_from_ascii2netcdf(on_offshore_dem_name, use_cache=True, verbose=True) 59 convert_dem_from_ascii2netcdf(on_offshore10_dem_name, use_cache=True, verbose=True) 60 convert_dem_from_ascii2netcdf(nsw_dem_name, use_cache=True, verbose=True) 59 61 60 62 #creates pts file for onshore DEM 61 dem2pts(on_offshore_dem_name, use_cache=True, verbose=True) 63 dem2pts(on_offshore10_dem_name, use_cache=True, verbose=True) 64 dem2pts(nsw_dem_name, 65 easting_min=project_slide.eastingmin_nsw, 66 easting_max=project_slide.eastingmax_nsw, 67 northing_min=project_slide.northingmin_nsw, 68 northing_max= project_slide.northingmax_nsw, 69 use_cache=True, verbose=True) 62 70 63 71 print 'create offshore' 64 G1 = Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')+\ 65 Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')+\ 66 Geospatial_data(file_name = project.offshore_dem_name3 + '.xya')+\ 67 Geospatial_data(file_name = project.offshore_dem_name4 + '.xya')+\ 68 Geospatial_data(file_name = project.offshore_dem_name5 + '.xya')+\ 69 Geospatial_data(file_name = project.offshore_dem_name6 + '.xya')+\ 70 Geospatial_data(file_name = project.offshore_dem_name7 + '.xya')+\ 71 Geospatial_data(file_name = project.offshore_dem_name8 + '.xya')+\ 72 Geospatial_data(file_name = project.offshore_dem_name9 + '.xya')+\ 73 Geospatial_data(file_name = project.offshore_dem_name10 + '.xya')+\ 74 Geospatial_data(file_name = project.offshore_dem_name11 + '.xya')+\ 75 Geospatial_data(file_name = project.offshore_dem_name12 + '.xya')+\ 76 Geospatial_data(file_name = project.offshore_dem_name13 + '.xya')+\ 77 Geospatial_data(file_name = project.offshore_dem_name14 + '.xya')+\ 78 Geospatial_data(file_name = project.offshore_dem_name15 + '.xya')+\ 79 Geospatial_data(file_name = project.offshore_dem_name16 + '.xya')+\ 80 Geospatial_data(file_name = project.offshore_dem_name17 + '.xya')+\ 81 Geospatial_data(file_name = project.offshore_dem_name18 + '.xya')+\ 82 Geospatial_data(file_name = project.offshore_dem_name19 + '.xya')+\ 83 Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')+\ 84 Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\ 85 Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\ 86 Geospatial_data(file_name = project.offshore_interp_dem_name + '.pts') 72 G11 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya') 73 G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\ 74 Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\ 75 Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\ 76 Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\ 77 Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\ 78 Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya') 87 79 print 'create onshore' 88 G2 = Geospatial_data(file_name = project.on_offshore_dem_name + '.pts') 80 G2 = Geospatial_data(file_name = project_slide.on_offshore10_dem_name + '.pts') 81 G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts') 89 82 print 'add' 90 G = G1 + G2 + G3 83 G = G11.clip(Geospatial(project_slide.poly_surveyclip)) +\ 84 G12.clip(Geospatial_data(project_slide.polyAll)) +\ 85 G2.clip(Geospatial_data(project_slide.poly_10mclip)) +\ 86 (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)).clip_outside(Geospatial_data(project_slide.poly_10mclip)) 91 87 print 'export points' 92 G.export_points_file(project .combined_dem_name + '.pts')93 G.export_points_file(project.combined_dem_name + '.xya')88 G.export_points_file(project_slide.combined_dem_name + '.pts') 89 #G.export_points_file(project_slide.combined_dem_name + '.xya') 94 90 95 91 #----------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.