INTRODUCCION

El problema de la determinacin de una red geodsica a lo largo de la superficie terrestre, arranca siempre de unos datos de una partida conocida. Con ellos, se resuelven determinados problemas sobre el elipsoide, que permiten establecer coordenadas de nuevos puntos o vrtices sobre l. Para empezar el presente tema, y antes de plantear los problemas que se van a estudiar, se presentaran algunas definiciones, teoremas y desarrollos que intervienen directamente en su resolucin.

EXCESO ESFERICO

= + + 180 .. Grados = . =
648000 . Radianes 648000 . ..... Segundos

Expresando el de un tringulo, en una esfera de radio R, expresado en funcin del rea de ese tringulo esfrico. = + + 180 = =
648000. .2 180. . 2

El valor de la superficie de un triangulo esfrico, expresado en funcin del valor de sus ngulos y de la superficie de la esfera + + 180 = . 720

TEOREMA DE LEGENDRE
Un

tringulo esfrico, cuyos lados sean muy pequeos comparados con el radio de la esfera que lo contiene, puede reemplazar para su clculo, por uno plano cuyos lados tengan la misma longitud que los del esfrico, y sus ngulos sean los de este, disminuidos en un tercio del exceso esfrico.

esfrico

= 3

= 3

= 3

TEOREMA DE GAUSS:

Se define curvatura total de una superficie cualquiera en un punto en el que sus radios de curvatura principales son 1 y 2 =
1 1 .2

En concreto, para que elementos del elipsoide, puedan aplicarse sobre la esfera, se debe verificar que sus curvaturas totales sean iguales. Si es la curvatura total en un punto del elipsoide, y es la curvatura total en otro punto de la esfera. = = .
1 .

1 2

y, por lo tanto:

ESFERA DE JACOBI:
Se define Esfera de Jacobi, a la esfera que es tangente a un elipsoide de revolucin. Consideremos el punto A de coordenadas (, ), y por ltimo, consideremos un punto B de la misma geodsica, de coordenadas (1 , 1 ), en el que el acimut sea 1 = 90.

(1 . 1 ) = (2 . 2 ) = ( . ) =

0. 0 = . = 1 . 90 = ,De donde 0 = . = 1

De forma anloga, pero ahora sobre la esfera de Jacobi, consideremos el arco de circulo mximo . Sea en un punto de este arco de crculo mximo, de coordenadas (1 , 1 ).

0. 0 = . = 1 . 90 = de donde 0 = . = 1 . . = . si = = . = .

=
2 . 2

. 2

. 2

2 . 2

. 2

Y otra relacin del tema 2, que da el valor del radio de curvatura de la elipse meridiana, junto con la recin obtenida:
=
(1 2 ) 3

2 (1 2 ) . 3

2 . 3

. 2 .

. .

= .

DESARROLLO DE WEINGARTENPUISEUX

DESARROLLO DE WEINGARTENPUISEUX

Sea A un punto de latitud sobre el elipsoide, cuyos radios de curvatura principales en ese punto son N y p. Supongamos que desde este punto parte una lnea geodsica con azimut Z, cuyo radio de curvatura en A es . Sea B un punto separado de A una distancia s, medida a lo largo de esa curva geodsica. Pues bien las coordenadas de B, expresadas en el sistema local de las direcciones principales (2) de A son:
s 3 x = . + 6 s 3 y = . + 6

2 s3 z= 2 2 + 2 4

Las coordenadas as obtenidas equivalen a las componentes del vector que tiene origen en A y el extremo B, expresadas en el mencionado sistema de referencia. De estos desarrollos es posible obtener las formulas inversas, que en un proceso iterativo, permiten hallar el valor de la distancia s y el acimut Z, en funcin de las coordenadas x, y z de B expresadas en el sistema local de las direcciones principales de A .Para obtenerlas, se comenzara por sumar los cuadrados de las dos primeras ecuaciones anteriores.
x2 + 2 = 2 . 2 +
s6 2 2 2 36

s4 2 3

2 . 2

s6 2 2 2 36

s4 2 3

x2

s4 2 3

2 +

s6 2 + 36 2

2 + 2

La expresin que aparece entre parntesis en el segundo trmino del segundo miembro, es la inversa del radio de curvatura dado por tanto:

4 6 2 2 s s 2 2 2 x + = 2 + + + 2 2 2 3 36

Despejando :

4 6 2 2 s s 2 2 2 = x + + 2 + + 2 2 2 3 36

La suma + es el mayor termino que da mayor peso al valor de ya que en los otros trminos , N y son muy grandes comparados con s. Por este motivo, se puede sustituir en la expresin anterior, el valor por la suma + .
x 2 + 2 x 2 + 2 x 2 + 2 1+ 2 2 3 36
2

2 2 + 2 2

En esta frmula, s esta expresado en funcin de la incgnita Z, por lo que se ha de buscar otra frmula que ligue ambas incgnitas. Esta se obtendr despejando Z en la primera y segunda ecuacin de los desarrollos.
y
3 6

= =

x 3 6

3 y 6 y (6 3 ) = = 3 x x (6 3 ) 6

De donde se obtiene Z en funcin de la longitud s: y (6 3 ) = x (6 3 )

Despus de realizar observaciones angulares y de distancia sobre el terreno, y tras su compensacin, se han de resolver el traslado de posiciones geogrficas sobre el elipsoide, para establecer una red geodsica en una determinada zona. Estos problemas se denominan problemas geodsicos principales, y son el problema directo, y el problema inverso. Existen diversos procedimientos para su resolucin, de los cuales, se estudiaran aqu el mtodo algebraico y el de las esferas auxiliares. Problema Directo Problema Inverso

METODO ALGEBRAICO
Problema Directo Como ya se sabe, los datos del problema son las coordenadas de un punto A ( , ), el acimut de la geodsica a otro punto B, y la distancia s que los separa. Las incgnitas son las coordenadas ( , ), del punto B, y el acimut inverso . Los pasos a seguir para resolver el problema son los que se sealan a continuacin, de acuerdo a la fig. en la que se representa el sistema local de las direcciones principales (2) con origen en A, y el sistema general cartesiano con origen en O.

Mediante los parmetros a y f del elipsoide de trabajo se calcularan los valores de la primera y segunda excentricidad del elipsoide, e y e, y los siguientes parmetros para la latitud , del punto A. Radio de curvatura del vertical primario. Radio de curvatura de la elipse meridiana. Radio de curvatura () de la normal segn el acimut Z. Mediante los desarrollos de WeiingartenPuiseux, se calculan las componentes del vector referidas al sistema de las direcciones principales de A:
1 = = 3 . + 6() 3 . + 6() 2 3 2 2 + 2() 4()

Se calculan las coordenadas cartesianas del punto A, a partir de ( , ), y h=0, mediante: = = 1 2

Se calculan las coordenadas cartesianas del punto B, que es el resultado de realizar la siguiente suma vectorial en el sistema general cartesiano:

= + = + . = + . 90 + .

= +

1 2

Calculo de las coordenadas geodsicas del punto B ( , ), a partir de las recin calculadas coordenadas cartesianas ( , , ). En la primera de ellas se tendr en cuenta que h=0, por lo que no ser necesario realizar el proceso itera
= 1 2
2 2 +

Se calcula el acimut de la geodsica que parte de A en el punto B, se empleara la frmula de Clairaut aplicada a los puntos A y B . = . = . = El valor de puede ser positivo o negativo > 0 > 0 > 0 < 0 < 900 = > 900 = 1800 < 2700 = 1800 > 2700 = 3600 +

Se calcula el acimut inverso mediante:

= 1800

Si se desea se puede hallar la convergencia de meridianos mediante:

=

PROBLEMA INVERSO:
los datos del problema son las coordenadas de dos puntos A y B Las incgnitas son los dos azimut reciproco ZAB y ZBA entre A y B y la distancia s que los separa De acuerdo al grafico en la que se representa el sistema local de la direcciones principales 2 con el origen en A

De acuerdo a que elipsoide de regencia se trabaja

Mediante los parmetros a y f del elipsoide de trabajo se calcula la primera excentricidad e del elipsoide , y el gran valor de la norma Nb, para la latitud de B de

se calcula el Radio de curvatura NA del vertical primario mediante la expresin y tambin el Radio de curvatura pA de la elipse meridiana

PROBLEMA INVERSO:
los datos del problema son las coordenadas de dos puntos A y B Las incgnitas son los dos azimut reciproco ZAB y ZBA entre A y B y la distancia s que los separa De acuerdo al grafico en la que se representa el sistema local de la direcciones principales 2 con el origen en A

De acuerdo a que elipsoide de regencia se trabaja

Mediante los parmetros a y f del elipsoide de trabajo se calcula la primera excentricidad e del elipsoide , y el gran valor de la norma Nb, para la latitud de B de

se calcula el Radio de curvatura NA del vertical primario mediante la expresin y tambin el Radio de curvatura pA de la elipse meridiana

Se calculan las componentes delos vectores OA Y OB respecto al sistema general cartesiano, o lo que es lo mismo , las coordenadas cartesianas de A (A, A Y B (B, B

se aplica el tensor de giro R correspondiente para que dicho vector de diferencia quede referido al sistema de las direcciones principales de A

una vez calculadas (xb,yb,zb) del punto B segn el sistema de las direcciones principales del punto A. Se pueden hallar las ZAB Y S de acuerdo a las expresiones a8 y a9
segn el mtodo iterativo

Pero para eso tenemos que tomas una valor de inicio lo cual es ZAB=arctan(yb/xb) estonces hallamos Rz(AB) y s

Entonces con este nuevo valor encontramos un nuevo azimut

El valor de Zab dado por a9 puedes ser positivo o negativo y en ambos caso se puede obtener dos valores en tener en cuenta la cantidad de k.Para elegir el valor verdadero para un clculo programado basara tener en cuenta los signo de xb y yb

6:Se calcula el azimut Z BA la formula a10 y se desea obtener tambien la convergencia de meridianas se utiliza la siguiente formula a11 a10: ZBA=ZB+-180 a11: Z=ZB-ZAB

ejemplo Dadas las coordenada del vrtice geodsico cabeza (A) y las del vertice cortijo (B) en el sistema wgs84 se desean calcula los azimut es directo e inverso as como las distancia que hay entre ellos .Los parmetros del elipsoide WGS84 a=6378.137m Y f=1/2982572235

Mediante los parmetros ay f del elipsoide de trabajo se calcula el valor de la primera excentricidad e , NA y NB

Mediante las formula a6 y a7 se calculan las coordenadas geocntricas de A Y B

3: SE calculan las coordenadas del vector AB en el sistema de las direcciones principales del punto A mediante la formula a8

4 : Se calcula el ZAB y S de acuerdo alas formulas a8 y a9 segNu el iterativo Para ello se ha de considerar primero el azimut ZAB y se calcula mediante: ZAB=arctan(yb/xb)=artan(9375908/8331107)=48376817069

Y con este azimut se calcula el primer valor para RZ(AB) CON LA FORMULA

Y CON RZ(AB) se puede calcular un primer valor para para la distancia s con la formula a8

y con este valor se obtiene un nuevo valor ZAB

SE e continua este proceso de calculo iterativo de RZAB Y ZAB hasta conseguir la convergencia deseada tanto en azimut como en distancia

el cuadrante corresponde al azimut ZAB vendr determinado por los signos xb y yb segunda siguiente expresin

9.-METODO DE LAS ESFERAS AUXILIARES

Se basa en el establecimiento de las relaciones entre los puntos del elipsoide y determinadas esferas, de forma que sea posible trasladar a ellas distintos problemas planteados sobre el elipsoide, y resolverlos mediante la trigonometra esfrica clsica.( mximo 60km)

9.1.-PROBLEMA DIRECTO
Se tiene: -Las coordenadas de un punto A (A,A) -El acimut ZAB de la geodsica a otro punto B -La distancia s que los separa Se busca: -Las coordenadas (B, B) del punto B -El acimut inverso ZBA

Trazar el arco de lnea geodsica BQ, perpendicular al meridiano de A. Llamemos a los arcos BQ=x y AQ=y coordenadas geodsicas ortogonales de B respecto de A
= (5 26)

Por la aplicacin del Teorema de Gauss en el punto A se calcula los valores x e y sobre la esfera tangente en A, de radio = , resolviendo el tringulo esfrico equivalente al elipsoidal ABQ. Por la aplicacin del Teorema de Legendre, se calculara el anterior triangulo esfrico como plano, disminuyendo los ngulos A, B y Q en la tercera parte del exceso esfrico.

Si se llama al ngulo B del tringulo plano se cumplir:

1 1 2 = 180 90 + = 90 + 3 3 3

Aplicando el teorema del seno al triangulo plano:

1 (90 3 ) = 2 (90 + 3 ) = 1 ( 3 )

= = 1 2 1 3 ( 3 ) ( 3 )

Despejando
= 2 ( 3 ) 1 3 = 1 ( 3 ) 1 3

2 = ( ) 3

1 = ( ) 3

Desarrollando por el coseno y seno de la diferencia, respectivamente: 2 2 = + 3 3

2 2 2 3 = + = 2 3 3 3 2 = (1 + ) 3 1 1 3

1 1 3

1 1 1 3 = + = 1 3 3 3 1 = (1 ) 3

Para determinar de hallar los valores de x e y, es necesario determinar el valor del exceso esfrico . Para ello, hay q calcular el rea del tringulo ABQ, mediante el tringulo plano correspondiente, despus de aproximarlo al recto en Q, y utilizando las expresiones

2 0 1 1 2 2 1 1 2 3 = = cos = = cos 1 2 2 3 3 2 0 3 La expresin (5 2) da el exceso esfrico en

funcin del rea T, recin calculada, y el radio R de la esfera local en A:

1 2 = 2 = = = cos 2 2 ( )

Donde y son los radios de curvatura principales del elipsoide elegido, en el punto A. Llevando este valor de , a las expresiones se obtiene los valores de x e y:
2 1 2 1 2 1+ = 1 + 2 3 2 3

1 1 2 = 1 + 3 2 1 2 = 1 2 6

2.- Una vez determinado el valor del arco y, se calcula la diferencia de latitudes que lo definen. Para ello, se utiliza la expresin que da la longitud del arco de meridiano:
2 = (1 + 2 2 ) 8

Por ser pequeo

= ,
Sustituyendo: = (1 +

2 2

, = , =
2 2 1 (1 + 2 ) 2 8

2 2 2 ) 8 2

De donde: =

Como 1 + 1/(1 ) cuando x es un infinitsimo, queda como valor para , en radianes: 2 2 = (1 2 ) 8 2

Esta expresin da una primera aproximacin de . Con este valor se podr calcular, tambin en primera aproximacin, la latitud del punto Q mediante: = + Con ella se puede calcular la latitud media entre A y Q para volver a calcular un valor de mas exacto, utilizando una esfera tangente en el punto medio de ambos puntos, para lo cual se volver a emplear la ecuacin pero introduciendo en ella en vez de , y en vez de .

+ 2

(1 2 ) 3

Por tanto, la siguiente expresin ser la que de la diferencia definitiva, con el signo correspondiente, y que habr que introducir, posteriormente, en la expresin =

3.- Calculo del arco (Fig.5 6). Para ello, se traza una esfera tangente al elipsoide en el paralelo correspondiente a la latitud de Q; y de radio la normal NQ para esa latitud. El Polo de esta esfera es PQ. La prolongacin del punto B a travs de su normal hacia esta esfera, determina el punto B. En esta esfera, se resuelve el tringulo PQQB para calcular la diferencia entre los arcos PQB y PQQ que da el valor del arco BT BT=QQ= . Para tener una idea del orden de magnitud de este valor , suponiendo que BQ sea igual a 50km, y =45, se obtiene =190m, y tal como se ha mencionado, se puede suponer situado, tanto en el elipsoide, como en la esfera tangente

Resolucion del tringulo PQQB

Esfera auxiliar para el clculo de

El tringulo PQQB, que es rectngulo en Q, se resuelve mediante aproximaciones a la resolucin de tringulos esfricos rectngulos (6). Los arcos y x vienen expresados en grados. Para expresarlos como magnitudes lineales en una esfera de radio R, se har:
2 = 22 Como el radio de la esfera es 2 = 2 4.- Una vez determinado el valor de , se calcula, al igual que con el arco y, la diferencia de latitudes que define su arco. Para ello, se volver a utilizar la expresin que da la longitud del arco de meridiano teniendo en cuenta las mismas consideraciones que las realizadas para el arco y. El valor de , en radianes, tendr el mismo signo que , siempre positivo. 2 2 = (1 2 ) 8 2 2 = 2

5.- Se calcula la Latitud de B. Como ya se ha mencionado, el valor de , ser positivo o negativo segn el punto B est ms al Norte o al Sur de A, respectivamente. Por ser siempre positivo, se ha de restar, Asi, una vez conocidas y convertidas a grados las diferencias de latitudes que definen los dos arcos e y, junto con la latitud de partida, se obtiene la latitud del punto B: = + 6.- Calculo de la Longitud. La diferencia de longitudes , entre los puntos A y B, se obtendr a partir del valor del arco BQ=x, tomado sobre una nueva esfera tangente al elipsoide, esta vez a lo largo del paralelo de B; sobre la cual, la diferencia de longitud, conserva al mismo valor. El radio de esta esfera, ser el de la gran normal correspondiente a la latitud del punto B. En esta esfera, se resuelve el tringulo esfrico PBBQ correspondiente a una figura similar a la Fig.5 6. Este tringulo ser rectngulo en Q, por aproximacin al episdico PBQ. Teniendo en cuenta estas consideraciones, se resuelve el tringulo rectngulo de la Fig.5 8, expresando todas las cantidades angulares en radianes: 90 = ( 2 ) 1

1 =

Calculo de la longitud 7.- Se calcula la Convergencia de Meridianos, para lo cual se va emplear la Esfera de Jacobi que como ya se sabe del apartado 5, su radio es a, y es tangente al elipsoide a lo largo del Ecuador. En esta esfera, se conservaran los acimuts y iguales a los del elipsoide, siempre que la relacin entre los puntos de ambas superficies sea tal que A y B en la esfera, homlogos de A y B en el elipsoide, tengan como latitud esfrica la latitud reducida de A y B. La distancia s se convierte en s. Segn lo estudiado en el apartado 5, la diferencia de longitudes entre A y B en el elipsoide, y entre A y B en la esfera, en general, no se conservaran. Teniendo en cuenta la definicin dada para la convergencia de meridianos en el Tema 2 mediante la :
= 180

Se resolver el tringulo esfrico PAB en la Esfera de Jacobi, Para ello, se emplearan las analogas de Neper. Veamos con que elementos de la esfera se corresponden estos angulos y arcos: + = + = + 360 = 180 + 360 = 180 = 360 2 = 90 2 + 1 = 90 + = 90 = 90 + = 180 ( + ) 2 2 : = 90 = = = 2 2

La esfera de Jacobi

Sustituyendo estos valores en la analoga de Neper, se obtendr: 2 tan 90 = = 2 cos(90 ) 2 2 2 2 Aunque para distancias cortas las latitudes reducidas por las geodsicas , y el incremento de longitud sobre la esfera, por el incremento sobre el elipsoide (Teorema de Dalby), se realizara el clculo de sin realizar esta sustitucin,. Para ello, ser necesario obtener las latitudes reducidas a partir de las latitudes geodsicas, y obtener el incremento a partir de . El valor de a emplear ser el valor medio correspondiente a la latitud media de A y B: 1 = = = 2 2 1 1

Con estas expresiones, se calcula fcilmente la convergencia de meridianos . El signo de , que va implcito en ella, depender del signo de y de , por tanto, en el hemisferio norte, si > sera > , y por el contrario, si < , ser < . El se toma positivo, si para ir de A a B, la logitud varia hacia el Este; en caso contrario ser negativo.

8.- Una vez calculada la convergencia de meridianos en posible hallar el acimut inverso mediante la relacin

= 180

= + + 180

9.2 PROBLEMA INVERSO

1.- Calculo de los acimuts y convergencia de meridianos Se utiliza la Esfera de Jacobi , tangente al elipsoide en el Ecuador. Se emplear la relacin existente entre los puntos A y B de la esfera con los A y B del elipsoide: Luego se ve que se conservarn los acimuts, del elipsoide a la esfera.

2.- Calculo de la distancia S entre los dos vrtices A y B.

Este calculo se puede realizar tomando como punto de partida, cualquiera de las coordenadas ortogonales geodsicas xe y. Luego utilizando el azimut y para utilizarlo como dato del problema . Se procede a calcular la longitud S mediante la coordenada ortogonal de S.

EJEMPLO:
Dadas las coordenadas del vrtice geodsico A y del vrtice obtenido B, en el sistema WGS84, se desea calcular los acimuts directo e inverso, as como la distancia que hay entre ellos. Los parmetros del elipsoide WGS84 son: a = 6 378.137m y f = 1/298, 257223563

A (Cabeza) = 3622 18 , 6557 = 0605 01 , 9186

B(Cortijo)

= 3626 48 , 77054 = 0602 45 , 42737

1.Calculo de los acimuts y convergencia de meridianos.

Para ello, se han de determinar las latitudes reducidas de A y B, y el en radianes, y previamente .

2 = 2 = 0.0066943800

cos =

= 0.80611347488

= 0.6332056823

= 0.8053599267

= 0.6345139308

+ 2

= 0.6338598065 = = 2 2 = 2 2
1

0.0006541242
= 0.0018292605

Se calcula la Convergencia de Meridianos y el ambas cantidades en grados:

= = 0.0005416984 2 2 2
= 0034347

2 = = 0.8875964422 2
= 4335 31 , 73

Una vez calculados y , se resolver el sistema de ecuaciones cuando >

= + 90 2 = + + 270 2

= 04822 36 , 54 = 22826 20 , 01

2. Para el clculo de la distancia del arco geodsico, se comenzaran calculando y , y .

= 1 2 2

= 0.9988221777 3 = 0.9964706932 = 6.385.658,171

= 1 2 2 (1 2 ) ( 2 2 + 1) =
= ( )

= 6.385.658,975

= 6.357.878,230

= 9.375 912

Se calcula ahora, para obtener ; y luego con , se obtendr la longitud de arco.

= 6.357.918,254

(1 2 ) = 3

2 = (1 + 2 cos 2 ) 8

= 8.326,035

2 = 2

= 5.083

Una vez determinados los arcos y , se puede hallar el valor del arco, siendo su valor positivo por tener el punto B mayor latitud que A. Finalmente, se calculara la distancia entre los puntos A y B. = + = 8.326.065 + 5.083
= 8.331,119

1 2 = 1+ 6

= 12.542,536