Radios de Curvatura

ECUADOR ELIPSOIDAL
o G
Cálculos
de
Posicionamiento
Geodésico
INSTITUTO NACIONAL DE ESTADISTICA

GEOGRAFIA E INFORMATICA
ECUADOR ELIPSOIDAL
Cálculos
de
Posicionamiento
Geodésico
INSTITUTO NACIONAL DE ESTADISTICA

GEOGRAFIA E INFORMATICA
NorA DEL TRADUCTOR
Este trabajo es parte del esfuerzo que DETENAL está haciendo
con el prop6sito de elevar el nivel de los conocimientos geodésicos
dentro y fuera de la Instituci6n. Debe pues, agradecerse la ､ｩｳｰｾ＠
sici6n y el apoyo brindado por las autoridades de DETENAL, parti-
cularmente de aquellas responsables del área de Geodesia, que al
facilitarnos medios y personal, hicieron posible que éstas notas ｶｩｾ＠
ran la luz del día.
Debe agradecerse profundamente la gentileza de los autores,
Dr. E. J. Krakiwsky y Dr. D. B. Thomsón, de la Universidad de
New Brunswick, al dar su autorización para la traducci6n y divulg.!
ci6n de su obra en español.
El excelente trabajo de mecanografiado estuvo a cargo de la
Srita. Concepción Vega Ché de la Oficina de Apoyo Básico. El Sr.
Julio Bueyes Oliva, de la Oficina de Apoyo Vertical, tuvo la respon-
sabilidad de trazar los diagramas y escribir las f6rmulas que ｡ｰｲｾ＠
cen en éstas notas.
M. C. Rafael Sosa Torres.

PREFACIO
El propósito de estas notas es dar la teoría y el uso de algu-

nos métodos de cálculo de posición geodésica de puntos sobre un elipsoi-
de de referencia y sobre el terreno. La justificacióJi por las primeras -
tres secciónes de estas notas de lectura, que tratan el problema clásico
del "Cálculo de posición geodésica sobre la superficie de un elipsoide",
no es fácil de determinar. Solamente puede ser establecido que la inten
ción ha sido producir un paquete independiente que contenga el desarrollo
completo de algunos métodos representativos que existen en la literatura.
La última sección es una introducción a los métodos de cálculo tridimen-
sionales y es ofrecida como una alternativa al método clásico. Muchos-
problemas y su respectiva solución son presentados.
El enfoque que se da aquí es realizar derivaciones completas,

esto es, alejado de la práctica de dar una lista de fórmulas a usarse en -
la solución de un problema. Se espera que este enfoque dará el lector
una apreciación para el fundamento en que están basadas las fórmulas y
al final, las fórmulas mismas.
Las notas se desarrollaron de las conferencias preparadas -

por E.J. Krakiwski y del trabajo de investigación realizado por D. B. --
Thomson en afios recientes en U. N. B. Los autores reconocen el uso -
de ideas contenidas en las conferencias de los profesores Urho A. Uotila
y Richard H. Rapp. del Departamento de Ciencia Geodésica de la Univer
sidad del Estado de Ohio, Columbus, Ohio. Otras fuentes usadas para :-
detalles importantes están referidas dentro del texto.
Los autores desean reconocer la contribución hecha por los

estudiantes de post grado de Ingeniería Topográfica de 1975, para el me-
joramiento de estas notas por el descubrimiento de errores topográficos.
INEGI. Cálculos de Posicionamiento Geodésico

Mr. C. Chamberlain es particularmente reconocido por su -
crítica constructiva y su asistencia en la preparación de este manuscrito
para su publicación.
E. J. Krakiwsky
D. B. Thomson
Febrero 14 de 19 7 4.
TABLA DE CONTENiDOS
Página
Prefacio Ｍｾ i
Lista de ilustraciones - -- - -- - - - - - --- - - - - - - --- -- - - - - -- -- - - - - - --- ii
Introduccion ------------------------------------------------- 11
Sección 1: Geometría Elipsoidal - - - --- -- - --- - - ---- --- - --- - - -- --- 12
1. - El Elipsoide de Révolucion --------------------------------- 12
1. 1. - Parámetros Elipsoidales ----------------------------- 12
1.2. - Radio de Curvatura -------:- ---- ------------ ---- --- --- 14
1. 2.1. - Radio de Curvatura del Meridiano ------------------- 15
1. 2. 2. - Radio de Curvatura del Primer Vertical-------------- 18
1. 2. 3. - Radio de Curvatura en cualquier Azimut -------------- 20
1. 3. - Curvas sobre la superficie de un Elipsoide ------------- 24
1. 3. 1. - La Sección Normal -------------------- ------------ 24
1.3.2. - La Geodésica ------------------------------------- 26
Sección n. - Reducción de Observaciones Geodésicas Terrestres --- 30

2. - Reducci6n a la Superficie del Elipsoide de Referencia ------ -- 30
2.1. Reducci6n de Direcci6nes Horizontales (o ángulos) -------- 31
2.1. 1. - Efectos Geométricos ------------------------------- 31
2.1. 2. - Efectos Gravimétricos ----------------------------- 35
2.2. - Distancias Zenitales --- ---- -------- ---- Ｍｾ 36
2.3.:" Distancias Espaciales -------------------------------- 36
2.4. - Reducción de Cantidades ｇ･ｯ､ｾｳｩ｣｡＠ calculadas al terreno 37
Sección 1II. - Cálculos de Posiciones Geodésicas sobre el ElipSOide

de Referencia. ----------------------------------- 37
3. - Fórmula de Puissant. - Líneas Cortas --------------------- 39
3. 1. - Introducci6n -- - ---- -------- ---- ----- ---------- ------- 39
3.2. - Problema directo ------------------------------------ 39
3. 3. - Problema InverSO - --- --- ----------- - ------- - --- -- ---- 42

3.4. - Resumen de Ecuaciones para la soluci6n de Problemas Di-
rectos e inversos usando la F6rmula de Puissant. ---- --- 45
3.5. - F6rmula de Gauss! - Latitud Media -------------------- 46
3.6. - Otras fórmulas de línea eorla-------:------------------- 46
4. - F6rmula de Bessel. - Líneas Largas ----------------.------ 47

4. '.- Introducci6n ----------------------------------------- 47
4.2. - Relaciones Fundamentales ------ ---------------------- 47
4. 3. - Soluci6n de la Integral Elíptica -- --- - -- -- -- -- -- - - - - - --- 52
4.4. - Probrema Directo ------------------ ------------------ 58
4.5. - Problema Inverso --------------------------------- --- 59
4.6. - Otras fórmulas para líneas largas - ------------------- 61
ii
TABL.A DE CONTENIDO (CONT)
Pigina
Secci6n IV. - Cálculo de Posiciones Geod6sieaaen Tres
ｉｩｪｾｯｮ･ｳＭ 61
5. - Problemas Directo e Inverso en Tred DlmeDSianes ------ 62
5.1. -Problema Directo ----------------------------- 62
5.2. -Problema InV1!lrBO ----------------.----------- 66
6. - Problemas de Intersecci6n en Tres Dimensiones ------- 67
6.1. - Inter. .cci6n Asimutal ------------------------ 67
6.2. - Intersecci6n en las m-.ncias Espaciales ------ 73
7. - A manera de concluswn ---------------------------- 75

Referencias ---------------------------------'------- 77
Ili
LISTA DE ｉｌｕｓｔｒａｃＰＱｾｅ＠
FIG. No. PAJINA

1 El Elipsoide de Revolución 13
2 Sección Normal Meridiana Mostrando el -

radio de curvatura del meridiano (M) • 15
3 Sección Normal al primer vertical mos -

trando la curv:atura del primer vertical (N) 16
4 Radio de curvatura del meridiano (M) 17
5 Radio de curvatura del primer vertical (N) 19
6 Sección normal en un Azimut Cualquiera C! 21
7 Indicatríz para la solución de Ra 21
8 Sección a través de PP' (o) para la solu -

ción de Ra 22
9 Solución de Z para solucionar Ro 23
10 Secciones con normales recípr,o.cas 24
11 Sección Triangular con normal recíproca 25
12 Separación Angular entre secciónes con Nor

mal Recíproca. - 26
13 Geodésica 27

14 Ecuación Diferencial de una Geodésica sobre
la superficie de un Elipsoide de Revolución. 29
15 Separación entre Sección Normal y Geodésica 3D
16 Corrección por Normal Oblicua 32
17 Corrección por Desviación de la Vertical 34
18 Reducción de Distancia Espacial 36
19 Fórmula de Puissant para Problema Directo. 39
iv
FIG. No. PAGINA
20 Fórmula de Puissant para Problema Inverso 44
21 Relaciones fundamentales para el desarrollo

de la fórmula de Bessel. 49
22 Esfera Reducida y Elipsoide. 49
23 Solución de 52
24 Solución de d1./ dA 57
25 Soluci6n de Longitud de Arco (J'" 59
26 Problema Directo (Geodésico Local) 63
27 Problema Directo (Astron6mico Local) 65
28 Vectores Unitarios en el Sistema Geodési-

co Local 68
29 Intersecci6n Azimutal en Tres Dimensiones. 70
30 Intersecci6n de Distancia Espacial en Tres

Dimensiones. 74
INTRODU CCION
Las primeras tres secciones de estas notas tratan del cálculo

de posiciones geodesicas sobre un elipsoide. En el primer capí-
tulo, se da una revisi6n a la geometri'a elipsoidal ya que en el desa
rroUo de f6rmulas posteriores puede ser muy útil. Común a to :-
dos los cálculos elipsoidales clúicos es la necesidad de reducir
las observaciones geodésicas al elipsoide. por lo tanto un capítU-
lo completo es empleado en éste t6pico.
Dos clásicos problemaa de cálculo de geodesia geométrica -

son tratados; Ellos son los llamados problemas geodésicos "Di -
recto'" e "Inverso n • ,"
Hay varios métodos que pueden adoptarse para la solucifut de

estos problemas. Generalmente están claslficados en término de
f6rmula de línea "Corta" "mediana" o "Larga". Cada uno de
ellos involucra apro:ld.maaiones diferentes que restringen la dis-
tancia entre estaciones sobre las que algunas f6rmulas son muy -
útiles para una exactitud dada.
La última secci6n de las notas trata de los cálculos de posi -

ciones geodésicas en tres dimensiones. Primero son desarrolla
dos los problemas directo e inverso y luego dos problemas. espe-:

ciales (ambos de intersecci6n de Azimut y de distancia espacial )
son tratados. Estas soluciones ofrecen una alternativa a los mé
todos clúicoe de c6leuloe de posici6n geodéaiea. -
11
SECCJOK 1: GEOMETRIA ELIPSOTDAL.
l. - El Elipsoide de Revolución.
Ya que un elipsoide de revolucion (Elipsoide de Re-

ierencia) es generalmente considerado como la mejor aproxima-
ción al tamafio y a la forma de la tierra, es usado como la super-
ficie sobre la cual se hacen los cálculos. Inmediatamente des -
pués ｾｳｴｵ､ｩ｡ｭｯ＠ muchas propiedades geométricas de un elipsoide
de revolucion que son de especial interés para los geodestas.
En particular, el radio de curvatura de puntos sobre la

superficie del elipsoide y algunas curvas sobre su superficie, son
descritos.
1. 1. - Parámetros Elipsoidales.
La figura 1 muestra un elipsoide de revoluciono Los pa

rámetros de un elipsoide de referencia que describen su tamat'l.o y
su forma son:
i) El semieje mayor (a)

ii) El semieje menor (b)
La ecuación de cualquier curva meridiana (Intersección

de un plano meridiano con la superficie del elipsoide) es: (Ver -
Fig. 1).
X2 Z2.
- +--=1------- 1
2
0 b2
La superficie de un elipsoide de revolucion está dada -

por:

Z2
+--=1-------10
b2
12
z
p
ｅｾＭ［Ｋ＠ｅｾｖ＠
p'
FIGURA 1
11 EL ELIPSOIDE DE REVOL ueION 11
Los puntos F y F r en la figura 1 son los focos de la

elipse meridiana que pasa por los puntos p. E' • p'. E. Los
focos son equidistantes del centro geométrico (ó) de la elipse.
Las distancias PF y PF' sonigual al Semi eje mayor (a). Es
ta información es usada ahora para ayudar a describir propie ::-
dades posteriores de un elipsoide.
El achatamiento elipsoidal (Polar) está dado por:
a-b
f=-a- ------2
otras dos importantes propiedades que son descritas

para una sección meridiana· de un elipsoide son: La primera -
excentricidad.
13
y la segunda excentricidad:
Como un ejemplo de las magnitudes de éstos paráme-

tros para un elipsoide de referencia geodésico, presentamos --
aqui los valores del Elipsoide de Clarke de 1866 que es usado en
el presente para la mayoría de los cálculos de posición geodési-
ca en Norte América. (Bomford 1971, p 450):
a= 6378206.4m
b= 6356583.8m
Usando 2
f= 0.00339006
El cual es dado a menudo en la forma de 1-, que en -

este caso es ';. = 294.97869 f
f
Usando 3 y 4 respectivamente, tenemos:
e2 = 0.00676865----
e,2 = 0.00681478----
Los cuatro parámetros a, b, e, ( o e' ) y f; Y las re-

laciones entre ellos son los principales usados en el desarrollo de
más fórmulas geodésicas.
1. 2. - Radio de Curvatura.
Sobre la superficie de un elipsoide un número infinito de

planos pueden dibujarse a través de un puntos sobre la superficie-
que contiene la normal en ese punto. Estos planos son conocidos
como "Planos Normales". Las curvas de intersecci6n de los pla-

nos normales y las superficies del elipsoide son llamadas "Seccio
nes Normales " • En cada punto hay dos secciones normales mutua -
mente perpendiculares cuyas curvaturas son máximas y mínimas-
y son llamadas las "secciones normales ｰｲｩｮ｣｡ｬ･ｳｾＧ＠ Estas seccio
nes principales son las "Secci6nes Normales del Meridiano" y deC
"primer vertical" y sus radios de curvatura son denotados por (M)
y (N) respectivamente (figuras 2 y 3). En la figura 2 puede verse
que el radio de curvatura del meridiano aumenta del ecuador al po-
lo y el radio de curvatura, del primer vertical se comporta similar
mente (Fig. 3). La razón de esto será visto pronto, una vez ｱｵｾ＠ :-
hayan sido desarrolladas las fórmulas para (M) y (N).
14
1.2.1. - Radio de Curvatura del Meridiano.
Consideramos una secci6n meridiana de un elipsoide

de revoluci6n (Figura 4) dado por:
El radio de curvatura de ésta curva en cualquier punto

"p" está dado por:
ｾＭ PLANO TANG&ITE
FIGURA 2
SECCION NORMAL MERIDIANA MOSTRANDO

EL RADIO DE CURVATURA DEL MERIDIANO

(M).
15
FIGURA 3
SECCION NORMAL DEL PRIMER VERTICAL

MOSTRANDO EL RADIO DE CURVATURA DEL
PRIMER VERTICAL (N).
(Ph1Ups. 1957. pp. 194-197)
M: [1+ IF)Z]"I_______ 5
./z
dT
En el caso de una elipse meridiana!

2
dz X b
-
dx
:--
Z
--------6
al'
dz
17: --¿-( Z-X"""""
.2 Z ...2
Z:X ) --------7
dlz bl
- - : - - - (Z+-o--)----7a
x2 b2
dx 2 , a 2 z2 Z 02
16 ,
De la figura 4, podemos ver también que la pendiente
de la tangente en "P" está dada por:
dz
,1..
tan (90 0 +.,.) = -:¡x:: - cot q, --- ------ 8
Jraalando 6 y 8 tenemos:
X .
-cotq,:-T 7 ｾ＠ ----------------9
o tanq,:: ｾ＠ ! -------------90
Sustituyendo J
Itz
bao (I-e')----------- 9b
en 9a; tenemos:
Zo: x( 1-.') ton q, ------- - - - 10
Entonces. después de sustituir ｾＩ＠ y (z) en 1 con-

9b Y 10 respectivamente, y de algunas simples manipulaci6-
nes nos resulta que:
z
90"+ x
FIGURA 4
RADIO DE CURVATURA DEL MERIDIANO (M)

Sustituyendo la expresi6n anterior para X en la ecuaci6n
10, da la f6rmula:
z= o(t-l} sen! ------- 12
{t - r sen'cp} ｾＲ＠
Finalmente, ｳｵｴ￭ｹ･ｮｾ＠ (X) y (Z) en (6) y (7a) y colocan

do estos valores en (5) para ｾｺ＠ y --}- la expresi6n para el radio de ::
curvatura del meridiano 8e ｣ｯｮｾ･ｲｴＺ＠ d .
o(t-r)
M =-':;":":'-=-'''---- - - - - - - 13
En la ecuaci6n (13) el único parámetro variable es la latitud

geodésica 4> por lo que en el ecuador ( cp= c:1" );
M=a(t-.z }------ 130

yen el polo (cp=90 0 )
M= o --- 13b
ＨｈｦＩｾｴ＠
El radio de curvaturadél meridiano aumenta en longitud a -
mettida que el punto sobre el meridiano se mueve del ecuador al polo.
1.2.2. - Radio de Curvatura del Primer Vertical.
De la figura 5. -

coscp=+ -------14
6;
X
ｎ］ｾＭ 140
18
z
PARALELO DE LATI
ｾＭｲｦＷｌ ___ x
FIGURA 5
RADIO DE CURVATURA DEL PRIMER VERTICAL (N)
Sustituyendo la expresi6n para X (11) en (14a), nos da la ex-

presi6n para el radio de curvatura del primer vertical:
ｎｾ＠
Como solamente el parámetro variable en (15) es el> ,N --

variará entonces con el> • Cuando el> '" O" (ecuador) ｾ＠ N 11 8.; Y cuando -
el> =90· (polos)
N"'al (1-."'1. 11 M ------- 15a
Una cantidad importante que es usada muy a menudo en cál -

culos de geodéáiá geométrica es el "Radio Medio Gaussiano de Curvatura" --
que está dado por:
R = fo --- ----16
19
En muchas instancias el radio medio es suficientemente pre
ciso para cálculos de posición. -
Otro radio de curvatura que puede ser necesitado de vez en

cuando, es el de un paralelo de latitud. Cualquier paralelo de latitud visto _
desde el polo norte del elipsoide (eje Z) describe un círculo. Este radio, _
comp puede verse en la Fig. 5. es igual a la coordenada X (en el sistema ,_
del plano meridiano X - Z). Entonces, de la ecuación (14a) el radio de cur-
vatura de un paralelo de latitud está dado por:
IR cp = N cos ｃｰｾＭ - ----- 17
Se ve fácilmente que cuando cp=OO (ecuador); RCP=N ; Por

10 tanto Rcp=a ya que N = a para cp=O") y en cualquier polo ( cp = 90°):
cos cp =O; Y el radio desaparece
1. 2. 3. - Radio de Curvatura en un Azimut Cualquiera.
Como se mc:lEtr6 en las secciónes 1. 2. 1 • Y 1. 2. 2. El radio

de curvatura máximo y mínimo de un punto cualquiera P, sobre la superficie
de un elipSoide de revoluci6n está en los planos meridianos y primer vertical
respectivamente.
En algunos casos, los cálculos geodésicos requieren el radio

de curvatura en otro plano que no sea el principal (fig. 6),. La secci6n no!.
mal en algún azimut a tiene un radio de curvatura en un punto cualquiera -
P, designado por Ra Este se resuelve usando el Teorema de Euler (Lipa-
chutz, 1969 pago 196) Y es llamaQo el 'tRadio de Curvatura de Euler".
En la fig. 6, el punto P en el cual es requerido el radio Ra ,

está mostrado sobre la sección normal PP: solamente una parte diferencial
de la curva de secci6n normal (ds) es mostrada, de manera que el azimut a

de ésta secci6n pequefia es equivalente al azimut de una sección normal de -
cualquier longitud.
El teorema de Euler se resuelve como sigue. En el punto -

P, dibujamos un plano tangente y paralelo a el, otro plano (fig. 7) que inter-
secta la superficie del elipSoide. El ru.timo plano visto a través de la normal
que pasa por P, forma una elipse en el plano BB 'donde el plano tangente in -
tersecta la superficie del elipsoide. Los elementos de ésta'\ndicatriz" son -
mostrados en la fig. 7. Si vemos este plano a través del punto P: en el azi-
mut a , la sección resultante es la Fig. 8, recordamos que la ecuación de -
una elipse es:
Xl yl
-al- +--=
b2
1---------1
20
z
ｲＭｾｹ＠
FIGURA 6
SECCION NORMAL EN UN AZIMUT CUALQUIERA
:;:=p== Ｇｾ＠ _ •
ｏｾｂＧ＠
A
B:
I I
I I
I
I y
I
I
I
n
ｲＭｾｸ＠
m
FIGURA 7
INDICATRIZ PARA LA SOLUCION DE Ra
21
p
p'
o
Ro
FIGURA 8
SECCION A LO LARGO DE PP' a)
PARA LA SOLUCION DE Ra
De la Fig. 7.
X=ds sen a
-Y=ds cos a------- 17
Entonces (1) se convierte;
Usando la Fig. 9, podemos escribi-r.
Sen. e=.l
e
- - - - - - 19

y
2- e
Sen. e=-=.2Ra_ _ - - - - - - - - 200
De lo que resulta:
c
z=-
2
----------21
2Ra
De manera que si pp'es W1a distancia diferencial muy pequena entonces

C ｾ＠ ds y podemos escribir.
. ds z
Z = - - ----------- 22
2Ra
22
• n2
Cuando a =: 0° ; (S) es igual a (n) y Z=2M----- 23
. 2
Cuando a = 90° (S) es igual a (m) y z = _m__ ------- --24
2N
Combinando (22) Y (23); y (22) Y (24) tenemos:
2
n = ｾＧｍＭＲＵ＠ dr
y
_..:ds=-2_.N _________ 26
Ra
p'
FIGURA 9
SOLUCION DE Z PARA LA SOLUCION DE
Ra
Sustituyendo n 2 y m2 en (18) da:
Ra sen 2 a + Ra cos2 a = 1- __________ 27

N M
Finalmente después de arreglar los términos de (27), tenemos la expresión
para el radio de curvatura de Euler.
M. N
Ra = ---=-----:---.r-- -- --- -- 28
M sen 2 a +Ncos 2 a
23
1. 3. - Curvas sobre la Superficie de un Elipsoide.
Hay dos curvas principales sobre la superficie de un Elip -

soide, que son de especial interés en geodesia geométrica. Ellas son la - -
"Curva de Secci6n Normal" y la "Curva Geodésica" descritas abajo.
1.3.1. - La Secci6n Normal.
En la Secci6n 1. 2 la sección normal fue definida como la --

línea de intersección de un plano normal (en un punto P) y la superficie del-
elipsoUie. Considerese dos puntos sobre la superficie de un elipaoide (Pl y
P2) que están sobre meridianos diferentes y sobre diferentes latitudes. La
Secci6n Normal de Pl a P 2 (Secci6n Normal Directa) no coincide con la - -
secci6n normal de P a Pl (Secci6n Normal Inversa). Figura 10.
2
El plano normal de la Secci6n Normal Directa conteniendo -
los puntos P 1 , n-1 y P J contiene la normal en PI yel plano normal inver-
so, P 2 n2 PI' contiene ｾ｡＠ normal en P 2 y el punto PI' Si las Secci6nes --
Normiles PI P 2 Y P2 PI fueran coincidentes, entonces las normales
PI n 1 y P2 n..2' en SUS respectivos planos meridianos podría intersectar
el eje menor en erm1smo purlto.
FIGURA 10
SECCIONES DE NORMAL RECIPROCA.
24
Puede mostrarse que el punto de intersección Zn de cu al -
quier sección normal elipsoidal intersecta al ej e menor en:
(Zacatov, 1 9 53 p 39 - 40).
Si dos puntos tienen longitudes diferentes y CPP < cpP2

(Fig. 10) entonces Zn, < Zn2 y las normales P, n, y ' P. n no
están en el nrlsmoplano. Ellas se dice que son "Normales ｓ･ｳｧ｡ｊｾＮ＠ Sin-
embargo, siCPPI' es igual a<pP2' las secciónes normales directas e inversas
son coincidentes. Para dos puntos sobre el núsmo meridiano, las normales
elipsoidales no intersectan en el núsmo punto sobre el eje menor. Ellos es-
tán, sin embargo, en el núsmo plano (el plano meridiano común), por 10 tan-
to, las secciónes normales PI P2 Y P 2 PI son coincidentes.
El resultado de lo predicho, es que sobre la superficeie del

elipsoide, la sección normal no da una línea única entre dos puntos. Por 10
tanto un triángulo elipsoidal no está definido unicamente por las secci6nes -
normales. En la Fig. 11, la secci6n normal directa de A hacia B, A a B;
no es coincidente con la secci6n normal inversa B b A. Por 10 tanto, el --
Azimut geodésico aA no se refiere a la misma curva de aBo Un proble-
ma ｳｩｭｾ｡ｲ＠ existe para los Azimutes A a e, B a C, etc.
Ahora veremos rápidamente a la magnitud de la separaci6n -

entre las secciones normales directa e inversa. En la figura 12 ésta separa-
ci6n es mostrada como el ánguloll.. La fórmula para la solución de ll. está
dada por: ( ZOkotov, 1953 P 5' )
FIGURA 11.
TRIANGULO DE SECCION NORMAL RECIPROCA.
25
z
ｾ＠ .... ",
" ,,
" ....
.... ,
, IIP
" .... ....'
,
....
FIGURA 12
SEPARACION ANGULAR ENTRE SECCIONES DE

NORMAL REeIPROCAS.
Donde. 4>,n= 4?1 + P,.z ---- ----- 31

2
y S
0-=
Nm
y Nm= NI + NI ----- ----- 31a.

2
Por ejemplo, una línea P1 P2 con 200 km. de longitud y --

condici6nes máximas ( 4+rF 0° y a..z= 45° ): 6. =0'.'36. Ya que mu-

chas poligonales-· '0' líneas de triangulaci6n son más cortas que ésta y -
debido a que la situación máxima no siempre ocurre, el valor de 6. es ge- .
neralmente una cantidad pequefta y en última instancia, prácticamente des -
preclable.
1.3.2.- LaGeodésica.
La geodésica o línea geodésica entre dos puntos cualquiera

sobre la superficie de un elipabide es la curva única sobre la superficie en-
tre los dos puntos. En un punto cualquiera a lo largo de la geodésica, el ra
dio principal del vector de curvatura es coincidente con la normal elipsoidái..
La geodésica (Fig. 13) entre dos puntos P 1 , P2 es la ､ｩｳｴｾ＠

cia menor sobre la superficie entre esos dos puntos.
26
L a posición de la geodésica respecto de las secciones norma -
les directas e inversas se muestra en la Fig •.. 13.
Para describir matematicamente la geodésica, desarrollaremos

las ecuaciones diferenciales para líneas geodésicas, sobre una superficie
de revolución; la diferencial básica geométrica requerida para esto puede
fundamentarse en Phillips 1957 y Lipschuts 1969 La ecuación ge-
neral para una superficie de revolucion puede ser expresada como:
FIGURA 13
GEOOESICA
F(X,Y,Zl: 0--------- 32
Las ecuaciones paramétricas para una geodésica sobre esta supe!.
ficie son:
Los cosenos directores de la normal a la sup erficie son:

aF aF OF
Coso Q:
1-'1
OX
o
. cos o :. OY
' 1-'2 o'
. cos {3.: aZo - - - - - - 34
I
• [aF
Don d e.o: F
(--1 + (--) 2 a Z o 1] ｾＲ＠
F
+ (--) --- - - - - - -35
ax ay az
Los cosenos directores de la norI!B.l principal a la curva (33) son:
. cos/3. =R.JÓ:......· cosf3. ］ｒｾＭＳＶ＠

d2X
cosf3 :R
"1 dS2' N2 dSt' NS ds"
27
Donde R es el radio principal de curvatura de la superficie. En
la definicion de la geodésica se estableció que en cualquier punto sobre la
curva, la normal a la superficie y el radio vector principal (normal princi
pal) están siendo coincidentes. Para satisfacer esto, igualamos ( 34 ) Y
( 36 ) lo que se reduce a:
OF OF OF
ｾ＠
OX ay az
=--¡:;-= /z -----------37
ds
2
7 dT
Como estamos tratando con un elipsoide de revolución, la super-
ficie la podemos"tepresentar por la ecuacion:
X' + y2 + f tZl=O --------- 38

Entonces:
OF OF
--=2X·--=2Y· --=HZI------39
OF,
Ox ' OY , OZ
Las cuales cuando se sustituyen en (37) tenemos:
2
d x /y
Y - - - X - - = 0--------- 40
.
ds! di'
Integrando (40) tenemos:

Ydx- Xdy = CcIs----- ----- 41
Donde C es la constante de integración.
En la Fig. 14, la lfnea PP' representa una parte diferencial de una geodé-
sica sobre la superficie del elipSoide. Teniendo las coordenadas cartesia-
nas de P ( X, Y,Z ) podemos calcular las coordenadas de pi ( X+dx¡ Y+dYi
Z+dz ) ya que ds es una distancia muy pequena. Las coordenadas de A
(Proy.ecci6n de Pi) (sobre el plano del. paralelo de latitud de P) son en-
tonces X+dx: Y+dy. Z • El radio de éste paralelo es simbolizado por

r. El area del triángtÜO CPA es:
Area CPA ·t(YdX-XdYI----------42
yel área del sector PP" Ces:
Area PP .. C • t r.ds. sen a ------- 43
Cuando ds es muy pequefto:
Area CPA ,. Area PP" C:
Entoncesl
28
y sustituyendo (41) en (44) tenemos:
r----------------------y
z
FIGURA 14
ECUACION DIFERENCIAL DE UNA GEODESICA SOBRE

LA SUPERFICIE DE UN ELIPSOIDE DE REVOLUCION.
Cds=r. sen a ds--------45
o r sen a = e -- - -- - - - - -- 46
Finalmente, sustituyendo (17) en (46) encontram<B que:
IN coscp sena =cl--------- 47
Para cualquier pilnto a lo largo de una geodésica sobre la super-
ficie de un elipsoide de revoluciono En cálculos de geodésia geométrica,
es necesario definir nuestros azimuts directo e inverso respecto a la mis-
ma curva sobre la superficie y no respecto a las dos secciones normales. -
Por lo tanto necesitamos la separación entre la sección normal y las cur -
vas geodésicas. La separación, definida aqui sin comprobación, está dada
por (Zacatov1953 p.p. 41-51) .
8= ｾ＠ -------------48
Donde 8 es el ángulo entre la sección de normal directa y la geo-
désica en cualquier punto, y t:. es el ángulo entre las secciones normales
recíprocas entre dos puntos. Desarrollos posteriores de esto y la aplica -
ción de correcciones apropiadas son dadas en 2.1.1.
Luego, la distancia S entre dos puntos sobre la superficie de un

elipsoide es diferente si uno usa sección normal mayor que la geodésica. La
diferencia está dada por (Zacatov 1953, p. 51).
29
La ,!'Se para W1a lÍnea de 600 kms. en longitud es aproximadamen-
te de 9 X 10m. lo que es obviamente despreciable para todos los prop6
sitos prácticos. -
ｾ＠
FIGURA 15
SEPARACION ENTRE LA SECCION NORMAL y LA

GEODESICA.
SECCION II: REDUCCION DE OBSERVACIONES GEODE-

SICAS TERRESTRES.
2 • - Reducci6n a la Superficie del Elipsoide de Referencia.
Mediciones geodésicas (direcciones terrestres, distancias, dis-

tancias zenitales) son hechas sobre la superficie de la tierra. Los cálcu-
los de las posiciones geodésicas son hechas sobre el elipsoide de referen -
cia. Por lo tanto las mediciones deben ser reducidas de la superficie de la
tierra al elipsoide de referencia. Cuando se reducen cantidades medidas,

hay dos conjW1tos de efectos que deben considerarse - efectos geomét'ricos
y los efectos de las variaciones en el campo de gravedad de la tierra.
Debe notarse que las reducciones desarrolladas aquí pueden apli-

carse de W1a menera inversa. Esto es, calculadas las cantidades geodésicas
elipsoidales (distancias, por ejemplo), pueden ser "Reducida" a la superficie
de la tierra (2.4.).
30
2.1. - Reducci6n de Direcciones Horizontales (o ángulos).
Cuando medimos direcciones sobre la superficie de la tierra, ni-

velamos el instrumento hasta asegurar que el eje vertical Sea coincidente -
con el vector de gravedad local. Sabemos que el vector de gravedad local y
la normal al elipsoide, generalmente no coinciden. Para referir las direc-
ciones a la normal elipsoidal, es necesaria una correcci6n por desviaci6n -
de la vertical.
Otras dos consideraciones son aquellas de la geometría elipsoidal.

Primero, las normales en dos puntos sobre un elipsoide son "desviadas" una
con respecto a la otra, asi cuando un punto está sobre el elipsoide, este pun
to no está en el mismo plano de la proyecci6n normal del punto sobre el eliE
soide.
La correcci6n asociada con este fen6menose llama correcci6n de

la normal desviada. Seguida:meate deseamos tener las direcci6nes geodési-
·cas y no las direcciones de las secciones, normales, por lo que se necesita
una correcci6n geodésica - secci6n normal.
2.1.1. - Efectos ｇ･ｯｭｾｴｲｩｳ＠
La Fig. 16 muestra la situaci6n sobre la superficie de la tierra -

de mediciones de direcci6nes, después .de que han sido quitados los efectos
de gravedad (2.1.2. ) En ésta figura, P1', es la estaci6n de tnedida, que es-
tá.sobre la normal P 1 n1' P'2 es visado a la altura h sobre el punto elipso.!.
dal P 2 • Si h2 = O, la direcci6n medida (mostrada agui oomo un azimut, - -
Le. a 12 ;' d12 + ZII. donde Z12 se asume aue es :un parámetro de orientaci6n co-
nocida). Será entre los planos P,Z", y P, p. ｾＢ＠ esto eS"'li1.\2' el azimut de la
secci6n normal directa. En la práctica h 1: O, por lo que la direcci6n medida
a12 debe ser corregida. La reducci6n por éste efecto, llamada la reducci6n
por'normal derivada" 6 por "altura del punto visado", debe ser aplicada.
De (29)
"1 "2 = a el (cp,-4>¡lcos cpm - - - - -- - - ｾ＠

y
(cp _cp 1 = S cos a l2 --- ---- -- 51
I I Mm
Donde
Mm = MI + M,
2
Tenemos:
-
""I=ae
I
2 S
-cosa
Mm 12
cos cp m-------52
Donde S es el arco de longitud P1 P2.

31
Ahora para encontrar la reducci6n ah procedemos como sigue:
Primero calculamos:
FIGURA 16 nz
CORRECCION POR NORMAL SESGADA

R"'=","I cos epi: o ･Ｒｾ＠ Mm
ｊＢｾＲ＠ --
cos a ,.... "'2
Á. -- --- 53
.Donde cpm ha sido sustituida por epz • Ya que la diferen-
cia nos da un efecto despreciable. Entonces, el ángulo en ｰｾ＠ está dado por:
I
oel S cos a'l tOS epi
do: ＭＮＺ［］Ａｾ
Mm' P,:1t
---- 54
Ahora si aproximamos la longitud Pz R con el semieje mayor a, (54)

se convierte:
do - e2 ＠ｾ Mm cos a 12 cosl Á. -
't"l
- - - - - 55
Ahora, usando (55) calculamos ｾ＠ P," como:
p ｾＢＺ＠ h el _
S_ cos a casi Á. - - - - - 550
2 2 I Mm 1I 't"l
Entonces para el triángulo PI ｐｉｾＧ＠ podemos escribir (consideran-

do un triángulo plano).
ｾ＠ I:>h _ Po e"
...:z.....L ______ 56
sen(a z,-1801 s
32
La cual finalmente nos da, después de algunos arreglos, la fórmu-
la final para la corrección de la normal sesgada.
ｾｉ＠ _ "( h2 2
on - p -
Mm-e cos
sen O1 2 O cos 2 A.
.12 't'2
1 --- --- 57
Cuando 4>z = 45° Y h2 = 200m y 1000m. 8h, es igual a 0.00 s" y 0.05,
respectivamente. Cl:>viiirilente habrá casos aonoe el efecto sea significati-
vo y deba ser tomado en cuenta. Esto es particularmente cierto en traba-
jos de cálculo de posiciones geodésicas de orden ｳｵｰ･ｲｩｯｄｾ＠
El segundo efecto geométrico a considerar en la reducción de me

diciones de direcci6n es aquel de la diferencia entre la secci6n normal, an-
la cual tenemos ahora reducidas nuestras mediciones, y la geodésica. Es-
ta corrección, la cual se deriva simplemente por la combinaci6n de las e _
cuaci6nes (30) y (48), se expresa como:
Cuando c/>m=O"; 0 12 =45° y S= 200 Km, 100 Km, 50Km

8g es 0.12"; 0.02", Y . 0.006" • Este efecto podría ser signi-
ficativo y debería tomarse en cuenta para trabajos geOdésicos •
.Algunos puntos finales de acuerdo a estos efectos geométricos son

descritos inmediatamente abajo.
1 • - En la ecuaci6n (57) la altura elipsoidal h puede ser sustitui-

da por la altura Ortométrica H sin un efecto significativo sobre 8h.
2 • - En casos extremos 8h y 8g seránde aproximadamente --

igual magnitud y de signo contrario. Ellas deberán· calcularse, sin
embargo, particularmente para ｣ｾｵｬｯｳ＠ de posici6n geodéSica ｰｲｾＺ＠
cisa.
3 • - Las ecuaciones (57) y (58) a menudo son expresadas en otra

formal' dando resultados Eiquivalentes,pero en las cuales pueden
incluirse aproximaciones posteriores. Como un ejemplo, la (57)
puede expresarse como (Bomford, 1971, No. 122) •
LI_8h_=_h_2_e_,2_sen
__ ＲＰ｟Ｑ［ｾｏｓ｣Ｏ＾ｭ＠ __ -,I _______ 59
y la (58) como (Bomford, 1971, p 124.).
I SKm 2sen 20
8g =+0.02S (\'00) cos 2 c/>m l ',
-----+f;O
33
FIGURA 17
CORRECCION POR DESVI.ACION DE LA VERTICAL.
34
2.1.2.- "Efectos Gravimétricos".
Un teodolito se nivela respecto al vector de gravedad local y no

a la normal elipsoidal. Una correcci6n para el ángulo (desviaci6n de la -
vertical) entre el vector de gravedad y la normal elipsoidal, es necesaria.
La Fig. 17 describe la correcci6n que debe ser aplicada. Este punto es -
tratado con profundidad en Vanicek, 1972, p. p. 164 - 166 • Nosotros
solamente establecemos aquí la f6rmula reducida como:
88 =-8cot ｚ］ＭＨｾ＠ 1
sen a-'TJ
It 1
cosa 12)cotZ ------ 61
Donde ｾ＠ es la componente meridiana de la desviaci6n de la ver

tical "l T¡, es la componente del primer vertical de la desviaci6n de la ve;
tical,:Z es la distancia cenital. El efecto de ésta reducci6n puede va :-
riar desde una cantidad insignificante (si 8!: O O si Z = O) a valo-
res de magnitud 2" - 3 11 cuando por ejemplo 8 =20 11 Y Z =ＸＰｾ＠
Para aplicar ésta correcci6n y la requerida en 2. 2, son reque-

ridas las desviaciones de la vertical en dada punto. Estas pueden obtener-
se de varias formas '(Vanicek ＱＹＷＬｾｖ｡ｮｩ￩･ｫ＠ 1972), pero dos métodos -
son usados comunmente llamados:, '
1. - Teniendo las coordenadas astron6micas ( ｾＮ＠ A ) en cada e,!

taci6n. lo cual puede ser una taréa dificil.
2. - Usando una contemporánea técnica de cálculo del geoide

(Vanicek y Merry 1-973) Y cálculando ｾ＠ y 'TJ en cada punto.
2.2. - Distancias Cenitales.-
El único efecto sobre una medici6n de distancia cenital está so-

bre las variaciones en el campo de gravedad. Esto es, las desviaciones de
la vertical. Como en 2. 1. 3., estableceremos solamente aquí la f6rmula -
reducida como:
ZR= Zm + Ｈｾｉ＠ cos a ll + 'TJI sen a lz ) - - - - - - - - 62
Donde Zm 'es el valor medido de la distanCia cenital.
Este tema es discutido en (Vanicek 1972, p. 170 Y Heiskanen -

and Moritz, 1967 p. 173 - 175) Y no será discutido aquí posteriormente.
35
2.3. - Distancias Espaciales.
En ésta secci6n trataremos la reducci6n de una distancia espacial

medida sobre la. superficie de la tierra a la superficie del elipsoide. Des-
pues de tener hechas varias correcciones instrumentales y atmosféricas a -
la distancia medida electrónicamente nos quedamos con una distancia eapa -
C
cial en línea recta "1" (Fig. 18). Esta distancia espacial es entonces redu-
cida al elipsoide. La reducci6n.es lograda cómo sigue :
Primero se calcula:
RI +Rt
= ｾ＠
- - - - - 63
2
Donde R¡ y ｾＮ＠ son los radios de curvatura de Euler ( eco 28).
Luego, del triángulo PI L 1>' 2, O Y por la ley de los cosenos tenema3:
Donde
Las cuales son las alturas elipsOidales y son iguales a la suma de

sus respectivas alturas ortogmétricas ( HI Y Hz ) Y sus alturas geoida-
les ( NI Y Nz ).
ｓｾｴｵｹ･ｮ､ｯＮ＠ｾ＠
cos ｾ＠ =1- sen! -2- ------- - 66
en (64) y reordenando los términos tenemos:
o
FIGURA 18
REDUCCION DE LA DISTANCIA ESPACIAL.
36
ｾ
Ｒ＠ 2 2 hl h2 2 \ji
= (h 2 - h) + 4 R (1 +-
1 R ,) ( 1+ - R ) sen - 2 - - - - - - - 67
Del triángulo P¡,P2'10' la ley de los cosenos y la fórmula de un

simi-ángulo tenemos:
ｾｯ＠ = 2R sen i 2 ---- -- 68
o ,1,
't'
ｾｏ＠
= 2 sen-1 ../;!:!.---------68a
2R
poniendo
(67) se convierte:
"
02 =uh
"2
+ (Ih l h2
+--}{I+--)
02
"0--------70
R R
y rearreglando:
ｾ＠
0=
[ t - b.h JiLJ ｾ＠
(1+:hL) (1+
2
2 --- -- - - 71
R R
Ahora:
I So: RI/I= 2R señ
l
i%- I ------ ----72
Entonces usando (71) y (72) podemos reducir una distancia es -

pacial a la superficie del elipsoide. Estas fórmulas son suficientemente
rigurosas para trabajos ｮｯｲｾ･ｳ＠ de geodesia (Thomson and Vanicek 1973).
Note qte para una reducción rigurosa de distancia la altura geoi-

dal N. es necesaria. Hay varios métodos de calcular N, uno de los cuales
está desarrolando por U.N.B. (Vanicek and Merry 1973).
No sé ha hecho mención aquí de las consideraciones de las lineas

base precisas. La razón de esta omisión es que las líneas base precisas no
están siendo medidas por nadie mas excepto para instrumentos de calibra -
ción. EDM. para los cuales no es necesario la reducción al elipsoide.
Finalmente, puede notarse que hay muchas fórmulas de reducción

de distancias en uso. algunas de las cuales han sido desarrolladas para elip-
soides de referencia especificos o regiones de paises.
2.4 • - Reducción al Terreno de Cantidades Geodésicas Calculadas.
La situación ocurre frecuentemente en la práctica donde las can-

tidades geodésicas calculadas, llamadas distancias y ángulos, pueden ser --
medidas sobre el terreno. Estas no pueden ser comparadas generalmente -
directament e con los valores calculados, ya que estos usualmente son
dados sobre la superficie del elipsoide de referencia, entonces ellos pueden
ser" reducidos'" al terreno.
Para reducir los ángulos requeridos uno procede como sigue:
37
Primero, calcule las direcciones (azimutes) entre los puntos involucrado,S
Luego usando las ecuaciones (57) (58) Y (61) calcular las cantidades Ｘｾ＠ ;
8é
80 y respectivamente.
Estas correcciones son entonces aplicadas a la dirección calcula-
da ai¡ ; con "signos opuestos" a los que fueran usados para la r¡duc -
ci6n al elipsoide para obtener la dirección que deberá ser medida ｡ｾ･＠
Obviamente uno no sería capaz de medir esta dirección (ó ángulo) exacta-
mente y la medida tomada tendrá algunas desviaciones normales. Un proce-
dimiento similar es usado para reducción de distancia. Un simple reorde -
namiento de términos en la ecuación (72) nos dá.
ｾｯ＠ = 2R sen ｾ＠ 2R -- ----- 720

Y similarmente (71) nos da:
2 (hl h2 ｾ＠
ｾ］＠ [ ｾｏ＠ｉＫｆｲｾｨＺＭｬ｝ＷＱＰ＠
Entonces podemos calcular la distancia espacial sobre el terreno

( \ ) dando la distancia elipsoidal So Rapidamente otra vez, como con la
dirección, podría ser.
38
SECCION III • - ｃａｌｐｾＢｏ＠ DE POSICIONES GEODESICAS SOBRE
EL ELIPSOIDE DE REFERENCIA.
3 • - Fórmula de Puissant - Líneas Cortas. -
3.1.- Introducción.
Estas fórmulas son llamadas así en honor del matematico, Fran-

cés a quien se acreditó su desarrollo. Su derivación está basada sobre una
aprmd.mación esférica. Estas fórmulas generalmente están consideradas
con precisión de 1 ppm. en 100 km. más aUa de lo cual ellas quiebran
hacia abajo rapidamente (40 p. p. m. en 250 kms. cuando 4> = 60°) (Bom-
ford, 1971, p. 134) Por lo tanto, decimos que la fórmula de Puissant es
una fórmula de linea "corta".
3.2 • - Problema Directo.
Dadas las cantidades geodésicas <p" A" S'2 y a'2 (Fig. 19), ｱｵｾ＠
remos calcular las cantidades 4>z A2 Y a 2, En la derivación, cal-
culamos primero <P2' Obtenemos, por la aproximación esférica, de la --
trigonometría esférica (ley de los cosenos):
sen A..
't'2 =sen 't',
A.. cos (PP
, 2) + cos 't',
A.. sen (PP)
, 2 cos a --- - - - 73
- .. S,z .J,..J,. .J,.
dado que se
P ero ｾ＠ Pz = -W;Y 't'z='t',+d't'; ya =a'2
estipuló que los fueridianos están en el mismo plano.
Entonces:
A.. 1 S'2
sen .J,.
('t', +d't' = sen .J,.
't" cos N:"" 't', sen S'2
+cos .J,. Ñ;-' cos a,z - ---- - 74
Se requiere ahora conseguir una expresión para d<p De la ecua
ción (74) podemos expresar la parte del lado izquierdo por:
sen(<p, +d4» = sen<p, CO$ d<p +
+ cos cP,
sen d<p --·--75

ELIPSOIDE
P¡ (</>, A,) F 1G U R A 19
FORMULA DE PUISSANT PARA PROBLEMA DIRECTO.
39
La fórmula anterior obviamente no nos dará la solución requerida
puesto que dep aparece en el lado derecho de la ecuación. Para comenzar
a solucionar este problema, usamos otra vez la aproximación esférica y ob
tendremos: -
SIZ
dep;" - - - cos all ----- ---- 81
NI
Sustituyendo (82) en (81) nos da:

. S'I s¡zl sfa ｓｾｉ＠ d.l. l
dep: -cos all-::::-:r tan4>t- -.- cosa12 +-- - cosZo.1tan epl+ ｾ＠ +--- 82
NI 2NI 6N I 2N,z 6
Del (82) anterior podemos obtener ahora una aproximaci6n más -
precisa para dep (despreciando los términos mayores que la segunda poten-
cia nombrada) •
. Sil ｾｉ＠ｓ 1
d.l.: -
ｾ＠ NI
cos a ll - - - tan.l. l (I-cos al)
2Nf ｾ＠ ,
+ --,.------ 83
La cual puede simplificarse como:
､Ｎｬ］ＡｌｃｏｓｴＭｾ｡ｮＧ﾿ｳ･＠ｾ＠ｾ＠ '1 Ｒｎｾ＠ｾ＠

a ｾ＠ +---------84
Elevando al cuadrado (84) y despreeialdo los términos mayores que la ter-
cera potencia nos da: .
S'
64"= ｾ｣Ｚｯｳｉ｡＠ｎｾ＠ 1
It
- ｾ＠
ｎｾ＠
COSClttsentalltan ep,+-------85
Y posteriormente.
dep':1: ｾ＠ｎｾ＠
cos'cl
It
+-------- 86
Finalmente sustituyendo (85) y (86) en (80) y reordenando los térmi -
nos nos da;
S SI S' S2
dep :1: .....!!.. cos a12 - -!!.a tan epi - ｾ＠ cos a12+ ｾ＠ cosa alz tan 4>1 -
NI
, 2NI 6N
,I 2N,

- ｾ＠ cosa senza tan!.I. + ｾ＠ cosla +---- --- 87
2Nr 12 I1 ｾｉ＠ｎｾ＠ 12
ReSUmiendo ｴｾｉＧｭｩｮｯｳ＠ nos da:
- ＠ｾ s' cos tL sen Z a - - -- - - - - - 88

Ｖｎｾ＠ -.z 12
40
Simplificaciones Il,layores son realizadas para obtener:
La cual cuando es colocada en (88) finalmente nos da:

ｾ＠ 5,2 ｓｾ＠ｓｾＲ＠
a,z sen2 ala
2
dep - ｾ＠ cos a,! - ｾ＠ tan epi sen a,¡ Ｖｎｾ＠ cos (1 +:3 tan! ep,) - - - - 90
La ecuación ( 90 ) no es una solución rigurosa puesto, que el-

radio de curvatura a través de la sección normal de Ff a P2 se esta
tomandó como un valor constante NI ,cuando de hecho, cambia cm
lalatitud (N=t.<epl y M=fe{ep). (ecuaciones (15) Y (13) re!.
ｰ･｣ｴｩｶ｡ｭｮｾＩＮｐｲ＠ tomar en cuenta este cambio en curvatura podemos
escribir.
dep= M1m [Lada derecho de 90] --------91
Donde:
M.+Ma
Mm:: - - 2 - -----------------92
En tanto no conozcamos epz debemos usar la aproximación.
En orden de cálcular Mz De ( 13 ) calculamos:.
A.. -en
dM
- - '- :: a (l-e Z ) (-:3/2)(I-e 2 senz't', A..
) (- 2e2 sen 't") A..
cos 't"- ---- - -- 94
dep
41
\ _ d'; SIZ ｓｾ＠
d 1\ + --- -_ ( - - -- - __ o) (sen O sec<,/>--- 103
6 6N: 6N: 12 2
_ SI2 ｓｾｺ＠ d';
ó dA- - - sen a l2" sec4>z- - - sen a sec 4> + - - - --- 1030
Nz Ｖｎｾ＠ 12 z 6
Ahora, de los dos primeros términos de (103 a) (despreciando los tér-

minos mayores a la tercera potencia).
d ｾ＠ = Ｖｎｾ＠Srz sen s a secs 4> + - - - - - - - - - - 104

IZ 2
Lo que nos da:
dA'=p'l ｾｳ･ｮｑ＠｛ Nz IZ
secA-. {1-L(I-sen2 Q sec 2 A-.)}]--------105
't'z Ｖｎｾ＠ IZ 't'2
La cual cuando es colocada en ( 101 ) da la solución para >"2 •
Aunque a ZI es tamb ién una parte del problema directo, la deri -

vación para su solución se da e.n la sección siguiente.
3.3. - Problema Inverso.
Aquí estamos dando las cantidades 4>1. >"1 de y 4>2' >"z de P2 r::,
(Fig. 20). Las cantidades requeridas son SI2. a IZ ya 21 comenzamos-
por determinar a21 • Usando una aproximación esférica.
P' ｾ＠ｾ＠ = 360° a 21 ------ - - 106

y
i2( P'plPZ +i: p'p p)=_I_(a +360-a )---------107
2 1 2 IZ 21

Donde da es el término que expresa la convergencia de los meridia-
nos entre los puntos P1 Y P2' Usando la figura 20, podemos escribir:
a 21 : 0;2 + 180" --------- 109
Y sustituyendo CJ{t Pllt .(líOS ＧｬＡＦｾ＠ｾｦＴ＠
Luego sustituyendo-a zl en (107) por (109 a):
42
La cual se reduce a;
2
ｾ＠ = M 3e sen p,cosp,
- - - - - - - 95
dep '(I-e 2 sen 2 .p,)
La cual al sustituirla en (92) nos da:
Mm = ｟ｍｾＬＫＨＭｌｉ＠ +.;. .d.;:.:M.--'IL.) = M + _d_M_,_ - - - - _____ 96

2 '2
dM d"+'"
Mm= M, + --'- ＨｾＩ＠ --- - --- 960
dep, 2p"
2
Mm=M,+2..M, e senp' COSPI ＨｾＩＭＹＷ＠
2 (I-ésen ep¡) pI'
De(97): _1_= _1_ [1_2- e25ef1cP,COSPI dcP" 1] -------98
Mm M, 2 (1- eZsen'-<f> ) p"
La cual cuando la colocamos en (91) nos da. el resultado final:
dep"=[ p"( S,ZCOS a,z ｓｾ＠ ton epi sen! a,z

MI 2 MI NI
Donde d4> 11 en el último término de (99) es calculado usando la e -

cuación (90) (multiplicada por p").
Finalmente, calculamos epI por:

tpz =tpl + d tp - - - - - - - - - - 100
La longitud de P2 puede calcularse de :
ｾＲ］＠ ).,+ ､ｾ＠ --- --------- --101
De la Fig. 19 usando una aproximación esférica, . la ley de los senos nos

da:
sen d ｾ＠
S'2
sen ｾ＠
o
sen d>.. = sen.§L sen al! sec epz- - - - - - - - 1020
Nz
Ahora, aproximando los términos senoidales sobre cada lado de --
(102 a) por una serie de Taylor, podemos escribir (despreciando los térmi-
nos mayores a la tercera potencia).
43
FIGURA 20
FORMULA DE PUISSANT PARA PROBLEMA INVERSO
_1- ＨｾｐＧ＠ + ｾ＠ P'F! P. ) = 90 _ ｾ＠

0
P. P- ---- --- 1100
2' I 2 2 I 2
Usando trigonometría esférica, la ley de las tangentes nos da:
tan (90
0_ ｾ＠ _
) - cot 2
A cos + [(90 °-epz) _(90 _epl)]
I
0 0
o --- - - - - 111
2 cosT [(90-ep2)+(90-ep,l]
La cual se reduce a [Invirtiendo ambos lados (111 )]
da cos (90 - - 2 -
0 ｣ｦｲＧｾＲ＠
dA
)
tan -2-= - - 1 - - = - - - - - tan -2----- ---- 112
o cos 2" (epi - epl)
I
tan da = Sen "2 (pi + Pz) ton ｾ＠ --- - - - - - 1120
2 dep 2

cos --
Enseguida ､･ｾ｡ｲｯｬｭｳ＠ los términos tangentes sobre ambos la-
dos de (112 a) que pueden ser expresados por (despreciando los térininos
mayores a la tercera potencia).
da
tan -=sen epm sec ｾＨＭ､ｾ＠ dA
+､ｾ＠ - - + ----)------113
2 2 2 24
y
ｴｯｮｾ］＠ -ºº-+ das + ------114
2 2 24
La cual nos da la ecuación final.
44
d a " = p,,[ dAsen't'm
A..
sec dA..
ｾ＠ + -dXS (sen't'm
A..
sec ｾ＠ d A.. - sen3q,m sec'( ｾ＠ dA..
) + ] ---- --115
<
2 12 2 2
Dondeq,mes la latitud media.
Sustituyendo da en (109 a) por (115) nos da la requerida 021 una vez

que tenemos una expresión para al2 • La solución para al2 es como si- -
gue. Tomando la ecuación (99) y reordenando términos, obtenemos:
y usando (105), un rearreglo de terminas da:
Ｖｎｾ＠
Ahora, dividiendo (117) entre (116) da, después de algunas manipulaciones -

de términos.
117
tan a 12= -1-16- --------- 118
Puesto que a l2 aparece sobre el lado derecho de (118) es necesa -

ria la iteración. Primero comenzando por obtener valores aproximados pa-
ra a l2 de (118) usando solamente el primer término en el numerador y de -
nominador y para Sl2 de (116) o (117) usando nuevamente solo el primer
término del lado derecho de las ecuaciónes. Valores mas precisos de a l2 Y
SI2 son obtenidos usando todos los términos en (118) y (116) o (117) res-
pectivamente. Se itera hasta que las variaciones en a l2 Y SI2 sean des -
preciables.
3.4. - :Resumen de Ecuaciones para lo solucion de los problemas

Directos e Inversos usando la fórmula de Puissant.
El siguiente, es un bosquejo de las etapas requeridas para la solu -

ción del problema directo usando la fórmula de Puissant.
1. - Calcule MI y N l usando (13) y (15) respectivamente.

2. - Calcule una dq," aproximada usando (90).
3. - Resuelva para dq," usando (99) y q,2 usando (100).
4. - Calcular N2 con (15).
5. - Resolver para dA" con (105) y A2 usando (101).
6. - Usando (115), calcular da" y finalmente a 21 con (109 a).
Similarmente podemos delinear las etapas requeridas para la solu -

ción del problema inverso como sigue:
45
1 . - Calcule MI COil (13) y ｾ｜｝Ｑ＠ Y N 2 usa1do (15);
2 . - Calcule al2 con (118).
3 . - Calcule da" C01 (115), luego aZI usaado (109 a).
4 . - Usando cualquiera (116) o (117) calcular S'Z •
3.5. - La Fórmula de la Latitud Media de Gauss.
Estas formulas se publicaron primero en Inglés en 1861. Están ba-

sadas sobre una aproximación esférica de la tierra y solamente podrá usarse
para puntos separados por me 10S de 40 kms. y con latitudes menores de 80 o
(Allan et al, 1968) • Las fórmulas son (Allan et al, 1968) ,
dQII:d),."sen4>m------ 119
d<f¡": ( S'ZCOS °m 1.p"-----120

Mm
d>..l1:p"( S'2sen a m )-------121

Nm cos<f¡m
Las similitudes de las fórmulas anteriores con las fórmulas de ｐｵｩｾ＠

ssant son fácilmente vistas comparando (119), (120) Y (121) con los ｰｲｩｭｾ＠
ros términos de (115), (99) Y (105) respectivamente.
Para resolver el problema directo con la fórmula de la latitud me -

dia puede ser usado un procedimiento iterativo. Primero, d<f¡" puede ser -
aproximado usando el azimut medido en lugar de am y MI puede ser usado
en lugar de Mm'
Luego una primera aproximación de <f¡2 se obtiene usando (100), en

tonces da se calcula vía (119). El proceso iterativo puede ser continuado
ahora usando valores aproximados sucesivos para d<f¡ ,da (entonces am
y <f¡m ) hasta que los límites deseados hayan sido alcanzados. Finalmente-
d ｾｉ＠ es calculada para obtener },,2

El problema inverso es calcular sin iteración ya que <f¡m es i.nmedia-
tamente aprovechable. Usando (119) da es calculado, luego de (121) divi-
dido por (120) uno obtiene, Tan Qm ,entonces a '2 Y a 21 finalmente la dis -
tancia S'2 puede calcularse con cualquiera de (120) O (121).
3.6.- Otras Fórmulas De Lín eas Cortas.
Hay muchas fórmulas de líneas cortas en uso. Algunas de ellas son

incluidas en (Bomford 1971 pp. 133, 139) Y son llamadas por nombres ta-
les como "Fórmula aproximáda de Clarke" (1pp. m en < 150 km) y--
"Fórmula aJlroxitnadlir' de Lilly" (15 ¡no en 1000 km). Todos estos tipos
de fórmulas directas o inversas (líneas cortas) están basadas sobre ｡ｰｲｯｸｩｭｾ＠
ciones esféricas y no son tan rigurosas como otras tales como la fórmula de
Bessel para linea larga, desarrollada en 4.
46
4. - Fórmulas de Bessel - Líneas Largas.
4.1.- Introducción.
Las fórmulas para los problemas geodésicos directos e inversos
desarrollados abajo han sido acreditados a Bessel (Jordon 1962). La deri ..
vación está basada en la Geodésica sobre el elipsoide. Este hecho distin -
gue las fórmulas de Bessel de las fórmulas que están basadas sobre una --
aproximación esférica (ej. la de Puissant) o también de las fórmulas basa-
das elipsoidalmente, pero que usan la curva de sección normal como funda-
mento para la derivación (ej. Robbins 1962). La exáctitud de las fórmulas
de Bessel no está limitada por la separación entre los dos puntos en cues -
tión 6 por la localización de los puntos sobre la tierra. La exáctitud está li-
mitada simplemente por el número de términos que uno desee retener en las
series desarrolladas de las diversas expresiones.
Las siguientes derivaciones comienzan por desarrollar las rela -
ciones entre los elementos correspondientes sobre la esfera y el elipsoide -
(no una aproximaci6n esférica, sino un tratamiento riguroso). La soluci6n de
una integral elíptica es entonces realizada. Finalmente son enunciados los -"'
problemas directos e inversos.

4.2. - Relaci6n Fundamental.
Comenzamos estableciendo algunas relaciones rigurosas entre los
parámetros sobre la esfera y los parámetros sobre el elipsoide. En la sec-
ción (1. 3. 2.) desarrollamos la propiedad básica de una geodésica (47) la que
47
sobre una esfera puede ser expresada como.
cos{3sena: cos {30--------- 122
Donde {3 es la latitud reducida (Krakiwsky arid Wells, 1971,
p.23) Y {30 es llamada la latitud reducida del "Punto de Giro" ( a-90°).
De la figura 21 a a sobre la esfera reducida es igual a a sobre el elip -
sOide, como están {3 en la esfera reducida y {3 en el elipsoide entonces
podemos escribir para ambas:

cos {3 sen a :: cos {30 - ---- - --- - 1220
Desarrollamos ahora alguna relación difer81cial con la ayuda de -
la Fig. 21b. De los triangulos en las figuras esféricas podemos escribir:
o do- cos a 12: od {3

--- -- ------- 123
y odo- sen al!: o cos{3' d A
Donde (a) es el radio de la esfera reducida (Fig. 22) Y (do- ) es -
el ángulo subtendido (en el origen de la esfera) por las normales en P y PI.
Similarmente de los triángulos de la figura elipsoidal podemos escribir:
ds cos a 12 : Md<p

y ds sen a 12 : N'cos <p'dl-- -- -- -- 124
Dividiendo.(124) entre (123) tenemos:
ｾＺ＠ Mdp : N' coscP'd' -:-------- 125

do- d{3 cos{3' dA
48
De la figura 22 y de la ecuación (17)
N'cosep'= o cos {3'- --- --- - 126
La cual al ser sustituida en (125) da:
ds dl
- - = 0 - - ------- 127
da- dA
POLO
ESFERA REDUCIDA
FIGURA 21a
RELACIONES FUNDAMENTALES PARA EL DESARROLLO

DE LAS FORMULAS DE BESSEL.
ｾ＠
p ｾ＠
f \
Q.
o cos (3'
'\ !
ｾ＠ .,. ｾ＠
-
FIGURA 22
ESFERA REDUCIDA Y ELIPSOIDE
49
FIGURA 21b.
RELACIONES FUNDAMENTALES PARA EL DESARROLLO DE

LAS FORMULAS DE BESSEL.
di , ds
o ｾ＠ =-0- d'<T ------1270
En la cual. de (125) tenemos:
-2.L =.!!... ｾＭｾ＠ -- - -- 127b

dX. a d{3
Recordando que (Krakiwsky and Wells 1971. p. 28) •
tan {3= ＨｉＭ･ＲＩｾ＠ tan</> ------- - - 128
Podemos diferenciar y obtener:
ｾ］Ｈｉ｟･ＲＩＱＯ＠ｾＭ 129
oos2{3 cos2</>
o

...!!5L = - ' - ｾＭ 1290
d {3 ＨｉＮＬ･ＲＩＱｾ＠ cos 2 {3
La cual nos da. cuando la sustituimos en (127b):
ｾ］＠ M ｾＭＧＳＰ＠
d X. a (l-e 2 )lk cos 2 {3
Para cualquier punto sobre el elipsoide.
Ahora queremos obtener:
2L= f ({3)
Comenzamos por expresar

dX.
a cos {3 = + cos </> ---- 131
2 2 -1'2
donde V = (I-e cos {3) - - - - - 132
50
y (la curvatura en el polo - ecuación (5 a».
ｃＺｾ＠ - - - - - - 133
b
Elevando al cuadrado (123) Y rearreglando los términos nos da:
d! o
---:-- - - --- -- 134
dA Ve (l-e 2 ) ｾＲ＠
donde
M : _C_ - - - - - - -- 135
V
Una reducción mayor de (134) usando (13.3), (3) y 131, finalmente da:
ｾＺ｟ｉＭｾＭＱＳＶ＠
dA VadO"
Antes de procesos posteriores, derivaremos (132) de (131).
coscp:(-{-)Veos/3------ 137
La que cuando la elevamos al cuadrado da: eos 2 cp = ｾｶＲ･ｯｳ＠ 2
/3 ---- 1370
0
o
Sustituyendo (137b) en (137)

2
V : [1 +e'2(I-ll yleas
2
/3]----- - 138
Lo que se reduce a:
V2 [1_eI2(1_ e2 ) cas2 /3] =1 -- -- - -- 138 o
Ahora de las ecuaciones (3) y (4)
2
(l-e,l) (l+e )= 1- - - - - - - 139
y
La que al ser sustituida en (138a) da:

I
V (I-e
2
col /3) = 1-------- 140
Retornando a nuestro problema entre manos,

Sustituyendo (140 a) en (136) obtenemos:
ｾ＠ :(I-e2cos2 /3)Y!--- - - 141

dA
y
ｾＺ＠ o (I-e2 COSE /31/1 -- - - - 142
dO"
respectivamente. -
51
4.3. - Solución de la Integral ElÍptica.
Enseguida resolvemos (141) y (142) Y 10 hacemos así como por

integración. Comenzamos por resolver (142) para obtener una solución -
para ds/ dcr. De la Fig. 23, usamos la ley de los senos de la trigonome -
tría esférica y obtenemos.
sen O,Z = sen '300

- - - - - 143
sen (soo-f301
o
cos/30 = sen altcos f31 - - - - - - - - - 1430
La propiedad fundamental de una geodésica y de un círculo prin
cipal. Posteriormente usando la regla de Napiers de partes circulares:-
FIGURA 23
SOL UCION DE ds

dT
cos Cl¡z = coto; ton f31 -- - - - - - 144
o tono; = tonB - - - - - - - - 1440

cos a lz
y otra relación requerida.

sen f32 =sen (0;+ a;. 1 senf30 --- - -- -- 145
Generalizamos (145) por propósitos de integración (entre los puntos PI Y P2
Fig. 23) como:
52
sen /3 = sen (o; ter) sen /30 -- - - ---- - 1450
De modo que si er es variable. medida desde PI- Note que

cualdo er= er2 ; /3= /32 y cua ldo er= o. ＯＳｾ＠ /3,
Describieldo (142) como:
ds=o ( I-e2 cos 2/3)'/Z der---------- 146
y resol vie ldo eltonces para cosz /3 de (145a) por:

cos 2/3= l-senZ{o;+er) senz /30 ---- - - - 147
en la que sustituimos O; = er y X=a; ter (una nueva variable para integra-
ción) e'ltolces dl(= der y reescribimos (147) como:
cosz /3= l-senZ Xsen l /30 ---------1470
La que finalmente da:
ds = o (l-e z +e 2senz /30 senz X ＩｾＲ＠ dx ---- 148
de (3) y (4)
2
e' y '_el = - - '/2- - ---- 149
I+e'! '_e
Los que cuando se sustituyen en (148) da:
ds =0 [ - - '-
l-e /2
+ ｾ＠ l+e,2
sen! /30 sen! xl'J.z dI( ------1480
o
ds=
Puesto que _b_= - - - - - - 150

Q
y obteniendo: k2=e,2 sen 2/30 - - - - - - - --- 15/
I
(149) se convierte finalmente: ds= ｢ＨｬＫｫＲｳ･ｮｘｾ､ｸ＠
------- 152
Esta expresión es integrada ahora y evaluada para nuestros parámetros

particulares. los que:
1( ］｡［Ｋｾ＠
x
S::b
f
=a;
(I+k sen Xl
2 2 \re
dI( - - - - - - - 153
En matemáticas ésta es conocida como una integral elíptica (<\bramowitz and

ｾ･ｧ￺ｮＬ＠ 1968, p_ 589)
Los límites sobre X {er+erTl son: ｏｾ･ｲｔ＠ -- - - 154
Entonces cuando: er =0; X= O; - ------- 1540
53
y cuandoCT =CTri X=o¡ tCT -------154b
r
Resolviendo la ecuación (153), sabemos que debido a que k 2 es pequefia
entonces. I
z 1 1 k6
lI+ k2sen 2 X1 2 = 1+ - k2 sen 2 X - - k 4 sen 4 X t - sen6 X ------ 155
2 8 ｾ＠
Usando las identidades trigonométricas:

s8n 2 X = .L ( 1- cos 2 Xl· sen4 X =- - - - etc - - - - - 156
2 '
Y sustituyendo en (155) da:
(ltk2sen2Xl\l2= [Ii-Jé.
4
- ｾＭ 64 k4 +----]+[- _,-
4
k2 t _,- k4
16
+ -----Jcos X-
. k4
- 64 cos4X+-------155a
Recolocando :
2
A= I +k- - - -3k ' t4 - - - - - - 157
4 64
B= _1- k2 - _1- k4 t - - - - 157 a

4 16
c= ｾ＠ t-----------157b
64
0:------- ---------157c
En (153) da:
s
-b- =A
fO¡ +CTdx-B
r fa¡+CT
cos 2Xdx-C
¡o¡tCTr
T
cos 4Xdx---------158
o¡ a¡ o¡
Antes de obtener la integración actual de (157) consideraremos la solución
de la integral general.

o¡
CTr t o: +0:
la¡ cos n xdx = -1. sen nx I
na¡, [
T ----- 159
= ｾ＠
]
sen (o¡ tCTr , - sen n o¡ ---- - - -- 159a
Otra sustitución da una forma mejor, nombrada:
sen nx - sen ny = 2cos1!...l X +v) sen JL \X-V' ------160
2 2
La que cuando las asociamos con nuestro problema tenemos:
X=o;+CTr , V=o¡ -------161

entonces
xtV = 2o¡ tCT -- ---- --- 161a
r
y
X-V= CT
r
Ahora, en (159a) el lado derecho se convierte:
54
Esta representa la integración de la distancia sobre el elipsoide
respecto a la distancia sobre la esfera.
Ahora regresamos nuestra atención a la solucion de

(141) . Reescribiendo (141) obtenemos:
､ｾ＠ =(I-e2 cos 2 (3) '1.2 dA - - - --- 1410
De la Fig. 24.
dA cos (3 = dO" sen 0'2 - - - - - -- 167

o
dA= sen 0'2 . dO- ------1670
cos (3
Aplicando la ley de los senos (Trigonometría Esférica).
0
sen ct '2 sen 90
= ____ ----168
ｳ･ｮＨＹＰｾｦＳＩ＠ sen (90°-(3)
o
cos (30
sen a,¡ - - - - - - - - 1680
cos (3
La que cuando es sustituida en (167a) da:
dA= cos {30 .,,0"-------- 167b

cos2 (3
Sustituyendo para dA en (141a) obtenemos:
R 2 2 (3 '/2 COS (30
d,l. =O-e cos ) - -2 _ . dO" ---- - - 169
cos (3
Enseguida tomamos ､ｾ＠ menos (16Th) lo que nos da:
- - - - ] dO" - - - --- 170

l¡, cos2 (3
Desarrollando ( l-e 2 cos 2 (3 ) 2 en una serie nos da:
2
= 1- -2e cos2 (3 - -8e- cos 4 (3 - -16
e
'1. 4 6
- cos· f3 ------171
(I-e 2 cos 2 (3)2
55
senn(O"+a: l-sennO" ］Ｒ｣ｯｳｾＨＰＢＫ｡Ｚ＠ lsen..!!-a: -------162
IT I 2 IT 2T
Ahora evaluando (158) obtenemos.

rO; +CTT
Ja: dx = CTT- - - - - - - - 163
I
rOj+CTT
ｊｾ＠ ces 2Xdx=cosl2a:+a:lsena: ---------- 1630
o¡ I T T
¡ o¡
a;+O"T 1
cos 4Xdx = -2 cos l4Oj+2o;. 1 sen 20;. - - - - 163b
etc.
Poniendo.
o;. =OZ-O¡---------164
Entonces:
20; +0;-= 20;+CT2-a¡ --------- 1640

o
y
o¡ +CT2 -----------164c
2
o
2CTm = 2a¡ +CTT - - - - - - - - - - 164d
Cuando sustituimos en (163 la solución para (158) es:
-Sb =Aa:T - B cosO" sena: - -

e o
cos 4CT sen 2a: - - cos 6CT sen 3a:T -------- 165
m T2 m T3 m
De (164) obtenemos una solución para CTT como:

S 8 e
o;. =Ab +A ces 2CTm sen CTT +"2A cos 4CTm sen 20;. ------ - - - 166
Donde: k2 3
4
A = 1 + - - - - - k + -----
+
4 64 •
B = -+- k
2
- k
4
+---------.
_1_ k4 + _____ _
C = 64 '
D=
5
E = 65536
56
FIGURA 24
SOLUCION DE ｾ＠
dA
La que cuando la dividimos por cos 2 f3 da:

1 el e 4 e 6
1 = -2- - - - - coszf3 - - cos 4 f3 - - - - - - - 171 a
cosf3 2 8 16
La ecuación (170) es ahora:
4
d,(=dA- cos f30 [ -2e2 + - e8 cos 2
eS
f3 + -16 cos· f3+---- - ] do- - --- 172
o
e22
[e4
dl=dA--cosf3o 1+-cos2 f3+-cos 4f3+---- ] do------172a
2
e8
4
Para la solución de (172a) sustituimos cos 2 f3 , cos· f3 , etc.

por cos 2
f3 = 1- sen f30 2
sen! X- - - - - - - 173
y
cos 4 f3 = 1- 2 sen! f30 sen 2 X + sen 4 f30 sen4 X -- - - - - - 1730
(X está definida en la pago 43), los que cuando son sustituidos en (172a) nos
da:
La expresión anteriDr es simplificada y adecuada para integración

en mucho de la misma manera como fué hecho para la solución de ds / do-.
Los resultados son como sigue. La diferencia de longitud sobre el elipsoi
de , está dada por: -
Lz
L= f d.l- ------- 175
L,
57
y sobre la esfera por:
A= ¡ A2 dA ________ _
1750
A,
Entonces:
L= A- 1-2 cos {30 [ Jo:ra¡, +o-T(A+B'cos2X +C' cos4X+----)dx ] -------176

Donde:
4
, e2
B=-sen 2 e
{30+-sen 2 e4
{3o---sen4
130-------1770
8 8 16
e4
I
C =-
64
sen 4 130 + - - - - - - 177 b
/
D =---------- ---
El resultado es dado entonces por:
L= A-.f..
2 T T
+ ｾ＠
cos /30 [AO: + B'seno: cos20:m2 +
sen 20:T cos40:m senJL
3'
30- COS 60: + -- --
m-
-J ---- -178
Ahora con todas las relaciones necesarias desarrolladas, retornamos nues

tra atenci6n a los problemas directos e inversos.
4.4. - Problema Directo.-
Recordamos que para el problema directo debemos conocer las

coordenadas geoáésicas ep, ,A, del primer punto P 1 Y la distancia geo-
désica S'2 yel azimut a'2 al otro punto P2' entonces resolvemos -
para ep2' A. de ｾ＠ °
y 2 , • Las etapas en la soluci6n son como sigue:
1. - Calcule la latitud reducida 13, usando (128).

2. - Calcule el Azimut de la geodésica en el ecuador que es:
sen a = sen a,. cos 13,--- --- --_. 1220
58
3. - ｃ｡ｬ｣ｾ･＠ el a;-co. esférico ｡ｾｲｯｸｩｭ､＠ 00 de (166) usando solamente
el pnmer ternnno, (e. g. 00 =--¡)A ) luego· calcule 0]+1 por:
o:1+1 =0-0+ -ª-cos20- sen m +-------

A m I
Donde la primera iteración CT¡ = 0-0 Y recordando que:
20-m = 2o¡ + o¡ -- ----

En la que 0- 1 es resuelta por (142a) ésta etapa se repite hasta decir
100n-0;1'"0.0000¡"
4. - Calcule f32 por (145) donde f30 es calculado usando (l43a).
5. - Calcule ep2 usando (128)
6. - Calcule la diferencia de longitud esférica A usando la ley de los senos
(Fig. 24), la que da:
sen A= sen o- sen a '2 _ _ _ __ _ _ _ 180
cos f3z
Luego usando A de (180). calcular a • cos 2 o-m. cos 4 <Tm ; cos 6 <Tm ;
Usando(184), (185), (185a) (185b) respectivamente; Usando (179) resolver pa
ra ( A-L ); esta etapa es entonces repetida, con L= A - ( A-L ) (186) has-=-
ta digamos. Il A-Ll-l A-LlI '::0.00001"; finalmente:
A2 = A, +L
7. - El azimut inverso es entonces calculado vía (186a) o (187a)

POLO
ECUADOR
FIGURA 25
SOLUCION DE LA LONGITUD DE ARCO cr:.
4.5.·- Problema Inverso. -
En éste problema estamos dando ｾ＠ ( epi A, ) y ｾ＠ ( epa A2 ) de

los cuales calculamos S,z ; alz Y a z, •
La primera etapa es calcular 73, y f3z (latitudes reducidas) usando
128. Luego de la esfera reducida (Fig. 25) podemos calcular la longitud
de arco ( o- =o-T ). usando la ley de los cosenos de la trigonometría esfé -
rica como.
59
cos O" : sen {3, sen {3z +cos {3, cos {3z cos A-- - - - - - 181
o sen 0": [( sen Acos {32 )2+ (sen {32 cos {3, - sen {3, cos {3z cos A lj - - ｾＭ ---- -- -- 1810
Ya que éste es un problema iterativo, (181) se resuelve primero usando

A :: L en la primera aproximación calculamos entonces:
sen A COS{32
sen a lz : ---------- 182
sen O"
Para calcular el azimut de la geodésica en el ecuador a, combinamos

(143a) y (182) 10 que nos da:
sen a lz cos{3, =senacosO o -·-------183

o sen a
sena: --------1830
'2 cos ft,
Lo que cuando es sustItuida en (182) da:
sen a :
cos {3, cos f3z sen A
-- - - - - 184
sen O"
Otra vez. sen a es solamente una primera aproximación ya que A'; L •
Luego calcule:
cos 2 o: : COS O" _ 2 sen {3, sen /3z --------185

m cos 2 a
cos 4am: 2 cos2 2 am- 1 - - - - - - - - - - - - - 1850
y cos 6 O" :o 4 cos3 2 O" - 3 cos 2 O" - - - - - - 185 b
m m m
Usamos entonces (179) para calcular ( A- L ). Después de completar esta

etapa, calculamos:
A=L+(A-L)--- - - -186

Y regresamos a (181) y recalculamos las cantidades CT,a, 20"m '
40"m , 60"m , usando (181). (184), (185). (185a) y (l85b) respectiva-
mente. Después de recalcular (A-U Usando (179) ensayamos 1< A-UI+,-
< A-L\.1S 0.00001". Cuando esta prueba pasa, i + I procedemos a ｣｡ｬｾ＠
lar a'2' a z, y 5'2 . • El azimut de adelante es calculado usando -
(l83a) que es reescrita aquí como:
sen a
- - -- - - - 186
cos {3,
y
sen a
cos {32
Alternativamente los azimutes pueden calcularse por:
60
sen >.. cos /32
ton Q'2= - - - - - - - - - - - - 187
sen ,8zCos ,8, - cos >.. sen ,8, cos,82
y
- - - - - - - - 187 o
Para completar el problema, la distancia S12 es calculada usando (165)
4.6. - "Otras F6rmulas para Lineas Largas".
Muchos métodos para la soluci6n de los problemas directos e in-

versos para puntos ampliamente separados sobre un elipsoide de referen
cia, están disponibles en la literatura. Como con las f6rmulas de líneas
"cortas" y "medias", ellas generalmente dan los nombres de sus descu
bridores. Dos de éstas, que han sido usadas por los autores, son los :-
métodos de Rainsfotd (Rainsford, 1955) y Sodano (Sodano 1963). Las
f6rmulas de Rainsford están desarrolladas sobre los mismos principios
que las de Bessel. La diferencia mayor es que los coeficientes de la di -
ferencia de longitud (179) están desarrollados en términos de (f) ya que -
･ｾｯｳ＠ convergen más rápidamente que cuando son dados en funci6n de -
(e. ). La diferencia principal entre el Método de Sodano y éstos de Be -
ssel y Rainsford es que tanto el problema directo como el inverso pueden
resolverse en una manera no iterativa.
SECCION IV. - CALCULO DE POSICIONES GEODESICAS EN

TRES DIMENSIONES.
La posici6n geodésica de un punto del terreno puede describirse -

matemáticamente en términos de una triada de coordenadas cartesianas
(X, Y, Z,) referidas a los sistemas de coordenadas terrestre promedio,
geodésico, geodésiCO local, o astron6mico local o por latitud geodésica
( 4> ), longitud ( A ) y altura elipsoidal ( h ) referidas a algún elipsoide de

referencia. En las secci6nes previas en las q'IE fuéron presentados los -
cálculos clásicos de posicionainiento bidimensional, las posiciones geodé
sicas fuéron descritas solamente por dos coordenadas, llamadas la lati :-
tud y la longitud geodesicas. La tercera componente. altura elipsoidal, -
fue usada solamente para la reducci6n de mediciones terrestres al eUpso.!.
de de referencia.
Los cálculos de posici6nes geodésicas en tres dimensiones difiere

del clási. co proceso bidimensional en dos formas significativas: La prime
ra es que éste Último se basa en la geometría elipsoidal mientras que eC
primero se basa sobre los principios Euclidianos Tridimensionales yem-
plea el algebra vectorial y matricial. La segunda es que el método clási-
co requiere el uso de distancias geodésicas y azimutes para cálculos rig.!!
61
rosos mientras que distancias espaciales en línea recta (cuerdas)y azimuts
tridimensionales de la sección normal son usados en cálculos tridimensio-
nales. Considerando el azimut aquí usado, deberá notar se que se refiere a
la sección normal que pasa a través de los puntos del terreno en cuestion -
y no de la sección que pasa a través de los puntos proyectados sobre del
elipsoide. En vista del tratamiento diferente de las observaciones en cálcu
lo de posición tridimensional, no se considera un capítulo especial para pre
sentarlas. En su lugar, son dadas explicaciónes a fondo, donde se requie :-
ran, dentro del contexto del desarrollo de los prcblemas directo, inverso, -
intersección azimutal e intersección en la distancia espacial.
5. - PROBLEMAS DIRECTO E INVERSO EN TRES DIMENSIONES.
5.1. - Problema directo.
El problema directo puede definirse como: Dadas las coordenadas

(Xi, Yi. Zi) o (4)., A. , h. ) de un punto (i) y la distancia espacial terres
tre, azimut y ánguÍo vértibal (o diferencia de altura) a un segundo punto :-
(j) calcular las coordenadas (Xj. Yj. Zj,) o ( 4>. , A. , h. ) • Dos casos -
J
del problema directo pueden surgir. dependienao de 4ueJ el azimut y el áng.!!
10 vertical sean referidos al sistema de coordenadas geodésico local (nor -
mal. elipsoidal) al astronómico local (vertical de la gravedad). Entonces-
denotando los azimutes y los ángulos verticales en el sistema geodésico lo-
cal por ( a ) y (a) y de la misma manera en el sistema astronómico local
por ( A ) Y ( V ) respectivamente. (Fig. 26).
El método mas simple de solución de problemas tridimensionales

es usar coordenadas cartesianas. Si las coordenadas que se requieren en
los cálculos son dadas por (4), A. h ) una simple transformación de coorde-
nadas (krakiwsky and Wells, 1971) da las coordenadas cartesianas. Simi-
larmente si el resultado requerido es el de latitud, longitud y altura elip -
soidal, entonces las coordenadas cartesianas son transformadas a (q>,A.,h).
después de que los cálculos de posición son completados (Krakiwsky and -

Wells, 1971).
El vector entre dos puntos terrestres en un sistema de coordena-

das geodésicas está dado por la expresión.
- - - - - - - - /88
62
Z6L
ZGL (NORMAL ELlPSOIDALl
VGL
(ESTE)
Jk \i l• Y
G
O'l
w
ECUADOR ELIPSOIDAL.
XG
a). - Sistema Geodésico y Sistema Geodésico Local. b). - Distancia Espacial ( dij )
Azimut ( a;1 ) y Angulo Vertical
(aij) en el Sistema Geodésico Local.
FIGURA 26
PROBLEMA DIRECTO ( GEODESICO LOCAL) .

Ahora el vector de posición de un punto (J) en el sistema geodésico
local en (i)(Fig. 26) está dada por:
ｾﾡ＠
(¡:ij IGL" djj
cos ajj
cos ajj
",a]
sen a jj - - - - - - - - 189
djj sen ajj
y (rij)G puede escribirse
(f¡j IG =R 3 (180-ATI R2(90-<pjl ｾＨｦﾡｪｊ＠ -------190
La matrÍz de reflexión P2 y dos matrices rotacionales R2 y R3

transforman el vector topocéntrico del sistema geodésico local al siste-
ma geodésico. El vector de posición del segundo punto (J) es obtenido -
por suma vectorial como:
(Fj IG = ni IG+ (rij)G - - - - - - - - - 191
Donde (rij) está dado por (190) y (ri)G es el vector de posición del
punto dado (O. Como se mencionó previamente, las coordenadas geodé-
sicas ( <Pi Aj h¡ ) pueden obtenerse via, una simple transformación de -
coordenadas.
El procedimiento para el cálculo del problema directo cuando el -

azimut y el ángulo vertical son dados en el sistema astronómico local -
(Fig. 27) es completamente análogo al descrito con respecto al sistema
geodésico local anterior. La única diferencia es en la expresión usada
para calcular el vector de posición topocéntrico Fij. En éste caso, está
dada por:
1
,
(F.. )
1) G
= Rs(180 -A-l
1
R (90 - Á.·)P (,f.. )
2 't'1 2 IJ A.L,
1-------- - 192

Donde <Pi y A i son la latitud y longitud astronómicas del punto
dado, y.
(r:. )
ｾｩﾡ＠
::: dij
COSV
ij
casv ij senA i) ］ａﾡｪｾ＠ --- - - - - 193
1) AL ld"
IJ SenVij
64
ZAL
VAL
ZAL
en
01
VAL
XII
a). - Sistemas Geodésico y Astronómico b). - Distancia Espacial (dij) Azimut (Aij)
Local. y ángulo vertical (Vij) en el Sistema
Astronómico Local.
FIGURA 27
PROBLEMA DIRECTO (ASTRONOMICO LOCAL)

Note que en éste caso (192) el vector de posición está girado directa-
mente del sistema astronómico local al sistema geodésico. Una trans-
formación alternativa es posible vía el sistema geodésico local usando -
la expresión:
(r.. ) =Rf180- A.)
1) G ; j 1
R (90-,.1.) P R (A .. - a. ) Rf-f ) R.("",) (r:')A L - - - - - - - 194
2 'tí 2 5 1) 1j ｾ＠ 1 l ' fI 1) .
En la expresión anterior (194), Aij Y aij son los azimutes astronó

mico y geodésico respectivamente y las cantidades ｾｩ＠ y 'T) i son las dos
componentes de la desviación de la vertical en el punto i.
5 • 2 • - Problema Inverso. -
En éste caso, las triadas de coordenadas ( <P, A, h ) Ó ( x, v, Z

están dadas por dos puntos en el terreno. Son requeridas la distancia -
espacial (dij) los azimutes directos (dij) e inverso (ají ) Y los ángillos -
verticales (aij) y (aji). .
Los vectores de posición de los dos puntos (i) y (j) en el sistema -

geodésico están dados por:
X-JI [Ni + hi) cos <Pi COS Ajl

('r¡)G = V¡ = (Ni+h¡) cos<p¡ SenA¡ --------195
[
Z¡ N¡(I-e2 )+h¡)sen <p¡
G
y

Primero, el vector diferencia Fij en el sistema geodésico está deter-
minado por:
(roo)
1) G
=('f.) -(r·):
) G 1GIl
X. ｸｾ＠
/-y. = ｛ａｘｾｪ＠ Ayl.
1)
-------- -197
G G U Zj Z¡ AZij
Enseguida, el vector diferencia anterior es girado al sistema de coor-

denadas geodésico local vía una expresión que es inversa de (190) y está
dada por:
66
A.hora, para determinar la distancia espacial, el azimut y el ángulo
vertical en (i) usamos las componentes del vector (rij) en las expre
siones: G. L -
[2
di) = 6X ij + 2 + 6Z2JI'z
6Yij - - - - - - - 199
ij
a ij = ton-1 [ 6Yij
6Xij
] -------200
y ajj = sen-1 [
6Zij
dij
] --------201
Las expresiones correspondientes para determinar el azimut a·· yel

ángulo vertical aji en el sistema geodésico local en ( j ) son: JI
(r.. )
)1
= p.2 R2 (A..-90)
GL '"t')
R(A.-180)
a )
(r.. ) - - - - - - - 202
)1 G
a =ton- 1r 6Yji J --------203

ji t 6Xji
.y
Oji = sen- I ｛ｾｊ＠
dij
-- ------ 204
6. - Problemas de Intersecci6n en Tres Dimensiones. -
El problema de determinar las coordenadas de un punto sobre un plano

usando una intersecci.6n de dos azimutes o distancias de dos puntos (coor
denados) conocidos. es un proceso directo (Faig. 1972) • Este tipo de
problemas no es tratado genralmente para cllculos sobre un elipsoide de
referencia. El problema de intersecci6n para la determinaci6n de las - -
coordenadas geodésicas ( ep, A ) puede ser tratado con harta simpleza
usando el álgebra vectorial. Son presentados aquí dos casos. cada una -
requiere en formaci6n similar a la que se requeriría para dos cálculos di
mensionales rigurosos. -
6.1. - Intersecci6n A.zimutal.-
El problema es definido como: Dadas las triadas de coordenadas ( - -

( ep.,A., h. ) Y ( epi ' Aj , h J ) para dos puntos en el terreno
(O y (h Y los azimutes de la ｳｾ｣ｩＶｮ＠ normal del terreno aik y a¡k de
los puntos conocidos al punto desconocido (k) calcular las coordenadas - -
geodésicas epk y Ak del punto desconocido k. Note que la altura -
elipsoidal aproximada h k es requerida para los cálculos.
Para comenzar la solucion es necesario defi-w-r un vector unitario en -

un azimut cualquiera. Este ｶ･｣ｴｾｲ＠ es denotado fa y se expresa en tér -
minos de los vectores unitarios fLx y ;y que son respectivamente las
direcciones Norte y Este del sistema geoa.ésico local. (Fig. 28)
67
A
TQ
A
/'-v
ｾＭ＠
ECUADOR ELIPSOIDAL
F I GU R A 2 8

VECTORES UNITARIOS EN EL SISTEMA GEODESICO
LOCAL.
68
Este está dado por la ecuación.
ｾ＠ｾ＠ ,
to = fLx cos o + fLy sen o - - - - - - - 205
donde:
ｾ＠ ,
Usando las expresiones para fLx y fLy' (205) puede ser reescri-
ta como:
[!; = ＭＺ［ｾ＠｛
tz
- G-L
Ｚｾ＠Ｚｾ＠
cosq, cos o
ＺｾＭ
ｾｉ＠
-G-L
208
Ahora un vector unitario perpendicular al azimut o es de-

finido por:
A A A
t =u cos(ot90l+1L sen(0+90l--------209
0+90 'x 'Y
Para resolver para q,k Y Ak , deben ser formuladas dos
ecuaciones donde aparezcan explícitamente estas dos cantidades.
Primero, dos productosescalaresson formados, cada uno involu-

cra un vector en un plano definido por un par de puntos en el terreno y el o-
rigen del sistema coordenado, y un segundo vector que está en un azimut a
90° de éste plano, (Fig. 29). Los dos productos básicos son:
(rk -ri l· ｾ｡ｩｫ＠ +90°=0 - - - --- - 210

_ _ fI
y (rk-rjl- t Ojk +90°=0----------211

donde:
ｕ
ｉＧｬｘｩｾ＠
I'lYik - - - - - 214
lIZik
69
ｾＭＫｔｴｾ＠
FIGURA 29 INEGI. Cálculos de Posicionamiento Geodésico

INTERSECCION AZIMUTAL EN TRES DIMENSIONES.
70
y
XiG
Ü 8Yjk
8Zjk
------- 215
En las ecuaciones (214) .y (215) las coordenadas para (i) y (j ) son -

consideradas constantes, mientras que aquellas para (k) son consideradas
dadas por tres funciones desconocidas (Krakiwsky and Wells, 1971).
_ ｾｘ＠ [ces {3k cos Ak + hk cos </>k cos Aj

rk = Yk = a cosf:3k senAk +hk cos</>k senAk - - - - - - 216
Zk b sen f:3k + hk sen </>k
Primeros términos de (216) dan las coordenadas de (k) sobre la su-

perficie del elipsoide (def"wido por los semiejes mayor y menor (a) y (b)-
respectivamente) en términos de la latitud reducida ( f:3k ) Y la longitud geo
désica (Ak ). Los segundos términos se toma en cuenta el hecho de -::
que el punto (k) sobre el terreno está localizado en una altura elipsoidal -
(hk ) sobre el elipsoide de referencia, y están expresados en términos de
la latitud geodésica </>k y longitud Ak
Ahora las ecuaciones (210) y (211) pueden reescribirse como:
t I = AX'k
I
t,
XI
+ tJ.Y.Ik t,
yl
+ b.Z'k
I
t.
ZI
= O - - - - - - 217
fe = b.Xj!I\J +8V,k I ' +b.Z jk tz' = O - - - - - 218
J
Las cantidades desconocidas en las ecuaciones ari1erioreJ son las -
coordenadas de k·y en términos de estas (217) y (218) son no lineales.-
La siguiente etapa en la soluci6n es aproximar las ecuaciones (217) y ｾ＠
(218) por una serie lineal de Taylor usando valores aproximados para la
latitud reducida y longitud denotadas por f:3°k y A"k respectivamente,
entonces:
° 0'1 t:J
t =f + - - d""k+ - - dAk +-----=0 ---------219
Ot, .
I I af:3k OAk
f =fO+ afz ､ｦＺＳｫＫｾａＭ］ｏ＠ ------ 220

e 2 af:3k oAk
Donde
fO=b.X· t + 8VO t + b.ZO t --------- 221
I ik xi ik yi ik zi
71
AXo = o cos f3 0
cos >..0 + hO cos,4..o cos >..0 - X - - - - - - 223
Ik k k k 't'k k i
AY¡'k = o cos f3 k sen >"°k + h'k cos epk sen >"°k - Y¡ - - - - - - 224
AZfk = b sen f3'k + h'k sen ep'k - Zi - - - - - - - - - - - - - - - 225
ｦＺＮｘｾ＠ = o cos f3 0 cos >..0 + h'" cos epo cos >..0 - X --- - - --- 226
¡k k k k k k j
f:.Yk = o cos ¡9;
sen Á-k + h; cos 4>k sen >"'k. - 227 '1 ------ -
r.zojk =bsen F3°+hosenepo
k k k
-Z --------------228
J
O '
ｏｾｫ＠ = txi 1- a sen f3k cos >"°k - hksen <Pk cas >"°k 1+,
+ tyi (- a sen f3k sen ';. .<>k - h'k sen ･ｰｾ＠ sen >"°k 1 +
+tzilbcos ¡9k
+h'k cOs<P'kI------------- 229
afl 1_ Qo \0 ° ,4..0 \0
O>"k - Ixi -o cos I-'k sen "k - hk ces 't'k sen "'k 1+
+ t yi 1acos f3k cos >"'k +h¡; cos <P¡; cos >"°k 1- - - - - - - 2 30
Ahora reescribiendo (229) y (230) como:
afl
af3 =txi Xf3+tyi Yf3+ tzl Zf3-------231
k
Ofl _
ｾ＠ - txi X>.. + tyi y>.. -------- - - -- 232
y
--º.h....:
O{3k t· Xf3+
X¡
t· Yt.:l+l· Z-------o--- 233
y¡ 1-' Z¡ f3
Ofz

ｾ＠ "Ixj X>..+ty¡Y>.. ----------234
Se notará que tomando derivadas parciales, la latitud geodésica <Pk
fue considerada sin6nimaa la latitud reducida f3 k • No hay pérdida en -
exactitud en cálculos subsecuentes debido a ésta consideración. Adicion-
almedl:e, un valor aproximado de hk que es dentro de 100m del valor - -
verdadero es Slficiente.
Las ecuaciones (235) Y (236) son resueltas en un proceso iterativo

hasta que las correcciones a f3 k y >"k sean despreciables ( < O, 0001 ").
El valor de la latitud geodésica <Pk es entonces resuelta por (Krakiwsky
and Wells, 1971): o
72
1-1 ep_k_=_to_ñ'_[_ ｾ＠ ton ,8J------ --- - 237
6.2. - Intersección en la Distancia Espacial.-
La deternúnación de la latitud geodésica (epk ) Y la longitud (Ak )

de un punto en el terreno usando dos distancias espaciales terrestres, -
se resuelve de una manera sinúlar a la usada para una intersección azi
mutal ( 6.1). Son dadas las dos triadas de coordenadas ( epi I A¡ I h¡ )"
Y ( epo I Ao I h o ) y dos distancias espaciales terrestres r¡k y r"k
de lo;} punths donocidos al punto desconocido (k). Además, una Jtura
elipsoidal aproximada h"k se requiere (Dentro de los 100 m del valor de
h
k es suficiente).
La llave para la solución es la formación de dos ecuaciones linea -
les que están expresadas en términos de los parámetros conocidos y des
conocidos (Fig. 30). Comenzamos con las relaciones. -
fl " [(Xk- X¡)2 +(Yk-Y¡)2HZk-Z¡)2J ｾＲ＠ -r¡k =0----------238

ｾＲ］＠ [(Xk-X j )2 ＫＨｙｫＭｾＩｚｪ＠
2 Ｒ｝ｾ＠
-rjk =0--------239
Donde ( Xk I Yk I Zk ) están dadas por (216). Las ecuaciones an-
teriores son no lineales en térllÚnos de,8k y Ak por lo que son apro:a;i
madas por una serie lineal °de Taylor usando valores aproximados pa
ra la latitud reducida,,8; y la lonltitud geodésica ａｾ＠ La forma
lineal de las ecuaciónes (238) y (239) están dadas por:
f = fO + ｾ＠ d,8 + ｾ＠ dA +---- =0 - ------ 240

I I a,8k k aAk k
y
Donde:
73
F 1 G U R A 30
INTERSECCION DE LA DISTANCIA ESPACIAL EN

TRES DLMENSIONES.
tlO = rik - rik - ------ 242
t; = rjk - rjk -- - -- - - - 243
ｾ＠
Of3k
= _ 1 - [lX O-X.l OXk +lY o -Y.l OYk +(ZO -y )
r¡k k I Of3k k I af3k k i
ｾｊ＠
Of3k
-------- 244

-ª..L = _1_ [lXO -X.) ｾ＠ +lYo-y.) OYk + lZo -y.) OZk ] --- - - --- 245
OAk rik k I OAk k I OAk k I OAk
-ª.L. = _ 1 - [lX O -X.l OXk +lYo-y.) OVk + lZo -Y.l OZk ] ---------246
Of3k rjk k 1 Of3k k 1 af3k k 1 Of3k
ｾ＠
OA2
= _1- [(X
rjk k- Xjl OXk
OAk
+ (yo -y, ) OYk + (ZO - Y.) OZk ] -- - __ - ___ 247
k 1 OAk k 1 OAk
Ahora los términos en las ecuaciones (244) - (247) son derivadas de

(216) y están dadas por:
74
OXk
- - :-0 sen 13° COSAo -he sencj>° cos AO = X - -- --- - 248
iJf3k k k k k k f3
-ª..xL =-o sen f3ksen A'k - h'k sen cj>'k sen A'k = Y{3 - - - - - - 249
iJ{3k
ｾ＠ O{3k
= b cos{3° +ho coscj>° =Zf3 --- -------- --- 250
k k k .
iJXk =_ o cos {30 sen AO - hO cos cj>0 senA" =X - - -- -- - - 251

iJAk k k k k k A
iJYk
- - = a cos {30 cos>e + hOcoscj>° cos AO =Y, - ------- - 252
iJAk k k k k k 1\
OZk
- - = 0-- - - - - - - - - - - - - - - - - - - - - - 253
iJAk
Como en el caso de la intersecci6n azimutal. ,la latitud geodésica,

4>k ' fué considerada sinonima con 1 a latitud reducida {3k •
Ahora (240) y (241) son reescritas para la solución como:
f. =fr++ [IlXik ｘｾａｙｩｫ＠ Y{3 +AZ ik l¡3]d{3k++[AXikXA+AYik YA] dAk --- - 254
ｲｾ＠ *
------255
Las correcciónes dB. y dA. se resuelven usando un procedimiento

iterativo. Cuando las correcciónes se convierten en despreciables -
( 0.0001") se obtienen los valores finales de B. y A. ,cj>. se obtiene
usando (237).
7. - A Manera de Conclusión.
A primera vista, parece que el método clásico de cálculos de posi-

ción geodésica sobre la superficie de un elipsoide de revolución debería
ser abandonado en favor del método tridirrensional.
Las fórmulas para este último son mas simples de derivar e imple
mentar y en el caso de los problemas directo e inverso, están dadas en -
forma definida. Además, si las coordenadas curvilíneas (problema direc
to) o la distancia elipsoidal y azimutes de sección normal (problema inve!.
so) son requeridos, fórmulas de transformación rigurosa están disponi -
bIes para obtenerlas (Krakiwsky and Wells 1971, Sección II).
La dificultad mayor para usar el método tridimensional yace en las

observables geodésicas o en su carencia. Esto es particularmente cierto
en el caso del problema directo o en cualquier problema donde es requeri-
do el ángulo vertical (90° distancia zenital). Debido a los problemas de -
refracción, la distancia zenital no puede obtenerse mejor que 2: 1" 10-
75
que sobre una línea de 10 kms. da una desviación estandar en altura de
10 cms. (Heiskanen and Moritz, 1967)· Este error obviamente -
podría afectar los cálculos de las coordenadas tridimensionales ( X, Y,
Z) ó ( cf>, x. , h ) de un punto requerido. El problema puede ser su -
perado con nivelación geométrica pero es improbable que éstas observa
ciones se hicieran siempre disponibles. -
Los dos problemas de intersección que han sido presentados mues-

tran como el método tridimensional puede usarse para resolver directa
mente las coordenadas curvilineas. Sería obvio que si fuera disponible
suficiente información observada ( P. ej, distancias triespaciales) los -
problemas podrían ser formulados y resueltos diarectamente en térmi"\
nos de las coordenadas cartesianas tridimensionalés.
Finalmente se notará que una cantidad equivalente de información -

observada se requiere para los métodos clásicos y tridimensional. La
diferencia principal es que para los cálculos elipsoidales (P. ej, proble-
ma directo) la altura elipsoidal no necesita, ser conocida en forma tan -
precisa como en los cálculos tridimensionales. Sin embargo, no impor
ta que método se use, transfornaciones rigurosas mostrarán que los re
sultados son equivalentes. Esto es, las coordenadas cartesianas (X, -
Y, Z ) darán un conjunto ( cf>, x., h ) en que la latitud geodésica (cf> )
y la longitud (X. ) son iguales a las obtenidas por cálculos clásicos.
Aún más las distancias espaciales y los azimutes de 'las secciones nor
mal del terreno obtenidos de cálculos tridimensionales (Proble :-
ma inverso) y reducidos rigurosamente al elipsoide de referencia... son i -
guales a las distancias elipsoidales y azimutes geodésicos obtenidos del
problema inverso resuelto sobre el elipsoide.
76
REFEREN CIAS
Abroncowitz, M. and l. Se:;un (5th ed., 1968). Randbook of Mathematical

Functions. Dover Publication Inc., ]';ew York.
Allan, A. L., J. R. Hollwey, and J. H. B. Maynes (1958). Practical -

Field Surveying and Com putations. W illiam Heinemann Ltd.,
Toronto.
Bomford, G. (3rd ed., 1971>. Geodesy. Oxford University Press,

London.
Faig, w. (1972) Advanced Surveying 1. Lecture Notes No. 26, Depart--

ment of Surveying Engineering. University of New Brunswick,
Fredericton.
Heiskanen , W. A. and H. Moritz (1967) Physical Geodesy. W. H;

Freeman and Com pan y , London.
Jordon, W. and O. Eggert (1962). Handbuch der Vernessengskunde, Ed.

IIl. English Translation, Army Map Service, Washington.
Krakiwsky, E. J. and D. E. Wells (1971). Coordinate Systems in Geode-

sy. ｌｾ｟｣ｴｵｲ･＠ Notes No. 16, Department of Surveying Enginee-
ringo University of New Brunswick, Fredericton.
Lipschutz, M. M. (1969). Differential Geometry. Shaum 's Outline

Series, McGraw-Hill Book Company, Toronto.
Phillips, H. B. (2nd. ed., 1957). Analytical ｇ･ｾｴｲｹ＠ and Calculus.

John Wiley & Sons Ine., l\'ew York.
Rainsford, H. F. (1955). Long Geodesics on the Ellipsoid. Bulletin

Geodesique, No. 37.
Robbins, A. R. (1962). Long Lines on the Spheroid. Empire Survey

Review, Vol. XVI, No. 125.
Sodano, E. M. (1963). Generall\'on-Iterative Solution oí the Inverse and

Direct Geodetic Problems. ｒ･ｳｅｉＮｩＺｬｾ＠ I',otes I':o. 11, l:, S. -
ArE:!Y Development Agency. Fort Belvoir, Virginia.
Thomson, D. B. and P. Vanieek (197,';). A I\ote on the Reduction of

Spatial Distances to the Ellipsoid. ｓｵｾｶ＠ Review,_ Vol._XXII,
r\o·_11.!h
Vanicek, P. (1971). Physical Geodesy 1. ｌ･｣ｴｵＡＺｾｎｯｳ＠ No. 21, Depart-
m ent oí &lrveying Engineering. Uni versity oí N ew Brunswick.
Fredericton.
Vanicek. P. (1972). Physical Geodesy JI. Lecture Notes No. 24, Depart-
ment oLSurveying Engineering. University oí New Brunswick.
Fredericton.
Vanicek. P. and C. Merry (1973). Determinatíon oí the Geoid from DefleE.

tions of the Vertical Using a Least-Squares Fitting Technique.
Bullet in Geodesi9ue, No. 109.
Zakatov. P. S. (1953). A Course in Higher Geodesy. Published for the-

National Science Foundation, Washington, D. c. by the Israel
Program of Scientific Translations, Jerusalem, 1962.
78
Segunda reimpresi6n: mayo de 1983.
500 Ejemplares.
Direcci6n General d e Geografía.

San Antonio Abad No. 124.
México, D. F.
ISBN 968-809-921-X
ﾡＺｲｾＴ＠
5PP
IN$lITUTO NACIONAL oe ESTADISTICA
GEOGAA'IA e INFORMATleA

Radios de Curvatura

Cargado por

Copyright:

Formatos disponibles

Radios de Curvatura

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Radios de Curvatura

Cargado por

Copyright:

Formatos disponibles

ECUADOR ELIPSOIDAL

INSTITUTO NACIONAL DE ESTADISTICA

INSTITUTO NACIONAL DE ESTADISTICA

Este trabajo es parte del esfuerzo que DETENAL está haciendo

con el prop6sito de elevar el nivel de los conocimientos geodésicos

dentro y fuera de la Instituci6n. Debe pues, agradecerse la ､ｩｳｰｾ＠

sici6n y el apoyo brindado por las autoridades de DETENAL, parti-

cularmente de aquellas responsables del área de Geodesia, que al

facilitarnos medios y personal, hicieron posible que éstas notas ｶｩｾ＠

ran la luz del día.

Debe agradecerse profundamente la gentileza de los autores,

Dr. E. J. Krakiwsky y Dr. D. B. Thomsón, de la Universidad de

New Brunswick, al dar su autorización para la traducci6n y divulg.!

ci6n de su obra en español.

El excelente trabajo de mecanografiado estuvo a cargo de la

Srita. Concepción Vega Ché de la Oficina de Apoyo Básico. El Sr.

Julio Bueyes Oliva, de la Oficina de Apoyo Vertical, tuvo la respon-

sabilidad de trazar los diagramas y escribir las f6rmulas que ｡ｰｲｾ＠

cen en éstas notas.

M. C. Rafael Sosa Torres.

El propósito de estas notas es dar la teoría y el uso de algu-

El enfoque que se da aquí es realizar derivaciones completas,

Las notas se desarrollaron de las conferencias preparadas -

Los autores desean reconocer la contribución hecha por los

INEGI. Cálculos de Posicionamiento Geodésico

Sección n. - Reducción de Observaciones Geodésicas Terrestres --- 30

Sección 1II. - Cálculos de Posiciones Geodésicas sobre el ElipSOide

INEGI. Cálculos de Posicionamiento Geodésico

4. - F6rmula de Bessel. - Líneas Largas ----------------.------ 47

7. - A manera de concluswn ---------------------------- 75

INEGI. Cálculos de Posicionamiento Geodésico

FIG. No. PAJINA

2 Sección Normal Meridiana Mostrando el -

3 Sección Normal al primer vertical mos -

4 Radio de curvatura del meridiano (M) 17

5 Radio de curvatura del primer vertical (N) 19

6 Sección normal en un Azimut Cualquiera C! 21

7 Indicatríz para la solución de Ra 21

8 Sección a través de PP' (o) para la solu -

9 Solución de Z para solucionar Ro 23

10 Secciones con normales recípr,o.cas 24

11 Sección Triangular con normal recíproca 25

12 Separación Angular entre secciónes con Nor

INEGI. Cálculos de Posicionamiento Geodésico

15 Separación entre Sección Normal y Geodésica 3D

16 Corrección por Normal Oblicua 32

17 Corrección por Desviación de la Vertical 34

18 Reducción de Distancia Espacial 36

19 Fórmula de Puissant para Problema Directo. 39

20 Fórmula de Puissant para Problema Inverso 44

21 Relaciones fundamentales para el desarrollo

22 Esfera Reducida y Elipsoide. 49

25 Soluci6n de Longitud de Arco (J'" 59

26 Problema Directo (Geodésico Local) 63

27 Problema Directo (Astron6mico Local) 65

28 Vectores Unitarios en el Sistema Geodési-

29 Intersecci6n Azimutal en Tres Dimensiones. 70

30 Intersecci6n de Distancia Espacial en Tres

Las primeras tres secciones de estas notas tratan del cálculo

Dos clásicos problemaa de cálculo de geodesia geométrica -