Ciao, Ho scoperto queste fantastiche mappe di dominio pubblico https://legacy.lib.utexas.edu/maps/ams/italy_50k/
Aprendo una mappa si legge [1] che la proiezione è una Lambert Conic Conforme, ed altri paramentri come lat_0, lon_0, Spheroid, x_0 e y_0 della mappa. Per una proiezione LCC [2] il parametro +lat_1 (primo parallelo standard) è richiesto ed è possibile impostare anche un secondo parallelo standard (+lat_2). In lista proj [3] ho letto che se il fattore di scala (+k) è 1 basta utilizzare solo il parametro +lat_1 che coincide con Origin (+lat_0), mentre se +k < 1 il cono è secante e quindi occorre usare entrambi i paramentri +lat_1 e +lat_2, ma per le mappe in oggetto come definire questi paramentri?? Ho fatto due test, uno con una mappa del SUD Palermo [1] ed ho usato questa stringa proj: +proj=lcc +lat_0=39.50 +lon_0=14 +lat_1=39.50 +x_0=700000 +y_0=600000 +ellps=bessel +units=m +no_defs ottenendo un ottimo risultato https://user-images.githubusercontent.com/7631137/73636590-81696d00-4666-11ea-9588-102c2b3fc6f2.png; seconda prova con una mappa del Nord, Milano: aggiungo tutto lo script utilizzato: curl -sL "http://legacy.lib.utexas.edu/maps/ams/italy_50k/txu-pclmaps-oclc-6540719-milano-45-iii.jpg";>./output.jpg gdal_translate -of GTiff -gcp 530.48 453.918 407000 568000 -gcp 2388.58 467.385 423000 567000 -gcp 2625.35 2433.37 424000 550000 -gcp 532.218 2433.37 406000 551000 ./output.jpg ./output_gcp.tif gdalwarp -r near -order 1 -co COMPRESS=PACKBITS -t_srs "+proj=lcc +lat_0=45.90 +lon_0=14 +lat_1=45.9 +x_0=800000 +y_0=601000 +ellps=bessel +units=m +no_defs +k_0=0.998992911" ./output_gcp.tif ./output_warped.tif ma osservo uno spostamente di circa 150 m https://i.imgur.com/Ulr3GE2.png cosa sbaglio?? come definire correttamente la stringa proj=lcc per mappe del sud e del nord italia?? grazie [1] https://user-images.githubusercontent.com/7631137/73595475-78b75080-4519-11ea-9c79-1968126c136d.png [2] https://proj.org/operations/projections/lcc.html [3] http://osgeo-org.1560.x6.nabble.com/PROJ-Help-in-setting-proj-parameters-for-an-old-map-td5428992.html#a5428995 ----- https://pigrecoinfinito.wordpress.com/ -- Sent from: http://gfoss-geographic-free-and-open-source-software-italian-mailing.3056002.n2.nabble.com/ _______________________________________________ Gfoss@lists.gfoss.it http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it. 764 iscritti al 23/08/2019
![http://legacy.lib.utexas.edu/maps/ams/italy_50k/txu-pclmaps-oclc-6540719-milano-45-iii.jpg"](http://legacy.lib.utexas.edu/maps/ams/italy_50k/txu-pclmaps-oclc-6540719-milano-45-iii.jpg"){kind=link}
