Zillow का एक सेट है शेपफ़ाइलें प्रमुख अमेरिकी शहरों के विभिन्न इलाकों के लिए। मैं जाँच करना चाहता था कि क्या कुछ इमारतें आर का उपयोग करके कुछ पड़ोस में मौजूद थीं:
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
मूल्य प्राप्त करते हैं। किसी भी विचार मुझे और क्या प्रयास करना चाहिए?