source: anuga_work/production/busselton/Arc_txt2feature_class.py @ 6201

Last change on this file since 6201 was 6201, checked in by kristy, 15 years ago

Updating scripts

File size: 2.8 KB
Line 
1# ---------------------------------------------------------------------------
2# txt2fc_field.py
3# Created on: Thu Nov 27 2008 01:48:26 PM
4#   (generated by ArcGIS/ModelBuilder)
5# ---------------------------------------------------------------------------
6
7# Import system modules
8import sys, string, os, arcgisscripting, glob, os.path
9
10# Create the Geoprocessor object
11gp = arcgisscripting.create()
12
13# Load required toolboxes...
14gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")
15gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
16gp.overwriteoutput = 1
17
18scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\busselton_tsunami_scenario\\"
19output_dir = "anuga\\outputs\\"
20
21time_dir1 = '20081211_154006_run_final_0.6_27255_alpha0.1_kvanputt'
22##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
23##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
24##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt'
25##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt'
26##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt'
27
28time_dirs = [time_dir1] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
29
30for time_dir in time_dirs:
31
32    # Local variables...
33    folder =  scenario_dir + output_dir + time_dir + '\\'
34    input_dir_name = folder+'busselton_runup.csv'
35    input_name = os.path.basename(input_dir_name)
36    raster_gdb = folder + "raster.gdb"
37    layer = 'layer'
38    feature_class_name = input_name[:-4]
39
40##    print 'Process: Create File GDB'
41##    gp.CreateFileGDB_management(folder, "raster")
42
43    print 'Process: Make XY Event Layer'
44    gp.MakeXYEventLayer_management(input_dir_name, "x", "y", layer, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM"
45                                   "['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT"
46                                   "['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting'"
47                                   ",500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER"
48                                   "['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]];IsHighPrecision")
49
50    print 'Process: Feature Class to Feature Class'
51    gp.FeatureClassToFeatureClass_conversion(layer, raster_gdb, "busselton_runup", "", "area area true true false 255 Text 0 0 ,First,#,busselton_runup.csv_Features,area,-1,-1;q q true true false 8 Double 0 0 ,First,#,busselton_runup.csv_Features,q,-1,-1;x x true true false 8 Double 0 0 ,First,#,busselton_runup.csv_Features,x,-1,-1;y y true true false 8 Double 0 0 ,First,#,busselton_runup.csv_Features,y,-1,-1", "")
52
53
54 
Note: See TracBrowser for help on using the repository browser.