Bagaimana seseorang bisa menggunakan R to
- membagi shapefile dalam kotak / sub-poligon 200 meter,
- plot kisi ini (termasuk nomor ID untuk setiap kotak) di atas peta asli di bawah, dan
- mengevaluasi di mana koordinat geografis spesifik persegi berada .
Saya seorang pemula dalam GIS dan ini mungkin pertanyaan mendasar, tapi saya belum menemukan tutorial tentang cara melakukan ini di R.
Apa yang saya lakukan sejauh ini adalah memuat shapefile NYC dan merencanakan beberapa koordinat geografis yang patut dicontoh.
Saya mencari contoh (kode R) bagaimana ini dengan data di bawah ini.
# 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")