con due laptop e vari software GIS
sistema:
- processore: AMD FX-7500 Radeon R7, 10 Compute Cores 4C + 6G 2,10 GHz
- RAM: 8,00 GB
- Tipo sistema: Sistema operativo a 64 bit, processore basato du x64
- storage: HDD
sistema:
- processore: Intel i7-7700 2,10 GHz
- RAM: 16,00 GB
- Tipo sistema: Sistema operativo a 64 bit, processore basato du x64
- storage: SSD
Confini delle unità amministrative a fini statistici al 1 gennaio 2018 -ISTAT - (anno 2018 - versione generalizzata) cartella Com01012018_g. NB: gli shapefile dell'ISTAT hanno geometrie non valide, la correzione è stata fatta tramite spatialite_gui - leggere blog post su pigrecoinfinito
Comuni ISTAT | nro vertici | senza duplicati | nro duplicati |
---|---|---|---|
7960 MultiPolygon | 853.700 | 432.355 | 421.345 |
--> torna su
- QGIS 2.18.24, 3.2.3, 3.3 master (05/10/2018)
- R 3.5.1 + RStudio 1.1.456 (nella console: install.packages("rgdal") e poi install.packages("GISTools") )
- SpatiaLite_gui 2.10
- PostgreSQL 9.3/PostGIS 2.2.3
- Mapshaper 0.4.94 (riga di comando - GNU/Linux - Ubuntu 16.4 in Win10 64b)
## per installarlo
sudo apt-get update
sudo apt-get install nodejs
sudo apt-get install npm
sudo npm install -g mapshaper
## per aggiornarlo:
sudo npm update -g mapshaper
- estrai vertici;
- elimina geometrie duplicate (alludo ai vertici);
- dissolvi per regione (campo "cod_reg");
- buffer -2 km (su dissolvi per regione);
- conta punti nel poligono;
- spatial join (tra vertici e regione - trasferire campo "cod_reg");
Descrizione: il dataset rappresenta un insieme di (#7.960) poligoni (Comuni italiani ISTAT generalizzati) da cui estrarre i vertici (#853.700): questi, in corrispondenza dei limiti adiacenti tra i Comuni, saranno duplicati e quindi con la seconda prova si eliminano queste geometrie (#432.355); la terza prova richiede di dissolvere, il dataset di partenza, usando il campo "cod_reg" ottenendo i limiti regionali ISTAT generalizzati (#20); la quarta prova richiede un semplice buffer, negativo, di 2 km sul dataset risultante dal dissolve; la quinta prova conteggia il numero di punti per ogni regione con buffer -2 km; infine, ultima prova, trasferire il campo "cod_reg" nel vettore vertici senza duplicati. (Il dataset e le prove sono state studiate per essere facilmente realizzate nel mio laptop LZ50)
--> torna su
È importante (per PC/Laptop datati, come il mio) chiudere tutti gli applicativi non interessati alla prova in modo che ogni risorsa sia disponibile per i test; inoltre, la prova va ripetuta due volte e il tempo impiegato espresso in secondi. L'output va sempre salvato in un file (shp, csv) o tabella (in caso di DBMS).
Per i DBMS SpatiaLite e PostGIS - creazione indici spaziali
--> torna su
Grazie a questi test è stato possibile 'aprire'
quattro ticket (segnalazioni) sul bugtracker di QGIS, le quattro segnalazioni sono state risolte e la prossima QGIS 3.4 sarà più performante (Grazie Nyall e Alessandro).
--> torna su
- Ludovico Frate per l'idea e gli script per R + RStudio;
- Andrea Borruso per avermi fatto conoscere mapshaper e consigliato i relativi script