DTM a 5 m della LombardiaDa quest'anno (2016) il geoportale della regione Lombardia rende "disponibile" il DTM (Digital Terrain Model) con risoluzione 5 m del territorio lombardo. Si tratta di un notevole miglioramento rispetto al DTM a 20 m del Geoportale Nazionale (> catalogo WCS) o ai DEM mondiali gratiuti a 30 m. Prima d'ora per avere un DTM così dettagliato era necessario interpolare curve di livello e punti quotati. I dati altimetrici derivano da tre fonti diverse. Per la maggior parte il DTM deriva dall'interpolazione delle curve di livello del DBT (DataBase Topografico). Alcune aree, in particolare lungo le aste fluviali, sono coperte dai dati Lidar. Infine per altre zone non coperte da Lidar e dove il DBT non è ancora terminato, i dati derivano dal DTM a 20 m. DownloadCon il recente sistema di download dati del Geoportale della Regione Lombardia è ora possibile scaricare il DTM. Andando alla pagina di download del DTM 5x5 - edizione 2015 possiamo scegliere singole provincie o comunità montane come aree da scaricare, oppure filtrare in base al comune. Individuata l'area di interesse è sufficiente cliccare sulla rispettiva freccia a destra per avviare il download (necessario consentire l'apertura di finestre pop-up). Visualizzazione via WMSIl DTM è consultabile via WMS (Web Map Service) con l'indirizzo http://www.cartografia.servizirl.it/arcgis/services/wms/DTM5_RL_wms/MapServer/WMSServer La pagina dei metadati è disponibile qui: DTM 5X5 (WMS) - Modello digitale del terreno (ed. 2015) In essa viena dichiarata la presenza dei seguenti livelli: - DTM5-Base; - DTM5-Esposizione; - DTM5-Pendenza; - DTM5-Ombreggiatura; - Perimetrazione fonti. Tuttavia caricando il servizio in QGIS (> WMS in QGIS) compaioni solo due livelli: "DTM 5X5" e "Perimetrazione fonti"; il primo è una colorazione altimetrica dei dati originari e il secondo una mappa che indica le fonti dei dati (Lidar, DBT o DTM20m). DTM 5X5 + legenda: Perimetrazione fonti + legenda: Le proprietà generali sono indicate nel file xml prodotto dalla funzione "GetCapabilities": http://www.cartografia.servizirl.it/arcgis/services/wms/DTM5_RL_wms/MapServer/WmsServer?version=1.3.0&service=WMS&request=GetCapabilities LICENZA Come indicato nei metadati la licenza (che condiziona l'utilizzo legale del prodotto) è CC-BY-NC-SA 3.0 Italia, ovvero una licenza Creative Commons che consente di copiare, mostrare, ridistribuire e anche modificare i dati a patto di citarne la fonte/paternità, mantenere la stessa licenza e non usare il prodotto per scopi commerciali. (> licenze dei dati dal Geoportale Lombardia) Interrogazione WMS (GetFeatureInfo)Tramite il servizio WMS è possibile interrogare il livello "DTM 5X5" in punti specifici per ottenere sia il valore derivato (Stretchedvalue) che il valore di quota originario (PixelValue). Questa interrogazione può essere fatta in QGIS con l'apposito strumento: oppure può essere inviata direttamente come fosse un indirizzo internet usando la funzione "GetFeatureInfo" ad esempio: http://www.cartografia.servizirl.it/arcgis/services/wms/DTM5_RL_wms/MapServer/WmsServer?service=wms&version=1.3.0&request=GetFeatureInfo&crs=EPSG:32632&bbox=530000,5085000,530005,5085005&width=1&height=1&styles=&query_layers=DTM%205X5&info_format=text/xml&x=0&y=0 come si può vedere l'indirizzo è composto da diversi parametri separati da "&". I possibili valori dei parametri sono indicati nello schema XML delle "Capabilities". http://www.cartografia.servizirl.it/arcgis/services/wms/DTM5_RL_wms/MapServer/WmsServer? è l'indirizzo di base a cui si aggiungono i seguenti parametri:
nel mio esempio ho indicato il sistema di riferimento EPSG:32632 (ovvero UTM32-WGS84) e i parametri bbox=530000,5085000,530005,5085005 width=1&height=1&x=0&y=0 info_format=text/xml per interrogare il punto in alto a sinistra dell'area bbox=530000,5085000,530005,5085005 (ovvero il punto di coordinate x=530000,y=5085005) ed avere il risultato i formato xml. La risposta all'interrogazione é: <?xml version="1.0" encoding="UTF-8"?> <FeatureInfoResponse xmlns:esri_wms="http://www.esri.com/wms" xmlns="http://www.esri.com/wms"> <FIELDS Stretchedvalue="210" PixelValue="1765.035522"></FIELDS> </FeatureInfoResponse> che visualizzato in Firefox appare: <FeatureInfoResponse> <FIELDS Stretchedvalue="210" PixelValue="1765.035522"/> </FeatureInfoResponse> Interrogando in DTM al di fuori dei confini regionali, ma all'interno dei limiti (BoundingBox) indicati nelle capabilities si ottiene: <?xml version="1.0" encoding="UTF-8"?> <FeatureInfoResponse xmlns:esri_wms="http://www.esri.com/wms" xmlns="http://www.esri.com/wms"> <FIELDS PixelValue="NoData"></FIELDS> </FeatureInfoResponse> Interrogando in DTM esternamenteai limiti (BoundingBox) indicati nelle capabilities si ottiene una risposta vuota: <?xml version="1.0" encoding="UTF-8"?> <FeatureInfoResponse xmlns:esri_wms="http://www.esri.com/wms" xmlns="http://www.esri.com/wms"></FeatureInfoResponse> Utilizzando uno script è possibile interrogare automaticamente una certa area ogni 5 m per ricavare le quote e ricostruire il DTM >>> script per scaricare/ricostruire i dati DTM 5 della Lombardia Pagina visitata volte |
|