1 | """ |
---|
2 | Script for building the elevation data to run a tsunami inundation scenario |
---|
3 | for busselton, WA, Australia. |
---|
4 | |
---|
5 | Input: elevation data from project_250m.py |
---|
6 | Output: pts file stored in project_250m.topographies_dir |
---|
7 | The run_busselton.py is reliant on the output of this script. |
---|
8 | |
---|
9 | """ |
---|
10 | |
---|
11 | #------------------------------------------------------------------------------ |
---|
12 | # Import necessary modules |
---|
13 | #------------------------------------------------------------------------------ |
---|
14 | |
---|
15 | # Standard modules |
---|
16 | from os import sep |
---|
17 | from os.path import dirname, basename |
---|
18 | from os import mkdir, access, F_OK |
---|
19 | from shutil import copy |
---|
20 | import time |
---|
21 | import sys |
---|
22 | |
---|
23 | # Related major packages |
---|
24 | from anuga.shallow_water import Domain |
---|
25 | from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts |
---|
26 | from anuga.geospatial_data.geospatial_data import * |
---|
27 | from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files |
---|
28 | |
---|
29 | # Application specific imports |
---|
30 | import project # Definition of file names and polygons |
---|
31 | |
---|
32 | #------------------------------------------------------------------------------ |
---|
33 | # Copy scripts to time stamped output directory and capture screen |
---|
34 | # output to file |
---|
35 | #------------------------------------------------------------------------------ |
---|
36 | |
---|
37 | copy_code_files(project_250m.output_build_time_dir,__file__, |
---|
38 | dirname(project_250m.__file__)+sep+ project_250m.__name__+'.py' ) |
---|
39 | |
---|
40 | start_screen_catcher(project_250m.output_build_time_dir) |
---|
41 | |
---|
42 | print 'USER: ', project_250m.user |
---|
43 | |
---|
44 | #------------------------------------------------------------------------------- |
---|
45 | # Preparation of topographic data |
---|
46 | # |
---|
47 | # Convert ASC 2 DEM 2 PTS using source data and store result in source data |
---|
48 | # Do for coarse and fine data |
---|
49 | # Fine pts file to be clipped to area of interest |
---|
50 | #------------------------------------------------------------------------------- |
---|
51 | print "project_250m.poly_all", project_250m.poly_all |
---|
52 | print "project_250m.combined_dir_name", project_250m.combined_dir_name |
---|
53 | |
---|
54 | # input elevation directory filenames |
---|
55 | onshore_in_dir_name = project_250m.onshore_in_dir_name |
---|
56 | |
---|
57 | # output elevation directory filenames |
---|
58 | onshore_dir_name = project_250m.onshore_dir_name |
---|
59 | |
---|
60 | # creates DEM from asc data |
---|
61 | print "creates DEMs from asc data" |
---|
62 | convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True) |
---|
63 | |
---|
64 | # creates pts file for onshore DEM |
---|
65 | print "creates pts file for onshore DEM" |
---|
66 | dem2pts(onshore_dir_name ,use_cache=True,verbose=True) |
---|
67 | |
---|
68 | # create onshore pts files |
---|
69 | print'create Geospatial data1 objects from topographies' |
---|
70 | G = Geospatial_data(file_name = onshore_dir_name + '.pts') |
---|
71 | |
---|
72 | #------------------------------------------------------------------------------- |
---|
73 | # Combine, clip and export dataset |
---|
74 | #------------------------------------------------------------------------------- |
---|
75 | |
---|
76 | print'clip combined geospatial object by bounding polygon' |
---|
77 | G_clipped = G.clip(project_250m.poly_all) |
---|
78 | |
---|
79 | print'export combined DEM file' |
---|
80 | if access(project_250m.topographies_dir,F_OK) == 0: |
---|
81 | mkdir (project_250m.topographies_dir) |
---|
82 | G_clipped.export_points_file(project_250m.combined_dir_name + '.pts') |
---|
83 | #G_clipped.export_points_file(project_250m.combined_dir_name + '.txt') #Use for comparision in ARC |
---|
84 | |
---|