공간 점 데이터 집합에서 응집력있는 공간 픽셀을 만드는 방법

Question

나는 공간 점 DF spo (불규칙한 모양의 관심 영역을 다룹니다)이 있습니다. crs 변환으로 인해 데이터가 일반 그리드에 있지 않습니다.

내 목표는 사전 정의 된 해상도와 관심 영역의 범위 (더 많은 공간 포인트 데이터가이 마스터 래스터에 매핑 됨)가있는 래스터입니다. 내가 관심의 내 영역 내에서의 NA의 물결 무늬를 얻을 실 거예요 그래서 내가

rasterize(spo, raster(ncol, nrow, extent, crs), spo$param) 

내가하는 방식으로 nrow 및 ncol을 조정해야 할 때

문제는 시작합니다. rasterize에는 보간 기능이 없으므로 사전 정의 된 (더 높은) 해상도를 사용할 수 없습니다. 이에 대한 해결책으로

, 나는 내가 관심 (단지 같은 가 library(raster); data(meuse.grid)에서을 meuse.grid)의 내 모든 지역을 커버하고, 마스터 그리드의 역할 공간 픽셀 DF spi의 어떤 종류를해야 할 수도 있습니다 생각했다. 그런 다음 데이터를 보간하는 데 사용할 수 있습니다.

idw(param~1,spo,spi) 

그리고 이에 따라 관심 분야를 원하는 해상도로 완전하게 표시하십시오. 그러나 SpatialPixelsDataFrame을 포인트 데이터에서 어떻게 생성 할 수 있습니까? 어떻게 구절 데이터 세트에서 meuse.grid을 생산하는 내보기에 그래서

는 질문에 귀결?

어쩌면 여기서 잘못된 접근 방식을 사용하고있을 수 있습니다. 다른 방법을 사용하여 내가 수행 한 작업을보다 쉽게 ​​달성 할 수 있다면 알려 주시기 바랍니다.

Answer 1

관심 영역의 경계를 정의하는 다각형이 있다면 (해야하는) 직선입니다. 한 가지 방법은 당신이 경우 ROI

이라고 불리는 관심의 영역을 정의하는 다각형이 있다고 가정하자 SpatialPoints, expand.grid 및

overlay 단지 래퍼 인 자체 geoR에서 polygrid 기능을 사용하는 것입니다 나는 현실에서

data(meuse.grid) 
coordinates(meuse.grid) = ~x+y 
x <- chull(meuse.grid@coords) 
borders <- meuse.grid@coords[c(x,x[1]),] 

ROI <- SpatialPolygons(list(Polygons(list(Polygon(borders)), ID = 'border'))) 

meuse.grid

에서 하나를 만들 것입니다 만 관심의 영역을 정의 다각형의 좌표를 필요 polygrid을 사용합니다.

당신이 당신이 재현 예를 제공하는 경우에 작동하는 것이 더 쉽습니다 polygrid

# get the bounding box for ROI an convert to a list 
bboxROI <- apply(bbox(ROI), 1, as.list) 
# create a sequence from min(x) to max(x) in each dimension 
seqs <- lapply(bboxROI, function(x) seq(x$min, x$max, by= 10)) 

# rename to xgrid and ygrid 
names(seqs) <- c('xgrid','ygrid') 

thegrid <- do.call(polygrid,c(seqs, borders = list(ROI@polygons[[1]]@Polygons[[1]]@coords)))