如何從R中的另一個shapefile獲取包含多邊形質心的shapefile的多邊形？

Question

我正在R中使用兩個shapefile，我試圖選擇其中一個包含另一個shp的質心的多邊形。

我已經能夠分別獲得每個文件的質心（附圖），但是我找不到完成上述任務的方法。在這個例子中，我們假設我只想得到裏面有藍色質心（來自shp2）的多邊形（shp1）。

example

謝謝！

Answer 1

你可以使用gCentroid()和gContains()從rgeos包：

library(raster) ## For data and functions used to make example SpatialPolygons objects 
library(rgeos) ## For topological operations on geometries 

## Make a couple of example SpatialPolygons objects, p1 & p2 
p1 <- shapefile(system.file("external/lux.shp", package="raster")) 
r <- raster(extent(p1)) 
r[] <- 1:10 
p2 <- rasterToPolygons(r, dissolve=TRUE) 

## Find centroids of p2 
cc <- gCentroid(p2, byid=TRUE) 

## Select Polygons in p1 that contain at least one of centroids from p2 
p3 <- p1[apply(gContains(p1, cc, byid=TRUE), 2, any),] 

## Plot to check that that worked 
ared <- adjustcolor("red", alpha=0.6) 
plot(p1) 
plot(p3, add=TRUE, col="wheat") 
plot(p2, add=TRUE, border=ared) 
points(cc, pch=16, col=ared)