R을 사용하여 위도, 경도 값 형식의 데이터를 래스터 파일로 변환하는 방법은 무엇입니까?

미국 대륙의 km 격자에 대한 데이터 값이 있습니다. 열은 "위도", "경도"및 "관측"입니다. 예 :

"lat"    "lon"     "yield"
 25.567  -120.347  3.6 
 25.832  -120.400  2.6
 26.097  -120.454  3.4
 26.363  -120.508  3.1
 26.630  -120.562  4.4

또는 R 데이터 프레임으로 :

mydata <- structure(list(lat = c(25.567, 25.832, 26.097, 26.363, 26.63), 
lon = c(-120.347, -120.4, -120.454, -120.508, -120.562), 
yield = c(3.6, 2.6, 3.4, 3.1, 4.4)), .Names = c("lat", 
"lon", "yield"), class = "data.frame", row.names = c(NA, -5L))

(전체 데이터 세트는 여기 에서 CSV로 다운로드 할 수 있습니다 )

데이터는 30km x 30km 그리드 ( Miguez et al 2012 ) 의 작물 모델 (출력 예정)에서 출력됩니다 .

지도 투영과 같은 GIS 관련 메타 데이터를 사용하여 래스터 파일로 변환하려면 어떻게해야합니까?

플랫폼과 소프트웨어에 독립적이기 때문에 파일은 텍스트 (ASCII?) 파일이 이상적입니다.

몇 가지 단계가 필요합니다.

1. 규칙적인 1km 그리드라고 말하면 위도는 규칙적이지 않습니다. 먼저 정규 그리드 좌표 시스템으로 변환하여 X 및 Y 값이 정기적으로 간격을 두어야합니다.

   에이. 열 x, y, 수율로 데이터 프레임으로 R로 읽습니다.

   pts = read.table("file.csv",......)

   비. sp 패키지를 사용하여 데이터 프레임을 SpatialPointsDataFrame으로 변환하십시오.

   library(sp)
   library(rgdal)
   coordinates(pts)=~x+y

   씨. 먼저 CRS가 무엇인지 말한 다음 일반 km 시스템으로 변환 한 다음 목적지로 spTransform하십시오.

   proj4string(pts)=CRS("+init=epsg:4326") # set it to lat-long
   pts = spTransform(pts,CRS("insert your proj4 string here"))

   디. R에게 이것이 격자로 표시되어 있다고 말합니다.

   gridded(pts) = TRUE

   이 시점에서 좌표가 규칙적인 격자에 있지 않으면 오류가 발생합니다.

2. 이제 래스터 패키지를 사용하여 래스터로 변환하고 CRS를 설정하십시오.

   r = raster(pts)
   projection(r) = CRS("insert your proj4 string here")

3. 이제 살펴보십시오 :

   plot(r)

4. 이제 래스터 패키지를 사용하여 geoTIFF 파일로 작성하십시오.

   writeRaster(r,"pts.tif")

이 geoTIFF는 모든 주요 GIS 패키지에서 읽을 수 있어야합니다. 여기서 명백히 누락 된 부분은 변환 할 proj4 문자열입니다. 이것은 아마도 일종의 UTM 참조 시스템 일 것입니다. 더 이상의 데이터없이 말하기가 어렵습니다 ...

질문에 마지막으로 답변 한 이후, rasterFromXYZ필요한 모든 단계 (CRS 문자열 지정 포함)를 캡슐화 하는 래스터 패키지 기능을 사용하면 훨씬 쉬운 솔루션이 존재합니다 .

library(raster)
rasterFromXYZ(mydata)