Konvertieren der WMS-Getmap-Anforderung in Geotiff

8

Ich greife auf einen WMS-Server zu, der Geotiff nicht als eine der Formatierungsfunktionen bereitstellt, aber PNG und KMZ bereitstellt.

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

Ich denke, dass es möglich sein sollte, ein kleines Shell-Skript zu schreiben, das entweder:

  1. Fordern Sie KMZ an und konvertieren Sie diese unter Verwendung der implizierten SRS- und Bounding-Box-Informationen in Geotiff

  2. Fordern Sie PNG an und verwenden Sie den WMS-Abfragebegrenzungsrahmen und SRS, um den Geotiff zu erstellen, möglicherweise mit gdal oder ähnlichem

Hat jemand so etwas gemacht oder gibt es einen besseren Weg?

Danke, Rich

Rich Signell
quelle

Antworten:

9

Versuche dies:

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

Tippfehler behoben

user2690011
quelle
Okay, geben Sie BBOX und CRS sowohl an die WMS-Anforderung als auch an den Befehl gdal_translate weiter. Das macht Sinn. Und funktioniert. Vielen Dank!
Rich Signell
2
lists.osgeo.org/pipermail/gdal-dev/2013-April/035951.html scheint darauf hinzudeuten , dass Sie dies direkt in gdal_translate tun können - ich habe es jedoch nicht versucht
Ian Turton
Deshalb liebe ich GSE!
Brad Nesom