Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/build_busselton.py

    r5829 r6062  
    7676print "creates DEMs from asc data"
    7777convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True)
     78convert_dem_from_ascii2netcdf(offshore_in_dir_name5, basename_out=offshore_dir_name5, use_cache=True, verbose=True)
    7879
    7980# creates pts file for onshore DEM
    8081print "creates pts file for onshore DEM"
    8182dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
     83dem2pts(offshore_dir_name5 ,use_cache=True,verbose=True)
    8284
    8385# create onshore pts files
     
    101103G_off4 = Geospatial_data(file_name = offshore_in_dir_name4)
    102104print'create Geospatial data9 objects from topographies'
    103 G_off5 = Geospatial_data(file_name = offshore_in_dir_name5)
     105G_off5 = Geospatial_data(file_name = offshore_dir_name5 + '.pts')
    104106
    105107
     
    109111
    110112print'add all geospatial objects'
    111 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5
     113G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4
    112114
    113115print'clip combined geospatial object by bounding polygon'
    114 G_clipped = G.clip(project.poly_all)
     116G_clip = G.clip_outside(project.poly_aoi1)
     117G_all = G_clip + G_off5
     118G_clipped = G_all.clip(project.poly_all)
    115119
    116120print'export combined DEM file'
     
    118122    mkdir (project.topographies_dir)
    119123G_clipped.export_points_file(project.combined_dir_name + '.pts')
    120 #G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
     124G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
    121125
Note: See TracChangeset for help on using the changeset viewer.