Conversión de solicitud de getmap WMS a Geotiff

8

Estoy accediendo a un servidor WMS que no proporciona geotiff como una de las capacidades de formato, pero sí proporciona PNG y KMZ.

http://tds.marine.rutgers.edu:8081/ncWMS/wms?SERVICE=WMS&REQUEST=GetCapabilities&VERSION=1.3.0

Estoy pensando que debería ser posible escribir un pequeño script de shell que:

  1. Solicite KMZ y conviértalo a geotiff, utilizando el SRS implícito y la información del cuadro delimitador

  2. Solicite PNG y use el cuadro delimitador de consultas WMS y SRS para crear el geotiff, quizás usando gdal o similar

¿Alguien ha hecho algo como esto, o hay una mejor manera?

Gracias rico

Rich Signell
fuente

Respuestas:

9

Prueba esto:

#!/bin/bash
ulx=-76.80
lrx=-67.85
lry=34.58
uly=41.71
crs=EPSG:4326
url="http://tds.marine.rutgers.edu:8081/ncWMS/wms?LAYERS=espresso_4dvar%2Ftemp&ELEVATION=-0.986111111111111&TIME=2013-08-20T12%3A00%3A00.000Z&TRANSPARENT=true&STYLES=boxfill%2Frainbow&CRS=$crs&COLORSCALERANGE=2.156%2C28.6&NUMCOLORBANDS=254&LOGSCALE=false&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&EXCEPTIONS=XML&FORMAT=image%2Fpng&BBOX=$ulx,$lry,$lrx,$uly&WIDTH=1024&HEIGHT=1024"
curl -o foo.png $url 
gdal_translate -a_srs ${crs} -a_ullr $ulx $uly $lrx $lry foo.png foo.tif

Se corrigió el error tipográfico

usuario2690011
fuente
Bien, entonces solo alimente el BBOX y el CRS tanto a la solicitud WMS como al comando gdal_translate. Eso tiene sentido. Y funciona. ¡Gracias!
Rich Signell
2
lists.osgeo.org/pipermail/gdal-dev/2013-April/035951.html parece sugerir que puede hacerlo directamente en gdal_translate. Sin embargo, no lo he intentado
Ian Turton
¡Es por eso que amo a GSE!
Brad Nesom el