[Gfoss] Help come impostare i parametri proj per una vecchia mappa

[pigreco]

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/