4 votos

proj.4 definición para la cuadrícula (obsoleta) "Italia Sur" de la Oficina de Guerra del Reino Unido

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?

3voto

Bernd Haug Puntos 836

Ampliando el comentario de mkennedy, y traduciendo términos genéricos de mapas a parámetros de proj.4:

Projection - Lambert Conical Orthomorphic  →    +proj=lcc
Ellipsoid: Bessel 1841                     →    +ellps=bessel
False Easting : 700000                     →    +x_0=700000
False Northing : 600000                    →    +y_0=600000
Central Meridian : 14.0°                   →    +lon_0=14 
Central Parallel : 39.5°                   →    +lat_0=39.5 +lat_1=39.5
Scale Factor : 0.99906                     →    +k_0=0.99906
                       (other proj.4 terms)     +units=m +no_defs  

lo que da lugar a una definición de trabajo de

+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

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X