[PROJ] Help in setting proj parameters for an old map

[aborruso]

Hi,
I think my proj definition is not correct and/or complete.

If I run in example this script:

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 532.292 2433.53 405000 551000 -gcp 2618.57
2316.82 424000 551000 -gcp 2388.9 467.338 423000 567000 -gcp 537.666 569.428
407000 567000 -gcp 933.404 1364.4 410000 560000 -gcp 2090.5 1299.92 420000
560000 -gcp 1462.72 517.946 415000 567000 -gcp 1576.18 2376.68 415000 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.90 +x_0=800000 +y_0=601000 +ellps=bessel
+units=m +no_defs +k_0=0.998992911" ./output_gcp.tif ./output_warped.tif


I have a displacement of about 150 meters (as you can see in this output
https://i.imgur.com/Ulr3GE2.png)

Am I doing something wrong?

Thank you



--
Sent from: http://osgeo-org.1560.x6.nabble.com/PROJ-4-f3840930.html