Transformación de CANOA (Venezuela) a WGS84 en gvSIG mediante una rejilla que simula el método Molodensky-Badekas
Puesto que gvSIG 1.9 no permite utilizar la transformación Molodensky-Badekas de 10 parámetros, una alternativa consiste en crear una rejilla o grilla de transformación en formato NTv2 utilizando los resultados de la transformación en un gran número de puntos de control repartidos por todo el territorio (en este caso, Venezuela).
Hermógenes Suárez y José Nicolás Pérez Duin han publicado en la lista de usuarios de gvSIG una hoja de cálculo con las fórmulas y los parámetros necesarios para realizar la transformación.
Para realizar la rejilla de transformación, partimos del rectángulo mínimo que contiene a todo el territorio de Venezuela:
Longitud: [-73.5, -59.75]
Latitud: [0.5, 12.25]
La distancia entre los nodos de la rejilla debe, por un lado, ser lo suficientemente pequeña como para proporcionar una precisión adecuada, pero por otro lado no debemos generar una rejilla con un número desorbitado de nodos. En el primer intento, generaremos una rejilla con una distancia entre nodos de 1/80 = 0,0125 grados (estamos hablando de longitud y latitud en el datum de partida, es decir, CANOA). Esto implica una rejilla con 1036041 nodos.
Utilizando las fórmulas de la hoja de cálculo y con la ayuda del programa awk para el procesamiento de tablas en formato ASCII, obtenemos una tabla de la forma:
LON_CANOA LAT_CANOA LON_WGS84 LAT_WGS84 ALT_WGS84
-73.500000000000 12.250000000000 -73.502055919672 12.246810195168 -12.209
-73.487500000000 12.250000000000 -73.489555608802 12.246810276383 -12.241
-73.475000000000 12.250000000000 -73.477055299826 12.246810358221 -12.273
...
-59.774999999997 0.500000000000 -59.776567401559 0.496812651416 -3.089
-59.762499999997 0.500000000000 -59.764066972853 0.496812748998 -3.111
-59.749999999997 0.500000000000 -59.751566544971 0.496812847490 -3.133
Como puede verse, por un problema de redondeos, algunas filas tienen un error de partida de 3 billonésimas de grado, lo cual es irrelevante, pues corresponde a una distancia inferior a una micra (micrómetro).
El siguiente paso consiste en obtener la diferencia entre los elementos de cada fila de la tabla, que es precisamente el cuerpo de la rejilla en formato NTv2 (ASCII). Después sólo habrá que añadir el encabezado adecuado y la rejilla estará lista. Esto se verá en el próximo artículo de esta bitácora.
Hermógenes Suárez y José Nicolás Pérez Duin han publicado en la lista de usuarios de gvSIG una hoja de cálculo con las fórmulas y los parámetros necesarios para realizar la transformación.
Para realizar la rejilla de transformación, partimos del rectángulo mínimo que contiene a todo el territorio de Venezuela:
Longitud: [-73.5, -59.75]
Latitud: [0.5, 12.25]
La distancia entre los nodos de la rejilla debe, por un lado, ser lo suficientemente pequeña como para proporcionar una precisión adecuada, pero por otro lado no debemos generar una rejilla con un número desorbitado de nodos. En el primer intento, generaremos una rejilla con una distancia entre nodos de 1/80 = 0,0125 grados (estamos hablando de longitud y latitud en el datum de partida, es decir, CANOA). Esto implica una rejilla con 1036041 nodos.
Utilizando las fórmulas de la hoja de cálculo y con la ayuda del programa awk para el procesamiento de tablas en formato ASCII, obtenemos una tabla de la forma:
LON_CANOA LAT_CANOA LON_WGS84 LAT_WGS84 ALT_WGS84
-73.500000000000 12.250000000000 -73.502055919672 12.246810195168 -12.209
-73.487500000000 12.250000000000 -73.489555608802 12.246810276383 -12.241
-73.475000000000 12.250000000000 -73.477055299826 12.246810358221 -12.273
...
-59.774999999997 0.500000000000 -59.776567401559 0.496812651416 -3.089
-59.762499999997 0.500000000000 -59.764066972853 0.496812748998 -3.111
-59.749999999997 0.500000000000 -59.751566544971 0.496812847490 -3.133
Como puede verse, por un problema de redondeos, algunas filas tienen un error de partida de 3 billonésimas de grado, lo cual es irrelevante, pues corresponde a una distancia inferior a una micra (micrómetro).
El siguiente paso consiste en obtener la diferencia entre los elementos de cada fila de la tabla, que es precisamente el cuerpo de la rejilla en formato NTv2 (ASCII). Después sólo habrá que añadir el encabezado adecuado y la rejilla estará lista. Esto se verá en el próximo artículo de esta bitácora.
En este artículo:
ReplyDeletehttp://joseguerreroa.wordpress.com/2010/09/30/conversion-de-coordenadas-utm-con-el-modelo-de-10-parametros/
programé en C++ el modelo de 10 parámetros de Molodensky-Badekas para ambos sentidos: REGVEN -> CANOA, CANOA -> REGVEN. Podrían echarle un vistazo a ver si produce los resultados esperados.
Que características tenía el punto Canoas para ser considerado Dátum en Venezuela
ReplyDelete