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.

Esta entrada fue publicada en GDAL, Linux, SIG. Guarda el enlace permanente.

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

  1. Pingback: Cómo cortar imágenes georreferenciadas con gdalwarp (GDAL/OGR) en Linux |

  2. fpugaga dijo:

    Gracias por los créditos🙂

    Yo también trabajo en GIS y la verdad es que la potencia de gdal/ogr es la leche

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s