Mengonversi permintaan getmap WMS ke Geotiff


8

Saya mengakses server WMS yang tidak menyediakan geotiff sebagai salah satu kemampuan format, tetapi menyediakan PNG dan KMZ.

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

Saya berpikir bahwa mungkin untuk menulis skrip shell kecil yang akan:

  1. Minta KMZ dan ubah itu menjadi geotiff, menggunakan SRS tersirat dan info kotak pembatas

  2. Minta PNG dan gunakan kotak pembatas kueri WMS dan SRS untuk membuat geotiff, mungkin menggunakan gdal atau serupa

Adakah yang melakukan sesuatu seperti ini, atau adakah cara yang lebih baik?

Terima kasih, Rich

Jawaban:


9

Coba ini:

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

Memperbaiki kesalahan ketik


Oke, jadi beri makan BBOX dan CRS baik untuk permintaan WMS dan ke perintah gdal_translate. Masuk akal. Dan bekerja. Terima kasih!
Rich Signell

2
lists.osgeo.org/pipermail/gdal-dev/2013-April/035951.html tampaknya menyarankan Anda dapat melakukannya langsung di gdal_translate - Saya belum mencobanya
Ian Turton

Inilah sebabnya saya suka GSE!
Brad Nesom
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.