#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