仅将ESRI shapefile的一部分读取到R中

我有一个很大的ESRI shapefile（> 8000个多边形，.shp文件> 32MB），我想在Windows 7上读入R。

我熟悉rgdal和readOGR（）。但是，对于这种大小的shapefile，将整个shapefile读取到内存中的过程非常缓慢。此外，我只对shapefile中包含的8,000个聚集的多边形（<100个）感兴趣。我需要执行的过程是高度重复的，即shapefile需要加载多次，因此目前非常缓慢。

因此，我想知道是否有可能仅读取shapefile的一部分以加速该过程。加载如此大量的数据然后丢弃其中的大部分似乎不切实际。

唯一想到的就是利用新的“ sfr”库及其关联的简单要素类作为子集。该SFR图书馆目前可在GitHub上，这里是一个教程让你开始。

我不认为您可以基于子集查询来流数据，但是，这将大大加快读取时间，并且应该是R空间类的未来。

另一种选择是将您的空间数据存储在数据库（例如PostGIS）中并从数据库中查询。这可能是处理大数据的极其有效的方法。

我认为使用QGIS之类的GIS并选择所需的功能将是一种不错的方法。
如果需要将结果合并回去，则可以使用不想处理的7900个多边形创建一个文件，然后处理100个多边形，并在处理后重新加入。

您可以使用org2ogr命令制作shapefile的子集，如此处所示：

ogr2ogr通过属性选择特征

不过，您需要弄清楚选择所需功能的条件。