dans-gdal-scripts: uso de gdal_get_projected_bounds

El script gdal_get_projected_bounds proyecta un polígono y devuelve el rectángulo que lo limita. Para averiguar las opciones del comando, en consola, gdal_get_projected_bounds. La sintaxis que usé, en mi caso, fue:

gdal_get_projected_bounds -s_wkt archivo_wkt_imagen -s_srs Source_SRS -t_srs target_SRS -report out.ppm

Toda esa información se averigua con la ejecución de gdal_trace_outline y gdal_list_borders de la imagen fuente. El primer comando es para producir el archivo *.wkt; algo así como:

gdal_trace_outline imagen_fuente -out-cs en -wkt-out archivo_wkt

El contenido de mi archivo *.wkt fue:

POLYGON ((600915 1034685,752565 1034685,752565 879945,600915 879945,600915 1034685))

Para filtrar el código SRS en la info de la imagen fuente:

gdal_list_corners imagen_fuente| grep 's_srs'

En mi caso, éste fue el que obtuve:

'+proj=utm +zone=19 +ellps=WGS84 +datum=WGS84 +units=m +no_defs '

Después de ejecutar gdal_get_projected_bounds con esta sintaxis particular:

gdal_get_projected_bounds -s_wkt outline.wkt -s_srs '+proj=utm +zone=19 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ' -t_srs '+proj=utm +zone=19 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ' -report out.ppm

La salida fue:

bounds:
  min_e: 600914.999999978695996
  min_n: 879944.999999901512638
  max_e: 752564.999997962964699
  max_n: 1034684.999999999417923

La imagen out.ppm producida tiene este aspecto:

About these ads
Esta entrada fue publicada en GDAL/OGR, Linux, SIG, Teledetección. Guarda el enlace permanente.

Una respuesta a dans-gdal-scripts: uso de gdal_get_projected_bounds

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

Deja un comentario

Introduce tus datos o haz clic en un icono para iniciar sesión:

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