Obtener información de capas vectoriales con librería GDAL/OGR

Supongamos que se tiene un shapefile de una zona que nos interesa delimitar y hacer un corte en un raster al estilo ERDAS. Lo primero que se tendría que hacer es obtener información sobre el área en la cual se va a efectuar el corte antes de usar el comando para esa tarea. Ello es posible si se tiene instalada en su sistema GNU/Linux la librería GDAL/OGR.

El comando para esa operación es ogrinfo (nemotécnicamente fácil de recordar). Para más información sobre lo que hace el comando: man ogrinfo. En la salida del manual se detalla que la opción -spat produce los valores xmin, ymin, xmax y ymax del área de interés y condiciona a ogrinfo a reportar sólo información dentro del rectángulo respectivo. En nuestro caso, se va a filtrar con grep (comando de Linux) la salida de ogrinfo (con la opción -al) para ‘Extend’; que produce sólo la información requerida de xmin, ymin, xmax y ymax del área de interés.

Lo probé con un shapefile del municipio Atures (Estado Amazonas), proyectado en UTM-WGS84.

ogrinfo -al atures_plg_UTMWGS84.shp|grep ‘Extent’

cuya salida fue:

Extent: (628262.654156, 584293.869679) – (723651.956960, 685190.120639)

Esta información es la requerida por el comando gdalwarp para realizar el “corte”. De eso escribiré en otra oportunidad.

Nota: En realidad la información solicitada es más concisa. El script de bash que hace lo requerido es:

#!/bin/bash
clear;ls *.shp;echo 
echo "Cual es el archivo del que quiere la informacion?";echo
read -p "archivo = ? " nombre;echo
echo "El nombre del archivo empleado es = " $nombre;echo
echo "La informacion solicitada es: ";echo
salida1=`ogrinfo -al $nombre | grep 'Extent'`
salida2=`echo $salida1 | grep -oE '[0-9.]*'`
echo $salida2

La salida, en mi sistema, fue:

atures_plg.shp		 autana_plg.shp
atures_plg_UTMWGS84.shp  autana_plg_UTMWGS84.shp

Cual es el archivo del que quiere la informacion?

archivo = ? atures_plg_UTMWGS84.shp

El nombre del archivo empleado es =  atures_plg_UTMWGS84.shp

La informacion solicitada es: 

628262.654156 584293.869679 723651.956960 685190.120639
zeito@Debian:~/Desktop/Proyectos/Amazonas/poligonos$ 

Fue realizado con el aporte de fpuga en Espacio Linux.