Changeset 6499 for anuga_work


Ignore:
Timestamp:
Mar 12, 2009, 1:40:40 PM (16 years ago)
Author:
kristy
Message:

reading multiple internal polygons from the one csv file

Location:
anuga_work/production/hobart_2009
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/hobart_2009/project.py

    r6457 r6499  
    5959# Used in build_elevation.py
    6060# Format for ascii grids, as produced in ArcGIS + a projection file
    61 ascii_grid_filenames = ['hob3_project',   # Topo
    62                         'grid_250m_project'] # Busselton Topo
     61ascii_grid_filenames = ['hob3_pro_extract', 'grid_250m_project', 'ldr_bbay',
     62                        'ldr_bellerive', 'ldr_brunyislandneck',
     63                        'ldr_calvertsbeach', 'ldr_carltonbeach',
     64                        'ldr_connelly', 'ldr_cremorne', 'ldr_dodge1',
     65                        'ldr_dodge2', 'ldr_kingston', 'ldr_lauderdale',
     66                        'ldr_linersfarne', 'ldr_oppossumbay', 'ldr_primrose',
     67                        'ldr_rosny', 'ldr_sandybay', 'ldr_selfspoint',
     68                        'ldr_sevenmile', 'ldr_southarm1', 'ldr_southarm2',
     69                        'ldr_sullivanscove', 'ldr_tranmere', 'tasdem_mask']
    6370
    6471# Format for point is x,y,elevation (with header)
     
    7582# Format for points easting,northing (no header)
    7683bounding_polygon_filename = 'bounding_polygon.csv'
    77 bounding_polygon_maxarea = 100000
     84bounding_polygon_maxarea = 1000000
    7885
    7986# INTERIOR REGIONS -  for designing the mesh
    8087# Used in run_model.py
    8188# Format for points easting,northing (no header)                   
    82 interior_regions_data = []
     89interior_regions_data = [['aos1.csv', 1500],
     90                         ['aos2.csv', 1500],
     91                         ['sw.csv', 30000]]
     92PriorityArea_filename = 'PriorityAreas.csv'
    8393
    8494# LAND - used to set the initial stage/water to be offcoast only
     
    237247mux_input = join(event_folder, mux_input_filename)
    238248
     249#Multiple polygons in one CSV file to make internal polygons
     250PriorityAreas = join(polygons_folder, PriorityArea_filename)
     251
     252
  • anuga_work/production/hobart_2009/setup_model.py

    r6457 r6499  
    1111from os.path import join, exists
    1212from anuga.utilities.polygon import read_polygon, number_mesh_triangles
    13 
     13from anuga.shallow_water.data_manager import csv2polygons
    1414import project
    1515
     
    157157# Create list of interior polygons with scaling factor
    158158project.interior_regions = []
     159
     160polygons, maxareas = csv2polygons(project.PriorityAreas)
     161print 'Creating %d internal polygons' % len(polygons)
     162#def create_polygon_function(polygons, geo_reference=None):
     163project.interior_regions = []
     164for i, key in enumerate(polygons):
     165    if i%100==0: print i
     166    poly = polygons[key]
     167    maxarea = float(maxareas[key])
     168    project.interior_regions.append([poly,
     169                                     maxarea*project.scale_factor])
     170
    159171for filename, maxarea in project.interior_regions_data:
    160172    polygon = read_polygon(join(project.polygons_folder, filename))
Note: See TracChangeset for help on using the changeset viewer.