Conversión masiva de coordenadas UTM con cambio de datum

Para implementar este procedimiento modifiqué la aplicación que se encuentra aquí:

Conversión de coordenadas UTM con cambio de datum

para que pudiese leer de un archivo datos.in los valores a transformar, conjuntamente con todos los parámetros necesarios para la transformación, y después volcara los resultados en un archivo de salida datos.out. La estructura del archivo de entrada que usé para probar la aplicación se encuentra a continuación:

8

3 30 n 448500.79 4377580.93 771.76
4 30 n 448610.60 4377788.16 697.90
1 20 n 184598.00 1141186.00 771.76
2 20 n 184405.91 1140826.20 746.83
3 30 n 448500.79 4377580.93   0.00
4 30 n 448610.60 4377788.16   0.00
1 20 n 184598.00 1141186.00   0.00
2 20 n 184405.91 1140826.20   0.00

Al inicio tenemos el número de registros a procesar seguido de éstos con la siguiente información por línea:

    1. Código de la transformación: 1 para Canoa-Regven, 2 para Regven-Canoa, 3 para WGS 84-ED50 y 4 para ED50-WGS 84.
    2. Huso.
    3. Código para la ubicación hemisférica (Norte =n, Sur = s; pueden usarse también mayúsculas).
    4. La longitud en metros (coordenada X) y la latitud en metros (coordenada Y).
    5. La altura elipsoidal (si la desconoce puede colocar cero).

El archivo datos.in lo procesé con la aplicación que se encuentra comprimida en:

UTM_UTM_MAS0.1.zip

obteniéndose el siguiente archivo datos.out.

            X(m)        Y(m)      Z(m)   huso
     1   448610.60  4377788.16  697.90    30
     2   448500.79  4377580.93  771.76    30
     3   184405.91  1140826.20  746.83    20
     4   184598.00  1141186.00  771.76    20
     5   448610.61  4377788.17  -73.86    30
     6   448500.78  4377580.92   73.86    30
     7   184405.89  1140826.16  -24.93    20
     8   184598.02  1141186.03   24.93    20

Al omitir la altura geoidal (establecerla en cero), los resultados encontrados permiten observar que las diferencias en los valores X, Y están por debajo de 5 cm.

Esta entrada fue publicada en Código C++, SIG, Transformar Coordenadas. Guarda el enlace permanente.

9 respuestas a Conversión masiva de coordenadas UTM con cambio de datum

  1. KND #1 dijo:

    José buenisimo!!!!

  2. Rafael dijo:

    Saludos José, estoy corriendo utm_utm_masiva en una plataforma Win 7 de 64b, y cuando abre el programa y presiono una tecla, se cierra! Te ha pasado este Bug? Saludos

    • No sé si será un bug pero en XP funciona bien. Lo que si te anticipo es que fue producido con un compilador de Visual C++ que compré hace mucho tiempo (plataforma de 32 b) y que todavía conservo para su uso en XP. No he programado con Vista, ni con Seven (pero sé que el compilador tenía un comportamiento allí como el que refieres), ni tengo esperanzas de hacerlo para el 8 y más del 90 % de mi tiempo trabajo con Linux. No obstante, si estás interesado, la programación está basada en este artículo:

      http://www.gabrielortiz.com/index.asp?Info=064

      Con Python podrías programarla más fácilmente porque tiene librerías para resolver las operaciones con matrices.

      Saludos

      • Rafael dijo:

        Muchas gracias por tu pronta respuesta José! Si… yo tambien en proceso de migración a Linux con QGis y Grass. Sin embargo mira que curioso, utilicé otro programa tuyo “bursa_wolf0.2” que imagino lo hiciste con la misma versión de softw. C++ y corrió bien en la misma plataforma. Hice ajustes para intentar resolver el problema y corriera en “modo XP” (esto lo permite Win7) y hace algo distinto, pero al final, al presionar cualquier tecla, se cierra! La verdad, no estoy interesado en programar, pero gracias por el link! Saludos

  3. Entonces si puede haber alguna violación de segmento y probablemente está en el álgebra de matrices donde algún puntero debe estar mal direccionado o me olvidé de algún delete en la asignación dinámica de memoria; si es que la usé porque apenas son matrices 3×3 y con operaciones de álgebra sencilla (multiplicación y suma de matrices y multiplicación de matrices por un escalar). Para averiguarlo tendría que ubicar el código fuente y para ello tengo que buscar en varios directorios.

    Saludos

    • Rafael dijo:

      Gracias amigo! Tómate tu tiempo! Ya resolví mi pequeña duda por otro lado! Si conoces la ubicación de tutoriales de QGis (trabajo con Lisboa), en especial para redefinir proyecciones (cambios de La Canoa a WGS84) me sería de mucha ayuda. Hace poco tiempo que estoy manejando el programa y debo prepararlo para un próximo curso de postgrado en mi universidad y decidí basar la enseñanza de GIS en software de domínio público. Muchas gracias otra vez.

      • Para eso no se necesita ningún tutorial. Si es un vectorial, en Guardar como… seleccionas como SRC de la capa destino la correspondiente a WGS84 (por ejemplo, UTM 19 N el EPSG es 32619). Si es un ráster entonces sería Raster->Proyecciones->Combar(Reproyectar) seleccionando las de origen y destino. Desde Canoa/UTM 19 N, el código EPSG es 24719 donde hasta UTM WGS84 19 N (equivalente a REGVEN) el EPSG sería 32619.

        Siempre es bueno, para corroborar (calibrar), tener un punto registrado con el GPS en alguna(s) localidades(s) conocidas. Usa UTM WGS84 para minimizar los errores porque las conversiones internas de los accesibles GPSs “para navegar” producen importantes desviaciones porque no usan modelos con 7 o 10 parámetros como yo empleé en mi aplicación.

  4. no me funciona y es seguido las intruciones de este tutorial en youtube http://www.youtube.com/watch?v=lazApIojG0w

  5. todavia me encuentro en Xp 32

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s