R: Como obter latitudes e longitudes de um RasterLayer?

Eu sou um iniciante absoluto de dados geográficos; portanto, perdoe-me se a pergunta não for apropriada.

Baixei dados do NCDC NARR e consegui carregar no R usando o rasterpacote. Gostaria de obter uma lista com latitude, longitude e valor. Entendo que rasterToPoints()deve fazer exatamente o que quero, no entanto, meus valores de latitude e longitude parecem estranhos:

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

Suponho que devo fazer algo com a projeção que atualmente é Lambert Conformal Conic (LCC). Aqui estão mais informações sobre a varredura.

> 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

O que devo fazer para obter valores reais de latitude e longitude nos EUA?

você realmente precisa reprojetar a varredura em uma projeção geográfica (graus decimais) usando "projectRaster" ou "spTransform". Observe também as definições de CRS sp que especificam a sequência de projeção desejada. O exemplo da ajuda para o "projectRaster" é bastante claro em como fazer isso.

Se você coagir seus dados raster em um objeto SpatialPointsDataFrame, use "spTransform" e puxe as coordenadas do slot @coordinates e as adicione ao data.frame no slot @data. Aqui está um exemplo de como isso seria.

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)

Devo observar que não é uma boa prática converter rasters em uma classe de objeto de vetor e nega as vantagens do pacote raster, fornecendo um processamento seguro de memória. Muitas vezes, é prudente pensar realmente no seu problema e avaliar se você o está abordando corretamente. Se o OP tivesse fornecido o contexto do motivo pelo qual eles precisam de coordenadas [x, y] para cada célula, a comunidade do fórum pode ter sido capaz de fornecer alternativas computacionais que manteriam o problema em um ambiente de varredura.

Obtenha as coordenadas dos centros celulares e crie um objeto espacial:

spts <- rasterToPoints(r, spatial = TRUE)

Transforme os pontos no seu alvo desejado:

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

Os valores já são copiados como colunas neste SpatialPointsDataFrame.

print(llpts)

Agora, para finalizar, obtenha um data.frame:

x <- as.data.frame(llpts)

Há uma implementação geral disso no pacote SGAT, veja a função lonlatFromCellaqui:

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

Parece que você tem coordenadas projetadas lá (não as coordenadas de latitude / longitude, também conhecidas como GCS). Provavelmente não estava claro para você que esse era o problema. Veja este post. Convertendo sistema de coordenadas geográficas em R