Problemi con la cartina fisica dell'italia dal portale cartografico nazionale

30 views
Skip to first unread message

Cristiano Longo

unread,
Aug 16, 2022, 4:16:57 AM8/16/22
to opendatasicilia
Salve a tutti, come ben sapete sono un po "niubbo" nel settore GIS.
Volevo scaricare la cartina fisica dell'Italia dal portale cartografico
nazionale.

In http://www.pcn.minambiente.it/mattm/servizio-wms/ è disponibile il
catalogo dei servizi wms. Poniamo "Cartografia di base - IGM 25.000".
Suppongo che 25000 sia la scala. Per ottenere le capabilities eseguo la
seguente chiamata.

http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/IGM_25000.map&service=wms&request=getCapabilities&version=1.3.0

Tra le capabilities è presente Request > GetMap.

Tra le capabilities inoltre sono elencati vari layers. Scelgo il layer
con nome CB.IGM25000.32 (il primo che trovo per il quale è definito
l'attributo name)

Effettuo la richiesta getMaps secondo le capability indicate: uso il crs
EPSG:4326 e il corrispondente bounding box (sempre elencato nelle
capabilities) 35.0343,6.33499,47.3109,19.8408, come formato image/png e
configuro altezza e larghezza in modo che siano proporzionate rispetto
al bounding box.

http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/IGM_25000.map&service=wms&request=getMap&version=1.3.0&layers=CB.IGM25000.32&crs=EPSG%3A4326&bbox=35.0343,6.33499,47.3109,19.8408&format=image%2Fpng&width=1228&height=1350


Il risultato è un rettangolo bianco. Questo vuol dire che qualcosa non
ha funzionato. Suggerimenti?

Grazie in anticipo

andy

unread,
Aug 17, 2022, 7:47:25 AM8/17/22
to opendat...@googlegroups.com

Caro Cris,che piacere incontrarti.

On Tue, 16 Aug 2022 at 10:16, Cristiano Longo cristia...@opendatahacklab.org wrote:

In http://www.pcn.minambiente.it/mattm/servizio-wms/ è disponibile il
catalogo dei servizi wms. Poniamo "Cartografia di base - IGM 25.000".
Suppongo che 25000 sia la scala. Per ottenere le capabilities eseguo la
seguente chiamata.

http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/IGM_25000.map&service=wms&request=getCapabilities&version=1.3.0

Tra le capabilities è presente Request > GetMap.

Tra le capabilities inoltre sono elencati vari layers. Scelgo il layer
con nome CB.IGM25000.32 (il primo che trovo per il quale è definito
l'attributo name)

Effettuo la richiesta getMaps secondo le capability indicate: uso il crs
EPSG:4326 e il corrispondente bounding box (sempre elencato nelle
capabilities) 35.0343,6.33499,47.3109,19.8408, come formato image/png e
configuro altezza e larghezza in modo che siano proporzionate rispetto
al bounding box.

http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/IGM_25000.map&service=wms&request=getMap&version=1.3.0&layers=CB.IGM25000.32&crs=EPSG%3A4326&bbox=35.0343,6.33499,47.3109,19.8408&format=image%2Fpng&width=1228&height=1350

Alcuni punti spero utili. Tra le cose che hai in GetCapabilities, hai

<Abstract>Carta topografica d'Italia alla scala 1:25.000. Risoluzione 1:25.000. Il servizio è visualizzabile solo a scale superiori a 1:100.000.</Abstract>

Dai numeri che vedo, hai scelto un’areale molto grande, da rappresentare in pochi pixel. Quindi siamo con una scala probabilmente molto più piccola di 1:100.000 e quindi non vedi nulla.
Il GetCapabilities te lo dà anche in forma machine readable

<MaxScaleDenominator>100000</MaxScaleDenominator> </Layer>

Poi occhio ai due layer presenti, che hanno a che fare il con il fuso 32 (che va da 6° di longitudine a 12°) e il fuso 33 (che va da 12° di longitudine a 18°).

Nella tua chiamata inoltre mi sembra che non rispetti le specifiche di bbox, che devono essere “minx,miny,maxx,maxy in units of the SRS”, mentre le tue sembrano miny,minx,maxy,max.

Qui un esempio di chiamata, nel fuso 33, usando un sistema di coordinate metrico (è più comodo per i calcoli di scala):

http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/IGM_25000.map&REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&LAYERS=CB.IGM25000.33&STYLES=&FORMAT=image/png&CRS=EPSG:32633&BBOX=311849,4212710,313849,4214710&WIDTH=2000&HEIGHT=2000

Devi scaricare tutta l’Italia? Per stamparla su carta? Per lavorarci su un PC?
Se mi dettagli cosa devi fare, posso provare a suggerirti una modalità di lavoro.

Buon pranzo

--
___________________

Andrea Borruso
website: https://medium.com/tantotanto
38° 7' 48" N, 13° 21' 9" E, EPSG:4326
___________________

"cercare e saper riconoscere chi e cosa,
 in mezzo all’inferno, non è inferno, 
e farlo durare, e dargli spazio"

Italo Calvino

Cristiano Longo

unread,
Aug 30, 2022, 3:54:36 AM8/30/22
to opendat...@googlegroups.com

Grazie Andrea ti rispondo in linea.

ma in realtà ho copiato da quanto riportato nelle capabilities.

<BoundingBox CRS="EPSG:4326" maxx="47.3109" maxy="19.8408" minx="35.0343" miny="6.33499"/>

Grazie, ancora non riesco a raccapezzarmi del perchè la tua chiamata funziona e la mia no. Ho provato anche con un bounding box più piccolo (sempre sistema di coordinate wsg84) ma niente.

Devi scaricare tutta l’Italia? Per stamparla su carta? Per lavorarci su un PC?

Volevo scaricarla sul PC per usarla come base per dei disegnini.

Se mi dettagli cosa devi fare, posso provare a suggerirti una modalità di lavoro.

Buon pranzo

--
___________________

Andrea Borruso
website: https://medium.com/tantotanto
38° 7' 48" N, 13° 21' 9" E, EPSG:4326
___________________

"cercare e saper riconoscere chi e cosa,
 in mezzo all’inferno, non è inferno, 
e farlo durare, e dargli spazio"

Italo Calvino
--
Sito: http://opendatasicilia.it
Facebook: https://www.facebook.com/groups/opendatasicilia/
twitter: http://twitter.com/opendatasicilia
Gruppo Telegram: https://t.me/opendatasicilia
---
Hai ricevuto questo messaggio perché sei iscritto al gruppo "opendatasicilia" di Google Gruppi.
Per annullare l'iscrizione a questo gruppo e non ricevere più le sue email, invia un'email a opendatasicil...@googlegroups.com.
Per visualizzare questa discussione sul Web, visita https://groups.google.com/d/msgid/opendatasicilia/CAHEdGZNTnxajj9g1XS0yCwB-JusyAs4HHroztyDZdypfxzgnjQ%40mail.gmail.com.

andy

unread,
Aug 30, 2022, 4:42:31 AM8/30/22
to opendat...@googlegroups.com
Ciao Cris,


On Tue, 30 Aug 2022 at 09:54, Cristiano Longo <cristia...@opendatahacklab.org> wrote:

Nella tua chiamata inoltre mi sembra che non rispetti le specifiche di bbox, che devono essere “minx,miny,maxx,maxy in units of the SRS”, mentre le tue sembrano miny,minx,maxy,max.

ma in realtà ho copiato da quanto riportato nelle capabilities.

<BoundingBox CRS="EPSG:4326" maxx="47.3109" maxy="19.8408" minx="35.0343" miny="6.33499"/>

e ti confermo che c'è qualcosa che mi torna.
Perché l'italia ricade (non esattamente) circa tra 12 e 18 gradi di longitudine, tra 12 e 18 di x. Mentre qui leggo un minx di "35.0343".
 

Devi scaricare tutta l’Italia? Per stamparla su carta? Per lavorarci su un PC?

Volevo scaricarla sul PC per usarla come base per dei disegnini.

Tutta l'Italia a quella scala, al 25.000 è un file di discrete dimensioni.
Quindi zommerai a Brescia e disegnerai qualcosa, poi a Canicattì, un po' in tutto il territorio?

Con che strumento farai i disegnini? 
Perché un programma di disegno potrebbe avere difficoltà o barriere grosse al disegno su questa enormità di pixel da scaricare.
Mentre un file ottimizzato a tasselli, e piramidato (un copia della mappa "renderizzata" a diverse scale) si gestisce con grande agitiltà.
Ma piramide e tasselli sono gestiti da software GIS.

Una soluzione easy, potrebbe essere proprio quella di usare soltanto i servizi WMS. Ti installi QGIS, carichi il layer come servizio, zoomi nel punto di tuo interesse e fai i disegnini.
Non sei costretto a scaricare nulla.

Ripeto, se dettagli meglio il vero caso d'uso, oltre "disegnini", è più facile suggerirti una soluzione giusta.

Ti allego un file QGIS 3, con "acceso" il layer WMS.

Saluti
cris.qgz
Reply all
Reply to author
Forward
0 new messages