Changeset 6901

Timestamp:

Apr 24, 2009, 4:16:26 PM (14 years ago)

Author:

myall

Message:

making rasters and images

Location:

anuga_work/production/australia_ph2

Files:

• Communities_model_tracking.xls (modified) (previous)
• bamaga/Arc_asc2raster_GDA94z50.py (added)
• broome/Arc_asc2raster_GDA94z50.py (added)
• dampier/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)
• derby/Arc_asc2raster_GDA94z50.py (added)

anuga_work/production/australia_ph2/dampier/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090324_140112_run_final_0_70385_Bf_all_kvanputt'
+time_dir1 = '20090415_124319_run_final_0_70844_1307_Tb__kvanputt'
+time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
+time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
 
-time_dirs = [time_dir1]
+time_dirs = [time_dir1, time_dir2, time_dir3]
 
 for time_dir in time_dirs:
-
     # Local variables...
     folder = scenario_dir + output_dir +  time_dir +'\\'
@@ -46 +47 @@
     #replication dictionary
     replicate = (('dampier', ''),('_', ''),('max','_M'),
-                 ('CBD', 'CDB'),
+                 ('CBD', 'CDB'),('All',''),
                  ('depth','_depth'),('speed', '_speed'),
                  ('elevation', '_ele_'), ('stage','_stage'))
@@ -68 +69 @@
         gp.ASCIIToRaster_conversion(infile, output_DEM, "FLOAT")
 
-        print 'Process: Define Projection' 
+        print 'Process: Define Projection'
+        #GDA_1994_MGA_Zone_54
         gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
                                                    ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"