Changeset 6280

Timestamp:

Feb 5, 2009, 3:06:41 PM (16 years ago)

Author:

ole

Message:

Notes on combining Geospatial data

File:

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

anuga_work/production/busselton/standardised_version/build_busselton.py

@@ -30 +30 @@
 # output to file
 #------------------------------------------------------------------------------
-copy_code_files(project.output_build,__file__, 
-                os.path.dirname(project.__file__)+os.sep+\
-                project.__name__+'.py' )
-start_screen_catcher(project.output_build)
+#copy_code_files(project.output_build,__file__, 
+#                os.path.dirname(project.__file__)+os.sep+\
+#                project.__name__+'.py' )
+#start_screen_catcher(project.output_build)
 
 
@@ -72 +72 @@
 
 print 'Add geospatial objects' # except', project.offshore_name5
-G = geospatial_data
-##None
-##for key in geospatial_data:
-##    if key != project.offshore_name5:
-##        G += geospatial_data[key]
+G = None
+for key in geospatial_data:
+    #if key == project.point_filenames[4] or key == project.ascii_filenames[1]:
+    #    # Skip these files for now
+    #    continue
+    
+    G += geospatial_data[key]
         
 print 'Clip combined geospatial data'
 ##G_clip = G.clip_outside(project.poly_aoi1)
-##G_all = G_clip + geospatial_data[project.offshore_name5]   # G_off5
+##G_all = G_clip + geospatial_data[project.point_filenames[4]] 
 #G_clipped = G_all.clip(project.poly_all)
 G_clip = G.clip(project.bounding_polygon)