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

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


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
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.