Changeset 6273


Ignore:
Timestamp:
Feb 4, 2009, 2:15:57 PM (10 years ago)
Author:
kristy
Message:

updating

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/patong/asc2raster_WGS84z47N.py

    r6066 r6273  
    2626
    2727
    28 time_dir1 = '20081126_141801_run_final_0.8_kvanputt'
     28time_dir1 = '20090130_135044_run_final_0.8_buildings_onielsen'
     29time_dir2 = '20090123_135355_run_final_0.8_buildings_onielsen'
    2930
    30 time_dirs = [time_dir1]
     31time_dirs = [time_dir1] #, time_dir2]
    3132
    3233for time_dir in time_dirs:
     
    3536    folder = scenario_dir + output_dir + time_dir + '\\'
    3637    raster_gbd = folder + 'raster.gdb'
    37    # land = scenario_dir + "map_work\\Patong.gdb\\poly_all_cut"
     38    land = scenario_dir + "map_work\\patong.gdb\\initial_conditions_Dissolve"
    3839   
    39     print 'Process: Create File GDB'
    40     gp.CreateFileGDB_management(folder, "raster")
     40##    print 'Process: Create File GDB'
     41##    gp.CreateFileGDB_management(folder, "raster")
    4142
    4243    gp.Workspace = raster_gbd
     
    4950
    5051    generate_filename = []
    51     input_ascii = glob.glob(folder + '*.asc')
     52    input_ascii = glob.glob(folder + '*speed.asc')
    5253
    5354    for infile in input_ascii:
     
    6162        generate_filename.append(output_DEM)
    6263        print 'Output DEM ',output_DEM
    63         output_extract = output_DEM + 'E'
    64         print 'Output Extract ',output_extract
    6564       
    6665        print 'Process: ASCII to Raster'
     
    7574                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    7675
    77 
    78 
     76##        output_extract = output_DEM + 'E'
     77##        print 'Output Extract ',output_extract
    7978##        print 'Process: Extract by Mask'
    8079##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
Note: See TracChangeset for help on using the changeset viewer.