我有SpatialPolygonsDataFrame
的列表1,其中包含179个对象(图块)。我还有另一个包含4个对象(tile)的SpatialPolygonsDataFrame
列表2。我想基于列表2中的对象(图块)对列表1进行子集化。基本上,列表2中的四个图块也包含在179个图块中的列表1中。我尝试了下面的代码行,但它不起作用:
Subset_Forest_2000_CA<-lapply(list_shp_Forest_2000_CA,
function(x, y=list_shp_Deforested_2000_CA){subset(x, x[x@data$Tile]==y[y@data$Tile])})
有人可以帮我解决这个问题吗?对不起我无法分享数据,因为它们太大了。如果您还需要澄清,请告诉我。