Conversion d'une requête getmap WMS en Geotiff

8

J'accède à un serveur WMS qui ne fournit pas de géotiff comme l'une des capacités de format, mais fournit PNG et KMZ.

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

Je pense qu'il devrait être possible d'écrire un petit script shell qui soit:

  1. Demandez KMZ et convertissez-le en géotiff, en utilisant le SRS implicite et les informations de la boîte englobante

  2. Demander PNG et utiliser la boîte englobante de requête WMS et SRS pour créer le géotiff, peut-être en utilisant gdal ou similaire

Quelqu'un at-il fait quelque chose comme ça, ou y a-t-il une meilleure façon?

Merci, Rich

Rich Signell
la source

Réponses:

9

Essaye ça:

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

Correction de la faute de frappe

user2690011
la source
D'accord, il suffit donc d'alimenter le BBOX et le CRS à la fois à la demande WMS et à la commande gdal_translate. C'est logique. Et ça marche. Merci!
Rich Signell
2
lists.osgeo.org/pipermail/gdal-dev/2013-April/035951.html semble suggérer que vous pouvez le faire directement dans gdal_translate - je n'ai pas essayé cependant
Ian Turton
C'est pourquoi j'aime GSE!
Brad Nesom