Estoy ayudando a un amigo a convertir las coordenadas del diario de guerra de su difunto padre a WGS84. Tengo una definición proj.4 que parece funcionar en QGIS, pero si transformo un mapa escaneado con gdalwarp, ni se acerca. ¿Qué estoy haciendo mal con mi definición CRS, por favor?
Las coordenadas se indicaron en el Sistema británico modificado que es una variante de la cuadrícula OSGB, pero ampliada a varias zonas de Europa. Los parámetros de proyección (que me facilitó el autor del sitio enlazado) son los siguientes:
Projection - Lambert Conical Orthomorphic
Ellipsoid: Bessel 1841
False Easting : 700000
False Northing : 600000
Central Meridian : 14.0°
Central Parallel : 39.5°
Scale Factor : 0.99906
Todo, excepto el factor de escala, puede confirmarse a partir de 100k Índice de mapas topográficos de Italia de la Segunda Guerra Mundial en la Universidad McMaster. A partir de ahí, desarrollé una definición de proj.4:
+proj=lcc +lat_1=39.5 +lon_0=14 +k_0=0.99906 +x_0=700000 +y_0=600000 +ellps=bessel +units=m +no_defs
Tras convertir las coordenadas alfanuméricas en coordenadas de cuadrícula y utilizar la definición anterior, obtengo los valores esperados para algunas ubicaciones de muestra:
- rN175845 717500, 784500 41.1629 °N, 14.2086 °E
- rN138862 713800, 786200 41.1783 °N, 14.1646 °E
Sin embargo, cuando intento utilizar gdalwarp para traducir una hoja de mapa georreferenciada a WGS84, termino con latitudes en torno a 3-4 °N. La línea de comandos que he utilizado es:
gdalwarp -s_srs "+proj=lcc +lat_1=39.5 +lon_0=14 +k_0=0.99906 +x_0=700000 +y_0=600000 +ellps=bessel +units=m +no_defs" -t_srs EPSG:4326 -r cubic -of GTiff in_modified.tif out-4326.tif
gdalinfo informa de lo siguiente para el archivo de entrada:
Size is 6641, 5162
Coordinate System is `'
Origin = (694095.553313978714868,954540.625156613648869)
Pixel Size = (8.495191956798934,-8.495191956798934)
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 694095.553, 954540.625)
Lower Left ( 694095.553, 910688.444)
Upper Right ( 750512.123, 954540.625)
Lower Right ( 750512.123, 910688.444)
Center ( 722303.838, 932614.535)
Band 1 Block=6641x1 Type=Byte, ColorInterp=Red
Band 2 Block=6641x1 Type=Byte, ColorInterp=Green
Band 3 Block=6641x1 Type=Byte, ColorInterp=Blue
(las coordenadas de la muestra están en otra hoja de mapa, por cierto).
Actualización una definición de trabajo de proj.4, basada en el comentario de mkennedy más abajo, es:
+proj=lcc +lat_0=39.5 +lat_1=39.5 +lon_0=14 +k_0=0.99906 +x_0=700000 +y_0=600000 +ellps=bessel +units=m +no_defs
3 votos
Prueba también a poner lat_0 a 39,5 (y quizá lat_2). El hecho de que obtengas valores de latitud muy bajos me hace pensar que la latitud de origen se está poniendo a cero por defecto.
0 votos
Sí, eso es, gracias. ¡Lat_0 era todo lo que necesitaba!
0 votos
¿Por qué no da su actualización como respuesta?
0 votos
¿Ayudaría? Lo contestó Mkennedy, no yo.
0 votos
Sin respuesta, la pregunta aparecerá periódicamente por el
Community
usuario, y no podemos marcar una nueva pregunta como duplicada de esta porque no tiene respuesta.0 votos
Aunque ahora he añadido una respuesta, creo que la cuestión que planteas más arriba es más un problema de diseño de SE que algo que los usuarios tengan que solucionar.