Leer solo parte de un archivo de forma ESRI en R

Tengo un archivo de forma ESRI muy grande (> 8000 polígonos, el archivo .shp es> 32 MB), que quiero leer en R en Windows 7.

Estoy familiarizado con rgdal y readOGR (). Sin embargo, con este tamaño de shapefile, leer todo el shapefile en la memoria es un proceso muy lento. Además, solo estoy interesado en algunos polígonos agrupados (<100) de los 8,000 que están contenidos en el archivo de forma. El proceso para el que necesito hacer esto es muy repetitivo, es decir, el archivo de forma debe cargarse muchas veces, por lo que actualmente es muy lento.

Por lo tanto, me pregunto si es posible leer solo una parte del archivo de forma para acelerar el proceso. Parece poco práctico cargar una cantidad tan enorme de datos y luego descartar la mayor parte.

Lo único que viene a la mente es aprovechar la nueva biblioteca "sfr" y su clase de entidad simple asociada para el subconjunto. La biblioteca sfr está actualmente disponible en GitHub y aquí hay un tutorial para comenzar.

No creo que se puede transmitir datos en base a una consulta subconjunto sin embargo, esto haría a la velocidad de arrendamiento hasta los tiempos de lectura considerablemente y es supuestamente el futuro de las clases espaciales R .

Otra opción sería almacenar sus datos espaciales en una base de datos (por ejemplo, PostGIS) y consultar desde la base de datos. Esta puede ser una forma extremadamente eficiente de manejar grandes datos.

Creo que usar un SIG como QGIS y seleccionar las características que necesita sería un buen enfoque.
Si necesita fusionar los resultados, puede crear un archivo con los 7900 polígonos que no desea procesar, luego procesar los 100 polígonos y volver a unirse después del procesamiento.

Puede hacer un subconjunto de un archivo de forma usando el org2ogrcomando, como se ejemplifica aquí:

ogr2ogr seleccionando características por atributos

Sin embargo, necesitarás determinar una condición para seleccionar las características que deseas.