R: Comment obtenir les latitudes et longitudes d'un RasterLayer?

Je suis un débutant absolu des données géographiques, alors s'il vous plaît, pardonnez-moi si la question n'est pas appropriée.

J'ai téléchargé des données depuis NCDC NARR et j'ai réussi à charger dans R à l'aide du rasterpackage. Je voudrais obtenir une liste avec latitude, longitude et valeur. Je comprends que cela rasterToPoints()devrait faire exactement ce que je veux, cependant, mes valeurs de latitude et de longitude semblent étranges:

r <- raster(myfile)
data_matrix <- rasterToPoints(r)
head(data_matrix)
            x       y value
[1,] -5405401 4347242    70
[2,] -5372938 4347242    88
[3,] -5340475 4347242    76
[4,] -5308012 4347242    85
[5,] -5275549 4347242    87
[6,] -5243086 4347242    88

Je suppose que je devrais faire quelque chose avec la projection qui est actuellement la conique conforme de Lambert (LCC). Voici d'autres informations sur le raster.

> r
class       : RasterLayer 
dimensions  : 277, 349, 96673  (nrow, ncol, ncell)
resolution  : 32463, 32463  (x, y)
extent      : -5648874, 5680713, -4628777, 4363474  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=lcc +lat_1=50 +lat_2=50 +lat_0=50 +lon_0=-107 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs 
data source : mypath-to-file
names       : value

Que dois-je faire pour obtenir des valeurs réelles de latitude et de longitude américaines?

vous devez réellement reprojeter le raster dans une projection géographique (degrés décimaux) à l'aide de "projectRaster" ou "spTransform". Regardez également les définitions de CRS sp qui spécifient la chaîne de projection souhaitée. L'exemple dans l'aide du "projectRaster" est assez clair sur la façon de procéder.

Si vous contraignez vos données raster dans un objet SpatialPointsDataFrame, vous utiliserez alors "spTransform" et tirerez les coordonnées de l'emplacement @coordinates et les ajouterez au data.frame dans l'emplacement @data. Voici un exemple de ce à quoi cela ressemblerait.

library(raster)
library(rgdal) # for spTransform

# Create data
r <- raster(ncols=100, nrows=100)
  r[] <- runif(ncell(r))
  crs(r) <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +ellps=WGS84"
  projection(r)

# Convert raster to SpatialPointsDataFrame
r.pts <- rasterToPoints(r, spatial=TRUE)
  proj4string(r.pts)

# reproject sp object
geo.prj <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0" 
r.pts <- spTransform(r.pts, CRS(geo.prj)) 
  proj4string(r.pts)

# Assign coordinates to @data slot, display first 6 rows of data.frame
r.pts@data <- data.frame(r.pts@data, long=coordinates(r.pts)[,1],
                         lat=coordinates(r.pts)[,2])                         
head(r.pts@data)

Je dois noter que ce n'est pas une bonne pratique de convertir des rasters en une classe d'objets vectoriels et que cela annule les avantages du package raster offrant un traitement sûr de la mémoire. Il est souvent prudent de vraiment réfléchir à votre problème et d'évaluer si vous l'abordez correctement. Si l'OP avait fourni un contexte expliquant pourquoi ils avaient besoin des coordonnées [x, y] pour chaque cellule, la communauté du forum aurait pu être en mesure de fournir des alternatives de calcul qui permettraient de conserver le problème dans un environnement raster.

Obtenez les coordonnées des centres cellulaires et créez un objet spatial:

spts <- rasterToPoints(r, spatial = TRUE)

Transformez les points en cible souhaitée:

library(rgdal)
llprj <-  "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0"
llpts <- spTransform(spts, CRS(llprj))

Les valeurs sont déjà copiées sous forme de colonnes sur ce SpatialPointsDataFrame.

print(llpts)

Maintenant, pour terminer, obtenez un data.frame:

x <- as.data.frame(llpts)

Il y a une implémentation générale de cela dans le package SGAT, voir la fonction lonlatFromCellici:

https://github.com/SWotherspoon/SGAT/blob/master/R/Raster.R

Il semble que vous ayez des coordonnées projetées là-bas (pas Latitude / Longitude alias coordonnées GCS). Ce n'était probablement pas clair pour vous que c'était le problème. Voir cet article. Conversion du système de coordonnées géographiques dans R