Groenlanda grosimea gheții | R-bloggeri

drop_na(valoare) |> pivot_wider(values_from = value, names_from = name) |> arrange(desc(row)) |> select(-row) |> as.matrix() |> rast(crs = „+proj=stere) +lat_0=90 +lat_ts=71 +lon_0=-39 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +type=crs”) ext(thick) = c(-800000, 700000, -3400000, -600000) Harta Iată harta brută într-o proiecție stereografică polară: gros |> plot(principal = ” Grosimea gheții din Groenlanda”, col = map.pal(„magma”)) Figura 1: Grosimea gheții din Groenlanda într-un proiecție stereografică polară. date: Bamber J., 2021. NASA National Snow and Ice Data Center Și pe o hartă interactivă după reproiectare: bibliotecă (pliant) # rezoluția nativă este de 5 km, deci la 66° N este aproximativ 0,1 grade # -75 – -15 = 65; 65 / 0,1 = 650 pixeli lățime thick_wgs84 project(rast(nrows = 250, ncols = 650, xmin = -75, xmax = -10, ymin = 60, ymax = 85, crs = „EPSG:4326”)) |> subst (x = _, 0, NA) range_m addLegend(pal = pal_rev, title = „GreenlandIce grosime (m)”, values ​​= range_m, labFormat = labelFormat(transform = function(x) sort(x + 500, descreasing = TRUE))) Figura 2: Grosimea gheții din Groenlanda Referințe Bamber, JL, RL Layberry, și SP Gogineni. 2001. „Un nou set de date despre grosimea și stratul de gheață pentru stratul de gheață din Groenlanda: 1. Măsurare, reducerea datelor și erori.” Journal of Geophysical Research: Atmospheres 106 (D24): 33773–80. https://doi.org/10.1029/2001JD900054. Bamber, Jonathan. 2001. „Greenland 5 Km DEM, Ice Thickness, and Bedrock Elevation Grids, Version 1.” Centrul național de date despre zăpadă și gheață NASA Centrul de arhivă activ distribuit. https://doi.org/10.5067/01A10Z9BM7KP. Layberry, RL și JL Bamber. 2001. „Un nou set de date despre grosimea și stratul de gheață pentru stratul de gheață din Groenlanda: 2. Relația dintre dinamică și topografia bazală.” Journal of Geophysical Research: Atmospheres 106 (D24): 33781–88. https://doi.org/10.1029/2001JD900053. ” />

(Acest articol a fost publicat pentru prima dată pe r.iresmi.netși cu amabilitate a contribuit la R-bloggeri). (Puteți raporta problema legată de conținutul acestei pagini aici)

---

Doriți să vă distribuiți conținutul pe R-bloggeri? dați clic aici dacă aveți un blog, sau aici dacă nu aveți.

Ziua 11 din 30DayMapChallenge: «Arctic» (anterior).

Vom folosi Groenlanda 5 km DEM, grosimea gheții și grile de ridicare a rocii de bază (J. Bamber 2001) din JL Bamber, Layberry și Gogineni (2001) şi Layberry și Bamber (2001). Descărcați aici (după înregistrare).

Datele necesită unele discuții, deoarece formatul nu este simplu: este un fișier ASCII cu lățime fixă ​​înfășurat (consultați ghidul utilizatorului). Trebuie să facem un rând din fiecare 31 de linii ale fișierului, să inversăm ordinea liniilor și să dăm proiecția și întinderea corecte.

library(terra)
library(readr)
library(dplyr)
library(tidyr)

thick <- read_fwf("thick_5km_corrected",
                  guess_max = 1e4) |> 
  mutate(row = ceiling(row_number() / 31)) |> 
  group_by(row) |> 
  group_modify(~ as_tibble(as.vector(t(as.matrix(.x))))) |> 
  ungroup() |> 
  mutate(name = rep(paste0("x", 1:310), 561)) |> 
  drop_na(value) |> 
  pivot_wider(values_from = value,
              names_from = name) |>
  arrange(desc(row)) |> 
  select(-row) |> 
  as.matrix() |> 
  rast(crs = "+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +type=crs") 

ext(thick) = c(-800000, 700000, -3400000, -600000)

Iată harta brută într-o proiecție stereografică polară:

thick |> 
  plot(main = "Greenland Ice thickness",
       col = map.pal("magma"))

Și pe o hartă interactivă după reproiectare:

library(leaflet)

# native resolution is 5 km, so at 66° N it's about 0.1 degree
# -75 - -15 = 65 ; 65 / 0.1 = 650 pixels wide
thick_wgs84 <- thick |>
  project(rast(nrows = 250, ncols = 650, 
               xmin = -75, xmax = -10, 
               ymin = 60, ymax = 85,
               crs = "EPSG:4326")) |> 
  subst(x = _, 0, NA)

range_m <- c(0, max(values(thick_wgs84), na.rm = TRUE))
pal <- colorNumeric("magma", domain = range_m, na.color = "#ffffff00")
# fix the legend order:
pal_rev <- colorNumeric("magma", domain = range_m, na.color = "#ffffff00", reverse = TRUE)

leaflet() |> 
  addTiles() |> 
  addRasterImage(thick_wgs84, 
                 colors = pal,
                 opacity = 0.5,
                 attribution = "Bamber, 2021 (NASA National Snow and Ice Data Center)") |> 
  addLegend(pal = pal_rev,
            title = "Greenland
Ice thickness (m)",
            values = range_m,
            labFormat = labelFormat(transform = function(x) sort(x + 500, decreasing = TRUE)))