Glaciers thiness in the Hunza basin Rcode data_path<-'../SWE/shp/' dem <- raster(paste0(data_path, "clipped_dem.tif")) crs(dem) rgi <- st_read( file.path(data_path, "gl_wgs84.shp"))|> st_transform(crs = crs(dem)) load("glacier_thickness_extractedRGIreg13.RData") basin <- st_read(paste0(data_path, "hunza_basin.shp"), quiet = TRUE) rgi02 <- rgi |> left_join(glacier_thickness |> dplyr::select(RGIId, weighted_mean), by = "RGIId") |> rename(thickness_m = weighted_mean) |> mutate(volum_m3 = Area_m2 * thickness_m) plot(rgi02) tmap_mode("view") tm_shape(dem, name = "DEM") + tm_raster(n = 6, palette = terrain.colors(6), alpha = 0.8, legend.show = TRUE, title = "Elevation (masl)") + tm_shape(glaciers_hugonnet |> dplyr::filter(period == 20), name = "Hugonnet") + tm_fill(col = "dmdtda", n = 6, palette = "RdBu", midpoint = 0, legend.show = TRUE, title = "Glacier thinning\n(m weq/a)") + tm_shape(rgi, name = "RGI v6.0") + tm_borders(col = "light blue", lwd = 0.4) + tm_shape(basin, name = "Basin outline") + tm_borders(col = "black", lwd = 0.6)

Unsupervised Learning - K-Means & Principal component analysis (PCA)

Data 605

14082426_IrvinYeeHin.Ng_TRM32931