#DataDirectory = "/mnt/SCRATCH/Analyses/HydrogeomorphPaper/peak_flood_maps/boscastle/erode_diff/"
DataDirectory = "/mnt/SCRATCH/Analyses/HydrogeomorphPaper/peak_flood_maps/boscastle/erode_diff/test_raster_diff_func/"
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)",
#trunc_cmap = LSDP.colours.truncate_colormap("Blues", 0.4, 1.0)
# Option for getting a discrete colour map
#discreet_cmap = mplext.colours.discrete_colourmap(8, "Blues")
#discreet_cmap = mplext.colours.cmap_discretize(8, trunc_cmap)
# Non-linear colourmap
levels = [-2.5, -1.5 - 0.5, -0.25, 0, 0.25, 0.5, 1.5, 2.5]
#levels = levels[levels <= tmax]
print(levels)
levels.sort()
print(levels)
cmap_lin = cm.jet
nonlincmap = LSDP.colours.nonlinear_colourmap(cmap_lin, levels)
print(type(nonlincmap))
#DataDirectory = "/run/media/dav/SHETLAND/Analyses/Ryedale_storms_simulation/Gridded/DetachLim/"
#DataDirectory = "/mnt/SCRATCH/Analyses/HydrogeomorphPaper/peak_flood_maps/boscastle/erode_diff/"
DataDirectory = "/mnt/SCRATCH/Analyses/HydrogeomorphPaper/erode_diff/test_raster_diff_func/"
#DataDirectory = "/run/media/dav/SHETLAND/Analyses/HydrogeomorphPaper/erode_diff/test_raster_diff_func/"
LSDP.MultiDrapeErodeDiffMaps(DataDirectory,
"BoscastleElevations0.asc",
"Boscastle*.bil",
cmap=cmap_lin,
drape_min_threshold=-2.5,
drape_max_threshold=2.5,
cbar_label="DEM difference (m)",
middle_mask_range=(-0.01, 0.01))
#
# return NewRaster
#DataDirectory = "/run/media/dav/SHETLAND/Analyses/Ryedale_storms_simulation/Gridded/DetachLim/"
#DataDirectory = "/mnt/SCRATCH/Analyses/HydrogeomorphPaper/peak_flood_maps/boscastle/erode_diff/"
DataDirectory = "/mnt/SCRATCH/Analyses/HydrogeomorphPaper/erode_diff/test_raster_diff_func/"
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, "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
