কেউ কীভাবে আর ব্যবহার করতে পারে
- 200 মিটার স্কোয়ার / উপ-বহুভুজগুলিতে একটি শেফফিল বিভক্ত করুন ,
- নীচের মূল মানচিত্রের উপরে এই গ্রিডটি (প্রতিটি বর্গক্ষেত্রের জন্য ID নম্বর সহ) প্লট করুন এবং
- বর্গ নির্দিষ্ট নির্দিষ্ট ভৌগলিক স্থানাঙ্কে অবস্থিত যা মূল্যায়ন ।
আমি জিআইএসে একটি শিক্ষানবিস এবং এটি সম্ভবত একটি প্রাথমিক প্রশ্ন, তবে আর কীভাবে এটি করা যায় তার কোনও টিউটোরিয়াল আমি পাইনি
আমি এখন পর্যন্ত যা করেছি তা হ'ল এনওয়াইসির একটি শেফফিল লোড করা এবং কিছু অনুকরণীয় ভৌগলিক স্থানাঙ্কের প্লট করা।
আমি নীচের ডেটা সহ এটি কীভাবে একটি উদাহরণ (আর কোড) সন্ধান করছি।
# 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")