source: anuga_work/production/australia_ph2/sydney/build_model.py @ 6288

Last change on this file since 6288 was 6288, checked in by jgriffin, 15 years ago

new format scripts added, old ones deleted

File size: 3.7 KB
Line 
1"""Build the elevation data to run a tsunami inundation scenario
2for Busselton, WA, Australia.
3
4Input: elevation data from project.py
5Output: pts file stored in project.topographies_dir
6The run_busselton.py is reliant on the output of this script.
7
8"""
9
10#------------------------------------------------------------------------------
11# Import necessary modules
12#------------------------------------------------------------------------------
13
14# Standard modules
15import os
16
17# Related major packages
18from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf
19from anuga.shallow_water.data_manager import dem2pts
20from anuga.shallow_water.data_manager import start_screen_catcher
21from anuga.shallow_water.data_manager import copy_code_files
22from anuga.geospatial_data.geospatial_data import Geospatial_data
23
24# Application specific imports
25import project   # Definition of file names and polygons
26
27
28#------------------------------------------------------------------------------
29# Copy scripts to time stamped output directory and capture screen
30# output to file
31#------------------------------------------------------------------------------
32copy_code_files(project.output_build,__file__, 
33                os.path.dirname(project.__file__)+os.sep+\
34                project.__name__+'.py' )
35start_screen_catcher(project.output_build)
36
37
38#------------------------------------------------------------------------------
39# Preparation of topographic data
40#
41# Convert ASC 2 DEM 2 PTS using source data and store result in source data
42# Do for coarse and fine data
43# Fine pts file to be clipped to area of interest
44#------------------------------------------------------------------------------
45
46print 'project.bounding_polygon', project.bounding_polygon
47print 'project.combined_elevation_basename', project.combined_elevation_basename
48
49# Create Geospatial data from ASCII files
50geospatial_data = {}
51##for filename in project.ascii_grid_filenames:
52##    absolute_filename = project.topographies_folder + os.sep + filename
53##    convert_dem_from_ascii2netcdf(absolute_filename,
54##                                  basename_out=absolute_filename,
55##                                  use_cache=True,
56##                                  verbose=True)
57##    dem2pts(absolute_filename, use_cache=True, verbose=True)
58##
59##    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
60##                                                verbose=True)
61
62# Create Geospatial data from TXT files
63for filename in project.point_filenames:
64    print filename
65    absolute_filename = project.topographies_folder  + os.sep + str(filename)
66    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
67                                                verbose=True)
68
69
70#-------------------------------------------------------------------------------
71# Combine, clip and export dataset
72#-------------------------------------------------------------------------------
73
74print 'Add geospatial objects' # except', project.offshore_name5
75G = None
76for key in geospatial_data:
77    #if key == project.point_filenames[4] or key == project.ascii_filenames[1]:
78    #    # Skip these files for now
79    #    continue
80   
81    G += geospatial_data[key]
82       
83print 'Clip combined geospatial data'
84##G_clip = G.clip_outside(project.poly_aoi1)
85##G_all = G_clip + geospatial_data[project.point_filenames[4]]
86#G_clipped = G_all.clip(project.poly_all)
87G_clip = G.clip(project.bounding_polygon)
88
89
90print 'Export combined DEM file'
91G_clip.export_points_file(project.combined_elevation + '.pts')
92print 'Do txt version too'
93# Use for comparision in ARC
94G_clip.export_points_file(project.combined_elevation + '.txt')
95
Note: See TracBrowser for help on using the repository browser.