Changeset 6062 for anuga_work/production/busselton/build_busselton.py

Timestamp:

Dec 11, 2008, 4:10:06 PM (16 years ago)

Author:

kristy

Message:

updated scripts for new data in bunbury and marina in busselton

File:

• anuga_work/production/busselton/build_busselton.py (modified) (4 diffs)

anuga_work/production/busselton/build_busselton.py

@@ -76 +76 @@
 print "creates DEMs from asc data"
 convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True)
+convert_dem_from_ascii2netcdf(offshore_in_dir_name5, basename_out=offshore_dir_name5, use_cache=True, verbose=True)
 
 # creates pts file for onshore DEM
 print "creates pts file for onshore DEM"
 dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
+dem2pts(offshore_dir_name5 ,use_cache=True,verbose=True)
 
 # create onshore pts files
@@ -101 +103 @@
 G_off4 = Geospatial_data(file_name = offshore_in_dir_name4)
 print'create Geospatial data9 objects from topographies'
-G_off5 = Geospatial_data(file_name = offshore_in_dir_name5)
+G_off5 = Geospatial_data(file_name = offshore_dir_name5 + '.pts')
 
 
@@ -109 +111 @@
 
 print'add all geospatial objects'
-G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5 
+G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4
 
 print'clip combined geospatial object by bounding polygon'
-G_clipped = G.clip(project.poly_all)
+G_clip = G.clip_outside(project.poly_aoi1)
+G_all = G_clip + G_off5
+G_clipped = G_all.clip(project.poly_all)
 
 print'export combined DEM file'
@@ -118 +122 @@
     mkdir (project.topographies_dir)
 G_clipped.export_points_file(project.combined_dir_name + '.pts')
-#G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
+G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC