Changeset 6971


Ignore:
Timestamp:
May 5, 2009, 5:45:54 PM (11 years ago)
Author:
myall
Message:

converting ascii to raster and making images

Location:
anuga_work/production/australia_ph2
Files:
3 added
9 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/broome/Arc_asc2raster_GDA94z50.py

    r6901 r6971  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090408_152311_run_final_0_27337_1661_Tb__mhingee'
    28 time_dir2 = '20090409_070404_run_final_0_70336_1661_Tb__mhingee'
     27##time_dir1 = '20090408_152311_run_final_0_27337_1661_Tb__mhingee'
     28##time_dir2 = '20090409_070404_run_final_0_70336_1661_Tb__mhingee'
    2929time_dir3 = '20090409_224602_run_final_0_70869_1661_Tb__mhingee'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dir1 = '20090421_174625_run_final_0_70869_1661_Tb_internal_mhingee'
     32
     33time_dirs = [time_dir1, time_dir3]#, time_dir3]
    3234
    3335for time_dir in time_dirs:
     
    3840##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    3941
    40     print 'Process: Create File GDB'
    41     gp.CreateFileGDB_management(folder, "raster")
     42##    print 'Process: Create File GDB'
     43##    gp.CreateFileGDB_management(folder, "raster")
    4244
    4345    gp.Workspace = raster_gbd
     
    5254
    5355    generate_filename = []
    54     input_ascii = glob.glob(folder + '*max.asc')
     56    input_ascii = glob.glob(folder + '*elevation.asc')
    5557
    5658    for infile in input_ascii:
  • anuga_work/production/australia_ph2/carnarvon/project.py

    r6949 r6971  
    5151
    5252
    53 internal_polygon = False
     53internal_polygon = True
    5454
    5555#-------------------------------------------------------------------------------
  • anuga_work/production/australia_ph2/dampier/Arc_asc2raster_GDA94z50.py

    r6901 r6971  
    2626
    2727time_dir1 = '20090415_124319_run_final_0_70844_1307_Tb__kvanputt'
    28 time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
    29 time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
     28##time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
     29##time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
     32
     33time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3234
    3335for time_dir in time_dirs:
     
    3840##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    3941
    40     print 'Process: Create File GDB'
    41     gp.CreateFileGDB_management(folder, "raster")
     42##    print 'Process: Create File GDB'
     43##    gp.CreateFileGDB_management(folder, "raster")
    4244
    4345    gp.Workspace = raster_gbd
     
    5254
    5355    generate_filename = []
    54     input_ascii = glob.glob(folder + '*max.asc')
     56    input_ascii = glob.glob(folder + '*elevation.asc')
    5557
    5658    for infile in input_ascii:
  • anuga_work/production/australia_ph2/derby/Arc_asc2raster_GDA94z50.py

    r6901 r6971  
    2626
    2727time_dir1 = '20090408_171910_run_final_0_27268_1708_Tb__kvanputt'
    28 time_dir2 = '20090410_070302_run_final_0_70437_1708_Tb__kvanputt'
    29 time_dir3 = '20090411_153927_run_final_0_70862_1708_Tb__kvanputt'
     28##time_dir2 = '20090410_070302_run_final_0_70437_1708_Tb__kvanputt'
     29##time_dir3 = '20090411_153927_run_final_0_70862_1708_Tb__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dir2 = '20090421_175503_run_final_0_27268_1708_Tb_internal_mhingee'
     32
     33time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3234
    3335for time_dir in time_dirs:
     
    3840##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    3941
    40     print 'Process: Create File GDB'
    41     gp.CreateFileGDB_management(folder, "raster")
     42##    print 'Process: Create File GDB'
     43##    gp.CreateFileGDB_management(folder, "raster")
    4244
    4345    gp.Workspace = raster_gbd
     
    5254
    5355    generate_filename = []
    54     input_ascii = glob.glob(folder + '*max.asc')
     56    input_ascii = glob.glob(folder + '*elevation.asc')
    5557
    5658    for infile in input_ascii:
  • anuga_work/production/australia_ph2/geraldton/Arc_asc2raster_GDA94z50.py

    r6883 r6971  
    2626
    2727time_dir1 = '20090414_114540_run_final_0_64475_576_Tb__kvanputt'
    28 time_dir2 = '20090415_014342_run_final_0_70921_576_Tb__kvanputt'
    29 time_dir3 = '20090415_153319_run_final_0_26985_576_Tb__kvanputt'
     28##time_dir2 = '20090415_014342_run_final_0_70921_576_Tb__kvanputt'
     29##time_dir3 = '20090415_153319_run_final_0_26985_576_Tb__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dir2 = '20090421_172030_run_final_0_64475_576_Tb_internal_mhingee'
     32
     33time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3234
    3335for time_dir in time_dirs:
     
    3840##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    3941
    40     print 'Process: Create File GDB'
    41     gp.CreateFileGDB_management(folder, "raster")
     42##    print 'Process: Create File GDB'
     43##    gp.CreateFileGDB_management(folder, "raster")
    4244
    4345    gp.Workspace = raster_gbd
     
    5254
    5355    generate_filename = []
    54     input_ascii = glob.glob(folder + '*max.asc')
     56    input_ascii = glob.glob(folder + '*elevation.asc')
    5557
    5658    for infile in input_ascii:
  • anuga_work/production/australia_ph2/melbourne_east/Arc_asc2raster_GDA94z50.py

    r6959 r6971  
    3434time_dir7 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
    3535
    36 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir7]
    37 ##time_dirs = [time_dir4]
     36##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir7]
     37time_dirs = [time_dir1, time_dir7]
    3838   
    3939for time_dir in time_dirs:
     
    4444##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4545
    46     print 'Process: Create File GDB'
    47     gp.CreateFileGDB_management(folder, "raster")
     46##    print 'Process: Create File GDB'
     47##    gp.CreateFileGDB_management(folder, "raster")
    4848
    4949    gp.Workspace = raster_gbd
     
    5858
    5959    generate_filename = []
    60     input_ascii = glob.glob(folder + '*max.asc')
     60    input_ascii = glob.glob(folder + '*elevation.asc')
    6161
    6262    for infile in input_ascii:
  • anuga_work/production/australia_ph2/onslow/Arc_asc2raster_GDA94z50.py

    r6883 r6971  
    2626
    2727time_dir1 = '20090414_121216_run_final_0_27279_711_Tb__kvanputt'
    28 time_dir2 = '20090415_001215_run_final_0_70371_711_Tb__kvanputt'
    29 time_dir3 = '20090415_114548_run_final_0_27168_711_Tb__kvanputt'
     28##time_dir2 = '20090415_001215_run_final_0_70371_711_Tb__kvanputt'
     29##time_dir3 = '20090415_114548_run_final_0_27168_711_Tb__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dir2 = '20090421_172702_run_final_0_27279_711_Tb_internal_mhingee'
     32
     33time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3234
    3335for time_dir in time_dirs:
     
    3840##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    3941
    40     print 'Process: Create File GDB'
    41     gp.CreateFileGDB_management(folder, "raster")
     42##    print 'Process: Create File GDB'
     43##    gp.CreateFileGDB_management(folder, "raster")
    4244
    4345    gp.Workspace = raster_gbd
     
    5254
    5355    generate_filename = []
    54     input_ascii = glob.glob(folder + '*max.asc')
     56    input_ascii = glob.glob(folder + '*elevation.asc')
    5557
    5658    for infile in input_ascii:
  • anuga_work/production/australia_ph2/strahan/Arc_asc2raster_GDA94z50.py

    r6959 r6971  
    3232
    3333##time_dirs = [time_dir1, time_dir2, time_dir3]
    34 time_dirs = [time_dir4]
     34time_dirs = [time_dir1]
    3535   
    3636for time_dir in time_dirs:
     
    4141##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4242
    43     print 'Process: Create File GDB'
    44     gp.CreateFileGDB_management(folder, "raster")
     43##    print 'Process: Create File GDB'
     44##    gp.CreateFileGDB_management(folder, "raster")
    4545
    4646    gp.Workspace = raster_gbd
     
    5555
    5656    generate_filename = []
    57     input_ascii = glob.glob(folder + '*.asc')
     57    input_ascii = glob.glob(folder + '*elevation.asc')
    5858
    5959    for infile in input_ascii:
Note: See TracChangeset for help on using the changeset viewer.