R: Bir RasterLayer'dan enlem ve boylam nasıl alınır?

Ben coğrafi verilerin mutlak bir başlangıç, bu yüzden lütfen soru uygun değilse beni affet.

NCDC NARR'dan veri indirdim ve rasterpaketi kullanarak R'ye yüklemeyi başardım . Enlem, boylam ve değer içeren bir liste almak istiyorum. Bunun rasterToPoints()tam olarak istediğimi yapması gerektiğini anlıyorum , ancak enlem ve boylam değerlerim garip görünüyor:

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

Sanırım şu anda Lambert Konformal Konik (LCC) olan projeksiyon ile bir şeyler yapmalıyım. Raster hakkında daha fazla bilgi.

> 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

Gerçek ABD enlem ve boylam değerlerini almak için ne yapmalıyım?

gerçekten "projectRaster" veya "spTransform" kullanarak raster coğrafi (ondalık derece) bir projeksiyon yeniden projeksiyon gerekir. Ayrıca, istediğiniz projeksiyon dizesini belirten CRS sp tanımlarına da bakın. Bunun nasıl yapılacağı konusunda "projectRaster" yardımındaki örnek oldukça açıktır.

Raster verilerinizi bir SpatialPointsDataFrame nesnesine zorlarsanız, "spTransform" kullanır ve koordinatları @coordinates yuvasından alıp @data yuvasındaki data.frame öğesine eklersiniz. İşte neye benzeyeceğine bir örnek.

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)

Ben rasterleri bir vektör nesnesi sınıfına dönüştürmek iyi bir uygulama değildir ve bellek güvenli işleme sağlayan raster paketinin avantajlarını reddetmek gerekir. Sorununuzu gerçekten düşünmek ve doğru bir şekilde yaklaşıp yaklaşmadığınızı değerlendirmek genellikle ihtiyatlı olur. Eğer OP her hücre için neden [x, y] koordinatlarına ihtiyaç duyduklarına dair bir bağlam sağlamışsa, forum topluluğu sorunu raster ortamında tutacak hesaplama alternatifleri sunmuş olabilir.

Hücre merkezlerinin koordinatlarını alın ve bir Mekansal nesne oluşturun:

spts <- rasterToPoints(r, spatial = TRUE)

Noktaları istediğiniz hedefe dönüştürün:

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

Değerler zaten bu SpatialPointsDataFrame öğesinde sütunlar olarak kopyalandı.

print(llpts)

Şimdi bitirmek için bir veri alın. Çerçeve:

x <- as.data.frame(llpts)

SGAT paketinde bunun genel bir uygulaması var, buradaki işleve bakın lonlatFromCell:

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

Orada Öngörülen Koordinatlar olduğu anlaşılıyor (Enlem / Boylam, yani GCS Koordinatları değil). Sorunun bu olduğu muhtemelen açık değildi. Bu gönderiye bakın. Coğrafi koordinat sistemini R cinsinden dönüştürme

    library(raster)
r<-pt$cpc.pr[[1]] #my raster taken from a first layer of a stack
    rlon<-rlat<-r #copy r to rlon and rlat rasters [1]][1]which will contain the longitude and latitude
    xy<-xyFromCell(r,1:length(r)) #matrix of logitudes (x) and latitudes(y)
    rlon[]<-xy[,1] #raster of longitudes
    rlat[]<-xy[,2] #raster of latitides
    par(mfrow=c(1,2))
    image(rlon,main="longitudes")
    image(rlat,main="latitudes")