Leer sólo una parte de un shapefile de ESRI en R

Question

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

Estoy familiarizado con rgdal y readOGR(). Sin embargo, con este tamaño de shapefile, la lectura de todo el shapefile en la memoria es un proceso muy lento. Además, sólo me interesan unos pocos polígonos agrupados (<100) de los 8.000 que contiene el shapefile. El proceso para el que necesito hacer esto es altamente repetitivo, es decir, el shapefile necesita ser cargado muchas veces, por lo que actualmente es muy lento.

Por ello, me pregunto si es posible leer sólo una parte del shapefile para acelerar el proceso. No parece práctico cargar una cantidad tan grande de datos y luego descartar la mayor parte.

Answer 1

Lo único que se me ocurre es aprovechar la nueva biblioteca "sfr" y su clase de característica simple asociada para el subconjunto. El biblioteca sfr está actualmente disponible en GitHub y aquí hay un tutorial para empezar.

No creo que se puedan transmitir datos basados en una consulta de subconjunto, sin embargo, esto al menos aceleraría los tiempos de lectura considerablemente y es supuestamente el futuro de las clases espaciales de R .

Otra opción sería almacenar los datos espaciales en una base de datos (por ejemplo, PostGIS) y consultarlos desde la base de datos. Esta puede ser una forma extremadamente eficiente de tratar con datos de gran tamaño.

Answer 2

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

Answer 3

Se puede hacer un subconjunto de un shapefile utilizando la función org2ogr como se ejemplifica aquí:

ogr2ogr selección de características por atributos

Sin embargo, tendrás que establecer una condición para seleccionar las características que desees.