기울기 임계 값 내에있는 도로망으로부터 거리를 멀리 두십시오.

Question

도로, 경사 및 식물로 제한되어 있다고 가정 할 때 도로를 얼마나 멀리 벗어날 수 있는지 찾는 노력을 시작했습니다. 문제를 순차적으로 다뤄 봤지만 해결할 수있는 것보다 많은 문제가 있습니다. 그래서이 게시물에 모든 것을 추가 할 것입니다. 내 첫 번째 시도는, 내가 사진을 다음 있습니다 @ 크리스토퍼-스테판에 find all the terrain that is <= 5 degrees.

덕분이었다 그래서

내 차량은, 5도 경사보다 더에서 구동 할 수 있습니다.

이전 작업이나 내가하고자하는 것을 처리하는 게시물을 찾을 수없는 것 같습니다. 그래서 설명해 드리겠습니다.

내가 차에서 길에있을 때,이 경우 팔라완, 나는 경사가 5도 (매우 중요)가되기 전에 내가 오프로드를 얼마나 멀리 운전할 수 있는지보고자한다. 정글에 빠지다 (중요하지 않음).

나는 5도 이내의 경사를 가지고 있다고 생각하고 있는데, 나는 길을 계획하고 5도 미만의 인접한 상자를 그늘지게한다. 경사도가 5 이상이면 다른 색상으로 음영 처리하십시오. (출력은이 그래픽 임).

그러면 음영 처리 된 영역 내에있는 영역의 비율을 알아야합니다 (p를 찾으면 다른 하나는 1-p입니다).

사진에 대한 코드는 당신이 http://philgis.org/country-basemaps/roads에서 Shape 파일로 도로 네트워크를 얻고 levelplot에 추가 할 수 있습니다

library(raster) 
library(rasterVis) 
elevation <- getData("alt", country = "PHL") 
x <- terrain(elevation, opt = c("slope", "aspect"), unit = "degrees") 
e <- drawExtent(show = TRUE) 
gewataSub <- crop(x, e) 
m <- c(0, 5, 0, 5, maxValue(gewataSub$slope), 1) 
rclmat <- matrix(m, ncol = 3, byrow = TRUE) 
rc <- reclassify(gewataSub$slope, rclmat) 
levelplot(rc, margin = F, col.regions = c("wheat", "gray"), colorkey = list(at = c(0, 1, 2), labels = list(at = c(0.5, 1.5), labels = c("<= 5", "> 5")))) 

감사

Answer 1

이하입니다. 거기에서 나는 당신의 예상 된 결과가 정확히 무엇인지 이해하지 못합니다.

library(raster) 
library(rasterVis) 
library(rgdal) 

elevation <- getData("alt", country = "PHL") 
x <- terrain(elevation, opt = c("slope", "aspect"), unit = "degrees") 
plot(x$slope) 

e <- drawExtent(show = TRUE) 
gewataSub <- crop(x, e) 
plot(gewataSub$slope, 1) 

m <- c(0, 5, 0, 5, maxValue(gewataSub$slope), 1) 
rclmat <- matrix(m, ncol = 3, byrow = TRUE) 
rc <- reclassify(gewataSub$slope, rclmat) 

roads = readOGR("/Users/christopherstephan/Downloads", layer = "roads") 

levelplot(
    rc, 
    margin = F, 
    col.regions = c("wheat", "gray"), 
    colorkey = list(at = c(0, 1, 2), labels = list(
    at = c(0.5, 1.5), labels = c("<= 5", "> 5") 
)) 
) + layer(sp.polygons(roads))