as a series of rasters and overlays of different data, e.g. Hillshade overlaid by erosion/deposition amounts, or water depth or d50 size, etc. @author: dav """ import LSDPlottingTools as lsdplt import LSDPlottingTools.LSDMap_GDALIO as mapio # fun trivia: this is Welsh for 'mapping'. import numpy as np """Plots a hillshade overlain with the erossion or deposition amount""" folder = "/mnt/DATA/DATA/VIRTUALBOX_SHARED/HydrogeomorphPaper/BOSCASTLE/PaperSimulations/Radar1km/TransLim/" laptop_folder = "" filename = "boscastle5m_bedrock_fill.asc" drapename = "waterdepth2400.txt" # Create the drape array from one of the Catchment model output rasters drape_array = mapio.ReadRasterArrayBlocks(folder + drapename) # Optional: A lot of the output rasters contain very small values for certain # things like water depth or elevation difference, so you can mask this below: low_values_index = drape_array < 0.005 drape_array[low_values_index] = np.nan # Alternatively you can just pass the name of a secondary drape raster instead # of an array below: The function takes either a string or a ready-extracted array. lsdplt.DrapedOverHillshade(folder + filename, drape_array, clim_val=(0,400), \ drape_cmap='Blues', colorbarlabel='Elevation in meters'\ ,ShowColorbar=True, ShowDrapeColorbar=True, drape_cbarlabel = "Water depth (m)")
filename = DataDirectory + "Elevations0.asc" drapename = DataDirectory + "WaterDepths2880.asc" #LSDP.RasterDifference(DataDirectory + "Elevations0.asc", # DataDirectory + "Elevations4200.asc", # raster_band=1, # OutFileName="TestDIFF.asc", # OutFileType="AAIGrid") LSDP.MultiDrapeErodeDiffMaps(DataDirectory, "Elevations0.asc", "Ryedale*.bil", "plasma", drape_min_threshold=-4.5, drape_max_threshold=-0.05, cbar_label="DEM difference (m)") # Create the drape array from one of the Catchment model output rasters #drape_array = LSDP.ReadRasterArrayBlocks(drapename) # Optional: A lot of the output rasters contain very small values for certain # things like water depth or elevation difference, so you can mask this below: #low_values_index = drape_array < 0.005 #drape_array[low_values_index] = np.nan """ LSDP.DrapedOverHillshade(filename,drape_array,clim_val=(0,400), \ drape_cmap=trunc_cmap, colorbarlabel='Elevation in meters',\ ShowColorbar=True, ShowDrapeColorbar=True, drape_cbarlabel = "Water depth (m)", drape_alpha=1.0) """
drapename = DataDirectory + "WaterDepths2880.asc" #LSDP.RasterDifference(DataDirectory + "Elevations0.asc", # DataDirectory + "Elevations4200.asc", # raster_band=1, # OutFileName="TestDIFF.asc", # OutFileType="AAIGrid") """ LSDP.MultiDrapeErodeDiffMaps(DataDirectory, "BoscastleElevations0.asc", "BoscastleElevDiff*.bil", "seismic", drape_min_threshold= -4.5, drape_max_threshold= 4.5, middle_mask_range=(-0.02,0.02), cbar_label = "Erosion/Deposition (m)") """ # Create the drape array from one of the Catchment model output rasters #drape_array = LSDP.ReadRasterArrayBlocks(drapename) # Optional: A lot of the output rasters contain very small values for certain # things like water depth or elevation difference, so you can mask this below: #low_values_index = drape_array < 0.005 #drape_array[low_values_index] = np.nan LSDP.DrapedOverHillshade(filename,drape_array,clim_val=(0,400), \ drape_cmap=trunc_cmap, colorbarlabel='Elevation in meters',\ ShowColorbar=True, ShowDrapeColorbar=True, drape_cbarlabel = "Water depth (m)", drape_alpha=1.0)