source: anuga_work/production/perth/Arc_txt2feature_class_res.py @ 6929

# ---------------------------------------------------------------------------
# txt2fc_res.py
# Created on: Mon Dec 01 2008 11:36:17 AM
#   (generated by ArcGIS/ModelBuilder)
# ---------------------------------------------------------------------------

# Import system modules
import sys, string, os, arcgisscripting
from os.path import join

# Create the Geoprocessor object
gp = arcgisscripting.create()

# Load required toolboxes...
gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")
gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
gp.overwriteoutput = 1

scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\perth_tsunami_scenario\\"
output_dir = "anuga\\outputs\\"

time_dir1 = '20090324_135427_run_final_0.6_27255_Bt_kvanputt'
time_dir2 = '20090324_154957_run_final_0.6_68693_Bt_kvanputt'
time_dir3 = '20090324_155028_run_final_0.6_27283_Bt_kvanputt'
time_dir4 = '20090327_155353_run_final_0_27283_Bt_kvanputt'
time_dir5 = '20090327_155659_run_final_0_27255_Bt_kvanputt'
time_dir6 = '20090327_155803_run_final_0_68693_Bt_kvanputt'

time_dirs = [time_dir3, time_dir4, time_dir5, time_dir6] 

def replace_in_file(filename, instr, repstr):
    '''Replace one string with another in a file.'''
    
    # get the file contents
    fd = open(filename, 'r')
    lines = fd.readlines()
    fd.close()

    # convert list to one string
    newlines = ''.join(lines)
    
    # replace 'instr' with 'repstr'
    newlines = newlines.replace(instr, repstr)
    
    # write changed lines back to file
    fd = open(filename, 'w')
    fd.write(newlines)
    fd.close()

for time_dir in time_dirs:
    print 'time_dir = ', time_dir
    
    # Local variables...
    folder =  scenario_dir + output_dir + time_dir + '\\'
    input_dir_name = folder+'perth_res.csv'
    raster_gdb = folder + "raster.gdb"
    layer = 'layer'

    print 'replace _$ with nothing'
    replace_in_file(input_dir_name, '_$', '')
    
##    print 'Process: Feature Class to Feature Class'
##    gp.FeatureClassToFeatureClass_conversion(layer, raster_gdb, "perth_res", "", "OBJECTID OBJECTID true true false 4 Long 0 10 ,First,#,perth_res.csv_Features,OBJECTID,-1,-1;LONGITUDE LONGITUDE true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,LONGITUDE,-1,-1;LATITUDE LATITUDE true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,LATITUDE,-1,-1;STRUCT_INUNDATED STRUCT_INUNDATED true true false 4 Long 0 10 ,First,#,perth_res.csv_Features,STRUCT_INUNDATED,-1,-1;STRUCT_DAMAGE_fraction STRUCT_DAMAGE_fraction true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,STRUCT_DAMAGE_fraction,-1,-1;COLLAPSE_CSV_INFO COLLAPSE_CSV_INFO true true false 255 Text 0 0 ,First,#,perth_res.csv_Features,COLLAPSE_CSV_INFO,-1,-1;MAX_DEPTH_m MAX_DEPTH_m true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,MAX_DEPTH_m,-1,-1;STRUCT_COLLAPSED STRUCT_COLLAPSED true true false 255 Text 0 0 ,First,#,perth_res.csv_Features,STRUCT_COLLAPSED,-1,-1;STRUCT_LOSS STRUCT_LOSS true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,STRUCT_LOSS,-1,-1;CONTENTS_DAMAGE_fraction CONTENTS_DAMAGE_fraction true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,CONTENTS_DAMAGE_fraction,-1,-1;CONTENTS_LOSS CONTENTS_LOSS true true false 8 Double 0 0 ,First,#,perth_res.csv_Features,CONTENTS_LOSS,-1,-1", "")
##