3
ตรวจสอบว่าจุดตกอยู่ในรูปหลายเหลี่ยม Shapefile
Zillow มีชุดรูปร่างที่แตกต่างกันสำหรับละแวกใกล้เคียงของเมืองสำคัญ ๆ ในสหรัฐอเมริกา ฉันต้องการตรวจสอบว่ามีอาคารบางแห่งอยู่ในละแวกใกล้เคียงโดยใช้ R: 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 …