2017-12-11 24 views
0

위도 (y)와 경도 (y)로 태그가 지정된 데이터 세트가 있습니다.data.frame을 SpatialPixelsDataFrame으로 변환하는 방법 R

x   y   pH  EC  RSC SAR 
644602.7184 3295635.709 7.82 1.069 6.8 2.612316328 
644219.281 3293883.389 7.45 2.27 6 2.40812669 
644368.6088 3295296.452 7.41 1.17 7.2 1.240347346 
646050.5728 3297548.156 7.8  0.971 4.5 1.589470393 
644641.7003 3295115.255 7.87 0.708 4.4 0.875069715 
644764.8507 3294939.506 8  0.634 4.1 2.154375393 
644885.5937 3295445.407 7.83 1.13 7.1 0.951290419 

SpatialPixelsDataFrame으로 변환하고 싶습니다. 다음 코드를 사용하고 있습니다 :

data <- read.csv("Try_data.csv", header = T) 
coordinates(data) <- ~x+y 
proj4string(data) <- CRS("+proj=utm +zone=43 ellps=WGS84") 

grid = spsample(data, type = "regular", cellsize = c(29.11785674, 29.11785674)) 
gridded(grid) = TRUE 
grid2 <- as(grid, "SpatialPixelsDataFrame") 

그러나 grid2 내에서 data.frame에는 0obs가 있습니다. 0 개의 변수로 구성됩니다. pH, EC, RSC, SAR 파라미터가 grid2 내의 data.frame에 있어야합니다. 그것을하는 방법?

+0

이 질문은 다음에서 또한 요청되었습니다 : https://gis.stackexchange.com/questions/264892/how-to-convert-data-frame-to-spatialpixelsdataframe- in-r – RobertH

답변

0

예제 데이터를 스크립트로 제공해주십시오. 불규칙한 점 데이터로부터 정규 래스터를 만들고 싶습니다. raster::rasterize과 같은 함수를 사용하거나 보간법을 사용할 수 있습니다. 예 : raster::interpolate 및 gstat와 같은 패키지의 메소드

+0

보간을 사용하여 data.frame을 SpatialPixelsDataFrame으로 변환했습니다. –