10-04-02 [Link]- pag.

1 de 28

GEODESIA GEOMÉTRICA
RENÉ ZEPEDA G. - AGOSTO 2001

ESTOS APUNTES SON PROVISORIOS, ESTÁN SUJETOS A REVISIÓN Y CAMBIOS

Petr Vaniceck: “... la llave del conocimiento reside en el dominio de sus conceptos ...”

INTRODUCCIÓN

La palabra Geodesia, tiene como origen la palabra griega Geodaisia, que significa “división de tierras”
(geo = Tierra y daisia = dividido)

Por lo expuesto en diversas enciclopedias y diccionarios, Geodesia puede ser definida como una
ciencia cuyo objetivo es determinar la forma de la Tierra y calcular sus dimensiones. Posee dos
campos de pesquisa: uno teórico, que examina la configuración de la Tierra en su conjunto,
considerando los factores internos y externos que la determinan; otro más práctico, que partiendo de la
elaboración de los dados teóricos, prepara las soluciones apropiadas para representación cartográfica
de la superficie terrestre.

Camil Gemael (1981), “el objetivo de la Geodesia es la determinación de la forma y de las

dimensiones de la Tierra. Encontrándose dividida en: Geométrica, Física y Celeste. la Geodesia
Geométrica es ejecutada a través de la medición de los ángulos y/o distancias, proporcionando el
cálculo de las coordenadas elipsoidales de un punto P, de la superficie física de la Tierra, sobre el
modelo de referencia. La Física se preocupa con el estudio de la gravedad y sus aplicaciones
geodésicas, mientras que la Celeste permite la determinación de la posición relativa o absoluta de
puntos de la superficie terrestre.”

Torge (1980), “la Geodesia puede ser dividida en global (global geodesy), de levantamientos (geodetic
survey) y de levantamientos planos (plane surveying). La Geodesia Global es responsable por la
determinación de la figura de la Tierra, incluido el campo de gravedad externa. La Geodesia de
Levantamientos es responsable por la definición de las redes nacionales establecidas en los países. El
Levantamiento Plano (levantamiento topográfico, catastral, etc.) es responsable por el detalle de la
 10-04-02 [Link]- pag. 2 de 28

superficie; siendo o plano horizontal, normalmente, superficie de referencia en este caso. La

integración entre ellas, se da de la siguiente forma: La Geodesia Global es responsable por la
definición de los parámetros que determinan la forma de la tierra y, en consecuencia, de las medidas
realizadas o conducidas, en su superficie. El Levantamiento Plano, a su vez, hace uso de los puntos de
control establecidos por el Levantamiento Geodésico, y que son utilizados para el mapeo nacional y en
los catastros estatales.”

Langley (1992), dice que “es la ciencia que determina el tamaño y la forma de la Tierra, incluido su
campo de gravedad, en cuatro dimensiones, espacio-tiempo”. Considerando esta definición, es
necesario por lo tanto, definir el sistema de coordenadas a ser adoptado, describir el campo de
gravedad global y estudiar las variaciones temporales de posiciones, sean ellas por causas naturales o
no. Aún en esta línea, Langley (1992) afirma que “cuando la Geodesia se refiere la superficie de la
Tierra, significa referirse la superficie equipotencial. Aunque existan diversas superficies
equipotenciales, apenas una tiene especial significado, es aquella que más se aproxima de la superficie
de los océanos en reposo, y no a su nivel, cuando eles se prolongan bajo os continentes y se
encuentran libres de los efectos de las mareas, ondas, vientos, corrientes, etc. Esta superficie es
llamada geoide.”

La forma real de la Tierra, según Torge (1991), “es el geoide, definido como la superficie
equipotencial, que en cualquier lugar es perpendicular a la vertical dada por un hilo de plomo y que
coincide con o nivel medio no perturbado de los mares.” Como el geoide es una superficie irregular,
no puede ser matemáticamente definida; es importante no solo para la investigación científica así
también, para diversas actividades cotidianas. A través del mapeo del geoide, se puede verificar la
estructura de la costra terrestre y acompañar la evolución de la tectónica de placas; en actividades
cotidianas, su uso pode ser comprobado la través de las altitudes referidas a la superficie, altitudes
ortométricas, que normalmente son empleadas en el mapeo topográfico.

Como el geoide es de difícil representación, la forma de la Tierra ha sido matemáticamente definida

por un elipsoide de revolución; que es la figura geométrica que más se aproxima a la forma real de la
tierra: achatada en los polos y alargada en el Ecuador. La superficie elipsoidal es conveniente como
referencia y facilita las operaciones matemáticas. En razón de ello, esta es la superficie de referencia
mas ampliamente empleada en levantamientos y mapeos; pues por ser una superficie
matemáticamente desarrollada, es largamente utilizada en proyecciones cartográficas y en el
establecimiento de coordenadas horizontales de las redes geodésicas, permitiendo la ejecución de
cálculos diversos, con una precisión necesaria para la cartografía de grandes áreas.
 10-04-02 [Link]- pag. 3 de 28

El elipsoide (superficie elipsoidal) es menos usado como superficie de referencia para las coordenadas
verticales (altitudes), ya que no refleja una superficie física de nivel, pero sí una superficie geométrica.
“La determinación de altitud, respecto al nivel medio de los mares es la operación normalmente
ejecutada en nivelación. Esta altitud pode ser interpretada como la altura ortométrica, ya que o geoide
es muy próximo al Nivel Medio de los Mares (NMM)”.

Bomford (1975) destaca que “cabe la Geodesia establecer una red de puntos que pueda ser empleada
para control de los levantamientos efectuados”.

Objetivos de la geodesia

Práctico: entregar referencias precisas para el control de levantamientos topográficos

Determinación de la forma y dimensiones de la Tierra (y otros cuerpos celestes)

Definiciones / conceptos:

Ciencia de medir o levantar la Tierra o parte de ella

Ciencia que determina la figura geométrica de la Tierra y su interrelación con puntos

seleccionados en su superficie

Hosmer: Ciencia que trata de las investigaciones de la forma y dimensiones de la superficie de la

Tierra

Zakatov: Estudio de la figura (forma y medidas) y del campo gravitacional exterior de la Tierra

Comittee on Geodesy – EEUU: Es el ramo de la metemática aplicada que determina, por medio de
observaciones y mediciones, la exacta posición de puntos, figuras y áreas de grandes porciones de
la superficie terrestre, la forma y tamaño de la Tierra y las variaciones de la gravedad terrestre

National Research Council – Canadá: Es la disciplina que lida con mediciones y representación de
la Tierra, incluyendo su campo de gravedad, con variaciones en el espacio-tiempo

Comunidad Europea: todas las actividades de evaluación, manejo de tierras, prueba de suelos,
cartografía, levantamientos subterráneos, mapeo nacional, levantamiento de limites y SIG.

OSU: Geodesia es una ciencia interdisciplinar la cual usa mediciones espaciales, aéreas y terrestres
para estudiar la forma y tamaño de la Tierra, los planetas y sus satélites, y sus cambios; para
 10-04-02 [Link]- pag. 4 de 28

determinar precisamente posición y velocidad de puntos y objetos en la superficie u órbita de los

planetas, dentro de un sistema de referencia terrestre y aplicar esos conocimientos a una variedad
de aplicaciones científicas y de ingeniería, usando herramientas de la matemática, física,
astronomía y computación.

Para lograr su objetivo puede valerse de operaciones geométricas realizadas sobre la superficie
terrestre (medidas angulares y de distancias) asociadas a determinaciones astronómicas y
gravimétricas; o más modernamente efectuadas sobre satélites artificiales.

Áreas de la Geodesia: Geodesia Geométrica, Geodesia Física y Geodesia Celeste (incluye la Geodesia
Satelital).

Bajo otro punto de vista: Geodesia Teórica: que estudia el elipsoide y el geoide (y su relación);
Geodesia Aplicada: descripción de la superficie terrestre.

Problemas típicos:

Conocimiento de la forma y dimensión de la Tierra;

Estudio del elipsoide como superficie de referencia;

Resolver problemas geométricos: métodos, fórmulas, aproximaciones;

Representar el elipsoide de acuerdo a sistemas de proyección (cartografía y topografía);

Estudio de redes geodésicas nacionales o globales;

Estudio del campo gravitacional de la Tierra (fuerza de la gravedad y desvío de la vertical);

Conocimiento del geoide (mapas geoidales);

Determinación de alturas y del NMM;

Estudio del movimiento de las placas terrestres;

Procedimientos de terreno para apoyar trabajos de levantamientos;

Establecer referenciales para proyectos de ingeniería.

 10-04-02 [Link]- pag. 5 de 28

Historia de la Geodesia

2400aC: Mapa más antiguo - valle del río Eufrates;

1333aC – 1300aC: catastro del valle del Nilo – Ramses II;

Pitágoras (580 – 500aC) fue el primero a suponer la Tierra como esférica

Aritósteles (384 – 322aC) observó el contorno circular de la sobra de la Tierra proyectada en la

Luna durante los eclipses; estimó el diámetro de la esfera terrestre en 400.000 estadios (84.000 a
63.000km, dependiendo de la conversión).

Arquímedes (287 – 212aC) calculó en 300.000 estadios (63.000 a 47.000km) usando diferente
longitud de estadio

Eratóstenes (276 – 194aC) filósofo y matemático, director de la biblioteca de Alejandría.

Observó en Syene (margen derecha del Nilo) que el Sol cruzaba el meridiano en el cenit y en
Alejandría el Sol causaba una sombra de 1/50 de circunferencia (7º12’). La distancia entre las dos
ciudades es de 5.000 estadios (medido por los “geomensores” reales en días a camello).

Suponiendo (errado) que ambas ciudades están en el mismo meridiano:

∆Z = ∆φ ; d = 5.000 estadios ; π = 256/81 = 3,16 y 1 estadio ≈ 157m

1/50 circulo ----- 5.000 estadios

1 circulo ----- 2πR

R = (5.000 x 50) / 2π = 39.556,96 estadios = 6.210km (error < 2%)

Poseidonius (135 – 50aC). Un siglo después a través de la distancia Alejandría – Rodas recalculó
el radio de la Tierra usando la estrella Canopus, obteniendo un valor semejante.

I-Hing (China siglo IIXaC), monje budista matemático y astrónomo midió un arco de 11.440li (1 li
≈442m) resultando P ≈ 56.700km y R ≈ 9.000km

Edad Media. obscuridad: prohibidos Copérnico, Kepler, Galileo, etc;

Jean Fernel (Francia). En 1525 se midió el arco de 1º entre París y Amiens usando las
revoluciones de una rueda. Fue obtenido 56.746 toesas (= 110.600 metros)
 10-04-02 [Link]- pag. 6 de 28

Jean Picard (Francia 1620 – 1682) introdujo el telescópico para observar alturas de estrellas y
ángulos en la traingulación; midió dos bases con reglas de madera. Calculó que 1º = 57.060 toesas
(1 toesa ≈ 1,95m) 111.210m, R=6.372km

Isaac Newton (Inglaterra 1642 – 1727) se valió de los resultados de Picard para sus estudios sobre
gravitación, considerando la Tierra achatada en los polos

Giovanni Cassini (Francia 1625 – 1712) concluyó que la longitud de un arco de meridiano
disminuye con el aumento de la latitud: achatada en los polos> Demarca el inicio de la Geodesia
Moderna.

1735 con el auspicio de la Academia de París organizan dos mediciones de arcos de 1º: a Perú
(hoy Ecuador) con Pierre Bouguer, La Condamine y Godin resultando 110.613m y, a Laponia con
Clairaut, Maupertuis y Camus, resultando 111.948m. Se adoptó el elipsoide de revolución; a=
6.376,45km y b= 6.355,88km

1790 se crea el metro

1924 la Asamblea General de la Asociación de Geodesia de la Unión de Geodesia y Geofísica

Internacional (UGGI) realizada en Madrid resolvió adoptar el elipsoide de Hayford como de
Referencia Internacional

1953 El IAGS (Servicio Geodésico Interamericano) terminó la triangulación desde México hasta
el sur de Chile

1956 se recomienda para América del Sur el elipsoide de Referencia Internacional. Se adopta el
PSAD56, con punto datum La Canoa (Venezuela) con deflexión de la vertical igual a cero.

1969 la UGGI recomienda para América del Sur el elipsoide de referencia 1967. Lleva a la
definición del SAD69 con elipsoide denominado UGGI67, con punto datum en Chua (Brasil)

1995 se efectúa la primera campaña del proyecto SIRGAS (Sistema de Referencia Geocéntrico
para América del Sur). Red científica medida con GPS, referida a ITRF95,4 que en la práctica es
igual a WGS84. Segunda campaña en 2000.

2003?: Sistema de referencia para Chile, SIRGAS2000.

 10-04-02 [Link]- pag. 7 de 28

UNIDADES DE MEDIDA

¿Que distancia existe entre dos puntos? Depende de las unidades, donde es adoptada y en que época

No hay equivalencia exacta para el estadio en la medición entre Alejandría y Syene. Historiadores
evalúan entre 157,5 y 210m aproximadamente.

Angulares

Egipto: se pensaba que el Sol giraba la Tierra en 360 días. 1 día = 1 grado

GON: circunferencia en 400 GON

Lineales

Estadio

Legua: origen en Galia, Francia. 1 legua = 1,5 millas = 1500 pasos

Milla: (mil) de los militares romanos; 1000 pasos; 1 paso = 5 pies romanos.

Milla marítima: distancia entre dos puntos en la misma longitud y separados por 1’ en latitud (1.852m)

Palma: mayor distancia entre el pulgar y el meñique

Pulgada: Segunda falange del pulgar

Pié: Inglaterra y EEUU; 12 pulgadas

Vara: trozo de madera con 5 palmas de longitud o 16 pies (Inglaterra); el juez poseía la “vara” legal.
En Inglaterra: “ ... a la salida de la iglesia, después del oficio religioso, dieciséis hombres tomados al
azar entre los fieles, altos y bajos, se colocarán en línea recta con sus respectivos pies izquierdos, unos
enseguida de los otros ...”

Metro:

1790 a petición de Luis XVI se introdujo la unidad metro y se depositó una barra de platino.
Un tercio de esta representa 1 pié que será dividido en 10 pulgadas.

1791 Comisión de Pesos y Medidas (Francia), se adopta el método de “fracción de un

cuadrante de meridiano”.

1795: 1 metro = cuadragésima millonésima parte del meridiano terrestre. 1 cuadrante =

5.130.740 toesas; 1 metro = 443,2959 líneas; 1 toesa (Perú) = 6 pies = 72 pulgadas = 864
líneas = 1,949 metros; 1m = 0,5130740 toesas.
 10-04-02 [Link]- pag. 8 de 28

1789: fueron fabricadas 4 barras bimetálicas en capas de cobre y platino

1870: primera tentativa internacional con la creación del “Bureau Internationale des Poids et
Mésures”

1889: 30 copias fueron hechas y distribuidas a diferentes países

1890: surge el patrón natural en función de la longitud de onda de la radiación cadmio rojo

1960: redefinición como la longitud de onda de la luz – 1.650.763,73λ del gas cripton-86 en el
vacío. Precisión 4 partes en 109

1983: La Conferencia General de Pesos y Medidas en París redefine en función del tiempo. La
longitud que viaja la luz en el vacío durante 1/299.792.458 segundos. Precisión 1 parte en 1010

Notación: Múltiplos y submúltiplos:

Metro: m 10-6 : µ (micro)
Gramo: g 10-3 : m (mili)
Segundo: s 103 : k (kilo)
106 : M (mega)

1 metro = 36,3636 pulgadas

1 vara Chile = 0,835 metros
1 milla terrestre = 1.609,31 metros
1 milla marítima = 1.851,85 metros
1 legua marítima = 5.555,55 metros
1 legua métrica = 5.500 metros
1 pié ingles = 0,30479 metros
1 yarda = 0,91438 metros

PROBLEMAS GEODESICOS Y TOPOGRÁFICOS PRÁCTICOS.

Transportar coordenadas: problema directo (α,D) (X,Y)

Calcular distancia y acimut: problema inverso (X,Y) (α,D)

Determinar forma de la Tierra (geoide)

¿?
 10-04-02 [Link]- pag. 9 de 28

X B = X C + dBC sen α dBC = ∆X 2 + ∆Y 2

YB = YC + dBC cos α ∆X
α = arctg
∆Y

a b
En triangulación: Ley de los senos =
sen  sen B̂

∆h = di ⋅ cos DZ + i − s
∆h = di ⋅ sen α + i − s
∆h AB = ∑ atrás − ∑ adelante
∆h = dh ⋅ cot DZ + i − s
∆h = dh ⋅ tgα + i − s

En geodesia las proyecciones de las distancias son más grandes y son referidas al elipsoide, por lo
tanto necesitan algunas reducciones.

COORDENADAS ASTRONÓMICAS

Todos los cuerpos en la Tierra están sujetos al campo gravitacional, resultante de la fuerza de
atracción ejercida por la Tierra y la fuerza centrífuga.

Superficies equipotenciales W=constante, denominados “geopes”.

 10-04-02 [Link]- pag. 10 de 28

Líneas de fuerza perpendiculares a los geopes: líneas de fuerza de campo = verticales, representa la
dirección del vector gravedad (eje de plomo o eje principal del teodolito).

Latitud astronómica: ángulo entre la vertical e su proyección ecuatorial.

Meridiano astronómico: plano vertical paralelo al eje de rotación terrestre

Longitud astronómica: ángulo diedro entre el y el meridiano astronómico y el meridiano medio

astronómico de Greenwich (origen).

Por consecuencia del movimiento de los polos terrestres que alteran el eje de rotación y
consecuentemente del ecuador, las coordenadas astronómicas son función del tiempo. Deben ser
reducidas a una misma época.

SUPERFICIES DE REFERENCIA

En geodesia se trabajan en 3 superficies

1. Superficie física terrestre: donde se realizan las operaciones de medida

2. Superficie del modelo geométrico de referencia, elipsoide de revolución: donde se realizan los
cálculos geodésicos

3. Geoide, superficie que representa la forma real de la Tierra en función de su campo gravitacional;
es una superficie equipotencial; un geope que más se aproxima al Nivel Medio del Mar (NMM);
coincide con la superficie de los océanos en reposo extendida idealmente sobre los continentes; es
una superficie “horizontal”; es el origen para las altitudes o altura ortométrica (distancia por la
vertical de un punto al geoide). Se obtiene por nivelación geométrica asociada a gravimetría.

Uno de los problemas geodésicos más importantes y complejos es la determinación de la separación

entre geoide y elipsoide (ondulación geoidal)
 10-04-02 [Link]- pag. 11 de 28

z de
Ge
lipsoi P
óid E
e SUPERFICIE
H TERRESTRE

GEOIDE
y N h
x
ELIPSOIDE
h=H+N

s
icale
vert

sup
equip erficie
otenc
ial
ión
trac
cen sa
con de ma

(tarea: investigar las 2 superficies de referencia, elipsoide y geoide)

GEOMETRÍA DEL ELIPSOIDE

P
Q
Elipsoide de revolución: cuerpo geométrico a
generado por la rotación de una elipse
alrededor del eje menor, el eje menor a

coincide con el eje polar terrestre.

o φ
d d
FQ + F´Q = constante = 2·a F F´

En el elipsoide tri-axial: a=c=b esfera

c=b elipsoide de revolución

El elipsoide de revolución es la “forma matemática de la Tierra”, donde se realizan los cálculos

x 2 + y2 z2
Ecuación del elipsoide de revolución: =1 +
a2 b2
Excentricidad es la distancia focal expresada en términos del semi eje mayor (a)
 10-04-02 [Link]- pag. 12 de 28

FO d
e= =
a a
2 d 2 a 2 − b2 b2
e = = = 1−
a2 a2 a2

Parámetros:
a−b b 1
Achatamiento (α): f= =1− α=
a a f

a 2 a 2 − b2 b2
1 excentricidad (e): e = = 1−
a2 a2

a 2 − b2 a2
2a excentricidad( e’): e' 2 = = −1
b2 b2

Otras relaciones:
e' 2 e2
e2 = 2 ⋅ f − f 2 e2 = e' 2 =
1 + e' 2 1 − e2

b
a= = b ⋅ 1 + e' 2 b = a ⋅ (1 − f ) = a ⋅ 1 − e 2 b 2 = a 2 ⋅ (1 − e 2 )
1− f

COORDENADAS GEODÉSICAS

Eje Polar
Basadas en un elipsoide de revolución P1
Sección 1er
Z
generado por una elipse girada en h Vertical
Meridiano Origen
torno al eje polar; es el modelo P
Sección
b
matemático de la Tierra. Meridiana

Sección Normal: sección que

a o R
φ Z1
Y
contiene la normal al elipsoide en
Q x1
P λ Plano Ecuatorial

Sección Meridiana: sección y1

X
normal particular, contiene el eje
menor (polar)
Sección 1º vertical: perpendicular a la sección meridiana en P
 10-04-02 [Link]- pag. 13 de 28

Gran Normal: segmento PQ de la normal; desde P hasta el eje polar

Pequeña Normal: segmento PR, hasta el plano ecuatorial
Meridiano Geodésico: intersección de la sección meridiana con el elipsoide
Paralelo Geodésico: intersección de un plano paralelo al ecuador y el elipsoide, es un círculo
Latitud Geodésica: ángulo formado por la normal en P y su proyección en el ecuador; (-) al sur
del ecuador; varía de +90º a -90º
Longitud Geodésica: ángulo formado entre el meridiano origen y la sección meridiana en P; (-) al
este de Greenwich; varía 0º a 360º o a +/-180º
Altura Geométrica o Elipsoidica: distancia por la normal entre el elipsoide (P) y el punto P1

Desvío de la Vertical (δ): ángulo entre la vertical local (en P1) y la normal al elipsoide; ayuda a
transformar grandezas astronómicas a geodésicas:

Componente meridiana ξ = φa – φ

Componente 1º vertical η = (λa – λ) cos φ = (Aa – A) cot φ

Ecuación de Laplace: A = Aa – (λa – λ) sen φ

Usada en astronomía geodésica para orientar redes geodésicas. En vértices de triangulación que se
realizan determinaciones astronómicas de azimut y longitud, se llaman “puntos de Laplace”

PSAD-56 SAD-69 WGS-84

Elipsoide Internacional 24 (Hayford) UGGI67 (GRS-67) GRS-80
a 6378388 6378160 6378137
1/f 297 298.25 298.257223563
b 6356911.946 6356774.719 6356752.314
e^2 0.00672267002 0.00669454185 0.00669437999
e´^2 0.00676817020 0.00673966080 0.00673949674

(tarea: investigar los sistemas PSAD56, SAD69 y WGS84)

 10-04-02 [Link]- pag. 14 de 28

RADIOS DE CURVATURA DE SECCIONES NORMALES

Gran Normal (N): distancia normal al elipsoide entre el punto (M) y la intersección con eje Z (H)

Pequeña Normal (N’): distancia normal al elipsoide entre el punto (M) y la intersección con ecuador

∆τ var iación de dirección de la tan gente τ

Curvatura = =
∆S var iación del arco S
∆τ 1
lim =K y Radio de curvatura R=
∆S → 0 ∆S K
dz d 2z
pero τ= y ∆τ =
dx dx 2
3
 dz 2  2
1 + ( dx ) 
Para una curva plana z = f(x), el radio de curvatura es: R = (*)
d 2z
dx 2
x2 z2
Elipse meridiana: + =1 /·a 2 b 2
2 2
a b
b2x 2 + a 2z 2 = a 2b2
b2 + x 2 a 2b2
z2 + =
a2 a2
b2
pero = (1 − e 2 )
2
a
z 2 + (1 − e 2 ) x 2 = a 2 (1 − e 2 )

2z·dz + (1 − e 2 )2 x·dx = 0
derivando: z·dz + (1 − e 2 ) x·dx = 0
dz x
donde : = −(1 − e 2 )
dx z
 10-04-02 [Link]- pag. 15 de 28

dz
diferenciando :
dx
dz
z−x
d 2z dx
= −(1 − e 2 )
2
dx z2
dz
reemplazando
dx
x2
2 z + (1 − e 2 ) 2 2 2
d z
= −(1 − e 2 ) z = −(1 − e 2 ) z + (1 − e ) x
2
dx z2 z3
pero : z 2 + (1 − e 2 ) x 2 = a 2 (1 − e 2 )
d 2z 2 a
2
= −(1 − e )
dx 2 z3

dz d 2z
reemplazando y en R :
dx dx 2
a ⋅ (1 − e 2 )
Designando por M el radio de curvatura M= 3
(1 − e 2 ⋅ sen 2φ) 2
En la elipse:
2 a2
después de algunos reemplazos, finalmente x =
1 + (1 − e 2 ) ⋅ tg 2 φ

a 2 ⋅ cos 2 φ a ⋅ cos φ a ⋅ (1 − e 2 ) ⋅ senφ

x2 = x = y z=
1 − e 2 ⋅ sen 2 φ 1 − e 2 ⋅ sen 2 φ 1 − e 2 ⋅ sen 2 φ

De la figura:

x = N cosφ y z = N’ senφ

a a ⋅ (1 − e 2 )
N= N' =
1 − e 2 ⋅ sen 2 φ 1 − e 2 ⋅ sen 2 φ
N
z=N´sen φ

φ Nsen φ
N' = N(1 − e 2 )
φ
x=Ncos φ
 10-04-02 [Link]- pag. 16 de 28

Secciones principales (para un punto):

Sección meridiana, radio de curvatura mínimo
Sección 1o vertical (acimut 90º), radio de curvatura máximo

Radio de curvatura de la sección meridiana (M):

a ⋅ (1 − e 2 )
M= 3
(1 − e 2 ⋅ sen 2φ) 2

Radio de curvatura de una sección normal cualquiera con acimut α (Rα):

Teorema de Euler:
1 cos 2 α sen 2 α Rα =
N⋅M
= +
Rα M N N ⋅ cos α + M ⋅ sen 2 α
2

Radio de curvatura del 1o vertical - α = 90º R 1o = N

b
Radio medio de curvatura (Ro) Ro = M ⋅ N =
1 − e ⋅ sen 2 φ
2

a2
En los polos: M P = N P = R oP =
b

b2
En el ecuador: ME = NE = a R oE = b
a

Radio de un paralelo (r): r = N ⋅ cosφ

r tiene valor máximo en el ecuador (=a) y nulo en los polos

 10-04-02 [Link]- pag. 17 de 28

LATITUDES GEOCÉNTRICA Y REDUCIDA

M’
En los problemas prácticos de la geodesia P
interfiere solo la latitud geodésica, pero en M
aspectos teóricos son útiles otros dos tipos de b ρ
latitud: z
Latitud geocéntrica (ψ): ángulo entre el a O µ ψ ϕ
Q Q’
radio vector de un punto M con su x
proyección en el ecuador;
Latitud reducida (µ): ángulo formado por el
radio (MO) y su proyección en el ecuador; H
MO formado por la prolongación de la P’ ϕ : latitud geodésica
ordenada en M, hasta la circunferencia ψ : latitud geocéntrica
circunscrita de radio “a”. µ : latitud reducida

Z = N’ senφ = x tgψ

N’= N(1-e2) y x = N cosφ

N(1-e2) senφ = N cosφ tgψ

x = a cosµ y z = b senµ

z b ⋅ senµ b
tgψ = = = ⋅ tg ⋅ µ = (1 − e 2 ) ⋅ tgµ tgψ = (1 − e 2 ) ⋅ tgφ
x a ⋅ cos µ a

tgµ = (1 − e 2 ) ⋅ tgφ

VALORES NUMERICOS

1. El cuadrado de la excentricidad es un valor del mismo orden que el achatamiento:

e2 = 0,006 ....

2. Dos elipsoides con el mismo achatamiento y 1 metro de diferencia en “a”, resulta en 3,5mm de
diferencia en “b”

3. Un arco de meridiano de 1º varía entre 110.567m (ecuador) y 111.699m (polos); variación de

~1.152m; los marinos simplifican para 1’ = 1.852m = 1 milla marítima
 10-04-02 [Link]- pag. 18 de 28

4. Los radios de la sección del 1º vertical son, para un punto, siempre más grandes; N ≥ a ; M ≤ a
2
5. Un cuadrado elipsóidico de 1º x 1º en el ecuador tiene área aproximada de 12.300km y en la

región central de Chile es aproximadamente 10.300km2.

LONGITUD DE UN ARCO DE ELIPSE MERIDIANA

[Geodesia Geométrica, DMA 1982, Richard Rapp]

S
Para el caso de un arco circular: S = R·α
α
R

Arco MM’ de la elipse meridiana . Radio de curvatura no varía.

y Z
P’ds ∆λ
P
C
b ∆φ B
D
A

a O φ x a o Y
Q
M

b
X

M: radio de curvatura de la sección meridiana

ds = M ⋅ dφ ds = dx 2 + dy 2

φ2 2 φ2 2 2 −3 / 2
∫φ1 ds = a ⋅ (1 − e ) ⋅ ∫φ1 (1 − e ⋅ sen φ) dφ

Haciendo: W = (1 − e 2 ⋅ sen 2φ)

φ2 2 φ2 1
∫φ1 ds = a ⋅ (1 − e ) ⋅ ∫φ1 dφ
W3

Usando el desarrollo en serie de McLaurin:

1 3 2 2 15 4 4 35 6 6 315 8 8
= 1 + e sen φ + e sen φ + e sen φ + e sen φ + ...
W3 2 8 16 128

Se reemplazan las potencias de senφ por ángulos múltiples:

 10-04-02 [Link]- pag. 19 de 28

1 1
sen 2φ = − cos 2φ
2 2
3 1 1
sen 4φ = − cos 2φ + cos 4φ
8 2 8
15 15 3 1
sen 6φ = − cos 2φ + cos 4φ − cos 6φ
16 32 16 32
35 7 7 1 1
sen 8φ = − cos 2φ + cos 4φ − cos 6φ + cos 8φ
128 16 32 16 128

1 1
s = a ⋅ (1 − e 2 ) ⋅ [A ⋅ (φ 2 − φ1) −
⋅ B ⋅ (sen 2φ 2 − sen 2φ1 ) + ⋅ C ⋅ (sen 4φ 2 − sen 4φ1 ) −
2 4
1 1 1
− ⋅ D ⋅ (sen 6φ 2 − sen 6φ1 ) + ⋅ E ⋅ (sen8φ 2 − sen8φ1 ) − ⋅ F ⋅ (sen10φ 2 − sen10φ1 )] + ...
6 8 10

3 2 45 4 176 6 11025 8 43659 10

A = 1+ ⋅e + ⋅e + ⋅e ⋅e + ⋅ e + ...
4 64 256 16384 65536
3 15 4 525 6 2205 8 72765 10
B = ⋅ e2 + ⋅e + ⋅e + ⋅e + ⋅ e + ...
4 16 512 2048 65536
15 4 105 6 2205 8 10395 10
C= ⋅e + ⋅e + ⋅e + ⋅ e + ...
64 256 4096 16384
35 6 315 8 31185 10
D= ⋅e + ⋅e + ⋅ e + ...
512 2048 131072
315 3465 10
E= ⋅ e8 + ⋅ e + ...
16384 65536
693
F= ⋅ e10 + ...
131072

Para el cuadrante meridiano: φ1 = 0º ; φ2 = 90º s = a(1-e2) A π/2

Para SAD69 s = 10.002.001,23m

[Rapp]

Zakatov en 1962:

[ 1 3 3 3 15 1
s = a ⋅ ∆φ ⋅ 1− e2 ⋅ ( + ⋅ cos2φm ) − e4 ⋅ ( + ⋅ cos2φm − ⋅ cos4φm ) + ⋅ e2 ⋅ ∆φ2 ⋅ cos2φm ]
4 4 64 16 64 8

Se considera exacta para líneas hasta 600km; con restricciones se puede usar hasta 2000km.

Zakatov simplificada:

1
Mm = radio de curvatura de la latitud media. s = M m ⋅ ∆φ ⋅ [1 + ⋅ e 2 ⋅ ∆φ2 ⋅ cos 2φm ]
8

Precisión 1mm hasta 400km.

Hasta 45km: s = M m ⋅ ∆φ
 10-04-02 [Link]- pag. 20 de 28

LONGITUD DE UN ARCO PARALELO

Puntos de longitudes λ1 y λ2 en el mismo paralelo, sea L el arco:

r = N cosφ

L = r ⋅ ∆λ = N ⋅ cos φ ⋅ ∆λ

(Calcular la distancia por el paraleloo desde el Meridiano Greenwich a Santiago)

(tarea: calcular y graficar 1” de arco meridiano y paralelo para diferentes latitudes en Chile)

ÁREA DE UN CUADRILÁTERO ELIPSOIDICO

Considerar el área en el elipsoide limitada por meridianos y paralelos conocidos (dφ y dλ).

AB = CD = M dφ

AD = BC = N cosφ dλ

Ärea diferencial: dA = AB * AD = M N cosφ dφ dλ

φ2 λ 2 φ2
A = ∫ ∫ M ⋅ N ⋅ cos ϕ ⋅ dϕ ⋅ dλ = (λ 2 − λ1 ) ∫ M ⋅ N ⋅ cos ϕ ⋅ dϕ
φ1 λ1 φ1

Área de la zona elipsoidica (dφ x 2π)

A ' = 1 + 1 e 2 + 3 e 4 + 5 e 6 + 35 e8 + 63 e10
2 8 16 128 256
B' = 1 e 2 + 3 e 4 + 3 e 6 + 35 e8 + 45 e10
6 16 16 192 256
C' = 3 e 4 + 1 e 6 + 5 e8 + 45 e10
80 16 64 512
1 6 5 8 15 10
D' = e + e + e
112 256 512
E ' = 5 e8 + 3 e10
2304 512

φ +φ
∆φ = φ 2 − φ1 y φm = 2 1
2

Área del cuadrilátero elipsoidico (dφ x dλ)

A = 2 ⋅ b 2 ⋅ ∆ λ[ A '⋅sen ∆ φ ⋅ cos φ m − B '⋅sen 3∆ φ ⋅ cos 3φ m + C '⋅sen 5 ∆ φ ⋅ cos 5φ − − ...]

 10-04-02 [Link]- pag. 21 de 28

φ +φ
∆φ = φ 2 − φ1 ; ∆λ = λ 2 − λ1 yφ m = 2 1
2

APROXIMACIÓN ESFÉRICA

En ciertos problemas la aproximación esférica (considerar la Tierra como esfera) puede ser suficiente,
para triángulos geodésicos pequeños.

Se adopta una familia de esferas con radios entre b2/a y a2/b, que son los radios medio de curvatura
en el ecuador y en los polos, respectivamente.

A cada triángulo corresponde un radio R 0 = M ⋅ N calculado en función de la latitud media del

triángulo.

2⋅a + b
Radio de esfera con media aritmética de los 3 ejes: R=
3

Radio de una esfera de igual área que el elipsoide (RA):

AE 2 3 4 5
RA = = b ⋅ (1 + ⋅ e 2 + ⋅ e 4 + ⋅ e 6 + ⋅ e8 + ...)
4π 3 5 7 9

Radio de una esfera de igual volumen que el elipsoide (RV):

4 4
VESFERA = ⋅ π ⋅ R 3V VELIPSOIDE = ⋅π⋅a2 ⋅b
3 3

3
R V = a 2 ⋅ b = a 6 (1 − e 2 )

CURVAS SOBRE EL ELIPSOIDE

SECCIONES NORMALES RECÍPROCAS

Sección normal directa respecto al punto “A”: sección normal en A que contiene el punto.

Sección normal recíproca respecto al punto “A”: sección normal en C que contiene el punto.

En general, para puntos distintos, las normales en A y C no son coplanares secciones normales
directa e inversa no son coplanares “camino” normal A-C ≠ “camino” normal C-A.

Coplanares solo si los puntos están en la misma latitud o misma longitud.

 10-04-02 [Link]- pag. 22 de 28

Se fuera posible calar con un teodolito, instalado en el elipsoide según la normal, los planos de
observación A-C es diferente a C-A, o sea, son diferentes direcciones.

Para punto más al sur curva directa más al sur.

Secciones normales no definen un triángulo geodésico.

El mejor camino entre los dos puntos es una curva, generalmente reversa, comprendida entre los
planos directo y recíproco, denominada línea geodésica.

SEPARACIÓN ENTRE SECCIONES NORMALES RECÍPROCAS

Z N’

C
β’ N
B
β
b

D
o A
a Y φ
φC
φA
φ’
X

Sección normal C-A

C
Sección normal A-C

Ángulo entre
planos normales

cuerda A-C
A

En la práctica interesan las diferencias en distancia y acimut entre secciones normales recíprocas.

Ángulo auxiliar (β) [Gemael]

e 2 ⋅ ( N1 ⋅ senφ1 − N ⋅ senφ) ⋅ cos φ

tgβ =
N + e 2 ⋅ ( N1 ⋅ senφ1 − N ⋅ senφ) ⋅ senφ

Ángulo ortogonal (V):

 10-04-02 [Link]- pag. 23 de 28

[Gemael] V = β ⋅ senA

A: acimut

Z: ángulo cenital

1 2 S
[Rapp]: V= ⋅ e ( ) ⋅ cos 2 φ m ⋅ sen 2A
2 N1

S: distancia geodésica

A: acimut

Para S = 100km; φm = 45º; A = 45º : V = 6” (valor máximo en A = 45º)

Separación acimutal (θ):

N B

A
A’
S

2θ/3
θ/3

A1

Ángulo en el plano tangente (horizontal) en N

[Gemael]: θ = β ⋅ senA ⋅ cot gZ

e2 S 2
[Rapp]: θ= ⋅ ( ) ⋅ cos 2 φ m ⋅ sen 2A
4 N1

e2 S 2 tgφ1 S
[Jordan]: θ= ⋅ ( ) ⋅ cos 2 φ1 ⋅ senA12 ⋅ (cos A12 − ⋅ )
2 N1 2 N1

Para φm = 0º y A = 45º
S 200km 100km 50km
θ” 0,36” 0,09” 0,023”
 10-04-02 [Link]- pag. 24 de 28

Para φm = 52º y A = 45º

S 150km 100km 30km
θ” 0,057” 0,032” 0,003”

En la práctica se hacen correcciones a distancias > 30km

SEPARACIÓN ENTRE ARCOS

En el punto medio entre A y C, la separación “L” será máxima:

e 2 ⋅ S3 ⋅ cos 2 φ ⋅ sen 2A
[Gemael]: L=
16 ⋅ N 2

Para φm = 45º y A = 45º

S 200km 100km 50km
L máximo 0,050m 0,006m 0,0008m
Para φm = 52º y A = 45º
S 150km 100km 30km

L máximo 0,013m 0,0038m 0,0001m

Corrección del acimut debido a la altura del punto observado

Las direcciones se miden entre puntos sobre la superficie terrestre, sin embargo los cálculos se
efectúan sobre la superficie del elipsoide, por lo tanto existe influencia de la altura del punto visado en
el acimut calculado.

Desde A se cala B, a altura h.

Acimut deseado: A;

Acimut observado: Ah

Puesto que el elipsoide es achatado, se debe considerar la diferencia (A – Ah)

 10-04-02 [Link]- pag. 25 de 28

Ah A b’
b δ2
A

S
[Rapp] δ2 = ⋅ e 2 ⋅ cos 2 φ m ⋅ cos A AB
Mm

h
δ" = A − Ah = ⋅ e2 ⋅ cos2 φm ⋅ sen2AAB ⋅ ρ"
2Mm

La corrección no depende de la distancia entre los puntos. Esta corrección se aplica solamente
cuando el punto calado (B) está en altura, independiente si el punto origen (A) de las visadas está en
altura.

(tarea: hacer gráfico de corrección por altura, para diferentes latitudes)

LÍNEA GEODÉSICA

B meridiano
normal B-A
C

tangente
B

Geodésica entre
AyB
A
A normal A-B
 10-04-02 [Link]- pag. 26 de 28

Para obtener un único triángulo elipsóidico, los vértices deben estar conectados por líneas
geodésicas.

Línea geodésica, yacente a una superficie, es la que en todos sus puntos el plano osculador es normal a
la superficie

Es única entre dos puntos

Es la distancia más corta sobre la superficie

En el plano es una recta

En una esfera es un arco de círculo máximo

En el elipsoide es reversa (curvatura espacial) ; no es plana

Propiedad importante (fig. 12): la normal principal de la geodésica coincide, en cualquier punto, con la
normal del elipsoide. La normal (principal) está contenida dentro del plano osculador que pasa por
tres puntos infinitamente cercanos en la curva. La sección normal no tiene esta propiedad.

(¿Paralelos y meridianos son líneas geodésicas?)

característica: r senA = constante

Correcciónn ángulo geodésica – sección normal.

Se necesita transformar el acimut de una sección normal en el acimut correspondiente de la línea

geodésica, ya que esta representa sin ambigüedad el lado de un triángulo geodésico.

La línea geodésica divide el ángulo “θ” de las secciones normales recíprocas (excepto en los casos de
acimut 90º o 270º), en razón 1:2, estando siempre más cerca de la sección normal directa. Designando
por “τ” la corrección:
 10-04-02 [Link]- pag. 27 de 28

N B
r = N cos φ
A N cosφ sen A = cte = K
A’
S
Ecuador A = 90º K=N
Meridiano A = 0 K=0

2θ/3
θ/3

A1

θ" e 2 ⋅ S2 S ⋅ sen 2φ ⋅ senA

τ" = = A'− A = ⋅ [cos 2 φ ⋅ sen 2A − ] ⋅ ρ"
3 12 N 2 2N

" θ" e 2 ⋅ S2
τ = = A'−A = ⋅ cos 2 φ ⋅ sen 2A ⋅ ρ"
3 2
12 N

(tarea: calcular la reducción anterior para diferentes distancias en azimut 45º)

Si los dos puntos están sobre el mismo meridiano, solo existe una sección normal entre ellos, la
geodésica es el meridiano (τ = 0). Pero .....
Si los dos puntos están sobre el mismo paralelo, solo existe una sección normal entre ellos, sin
embargo la geodésica NO coincide con esta sección normal. La geodésica estará fuera de las
secciones normales. La diferencia no influye en la distancia.

Corrección a la deflexión de la vertical.

(cenit geodésico)
Z ξ
e meridiano
η
Z’

Debido a que las mediciones no se efectúan sobre el elipsoide y sí sobre la Tierra verdadera, bajo
influencia del campo gravitacional, los ángulos se miden en la horizontal local (perpendicular a la
vertical) y deben ser llevados al plano perpendicular a la normal del elipsoide.
La diferencia se llama deflexión de la vertical, con componentes ξ (componente meridional) y η
(componente en el 1er vertical).

1
δ3 = − (ξsenAAB − η cos A AB ) ⋅ ⋅ ρ"
tgZ
 10-04-02 [Link]- pag. 28 de 28

Z: distancia cenital del punto observado = (90-φ)

1
δ3 = − (ξsenAAB − ηcosAAB) ⋅ ⋅ ρ"
tg(90− φ)

Esta corrección es normalmente muy pequeña

Corrección total:

Ac =Ao +δ+τ+δ3
 10-04-02 [Link]- pag. 1 de 11

COORDENADAS CARTESIANAS.

Comenzaron a tener uso más amplio con la geodesia satelital y los sistemas de referencia globales.

Los cálculos 3-dimensionales se facilitan, pero tiene el inconveniente que es no es apto a la

cartografía.

h
P1
Z
h
P
b
N N´
(N + h)
a o R
φ Z1
Y (N + h)senφ (N´+h)senφ
Q X1
λ φ
Y1
X φ
(N + h) cos φ

Relación entre coordenadas geodésicas y cartesianas.

d = ( N + h ) cos φ = X 2 + Y 2
X = d ⋅ cos λ
Y = d ⋅ senλ
N' = N ⋅ (1 − e 2 )

N: Gran Normal
N’: Pequeña Normal

X = (N + h) ⋅ cosφ⋅ cosλ
Y = (N + h) ⋅ cosφ⋅ senλ
Z = ( N ⋅ (1 − e2 ) + h) ⋅senφ

Fórmulas directas de Bowring.

Z1 + b ⋅ e'2 ⋅sen 3 ψ a 2 − b2
φ = arctan[ ] 2 a exentricidad : e'2 =
d − a ⋅ e 2 ⋅ cos3 ψ b2
Y valores auxiliares :
λ = arctan[ 1 ]
X1 d = (X 2
+ Y 2)
d
h = −N a⋅Z
cosφ ψ = arctan[ ]
b⋅d
 10-04-02 [Link]- pag. 2 de 11

SISTEMAS DE REFERENCIA CONVENCIONAL

Sistema de Referencia Celeste Convencional (CCRS): Eje Xc apunta al equinoccio vernal medio de las
12h del 1º de enero de 2000 (día Juliano 2451545,0 – J2000); eje Zc apunta en la dirección del polo
norte celeste medio de la misma época; eje Yc completa el sistema dextrógiro.

Sistema de Referencia Terrestre Convencional (CTRS):

Geocéntrico; centro de masa de la Tierra (Tierra y atmósfera) coincide con el origen
Fijo a la Tierra - ECEF
Orientación dada por el BIH (Bureau International de L’Heure) en la época 1984,0
Sin rotación
Eje Z en la dirección del polo terrestre convencional (CTP); eje X en la dirección del meridiano medio
de Greenwich . Se recomienda usar el elipsoide GRS80. El CTRS es definido como ITRF
(International Terrestrial Reference Frame) el cual es mantenido por el IRES (International Earth
Rotation Service)

La transformación entre CCRS y CTRS se efectúa usando rotaciones que consideran precesión,
nutación rotación y orientación de la Tierra (incluyendo el movimiento del polo)
Precesión: movimiento secular cónico del eje de rotación respecto a la eclíptica
Nutación: movimiento del eje de rotación respecto del eje de la figura; es parte del movimiento del
polo

(tarea: investigar el sistema ITRF)

SISTEMAS DE REFERENCIA.

Sistema Geodésico: adopta un elipsoide de referencia fijado espacialmente respecto al cuerpo

terrestre.

Los sistemas de referencia continentales o nacionales no son geocéntricos y a veces no paralelos al

CTS.
 10-04-02 [Link]- pag. 3 de 11

Un sistema de referencia puede ser :

DEFINIDO: sistema abstracto

- IDEAL: para no ser implementado

- CONVENCIONAL: asociado a la implementación (constantes y modelos físicos), ejemplo es

el ITRF

REALIZADO: asume características físicas. Un sistema realizado es una Red de Referencia.

La realización no siempre corresponde a la definición. La realización depende de las técnicas

utilizadas. Ejemplo, SAD-69 es definido de forma única y realizado de formas diferentes.

Se puede definir un sistema local pero al realizarlo se introducen errores. Al relacionar estos pueden
aparecer rotaciones producto de las deformaciones naturales de la realización.

TRANSFORMACIÓN ENTRE SISTEMAS GEODÉSICOS.

La transformación de sistemas se puede realizar de varias formas.

Parámetros geocéntricos

- 3 parámetros: solo translación (sistemas paralelos)

- 6 parámetros: translación (3) y rotación (3)

- 7 parámetros: translación, rotación y escala

- n parámetros: distorsiones

Modelos diferenciales: ej. Rapp, Veis, Molodenskii, simplificadas de Molodenskii, Vincenty.

Modelos Cartesianos (o matriciales): más fáciles de “comprender”

 10-04-02 [Link]- pag. 4 de 11

Ángulos de Euler.

Z1 Z2 Z1 Z3 Z2
ω
ε ψ
Y3
Y2
ψ Y3
ω Y1 ε Y2
X1 ω ε
X2 X2 ψ
X2 X3

X 1  X 
 Y  = R (ω) ⋅  Y 
 1 3  
 Z1   Z 

X 2  X 1  X 
 Y  = R (ε) ⋅  Y  = R (ε) ⋅ R (ω) ⋅ Y 
 2 1  1 1 3  
 Z 2   Z1   Z 

X 3  X 2  X 
 Y  = R (ψ ) ⋅  Y  = R (ψ ) ⋅ R (ε) ⋅ R (ω) ⋅  Y 
 3 2  2 2 1 3  
 Z 3   Z 2   Z 

 cos ω sen ω 0 1 0 0  cos ψ 0 − sen ψ 

R 3 (ω) = − sen ω cos ω 0 R1(ε) = 0 cos ε sen ε R 2 (ψ ) =  0 1 0 
 0 0 1 0 − sen ε cos ε  sen ψ 0 cos ψ 

cos ω ⋅ cos ψ − sen ω ⋅ sen ε ⋅ sen ψ sen ω ⋅ cos ψ + cos ω ⋅ sen ε ⋅ sen ψ − cos ε ⋅ sen ψ 
R 3 ( ω) ⋅ R1( ε) ⋅ R 2 ( ψ ) =  − sen ω ⋅ cos ε cos ω ⋅ cos ε sen ε 
 cos ω ⋅ sen ψ + sen ω ⋅ sen ε ⋅ cos ψ sen ω ⋅ sen ψ − cos ω ⋅ sen ε ⋅ cos ψ cos ε ⋅ cos ψ 
 

CONSIDERANDO QUE LAS ROTACIONES SON PEQUEÑAS:

sen(α) = α ; cos(α) = 1 ; sen(α)∗sen(α) = 0

 1 ω − ψ

E = − ω 1 ε 
 ψ − ε 1 

X 3  TX   1 ω − ψ  X1 
 Y  = TY  + k ⋅ − ω 1 ε   Y1 
 3   
 Z 3   TZ   ψ − ε 1   Z1 

Generalizando:
 10-04-02 [Link]- pag. 5 de 11

X 3  X 
 Y  = Y  y k = 1 + ∆L
 3  
 Z 3   Z 

X = ∆X + X + ω ⋅ Y − ψ ⋅ Z + X ⋅ ∆L
Y = ∆Y + Y − ω ⋅ X + ε ⋅ Z + Y ⋅ ∆L
Z = ∆Z + Z + ψ ⋅ X − ε ⋅ Y + Z ⋅ ∆L

Este corresponde al modelo completo de 7 parámetros:

3 traslaciones de origen (TX, TY, TZ)

3 rotaciones (ω, ε, ψ)

escala (k)

Afortunadamente los sistemas geodésicos son (generalmente) definidos y realizados (casi) paralelos al
los sistemas convencionales; se utilizan mediciones modernas que no introducen escala. Se eliminan
las rotaciones y escala.

X 3  TX  X1 
 Y  = TY  +  Y 
 3     1
 Z 3   TZ   Z1 

SISTEMAS GLOBALES DE REFERENCIA.

El posicionamiento con GPS, así como su homólogo ruso GLONASS (GLobal NAvigation
Satellite System), requiere sistemas de referencia bien definidos y consistentes, globales y
geocéntricos, esto implica que consideran todo el globo terrestre y tienen su origen en el centro de
masa de la Tierra.

Sistema ITRF.

El Sistema de Referencia Terrestre Internacional – ITRF (International Terrestrial

Reference Frame), materializa un sistema global de carácter científico establecido por el Servicio
Internacional de Rotación Terrestre - IERS (International Earth Rotation Service) y está
materializado por redes geodésicas continentales implantadas a través de técnicas geodésicas
espaciales modernas.
 10-04-02 [Link]- pag. 6 de 11

Debido a la precisión alcanzada en la implantación y a los movimientos tectónicos

sufridos en la corteza terrestre, las coordenadas asignadas a las estaciones deben ser reducidas a una
época de referencia común t0. Significa esto la puesta en práctica de la geodesia global 4D (tetra-
dimensional), donde a las coordenadas geocéntricas 3D les son asignadas sus variaciones o
velocidades, o sea, las coordenadas pasan a tener validez respecto a una determinada época.

SISTEMA WGS-84.

El Sistema Geodésico Mundial 1984 – WGS-84 (WorldGeodetic System 1984), es el

sistema de referencia para el GPS y compatible con un ITRF básicamente bajo los siguientes aspectos:

Posición: geocéntrico, con origen en el centro de masa de la Tierra, incluyendo océanos y

atmósfera;

Orientación:

- eje Z en la dirección del Polo de Referencia IERS;

- eje X en la intersección del Meridiano de Referencia IERS y el plano ecuatorial;

- eje Y completa el sistema ortogonal dextrógiro (sentido mano derecha).

Al sistema cartesiano se asigna un elipsoide denominado también de WGS-84. Este

elipsoide posee los parámetros del Sistema Geodésico de Referencia 1980 – GRS-80.

Recientes refinamientos del WGS-84 han llevado a la realización del denominado WGS-
84 (G730*) en concordancia con el ITRF92. Los parámetros del WG-S84 son:

Semieje mayor: a = 6 378 137m

Achatamiento: f =1 / 298,257 223 563

Velocidad angular de la Tierra: ω = 7 292 115 *10-11 rad/s

Constante gravitacional: µ = 3 986 004,418 *108 m3/s2

(*) G730 define la fecha de implantación del sistema, corresponde a la semana GPS 730 (junio 1994)

SISTEMA SIRGAS.

La comunidad geodésica de América del Sur ha desarrollado un proyecto, aún en

ejecución, denominado SIRGAS (SIstema de Referencia Geocéntrico para América del Sur), con
 10-04-02 [Link]- pag. 7 de 11

fines de adoptar, para el continente, una red de referencia de precisión compatible con las técnicas
modernas de posicionamiento, principalmente GPS. SIRGAS adopta como sistema de referencia el
ITRF y elipsoide del WGS-84, que en la práctica es idéntico al WGS-84. El órgano representativo de
Chile ante el Proyecto SIRGAS, es el Instituto Geográfico Militar – IGM, responsable por la Red
Geodésica Nacional - RGN.

La materialización de SIRGAS se inició en mayo de 1995 y fueron ocupadas

simultáneamente, durante 10 días, 58 estaciones en todo el continente, siendo 8 en Chile. El resultado
de la solución final definió coordenadas referidas a ITRF95,4 , con precisión de orden científica. Una
segunda campaña de SIRGAS entregará coordenadas relacionadas a otra época de referencia y las
variaciones temporales de las coordenadas geocéntricas para cada estación.

La adopción de SIRGAS por parte de varios países sudamericanos como referencia para
los sistemas geodésicos nacionales, refleja la tendencia global de compatibilizar éstos a las tecnologías
modernas y Chile en el futuro no debe estar fuera de ese contexto.

(tarea: investigar SIRGAS)

SISTEMAS PSAD-56 Y SAD-69.

En las décadas de los años cincuenta y sesenta, para fines geodésicos y cartográficos
fueron definidos los sistemas de referencia sudamericanos Datum Provisorio Sudamericano 1956 –
PSAD-56, con su vértice de origen en La Canoa, Venezuela y Datum Sudamericano 1969 – SAD-69,
con origen en Chua, Brasil.

En Chile, el Instituto Geográfico Militar (IGM) implementó el PSAD-56 como sistema de

referencia oficial para el territorio nacional desde el extremo norte hasta la latitud 43º 30‘ Sur, lo que
coincide aproximadamente con el límite entre las regiones X y XI. En el sur de Chile se usa el SAD-
69 como referencia cartográfica, como también el datum Hito XVIII en el extremo sur de la XII
Región. La cartografía sistemática escala 1/50 000 editada por el IGM es referida a los datums PSAD-
56, SAD-69 e Hito XVIII, en cada región correspondiente. Las cartas escala 1/25 000 son referidas al
SAD-69.
 10-04-02 [Link]- pag. 8 de 11

Datum Elipsoide Semi-eje mayor (a) Achatamiento (1/f)

PSAD-56 Internacional (Hayford) 6 378 388m 297

SAD-69 SAD-69 (UGGI-67) 6 378 160m 298,25

Hito XVIII Internacional (Hayford) 6 378 388m 297

WGS-84 GRS-80 6 378 137m 298,257 223 563

Transformación de Sistemas.

La transformación de coordenadas respecto a diferentes sistemas es de fundamental

relevancia al compatibilizar sistemas de referencia, especialmente en posicionamiento por GPS, donde
esta fase de cálculo se realiza de forma automática por los programas que acompañan a los equipos y
no siempre es clara la metodología ni los valores utilizados, llevando en algunos casos, al usuario a
realizar una transformación poco rigurosa o incorrecta.

La figura ejemplifica la relación entre el datum GPS, WGS-84 y el datum SAD-69. Por
definición, ellos son considerados paralelos, habiendo en este caso solo translación tridimensional
entre sus orígenes, aunque en su materialización puedan existir deformaciones que produzcan
rotaciones. Por ese motivo, el trío de valores correspondientes a tal translación, se denominan
“parámetros de transformación” – PT entre datums, a saber ∆X, ∆Y y ∆Z, los que deben ser
adicionados (considerando su signo) a las coordenadas cartesianas (X, Y, Z) del punto a ser
transformado.

ZSAD-69

ZWGS-84
YSAD-69

XSAD-69
TZ YWGS-84
TX
TY
XWGS-84

Entre los diversos enfoques para la transformación de coordenadas, las formas más usadas
de aplicar los PT, son: las Ecuaciones Diferenciales de Molodensky y los Modelos Cartesianos,
 10-04-02 [Link]- pag. 9 de 11

aunque también existe el modelo de “Regresión Múltiple”, basados en desarrollo polinomial. A

continuación se muestran los dos primeros modelos.

Ecuaciones Diferenciales Simplificadas de Molodensky.

Este modelo posee la particularidad de transformar coordenadas, del primer datum al

segundo, en un solo modelo de ecuaciones, donde las coordenadas transformadas al 2o datum son
dadas por: φ2 = φ1 + ∆φ , λ2 = λ1 + ∆λ y h2 = h1 + ∆h

1
∆φ = [(a1 ⋅ ∆f + f1 ⋅ ∆a) ⋅ sen2φ1 − TX ⋅ senφ1 ⋅ cosλ1 − TY ⋅ senφ1 ⋅ senλ1 + TZ ⋅ cosφ1 ]
M1
1
∆λ = [−TX ⋅ senλ1 + TY ⋅ cosλ1 ]
N1 ⋅ cosφ1
∆h = TX ⋅ cosφ1 ⋅ cosλ1 + TY ⋅ cosφ1 ⋅ senλ1 + TZ ⋅ senφ1 + (a 1 ⋅ ∆f + f1 ⋅ ∆a) ⋅ sen 2 φ1 − ∆a

∆φº = ∆φ ρº ∆λº = ∆λ ρº ρº = 180/π

2
a 1 ⋅ (1 − e1 )
Con: M1 = ; e2 = 2 ⋅ f − f 2
2 2 3
(1 − e1 ⋅ sen φ1 )

a1 , f1 : parámetros del primer elipsoide

a2 , f2 : parámetros del segundo elipsoide

∆a = a2 – a1 , ∆f = f2 – f1

TX, TY, TZ : parámetros de translación entre los datums

Modelo Cartesiano.

La aplicación de las ecuaciones trigonométricas de transformación de coordenadas

geodésicas a cartesianas, y viceversa, comenzaron a emplearse preferentemente con el advenimiento
de las computadoras y aunque siendo este método de mayor número fases analíticas, es de más fácil
visualización en cuanto a su concepto.
 10-04-02 [Link]- pag. 10 de 11

Sist. Geodésico 1 Sist. Geodésico 1

(φ, λ, h)1 (X, Y, Z)1

MOLODENSKII CARTESIANO
Sist. Geodésico 2 Sist. Geodésico 2

(φ, λ, h)2 (X, Y, Z)2

Este se basa en la conversión a coordenadas cartesianas y aplicación por separado de los

PT, de acuerdo a las siguientes etapas:

Convertir las coordenadas geodésicas en el primer sistema a coordenadas cartesianas en el mismo

sistema: (φ,λ,h)1 (X,Y,Z)

Aplicar los PT a las coordenadas cartesianas, trasladando el origen del sistema al segundo
sistema: (X,Y,Z)1 + (TX,TY,TZ) = (X,Y,Z)2

Convertir las coordenadas cartesianas del segundo sistema a coordenadas geodésicas: (X,Y,Z)2
(φ,λ,h)2

Parámetros de Transformación.

Considerando los sistemas geodésicos materializados según la región geográfica de que se

trate, los programas utilizados y la literatura técnica consultada indican diferentes valores para los PT.
Se debe prestar especial atención a los valores escogidos debido a que no existen valores únicos. Es
esencial aplicar los mismos PT a coordenadas dentro de un mismo proyecto, con el fin de evitar
deformaciones sistemáticas. Alternativamente se pueden calcular PT particulares con validez para una
área específica, en base a puntos comunes con coordenadas conocidas en los dos sistemas geodésicos,
o datums, considerados. Se destaca que esta alternativa requiere extremo cuidado, ya que, debido a la
imprecisión en las coordenadas, los valores de PT calculados serán igualmente imprecisos.

Se muestra a título indicativo los valores de PT entre diferentes datums, calculados y

difundidos por la Agencia Nacional Estadounidense de Imágenes y Mapas – NIMA:
 10-04-02 [Link]- pag. 11 de 11

TRANSFORMACIÓN VALORES [m] OBSERVACIÓN

TX = -270 ± 25 Válidos para Chile, al norte del
PSAD-56 WGS-84 TY = +183 ± 25 paralelo 19ºS aproximadamente;

TZ = -390 ± 25 calculados con 1 estación

TX = -305 ± 20 Válidos para Chile, al sur del

PSAD-56 WGS-84 TY = +243 ± 20 paralelo 19ºS aproximadamente;

TZ = -442 ± 20 calculados con 3 estaciones

TX = +16 ± 25 Válidos para Chile, al sur del

Hito XVIII WGS-84 TY = +196 ± 25 paralelo 53ºS aproximadamente;

TZ = +93 ± 25 calculados con 2 estaciones

TX = -75 ± 15 Válidos para todo Chile,

SAD-69 WGS-84 TY = -1 ± 8 calculados con 9 estaciones

TZ = -44 ± 11

Especial atención merecen las precisiones asociadas a los parámetros de translación. Programas que
acompañan los equipos GPS, pueden tener incorporados diversos valores de PT, aplicados según
diferentes modelos de transformación de sistemas, incluyendo una opción para que el usuario imponga
valores de PT y modelo de transformación, según su propio criterio.
(tarea: investigar como se comporta la diferencia de traslación en las coordenadas horizontales,
en metros)
 10-04-02 [Link]- pag. 1 de 15

TRIANGULACIÓN

La materialización del sistema cartográfico se efectúa por un proceso llamado de “triangulación”

(proviene de la época en que se formaban solo triángulos adyacentes) y los puntos llamados “vértices
geodésicos”. La marca física se llama marco o monolito y representan geométricamente las visadas
realizadas con teodolito. Las medidas son hechas sobre el terreno.

Imposiciones iniciales:

Posición; punto origen (o datum) de coordenadas geodésicas conocidas, amarración al elipsoide

(evita translación);

Orientación; conocer un acimut de partida (evita rotación);

Escala; base geodésica inicial medida y reducida, significa imponer una escala.

Con las 4 imposiciones iniciales se puede proyectar la triangulación sobre el elipsoide y se pueden
transportar coordenadas.

El problema del datum

La mayoría de las triangulaciones geodésicas se caracterizan por la imposición inicial:

ξ0 = η0 = N0 = 0

Desvío de la vertical = 0; ξ: componente meridiana; η: componente 1o vertical

Coinciden en el punto origen la normal y la vertical y, elipsoide con geoide coincidencia entre
coordenadas geodésicas y astronómicas.

φ0 = φa ; λ0 = λa ; A0 = Aa

De esta manera fueron definidas 3 de las imposiciones iniciales, por determinaciones astronómicas. La
escala se define “midiendo” una base.

PSAD-56. Con origen en el vértice La Canoa (Venezuela); elipsoide Hayford (Internacional).

Hito XVIII = PSC-63 (Provisional Southern Chile 63). Origen en el Hito XVIII; elipsoide
Internacional.

SAD-69. Origen en vértice Chuá (Brasil); elipsoide SAD-69 (GRS-67 o UGGI-67).

 10-04-02 [Link]- pag. 2 de 15

φ0 = 19º 45’ 41,6527”S ; λ0 = 48º 06’ 04,0639”W ; Aa = 271º 30’ 04,05”

φa = 41,34” S ; λa = 07,80” W ; Aa = 05,42”

ξ0 = 0,31” ; η0 = 3,59” ; N0 = 0

Paralelo al sistema terrestre medio.

1980: UGGI definió el GRS-80 adoptado por Transit y después por GPS.

Cálculo de la triangulación.

Operaciones astronómicas: puntos de Laplace (φ y A)

Operaciones geodésicas:

- Medida de base geodésica inicial y de verificación.

- Medida de ángulos o direcciones horizontales.

- Medida de ángulos verticales.

La triangulación (y cualquier medición) está sujeta a errores observación (accidentales) figuras no

cierran coordenadas dependen del “camino” utilizado. Después de eliminar los errores
sistemáticos, se procede al ajuste por el Método de los Mínimos Cuadrados (MMC) solución única
coordenadas únicas.

Precisión nominal ..... 1/100.000 ó 10ppm.

Los ángulos verticales se usan en la nivelación trigonométrica, no tan precisa como nivelación
geométrica. Esta es conducida en redes separadas a la triangulación.

Rigidez.

La figura básica de la triangulación es el cuadrilátero completo, con lados y diagonales visados en

ambas direcciones.

Rigidez es un concepto precursor de las técnicas modernas de optimización basadas en elipses de

errores. En la práctica no se usa mas. Consiste, a partir de la ley de propagación de errores, en el
error probable del logaritmo de un lado, calculado por la solución de uno o varios triángulos.
 10-04-02 [Link]- pag. 3 de 15

Medidas de Ángulos Horizontales.

Repetición y reiteración.

Repetición: necesita teodolito repetidor: arrastrar el limbo horizontal. Aprovecha para multiplicar de
forma mecánica un ángulo observado.

Reiteración. Se usa diferentes partes del limbo horizontal.

B
L1
L2 α

Método de las direcciones. El observador cala las direcciones 1, 2, 3 y 4, todas al mismo origen; pasa
a transito y vuelve calando a 4, 3, 2 y 1. Este conjunto de mediciones es una serie. En la Segunda
serie reinicia en transito y vuelve en directa.
1

1a serie Ida directa 2a serie Ida transito

Vuelta transito Vuelta directa

Según normas internacionales, un ángulo no se debe desviar más de 4” de la media.

Método de Schreiber. Se inicia como el método e las direcciones, 1a serie; se abandona el 1er punto
calado y se hace otra serie; se abandona el 2o punto calado y se hace otra serie, así hasta el penúltimo
punto. Se usa el método de Mínimos Cuadrados para calcular cada ángulo.
1

Reconocimiento. Objetivo: adoptar una figura prefijada; asegurar intervisibilidad; lugares accesibles y
permanentes; evitar refracción lateral.
 10-04-02 [Link]- pag. 4 de 15

Materialización. Marcos o monolitos de concreto con identificación y marca de centrado. Mediciones

de 1er orden se efectúan normalmente en la noche, por causa de menos refracción y visibilidad. Se
necesitan lámparas eléctricas.

Centrado. En caso de no haber visibilidad se construye una torre de madera o metal. Las
prefabricadas usadas por el DMA se llaman Torres Bilby de hasta 38m; consiste en dos torres
concéntricas independientes. Plomo óptico ni cualquier otro sirve para instalar el teodolito; se usa un
colimador vertical.

MEDICIONES DE BASES GEODÉSICAS.

Distanciómetros mecánicos. Hasta fines del siglo 19 eran medidas con reglas rígidas bimetálicas. En
1885 en Suecia se inventó la liga Invar (acero con 35% níquel) con bajo coeficiente de dilatación.

Coef. dilatación acero ~ 0,000 01

Coef. dilatación invar ~ 0,000 000 04 (250 X); más blando y menos elástico que acero. Se usó
hasta unos 40 años atrás.

Hilos (o cinta) de invar eran de 24 o 50m y poseían en los extremos pequeñas reglas graduadas de 8
centímetros . Se debía estacar cada 24 o 50m ±15cm. Se aplican tensores de 10kg en los extremos. Se
hace nivelación geométrica en las estacas para reducir al horizonte más correcciones de temperatura,
flecha, variación de tensión, etc. Después era reducida al geoide. .... semanas de trabajo para medir
una base.

Distanciómetros electrónicos. En 1947 se inventó en Suecia el Geodímeter (Geodetic Distance Meter)

basado en onda luminosas; en 1957 en Sudáfrica se inventó el Telurómetro (microondas), el más
popular en geodesia.

REDUCCIONES GEOMÉTRICAS.

Existe una reducción meteorológica que depende de temperatura, presión y humedad, medidos durante
la operación en terreno.

Para convertir la distancia electrónica en distancia geodésica se deben efectuar algunas correcciones
geométricas.
 10-04-02 [Link]- pag. 5 de 15

a- Reducción al horizonte. La distancia electrónica (De) es inclinada y la distancia horizontal (Dh)

es:

Dh = De * cosα ; α = ángulo vertical respecto al horizonte, o

Dh2 = De2 - ∆H2 ; ∆H = desnivel

Siendo la corrección Ch = Dh – De = (De2 - ∆H2)1/2 – De

b- Reducción al elipsoide. Designando por H la altura geométrica (elipsóidica) de la base, o del lado
de la poligonal; por Dp el lado proyectado en el elipsoide:
Dh

Dp H

Dh R + H
=
Dp R
Dh ⋅ R
Dp =
R+H

sea la corrección Ce = Dh – Dp

H H2
CP = Dh ⋅ ( − 2 + ...)
R R

R: radio de la sección normal del acimut. Fórmula de Euler. En ciertos casos se puede tomar
un radio medio.

(notar semejanza con la corrección al NMM, deducida por otro camino)

1 sen2 A cos2 A
= +
RA N M

Recordar que H = h + N; cuando N es desconocido (o era) los geodestas reducían las bases al geoide,
no obstante los cálculos de la triangulación eran realizadas en el elipsoide.
 10-04-02 [Link]- pag. 6 de 15

c- Reducción al arco. Esta corrección es aditiva (contrario a las anteriores)

DG
DP

α
DG = R ⋅ α
α α α3
DP = 2R ⋅ sen = 2R ⋅ ( − + ...)
2 2 48

R ⋅ α3 DP 3
CG = DG − DP = =
24 24R 2

(importante: estas son las 3 reducciones significativas a las distancias)

TRIÁNGULO GEODÉSICO.

El objetivo inmediato de la geodesia es determinar coordenadas de puntos en la triangulación,

trilateración y poligonales de precisión. Las coordenadas son transportadas vértice a vértice.
Significa que todos los lados deben ser conocidos y surge el problema de resolver los triángulos
geodésicos.

No se pueden utilizar los valores observados, pues las coordenadas dependerían del camino utilizado.
Debe haber solución “única”, para eso los valores deben ser ajustados. Aquí no se abordará el ajuste
por el Método de los Mínimos Cuadrados (MMC).

El ajuste presupone una resolución preliminar de los triángulos geodésicos.

El modelo geométrico adoptado es el elipsoide, por lo tanto los triángulos geodésicos son elipsóidicos
y no esféricos. Pero como el elipsoide tiene poca excentricidad, se puede asumir como esféricos,
siempre que se les atribuya una esfera correspondiente de radio igual a el radio medio R0 en función
de la latitud del centro de gravedad del triángulo, o latitud media.

R0 = M ⋅ N

En un triángulo esférico, la suma de los ángulos interiores es mayor que 180º (excluyendo errores de
medición)
 10-04-02 [Link]- pag. 7 de 15

EXCESO ESFÉRICO

Es el valor que excede dos ángulos rectos.

ε = A+B+C-180º ó =A+B+C-200g

TEOREMA DE LEGENDRE. Se conservan los lados y varían los ángulos.

A A’

γ β
γ’ β'
S
S’

B’ C’
B
α C α'

Sean dos triángulos, esférico y plano, cuyos lados correspondientes son iguales en longitud:

α = α’ ; β = β’ ; γ = γ’

Condiciones del teorema de Legendre:

1. Los dos triángulos tienen la misma área (S = S´)

2. Los ángulos del triángulo esférico son iguales a los correspondientes del triángulo plano, más 1/3
del exceso esférico (ε) (A-A´=ε/3)

En la práctica:

Se miden triángulos elipsóidicos.

Se suponen esféricos

Se calcula el exceso esférico se “transforman” en planos

Se calculan los lados

Se transportan coordenadas

S'
ε= S´: sup erficie triángulo esférico
R2
S = 21 ⋅ β ⋅ α ⋅ senA´ R = MN

β ⋅ α ⋅ senA '
ε" = ρ"
2MN
 10-04-02 [Link]- pag. 8 de 15

1
factor de exceso esférico : m = ρ" ; ε " = m ⋅ β ⋅ α ⋅ sen A '
2MN

Cálculo provisorio del triángulo.

El cálculo definitivo de los triángulos es después del ajuste de los ángulos y este presupone conocer el
exceso esférico, que a su vez exige el cálculo preliminar de los triángulos.

Recordando que cuando los ángulos de un triángulo tienen el mismo peso, son corregidos en 1/3 del
error de cierre (E). Por Legendre, los ángulos planos corregidos se obtienen de los ángulos esféricos
“observados”.

ε+E A + B + C − 180 o
A' = A − =A− E : error de cierre
3 3
ε+E A + B + C − 180o
B' = B − =B−
3 3
ε+E A + B + C − 180o
C' = C − =C−
3 3

Resuelto el cálculo provisorio (ángulos A´,B´ y C´), se determina el exceso esférico (ε),
posteriormente se ajustan los triángulos y finalmente se transportan las coordenadas.

Reducción de valores observados.

Obviamente los valores brutos no pueden ser considerados en los cálculos geodésicos. Ellos vienen
con errores de medición debido al observador, al equipo y a los efectos ambientales. Después de la
eliminación de los errores sistemáticos, las observaciones son ajustadas.

Como las observaciones son realizadas en la superficie física de la Tierra y los cálculos son en el
elipsoide, ellas deben ser reducidas.

1 Ángulos horizontales.

1.1 reducción geométrica: efecto de la altura de la mira y ángulo sección normal - geodésica

1.2 reducción física: medidas sobre la vertical y cálculos en el elipsoide: corrección del desvío de la
vertical

2 distancias.
 10-04-02 [Link]- pag. 9 de 15

2.1 Reducciones geométricas: al horizonte; al elipsoide; al arco

3 Ángulos verticales. Se verá en nivelación trigonométrica.

MÉTODO DE LOS ADITAMENTOS (O AGREGACIONES). Se mantienen fijos los ángulos y se

modifican los lados

Permite resolver triángulos geodésicos como triángulos planos

sen A α'
Del triángulo plano: =
sen B β'

α3 β3 χ3
α' = α − β' = β − γ' = γ −
6R 2 6R 2 6R 2

aditamentos que convierte el triángulo esférico en plano

Este procedimiento es poco usado, pero tiene la ventaja que de una serie de triángulos solo el primer
lado es convertido a plano; de allí en adelante se resuelven todos como planos.

TRANSPORTE DE COORDENADAS

N
N
A12
A23
2 A21
S12
1
α

S23

4
3

Problema directo de la geodesia.

Dadas las coordenadas de un punto P1 en el elipsoide, la distancia y acimut a un punto P2, determinar
las coordenadas de P2.

(φ1, λ1, S1-2, A1-2) calcular (φ2, λ2, A21)

φ2 = φ1 + ∆φ ; λ2 = λ1 + ∆λ ; A21 = A12 + ∆A + 180º

Este es el objetivo inmediato de la geodesia (transporte de coordenadas).

 10-04-02 [Link]- pag. 10 de 15

Problema inverso de la geodesia.

Dadas las coordenadas de dos puntos P1 y P2 en el elipsoide, determinar la distancia y los azimutes
entre P1 y P2.

Z dλ
N cosφ

N cosφ dλ

M dφ dS

(φ1, λ1, φ2, λ2) calcular (S1-2, A1-2, A2-1)

Se pueden considerar dos hipótesis:

Distancias pequeñas < 50km, fórmulas son equivalentes en precisión

Distancias grandes > 50km, fórmulas pierden precisión respecto a la distancia y deben ser tratados
con más cuidado.

M ⋅ dφ = cos A ⋅ ds
N ⋅ cos φ ⋅ dλ = sen A ⋅ ds
dA tgφ
= sen A
ds N

desarrollando en serie:

dA − sen φ dλ dλ dλ
γ= = − sen φ
ds ds ds
 10-04-02 [Link]- pag. 11 de 15

FÓRMULAS DE PUISSANT (PROBLEMA DIRECTO).

P’

P
’
-φ 1
90
e

2’
oid
ra

φ
fe

ps

90-
i
es

el

P2

A12 s φ2

P1
φ1

Ecuador elipsoide φ

Ecuador esfera φ H

Las fórmulas derivadas por el geodesta francés Puissant son las más utilizadas en líneas hasta 80km
[Rapp].

Latitud

Se utiliza la esfera de Jacobi en el teorema de Dalby para el transporte de latitud. La esfera auxiliar es
tangente al elipsoide en el 1er vértice.

P1 es el vértice conocido y B2 es el vértice a determinar.

PP1 y PP2 son los meridianos geodésicos; P’P1 y P’P2 son los respectivos esféricos.

Elipsoide y esfera tienen en común el paralelo (φ) de P1, su radio coincide con la Gran Normal (N’) en
ese punto, su centro es en (o); significa que en el vértice P1 las latitudes esférica y geodésica son
iguales. El acimut y las distancia son iguales en la esfera y en el elipsoide.

Fórmula de los 4 elementos aplicado al lado P’P1 del triángulo esférico y el valor angular (σ) del arco
P1P2 , o sea σ = s/N

s s
sen φ2 = sen φ1 ⋅ cos + cos φ1 ⋅ sen ⋅ cos A12
N1 N1

s s
sen(φ1 + ∆φ' ) = sen φ1 ⋅ cos + cos φ1 ⋅ sen ⋅ cos A12
N1 N1
 10-04-02 [Link]- pag. 12 de 15

2 3 3
s s s ∆φ'
∆φ' = ⋅ cos A12 − 2
⋅ tgφ1 ⋅ sen2 A12 − 3
⋅ cos A12 +
N1 2N1 6N1 6

Una mejor aproximación en la iteración es:

s s2
∆φ' = ⋅ cos A12 − 2
⋅ tgφ1 ⋅ sen2 A12
N1 2N1

s s2 2 s3
∆φ' = ⋅ cos A12 − 2
⋅ tg φ1 ⋅ sen A 12 − 3
⋅ cos A12 ⋅ sen2 A12 ⋅ (1 + 3tg2φ1 ) + ...
N1 2N1 6N1

Hasta ahora se trabajó en una esfera de radio N y debe cambiar ∆φ’ (medido en la esfera de radio N1)
por ∆φ a lo largo del meridiano. Para eso se asume que N1 ∆φ’ en la esfera es igual a la distancia
correspondiente en el elipsoide. Permitiendo que Mm sea el radio meridiano de curvatura en la latitud
media.

N
N ⋅ ∆φ' = Mn ⋅ ∆φ ⇒ ∆φ = ∆φ'⋅ 00
Mm

δφ: diferencia aproximada en latitud entre P1 y P2

∆φ = B ⋅ s ⋅ cos A12 − C ⋅ s2 ⋅ sen2 A12 − E ⋅ h ⋅ s2 ⋅ sen2 A12 − D ⋅ (δφ)2

Con:

s s2 s3
δφ = ⋅ cos A12 − ⋅ sen2 A12 ⋅ tgφ1 − 2
⋅ sen2 A12 ⋅ cos A12 ⋅ (1 + 3tg2φ1 )
M1 2M1N1 6M1N1
1 tgφ1 3c 2 ⋅ sen φ1 ⋅ cos φ1
B= C= D=
M1 2M1N1 2(1 − e 2 ⋅ sen2 φ1 )
1 + 3tg2φ1 s ⋅ cos A12
E= 2
h=
6N1 M1

φ2 = φ1 + ∆φ
 10-04-02 [Link]- pag. 13 de 15

Esta serie es rápidamente convergente. Los términos C y D son negativos en el hemisferio sur.
Antiguamente (sin computadores) los cálculos eran realizados con tablas de logaritmos de B, C, D, y
E.

Longitud.

Analogía de los senos al triángulo elipsóidico, supuesto esférico P1P2P’ (de radio N) con φ2 conocido.

sen ∆λ sen A12 s sen A12 s

= ⇒ sen ∆λ = sen ⋅ ⇒ sen ∆λ = sen ⋅ sen A12 ⋅ sec φ2
s cos φ2 N2 cos φ2 N2
sen
N2

Desarrollado en serie:

s s2
∆λ = ⋅ sen A12 ⋅ sec φ2 [1 − 2
⋅ (1 − sen2 A12 ⋅ sec 2 φ2 )]
N2 6N2

Contra acimut

Usando las analogías de Napier:

c b

B C

cos 21 (b − c )
tg 21 (B + C) = ⋅ cot 21 A
cos 21 (b + c)

B = A12 ; b = 90º-φ2

C = 360º-A21 ; c = 90º-φ1

A = ∆λ

A21 = A12 + ∆A + 180º

∆A : convergencia meridiana
 10-04-02 [Link]- pag. 14 de 15

3
∆λ
∆A = ∆λ ⋅ sen φm ⋅ sec 21 ∆φ + ⋅ (sen φm ⋅ sec 21 ∆φ − sen3 φm ⋅ sec 3 21 ∆φ
12

Resumen

Calcular M y N

Calcular ∆φ’ aproximado

Calcular ∆φ y φ2

Calcular N2

Calcular ∆λ y λ2

Calcular ∆Α y A21

PROBLEMA INVERSO.

(φ1, λ1, φ2, λ2) calcular (S1-2, A1-2, A2-1)

usando la ecuación directa de la longitud de Puissant.

N2 ⋅ ∆λ ⋅ cos φ2
s ⋅ sen A12 = 2
=X (1)
s 2 2
[1 − 2
⋅ (1 − sen A12 ⋅ sec φ2 )]
6N2

1 (2)
s ⋅ cos A12 = ⋅ [ ∆φ + C ⋅ s 2 ⋅ sen2 A12 − E ⋅ h ⋅ s2 sen2 A12 + D ⋅ (δφ)2 ]
B
1
s ⋅ cos A12 = ⋅ [ ∆φ + C ⋅ X 2 − E ⋅ h ⋅ X 2 + D ⋅ (δφ)2 ]
B

s ⋅ sen A12 (3)

tgA12 =
s ⋅ cos A12

Estas fórmulas son iterativas, acimut (A12) y distancia (s) están presentes en las ecuaciones.

Calcular N1 , N2 y M1

Usando solo el denominador de (1) calcular (s*senA12), reemplazar en (2) y calcular (s*cosA12)

Calcular (tgA12) con (3), resulta A12 aproximado

Usando (1) o (2) calcular “s”

Realizar iteraciones.
 10-04-02 [Link]- pag. 15 de 15

Fórmulas no iterativas usando la latitud media de Gauss (hasta ~40km).

∆λ
cos φm ⋅ sen
∆A 2
(4)
tg( A12 + )=
2 M ∆φ
sen[ m ⋅ ∆φ] ⋅ cos
2Nm 2

∆λ
cos φm ⋅ sen
s 2 (5)
sen =
2Nm sen( A + ∆λ )
12
2

El lado derecho de (4) es conocido y ∆A es conocido

Con la (4) se encuentra el término (A12+∆A/2) y se calcula el acimut A12

Con la (5) se calcula “s”

 10-04-02 [Link]- pag. 1 de 5

ALTIMETRÍA.

SUPERFICIE
H TERRESTRE

GEOIDE
N h

ELIPSOIDE
h=H+N

La superficie de referencia altimétrica es el geoide, el cual está definido como la superficie

equipotencial (de igual valor de atracción gravitacional) que coincide con la superficie de los océanos
en reposo, extendida sobre los continentes, su denominación más común es Nivel Medio del Mar –
NMM. La altura sobre el geoide (o sobre el NMM) se denomina “altura ortométrica”, también
referida como altitud o elevación.

La “altura ortométrica” se definie como la distancia vertical desde el geoide a un punto en

la superficie de la Tierra. La altura elipsóidica se mide por la normal al punto en la superficie
terrestre, como muestra la figura. Para fines prácticos ellas se consideran colineales, aunque
rigurosamente no lo son. La relación entre la superficie elipsoidal y la superficie del geoide está dada
por la “ondulación geoidal” , designada por “N”, ella representa en un punto la altura del geoide
respecto al elipsoide. El conocimiento de este valor es necesario para la reducción de alturas
elipsóidicas a alturas sobre el NMM, de acuerdo a la expresión: h = H + N. La altura elipsóidica sólo
interesa en posicionamiento con GPS.

El tratamiento matemático del geoide es un problema complejo que se resuelve

puntualmente y el usuario debe recurrir a modelos geoidales. Los modelos existentes están presentes
en algunos programas computacionales de procesamiento GPS o se puede recurrir externamente a
modelos continentales o aún, modelos globales modernos como el EGM96 (Earth Gravity Model
1996), para obtener valores de ondulación. El EGM96 es de uso público y está a disposición un
programa de extracción automática, con su respectivo banco de datos, en el “web site” de la NASA
([Link] o [Link]
84/[Link]). Aún así, los usuarios deben estar atentos a nuevos modelos globales, continentales o
regionales mejor adaptados a nuestro país.

(tarea: investigar EGM96)

 10-04-02 [Link]- pag. 2 de 5

SISTEMA TOPOCÉNTRICO O SISTEMA LOCAL DE COORDENADAS.

y
N
r
x
o P1
Z
h

o φ Y
λ

La situación correspondiente a un Plano Topográfico Local – PTL puede ser ilustrada por
la figura, donde el área que se desea representar es una porción de la superficie terrestre. En este caso
la figura matemática de la Tierra pasa a estar representada por una esfera de radio igual al radio medio
del elipsoide. El PTL es un sistema de representación plano-rectangular, horizontal topocéntrico para
fines topográficos, con un origen de coordenadas geodésicas conocidas.

B
Terreno
SH
∆H
CM A
Ω
CA Plano
SP ∆S
Topográfico
Local
HPTL
NMM

Para que el PTL pueda ser considerado “georreferenciado”, este debe estar ligado
geométricamente al sistema elipsóidico, en posición, orientación y altura del plano horizontal.
Generalmente esta ligazón está dada por:
 10-04-02 [Link]- pag. 3 de 5

Posición: dada por las coordenadas geodésicas (φ0,λ0) del origen del sistema rectangular local
(X0,Y0);

Orientación: dada por el acimut geodésico (α0) del eje “Y”, ordenada del sistema local;

Altura del plano de referencia: dada por la altura (h0) del plano de referencia.

Aún cuando el área representada en el PTL sea suficientemente pequeña como para no
considerar la curvatura terrestre, el efecto de la altura del PTL y de los desniveles requiere un factor de
corrección para las distancias horizontales. En efecto, según la figura, la distancia horizontal media
entre los puntos A y B del terreno (SH) se proyecta en el plano del horizonte del PTL como SP; este
plano está a una altura HPTL. La magnitud del factor de reducción de distancias horizontales (∆S) a ser
representadas en el PTL, depende de la altura media al PTL (CM), de la distancia horizontal entre los
puntos y en menor medida del valor de HPTL. Esto implica que la adopción de un plano medio
involucra criterios que relacionan el área del proyecto (distancias máximas) y desnivel del terreno
(desniveles máximos).

La distancia reducida al plano topográfico (SP) se calcula por la siguiente expresión:

SP = SH - ∆S

El factor de reducción ∆S puede calcularse a partir de:

R + HPTL + CM CM
=
SH ∆S
SH ⋅ CM
∆S =
R + HPTL + CM

Sean:

R: radio medio, puede usarse un valor aproximado, por ejemplo,

R = 6 378 000m

SH: distancia horizontal media entre A y B

Sp: distancia reducida al PTL

∆S: corrección a la distancia horizontal

HPTL: altura media del PTL respecto al NMM

CM: cota media (ente A y B) respecto al PTL

 10-04-02 [Link]- pag. 4 de 5

De la proporcionalidad de triángulos:

Considerando que (HPTL + CM) es de magnitud pequeña frente a R, la expresión se escribe:

SH ⋅ CM
∆S =
R

Donde:

CM: cota media;

SH: distancia horizontal proyectada sobre el plano medio;

R: radio medio de curvatura, puede ser adoptado un radio medio sin perjuicio de la precisión;

HPTL: altura del PTL.

Esta corrección puede ser expresada por un factor proporcional como factor de escala, en mm/km o
ppm.

R+H
KH =
R

Para H>0 KH>0

distancia elipsóidica x KH = distancia terreno

Para ilustrar la magnitud de ∆S, la figura 1 contiene un gráfico que muestra para diferentes SH, la
relación entre el CM y la reducción ∆S. Por ejemplo, para una distancia de 20km y un desnivel de
500m al PTL, la reducción de SH a SP será de 1,568m, que equivale a una corrección de 78,5mm/km
(1/12.740). Por otra parte el gráfico en la figura 2 , muestra la corrección proporcional (en mm/km)
respecto al desnivel. Por ejemplo, a una altura de 300m se debe aplicar a las SH una corrección de
47mm/km, para ser consideradas como distancias en el PTL.

En la práctica, para proyectos viales en áreas con gran desnivel, existen dos medios de
compatibilizar distancias proyectadas en el PTL con distancias horizontales de terreno:

Utilizar solo un PTL. En este caso en las fases de levantamiento y replanteo se deben reducir las
distancias horizontales al PTL. Esto implica que para el cálculo del levantamiento y del replanteo
habrá que considerarse la reducción de acuerdo a intervalos de desniveles, determinados por la
fórmula (3) citada o por el gráfico de la figura 2.
 10-04-02 [Link]- pag. 5 de 5

Utilizar más de un PTL. En esta opción no se aplica la reducción facilitándose las operaciones de
levantamiento y replanteo, pero por otro lado, se deben definir tantos PLT como el desnivel del
área de trabajo lo requiera, de forma que las distancias proyectadas en el PTL no difieran de las
distancias horizontales de terreno, más que la tolerancia horizontal permitida, a fin que en la fase de
replanteo no sea necesario la imposición de un factor de escala a las distancias.

10,0
SH = 50km SH = 40km
9,0
SH = 30km
SH: Distancia horizontal
8,0
Corrección dSH [m]

7,0
SH = 20km
6,0

5,0

4,0
SH = 10km
3,0

2,0 SH = 5km
1,0 SH = 2km
SH = 1km
0,0
0 250 500 750 1000 1250 1500 1750 2000
Altura Plano Medio - Hm [m]

500
450
400
350
Desnivel [m]

300
250
200
150
100
50
0
0 10 20 30 40 50 60 70 80
Corrección [mm/km]