Scaricare dati da una sorgente WCS con GDAL

Il WCS, Web Coverage Service, “è uno standard Open Geospatial Consortium che definisce un’interfaccia per lo scambio dei dati geospaziali sul web” (da Wikipedia).

È un servizio molto più ricco del più “noto” WMS:

Unlike OGC Web Map Service (WMS), which portrays spatial data to return static maps (rendered as pictures by the server), the Web Coverage Service provides available data together with their detailed descriptions; defines a rich syntax for requests against these data; and returns data with its original semantics (instead of pictures) which may be interpreted, extrapolated, etc., and not just portrayed.

Tramite WCS vengono distribuiti ad esempio dati relativi a modelli digitali del terreno. Questo breve post per mostrare un esempio di accesso e download di dati, utilizzando la “riga di comando” e in particolare l’utility gdal_translate della libreria GDAL.

Come dato sorgente il “MDT (Modello Digitale del Terreno) — Volo anno 2012–2013 ETRS89 — scala 1:10000” della Regione Siciliana.
Con l’URL seguente — con una chiamata GetCapabilities - si hanno indietro le caratteristiche del suddetto servizio:

http://map.sitr.regione.sicilia.it/gis/services/DEM/MDT_2012_2013/MapServer/WCSServer?request=GetCapabilities&service=WCS

Tra queste ad esempio gli spatial reference system supportati:

<SupportedCRS>urn:ogc:def:crs:EPSG::4326</SupportedCRS>
<SupportedCRS>urn:ogc:def:crs:EPSG::25833</SupportedCRS>

Per fare il download dei dati con gdal_translate, verrà usato come input un file XML con questa struttura (ServiceURL è l'URL del servizio):

<WCS_GDAL>
<ServiceURL>http://map.sitr.regione.sicilia.it/gis/services/DEM/MDT_2012_2013/MapServer/WCSServer?</ServiceURL>
<CoverageName>1</CoverageName>
</WCS_GDAL>

Il comando sarà qualcosa come:

gdal_translate -projwin_srs EPSG:3004 -projwin 2463000 4157600 2463800 4156400 -outsize 400 600 input.xml output.tif

Alcune note:

  • -projwin_srs EPSG:3004 per indicare a GDAL che le coordinate dell'area di mio interesse sono espresse con EPSG=3004 (è solo un esempio, si possono usare quelle native o altre);
  • -projwin 2463000 4157600 2463800 4156400, per indicare le coordinate del vertice in alto a sinistra (x e y) e di quello in basso a destra (x e y), dell'area di cui voglio fare il download;
  • -outsize 400 600 per definire in pixel le dimensioni del file di output. Date le sottostanti coordinate dei vertici, e data la risoluzione di 2 metri, il file di input avrà una dimensione di 400 x 600.

Volevo un’immagine per la cover del post — quella usata in alto — di 1920x762 pixel. Ho modificato così il comando (coordinate e risoluzione):

gdal_translate -projwin_srs EPSG:3004 -projwin 2463000 4157600 2466840 4156076 -outsize 1920 762 input.xml output.tif

Evviva GDAL!


Aggiornamento: Salvatore Fiandaca ha scritto, stimolato da questo articolo, il post “QGIS VS CLI : Scaricare dati da una sorgente WCS