R: Bagaimana cara mendapatkan garis lintang dan bujur dari RasterLayer?

Saya seorang pemula mutlak data geografis, jadi tolong, maafkan saya jika pertanyaannya tidak tepat.

Saya mengunduh data dari NCDC NARR dan berhasil memuat ke R menggunakan rasterpaket. Saya ingin mendapatkan daftar dengan garis lintang, bujur dan nilai. Saya mengerti bahwa rasterToPoints()harus melakukan apa yang saya inginkan, namun, nilai garis lintang dan bujur saya terlihat aneh:

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

Saya kira saya harus melakukan sesuatu dengan proyeksi yang saat ini Lambert Conformal Conic (LCC). Berikut ini info lebih lanjut tentang 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

Apa yang harus saya lakukan untuk mendapatkan nilai garis lintang dan bujur AS yang sebenarnya?

Anda harus memproyeksi ulang raster ke dalam proyeksi geografis (derajat desimal) menggunakan "projectRaster" atau "spTransform". Lihat juga definisi CRS sp yang menentukan string proyeksi yang Anda inginkan. Contoh dalam bantuan untuk "projectRaster" cukup jelas dalam cara melakukan ini.

Jika Anda memaksa data raster Anda ke objek SpatialPointsDataFrame maka Anda akan menggunakan "spTransform" dan menarik koordinat dari slot @coordinates dan menambahkannya ke data.frame di slot @data. Berikut ini contohnya.

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)

Saya harus mencatat bahwa itu bukan praktik yang baik untuk mengkonversi raster ke kelas objek vektor dan meniadakan keuntungan dari paket raster yang menyediakan pemrosesan yang aman memori. Seringkali bijaksana untuk benar-benar memikirkan masalah Anda dan menilai apakah Anda mendekatinya dengan benar. Jika OP telah memberikan konteks mengapa mereka membutuhkan koordinat [x, y] untuk setiap sel, komunitas forum mungkin dapat memberikan alternatif komputasi yang akan menjaga masalah di lingkungan raster.

Dapatkan koordinat pusat sel dan buat objek Spasial:

spts <- rasterToPoints(r, spatial = TRUE)

Ubah poin ke target yang Anda inginkan:

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

Nilai-nilai sudah disalin sebagai kolom pada SpatialPointsDataFrame ini.

print(llpts)

Sekarang untuk menyelesaikan, dapatkan data.frame:

x <- as.data.frame(llpts)

Ada implementasi umum dari ini dalam paket SGAT, lihat fungsinya di lonlatFromCellsini:

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

Tampaknya Anda memiliki Koordinat yang Diproyeksikan di sana (bukan Latitude / Longitude alias Koordinat GCS). Mungkin tidak jelas bagi Anda bahwa itulah masalahnya. Lihat posting ini. Konversi sistem koordinat geografis di R