Comment peut-on utiliser R pour
- diviser un fichier de formes en carrés / sous-polygones de 200 mètres,
- tracer cette grille (y compris les numéros d’identification pour chaque carré) sur la carte originale ci-dessous, et
- Évaluer dans quel carré se trouvent les coordonnées géographiques spécifiques .
Je suis un débutant en système d’information géographique et c’est peut-être une question fondamentale, mais je n’ai pas trouvé de didacticiel expliquant comment procéder dans R.
Jusqu'à présent, j'ai chargé un fichier de formes de NYC et tracé des coordonnées géographiques exemplaires.
Je cherche un exemple (code R) comment procéder avec les données ci-dessous.
# 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")