Changeset 6062 for anuga_work/production/busselton/build_busselton.py
- Timestamp:
- Dec 11, 2008, 4:10:06 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/busselton/build_busselton.py
r5829 r6062 76 76 print "creates DEMs from asc data" 77 77 convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True) 78 convert_dem_from_ascii2netcdf(offshore_in_dir_name5, basename_out=offshore_dir_name5, use_cache=True, verbose=True) 78 79 79 80 # creates pts file for onshore DEM 80 81 print "creates pts file for onshore DEM" 81 82 dem2pts(onshore_dir_name ,use_cache=True,verbose=True) 83 dem2pts(offshore_dir_name5 ,use_cache=True,verbose=True) 82 84 83 85 # create onshore pts files … … 101 103 G_off4 = Geospatial_data(file_name = offshore_in_dir_name4) 102 104 print'create Geospatial data9 objects from topographies' 103 G_off5 = Geospatial_data(file_name = offshore_ in_dir_name5)105 G_off5 = Geospatial_data(file_name = offshore_dir_name5 + '.pts') 104 106 105 107 … … 109 111 110 112 print'add all geospatial objects' 111 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5113 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 112 114 113 115 print'clip combined geospatial object by bounding polygon' 114 G_clipped = G.clip(project.poly_all) 116 G_clip = G.clip_outside(project.poly_aoi1) 117 G_all = G_clip + G_off5 118 G_clipped = G_all.clip(project.poly_all) 115 119 116 120 print'export combined DEM file' … … 118 122 mkdir (project.topographies_dir) 119 123 G_clipped.export_points_file(project.combined_dir_name + '.pts') 120 #G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC124 G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC 121 125
Note: See TracChangeset
for help on using the changeset viewer.