6
วิธีการซ้อนโพลีกอนทับ SpatialPointsDataFrame และเก็บรักษาข้อมูล SPDF ไว้อย่างไร
ฉันมีSpatialPointsDataFrameข้อมูลเพิ่มเติมบางอย่างด้วย ฉันต้องการที่จะดึงจุดเหล่านั้นเข้าไปในรูปหลายเหลี่ยมและในเวลาเดียวกันให้เก็บรักษาSPDFวัตถุและข้อมูลที่เกี่ยวข้อง จนถึงตอนนี้ฉันมีโชคเล็กน้อยและใช้วิธีจับคู่และรวมกันผ่าน ID ทั่วไป แต่วิธีนี้ใช้ได้เพราะฉันมีข้อมูล gridded กับ IDS แต่ละรายการเท่านั้น นี่คือตัวอย่างย่อ ๆ ฉันกำลังมองหาจุดที่อยู่ภายในจัตุรัสแดง library(sp) set.seed(357) pts <- data.frame(x = rnorm(100), y = rnorm(100), var1 = runif(100), var2 = sample(letters, 100, replace = TRUE)) coordinates(pts) <- ~ x + y class(pts) plot(pts) axis(1); axis(2) ply <- matrix(c(-1,-1, 1,-1, 1,1, -1,1, -1,-1), …