UNIVERSIDAD TECNOLOGICA DE CHILE

INSTITUTO PROFESIONAL
CENTRO DE FORMACION TECNICA

APUNTES DE CLASES
GEODESIA I

DOCENTE: CAROLINA GONZALEZ VERA

INGENIERO GEOMENSOR
AREA DE CONSTRUCCION
INACAP LA SERENA

2015
 GEODESIA I

La asignatura de Geodesia I es una asignatura del área de la Geo-ciencia que pretende entregar a los
estudiantes las herramientas y técnicas de la Geodesia Geométrica para la determinación de sistemas
de referencias terrestres y modelos cartográficos, que les permitan desarrollar proyectos cartográficos,
viales fotogramétricos y mineros, desde su especialidad.

UNIDAD I: INTRODUCCIÓN A LA GEODESIA

La geodesia estudia la división y forma de la tierra, para determinar las posiciones y la dinámica
de la tierra y como esta influye en las mediciones que hacemos. La geodesia es el estudio y medición del
campo gravitatorio del campo terrestre. La geodesia es una ciencia y en ella están integradas la física,
matemática y la geometría.

Lo primero de la geodesia es determinar la forma de la tierra, esto con el fin de garantizar una
medición precisa y una determinación precisa de nuestra posición sobre ella.

La determinación de la forma de la tierra tiene como fin dar un marco de referencia más preciso
para establecer la posición de puntos sobre la superficie terrestre con una alta precisión.

Desde el punto de vista del objetivo de estudio, se puede establecer una división de la geodesia en
diferentes especialidades, aunque cualquier trabajo geodésico requiere la intervención de varias de
estas subdivisiones:

Geodesia Geométrica: Determinación de la forma y dimensiones de la Tierra en su aspecto

geométrico, lo cual incluye fundamentalmente la determinación de coordenadas de puntos en su
superficie.

Geodesia Física: Estudio del campo gravitatorio de la Tierra y sus variaciones, mareas (oceánicas y
terrestres) y su relación con el concepto de altitud.

Astronomía Geodésica: Determinación de coordenadas en la superficie terrestre a partir de

mediciones a los astros.

Geodesia Espacial: Determinación de coordenadas a partir de mediciones efectuadas a satélites

artificiales (GNSS, VLBI, SLR, DORIS) y relación con la definición de sistemas de referencia.

2
1.1 CONCEPTOS BÁSICOS

 Geoide: Superficie Gravitacional equipotencial localizada aproximadamente al nivel medio del

mar (N.M.M), la cual es perpendicular en todo punto a la dirección de la gravedad.

 Elipsoide: Superficie matemática obtenida al ser girar un elipse alrededor del eje polar de la
tierra. Las dimensiones de la elipse se seleccionan al obtener un buen ajuste del elipsoide con el
geoide.

Parámetros elipsoidales

a: semieje mayor b: semieje menor

f: achatamiento e: excentricidad e´: 2da excentricidad

𝑏 √𝑎2 −𝑏2 √𝑎2 −𝑏2
𝑓 =1−𝑎 𝑒= 𝑒 = √2𝑓 − 𝑓 2 𝑒´ =
𝑎 𝑏

 Altura Elipsoidal (h): Altura medida desde una superficie matemática de referencia adoptada
(elipsoide) hasta un punto cualquiera en la superficie topografía.

3
 Altura Ortometrica (H): Altura medida desde el geoide hasta un punto cualquiera en la superficie
topográfica.

 Ondulación Geoidal (N): Diferencia entre la altura elipsoidal y la altura ortometrica.

h= H + N N= h – H

Tipos de Alturas

 Deflexión de la Vertical: Para cualquier punto sobre el terreno, es el ángulo medido entre la
vertical (dirección de la gravedad) y la normal del elipsoide.

4
1.2 SISTEMA DE REFERENCIA

Es una definición conceptual de teorías, hipótesis y constantes que permiten situar una tripleta de
ejes considerados en el espacio definiendo su origen y orientación.

 Sistema de Referencia Terrestre Geocéntrica Convencional : Se utilizan ejes de coordenadas fijas

sobre la tierra, es decir, que giran con ella, de manera que las coordenadas de un punto en
principio sean siempre las mismas.

Origen: Se ubican en el centro de la masa terrestre incluyendo la atmosfera y los océanos.

Eje Z: Coincide con el eje de rotación de la tierra.

Plano Meridiano: Pasa por cualquier punto de la tierra y contiene el eje de rotación.

Plano Ecuador: Es perpendicular al eje de rotación y pasa por el geocentro.

Eje X: Se sitúa en la dirección del plano meridiano que pasa por Greenwich y contenido en el
plano ecuador. El meridiano de Greenwich es considerado internacionalmente como meridiano
de origen o cero desde 1884.

5
 Eje Y: Contenido en el plano ecuador perpendicular al eje x su sentido será tan que los tres ejes
forman una tripleta de dextrógira.

 Sistema de Coordenadas: Es la parametrización de las coordenadas de los puntos que forman el

marco de referencia (materialización de un sistema de referencia), existen infinitos sistemas de
coordenadas por ejemplo: Coordenadas cartesianas, Geográficas, esféricas, astronómicas, etc.

 Coordenadas Cartesianas y Esféricas: Cualquier punto en la superficie terrestre presentan

coordenadas (X, Y, Z) en la tripleta de ejes coordenadas siendo estos puntos los que constituyen
el marco de referencia.

6
En cualquier caso es habitual el uso de coordenadas esféricas para la parametrización de los puntos
sobre la tierra (r, ) donde r Es la distancia radial al geocentro,  Es la latitud geocéntrica y es el ángulo
entre r y el plano ecuador, finalmente Es la longitud geocéntrica y es el ángulo entre el plano
meridiano y el punto de cálculo medido en el plano ecuador.

La relación entre las coordenadas esféricas y cartesianas es un producto de matrices

correspondientes a la parametrización regular de la esfera.

𝑥 𝑐𝑜𝑠𝜌 𝑐𝑜𝑠λ
(𝑦) = 𝑟 (𝑠𝑒𝑛λ cosρ)
𝑧 𝑠𝑒𝑛𝜌

Para la correcta definición de la distancia radial r, se debe introducir una superficie de referencia de
aproximación a la forma real de la tierra.

 Sistema de Coordenadas Geográficas: La primera aproximación la tierra es una esfera

homogénea de radio R.

7
El eje de rotación corta la superficie terrestre en dos puntos, polo norte (PN) y polo sur (PS). El polo
norte es aquel desde que si se observa la tierra hacia su interior, está rota en el sentido contrario a la
agujas del reloj.

La vertical de cualquier punto sobre la superficie terrestre pasa por el centro de la tierra. La posición
de un punto P de la superficie de la tierra se determina por dos coordenadas, Latitud Geográfica y
Longitud Geográfica (φ, λ)

Latitud Geográfica (𝛗): Angulo que porfa la vertical del punto con el P con el plano del ecuador. La
latitud varía de 0° a 90° en el hemisferio norte y de 0° a 90° en el hemisferio sur.

Longitud Geográfica (𝛌): Angulo que forma el meridiano del unto con el meridiano de origen. La
longitud varía de 0° a 180° en el hemisferio oriental y de 0° a 180° en el hemisferio occidental.

8
 Sistema de coordenadas Geodésicas y Geocéntricas: La segunda aproximación de la tierra es un
elipsoide en revolución definido por su semieje mayor (a), su semieje menor (b) y su
achatamiento (f).
El centro del elipsoide coincide con el centro del sistema de referencia, es decir, con el geocentro
y el semieje menor se hace coincidir con el eje de rotación quedando así el sistema de
coordenadas geodésicas.
La vertical geodésica es un punto de la superficie del elipsoide coincide con la dirección del vector
normal del elipsoide en dicho punto y por tanto no pasa por el centro del elipsoide. Así las
coordenadas geodésicas son:

Latitud Geodésicas: Angulo que forma la vertical geodésica del punto con el plano ecuador.

Longitud Geodésica: Angulo formado por el meridiano geodésico del punto de cálculo y el
meridiano geodésico de origen.

Introducimos en este momentos las coordenadas geocéntricas ya que en el elipsoide no

coinciden con las geodésicas, en este caso la longitud geocéntrica será igual a la geodésica, pero
la latitud geocéntrica (B), será el ángulo entre la línea que une el centro del elipsoide y un punto
sobre el elipsoide y plano ecuador geodésico.

9
1.3 Datum – Sistema De Referencia Geodésico

Se define Datum como el conjunto de parámetros que define la posición de un elipsoide con respecto a
la tierra. Para determinarlo nos hace falta conocer la geometría del elipsoide y su posición con respecto
al geocentro ΔX, ΔY, ΔZ (coordenadas de centro del elipsoide respecto al geocentro), su orientación R1,
R2, R3 (orientación de los ejes del elipsoide) y la escala K.

Si de ΔX=ΔY=Δ

Z=0, el Datum geodésico se llama global o absoluto.

Los sistemas de referencia locales corresponden a una estación astronómica como por ejemplo SAP 56
y SAD 69. En el caso de los sistemas geocéntricos y el punto Datum, es el centro de masa de la tierra
como por ejemplo WGS 84.

En nuestro país el IGM (Instituto Geográfico Militar) se ha basado en dos Datum sudamericanos, el
primer Datum provisorio sudamericano año 1956 (PSAD 56) referido al elipsoide internación HAYFORD,
con punto fundamental en la Canoa Venezuela. El segundo Datum sudamericano de año 1969 con punto
fundamental en Chua, Brasil (SAD69).

El año 2003 el IGM oficializa un nuevo Datum denominado sistema de referencia geocéntrica para las
Américas (SIRGAS) utilizando el elipsoide WGS 84.

1.3 Calculo de Coordenadas

 Calculo de Coordenadas Cilíndricas

10
 𝑥
𝑠𝑒𝑛𝜆 = = 𝑦 = 𝑠𝑒𝑛𝜆
𝑟
𝑥 = 𝑟𝑐𝑜𝑠𝜆

𝑧=𝑧

 Coordenadas Cartográficas v/s Coordenadas Esféricas

𝑥 = 𝜌𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜆

𝑦 = 𝜌𝑠𝑒𝑛𝜑𝑠𝑒𝑛2

𝑧 = 𝜌𝑠𝑒𝑛𝜑

 Calculo para un elipsoide en revolución

Se parametriza la superficie del elipsoide de revolución con las coordenadas latitud y longitud (φ, λ)

11
 a= semieje mayor

b= semieje menor

f= achatamiento

𝑒 2 = Primera excentricidad

𝑒 ´2 = segunda excentricidad

𝑎−𝑏 𝑎2 −𝑏 2 𝑎2 −𝑏2
𝑓= 𝑒2 = 𝑒 ´2 =
𝑎 𝑎2 𝑏2

 Calculo de Coordenadas Cartesianas a Geodésicas

Relación de coordenadas cartesianas tridimensionales de un punto del elipsoide en función a

coordenadas latitud y longitud geodésicas. Dado el elipsoide de meridiano en el plano X, Z, se considera
la circunferencia inscrita de radio b y radio a. Se utilizan como paso intermedio una latitud reducida (µ),
sea:

X= a*cos(µ) Z= b*sen(µ)
1
cos(φ) 𝑠𝑒𝑛φ(1−𝑒 2 ) ⁄2
Cos(µ) = 1 Sen(µ)=
(1−𝑒 2 𝑠𝑒𝑛2 φ) ⁄2 √1−𝑒 2 𝑠𝑒𝑛2 φ

12
Al reemplazar los valores de cos(µ) y sen(µ) en las ecuaciones X, Z, obtenemos lo siguiente:
1
a cos(φ) 𝑏 𝑠𝑒𝑛φ(1−𝑒 2 ) ⁄2
X= 1 Z=
(1−𝑒 2 𝑠𝑒𝑛2 φ) ⁄2 √1−𝑒 2 𝑠𝑒𝑛2 φ

Ejercicio 1.1 Realice las siguientes demostración

cos(φ)
Cos(µ) = 1
(1−𝑒 2 𝑠𝑒𝑛2 φ) ⁄2

1
𝑠𝑒𝑛φ(1−𝑒 2 ) ⁄2
Sen(µ)=
√1−𝑒 2 𝑠𝑒𝑛2 φ

Luego es directo llegar a la parametrización del elipsoide de revolución introduciendo las coordenadas
geodésicas denominadas:
𝑎
𝑣= 1⁄
(1 − 𝑒 2 𝑠𝑒𝑛2 φ) 2

Gran normal o normal principal, es el radio de curvatura perpendicular al meridiano.

1⁄ 𝑎
(1 − 𝑒 2 𝑠𝑒𝑛2 φ) 2 =
𝑣

1−(𝑎2 −𝑏2 ) 𝑎2 −(𝑎2 −𝑏2 ) 𝑏2 𝑏

√1 − 𝑒 2 = √ = √ = √ =
𝑎2 𝑎2 𝑎2 𝑎

1 ⁄1
𝑏∗𝑠𝑒𝑛φ(1−𝑒 2 ) ⁄2 𝑏∗𝑠𝑒𝑛φ(1−𝑒 2 ) 2 ∗𝑣 𝑏𝑠𝑒𝑛φ 𝑏𝑣 𝑏2
Z= = Z= = ∗ = ∗ 𝑣𝑠𝑒𝑛φ
√1−𝑒 2 𝑠𝑒𝑛2 φ 𝑎 𝑎 𝑎 𝑎2

Por lo tanto Z= 𝑣(1 − 𝑒 2)sen φ

X=a cos µ * cos λ

Y=a cos µ * cos λ

13
Z= b sen µ

 Coordenadas Cartográficas v/s Coordenadas Geodésicas sin considerar la altura elipsoidal

X= v cos φ * cos λ

Y= v cos φ * sen λ

Z= 𝑣(1 − 𝑒 2 )𝑠𝑒𝑛 φ

 Coordenadas Cartografías v/s Coordenadas Geodésicas considerando altura elipsoidal (h)

X= (v+h) cos φ * cos λ

Y= (v+h) cos φ * sen λ

Z= (𝑣(1 − 𝑒 2 ) + ℎ)𝑠𝑒𝑛 φ

 Coordenadas Geodésicas a Coordenadas Cartesianas

(x,y,z) (φ, λ,h)

P=√𝑥 2 + 𝑦 2
𝑧∗𝑎
𝜃=𝑎𝑟𝑐𝑡𝑔 (𝑃∗𝑏)

𝑧+𝑒 ´2 𝑏𝑠𝑒𝑛3 𝜃
ϕ=𝑎𝑟𝑐𝑡𝑔 𝑝−𝑒 2 𝑎𝑐𝑜𝑠3 𝜃

𝑦
λ = arctg
𝑥
𝑝
ℎ= −𝑣
𝑐𝑜𝑠𝜑

14
Ejercicio 3.I

De las siguientes coordenadas geodésicas a partir del Datum WGS 84 calcule las coordenadas
cartesianas.

𝜑 = 36° 42´ 7,0337"

λ = −2° 50´ 26,7739"

h= 67,593

Semieje mayor= 6378137

Semieje menor= 6356752,3

15
UNIDAD II: RADIO DE CURVATURA Y REDUCCIÓN DE LAS OBSERVACIONES

2.1 Radio de Curvatura Principal

Todos los planos que pasan por la normal a un punto de la elipsoide cortan este en varias secciones
normales, que en su mayoría son elipses, de las infinitas secciones normales que pueden trazarse de
punto M del elipsoide, hay que señalar dos secciones principales, la sección meridiana y la sección
normal.

1° sección meridiana: Corresponde a un mínimo radio de curvatura y se denomina con la letra 𝜌 y su

expresión es:

𝑎(1 − 𝑒 2 )
𝜌= 3⁄
(1 − 𝑒 2 𝑠𝑒𝑛2 𝜑) 2

2° sección normal: Corresponde al máximo radio de curvatura y está comprendida por el segmento KM
de la normal al elipsoide se denomina gran normal o vertical primaria, su expresión es:

𝑎
𝑣= 1⁄
(1 − 𝑒 2 𝑠𝑒𝑛2 𝜑) 2

Al observar las expresiones de los radios de curvatura principalmente podemos identificar que la única
variable independiente que se utiliza es la latitud (𝜑) es decir, ambos radios permanecen constantes a
lo largo del mismo paralelo, siendo la curvatura en un punto M del elipsoide variable en la sección normal
que se considera, por lo que tomaremos como radio de curvatura del elipsoide en un punto considerado
como la media aritmética de los infinitos radios de curvatura correspondientes a las infinitas secciones
normales que pasan por el punto M. Con lo anterior podemos decir que el radio medio será igual a la
media geométrica de los radios de curvatura principal, se expresión queda dada por:

16
 𝑅𝑚 = √𝜌 ∗ 𝑣

 Los radios de curvatura principales iguales en los polos y en el ecuador valen:

𝑏2
V 0= a 𝜌 = 𝑎2

Ejercicio 2.1

Determine cada punto de los siguientes radios de curvatura.

a) Radio de curvatura de la elipse meridiana

b) Radio de curvatura del primer vertical o gran normal
c) Radio de curvatura medio

Punto 𝜑 λ h

FARO 43° 19´ 20,9406” -2° 0´ 20,9010” 188,351

CARCHE 38° 25´ 44,0590” -1° 09´ 46,0052” 1343,765

ELIPSOIDE DE REFERENCIA: HAYFORD

2.2 Valor Lineal De Los Arcos De Paralelo Y De Meridiano

Para determinar el valor lineal de un arco de paralelo, se calcula previamente el radio de este. El radio
del paralelo que pasa por un punto P es OP´ cateto del rectángulo OPA, en el que PQ=V es a gran normal
según se deduce en la siguiente figura:

17
Valor Lineal de un arco de paralelo:

𝑅𝑝 = 𝑉𝑠𝑒𝑛(90° − 𝜑𝑝) = 𝑉𝑐𝑜𝑠𝜑𝑝

Siendo 𝜑𝑝 la latitud del punto P

El valor del arco de paralelo será la diferencia de longitudes de sus extremos P y R, es decir, λP − λR. Si
tenemos en cuenta que la longitud total del paralelo será:

L= 2𝜋* Rp

2𝜋𝑅𝑝(λP − λR)" 𝑅𝑝 (λP − λR)"

𝑃𝑅 = =
360 ∗ 60 ∗ 60 𝑟

λP − λR Es la diferencia de longitudes expresadas en segundos y r es el número de segundos de un

radian (r=206,265).

Por lo tanto, S = longitud de un arco de paralelo es:

S = Rpx Δλ Rp= Vcos𝜑

18
2.3 Longitud De Un Arco De Meridiano

La primera cuadrática fundamental de una superficie del elipsoide de reducción parametrizado por la
latitud y longitud geodésica responde a la expresión:

𝐼(𝑑𝜑, 𝑑λ) = d𝑠 2 = 𝜑 2 𝑑𝜑 2 + 𝑣 2 𝑐𝑜𝑠 2 𝜑𝑑λ2

Por lo tanto, la longitud de una curva cualquiera sobre la superficie del elipsoide de revolución de un
punto 1 a un punto 2 se determina por:

λ2 𝜑2

𝑆 = ∫ ∫ √𝜌2 𝑑𝜑 2 + 𝑣 2 𝑐𝑜𝑠 2 𝜑𝑑λ2

λ1 𝜑1

Si se particulariza en el caso de que la curva será un arco meridiano, es decir, una de las curvas
paramétricas, la expresión se simplifica a:

𝜑2

𝑆 = ∫ 𝜌𝑑𝜑
𝜑1

Por lo tanto, la expresión permite obtener la longitud de un arco de meridiano resultante como:

𝜑2
𝜑2 𝑑𝜑
𝐿. 𝑎. 𝑚𝜑1 = 𝑎(1 − 𝑒 2 ) ∫ 3⁄
𝜑1
(1 − 𝑒 2 𝑠𝑒𝑛2 𝜑) 2

Es evidente que la longitud de un arco de meridiano desde un paralelo 𝜑, a otro 𝑒2 , se puede obtener a
partir de las longitudes de los arcos de meridiano que van del ecuador a los paralelos anteriores.

Se define la longitud del arco de meridiano de un determinado paralelo de latitud 𝜑 como la solución
integral definida por:
𝜑
𝑑𝜑
𝐿. 𝑎. 𝑚(𝜑) = 𝑎(1 − 𝑒 2 ) ∫ 3⁄
0
(1 − 𝑒 2 𝑠𝑒𝑛2 𝜑) 2

1
Sea 𝑓(𝜑) = 3 𝑥 = 𝑠𝑒𝑛𝜑
(1−𝑒 2 𝑠𝑒𝑛2 𝜑) ⁄2

19
 1
𝑓(𝑥) = 3⁄
(1 − 𝑒 2 𝑠𝑒𝑛2 𝜑) 2

Desarrollando en serie de Taylor nos queda:

3 15 35 6 6
𝑓(𝑥) = 1 + 𝑒 2 𝑥 2 + 𝑒 4 𝑥 4 + 𝑒 𝑥 … … . 𝑅(𝑥12 )
2 8 16

Si se realiza el cambio de variable resulta:

𝜑

𝐿. 𝑎. 𝑚 = ∫ 𝜌𝑑𝜑 = 𝑎(1 − 𝑒 2 )(𝑔1 + 𝑔2 + 𝑔3 + 𝑔4)

0

Siendo: 𝑔1 = 𝜑

3 2 −1 1
𝑔2 = 𝑒 ( 2 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜑 + 2 𝜑)
2

15 4 −1 3 3
𝑔3 = 𝑒 ( 4 𝑠𝑒𝑛3 𝜑𝑐𝑜𝑠𝜑 − 8 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜑 + 8 𝜑)
8

35 6 −1 5 5 5
𝑔4 = 𝑒 ( 6 𝑠𝑒𝑛5 𝜑𝑐𝑜𝑠𝜑 − 24 𝑠𝑒𝑛3 𝜑𝑐𝑜𝑠𝜑 − 16 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜑 + 16 𝜑
16

Por lo tanto, esta fórmula nos da la longitud de un arco meridiano desde cada punto al ecuador en el
meridiano sobre el que tratamos.

20
Ejercicio 2.2

Determine la distancia correspondiente al arco de paralelo del punto 1 al resto, en el paralelo del punto
1 (WGS 84)

Punto 𝜑 λ A.H. Elipsoidal

1 36,42070337 -2.50267739 67,593

2 40,49007808 -5,36526511 1066,000

3 43,33158796 -8,14256201 364,334

4 41,03343188 1,10121658 129,358

5 43,19170007 -2,00344814 255,681

6 38,25396946 -1,09503671 1424,348

a= 6378137,0

b= 6356752,3142

R1= Vcos𝜑

S1=∆ λ ∗ Rp

Ejercicio 2.3

Determina la distancia correspondiente al arco de meridiano del punto 1 al resto, en el meridiano del
punto 1 (WGS 84).

21
2.4 Sección Normal Y Línea Geodésica Sobre El Elipsoide

Sea un punto sobre la superficie del elipsoide y sea la dirección de la normal el elipsoide en P. Para
cualquier otro punto del elipsoide existe una sección normal que la contiene.

Para un punto Q, se define el plano normal de P a Q como el que contiene a la normal en P y al punto Q
sobre el elipsoide. Se ha definido un plano mediante una recta y un punto. Se define la sección normal
de P a Q como la curva de intersección del plano normal P a Q y la superficie del elipsoide en revolución.

Es evidente que la normal Q a P no coincidirán con el P a Q debido a que el segundo se define por la
normal de Q a P no coincide con la de P a Q.

2.5 Línea Geodésica

El problema que se presenta en el apartado anterior sobre la duplicidad sobre la línea que tiene dos
puntos sobre el elipsoide, sección normal directa y reciproca se resuelve en la línea geodésica.

La línea geodésica entre dos puntos se define como la curva de menos longitud que los une. En un punto
sería una recta, en una superficie esférica sería un círculo máximo y en un elipsoide sería una curva
definida por tres ecuaciones diferentes.

𝑑𝑦 𝑐𝑜𝑠𝑎 𝑑λ 𝑠𝑒𝑛 𝑎 𝑑𝑎 𝑠𝑒𝑛 𝑎

> > ∗ 𝑠𝑒𝑐𝜑 > ∗ 𝑡𝑔𝜑
𝑑𝑠 𝑝 𝑑𝑠 𝑣 𝑑𝑠 𝑣

22
Este sistema de ecuaciones tiene un significado muy importante dado que es el sistema de partida para
resolver los problemas geodésicas directo e inverso que se estudiaran en la unidad III de este curso.

2.6 Medida De Distancia Sobre Una Superficie

En la geometría del plano euclideo la distancia entre dos puntos de coordenadas (x,y) y (x+dx, y+dx)
viene dado por el vector que los une √𝑑𝑥 2 + 𝑑𝑦 2.

Se plantea determinar la distancia entre dos puntos de una superficie S. se considera que un entorno de
un punto P de S suficientemente pequeño, la superficie coincidirá con el p y tangente a ella en P, y por
tanto se puede trasladar las propiedades de la superficie a las del plano.

Sean los puntos X y X+DX, la distancia euclideo viene dada por:

𝑑𝑠 2 = 𝑑𝑥 ∗ 𝑑𝑥 = (𝑥𝜇 ∗ 𝑑𝜇 + 𝑥𝑣 ∗ 𝑑𝑣)(𝑥𝜇 ∗ 𝑑𝜇 + 𝑥𝜇 ∗ 𝑑𝑣)

𝑑𝑠 2 = 𝐸𝑑𝜇 2 + 2𝐹𝑑𝜇𝑑𝑣 + 𝐺𝑑𝑣 2

Siendo: E= x𝜇 * x𝜇

F= x𝜇* xv Primero formula cuadrática fundamental

G= xv* xv

23
Consideramos 𝜇 = 𝜇(𝑡) 𝑦 𝑣 = (𝑡), una curva sobre S y X (t+dt) dos puntos infinitamente próximos a
dicha curva. Entonces la distancia entre estos dos puntos medidos a lo largo de esta curva viene dado
por:

𝑑𝑦 2 𝑑𝑦 𝑑𝑣 2
𝑑𝑠 = √𝐸 ( ) + 2𝐹 ( ∗ ) ∗ 𝑑𝑡
𝑑𝑡 𝑑𝑡 𝑑𝑡

A continuación se determina la primera forma cuadrática fundamental para un elipsoide en revolución

retornando la expresión de los vectores tangentes a las curvas paramétricas.

𝐸 = 𝑥𝜑 ∗ 𝑥𝜑 = 𝜌2 𝑠𝑒𝑛2 𝜑𝑐𝑜𝑠 2 λ + 𝑝2 𝑠𝑒𝑛2 𝜑𝑠𝑒𝑛2 λ + 𝑝2 𝑐𝑜𝑠 2 𝜑 = 𝜌2

𝐹 = 𝑥𝜑 ∗ 𝑥λ = ρsenφcosλ ∗ vcosφsenλ − ρsenφsenλ ∗ vcosφcosλ = 0

𝐺 = 𝑥λ ∗ xλ = 𝑣 2 𝑐𝑜𝑠 2 𝜑𝑠𝑒𝑛2 λ + 𝑣 2 𝑐𝑜𝑠 2 𝜑𝑐𝑜𝑠 2 λ = 𝑣 2 𝑐𝑜𝑠 2 𝜑

Por lo tanto, 𝐼(𝑑𝜑, 𝑑𝑥) = 𝜌2 𝑑𝜑 2 + 𝑣 2 𝑐𝑜𝑠 2 𝜑𝑑λ2

Ahora se evalúa esta expresión sobre las líneas productivas.

𝑑λ = λ2 − λ1

𝑑𝜑 = 𝜑2 − 𝜑1

𝜑 𝜑1−𝜑2
𝑚=
2

𝑎(1 − 𝑒 2 )
𝜌𝑚 = 3⁄
(1 − 𝑒 2 𝑠𝑒𝑛2𝜑𝑚 ) 2

𝑎
𝑣𝑚 = 1⁄
(1 − 𝑒 2 𝑠𝑒𝑛2 𝜑𝑚 ) 2

Por lo tanto la distancia sobre un elipsoide entre dos untos queda definido por:

𝑑λ
𝑑𝑠 = √𝑣𝑚2 𝑐𝑜𝑠 2 𝜑𝑚 𝑑λ2 + 𝜌𝑚2 𝑑𝜑 2 𝑐𝑜𝑠 2 ( )
2

24
Ejercicio 2.4

En base a los datos y resultados obtenidos en el ejercicio anterior, calcule la distancia (ds):

Dato: Punto 1 fijo

Ejercicio 2.5

Con respecto al Datum WGS 84:

a) Determinar el radio de curvatura de la elipsoide de meridiano 𝜌

b) Determinar el radio de curvatura del primer vertical v
c) Determinar el radio medio de curvatura

Punto 𝜑 𝑣 𝑅𝑚

1 6358230,968 6385776,262 6371988,730

2 6362717,517 6387277,904 6374985,883

3 6365764,237 6388297,235 6377020,784

4 6362986,162 6387367,796 6375165,323

5 6365503,837 6388210,126 6376846,875

6 6360094,240 6386399,982 6373233,539

25
UNIDAD III: DETERMINACIÓN DE COORDENADAS GEODESICAS

Los problemas geodésicos principales se refieren en concreto al tema del transporte de coordenadas de
los puntos sobre el elipsoide, es decir, al cálculo de las coordenadas de un punto a partir de las
coordenadas conocidas de otro vértice, o bien al cálculo de la distancia y de los azimut correspondientes
entre dos puntos de coordenadas conocidas. En el primer caso el problema geodésico se conoce como
problema directo y se enuncia como sigue:

Dadas las coordenadas geodésicas ϕ y λ de un punto sobre el elipsoide, el azimut geodésico  de la
dirección a un segundo punto y la distancia geodésica S entre 1 y 2, se pretenden obtener la longitud y
latitud geodésicas del punto 2 así como el azimut inverso.

Por su parte, el segundo caso constituye el problema inverso cuyo enunciado resulta ser:

Conocidas las coordenadas geodésicas de dos puntos 1 y 2 sobre la superficie del elipsoide, se han de
calcular la distancia geodésica s entre ambos y los azimut de la dirección en los extremos de la línea
geodésica que los une.

26
Estos dos problemas, que se complementan con el cálculo de la convergencia de meridianos entre los
puntos 1 y 2, son la base para el establecimiento posterior de la red geodésica.

Ambos se pueden resolver por distintos procedimientos desde el punto de vista de la Geodesia clásica,
si bien cada uno de los métodos se caracteriza por el margen de distancias en que puede utilizarse al
objeto de obtener resultados aceptables en precisión. Atendiendo a la magnitud de la distancia de
separación entre los puntos 1 y 2, clasificaremos los procedimientos existentes según la siguiente
subdivisión:

– Métodos para distancias cortas. Empleo hasta 30 ó 40 kilómetros como máximo.

– Métodos para medias distancias. Utilizables hasta los 250 ó 300 kilómetros.
– Métodos para distancias grandes. Pueden emplearse hasta los 5000 kilómetros.
– Métodos para distancias muy grandes. Útiles hasta 20000 kilómetros y más.

La clasificación expuesta es típica de la Geodesia Geométrica. Sin embargo, con la aparición de las
tecnologías que dieron lugar a la Geodesia espacial, el estudio de los métodos clásicos diseñados para
transporte de coordenadas sobre el elipsoide tiene más interés teórico que práctico. Por tal razón
prescindiremos de tratar aquí los procedimientos elaborados para distancias muy grandes, limitando
nuestro estudio al de aquellos métodos aplicables hasta distancias medias o grandes que aún pudieran
ser de utilidad bajo diversas circunstancias.

Para el desarrollo de estos problemas clásicos en Geodesia se estudiara el método de distancias cortas
mediante la fórmula de las ecuaciones de Puissant.

27
3.1 FÓRMULAS DE PUISSANT.

Las ecuaciones de Puissant, comúnmente usadas, son válidas para bases cuya longitud no sea superior
a 80 km, debido al hecho de que estas fueron deducidas con simplificaciones para un modelo esférico,
se desarrollaron en series matemáticas y adoptan un radio medio terrestre para la región levantada. La
deducción de las ecuaciones de Puissant puede ser encontrada en HOSMER (1946) y en GEMAEL (1959).

3.1.1. Problema Directo.

Factores constantes utilizados para coordenadas del punto 2:

1 3
(1−𝑒 2 𝑠𝑒𝑛2 𝜑2 ) ⁄2 ⁄
(1−𝑒 2 𝑠𝑒𝑛2 𝜑1 ) 2
2
(1−𝑒 2 𝑠𝑒𝑛2 𝜑1 ) 𝑡𝑔𝜑1
𝑇= 𝐵= 𝐶=
𝑎𝑠𝑒𝑛1" 𝑎(1−𝑒 2 )𝑠𝑒𝑛1" 2𝑎2 (1−𝑒 2 )𝑠𝑒𝑛1"

3⁄ 𝑒 2 𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑 𝑠𝑒𝑛1"

2 1 1
𝐷= 2 2
(1−𝑒 𝑠𝑒𝑛 𝜑1 )

∆𝜑" = −(𝐵 ∗ 𝑆 ∗ 𝑐𝑜𝑠𝐴12 + 𝐶 ∗ 𝑆 2 ∗ 𝑠𝑒𝑛2 𝐴12 + 𝐷(𝐵 ∗ 𝑆 ∗ 𝑐𝑜𝑠𝐴12 )2 )

𝜑2 = 𝜑1 + ∆𝜑
𝑇 ∗ 𝑆 ∗ 𝑠𝑒𝑛𝐴12
∆λ" =
𝑐𝑜𝑠𝜑2

λ2 = λ1 + ∆λ
Azimut 21

𝜑1 + 𝜑2
𝜑𝑚 =
2
1
𝐹= ∗ 𝑠𝑒𝑛 𝜑𝑚 𝑐𝑜𝑠 2 𝜑𝑚 𝑠𝑒𝑛2 1"
12
Δ𝜑
γ" = Δλ" ∗ sen𝜑𝑚 ∗ 𝑠𝑒𝑐 ( ) + (Δλ")3 ∗ 𝐹 Convergencia meridiana
2

𝐴21 = 𝐴12 + 𝛾 + 180°

28
3.1.2. Problema Inverso.

𝑁𝑚 ∗ Δλ ∗ cosφm 𝛾"
𝐴1−2 = 𝑎𝑟𝑐𝑡𝑔 [ ]−
Δλ 2
𝜑𝑚 ∗ (−Δ𝜑)cos ( )
2

Δ𝜑
𝛾" = Δλ" ∗ sen𝜑𝑚 ∗ 𝑠𝑒𝑐 ( ) + (Δλ")3 ∗ 𝐹
2

1
𝐹= ∗ 𝑠𝑒𝑛 𝜑1 ∗ 𝑐𝑜𝑠 2 𝜑1 𝑠𝑒𝑛2 1"
12

29
3.2 GRAFICO

𝜃1 𝑦 𝜃2 = Azimut Geodésico (A1 y A2 ecuación Puissant) Se mide desde el Sur de la proyección de NV

(Norte verdadero o astronómico) en sentido de los punteros del reloj.

𝑇1 𝑦 𝑇2 = Azimut geodésico proyectado. Es una magnitud angular real que se diferencia del azimut
geodésico verdadero 𝜃.

𝐶1 𝑦 𝐶2 = convergencia.

𝑡1 𝑦 𝑡2 = Azimut plano. Si sobre el plano de la proyección se unen los puntos 𝑃1 𝑦 𝑃2 con una recta. Se
determina la línea de cuadricula para obtener el azimut plano t, el cual parte en el Nc (norte cuadricula)
que une los puntos 𝑃1 𝑦 𝑃2 .

30
UNIDAD IV: PROYECCIÓN TRANSVERSAL DE MERCATOR

COORDENADAS UTM

4.1 Sistema de Coordenadas UTM

El sistema de coordenadas universal transversal de Mercator, es un sistema de coordenadas basado en

la proyección cartográfica transversa de Mercator, que se construye como la proyección de Mercator
normal, pero en vez de hacerla tangente al Ecuador, se la hace tangente a un meridiano.

A diferencia del sistema de coordenadas geográficas, expresadas en longitud y latitud, las magnitudes
en el sistema UTM se expresan en metros únicamente al nivel del mar, que es la base de la proyección
del elipsoide de referencia.

El sistema de coordenadas UTM fue desarrollado en la década de 1940. El sistema se basó en un modelo
elipsoidal de la Tierra. En su inicio se usó el elipsoide de Clarke de 1866 para el caso de Estados Unidos.
Para el resto del mundo –incluidos Alaska y Hawái– se usó el Elipsoide Internacional. Actualmente se usa
el elipsoide WGS84 como modelo de base para el sistema de coordenadas UTM.

El sistema UTM implica el uso de escalas no lineales para las coordenadas X e Y (longitud y latitud
cartográficas) para asegurar que el mapa proyectado resulte conforme. las coordenadas utm, son
también reconocidas como coordenadas planas.

31
4.2 Proyección Transversal de Mercator

La UTM es una proyección cilíndrica conforme. El factor de escala en la dirección del paralelo y en la
dirección del meridiano son iguales (h = k). Las líneas loxodrómicas (línea que une dos puntos cualesquiera
de la superficie terrestre cortando a todos los meridianos con el mismo ángulo) se representan como líneas
rectas sobre el mapa. Los meridianos se proyectan sobre el plano con una separación proporcional a la
del modelo, así hay equidistancia entre ellos. Sin embargo los paralelos se van separando a medida que
nos alejamos del Ecuador, por lo que al llegar al polo las deformaciones serán infinitas. Por eso sólo se
representa la región entre los paralelos 84ºN y 80ºS. Además es una proyección compuesta; la esfera se
representa en trozos, no entera. Para ello se divide la Tierra en husos de 6º de longitud cada uno.

La proyección UTM tiene la ventaja de que ningún punto está demasiado alejado del meridiano central
de su zona, por lo que las distorsiones son pequeñas. Pero esto se consigue al coste de la discontinuidad:
un punto en el límite de la zona se proyecta en coordenadas distintas propias de cada Huso.

Para evitar estas discontinuidades, a veces se extienden las zonas, para que el meridiano tangente sea
el mismo. Esto permite mapas continuos casi compatibles con los estándar. Sin embargo, en los límites
de esas zonas, las distorsiones son mayores que en las zonas estándar.

4.3 Husos

Se divide la Tierra en 60 husos de 6º de longitud, la zona de proyección de la UTM se define entre los
paralelos 80º S y 84º N. Cada huso se numera con un número entre el 1 y el 60, estando el primer huso
limitado entre las longitudes 180° y 174° W y centrado en el meridiano 177º W. Cada huso tiene asignado
un meridiano central, que es donde se sitúa el origen de coordenadas, junto con el ecuador. Los husos
se numeran en orden ascendente hacia el este. En el sistema de coordenadas geográfico las longitudes
se representan tradicionalmente con valores que van desde los -180º hasta casi 180º (intervalo -180º →
0º → 180º); el valor de longitud 180º se corresponde con el valor -180º, pues ambos son el mismo.

4.3.1 Husos en Chile

- Chile continental: huso 18 y 19

- Chile insular: huso 12, 13, 17
- Territorio antártico: 16, 17, 18, 19, 20, 21 y 22.

32
4.4 Bandas

Se divide la Tierra en 20 bandas de 8º Grados de Latitud, que se denominan con letras desde la C hasta
la X excluyendo las letras "I" y "O", por su parecido con los números uno (1) y cero (0), respectivamente.
Puesto que es un sistema norteamericano (estadounidense), tampoco se utiliza la letra "Ñ". La zona C
coincide con el intervalo de latitudes que va desde 80º Sur (o -80º latitud) hasta 72º S (o -72º latitud).
Las bandas polares no están consideradas en este sistema de referencia. Para definir un punto en
cualquiera de los polos, se usa el sistema de coordenadas UPS (proyección estereografica polar). Si una
banda tiene una letra igual o mayor que la N, la banda está en el hemisferio norte, mientras que está en
el sur si su letra es menor que la "N".

4.5 Geometría del Huso

En la figura se ilustra el resultado de proyectar el huso 31 según su meridiano central (3º E). Como se ve,
éste lo divide en dos partes iguales. Esto permite establecer dos ejes cartesianos X e Y sobre el huso, de
tal manera que el eje X es el ecuador y el eje Y el meridiano central. Estos ejes cartesianos permiten,

33
pues, determinar puntos sobre el huso haciendo uso de dos coordenadas rectangulares X e Y, que se
denominan coordenadas UTM.

El origen del sistema de coordenadas UTM se encuentra en la intersección del Ecuador con el meridiano
central del huso. Cada huso, pues, posee su propio origen de coordenadas.

La idea de las coordenadas UTM es que sus dos valores X e Y siempre sean positivos. Por ello no se han
elegido las coordenadas X=0 e Y=0 para el origen.

El meridiano central tiene un valor de cuadricula de 500.000 mts. creciente hacia el este y decreciente
al oeste y el ecuador con un valor de cuadricula de 0 mts. Para el hemisferio norte y en aumento hacia
el polo norte y de 10.000.000 mts para el caso del hemisferio sur disminuyendo hacia el polo sur.

4.6 Transformación De Coordenadas Geodésicas A Coordenadas UTM -Método Coticchia Surace

Coordenadas UTM 𝑥: 𝜉 ∗ 𝑉 (1 − 𝜍) + 500.000

𝑦: 𝜂 ∗ 𝑉(1 − 𝜍) + 𝑏0

𝑎2 √𝑎2 −𝑏2 1 1+𝐴

𝑐= 𝑒´ = 𝐴 = 𝑐𝑜𝑠𝜑𝑠𝑒𝑛Δ λ 𝜉= 𝑙𝑛 [1−𝐴]
𝑏 𝑏 2

𝑡𝑔𝜑 𝐶 𝑒´2
𝜂 = 𝑎𝑟𝑐𝑡𝑔 [𝑐𝑜𝑠Δλ] 𝑉= 1 ∗ 0,9996 𝜍= ∗ 𝜉 2 ∗ 𝑐𝑜𝑠 2 𝜑
(1−𝑒´2 𝑐𝑜𝑠2 𝜑) ⁄2 2

𝐴
𝐴1 = 𝑠𝑒𝑛(2𝜑) 𝐴2 = 𝐴1 ∗ 𝑐𝑜𝑠 2 𝜑 𝐽2 = 𝜑 + ( 21 )

(3∗𝐽2 )+𝐴2 (5∗𝐽4 )+(𝐴2 ∗𝑐𝑜𝑠2 𝜑) 3 5

𝐽4 = 𝐽6 = 𝛼 = 4 ∗ 𝑒 ´2 𝛽 = 3 ∗ 𝛼2
4 3

35
𝛾 = 27 ∗ 𝛼 3 𝛽0 = 0,9996 ∗ 𝐶(𝜑 − 𝛼𝐽2 + 𝛽𝐽4 − 𝛾𝐽6 )

34
4.7 Transformación De Coordenadas UTM A Geodésicas – Modelo Coticchia Surace

λ = ∆λ + λ0

3
𝜑 = 𝜑´ + [1 + (𝑒 +2 ∗ 𝑐𝑜𝑠 2 𝜑´ − 𝑒 ´2 ∗ 𝑠𝑒𝑛𝜑´𝑐𝑜𝑠𝜑´(𝜏 − 𝜑´))] (𝜏 − 𝜑´)
2

𝑣 𝐶 𝑥
𝜑´ = 6366197,723∗0,9996 𝑣= 1⁄ ∗ 0,9996 𝑎=𝑣
(𝑒 ´2 𝑐𝑜𝑠2 𝜑´) 2

𝐴1 (3∗𝐽2 )+𝐴2
𝐴1 = 𝑠𝑒𝑛(𝜑´) 𝐴2 = 𝐴1 ∗ 𝑐𝑜𝑠 2 𝜑´ 𝐽2 = 𝜑´ + 𝐽4 =
2 4

5∗𝐽4 +𝐴2 𝑐𝑜𝑠2 𝜑´ 3 5 35

𝐽6 = 𝛼 = 4 ∗ 𝑒 ´2 𝛽 = 3 ∗ 𝛼2 𝛾 = 27 ∗ 𝛼 3
3

𝑦−𝐵0
𝛽0 = 0,9996 ∗ 𝐶(𝛾´ − 𝛼𝐽2 + 𝛽𝐽4 − 𝛾𝐽6 𝑏= 𝑣

𝑒 ´2 ∗𝑎2 𝜁
𝜁= ∗ 𝑐𝑜𝑠 2 𝜑´ 𝜉 = 𝑎 ∗ [1 − 3] 𝜂 = 𝑏 ∗ (1 − 𝜁) + 𝜑´
2

𝑠𝑒𝑛 𝜂∗𝜉
𝑠𝑒𝑛(𝜂 ∗ 𝜉) = 𝑒 𝜉 − 𝑒 (−𝜉) Δλ = arctg ( ) 𝜏 = 𝑎𝑟𝑐𝑡𝑔(cos(Δλ) tg η)
cos 𝜂

35
Ejercicio 4.1

Basándose en los datos realice una transformación de coordenadas geodésicas a coordenadas UTM

a= 6378137

b= 6356752,314

𝜑 = 23°33´32,95" − 𝑆

λ = 70°24´10,06" − W

Ejercicio 4.2

Basándose en los datos realice una transformación de coordenadas geodésicas a coordenadas UTM

a= 6378137

b= 6356752,314

𝜑 = 23°34´17,17" − 𝑁

λ = 70°23´33,75" − W

36
BIBLIOGRAFIA

 Manual de Normas Técnicas de Mensura, 2010 - Ministerio de Bienes Nacionales

 Sistema y Marco de Referencia Terrestre. Sistemas de Coordenadas, 2010, Ángel Martín Furones,
Universidad Politécnica de Valencia
 Apuntes de clases Ingeniería Civil en Geomática, Patricio Alcota Aguirre, Universidad de
Antofagasta
 Mena, J.B., 2008. Geodesia Superior. Texto en dos tomos. En proceso de publicación por el
Instituto Geográfico Nacional. Madrid.
 Múgica, F. 1971. El ordenador y el análisis numérico. Boletín del Servicio Geográfico del Ejército,
1971.
 Zakatov, P.S., 1997. Curso de Geodesia Superior. Ed. Mir. Moscú.
 Geodesia y Cartografía Matemática, David Hernández, Universidad Politécnica De Valencia
 www.wikipedia.org

37