Changeset 5811
- Timestamp:
- Oct 2, 2008, 4:07:08 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/Broome_2008/build_broome.py
r5810 r5811 92 92 #creates pts file for onshore DEM 93 93 print "creates pts file for onshore DEM" 94 dem2pts(o ffshore_dir_name ,use_cache=True,verbose=True)94 dem2pts(onshore_dir_name ,use_cache=True,verbose=True) 95 95 dem2pts(offshore_dir_name1 ,use_cache=True,verbose=True) 96 96 dem2pts(offshore_dir_name2 ,use_cache=True,verbose=True) … … 120 120 G_off = Geospatial_data(file_name = offshore_in_dir_name) 121 121 print'create Geospatial data8 objects from topographies' 122 G_off1 = Geospatial_data(file_name = offshore_in_dir_name1 )122 G_off1 = Geospatial_data(file_name = offshore_in_dir_name1 + '.pts') 123 123 print'create Geospatial data9 objects from topographies' 124 G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 )124 G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.pts') 125 125 print'create Geospatial data10 objects from topographies' 126 G_off3 = Geospatial_data(file_name = offshore_in_dir_name3 )126 G_off3 = Geospatial_data(file_name = offshore_in_dir_name3 + '.pts') 127 127 print'create Geospatial data10 objects from topographies' 128 G_off 3 = Geospatial_data(file_name = offshore_in_dir_name4)128 G_off4 = Geospatial_data(file_name = offshore_in_dir_name4 + '.pts') 129 129 print'create Geospatial data10 objects from topographies' 130 G_off 3 = Geospatial_data(file_name = offshore_in_dir_name5)130 G_off5 = Geospatial_data(file_name = offshore_in_dir_name5 + '.pts') 131 131 print'create Geospatial data10 objects from topographies' 132 G_off 3 = Geospatial_data(file_name = offshore_in_dir_name6)132 G_off6 = Geospatial_data(file_name = offshore_in_dir_name6 + '.pts') 133 133 print'add all geospatial objects' 134 G = G1 + G_coast + G_off + G_off1 + G_off2 + G_off3 134 G = G1 + G_coast + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5 + G_off6 135 135 print'clip combined geospatial object by bounding polygon' 136 136 G_clipped = G.clip(project.poly_all)
Note: See TracChangeset
for help on using the changeset viewer.