WMS getmapリクエストをGeotiffに変換する

フォーマット機能の1つとしてgeotiffを提供しないが、PNGおよびKMZを提供するWMSサーバーにアクセスしています。

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

私は次のいずれかになる小さなシェルスクリプトを書くことが可能であるべきだと考えています：

暗黙のSRSおよび境界ボックス情報を使用して、KMZを要求し、それをgeotiffに変換します。
PNGを要求し、WMSクエリバウンディングボックスとSRSを使用して、おそらくgdalなどを使用してgeotiffを作成します。

誰かがこのようなことをしましたか、それとももっと良い方法がありますか？

ありがとう、リッチ

gdal wms shell

— リッチシネル
ソース

これを試して：

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

タイプミスを修正しました

— user2690011
ソース

では、BBOXとCRSをWMSリクエストとgdal_translateコマンドの両方にフィードします。意味がある。そして動作します。ありがとう！

— リッチシグネル2013

lists.osgeo.org/pipermail/gdal-dev/2013-April/035951.htmlは、gdal_translateで直接実行できることを示唆しているようです-私はまだ試していません

— Ian Turton

これが、GSEが大好きな理由です。

— Brad Nesom 2013