জিলোতে মার্কিন যুক্তরাষ্ট্রের প্রধান শহরগুলির বিভিন্ন পাড়ার জন্য শেফফিলের একটি সেট রয়েছে । আমি পরীক্ষা করে দেখতে চেয়েছিলাম যে নির্দিষ্ট কিছু আশেপাশে কোনও বিল্ডিং উপস্থিত ছিল কিনা:
library(rgeos)
library(sp)
library(rgdal)
df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152),
Longitude = c(-122.332271,-122.353985,-122.331639),
names = c("Safeco Field", "Key Arena", "Century Link"))
coordinates(df) <- ~ Latitude + Longitude
wa.map <- readOGR("ZillowNeighborhoods-WA.shp", layer="ZillowNeighborhoods-WA")
sodo <- wa.map[wa.map$CITY == "Seattle" & wa.map$NAME == "Industrial District", ]
আমি কোনও সমস্যা ছাড়াই প্লট করতে পারি
plot(sodo)
points(df$Latitude ~ df$Longitude, col = "red", cex = 1)
আমি proj4
স্ট্রিংটি শেফফাইল থেকে আমার ডেটা ফ্রেমে মেলে
CRSobj <- CRS("+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 ")
df@proj4string <- CRSobj
over(df, sodo)
এটি কেবল আমাকে NA
মানগুলির একটি গুচ্ছ দেয় । আমি এই উত্তর চেষ্টা করেছি
spp <- SpatialPoints(df)
spp@proj4string <- CRSobj
over(spp, sodo)
তবে এখনও শুধুমাত্র NA
মান পেতে । কোন ধারণা আমার আর কি চেষ্টা করা উচিত?