如何用R中的另一个多边形shapefile剪切多边形shapefile?
问题描述:
我有两个多边形shape文件,我想剪辑一个。我在谷歌搜索,但我可以找到只有边界框裁剪或多边形裁剪点,这不是我所需要的。 除了R(http://rosettacode.org/wiki/Sutherland-Hodgman_polygon_clipping#Python)外,我还在其他编程语言中找到了一些东西。 你能帮我吗?如何用R中的另一个多边形shapefile剪切多边形shapefile?
感谢 蒂亚戈
答
削波的另一种方法是让一个选择像"polygontoclip"["templatepolygon", ]
,发现这个对点(http://robinlovelace.net/r/2014/07/29/clipping-with-r.html),而且还与多边形。
尝试使用'rgeos :: gIntersection()'作为示例(例如)[这个问题]的接受答案(http://stackoverflow.com/questions/13982773/crop-for-spatialpolygonsdataframe/13986029# 13986029)。或者,'raster :: intersect()'也应该这样做。 –
@ JoshO'Brien感谢您的回答。它帮助我获得几乎完全我想要的东西。我仍然只有一个小问题。当我做'rgeos :: gIntersection()'时,我丢失了存储在属性表中的所有行。有没有一种方法可以保持它与相交的相对应? – Tiago
也许试试'raster :: intersect()',它被明确地设计来保持关联的data.frame属性。 (而且我偶尔说“也许”的唯一原因是,它不如'rgeos :: gIntersection()'好,例如,上面链接的答案中的数据,但它通常会工作很好。)希望有所帮助! –