Suite

Extraction des informations de point du raster dans R


J'aimerais travailler avec les données climatiques disponibles à partir d'ici et extraire les valeurs de certains emplacements (pour lesquels j'ai lat et lon).

Voilà ce que j'ai fait jusqu'à présent :

library(SDMTools) library(raster) test <- raster(read.asc.gz("C:/… /grids_germany_annual_frost_days_1951_17.asc.gz")) > test class : RasterLayer dimensions : 866, 654, 566364 (nrow, ncol, ncell) résolution : 1000, 1000 (x, y) étendue : 3280415, 3934415, 5237501, 6103501 (xmin, xmax, ymin, ymax) coord. réf. : NA source de données : en mémoire noms : valeurs de couche : 13, 288 (min, max)

Maintenant, j'ai pensé que je cherchais peut-être leextraitfonction mais :

> extrait(test, c(3634415,5603501)) [,1] [1,] NA [2,] NA

Ajout duamortirl'argument donne :

extract(test, c(3634415,5603501), buffer=1000) Erreur dans .cellValues(x, y,… ) : argument non utilisé (buffer = 1000)

Pourriez-vous m'expliquer ce que je fais mal et comment je peux extraire des valeurs du raster ?


leouiargumenter àextraitdoit être 2D, c'est-à-dire une matrice ou un objet Spatial* (ou autre).

Tonouila requête n'est qu'un vecteur numérique atomique et est donc interprétée comme deux numéros de cellule (c'est-à-dire des indices), c'est pourquoi vous obtenez deux valeurs manquantes car elles sont hors des limites (cela déclenche un débat sur la question de savoir si ce type de hors limites devrait déclencher un autre type de "manquant").

Faites ceci :

extraire(test, matrice(c(3634415,5603501), ncol = 2))

Cela peut sembler être un comportement obscur localisé àraster, mais c'est la même chose avec d'autres fonctions R standard commeterrain. Essentiellement, votre vecteur sans dimension est traité de manière "vectorisée". Si vous donnez à plot un vecteur non dimensionnel (c'est vraiment ne pas 1-D dans le sens d'un tableau, c'est un épouvantail dans R et d'autres langages), il ajoute en fait son propre index pour l'axe X.

Je me demande ce que font l'extrait et le tracé avec de vrais tableaux 1D avec des dimensions dégénérées…


Voir la vidéo: 002 - Données Géographiques Raster u0026 Vecteur (Octobre 2021).