그리드값을 폴리곤과 중첩

[Audrey]

안녕하세요. R을 이용해 분석 중 사소한 난관에 부딪혀 다시 질문올리게 되었습니다.

이번 과정의 목적은 R 에서 생성한 그리드셀값과 특정 구역.shp 파일을 중첩시켜 (중첩시키면 그리드셀값이 shp파일에 들어가겠죠......................................................? 잘모르면서 일단 해보기, 하하하;)

shp 파일을 밖으로 꺼내오려고 합니다.

R에서 그리드셀을 생성해 값을 생성한후 중첩시켰습니다.

## Perform Prediction

p <- predict(g, model = v.fit, newdata = grd)

## Intersect polygon with grid

library(raster)

district <- readShapePoly("---/DONG_123.shp")
gridded(p) = TRUE

grid.raster <- raster(p)

rr <- disaggregate(raster(grid.raster), 10)
pr <- rasterize(district, rr, -1, progress = 'text')
ra <- aggregate(pr, 10, sum, progress = 'text')

plot(ra)
plot(district, add = TRUE)

첫번째 그림이 그리드 파일이고, 두번째가 중첩한 이미지인데 제대로 중첩된거 같지 않습니다.

shp 파일 값으로 그리드 파일을 생성했기 때문에 크기가 같다고 보아야하는데 중첩한 결과를 보면... ????

그리드 파일에 shp를 중첩시켜, 최종적으로 그리드값이 포함된 shp파일을 저장할 수 있게

조언부탁드립니다.

[강성찬]

제가 이것은 해 본 적이 없어 확실히 답변드리기 어렵지만 중첩을 하려면 먼저 좌표계를 일치시켜야 할 겁니다.예를 들어 가로, 세로 범위를 좌표로 잡고 격자점 같은 걸 생성해서 해당 점에 예측값을 대입하는 식으로 해야 할 텐데 저도 자신이 없네요.

[강성찬]

http://spatial-analyst.net/book/download?size=thumbnail 여기서 pdf파일을 다운받으셔서 117페이지를 한번 참고하시면 좋을 것 같습니다.

[Audrey]

매번 도움주셔서 감사합니다. 다른 소프트웨어를 이용하면 클릭 몇번으로 쉽게 할 수 있지만 R이 계산을 더 유연하게 잘해줘서 욕심내는 중이에요. pdf 파일 참고해서 더 고민해보겠습니다^^