- Timestamp:
- May 6, 2009, 4:06:54 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/maningrida/Arc_asc2raster_GDA94z53.py
r6969 r6985 30 30 time_dir4 = '20090422_115416_run_final_0_64970_1878_Tb_internal_mhingee' 31 31 32 ##time_dirs = [time_dir2, time_dir3, time_dir4] 33 time_dirs = [time_dir4] 32 34 33 events = [[time_dir1,0.430643],[time_dir2,0.462943],[time_dir3,0.418183],[time_dir4,0.418183]] 34 35 for event in events: 36 ##for time_dir in time_dirs: 37 time_dir = event[0] 38 max_wave = event[1] 35 for time_dir in time_dirs: 39 36 # Local variables... 40 37 folder = scenario_dir + output_dir + time_dir +'\\' … … 44 41 45 42 46 print 'Process: Create File GDB'47 gp.CreateFileGDB_management(folder, "raster")43 ## print 'Process: Create File GDB' 44 ## gp.CreateFileGDB_management(folder, "raster") 48 45 49 46 gp.Workspace = raster_gbd … … 52 49 53 50 #replication dictionary 54 replicate = (('maningrida', ''),('_', ''), 55 ('depth','_depth'),('speed', '_speed'), 51 replicate = (('maningrida', ''),('_', ''),('All',''), 52 ('depth','_depth'),('speed', '_speed'),('max','_M'), 56 53 ('elevation', '_ele_'), ('stage','_stage')) 57 54 58 55 generate_filename = [] 59 input_ascii = glob.glob(folder + '* .asc')56 input_ascii = glob.glob(folder + '*elevation.asc') 60 57 61 58 for infile in input_ascii: … … 75 72 76 73 print 'Process: Define Projection' 77 gp.DefineProjection_management(output_DEM, "PROJCS[' GDA_1994_MGA_Zone_53',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"78 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',13 5.0],PARAMETER['Scale_Factor',0.9996]"74 gp.DefineProjection_management(output_DEM, "PROJCS['CM_133',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" 75 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',133.0],PARAMETER['Scale_Factor',0.9996]" 79 76 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 80 77 ## output_extract = output_DEM + '_E' … … 82 79 ## print 'Process: Extract by Mask' 83 80 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 84 # do this bit only if there are only stage asc files85 div_file = output_DEM.replace('stage','div_stage')86 print 'divide', output_DEM,' by', max_wave,' and call',div_file87 gp.Divide_sa(output_DEM,max_wave,div_file)88 81 89 82 90 83 84
Note: See TracChangeset
for help on using the changeset viewer.