def PalumboPlots(): DataDirectory = "T://analysis_for_papers//Cosmo_paper//Palumbo_shield_plots//" Filename = "SpawnedBasin_07C13-(Q8).bil" Drapename = "SpawnedBasin_07C13-(Q8)_BASINS.bil" ThisFile = DataDirectory+Filename DrapeFile = DataDirectory+Drapename LSDP.BasicDrapedPlotGridPlot(ThisFile, DrapeFile, thiscmap='terrain',drape_cmap='gray', colorbarlabel='Elevation in meters',clim_val = (2000,5000), drape_alpha = 0.4,FigFileName = 'Basin.pdf',FigFormat = 'pdf')
def TestNewMappingTools(): #DataDirectory = "C://basin_data//Chile//lat26p0//" #Filename = "Site_lat26p0_UTM19_DEM_FILL.bil" #DataDirectory = "C://basin_data//Model_results//June2015_Results//HighK//" #DataDirectory = "T://test_clone//topodata//" #DataDirectory = "T://analysis_for_papers//Cosmo_paper//Tibet//for_plotting//" DataDirectory = "C://basin_data//CosmoPaper//DEMs//" #Filename = "SanBern.bil" Filename = "SpawnedBasin_07C13-(Q8)_SH.bil" #Filename = "CRNvariable_long_0_0_1var128.asc" #DrapeFileName = "CRNvariable_long_0_0_1var128_erosion.asc" DrapeFileName = "SpawnedBasin_07C13-(Q8)_BASINS.bil" DEMName = "SpawnedBasin_07C13-(Q8).bil" ThisFile = DataDirectory + Filename DrapeFile = DataDirectory + DrapeFileName DEMFile = DataDirectory + DEMName Shield2 = DataDirectory + "SpawnedBasin_Q8_phi30.bil" FigFormat = 'svg' FigFileN = 'Elevation.svg' FigFileName = DataDirectory + FigFileN ShieldFigName = DataDirectory + 'Shielding.svg' NDV, xsize, ysize, GeoT, Projection, DataType = LSDP.GetGeoInfo(ThisFile) print "NDV: " + str(NDV) print "xsize: " + str(xsize) print "ysize: " + str(ysize) print "GeoT: " + str(GeoT) print "Projection: " + str(Projection) print "DataType: " + str(DataType) CellSize, XMin, XMax, YMin, YMax = LSDP.GetUTMMaxMin(ThisFile) print "CellSize: " + str(CellSize) print "XMin: " + str(XMin) print "XMax: " + str(XMax) print "YMin: " + str(YMin) print "YMax: " + str(YMax) tcmap = 'jet' drape_cmap = 'gray' shield_map = 'summer' drape_alpha = 0.4 tcmapcolorbarlabel = 'Topographic shielding' clim_val = (0.72, 1) auto_clim = (1000, 4500) #LSDP.BasicDensityPlot(ThisFile,tcmap,tcmapcolorbarlabel,clim_val) #LSDP.DrapedPlot(ThisFile,DrapeFile) #LSDP.BasicDensityPlotGridPlot(ThisFile,tcmap,tcmapcolorbarlabel) #Plot the basin over the elevation cmap_label = "Elevation (m)" LSDP.BasicDrapedPlotGridPlot(DEMFile, DrapeFile, tcmap, drape_cmap, cmap_label, auto_clim, drape_alpha, FigFileName, FigFormat) # Now plot the two shielding rasters LSDP.BasicDensityPlotGridPlot(ThisFile, shield_map, tcmapcolorbarlabel, clim_val, ShieldFigName, FigFormat) FigFormat = 'svg' SNFileN = 'Shield30.svg' ShieldFigName2 = DataDirectory + SNFileN # Now plot the two shielding rasters LSDP.BasicDensityPlotGridPlot(Shield2, shield_map, tcmapcolorbarlabel, clim_val, ShieldFigName2, FigFormat)