R : RasterLayer에서 위도와 경도를 얻는 방법은 무엇입니까?

나는 지리 데이터의 절대 초보자이므로 질문이 적절하지 않은 경우 용서하십시오.

NCDC NARR에서 데이터를 다운로드하고 raster패키지를 사용하여 R에로드 할 수있었습니다 . 위도, 경도 및 가치가있는 목록을 얻고 싶습니다. 나는 rasterToPoints()내가 원하는 것을 정확하게해야한다는 것을 이해 하지만 위도 및 경도 값은 이상하게 보입니다.

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

현재 Lambert Conformal Conic (LCC) 인 투영으로 무언가를해야한다고 가정합니다. 래스터에 대한 추가 정보는 다음과 같습니다.

> 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

실제 미국 위도 및 경도 값을 얻으려면 어떻게해야합니까?

실제로 "projectRaster"또는 "spTransform"을 사용하여 래스터를 지리적 (십진도) 투영으로 다시 투영해야합니다. 또한 원하는 투영 문자열을 지정하는 CRS sp 정의를 살펴보십시오. "projectRaster"에 대한 도움말의 예제는이를 수행하는 방법에 대해 매우 분명합니다.

래스터 데이터를 SpatialPointsDataFrame 객체로 강제 변환하는 경우 "spTransform"을 사용하여 @ 좌표 슬롯에서 좌표를 가져와 @data 슬롯의 data.frame에 추가합니다. 다음은 그 모습에 대한 예입니다.

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)

래스터를 벡터 객체 클래스로 변환하고 메모리 안전 처리를 제공하는 래스터 패키지의 장점을 무효화하는 것은 좋은 방법이 아닙니다. 문제에 대해 실제로 생각하고 올바르게 접근하고 있는지 평가하는 것이 현명합니다. OP가 모든 셀에 대해 [x, y] 좌표가 필요한 이유에 대한 컨텍스트를 제공 한 경우 포럼 커뮤니티는 래스터 환경에서 문제를 유지할 계산 대안을 제공 할 수있었습니다.

셀 중심의 좌표를 가져 와서 공간 개체를 만듭니다.

spts <- rasterToPoints(r, spatial = TRUE)

포인트를 원하는 대상으로 변환하십시오.

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

이 SpatialPointsDataFrame에서 값이 이미 열로 복사되었습니다.

print(llpts)

이제 완료하려면 data.frame을 얻으십시오.

x <- as.data.frame(llpts)

SGAT 패키지에는 일반적인 구현이 lonlatFromCell있습니다. 여기의 기능을 참조 하십시오.

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

거기에 투영 좌표가있는 것으로 보입니다 (위도 / 경도가 일명 GCS 좌표가 아님). 그게 문제라는 것이 분명하지 않았을 것입니다. 이 게시물을 참조하십시오. R에서 지리 좌표계 변환