NOTAS DE CLASE GEODESIA

GEOMÉTRICA LATEX

ANDRÉS CARDENAS CONTRERAS

UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS

28 de noviembre de 2014

1
 Notas De Clase Geodesia Geométrica Prof:Andres Cardenas

Índice
INTRODUCCIÓN 4

1. CONCEPTOS 5
1.1. Geodesia Geométrica . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.1. Geoide . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.2. Elipsoide . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.1.3. Elipse Meridiana . . . . . . . . . . . . . . . . . . . . . . 7

2. SISTEMAS DE COORDENADAS 9
2.1. Geodésicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2. Geocéntricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.3. Paramétricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.4. Relaciones entre los sistemas de coordenadas . . . . . . . . . . . 14
2.5. Problema inverso de coordenadas rectangulares . . . . . . . . . 15

3. NIVELACIÓN 19
3.1. Trigonométrica . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.1.1. Reciproca . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.1.2. No Reciproca . . . . . . . . . . . . . . . . . . . . . . . . 22
3.2. Trisección . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.3. Bisección . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.4. Diferencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28

4. LONGITUD DE ARCO 31
4.1. Longitud de arco meridiano . . . . . . . . . . . . . . . . . . . . 31
4.2. Longitud de arco paralelo . . . . . . . . . . . . . . . . . . . . . 35
4.3. Área de un cuadrilátero sobre el elipsoide . . . . . . . . . . . . 36

5. RADIOS DE CURVATURA 36
5.1. Radio de curvatura de la sección meridiana (ρm ) . . . . . . . . . 37
5.2. Radio de curvatura de un paralelo . . . . . . . . . . . . . . . . . 39
5.3. Radio de curvatura en el primer vertical . . . . . . . . . . . . . 39
5.4. Radio de curvatura de una sección normal en un azimut cualquiera 40
5.5. Valor medio curvatura . . . . . . . . . . . . . . . . . . . . . . . 40

6. PROBLEMA GEODÉSICO 43
6.1. Secciones Normales . . . . . . . . . . . . . . . . . . . . . . . . . 43
6.2. Lı́nea Geodésica . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
6.3. Método Puissant . . . . . . . . . . . . . . . . . . . . . . . . . . 49
6.3.1. Método directo . . . . . . . . . . . . . . . . . . . . . . . 49
6.3.2. Método inverso . . . . . . . . . . . . . . . . . . . . . . . 56
 Notas De Clase Geodesia Geométrica Prof:Andres Cardenas

Índice de figuras
1. Figura 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2. Figura 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3. Figura 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
4. Figura 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
5. Figura 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
6. Figura 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
7. Figura 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
8. Figura 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
9. Figura 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
10. Figura 10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
11. Figura 11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
12. Figura 12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
13. Figura 13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
14. Figura 14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
15. Figura 15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
16. Figura 16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
17. Figura 17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
18. Figura 18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
19. Figura 19 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
20. Figura 20 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
21. Figura 21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
22. Figura 22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
23. Figura 23 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
24. Figura 24 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
25. Figura 25 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
 INTRODUCCIÓN

Las siguientes notas de clase, son producto del curso de Geodesia Geométrica
que imparte el docente Andrés Cárdenas Contreras en la Universidad Distri-
tal Francisco José de Caldas en el proyecto curricular de Ingenierı́a Catastral
y Geodesia; para dicho espacio académico es necesario manejar algunos con-
ceptos y metodologı́as básicas de Topografı́a, Astronomı́a Esférica, además de
la implementación de otras áreas como la matemática, la fı́sica y las ciencias
computacionales.

Esta recopilación de las notas de clase tienen como propósito brindar a los
estudiantes de Geodesia Geométrica una herramienta que proporcione los con-
ceptos básicos y demás elementos que sean necesarios para el desarrollo e
interpretación de procedimientos que aquı́ se expongan; por tanto estas notas
tienen un fin completamente pedagógico, ya que serán un gran apoyo a lo largo
de tal curso.

Para efectos de digitalización se usó el procesador de texto cientı́fico LATEXque

se compila sobre TEX, facilitando ası́ la escritura de posteriores expresiones
matemáticas, y la inserción de gráficas y figuras, los colaboradores en el pro-
ceso de edición fueron Jhon Alexander Galindo Ambuila y Danna Juliette
Castañeda Franco, estudiantes de geodesia geométrica en el segundo semestre
del año 2014

4
1. CONCEPTOS
1.1. Geodesia Geométrica
El objeto de estudio de la geodesia es determinar las dimensiones y forma de
la tierra, que corresponden a dimensiones en el espacio no sólo del sistema tie-
rra. Para lograr este propósito se requiere definir las ecuaciones matemáticas
de la forma geométrica (geodesia geométrica), las relaciones fı́sico - matemáti-
cas (geodesia fı́sica) correspondientes a un modelo geopotencial, las actuales
técnicas de posicionamiento GNSS (Global Navegation Satellite System) y el
estudio de la precisión de los datos obtenidos (ajustes geodésicos).
Esta ciencia se encarga de determinar las posiciones de los puntos de la su-
perficie de la tierra, esto significa los valores numéricos y las coordenadas que
tienen asignados estos puntos, pero estas coordenadas deben estar asociadas
a un sistema de referencia. Lo que implica involucrar un problema que busca
encontrar el modelo de la forma de la tierra, según sea su interés de represen-
tación: Plana, Esférica, Elipsoidal o Geoidal, La geodesia geométrica involucra
las relaciones geométricas del elipsoide de rotación. Su objeto de estudio tiene
dos aspectos principales:
La consideración que la forma de la tierra es un elipsoide de revolución
biaxial.
Es posible representar la superficie de la tierra, en mapas cartográficos,
a través de cartografı́a matemática fundamentada en el elipsoide.

1.1.1. Geoide
Corresponde a una superficie equipotencial del campo gravitatorio de la tierra.
Dicha superficie fı́sica real, puede considerarse coincidente con la superficie del
agua en reposo de los océanos idealmente extendida bajo los continentes y
sobre la cual la gravedad en todos sus puntos es normal a ella.

Figura 1: Geoide. Tomado de Wikipedia

5
1.1.2. Elipsoide
EL elipsoide es un modelo matemático de la forma de la tierra, se asumió ası́,
considerando el achatamiento de los polos y corresponde a la rotación de una
elipse alrededor del eje z. Si el modelo se observa desde un punto de fuga
normal al plano xy, es decir en la dirección de z se tiene una circunferencia de
radio a.

Figura 2: Elipsoide

El Elipsoide es una superficie Gaussiana cerrada, No Euclidiana, que posee dos

parámetros: semiejes y Achatamiento.

siendo:

a: semieje mayor
b: semieje menor
c: distancia del origen al foco

Elipsoide Biaxial:
Se define el elipsoide biaxial de revolución como la región del espacio en
la cual se cumple que:
a=c y b<a
Por tanto,
x2 y2 z2
a2 + a2 + b2 =1
Elipsoide triaxial:
a>c b<a y c 6= a

6
1.1.3. Elipse Meridiana
Considérese que una elipse es una figura cerrada, continua y derivable, en la
cual se cumple que todos los puntos, están distanciados de sus dos focos un
valor constante k = 2a.

Figura 3: Elipse meridiana

Siendo:

a: semieje mayor
b: semieje menor
c: distancia del origen al foco

Se define en R2 una Elipse centrada en el origen del sistema, con semiejes

paralelos a los ejes coordenados, que tendrá ecuación:

x2 z 2
+ 2 =1 (1.1)
a2 b

Sea c la distancia del centro a uno de los focos. Tenemos por definición, las
excentricidades de la elipse:

Primera,
c
e= a

Segunda,
c
e0 = b

7
De la relación que vemos en la Elipse tenemos que a2 = c2 + b2 , por tanto
definimos:

a2 − b 2
e2 = (1.2)
a2

a2 − b 2
e02 = (1.3)
b2

Y la relación con el aplanamiento de la elipse:

a−b
f= (1.4)
a

Y su equivalencia con:

e2 = 2f − f 2 (1.5)

Con dos de los parámetros mencionados, generalmente a y f , se define un elip-

soide; antes del desarrollo de los actuales sistemas de posicionamiento satelital
GNSS, se estableció en cada región o paı́s, según las caracterı́sticas propias,
diferentes parámetros para sus elipsoides, conocidos como dátum locales. En
Colombia, se utilizan los parámetros del Datum Bogotá, que corresponden
al elipsoide internacional y de los dos dátum geocéntricos más importantes
GRS80 y WGS84. En la tabla que se muestra a continuación se muestran los
valores.

Datum a = Semieje mayor b = Semieje menor Aplanamientof − 1

Internacional 6378388 6356911.9462 297
GRS 80 6378137 6356752.3141 298.257222101
WGS84 6378137 6356752.3141 298.2572223563

8
2. SISTEMAS DE COORDENADAS
Un sistema de coordenadas se define a partir de un meridiano y un paralelo
establecidos por convención internacional, estos son el meridiano de Greenwich
y el plano ecuatorial respectivamente; a partir de ellos se puede determinar la
posición de un punto en especı́fico. Cabe resaltar que para todos los sistemas
de coordenadas es posible definir su altura elipsoidal (h), esta corresponde a la
distancia geométrica entre la superficie del elipsoide y la posición de un punto
T en la superficie topográfica sobre el vector normal.
Para cualquier sistema de coordenadas, se define a la longitud (λ) como el
ángulo diedro entre los planos del primer meridiano (Greenwich) y el meridiano
que pasa por el punto de interés.

2.1. Geodésicas
Latitud Geodésica-Geográfica (ϕ): Se define como el ángulo diedro entre
el Ecuador y el vector normal al plano tangente.

Figura 4: Coordenadas geodésicas

Usando la Ecuación (1.1) podemos multiplicar por a2 b2 y luego obtener,

b2 x2 + a2 z 2 = a2 b2
−dx a2 z
2b2 xdx + 2a2 zdz = 0 ⇒ dz = b2 x
De la figura anterior podemos obtener la siguiente información:
−dx a2 z
tan ϕ = dz = b2 x
9
con lo cual,

b2
z = 2 x tan ϕ ⇒ z 2 = (1 − e2 )x2 tan2 ϕ
a
llevando este resultado a (1.1):
2 )x2 tan2 ϕ
x2
a2
+ (1−e b2

a2 cos2 ϕ
x2 =
(1−e2 )sin2 ϕ

a2 cos ϕ
x= 2
p (2.1)
(1 − e2 )sin2 ϕ

despejando x en z, obtenemos:

a(1 − e2 ) sin ϕ
z=p (2.2)
(1 − e2 )sin2 ϕ
sı́
p
ω= 1 − e2 sin2 ϕ

se define la gran normal (N),

a
N= ω

x = N cos ϕ (2.3)

z = N (1 − e2 ) sin ϕ (2.4)

10
2.2. Geocéntricas

Figura 5: Coordenadas geocéntricas

Latitud Geocéntrica (ω): Es el ángulo entre el vector que une a un punto

T y el origen con el ecuador.

ω 6= ϕM

Z
tan ω = Xm
m

Xm = Rg cos ω

ZM = Rg sin ω

llevando las cantidades mencionadas antes a (1.1):

(Rg cos ω)2 (Rg sin ω)2 1−sin ω

a2
+ b2
= Rg2 [
a2
+ sinb2ω ]
1−sin ω
= Rg2 [
a2
+ sinb2ω ]
luego
√
a 1 − e2
Rg = p (2.5)
1 − e2 sin2 ω

11
 Figura 6: Coordenadas geocéntricas

Zt
sin ω = R xz
g

Xt
cos ω = R xz
g

Xt = Rg xz cos ω

Zt = Rg xz sin ω

2.3. Paramétricas
Latitud Paramétrica (θ): Es el ángulo cuyo vértice se ubica en el centro de
una esfera que es tangente al elipsoide a lo largo del ecuador, entre el plano
ecuatorial y el radio del punto de interés, originado en la esfera por una linea
recta perpendicular al plano del ecuador, que es ademas un vector normal al
ecuador y que pasa por el punto dado y corta la circunferencia de radio a.

En la tabla que se encuentra a continuación, es posible resumir las relaciones

matemáticas de las coordenadas sobre la elipse meridiana y el elipsoide en
función de las latitudes.

12
 Figura 7: Coordenadas paramétricas

X
cos θ = at

Xt = a cos θ

Zt = b sin θ

TABLA LATITUDES DE LA ELIPSE MERIDIANA

Geodésicas Geocéntricas Paramétricas

Z = Rgxz sin ω
Z = N (1 − e2 ) sin φ √ Z = B sin θ
Rgxz = x2 + z 2

X = N cos φ X = Rgxz cos ω X = a cos θ

Z Z Z b
X
= tan φ(1 − e2 ) X
= tan ω X
= a
tan θ

13
 TABLA COORDENADAS SOBRE EL ELIPSOIDE LATITUDES

Geodésicas Geocéntricas Paramétricas

X = Rgxyz cos ω cos λ
X = N cos λ sin φ p X = a cos θ cos λ
Rgxyz = x2 + y 2 + z 2

Y = N sin λ sin φ Y = Rgxyz cos ω sin λ Y = a cos θ sin λ

Z = N (1 − e2 ) sin φ Z = Rgxyz sin ω Z = b sin θ

2.4. Relaciones entre los sistemas de coordenadas

Figura 8: Relación de los sistemas de coordenadas

De la Elipse Meridiana (fig. 8), tenemos υR = N cos ϕ y OR = Rg

Xm = Rg cos ω

Zm = Rg sin ω

14
Luego, por parte usando (2.3)(2.4) y la relación b2 = a2 (1 − e2 ). Y por otra
parte observando el Elipsoide por su parte superior, se puede establecer:
b
tan θ = (1 − e2 ) tan ϕ (2.6)
a

tan θ = (1 − f ) tan ϕ (2.7)

tan ϕ(1 − e2 ) = tan ω (2.8)

2.5. Problema inverso de coordenadas rectangulares

Conocidos ϕ, λ, h → x,y,z

Método 1
Dadas las coordenadas de un punto (p):(x,y,z) y los elementos del elipsoide
en mención: semieje mayor (a) y achatamiento (f). Encontrar las coordenadas
geodésicas: latitud(),longitud() y la altura del punto () sobre el elipsoide. A
partir del gráfico del elipsoide y de la elipse meridiana, se pueden encontrar
las siguientes relaciones:
Yp
tan λ = X
p

 
−1 Yp
λ = tan
Xp
p 2
υp = xp + yp2

pp0
sin ϕ = N +h

pp0 = (N + h) sin ϕ (2.9)

√ √
υp x2p +yp2 x2p +yp2
cos ϕ = N +h = N +h =⇒ N + h = cos ϕ
p
x2p + yp2
h= −N (2.10)
cos ϕ
pp0 q
tan ϕ = =⇒ pp0 = υp tan ϕ = tan ϕ x2p + yp2 (2.11)
υp
igualando (2.10) y (2.12):
p
tan ϕ x2p + yp2 = (N + h) sin ϕ

15
 zp +e2 N sin ϕ zp + e2 N sin ϕ
 
tan ϕ = √ =⇒ ϕ = tan −1 p 2
x2p +yp2 xp + yp2

e2 N sin ϕ
 
z

ϕ = tan −1 p 2p 2 1+
xp + y p zp
Si h = 0 tenemos en (2.11) que:

√
x2p +yp2
N= cos ϕ
  
2
ϕ0 = tan−1 √ z2p 2
1 + e N sin ϕ
zp
xp +yp

√ a
N1 =
1−e2 sin2 ϕ0

Ahora de manera general, iterar el procedimiento:

√ a
Nn =
1−e2 sin2 ϕn−1
 
zp +e2 Nn sin ϕn−1
ϕn = tan−1 √
x2p +yp2

Para n = 1,2,3,...
Hasta que |ϕn − ϕn−1 | ≤ 10−9 Rad
p 2
xp + yp2
h= − Nn (2.12)
cos ϕ

ϕ = ϕn (2.13)
Método 2
Usando la misma elipse meridiana usada para el método anterior, tenemos:

υp
cos ϕ = N +h

υp = (N + h) cos ϕ (2.14)

zp = [N (1 − e2 ) + h] sin ϕ (2.15)

16
Dividiendo (2.15) entre (2.14):

zp N (1−e2 )+h
υp = tan ϕ N +h

zp N +h
tan ϕ = υ
p N (1−e2 )+h
n o
zp N +h
ϕ = tan−1 υp N (1−e2 )+h

Luego hallamos una expresión para ϕ

  −1 
−1 zp N e2
ϕ = tan 1−
υp N +h
Ahora para hallar φ, Se usa la expresión
Yp
tan λ = X
p

tal cual se realizo en el método anterior para hallar ϕ hay que hacer un proceso
iterativo, como se indica a continuación:
p 2 √
h0 = xp + yp2 + zp2 − ab y N0 = a
 −1 
N0 e2

−1 zp
ϕ0 = tan 1−
υp N0 + h0
luego a continuación la iteración, de la siguiente forma:
a
Nn = p (2.16)
1 − e sin2 ϕn − 1
2

p 2
xp + yp2
hn = − Nn (2.17)
cos ϕn − 1

 −1 
Nn e2

−1 zp
ϕn = tan 1−
υp Nn + hn

Para n = 1,2,3,... Hasta que |ϕn − ϕn−1 | ≤ 10−9 Rad En tal caso,

ϕ = ϕn

h = hn

17
 Cuadrante Signo x Signo y Hemisferio Localización
1 + + E 0 ≥ λ ≤ 90
2 - + E 90 ≥ λ ≤ 180
3 - - W 90 ≥ λ ≤ 180
4 + - W 0 ≥ λ ≤ 90

Ejemplo

Dado un punto cercano a la Universidad Distrital Francisco José de Caldas

con coordenadas:

Latitud ϕ = 4◦ 370 41,01”N

Longitud λ = 74◦ 30 57,21”W

Altura h = 2591mts

Calcular sus coordenadas sobre el Elipsoide Internacional, para dicho Elipsoide

tenemos:
a = 6378388 b = 6356911,9462

f = 0,003367003356 e = 0,00672267

obtenemos

X = 1746107,48284

Y = −6115937,23274

Z = 511413,699017

18
3. NIVELACIÓN
3.1. Trigonométrica
A partir de observaciones de ángulos cenitales podemos determinar las alturas
ortométricas desconocidas del (los) punto(s) de interés.

Figura 9: Nivelación trigonométrica

De la figura se tiene que:

Z1 y Z2 Ángulos cenitales observados

P1 y P2 Estaciones geodésicas
H1 y H2 Alturas ortométricas
∆ξ Ángulo de refracción atmosférica
S Distancia geodésica entre las estaciones P1 y P2

Del triángulo P1 P2 L se tiene que:

◦ θ ◦
α + β + 90 + 2
= 180

19
De la figura también se tiene que:
◦
Z2 + ∆ξ + β = 180
◦ θ ◦
Z1 + ∆ξ + α + 90 − 2
= 180
◦ θ ◦
Z2 + ∆ξ + α + 90 − 2
= 180
◦ θ ◦
−Z1 − ∆ξ − α + 90 + 2
= −180

Z2 − Z1 − 2α = 0
Z2 −Z1
2 =α
Del triángulo P1 M O se tiene:

sen 2θ = ρP+H
1M


α 1

P1 M = (ρα + H1 ) sen 2θ

P1 L = 2P0 M

Utilizando el teorema de los senos en el triángulo P1 P2 L

sen β sen α
P1 L = H2 −H1
sen α
H2 − H1 = P1 L( sen β )

β = 90 − (α+ 2θ )
◦

◦ ◦ ◦
sen(90 − x) = sen(90 ) cos x − sen x cos(90 )
◦
sen β = sen[90 − (α + 2θ )]

= cos[α + 2θ ]
Z2 −Z1


θ

sen 2
H2 − H1 = 2(ρα + H1 ) sen 2
 Z2 −Z1


cos 2 + θ2

Desarrollando el denominador y dividiendo por cos (Z2 −Z

2
1)
cos 2θ se tiene que:

Z2 −Z1
θ
sen sen 2 2



Z2 −Z1
cos θ cos
H2 − H1 = 2(ρα + H1 )
2
2



Z2 −Z1 θ Z2 −Z1 θ
cos 2 cos sen sen

2
− 2

2


Z2 −Z1 θ Z2 −Z1 θ
cos 2 cos 2 cos 2 cos 2

20
 Z2 −Z1
tan( θ2 ) tan(
2 )
H2 − H1 = 2(ρα + H1 ) Z2 −Z1
1−tan( 2 ) tan( θ2 )

Desarrollando en series tan( 2θ ) y tomando sólo los dos primeros términos, se

tiene:
3
tan( 2θ ) = 2θ + 24
θ

θ = ρS
α
3
tan( 2θ ) = 2ρS + 24ρ
S
α α

S = Distancia geodésica que se reduce al elipsoide

ρ = Distancia (metros)
θ = Ángulo
S 3 tan( Z2 −Z1 )
( 2ρSα )+ 24ρ α 2
H2 − H1 = 2(ρα + H1 )
S 3 ) tan( Z2 −Z1 )
1−( 2ρSα + 24ρ α 2
Aplicando el teorema del binomio para el denominador y tomando el primer
término
 2
  1+S tan Z2 −Z1

H2 − H1 = 2(ρα + H1 ) 2ρSα 1 + 12ρS
tan Z2 −Z


2
1
2 2 2ρα
α

   1+S tan Z2 −Z1


Z2 −Z1 H1 S2
H2 − H1 = S tan 2
1 + ρα 1 + 12ρ2 2ρα
2
α

Se definen:
H1
A=1+ ρα
Corrección por elevación de la estación conocida
Z2 −Z1
B = 1 + 2ρSα tan


2
Corrección por diferencia en elevación
S2
C =1+ 12ρ2α
Corrección por distancia entre estaciones

3.1.1. Reciproca
Partiendo de las anteriores definiciones se deduce la fórmula de la nivelación
trigonométrica reciproca
H2 − H1 = S tan Z2 −Z


2
1
ABC
Adicionalmente, se aplica la corrección por refracción del ángulo cenital
 
00 i−o
C = S sen(0◦ 00 1” )

Siendo:
S: Distancia entre estaciones
i: Altura instrumental
o: Altura de la torre

21
3.1.2. No Reciproca
Adicional a las anteriores deducciones generales, se define k como k = 0,5 − m
y m = 0,071.
Posteriormente, se tiene un nuevo valor para la corrección por diferencia de
elevación (mencionada anteriormente como B)
◦
S tan(90 +k−Z1 )
B1 = 1 + 2ρα

Obteniendo ası́ una fórmula diferente para la nivelación trigonométrica no

reciproca.
◦
H2 − H1 = S tan(90 + k − Z1 )AB1 C

Igualmente, se define una corrección por refracción al ángulo cenital, esta es:
00
C = −i − o

Con i como la altura instrumental y o la altura de la torre.

3.2. Trisección
Corresponde a una técnica cuyo propósito es determinar las coordenadas geodési-
cas de un punto, conocidas las coordenadas geodésicas de tres puntos.
Conocidos (X,Y) → (N ,E).
Además obtener H altura ortométrica y h altura elipsoidal.
Esto siempre y cuando A, B y C sean intervisibles.
De la figura se tiene que:

pˆ1 = p̂ − pˆ2
sen(pˆ1 )
sen(pˆ2 ) = k

sen(p̂−pˆ2 ) 1
sen(pˆ2 ) = sen(pˆ2 ) sen p̂ cos pˆ2 − cos p̂ = tan1 pˆ sen p̂ − cos p̂
2

k = tan1 pˆ sen p̂ − cos p̂

2

sen p̂
 
pˆ2 = tan−1 k+cos p̂

ahora como
sen pˆ1 sen Â
AB
= pB

sen pˆ2 sen Ĉ

BC
= pB

22
 Figura 10: Trisección

se concluye que

sen  sen Ĉ
sen pˆ1 AB = sen pˆ2 BC

entonces
sen pˆ1 sen  BC
sen pˆ2 = sen Ĉ AB

se puede definir

k = sen  BC
sen Ĉ AB
◦
y como pˆ1 = p̂ − pˆ2 , B̂1 = 180 − (Â − p̂)
sen(p̂−pˆ2 ) 1
sen p2 = tan pˆ2 sen p̂ − cos p̂

Ap
= AB
sen pˆ1 → Ap = sen B̂1
sen p̂ AB
sen B̂1

23
 αAp̂ = αAB − Â
luego con
Ap̂ sen αAp̂ = ∆EAp
Ap̂ cos αAp̂ = ∆NAp
Las coordenadas de p serán:
Np = NA + Ap sen αAp̂
Ep = EA + Ap cos αAp̂
Método 2
Nota: No se cumple cuando es un cuadrado
h i h i
EC −EA −1 EB −EA
â = tan−1 N C −NA
− tan NB −NA
h i h i
EA −EB −1 EC −EB
b̂ = tan−1 N A −NB
− tan NC −NB
h i h i
−1 EB −EC −1 EA −EC
ĉ = tan NB −NC
− tan NA −NC

1
K1 = cot â−cot α̂
1
K2 =
cot b̂−cot β̂
1
K3 = cot ĉ−cot γ̂

Ep = K1 EAK+K2 EB +K3 EC
+K +K
1 2 3

Np = K1 NAK+K2 NB +K3 NC
+K +K
1 2 3

Procedimiento general:
a) Calcular las coordenadas cartesianas locales, conocidas las coordenadas
geodésicas.
b) Medir direcciones horizontales
c) Medir direcciones verticales
d) Calcular las direcciones definitivas
e) Calcular los ángulos de la figura α, β, γ, τ , η, δ
f) Calcular los ángulos verticales definitivos
g) Calcular Np y Ep
h) Determinar la altura ortométrica

24
 Figura 11: Métodos de trisección

3.3. Bisección
Consiste en determinar Np , Ep conocidos dos vértices A y B de coordenadas
NA , NB y EA , EB , y el ángulo entre Ap y Bp.
Distancia AB
p
AB = (∆NAB )2 + (∆EAB )2

∆NAB = NB − NA

∆EAB = EB − EA

Para calcular el azimut de AB, procedemos a:

h i
∆EAB
αAB = tan−1 ∆N AB

25
 Figura 12: Ángulos a calcular en la trisección

Del triángulo ABP utilizando la ley de senos se tiene que:

sen β sen ρ sen α
AP
= AB
= PB

sen β
AP = sen ρ AB

αAP = αAB − α

sen αAP = ∆EAP ∴ cos αAP = ∆NAP

AP AP

∆NAP = AP cos αAP

∆EAP = AP sen αAP

Sabemos además que

◦
ρ = 180 − (α + β)

26
 Figura 13: Bisección

entonces
sen βAB
AP = sen(180◦ −(α+β))

Las coordenadas de p serán:

Np = NA + ∆NAP

Ep = EA + ∆EAP

También es posible determinar los errores por ángulo y distancia que se expe-
rimentan en la bisección.

1) Por ángulo
◦
ρ = 180 − (α + β)

 
∂ρ
= −1
∂α

27
 que es lo mismo que  
∂ρ
= −1
∂β

∂ p̂ 2 ∂ p̂ 2



∇p̂ 2 = ∇A 2 ∂ α̂ +∇B 2
∂ β̂

2) Por distancia

sen β
AP = AB sen ρ

sen β
AP = AB sen(180◦ −(α+β))

∂AP
2 ∂AP
2
∇AP 2 = ∇β̂ 2 +∇α̂ 2 ∂ α̂
∂ β̂

3.4. Diferencial
Ofrece mayor precisión y se conoce como nivelación geométrica, consiste en
determinar las diferencias de nivel que tienen los diferentes trayectos para ob-
tener las alturas ortométricas y los puntos de interés.

CARTERA DIRECCIONES HORIZONTALES

◦
Estación D/I PV ’ ” ” Promedio Ángulo
D 60 15 17 18 17.5
Antena 1 60◦ 15’ 18.25”
I 240 15 18 20 19
D 100 15 17 16 16.5
G1 Antena 2 100◦ 15’ 18.5”
I 280 15 20 21 20.5
D 120 16 17 20 18.5
Antena 3 120◦ 16’ 18”
I 300 16 20 15 17.5

28
 CARTERA DIRECCIONES VERTICALES
◦
Estación D/I PV ’ ” ” Promedio Ángulo
D 100 10 15 16 15.5
Antena 1 100◦ 10’ 14.75”
I 259 49 45 47 46
D 99 59 59 57 58
G1 Antena 2 101◦ 22’ 17”
I 257 15 23 25 24
D 101 15 15 16 15.5
Antena 3 101◦ 27’ 20.25”
I 258 20 32 38 35

Procedimiento:

a) Leer el horizontal G1-Antena 1 en directo, generar pasos lentos y leer nue-

vamente G1-Antena 1

b) Leer el vertical G1-Antena 1 en directo, generar pasos lentos y leer nueva-

mente G1-Antena 1

c) Repetir los pasos a y b, con los demás puntos a observar

d) Repetir los pasos a,b, y c, pero ahora con inversos

e) En la cartera de direcciones horizontales, calcular el promedio de los segun-

dos y doble segundos
◦
f) Calcular los ángulos horizontales con la fórmula ( D+I−180
2
)

g) En la cartera de ángulos verticales, promediar los segundos y doble segundos

◦
de las lecturas, y aplicar la fórmula ( 360 +D−I
2
) para el cálculo del ángulo
promedio

Estación PV V+ V- hi Cota
BN107 1.911518 102.061518 100.15
BN108 BN108 0.451718 0.911516 101.60172 101.150002
BN109 BN109 1.414516 1.111716 101.90452 100.490004
BN109A BN109A 1.717113 100.187407
BN110 BN110 0.371618 1.353516 100.922622 100.551004
BN111 BN111 1.717136 1.211617 101.428141 99.711005
BN111A BN111A 1.616113 99.812028
BN111B BN111B 1.371617 100.056524
BN112 BN112 1.454516 99.973625

29
Nivelación por el método de subes y bajas,

Estación Sube Baja Cota

BN107 100.15
BN108 1.000002 101.150002
BN109 0.659998 100.490004
BN109A 0.302597 100.187407
BN110 0.061 100.551004
BN111 0.839999 99.711005
BN111A 0.101023 99.812028
BN111B 0.345519 100.056524
BN112 0.26262 99.973625

También es posible determinar la distancia entre estaciones, el error de la mira

y el error de lectura.
E = error de Lectura = HSD − 2HM D + HID
K = error de la Mira = HM I − 301550 − HM D
D = Distancia = (HSD − HID)/1000

Siendo,
HSD Hilo superior derecho
HM D Hilo medio derecho
HID Hilo inferior derecho
HM I Hilo medio izquierdo

30
4. LONGITUD DE ARCO
4.1. Longitud de arco meridiano
Sea A0 un punto sobre la superficie meridiana que se mueve una cantidad ds de
longitud de arco hasta el punto A0 .Dicho incremento produce un crecimiento
ϕ , en la latitud geodesica ϕ

Figura 14: Longitud de arco meridiano

A partir de la figura se obtiene:

dϕ = ϕA − ϕB

d0 = ds
AA

Como:
sin dϕ
lı́m dϕ = 1
dϕ→0

Es decir sin dϕ ≈ dϕ entonces:

ds
sin dϕ = ρ → ds = ρdϕ
R
s = ρdϕ
R a(1−e2 )
s=
(1−e2 sin ϕ2 )3/2

31
 Z ϕ2
2 dϕ
s = a(1 − e ) (4.1)
ϕ1 (1 − e2 sin ϕ2 )3/2

Para el desarrollo de esta integrar hacemos una expansión en la serie binómica,

y se tomara una aproximación hasta el cuarto termino; lo anterior por que aun
que es integrable la función, no posee antiderivada, que es lo que nos interesa.

Recordemos antes, que es la serie binómica, y verla como un caso particular

de la Serie de Taylor.

Teorema de Bernstein Suponga que f y todas sus derivadas sean no nega-

tivas en un intervalo compacto [b, b + r]. Entonces si b < x < b + r, la Serie de
Taylor
∞
P f (n) (n)
n! (x − c)n
n=0

Converge hacia f (x)

A modo particular se puede obtener el desarrollo para la siguiente serie, cono-

cida como la serie binómica:
∞
(an) xn Si −1 < x < 1
P
(1 + x)a =
n=0

Donde a ∈ R, y además

a(a−1)...(a−n+1)
(an) = n!
El teorema de bernstein no es directamente aplicable en este caso; más aun se
puede argumentar que siendo

f (x) = (1 − x)−c

Donde c > 0yx < 1 .Entonces

f n (x) = c(c + 1)...(c + n − 1)(1 − x)−c−n

por lo tanto f n (x) ≥ 0 para cada n con tal que x < 1. aplicando el teorema de
bernstein con b = −1yr = 2, hallamos que f (x) tiene un desarrollo en serie de
potencias convergente para −1 < x < 1. por lo tanto f (x) tiene desarrollo en
serie de potencias en torno a 0,

32
 ∞
f (k) (0) n
n! x
P
f (x) =
n=0

k


k
Convergente para −1 < x < 1. Pero f (0) = −c (−1)k k!, luego

∞
1 k


=
P
(1+x)c −c (−1)k xk Si −1 < x < 1
n=0

sustituyendo c por −a y −x, encontramos que la ecuación anterior es valida

también para a < 0 y por integraciones sucesivas a todo a ∈ R. Es evidente
que si a ∈ R, por ejemplo a = m, entonces (nm ) = 0 para m > n. y la ecuación
anterior se reduce a una suma finita (Teorema del binomio).

ϕ2
a(1 − e2 )
Z
s=  3 · dϕ (4.2)
1 − e2 sin2 ϕ 2

ϕ1

− 3
1 − e2 sin2 ϕ 2 = 1 − 32 (−e2 sin2 ϕ) + − 32 − 52 − 2!
1




(−e2 sin2 ϕ)2


+ − 23 − 52 − 72 1





3! (−e2 sin2 ϕ)3 + ...

1+ 23 e2 sin2 ϕ + 15 4 35 6 6
8 e sin ϕ + 16 e sin ϕ + ...
4

Entonces, reemplazando en (4.2)

Z ϕ2
2
n 3 15 35 o
s = a(1 − e ) 1 + e2 sin2 ϕ + e4 sin4 ϕ + e6 sin6 ϕ + ... dϕ (4.3)
ϕ1 2 8 16

Z ϕ2
2
n 3 15 35 o
s ' a(1 − e ) 1 + e2 sin2 ϕ + e4 sin4 ϕ + e6 sin6 ϕ dϕ (4.4)
ϕ1 2 8 16

para el calculo de tal integral haremos, un breve repaso, y expresaremos poten-

cias de sin ϕ, como combinación lineal de sin nϕ y cos nϕ, para n = 1, 2, 3, 4, 5, 6, ...
; empezaremos con:

cos (a + b) = cos a cos b − sin a sin b (4.5)

cos (a − b) = cos a cos b + sin a sin b (4.6)

33
Sumando (4.5) y (4.6), obtenemos:
1
cos a cos b = (cos (a + b) + cos (a − b)) (4.7)
2
1
sin2 ϕ = (1 − cos 2ϕ) (4.8)
2
1
cos2 ϕ = (1 + cos 2ϕ) (4.9)
2
 
1 1
cos3 ϕ = cos ϕ (1 + cos 2ϕ) = {cos 2ϕ cos ϕ + cos ϕ} (4.10)
2 2
1 1
cos 2ϕ cos ϕ = cos 3ϕ + cos ϕ (4.11)
2 2
 2
1 1
cos ϕ = (1 + cos 2ϕ) = (1 + 2 cos 2ϕ + cos2 2ϕ)
4
2 4
1 1
= (1 + 2 cos 2ϕ + (1 + cos4ϕ))
4 2
1 1 1
= (1 + 2 cos 2ϕ + + + cos4ϕ)
4 2 2
3 1 1
= + cos 2ϕ + cos4ϕ)
8 2 8
sin ϕ = (1 − cos ϕ) = 1 − 2 cos2 ϕ + cos4 ϕ
4 2 2

   2
1 1
= 1 − 2 (1 + cos 2ϕ) + (1 + cos 2ϕ)
2 2
1 1 1 1 1 1
= + cos 2ϕ + cos2 2ϕ − cos 2ϕ = − cos 2ϕ + (1 + cos 4ϕ)
4 2 4 4 2 8
3 1 1
= − cos 2ϕ + cos 4ϕ
8 2 8
 3  
6 2 3 1 1 3 1 1
sin ϕ = (sin ϕ) = (1 − cos 2ϕ) = (1−cos 2ϕ) − cos 2ϕ + cos 4ϕ
2 2 8 2 8
 
1 3 1 1
= (1 − cos 2ϕ) − cos 2ϕ + cos 4ϕ
2 8 2 8
3 1 1 3 1 1
= − cos 2ϕ + cos 4ϕ − cos 2ϕ + cos 2ϕ cos 2ϕ − cos 4ϕ cos 2ϕ
16 4 16 16 4 16
3 7 1 1 1
= − cos 2ϕ + cos 4ϕ + (cos 4ϕ + 1) − (cos 6ϕ + cos 2ϕ)
16 16 16 8 32
5 15 3 1
= − cos 2ϕ + cos 4ϕ − cos 6ϕ
16 32 16 32

34
Por lo tanto,

3 15 35 3 45 175 6 n 3 2 15 4 525 6 o
1+ e2 sin2 ϕ+ e4 sin4 ϕ+ e6 sin6 ϕ = 1+ e2 + e4 + e− e+ e+ e
2 8 16 4 64 256 4 16 512
n 15 105 6 o 35 6
cos 2ϕ + e4 + e cos 4ϕ − e cos 6ϕ
64 256 512
Definiendo:
3 45 175 6
A = 1 + e2 + e4 + e (4.12)
4 64 256
3 15 525 6
B = e2 + e4 + e (4.13)
4 16 512
15 105 6
C = e4 + e (4.14)
64 256
35 6
D= e (4.15)
512
y reemplazando en (4.4) se tendrá:
Z ϕ2
2
s ' a(1 − e ) {A − B cos 2ϕ + C cos 4ϕ − D cos 6ϕ}dϕ (4.16)
ϕ1
 ϕ=ϕ2
2 1 1 1
s ' a(1 − e ) Aϕ − B cos 2ϕ + C cos 4ϕ − D cos 6ϕ (4.17)
2 4 6 ϕ=ϕ1

4.2. Longitud de arco paralelo

La longitud de arco L, sobre un mismo paralelo corresponde a un arco de
circunferencia (Elipsoide biaxial de revolucion sobre el eje z) Delimitado por
los meridianos que pasan por dos puntos cualquiera (RyS) que tienen en comun
la misma latitud. En el triangulo y la elipse de la figura, vemos que: 4λ =
dλ = λR − λS ;
Q
cos ϕ =→ Q = N cos ϕ
N
L L
sin dλ ' dλ → sin dλ = → dλ = → L = Qdλ
Q Q
L = N cos ϕdϕ
L = N cos ϕ 4 ϕ

35
 Figura 15: Longitud de arco paralelo

4.3. Área de un cuadrilátero sobre el elipsoide

Dado el cuadrilátero infinitesimal ABCD sobre el elipsoide, donde AD = BC

AB = DC = L = N cos ϕ 4 λ
A = ρ sin ϕ ' dϕ
Z λ2 Z ϕ2
A= N ρ cos ϕdϕdλ
λ1 ϕ1
entonces
 ϕ=ϕ2
2 1 sin ϕ 1 e sin ϕ + 1
A = b (λ2 − λ1 ) 2 + ln
2 1 − e sin ϕ 4e e sin ϕ − 1
2
ϕ=ϕ1
La situación anterior se aprecia en la siguiente figura,

5. RADIOS DE CURVATURA
Plano normal, es un plano que corta el elipsoide y contiene a la normal.
Sección normal, es una curva generada por el corte o intersección del elip-
soide con un plano normal.
Sección meridiana, es aquella que contiene un punto (M ) y los polos, es
decir, esta sección normal corresponde a una elipse que tiene los puntos P M
P”.

36
 Figura 16: Área de un cuadrilátero

Primer vertical, sección normal, perpendicular a la sección meridiana que

contiene al punto en estudio.
La curvatura (K ), se define como:
 3
dy 2 2
1 + ( dx )
K= d2 y
dx2

Para hallar tal valor, partimos de las expresiones para x y z. Y hallamos las
derivadas necesarias:
a cos ϕ a(1−e2 ) sin ϕ
x =√ y z =√
1−e2 sin2 ϕ 1−e2 sin2 ϕ

La curvatura, es la variación de la dirección del vector tangente sobre la elipse

meridiana. De la figura tenemos:
−dx dz
tan ϕ = → = − cot ϕ
dz dx

5.1. Radio de curvatura de la sección meridiana (ρm )

h  dy 2 i 23
1 + dx
K= d2 y
dx2

37
Para este cálculo, se usa la regla de la cadena, es decir:
dz dz dϕ
=
dx dϕ dx
dx a(1 − e2 ) sin ϕ
= 3
dϕ [1 − e2 sin2 ϕ] 2
dz a(1 − e2 ) cos ϕ
= 3
dϕ [1 − e2 sin2 ϕ] 2
d2 z d dϕ
2
= (− cot ϕ) = csc2 ϕ
dx dϕ dx
3 3
1 − e2 sin2 ϕ 2 1 − e2 sin2 ϕ 2
  
1
= =
sin2 ϕ a(1 − e2 ) sin ϕ a(1 − e2 ) sin3 ϕ
3  2  32
csc ϕ a(1 − e2 ) sin3 ϕ

1 + [− cot]2 2
ρm = 3 = 3
[1 − e2 sin2 ϕ 2 1 − e2 sin2 ϕ 2


a(1 − e2 ) sin3 ϕ
a(1 − e2 ) csc3 ϕ sin3 ϕ
= 3
1 − e2 sin2 ϕ 2


a(1 − e2 )
ρm =  3
1 − e2 sin2 ϕ 2
Esta fórmula también se puede deducir de otra forma, a partir de la figura 4,
tenemos:
r  dx 2
p p
2 2
s = (dx) + (dz) = (1 + dz = 1 + (− tan ϕ)2 dz
dz
p dz
= sec2 ϕdz = sec ϕdz =
cos ϕ
Sabiendo que sin dϕ y dϕ tienen un comportamiento asintotico alrededor de
cero, porque:

sin dϕ
lı́m dϕ = 1
dϕ→0

tenemos,
s s
sen dϕ = → dϕ ≈
ρm ρm
s 1
ρm = =s
dϕ dϕ

38
 h i
dz 1 dz 1 a(1 − e2 ) cos ϕ 1
ρm = cos ϕ dϕ = dϕ cos ϕ = 3
cos ϕ
[1 − e2 sin2 ϕ] 2

a(1 − e2 ) a(1 − e2 )
ρm = 3 =
[1 − e2 sin2 ϕ] 2 w3

Algunos resultados importantes de esta sección son:

a
lı́mπ N =
ϕ→ 1−f
2
a
lı́mπ ρm =
ϕ→ (1 − f )3
2

lı́m ρm = a(1 − e2 )
ϕ→0

lı́m N = a
ϕ→0

5.2. Radio de curvatura de un paralelo

Ahora para un punto (p), tenemos ρp el radio de curvatura de un paralelo. Es
evidente que para cualquier punto sobre la superficie del elipsoide:

ρp = xp = N cos ϕ

5.3. Radio de curvatura en el primer vertical

Sobre la sección paralela no se tiene la Normal. El radio de curvatura de la
sección normal en una superficie cualquiera es el radio de curvatura de la
sección oblicua que tenga la misma recta tangente dividida por el coseno del
ángulo que forman ambas secciones.
1
ρN = N cos ϕ =N
cos ϕ
Al hacer la relación:
a
p
N 1 − e2 sin2 ϕ 1 − e2 sin2 ϕ 1 − e2 + e2 − e2 sin2 ϕ
= = =
ρm a(1 − e2 ) 1 − e2 1 − e2
3
1 − e2 sin2 ϕ 2


1 − e2 + e2 (1 − sin2 ϕ) 1 − e2 + e2 cos2 ϕ e2 cos2 ϕ

= = = 1 +
1 − e2 1 − e2 1 − e2

39
Del resultado anterior, dado que:

N e2 cos2 ϕ
=1+ >1
ρm 1 − e2
concluimos que N > ρm siempre se cumple.

5.4. Radio de curvatura de una sección normal en un

azimut cualquiera
Es el radio de curvatura en cualquier dirección (α) con respecto a un observa-
dor.
Para hallar el radio de curvatura (r.c.) de una sección Normal en un azimut
cualquiera. se usa la fórmula de Euler, ya que usualmente N > M .

1 cos2 α sin2 α
= +
ρα ρ1 ρ2
Donde:
ρα : r.c. de una sección Normal en un azimut cualquiera.
ρ1 : r.c. de la sección meridiana ρm
ρ2 : r.c. de la sección Normal. Primer vertical N

1 cos2 α sin2 α cos2 α sin2 α

= + = +
ρα ρ1 ρ2 ρm N
N cos2 α + ρm sin2 α N cos2 α + M sin2 α
⇒ ≡
N ρm MN
MN
ρα =
N cos2 α + M sin2 α

5.5. Valor medio curvatura

En un punto de la superficie, se toma como el limite de la media aritmética de
los radios de curvatura de las secciones normales, cuando el número de estos
tiende a infinito se denota ρα .
A partir del teorema del valor medio para integrales se puede hallar dicho
valor. Este teorema se enuncia como:
Teorema del valor medio (para integrales)
Z b
f (x)dx = f (c)(b − a)
a
Dicho de otra manera, buscamos un valor x, que tenga la siguiente propiedad:
Z b
1
x= f (x)dx
b−a a

40
Para nuestro caso a = 0 y b = 2π y
N ρm
f (α) =
N cos2 α + ρm sin2 α
Z 2π
1 N ρm
ρα = dα
2π 0 N cos2 α + ρm sin2 α
Z π
2 2 1
= N ρm 2 2 dα
π 0 N cos α + ρm sin α

Z π 1
2 2
N cos2 α dα
= N ρm
π 0 ρm sin2 α
1+
N cos2 α
1
2
Z π
2 sec2 α
= N ρm N dα
π ρm
0 1+ tan2 α
N
Z π
2 2 sec2 α
= ρm 2 dα
π
p ρm
0 1+ tan α
N

Sustituyendo,
r s
ρm N
u= tan α → du = sec2 αdα
N ρm
Ası́
π
α=0→u=0 y α= =⇒ u → ∞
2
s Z ∞
2 N du
ρα = ρm
π ρ m 0 1 + u2
Z ∞
2 p du
ρα = ρm N ρ m
π 0 1 + u2
2 p u=t
= ρm N ρm · lı́m [tan−1 u]
π t→∞ u=0

2 p hπ i
ρα = ρm N ρm
π 2
p
ρα = N ρ m

√
ρm = M =⇒ ρα = NM

41
Se ha usado que:
π
lı́m tan−1 ϕ =
ϕ→∞ 2
Se puede definir otro radio para el elipsoide, el radio medio (Rm ) como:
a+b+c
Rm =
3
√
tomando a = c y b = a 1 − e2
√
a+b+c 2 a 1 − e2
Rm = = a+
3 3 3
Ejercicio
Mostrar que,
Z
cos ϕ 1 sin ϕ 1 e sin ϕ + 1
2 dϕ = 2 + ln
2
(1 − e sin ϕ) 2 2
2 1 − e sin 4e e sin ϕ − 1
Demostración: si se sustituye x = sin ϕ → dx = cos ϕdϕ se tiene que:

(1 − e2 sin2 )2 = (1 − e2 x2 )2 = 1 − 2e2 x2 + e4 x4
además,

1 1 1 1 1
= + + +
1− 2e2 x2 + e 4 x4 4(ex − 1) 4(ex − 1) 2 4(ex + 1) 4(ex + 1)2
Por tanto,
Z Z
cos ϕ dx
dϕ =
(1 − e2 sin2 ϕ)2 1− 2e2 x2 + e 4 x4

−dx
Z Z Z Z Z
dx dx dx dx
= + + +
1 − 2e x2 + e4 x4
2 4(ex − 1) 4(ex − 1) 2 4(ex + 1) 4(ex + 1)2
1 1 1 1 1 1
= − ln |ex − 1| − + ln |ex + 1| −
4 4e ex − 1  4e 4e ex + 1
1 ex + 1 2ex
= ln − 2 2
4e ex − 1 e x −1
Z  
cos ϕ 1 e sin ϕ + 1 2e sin ϕ
dϕ = ln − 2 2
(1 − e2 sin2 ϕ)2 4e e sin ϕ − 1 e sin ϕ − 1
1 sin ϕ 1 e sin ϕ + 1
= 2 + ln
2 1 − e2 sin ϕ 4e e sin ϕ − 1
Como caso particular,
ϕ2 h1 sin ϕ 1 e sin ϕ + 1 iϕ=ϕ2
Z
cos ϕ
dϕ = + ln
ϕ1 (1 − e2 sin2 ϕ)2 2 1 − e2 sin2 ϕ 4e e sin ϕ − 1 ϕ=ϕ1

42
6. PROBLEMA GEODÉSICO
6.1. Secciones Normales
Considerando que el elipsoide es una figura que tiene aplanamiento, todos
los puntos que lo conforman no son equidistantes al origen, por tal razón es
necesario definir las secciones normales, las cuales resultan de la intersección
entre el elipsoide y un plano que contiene la normal N .
De las infinitas secciones normales que se obtienen al realizar las intersecciones,
existen dos secciones normales principales, llamadas la sección normal (primer
vertical) y la sección meridiana; son principales porque tienen el mayor (N ) y
menor radio de curvatura (M ).
La sección meridiana contiene los polos, el punto T y el radio N; mientras que
el primer vertical es perpendicular a la sección anterior y contiene el punto T
y el radio N.
Una sección normal especı́fica desde el punto A al punto B está formada por
la intersección de un plano, que contiene la normal en el punto A que pasa a
través del punto B, con la superficie del elipsoide de referencia.
La sección normal A a B contiene los puntos A y B la Normal asociada a A y
no contiene NB , mientras que la sección B a A contiene los puntos A y B, NB ,
y no contiene NA .

Figura 17: Secciones normales

De la figura se ven las secciones normales de A a B y de B a A; dichas secciones

son diferentes por que la intersección de normales a diferentes latitudes con
el semieje menor no es la misma. La distancia OnA y OnB se puede calcular ası́:

Del triángulo NA AJ se tiene que:

43
 Figura 18:

AJ
sin ϕA = =⇒ AJ = NA sin ϕA
NA
Por definición ZA = NA (1 − e2 ) sin ϕA

OnA = AJ − ZA

OnA = NA sin ϕA − NA sin ϕA + e2 NA sin ϕA

OnA = e2 NA sin ϕA
Por analogı́a,
OnB = e2 NB sin ϕB
Si ϕA > ϕB =⇒ OnA > OnB .
Mientras más septentrional este el punto, más grande Onn .

Onn = f (ϕn )

entonces si, ϕA < ϕB es decir ϕA 6= ϕB se tiene que

OnA 6= OnB
Demostrándose que las secciones normales reciprocas no son coincidentes.

44
 Figura 19: Linea geodésica - Plano Osculante

En general, tener dos secciones normales entre dos puntos crea problemas en
los cómputos. Puede verse en la figura, que donde las lı́neas observadas se
indican para un triángulo en la superficie del elipsoide. Los ángulos medidos
son α1 , α2 , α3 ; para esta figura no se ha observado ninguna figura cerrada.

Figura 20:

45
6.2. Lı́nea Geodésica
Se define lı́nea geodésica como aquella lı́nea en la cual los infinitos puntos que
la componen son normales al plano tangente asociado a los puntos. Es además
la lı́nea que representa la distancia más corta entre dos puntos del espacio sin
importar su curvatura.
Ası́ por ejemplo:

En una geometrı́a euclidiana → Lı́nea recta (segmento)

Esfera (geometrı́a no euclidiana) → Arco de cı́rculo máximo

Elipsoide → Lı́nea geodésica, proyección en superficie del plano osculan-

te.

La lı́nea geodésica contiene el punto A, el punto B, la NA , la NB y las infinitas

normales de los infinitos puntos que la definen.
Plano Osculante: También llamado de contacto, es el plano intersección de
la linea geodésica y el elipsoide.

Figura 21: La lı́nea geodésica

N
tan ϕ =
PQ
P Q = N (tan ϕ)
De la figura:
ρm dϕ N cos ϕdλ N cos ϕdλ
cos α = sin α = tan α =
ds ds ρm dϕ

46
 N cos ϕdλ
dα =
N cot ϕ

ds2 = (ρm dϕ)2 + (N cos ϕdλ)2

De la anterior relación

dsdϕ2 = ρ2 + N 2 cos2 dϕdλ

Pero:
dλ ρ tan α
=
dϕ N cos ϕ
ds2
2
= ρ2 + ρ2 tan2 α =⇒ ds2 = ρ2 dϕ2 + ρ2 tan2 αdϕ2
dϕ
p
ds = ρdϕ 1 + tan2 α
Con lo que se obtiene que la magnitud de la lı́nea geodésica (longitud), es una
función de la latitud. Después del proceso algebraico en el que se debe obtener
el valor de S con la condición que sea mı́nima. Se obtiene:

N cos ϕ sin α = cte

Ecuación de Clarrut:
Comprobación,
1
n cos ϕ = pp = a cos ϕ(1 − e2 sin2 ϕ)− 2

dp −a sin ϕ 3
=p = −a cos ϕ12(1 − e2 sin2 ϕ)− 2 (−2e2 sin ϕ cos ϕ)
dϕ 2
1 − e sin ϕ2

−a sin ϕ 2 2 2 2
= 2 3 [1(1 − e sin ϕ) − cos ϕe ]
2
(1 − e sin ϕ) 2

a sin ϕ 2 2 2
=− 2 3 [1 − e (sin ϕ + cos ϕ)]
2
(1 − e sin ϕ) 2
a(1 − e2 ) sin ϕ
=− 3
(1 − e2 sin2 ϕ) 2
= −ρm sin ϕdϕ

De la ecuación tenemos que:

cos ϕ dα
dα = dλ = sin ϕdλ → sin ϕ =
cot ϕ dλ

47
Reemplazando en la ecuación anterior tenemos que:
dα
dp = −ρm dϕ
dλ
Pero
dϕ ds cos α
cosα = ρm =⇒ dϕ =
ds ρm
entonces,
dp = −ρm dαds cos αρm dλ
dαds cos α N cos ϕ pdλ
dp = − ds = dλ =
dλ sin α sin α
dα pdλ dp
dp = − cos α =⇒ dα = − tan α
dλ sin α p
p cos αdα + sin αdp = 0
Lo que implica que:

p1 sin α1 = p2 sin α2 = pn sin αn = K

p sin α = N cos ϕ sin α = K

◦
K es una constante. Para casos extremos, el Ecuador ϕ = 0, N = a, α = 90 ,
se tiene el máximo azimut.

48
6.3. Método Puissant
6.3.1. Método directo

Figura 22:

Objetivo: el problema geodésico directo busca resolver a partir de los datos

conocidos de un punto P cualquiera; los datos desconocidos asociados a este;
para otro punto cualquiera Q.
es decir conocidos ϕp , λP , αP Q , SP Q . Determinar: ϕQ , λQ , yαQ .

49
 Figura 23:

En el triangulo esférico anterior P1 tiene coordenadas (ϕ01 , λ01 ); P2 tiene coor-

denadas (ϕ01 , λ01 ).

Este problema ha sido resuelto asumiendo diferentes condiciones iniciales y una

de las mas sencillas es la de puissant. para deducir las ecuaciones necesarias
del problema directo. consideremos una esfera de radio N1 ; es decir, se hace
coincidir el elipsoide y la esfera en el punto P1 (como se observa en las figu-
ras); significa que este punto es tangente a lo largo del paralelo que pasa por la
latitud ϕ1 . en distancias cortas la esfera estará aproximadamente coincidente
con el punto 2: ademas se supone que el acimut y las distancias son iguales en
la esfera y el elipsoide.
P1 = 90−ϕ01 , Pd
Del triangulo esférico P1 P2 P , se tienen los arcos Pd P2 = 90−ϕ02 ;
Al medir estos arcos se tiene que ϕ1 = ϕ01 , aplicando la ley de coseno se tiene:

cos (90 − ϕ02 ) = cos (90 − ϕ01 ) cos P

[ 0
1 P2 + sin 90 − ϕ1 sin P1 P2 cos α12
[ (6.1)

sin ϕ02 = sin ϕ01 cos P

[ 0
1 P2 + cos ϕ1 sin P1 P2 cos α12
[ (6.2)
Como el problema resuelve con precision las lineas geodésicas de cortas dis-
tancias se tiene la siguiente relación:

50
 Figura 24:

ϕ02 = ϕ1 + 4ϕ

Donde 4ϕ es una cantidad pequeña en radianes; ademas para resolver el pro-

blema se debe hacer la relacion con el radio de la esfera N1 asi:

S12
el arco P
[ 1 P2 = . lógicamente S12 = S21 ; de tal forma que la ecuación (6.2)
N1
se convierte en:

   
S S
sin (4ϕ + ϕ1 ) = sin ϕ01 cos + cos ϕ01 cos α12 sin (6.3)
N1 N1
Como ϕ1 = ϕ01 entonces se tiene que:
   
S S
sin (4ϕ + ϕ1 ) = sin ϕ1 cos + cos ϕ1 cos α12 sin (6.4)
N1 N1
El problema se debe seguir desarrollando mediante una serie de Taylor para
sin 4ϕ + ϕ1 .

recordando que:

Dada una función f (x) según la definición de la serie se tiene que:

51
 (x − x0 )2 00 (x − x0 )3 000
f (x) = f (x0 ) + (x − x0 )f 0 (x0 ) + f (x0 ) + f (x0 ) + ...
2! 3!
∞
X (x − x0 )n f (n) (x0 )
f (x) =
n=0
k!
Si por ejemplo se tiene que:
f (x) = sin xy(x − x0 ) = hy(x = x0 )entoncessegunlaseriesetieneque.

h2 h3 h4
sin (x + h) = sinx + h cos x − sin x − cos x + sin x + ...
2 3! 4!
como se quiere desarrollar sin(ϕ1 + 4ϕ0 ) se tiene que:

4ϕ02 4ϕ02 4ϕ04

sin (ϕ1 + 4ϕ0 ) = sin ϕ1 +4ϕ0 cos ϕ1 − sin ϕ1 − cos ϕ1 + sin ϕ1 +...
2 3! 4!
(6.5)
S
Al considerar que es muy pequeño se tiene que:
N1

S S3
sin SN1 ' − (6.6)
N1 6N12
S2
cos SN1 ' 1 − (6.7)
2N12
AL reemplazar (6.6) y (6.7) en (6.4) y hacer la igualdad con la ecuacion (4.5)
se tiene que:

4ϕ02 4ϕ02 S2 S3
   
0 S
sin ϕ1 +4ϕ cos ϕ1 − sin ϕ1 − = sin ϕ1 1 − +cos ϕ1 − cos α12
2 3! 2N12 N1 6N12
(6.8)
Al despejar 4ϕ1 de la ecuación (6.8) se tiene:

S2 S3
   
sin ϕ1 cos ϕ1 S sin ϕ1 sin ϕ1
4ϕ1 = 1− 2
+ − 2
cos α12 − +
cos ϕ1 2N1 cos ϕ1 N1 6N1 cos ϕ1 cos ϕ1

4ϕ02 cos ϕ1 4ϕ03 S2 S3

   
S
= tan ϕ1 1 − + cos α12 − cos α12
2 cos ϕ1 6 2N12 N1 6N13

52
 4ϕ02 4ϕ02
   
− tan ϕ1 + tan ϕ1 + (6.9)
2 3!

S2 S3 4ϕ02 4ϕ03
 
0 S
4ϕ = cos α12− tan ϕ 1 + cos α12 + tan ϕ 1 + (6.10)
N1 2N12 6N13 2 6

Como se nota en la ecuación (6.10) 4ϕ0 aparece a los dos lados de la relacion,
el problema se debe resolver por aproximaciones sucesivas. si se considera
S cos α12
4ϕ0 =
N1
como la primera aproximación

S2 S3 S2 4ϕ03
 
0 S 2
4ϕ = cos α12 − tan ϕ 1 + cos α12 + tan ϕ1 cos α12 +
N1 2N12 6N13 2N12 6

S2 S3 S2 4ϕ03
 
S 2
= cos α12 − tan ϕ1 − cos α 12 + tan ϕ1 (1 − sin α12 ) +
N1 2N12 6N13 2N12 6

S3 S2 4ϕ03
 
S 2
= cos α12 − cos α 12 − tan ϕ1 sin α 12 + (6.11)
N1 6N13 2N12 6
De la ecuación anterior se nota que una mejor aproximación a 4ϕ0 es:

S2
 
0 S
4ϕ = cos α12 − tan ϕ1 sin2α1 2 (6.12)
N1 2N12
Sustituyendo la ecuación (6.12) en la ecuación (6.10) se tiene:

S2 S3
 
0 S
4ϕ = cos α12 − tan ϕ 1 − cos α12
N1 2N12 6N13
2
S2 tan ϕ1 4ϕ03

S 2
+ cos α12 − tan ϕ sin α 12 + (6.13)
N1 2N12 2 6
S2 S3
 
0 S
4ϕ = cos α12 − tan ϕ1 − cos α12
N1 2N12 6N13
S2 2 S2 2 S3
+ cos α 12 − cos α12 tan ϕ 1 − tan2 ϕ1 cos α12 sin2 α12
2N12 2N12 2N13
S4 3 4 4ϕ03
+ tan ϕ1 sin α 12 + (6.14)
N14 6
Despreciando el termino de la cuarta potencia

53
 S2 S3
 
0 S
4ϕ = cos α12 − tan ϕ1 − cos α12
N1 2N12 6N13

S2 2 S3 2 2 4ϕ03
− (1 − sin α 12 tan ϕ1 ) − tan ϕ1 cos α 12 sin α 12 + (6.15)
2N12 2N13 6

De la ecuación (6.10) se tiene:

4ϕ03 S3
= cos α12 (6.16)
6 6N13
Reemplazando la ecuación (6.16) en la ecuación (6.11)

S2 S3
 
0 S
4ϕ = cos α12 − tan ϕ1 − cos α12
N1 2N12 6N13

S2 2 S3 2 2 S3
+ (1 − sin α 12 tan ϕ1 ) − tan ϕ1 cos α 12 sin α 12 + cos α12 (6.17)
2N12 2N13 2N13

Si se considera cos3 α12 = cos α12 (1 − sin2 α12 ) se obtiene:

S S2
4ϕ0 = cos α12 − sin2 α12 tan ϕ1
N1 2N12
S3
− cos α12 sin2 α12 [1 + 3 tan2 α12 ] (6.18)
6N13
Recordar que 4ϕ0 seconsideroparalaesf eraderadioN1 ; sim embargo se debe
cambiar para 4ϕ medido a lo largo de un arco de meridiano. para hacerlo se
asume que la distancia N1 4ϕ0 en la esfera es igual que la distancia correspon-
diente en el elipsoide. Considerando que ρm (ϕm ) sea la longitud del radio de
curvatura en la latitud medida se tiene que:

N1 4ϕ0 = ρm (ϕm )4ϕ (6.19)

De la ecuación anterior (6.19), notese que no es posible determinar ρm (ϕm ) sin
lo cual, tampoco se puede determinar 4ϕ que es el valor buscado. entonces
aunque no se tenga la latitud del segundo punto ϕ2 ; se resuelve el problema
de hallar ρm mediante un desarrollo de ρm (ϕm ) alrededor del primer punto.
∂ρm 4ϕ
ρm (ϕm ) = ρm ϕ1 + + ... (6.20)
∂ϕm ϕ=ϕm 2

Al efectuar la diferenciación se tiene:

3  ρm ϕm e2 sin ϕ1 cos ϕ1 
ρm ϕ m = ρm ϕ 1 + 4ϕ (6.21)
2 1 − e2 sin2 ϕ1

54
Recordar que:

a(1 − e2 )
ρm ϕm =
(1 −e2 sin2 ϕ1 )3/2
Diferenciando con respecto a ϕ, y evaluando en ϕ1 se tiene:

∂ρm 3ae2 (1 − e2 ) sin ϕ cos ϕ

=
∂ϕ (1 −e2 sin2 ϕ)5/2
∂ρm 3ae2 (1 − e2 ) sin ϕ1 cos ϕ1
= (6.22)
∂ϕ ϕ=ϕ1 (1 −e2 sin2 ϕ1 )5/2
Resolviendo la ecuación (6.19) para 4ϕ y sustituyendo la ecuación (6.21) y
(6.22) en esta expresión, se tiene:

N1 4ϕ0 = ρm ϕm 4ϕ

3 e2 sin ϕ cos ϕ
N1 4ϕ0 = ρm ϕ1 + ρm ϕ1
 
2 4ϕ 4ϕ (6.23)
2 (1 − e2 sin ϕ)
Usando la ecuación (6.18) y haciendo ρm (ϕ1 ) = ρ se tiene para la ecuación
(6.23) que:
 3e2 sin ϕ1 cos ϕ 
4ϕ + 2 4ϕ 4ϕ
2(1 − e2 sin ϕ1 )
S3
− cos α12 sin2 α12 (1 + 3 tan2 ϕ1 ) (6.24)
6N12 ρ

S S2 S3
4ϕ = cos α12 − sin2 α12 tan ϕ1 − cos α12 sin2 α1 2(1 + 3 tan2 ϕ1 )
ρ 2N1 ρ 6N12

3e2 sin ϕ1 cos ϕ1

− 4ϕ2 (6.25)
2(1 − e2 sin2 ϕ1 )
Es decir:

4ϕ = δϕ − C(4ϕ)2 (6.26)
Donde

S S2
δϕ = cos α12 − sin2 α12 tan ϕ1
ρ 2N1 ρ
S3
− cos α12 sin2 α12 (1 + 3 tan2 ϕ1 ) (6.27)
6N12 ρ
y

55
 3e2 sin ϕ1 cos ϕ1
C= (6.28)
2(1 − e2 sin2 ϕ1 )
Con la consideración que (δϕ − 4ϕ) es pequeño, se puede dejar a (4ϕ)2 que
aparece en la ecuación e introduciendo las siguientes constantes determinadas
por Hosmer se tienen:
1
B= (6.29)
ρm (ϕ1 )
tan ϕ1
C= (6.30)
2ρm (ϕ1 )N1
3e2 sin ϕ1 cos ϕ1
D= (6.31)
2(1 − e2 sin2 ϕ1 )
1 + 3 tan2 ϕ1
E= (6.32)
6N12
S cos α12
h= (6.33)
ρm (ϕ1 )
Se tendrá que:

4ϕ” = S cos α12 B − S 2 sin2 α12 C − hS 2 sin2 α12 E − (δϕ)2 D (6.34)

Donde δϕ resulta de la ecuación (6.27) y/o por la suma de los tres primeros
términos. Entonces la latitud del segundo punto sera:

ϕ2 = ϕ1 + 4ϕ” (6.35)
4ϕ Es un valor en segundos.

6.3.2. Método inverso

el objetivo del problema geodésico inverso es determinar la distancia entre dos
puntos de los que se conocen las coordenadas geográficas, elipsoidales y los
azimuts entre puntos. Se usa actualmente para calcular los datos conocidos
(α12 , α21 , S) entre dos puntos georeferenciados con G.P.S los datos se obtienen
a partir del problema geodésico directo y como ya se menciono el objetivo es:

Conocidos: ϕ1 , λ1 , ϕ2 , λ2

Determinar: α1 , α2 , S

56
 S S2
(1 − sin2 α12 sec2 ϕ2 )
 
4λ = sin α12 sec ϕ2 1 − 2
N2 6N2
Despejando se tiene:
N2 4λ cos ϕ2
S sin α12 = S2
1− 6N22
(1 − sin2 α12 sec2 ϕ2 )
Ahora también se sabe que:

4ϕ” = S cos α12 B − S 2 sin2 α12 C − hS 2 sin2 α12 E − (δϕ)2 D

Despejando

4ϕ” + S 2 sin2 α12 C − hS 2 sin2 α12 E − (δϕ)2 D

S cos α12 = (6.36)
B

Figura 25:

En el triángulo esférico anterior P1 tiene coordenadas (ϕ01 , λ01 ); P2 tiene coor-

denadas (ϕ02 , λ02 )

De la figura se tiene:
x = S sin α12

57
y = S cos α12
En donde para el primer caso se determina como una condición.

S2
1− (1 − sin2 α12 sec2 ϕ2 ) = 1
6N22

S sin α12 = x = N2 4λ” cos ϕ2

1
S cos α12 = y = (4ϕ + x2 C + BEx2 + (δϕ)2 D)
B
De la ecuación anterior se tiene:

yB − BEx2 = 4ϕ + x2 C + (δϕ)2 D

4ϕ + x2 C + (δϕ)2 D
y= (6.37)
B(1 − Ex2 )
De la figura se tiene:
x
tanα12 =
y
De donde:
x
α12 = arctan
y
Para calcular ”S” se tiene: usando el teorema de pitagoras:

S 2 = x2 + y 2
De donde
p
S= x2 + y
Tambien
x
S=
sin α12
calculo del α21 :

α21 = α21 ± 180 + 4α

De donde
4ϕ
4α = 4λ” sin ϕm sec
2
Hay que resaltar que S debe ser una distancia restringida, S < 30Km

58