raster - extraire une ligne correspondant à un isotherm

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Laure Trudel
Messages : 24
Enregistré le : 26 Mar 2013, 11:02

raster - extraire une ligne correspondant à un isotherm

Messagepar Laure Trudel » 15 Fév 2019, 22:47

Bonjour,

Peut-être que ma réponse est déjà dans une discussion mais bien qu'ayant cherché je ne l'ai pas trouvée.

j'ai créée un rasterlayer correspond à la température de surface moyenne sur l'ensemble du globe sur plusieurs années. Maintenant je cherche à extraire un isotherm précis (21°C) de cette couche afin de couper/délimiter ensuite des trajectoires d'animaux selon cet isotherm.
Je suppose qu'il faut chercher à extraire une ligne. J'ai essayé de voir avec la fonction extract du package raster mais je n'ai pas réussi.
Avez-vous une idée?

voilà la structure du fichier - pour info la sst est en Kelvin soit je souhaiterai avoir l'isotherm de 294 K

> sst_mean
class : RasterLayer
dimensions : 720, 1440, 1036800 (nrow, ncol, ncell)
resolution : 0.25, 0.25 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source : C:\Users\user\AppData\Local\Temp\RtmpMfLA7R\raster\r_tmp_2019-02-15_205147_2736_12049.grd
names : layer
values : 271.3491, 303.3241 (min, max)

Laure Trudel
Messages : 24
Enregistré le : 26 Mar 2013, 11:02

Re: raster - extraire une ligne correspondant à un isotherm

Messagepar Laure Trudel » 28 Fév 2019, 18:45

Problème résolu !

Pour ceux que ça intéresse il faut utiliser la fonction cut

library(raster)
vol <- raster(volcano)


# bien definir l’argument breaks pour inclure ou exclure les valeurs souhaitées

vol160 <- cut(vol, breaks = c(0, 160, 200))



plot(rasterToContour(vol160, nlevels = 1))

plot(rasterToPolygons(vol160, dissolve = TRUE))


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité