วิธีหนึ่งสามารถใช้ Rเพื่อ
- แยก Shapefileใน 200 เมตรสี่เหลี่ยม / รูปหลายเหลี่ยมย่อย
- พล็อตกริดนี้ (รวมถึงหมายเลข ID สำหรับแต่ละช่อง) เหนือแผนที่ดั้งเดิมด้านล่างและ
- ประเมินว่าพิกัดทางภูมิศาสตร์เฉพาะของรูปสี่เหลี่ยมจัตุรัสตั้งอยู่ที่ไหน
ฉันเป็นผู้เริ่มต้นใน GIS และนี่อาจเป็นคำถามพื้นฐาน แต่ฉันไม่พบบทเรียนเกี่ยวกับวิธีการทำเช่นนี้ใน R
สิ่งที่ฉันทำไปแล้วคือการโหลดรูปร่างไฟล์ของ NYC และวางแผนพิกัดทางภูมิศาสตร์ที่เป็นแบบอย่าง
ฉันกำลังมองหาตัวอย่าง (รหัส 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")