Convertendo solicitação getmap WMS para Geotiff

8

Estou acessando um servidor WMS que não fornece geotiff como um dos recursos de formato, mas fornece PNG e KMZ.

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

Eu estou pensando que deveria ser possível escrever um pequeno script de shell que:

  1. Solicite o KMZ e converta-o em geotiff, usando as informações implícitas do SRS e da caixa delimitadora

  2. Solicite PNG e use a caixa delimitadora de consulta WMS e o SRS para criar o geotiff, talvez usando gdal ou similar

Alguém já fez algo assim, ou existe uma maneira melhor?

Obrigado, Rich

Rich Signell
fonte

Respostas:

9

Tente o seguinte:

#!/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

Corrigido o erro de digitação

user2690011
fonte
Ok, apenas alimente o BBOX e o CRS à solicitação do WMS e ao comando gdal_translate. Isso faz sentido. E funciona. Obrigado!
Ricos Signell
2
lists.osgeo.org/pipermail/gdal-dev/2013-April/035951.html parece sugerir que você pode fazê-lo diretamente no gdal_translate - Eu não tentei embora
Ian Turton
É por isso que eu amo GSE!
Brad Nesom