Estoy tratando de transformar algunos datos (ráster y vectoriales) de WGS84 a un croata sistema de Coordenadas definido por:
+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs
Esto funciona muy bien en CS2CS (Rel 4.8.0):
cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to +proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs -r
Y también grandes en OGR2OGR (versión 1.9.2):
ogr2ogr -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs" out.shp in.shp
Pero en GDALWARP (versión 1.9.2), los datos sale c.400m desplazamiento hacia el Oeste (creo que es posiblemente un error de rotación basado en la aparición):
gdalwarp -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs" in.tif out.tif
Un ejemplo de esto puede verse más abajo, donde yo digitalizados una capa vectorial en torno a las actuales datos ráster, transformado tanto el uso de los comandos de arriba, y luego cargado en QGIS:
Yo no puedo por la vida de averiguar por qué esto podría ser - lo es digno de nota que he intentado hacer esto con varios datasets Ráster a partir de una variedad de fuentes, y ellos siempre parecen estar de desplazamiento hacia el Oeste por la misma cantidad.
He intentado actualizar a GDAL 1.10.1, y esto no tiene ningún efecto sobre la salida.
¿Alguien ha visto esto antes / tienes alguna idea de por qué podría ser?
Gracias!
