어떻게 R 을 사용 하여
- shapefile을 200m2 정사각형 / 하위 다각형으로 분할
- 아래의 원래지도 위에이 격자 (각 사각형의 ID 번호 포함)를 그리고
- 어떤 제곱 특정 지리 좌표가 있는지 평가하십시오 .
나는 GIS의 초보자이며 아마도 기본적인 질문 일지 모르지만 R에서이를 수행하는 방법에 대한 자습서를 찾지 못했습니다.
내가 지금까지 한 일은 NYC의 shapefile을로드하고 예시적인 지리적 좌표를 플로팅하는 것입니다.
아래 데이터를 사용하여 예제 (R 코드)를 찾고 있습니다.
# Load packages
library(maptools)
# Download shapefile for NYC
# OLD URL (no longer working)
# shpurl <- "http://www.nyc.gov/html/dcp/download/bytes/nybb_13a.zip"
shpurl <- "https://www1.nyc.gov/assets/planning/download/zip/data-maps/open-data/nybb_13a.zip"
tmp <- tempfile(fileext=".zip")
download.file(shpurl, destfile=tmp)
files <- unzip(tmp, exdir=getwd())
# Load & plot shapefile
shp <- readShapePoly(files[grep(".shp$", files)])
plot(shp)
# Define coordinates
points_of_interest <- data.frame(y=c(919500, 959500, 1019500, 1049500, 1029500, 989500),
x =c(130600, 150600, 180600, 198000, 248000, 218000),
id =c("A"), stringsAsFactors=F)
# Plot coordinates
points(points_of_interest$y, points_of_interest$x, pch=19, col="red")