เซ็ตย่อย SpatialPolygonsDataFrame

ฉันมี SpatialPolygonsDataFrame ที่ฉันสร้างขึ้นโดยการอ่านในเชปไฟล์โดยใช้ readOGR ในแพ็คเกจ rgdal ฉันกำลังพยายามใช้มันเพื่อสร้างตารางสุ่มตัวอย่างโดยใช้ spsample ในแพ็คเกจ sp เพื่อใช้ในการประมาณค่าจากข้อมูลการสำรวจที่รวบรวมในพื้นที่ อย่างไรก็ตาม SpatialPointsDataFrame ครอบคลุมพื้นที่ขนาดใหญ่กว่าแบบสำรวจมาก และด้วยเหตุนี้ การแก้ไขจึงเป็นการทำนายค่าที่อยู่ห่างไกลจากจุดที่ดำเนินการสำรวจใดๆ ข้อมูลการสำรวจและเชปไฟล์ถูกฉายโดยใช้ proj4string เดียวกัน

ฉันต้องการสับเซ็ต SpatialPolygonsDataFrame โดยใช้พิกัดที่กำหนดโดยสถานีสำรวจ แต่ฉันไม่แน่ใจว่าค่าที่เกี่ยวข้องถูกจัดเก็บไว้ที่ใดในวัตถุ

ฉันเกรงว่าจะให้ข้อมูลที่เกี่ยวข้องไม่ได้เนื่องจากเชปไฟล์ไม่ได้โฮสต์ออนไลน์ อย่างไรก็ตาม ฉันจะยืมโค้ดบางส่วนจากคำตอบของ Paul Hiemstra ต่อ โพสต์สำหรับเนเธอร์แลนด์:

library(ggplot2)
library(sp)
library(automap)
library(rgdal)

#get the spatial data for the Netherlands
con <- url("http://gadm.org/data/rda/NLD_adm0.RData")
print(load(con))
close(con)

class(gadm)
bbox(gadm)
>         min      max
>r1  3.360782  7.29271
>r2 50.755165 53.55458

สมมติว่าการสำรวจได้ดำเนินการในพื้นที่นี้:

bbox(surveys)
    >         min      max
    >r1     4.000    7.000
    >r2    51.000   53.000

ฉันจะครอบตัดพื้นที่ของ SpatialPolygonsDataFrame นั้นได้อย่างไร

แก้ไข: คำถามนี้ ดูเหมือนจะตอบของฉัน ขออภัยที่ค้นหาไม่มากพอ (แม้ว่าความคิดเห็นจะทำให้ฉันรู้สึกได้ว่าควรหันไปใช้ rgeos ที่ไหน) อย่างไรก็ตาม gIntersection ทำให้ R ขัดข้อง...