dans-gdal-scripts: uso de gdal_list_corners

Otro de los dans-gdal-scripts es gdal_list_corners. Como su nombre sugiere, imprime información, entre otras, de las coordenadas del polígono en la cual se inscribe el raster en el lenguaje de marcado YAML. Su operación es similar a gdalinfo pero con YAML. No obstante, también produce información de las coordenadas en los puntos medios de cada lado que constituye el polígono (rectángulo) que circunscribe al raster, la proyección del raster en formato SRS (útil para ser usado en otros scripts), dimensiones de la imagen y de los píxeles (ancho x altura), etc.

A modo de ejemplo, usé gdal_list_corners para producir un script de bash que determina el área (m2/ha/km2) del raster con base en las dimensiones de la imagen y de los píxeles.

#!/bin/bash
clear;ls *.tif;echo
read -p "imagen fuente = ? " imagen_fuente
ancho=`gdal_list_corners $imagen_fuente|grep width|grep -oE [0-9]*`
alto=`gdal_list_corners $imagen_fuente|grep height|grep -oE [0-9]*`
var_tmp=`gdal_list_corners $imagen_fuente|grep res_meters|grep -oE [0-9.]*`
ancho_pixel=`echo $var_tmp|awk '{print $1}'`
alto_pixel=`echo $var_tmp|awk '{print $2}'`
area_m=$(echo "scale=10; $ancho*$alto*$ancho_pixel*$alto_pixel"|bc)
area_h=$(echo "scale=10; $area_m/10000.0"|bc)
area_k=$(echo "scale=10; $area_h/100.0"|bc)
echo "Area imagen (metros2) = " $area_m
echo "Area imagen (ha) = " $area_h
echo "Area imagen (km2) = " $area_k

Esta es la salida con una de las imágenes en el directorio en el cual lo ejecuté:

b12345script_esc.tif  b1234.tif  b1_esc.tif  b2_esc.tif  b321script_esc.tif  b321.tif  b3_esc.tif  b432.tif  b4_esc.tif  b5_esc.tif

imagen fuente = ? b321script_esc.tif
Area imagen (metros2) =  23466321000.000000000000000
Area imagen (ha) =  2346632.1000000000
Area imagen (km2) =  23466.3210000000
Esta entrada fue publicada en GDAL/OGR, SIG, Teledetección. Guarda el enlace permanente.

2 respuestas a dans-gdal-scripts: uso de gdal_list_corners

  1. Pingback: dans-gdal-scripts: uso de gdal_dem2rgb |

  2. Pingback: dans-gdal-scripts: uso de gdal_trace_outline |

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