Molodensky-Badekas para Atacama y Antofagasta

Con base en la referencia:

Parámetros de Transformación entre Sistemas de Referencia Geodésicos y Cartográficos

se extendió el programa que se encuentra aquí:

Conversión de coordenadas UTM con el modelo de 10 parámetros

para incluir transformaciones de datum PSAD56-SIRGAS y SIRGAS-PSAD56, usando el modelo de Molodesnsky-Badekas, para las regiones de Atacama y Antofagasta, Chile. Los parámetros usados se resumen a continuación:

	PSAD56 a SIRGAS	-> elipsoide de entrada = Hayford; elipsoide de salida = WGS84	

	Área R1-2				
	-10.71	-10.64	15.22	1.000009093	
	2065900	-5510100	-2441700		
	-305.257	278.485	-358.598		
					
	Área R1				
	-13.87	-11.73	19.68	1.000005214	
	2103800	-5615900	-2165300		
	-298.216	259.531	-363.822		
					
	Área R2				
	-10.15	-9.69	13.09	1.000010216	
	2063000	-5477800	-2524600		
	-306.954	284.475	-356.698		
					
	Área R2N				
	-13.40	-10.51	13.89	1.000008413	
	2097600	-5522000	-2401400		
	-303.299	275.337	-359.482		

La información anterior fue incluida en un programa escrito en C/C++, donde el copy/paste de la definición de parámetros dentro del programa principal se incluye a continuación con propósitos de verificación:

	if (selector == 3) {

//Molodensky-Badekas para PSAD56-SIRGAS-R1-2-cl 

		 Rx= -10.71; Ry= -10.64; Rz= 15.22; esc = 1+(0.000009093);

 		 Xp = 2065900; Yp = -5510100; Zp = -2441700;

		 D[0][0] = -305.257; D[1][0] = 278.485; D[2][0] = -358.598;

	// Elipsoide de WGS 84

		cout << "\nElipsoide de salida WGS 84\n";

		a = 6378137.0;  b = 6356752.3142;

//*************************************************************

	}

	if (selector == 4) {

//Molodensky-Badekas para SIRGAS-PSAD56-R1-2-cl 

		 Rx = 10.71; Ry= 10.64; Rz= -15.22; esc = 1+(-0.000009093);

 		 Xp = 2065900; Yp = -5510100; Zp = -2441700;

		 D[0][0] = 305.257; D[1][0] = -278.485; D[2][0] = 358.598;

	// Elipsoide de Hayford

		cout << "\nElipsoide de salida Hayford\n";

		a = 6378388.0, b = 6356911.946130;

//*************************************************************

	}

	if (selector == 5) {

//Molodensky-Badekas para PSAD56-SIRGAS-R1-cl

		 Rx= -13.87; Ry= -11.73; Rz= 19.68; esc = 1+(0.000005214);

 		 Xp = 2103800; Yp = -5615900; Zp = -2165300;

		 D[0][0] = -298.216; D[1][0] = 259.531; D[2][0] = -363.822;

	// Elipsoide de WGS 84

		cout << "\nElipsoide de salida WGS 84\n";

		a = 6378137.0;  b = 6356752.3142;

//*************************************************************

	}

	if (selector == 6) {

//Molodensky-Badekas para SIRGAS-PSAD56-R1-cl 

		 Rx= 13.87; Ry= 11.73; Rz= -19.68; esc = 1+(-0.000005214);

 		 Xp = 2103800; Yp = -5615900; Zp = -2165300;

		 D[0][0] = 298.216; D[1][0] = -259.531; D[2][0] = 363.822;

	// Elipsoide de Hayford

		cout << "\nElipsoide de salida Hayford\n";

		a = 6378388.0, b = 6356911.946130;

//*************************************************************

	}

	if (selector == 7) {

//Molodensky-Badekas para PSAD56-SIRGAS-R2-cl 

		 Rx= -10.15; Ry= -9.69; Rz= 13.09; esc = 1+(0.000010216);

 		 Xp = 2063000; Yp = -5477800; Zp = -2524600;

		 D[0][0] = -306.954; D[1][0] = 284.475; D[2][0] = -356.698;

	// Elipsoide de WGS 84

		cout << "\nElipsoide de salida WGS 84\n";

		a = 6378137.0;  b = 6356752.3142;

//*************************************************************

	}

	if (selector == 8 ) {

//Molodensky-Badekas para SIRGAS-PSAD56-R2-cl 

		 Rx= 10.15; Ry= 9.69; Rz= -13.09; esc = 1+(-0.000010216);

 		 Xp = 2063000; Yp = -5477800; Zp = -2524600;

		 D[0][0] = 306.954; D[1][0] = -284.475; D[2][0] = 356.698;

	// Elipsoide de Hayford

		cout << "\nElipsoide de salida Hayford\n";

		a = 6378388.0, b = 6356911.946130;

//*************************************************************

	}

	if (selector == 9) {

//Molodensky-Badekas para PSAD56-SIRGAS-R2N-cl 

		 Rx= -13.40; Ry= -10.51; Rz= 13.89; esc = 1+(0.000008413);

 		 Xp = 2097600; Yp = -5522000; Zp = -2401400;

		 D[0][0] = -303.299; D[1][0] = 275.337; D[2][0] = -359.482;

	// Elipsoide de WGS 84

		cout << "\nElipsoide de salida WGS 84\n";

		a = 6378137.0;  b = 6356752.3142;

//*************************************************************

	}

	if (selector == 10) {

//Molodensky-Badekas para SIRGAS-PSAD56-R2N-cl

		 Rx= 13.40; Ry= 10.51; Rz= -13.89; esc = 1+(-0.000008413);

 		 Xp = 2097600; Yp = -5522000; Zp = -2401400;

		 D[0][0] = 303.299; D[1][0] = -275.337; D[2][0] = 359.482;

	// Elipsoide de Hayford

		cout << "\nElipsoide de salida Hayford\n";

		a = 6378388.0, b = 6356911.946130;

//*************************************************************

	}

El ejecutable para Windows, en formato comprimido, puede ser bajado del siguiente link:

badekas0.3.zip

Estoy en procura de la información correspondiente a los vértices geodésicos para las regiones consideradas a fin de proceder a la validación respectiva.

Links relacionados:

Validación de transformación de datum (Molodensky-Badekas) para Antofagasta

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

6 respuestas a Molodensky-Badekas para Atacama y Antofagasta

  1. Pingback: Validación de transformación de datum (Molodensky-Badekas) para Antofagasta |

  2. angel garcia dijo:

    Hola, es estado probando tu convertidor de coordenadas badekaso.3 y tengo algunas dudas. Si me puedes ayudar con un ejemplo directo para santiago de chile. Por ejemplo yo tomo las coordenadas del google eart que, segun he averiguado estan en sirgas- wgs84, o regven. Yo quiero cambiarlas a coordenadas UTM psad 56.
    En el programa existen diferentes parametros que no se para que sirven por ejemplo el 4, 6,8 y 10 que son para cambiar de 84 a 56 pero nose su finalidad. Me gustaria saber si para santiago valparaiso y coquimbo cual es el que debo seleccionar, me podrias dar un ejemplo de como realizar la operación para pasar de coordemnadas del google en santiago que son 84 a utm 56. Muchas gracias por tu ayuda.

    • Regreso el 7 de Enero, Dios mediante. Después te contesto mejor. En principio, necesitas otro juego de parámetros diferente. Estos son válidos para el norte de Chile. Para Santiago, si consigo los parámetros sólo tendría que hacer una ligera modificación al programa para hacer el cálculo. A veces hay que googlear algo para conseguirlos y luego para validarlos con puntos que estén en ambas sistemas de coordenadas.

      Saludos

    • Bueno, no se si es tarde pero ya en casa te respondo que es necesario encontrar esos parámetros en algún documento colgado en la red y luego, para algunos puntos expresados en ambos sistemas aplicar la transformación con fines de verificación. El ejemplo de validación está aquí:

      https://joseguerreroa.wordpress.com/2010/10/15/validacion-de-transformacion-de-datum-molodensky-badekas-para-antofagasta/

      pero para Antofagasta (con precisiones mejores que un metro). Para Santiago, Valparaiso y Coquimbo habría que buscar, como te dije, sus parámetros respectivos. Invirtiendo algunas horas de búsqueda y algo de suerte (ya que podrían no estar) para hallarlos. Las palabras clave, en Google, serían molodensky badekas santiago valparaiso coquimbo.

      Saludos

  3. Julio dijo:

    Hola, yo tengo algunos puntos de la región de atacama pero al sacar lo 7 parámetros con los programas AllTrans y SEVENPAR, al aplicarlos en el ArcGis me sale corrida la información, no se si estoy insertando bien la información o que pasa, no se si me puedes ayudar.
    Saludos

    • En principio, el ejecutable que está aquí debería funcionar. Si me señalas cuales son los puntos podríamos verificarlo. Para Antofagasta si sé que funciona porque encontre unos puntos de referencia para validarlo. Podría ser que no estás usando adecuadamente la información.

      Saludos

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