Espiago 2017 Cartografía para La Geografía - Georreferenciación

CARTOGRAFÍA PARA LA GEOGRAFÍA
FUNDAMENTOS
Javier Espiago, Servicio de Cartografía de la Universidad Autónoma de Madrid.
“El cálculo del posicionamiento con exactitud repetible es el problema central para la
referencia geográfica de la información terrestre y la función principal de la geodesia”.
Organización Hidrográfica Internacional, (2006), Manual de Hidrografía, pág. 37.
- JULIO 2017 -
1. GEORREFERENCIACIÓN Y SISTEMAS DE COORDENADAS.
Sistema Internacional de Referencia Terrestre (ITRS - ITRF).

A) Definición del sistema (ITRS).
B) Determinación del sistema (ITRF).
C) Control del sistema ITRF.
Elección de una figura de ajuste.
Coordenadas geográficas.
Cálculo de distancias y direcciones en la tierra esférica.
Elección del radio de la tierra esférica.
Coordenadas geodésicas.
Geometría de la elipse.
Referencias geográficas sobre el elipsoide
Cálculo de distancias y direcciones en la tierra elipsoidal.
Evolución de las determinaciones de las superficies de referencia
Elipsoide equipotencial GRS80, Determinación de elipsoides de ajuste.
Sistemas ETRS89 y WGS84.
2. GEORREFERENCIACIÓN Y REGISTRO DE LA ALTIMETRÍA.
Figuras de referencia.
Geoide.
Elipsoide equipotencial medio terrestre o elipsoide de nivel.
Altitud geográfica y datum vertical.

Determinación satelital del campo gravitatorio.
Modelos de gravedad en la superficie terrestre, EGM2008.
Anomalías de la gravedad.
Tipología de altitudes. anomalías de la altitud.

Otras figuras de la tierra (teluroide. cuasigeoide).
Modificaciones regionales del EGM2008.

Datum vertical europeo (EVRS, European Vertical Reference System).
EGM08- REDNAP
JULIO 2017
CARTOGRAFÍA PARA LA GEOGRAFÍA. FUNDAMENTOS Profesor Javier ESPIAGO
1. GEORREFERENCIACIÓN Y SISTEMAS DE COORDENADAS.

Introducción. Sistemas de coordenadas utilizadas en la georreferenciación.
GEORREFERENCIACIÓN: atribución de coordenadas a los puntos de la tierra. Es anterior a cualquier
sistema de organización, análisis y representación de información geográfica.
Existen varios Sistemas Terrestres de Referencia por Coordenadas (CTRS) aunque se tiende a la
normalización internacional. En este caso, la determinación por coordenadas es inequívoca. De todas for-
mas, perviven numerosos sistemas locales que dan lugar a costosos y complejos problemas de conversión
de coordenadas resueltos por los geodestas. En lugar de la necesaria exigencia de una norma común, se
habla de interoperabilidad o capacidad para el intercambio de datos en diferentes sistemas que han de
describirse, nombrarse y codificarse. En la posición, la interoperabilidad no asegura la sencillez.
_______________
La atribución de coordenadas a los puntos de la tierra es un proceso complejo que podemos describir en las
siguientes fases:
1) definición de un sistema de coordenadas tridimensionales terrestres X, Y Z.
Coordenadas cartesianas según ISO 19111
2) su materialización (determinación práctica) para un momento temporal y
3) el control de los cambios que afectan al sistema adoptado, pues la tierra no es estática ni indefor-
mable.
Con estas tres fases se persigue el posicionamiento astronómico, la localización y la orientación de la tierra
en el Universo. También son aplicables al posicionamiento de cualquier cuerpo celeste, incluyendo la loca-
lización de los satélites artificiales utilizados en geografía. Es decir, las coordenadas cartesianas tridimen-
sionales se emplean en la determinación de puntos de control para la georreferenciación, en navegación
aérea y en la navegación espacial, como es el caso del seguimiento de satélites. Se emplean en navega-
ción de vehículos, que no trataremos aquí, y en geodesia dinámica, donde se necesita una alta precisión.
Es el punto de partida para la posterior observación de las localizaciones de entidades cartográficas.
Para que la georreferenciación tridimensional quede determinada con respecto a la superficie del cuerpo
sólido terrestre es necesaria, además:
4) la definición de un sistema de coordenadas geográficas sobre una figura geométrica que ajuste la
forma y el tamaño de la tierra.
Es necesario un datum o superficie de referencia sobre las que se situarán las nuevas coordenadas (,).
Sin ella, la georreferenciación sería imprecisa y ambigua. La esfera fue utilizada primero como superficie de
ajuste. Actualmente se usa un elipsoide. Son figuras geométricas regulares con las que se ajustan la forma
y la figura física de la tierra en una reduccción matemática que facilita los cálculos.
Georreferenciación y sistemas de coordenadas. 3

En geografía, importa conocer si las coordenadas de un punto lo sitúan en la superficie terrestre, o por
encima o debajo de ésta. Sería muy difícil con el sistema X, Y, Z. Además, el cálculo de la distancia entre p
y d no será la línea recta que une estos dos puntos sino el arco de una elipse. También es conveniente en
geografía la adopción de un sistema que facilite el control de las direcciones convencionales N, S, E y O.
Con las figuras de ajuste adoptaremos un sistema de coordenadas angulares (,) que facilitan los cálculos
de distancias, de direcciones-orientación y de áreas sobre la superficie de la tierra. Son las denominadas
coordenadas geográficas, longitud y latitud. Con las coordenadas cartesianas (X, Y, Z) y con las coor-
denadas geográficas angulares (,), todo punto puede referirse con dos sistemas diferentes. Uno lo sitúa
1
en el espacio el otro, respecto a la tierra .
y
p
a a
b
c
d
x
Las últimas fases del proceso son:

5) determinación de la posición de la figura de ajuste en el sistema de coordenadas tridimensionales X, Y,
Z. En la práctica, puede ser anterior a la obtención de coordenadas geográficas. Actualmente, hacemos
coincidir el centro X, Y, Z con el centro del elipsoide.
Meridiano cero Meridiano local
Ecuador
Coordenadas cartesianas y elipsoidales. ISO 19111
6) determinación de un datum geodésico para el cálculo de la altitud, llamado datum vertical. Puede ser
distinto al datum utilizado para las coordenadas (,) llamado, desacertadamente, datum horizontal.
7) adopción de un sistema de proyección o de transformación de las coordenadas tridimensionales en
coordenadas sobre un plano (lo estudiaremos en un apartado propio).
1
Para el Ordnance Survey (2002): Un sistema terrestre de coordenadas es un sistema de coordenadas planeado
(designed) para la descripción de la posición de objetos sobre la superficie de la tierra (on the land surface of the
Earth)”, pág. 3. También habla de TRS (Terrestrial Reference System).
Las posiciones de los puntos anclados sobre la superficie sólida de la tierra, tienen coordenadas que solo
sufren pequeñas variaciones en el tiempo debidas a efectos geofísicos (deformaciones tectónicas o de
mareas).
Sistema Internacional de Referencia Terrestre (ITRS - ITRF).

A) Definición teórica del sistema (ITRS).
La posición astronómica de la tierra y la de un punto cualquiera asociado a ella (situado en el interior de
esta, en la superficie o en la atmósfera) se pueden proporcionar en un sistema tridimensional con
coordenadas cartesianas o rectangulares X, Y, Z. En la georreferenciación normalizada internacionalmente
2
se utiliza el conocido como International Terrestrial Reference System (ITRS) . Es un sistema ideal o
convencional con el que adoptamos unas propiedades teóricas, matemáticas y físicas.
En su definición se acepta que:
1. la tierra tiene una rotación constante alrededor de un polo astronómico.
2. el eje Z del sistema de coordenadas es el eje de rotación de la tierra determinado en la dirección
señalada por el polo astronómico.
3. el origen es el centro de masas de la tierra incluyendo los océanos y la atmósfera. Por ello el sistema se
denomina geocéntrico.
4. podemos situar un plano que contenga el origen, el eje Z y un punto de la superficie de la tierra
adoptado convencionalmente. El eje perpendicular al eje Z, que pasa por el origen y está contenido en
ese plano define el eje X.
5. el eje Y se establece 90º al este del eje X.
La escala del sistema se basa en la velocidad de la luz. La unidad de longitud es el metro y la de tiempo, el
segundo.
Las coordenadas en este sistema serán coordenadas tridimensionales X, Y, Z geocéntricas.
GTRS. CTRS. El sistema geocéntrico conocido con las siglas GTRS (Geocentric Terrestrial Reference
System) es la manera actual para designar un CTRS (Conventional Terrestrial Reference System). Con el
término CTRS designamos ahora cualquier TRS (Terrestrial Reference System), aunque no esté convenido
internacionalmente (ITRS). Siempre habrá que especificar el origen, la escala y la orientación del sistema.
B) Determinación práctica del sistema (ITRF).

La concreción práctica (la realización o materialización) de este sistema teórico recibe el nombre de
3
International Terrestrial Reference Frame (en siglas ITRF). Ha conocido distintas versiones según el cálculo
de los valores que dan solución al ITRS. Con ellas, se intenta la obtención de una mayor precisión. También
se persigue la adaptación a las variaciones de la masa terrestre y a los cambios de posición del polo
terrestre. Las propiedades que se determinan son el origen, la escala, la orientación y la evolución temporal
del sistema.
1. La posición del polo astronómico es conocida como Conventional Terrestrial Pole (CTP) cuya
4
determinación práctica correspondía al Bureau International de l´Heure (BIH) . Como la rotación terrestre
está afectada por cambios en el eje, fenómeno conocido como migración de los polos, se exige situarlos
siempre con referencia a un momento temporal. Actualmente se usa el primer momento temporal de 1984
(1984,0).
2. Las coordenadas del polo terrestre se implantan por el International Earth Rotation and Reference
5
Systems Service (IERS) con precisión de 0.005“Osobre el CTP. Este polo es conocido como IERS
6
Reference Pole (IRP) siendo el utilizado en la georreferenciación .
2
Como sistema de coordenadas normalizado para todas las ciencias de la tierra de la International Union of Geodesy
and Geophysics (IUGG) y su International Association of Geodesy (IAG) (Asamblea General de Perugia, 2007,
Resolución 2). Se definen conceptos generales en la página del ITRF: http://itrf.ign.fr/general.php.
3
Traducidas a veces en nuestro país como “marcos ITRF” pese a la redundancia. En http://itrf.ensg.ign.fr/ITRF_
solutions/ se descargan las versiones del International Terrestrial Reference Frame ITRF92, ITRF93, ITRF94, ITRF96,
ITRF97, ITRF2000, ITRF2005, ITRF2008 y el ITRF2014. Esta última versión concreta elOorigen, la escala y la orien-
tación. Se describe esta realización en IERS Technical Notes, (nº 38): Z. Altamimi, P. Rebischung, L. Métivier, X.
Collilieux: Analysis and results of ITRF2014, 2017.
4
El sistema BIH se traspasó al IERS orientando el International Terrestrial Reference Frame (ITRF). No es correcta la
utilización0del acrónimo CIO (Conventional International Origin) para nombrar el ITRF pole.
5
Anteriormente International Earth Rotation Service que coincidía en sus siglas (IERS). Es un organismo de la
International Astronomical Union y de la International Union of Geodesy and Geophysics, sustituyendo al International
Polar Motion Service (IPMS) y a la Earth Rotation Section del anteriormente nombrado Bureau International de l'Heure
(BIH).
El movimiento de los polos terrestres presenta un componente anual de trayectoria elíptica y un

componente con un periodo cercano a 435 días de trayectoria circular (componente Chandler). Resulta un
movimiento en espiral que, en su proyección sobre un plano y para distintos periodos de tiempo, vemos en
la siguiente figura del IERS:
Movimiento del Polo 2001-2006 y desplazamiento lineal 1900-2006
https://www.iers.org/IERS/EN/Science/EarthRotation/PolarMotionPlot.html
La migración de los polos se interpreta por los desplazamientos y otros comportamientos de los compo-
nentes físicos de la tierra (materiales del núcleo, manto, corteza, agua subterránea, agua de los océanos y
componentes de la atmósfera). Por ejemplo, la época de deshielo primaveral contribuye al componente
anual lo que ocurre, así mismo, con los grandes desplazamientos de las masas de aire. Aunque no se ha
demostrado con datos numéricos, se piensa que también intervienen los grandes terremotos y los embalses
de agua. Todo ello bajo la acción de la fuerza de la gravedad debida al Sol, los planetas y la Luna y bajo la
acción de la propia rotación. Estas acciones, además de a la migración de los polos, dan lugar a cambios en
la forma de la tierra.
3. El centro de masas es el determinado en el Geocentric Reference System (GRS) de la International
Astronomical Union. Desde 1996 se utiliza un valor que a la masa de la superficie sólida añade la masa de
la atmósfera. También en este caso la posición varía temporalmente debido a la redistribución de la masa
terrestre. Esta variación es conocida como migración del geocentro.
4. Como punto de la superficie instituido convencionalmente para la definición del plano meridiano origen (0°
7 8
0' 0") se elige el Observatorio de Greenwich .
6
http://www.iers.org/
7
En la vertical del Observatorio, el ángulo formado con el plano ecuatorial que sirve de orientación al sistema es de 51°
28' 38". El IERS Reference Meridian (IRM) coincide con el establecido por el BIH como Zero Meridian (epoch 1984,0)
con una incertidumbre de 0,005 “.
8
Es el IERS Reference Meridian. Además de Greenwich (0° 0' 0"), se utilizaron anteriormente otros meridianos de
referencia. Son, según OpenGis: Berna, 7° 26' 22.5" E; Bogotá, 74° 4' 51.3" W; Bruselas, 4° 22' 4.71" E; Ferro (Hierro),
C) Control del sistema ITRF.

Se han elegido estaciones terrestres en las que se observan permanentemente sus coordenadas para
controlar las variaciones del origen y los ejes del sistema con lo que también se controla la orientación. Esa
revisión corresponde al IERS y a un organismo norteamericano (National Imagery and Mapping Agency,
NIMA). Para algunas regiones, como es el caso de Europa, otras instituciones colaboran en el control o la
9.
mejora del sistema
Los movimientos de las placas de la litosfera, que han llegado a ser del orden de 7 cm/año, son las
variaciones de mayor incidencia. Se han elegido 1.499 estaciones que pueden localizarse sobre placas
10
rígidas (high-quality sites) o sobre áreas de deformación . También se analizan los movimientos no lineales
como las variaciones estacionales (anuales o no), las deformaciones debidas a seísmos (124 estaciones).
La atracción de los cuerpos celestes da lugar a mareas oceánicas y a las denominadas mareas terrestres
pues la parte sólida también sufre variaciones que llegan a algunos decímetros en el área emergida. Los
desplazamientos conocidos como mareas terrestres pueden ser temporales o permanentes, aunque
únicamente han de tenerse en cuenta cuando queramos alcanzar precisiones centimétricas. La mayor parte
de los cálculos utilizados en Geodesia no los consideran. Se pueden resolver parcialmente estas
variaciones definiendo la forma promedio lo que está justificado dado que algunas de estas variaciones lo
son en ciclos periódicos (alrededor de dieciocho años para las mareas lunares). La subsidencia, las
variaciones de la orientación, de la longitud de los días y de las fuerzas gravitatorias se unen al estudio de la
Fig. 1.4
estructura profunda de la tierra y a la aplicación de modelos atmosféricos.
17° 40' 0" W; Yacarta 106° 48' 27.79" E; Lisboa, 9° 7' 54.862" W; Madrid, 3° 41' 16.58" W; Paris, 2° 20' 14.025"E; Roma,
12° 27' 8.4" E y Estocolmo, 18° 3' 29" E.
9
“Scientific applications critically depend on the ITRF and impose the most stringent requirements on terrestrial
reference frame accuracy and long-term stability. (…), the applications demanding the highest accuracy and long-term
stability were sea level, geodynamics from vertical land motion, and large-scale horizontal deformation, and decadal
satellite survey missions”, en BLEWIT, Geofrey (2015).
10
Actualmente, la versión es el ITRF2014 “an improved realization of the International Terrestrial Reference System
(ITRS) and is demonstrated to be of higher qualitythan the past ITRF versions”, Altamimi Z. y otros (2017) Analysis and
results of ITRF2014, IERS TN 38.pdf.
El control se lleva a cabo, actualmente, mediante cuatro técnicas de geodesia espacial. Son las
llamativamente denominadas infraestructuras. Primero, el Global Navigation Satellite System (GNSS), que
incluye el norteamericano Global Positioning System (GPS), el ruso GLObal Navigation Satellite System
11
(GLONASS), y, en un futuro, el sistema europeo Galileo . Segundo, los Satellite Laser Ranging (SLR) y el
Lunar Laser Ranging (LLR). En tercer lugar, la Very Long Baseline Interferometry (VLBI). Finalmente, la
12
Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) . Actualmente, tienen,
respectivamente, 19, 32, 39 y 22 años de observaciones. Con el uso de estas técnicas de Geodesia Espa-
cial se calcula la posición y la orientación de las estaciones de control. Después, el origen y la escala del
sistema.
11
http://cddis.gsfc.nasa.gov/Techniques/GNSS/GNSS_Overview.html.
12
Para una posible ampliación: SCHUH, H.; BEHREND, D. (2012), BEHREND, D. (2013).
Actualmente, la Geodesia Dinámica es, por estas técnicas (VLBI, SLR, GPS, DORIS) una Geodesia espacial que
consigue una precisión centimétrica e, incluso, milimétrica. ALTAMIMI, Z, COLLILIEUX, X and MÉTIVIER, L (2013).
Sobre la geodesia espacial: “With the launch of Sputnik 1 on October 4, 1957 (and shortly after of Sputnik 2) modern
space age began. Already these two satellites had a fundamental effect on geodesy. Almost instantaneously a large part
of 100 years of diligent geodetic work dedicated to the determination of the figure of the Earth became out-dated”.
RUMMEL, Reiner (2010).
Para la definición0de0los Marcos de Referencia Terrestre (ITRF), se combinan los resultados de diferentes
técnicas de observación (VLBI, SLR, GNSS, DORIS, gravimetría absoluta, etc.). Cuando en un mismo
observatorio coexisten más de una de estas técnicas de observación (colocación), es necesario conocer
dentro del Sistema de Referencia la posición relativa de las diferentes marcas geodésicas a las que se
refieren las observaciones. Esto es lo que se conoce como la determinación del enlace local (localtie).
Se determinan las coordenadas, el campo de gravedad terrestre y su rotación. Se consigue una elevada
precisión que se aplica al vulcanismo, al estudio de los terremotos, de las regiones tectónicas activas o a los
deslizamientos de tierras. Con el Low Earth Orbit (LEO) se estudian con detalle el nivel del mar, el glaciaris-
mo, los embalses de agua, el agua atmosférica, el movimiento terrestre y las variaciones del campo de gra-
vedad.
Si a las técnicas espaciales añadimos las antiguas técnicas “terrestres” como la nivelación, la fotograme-
tría, la gravimetría o las estaciones que miden el nivel del mar, se estudia entonces el sistema en sus impli-
caciones astronómicas, atmosféricas y oceánicas, de la corteza… Es, de momento, la medida geométrica
de la posición, de las deformaciones y de los cambios del nivel del mar. Después será el análisis físico y la
altitud, como veremos. Cuestiones geométricas como la migración de los polos y el estudio de la rotación
terrestre y la gravedad son los tres ámbitos de análisis y algunos de sus ejemplos que se muestran en el
13
siguiente cuadro del GGOS (Global Geodetic Observing Sistem) de la IAG (International Association of
Geodesy):
13
http://www.ggos.org/. Ver bibliografía y documentación sobre los enlaces de GGOS.

Elección de una figura de ajuste.

SUPERFICIE DE REFERENCIA TERRESTRE: superficie de ajuste a la forma y al tamaño de la tierra en la
que se basa un sistema de coordenadas geográficas.
ESFEROIDE: figura matemática próxima a la figura de la tierra en forma y tamaño. Generalmente, en la
georreferenciación se usa un elipsoide, aunque en algunas deducciones teóricas y en operaciones de
cálculo que no necesiten precisión se usa también la esfera.
Cuando restringimos el cálculo de la posición únicamente a los puntos de la superficie de la tierra, puede
utilizarse un sistema de coordenadas que, por ello, se denomina sistema de coordenadas geográficas,
latitud y longitud (,). La posición relativa de los puntos utiliza estas medidas angulares (,) sobre una
esfera o un elipsoide. A diferencia de las coordenadas X, Y, Z que determinan la posición en el espacio
tridimensional euclidiano, estos sistemas lo determinan sobre la superficie de una figura próxima a la del
cuerpo sólido terrestre. Es también una determinación tridimensional pero restringida a la posición de los
puntos en la superficie de los esferoides. Por ello, las dos medidas angulares (,) han de acompañarse
con el radio de la esfera de ajuste o con los parámetros que definen la forma y tamaño del elipsoide (eje
semi mayor, a y aplanamiento f).
Más adelante puede comprobarse la diversidad de figuras de la tierra adoptadas en la elaboración de
cartografía. Utilizando cada una de ellas, las coordenadas geográficas serán distintas. La diversidad
aumenta si además pensamos que una misma figura puede variar en la posición del centro y en la
orientación de su eje de rotación. En la XVII General Assembly de la International Union of Geodesy and
Geophysics (IUGG) celebrada en Canberra, Australia, en 1979 se adoptó el elipsoide GRS80. Para su
completa definición se utiliza un parámetro geométrico y tres parámetros físicos. Lo veremos más adelante.
Coordenadas geográficas.
Fig. pág. 50 Maling

La más sencilla entre las aproximaciones cartográficas a la figura de la tierra es la de una esfera que pres-
enta irregularidades del terreno, aunque, a determinadas escalas, estas irregularidades no son de gran
14
mag-nitud. Con un radio terrestre de aproximadamente 6.371 km., la altitud mayor de los continentes se
sitúa próxima a los 9 km. y la profundidad oceánica máxima es del orden de los 11 km. Para la
representación de la tierra a pequeñas escalas y para los cálculos que no exijan precisión es suficiente
considerarla como una esfera. Con ello se simplifican las operaciones cartográficas.
En la aproximación esférica, todos los puntos de la superficie terrestre son equidistantes de su centro. El
tamaño de la tierra vendrá dado por su radio (R) que se define como la línea recta que une el centro con
cada uno de los puntos. Todo plano que divida a la tierra en dos partes iguales, que intersecte a la esfera
pasando por su centro, define en la superficie de esta un círculo de radio R que es denominado círculo
máximo. Utilizando el eje de rotación de la tierra, definimos dos puntos en su intersección con la superficie
de la esfera a los que denominamos polos. El plano perpendicular al eje y que pasa por el centro define un
14
En GRS80 se utiliza, como veremos, una medida de 6 371.008, 7714 metros.
0

círculo máximo llamado ecuador y el conjunto de planos que contienen el eje y son perpendiculares al
ecuador definen entre los dos polos semicírculos denominados meridianos geográficos (también referidos
como meridianos astronómicos o meridianos verdaderos). El término meridiano se limita a un arco de
círculo máximo entre los dos polos y no al círculo máximo completo. Para referir este último habrá que
indicar un meridiano y el correspondiente antimeridiano.
Si elegimos un meridiano de referencia, al que denominaremos meridiano de longitud cero, la posición de
un punto puede medirse mediante el ángulo de su distancia al ecuador, al que llamamos latitud y re-
presentamos por , y el que mide su distancia al meridiano de longitud cero, ángulo al que llamamos
longitud y representamos por . Como todo ángulo, las coordenadas geográficas ( y ) pueden
15
expresarse en grados sexagesimales, grados centesimales o en radianes . Se conviene que la latitud
precede a la longitud. También que las expresiones + y - expresan la latitud norte y la latitud sur. Si,
respecto al meridiano origen, utilizamos las direcciones este y oeste teniendo en cuenta el sentido de la
rotación de la tierra, se conviene, así mismo, que + y - expresan la longitud este y la longitud oeste. Al
o
ángulo complementario a la latitud de un punto (90 -) se le denomina su colatitud y es convencionalmente
representado por .
El conjunto de puntos que tienen un mismo valor de  definen una circunferencia cuyo círculo es paralelo al
ecuador por lo que son denominados paralelos de latitud. Las líneas de longitud, el conjunto de puntos con
el mismo valor deO, son los meridianos. Todos los meridianos geográficos convergen en los polos. El
trazado a intervalos regulares de meridianos y paralelos es uno de los sistemas de mayor implantación para
mostrar la localización de los fenómenos geográficos. Se la conoce, usualmente, como malla geográfica.
Se llaman meridianos magnéticos a los definidos con las direcciones indicadas por las brújulas que varían
según el tiempo y, además, según las localizaciones terrestres. Los meridianos geográficos se representan
en la cartografía, convencionalmente, mediante un asterisco () o una estrella () y los meridianos
magnéticos mediante una punta de flecha (). El ángulo horizontal entre ambos meridianos es denominado
declinación magnética (δ).
Cálculo de distancias y direcciones en la tierra esférica.

La distancia entre dos puntos A y B de una circunferencia se mide como su distancia angular (z), el ángulo
que dichos dos puntos forman con el centro, o como la longitud del arco que definen los dos puntos. En
este último caso, la distancia vendrá dada por:
Longitud del arco: AB = R.z z en radianes
Si z se expresa en grados sexagesimales:
Longitud del arco: AB = 2R.zº/360º
A
R
B
O
z
0
15
La notación sexagesimal es la que usa una esfera que se supone dividida en 360 unidades llamadas grados
(simbolizados por º) que, a su vez, se dividen en 60 minutos (simbolizados mediante ') y estos, por su parte, en 60
segundos ("). La notación centesimal emplea una esfera dividida en 400 grados (g), 100 minutos (m) y 100 segundos
(s). La longitud de los 360º, o de los 400 g del arco de una circunferencia viene dada por 2r. Si tomamos un radio de
valor = 1, definimos el radián como el arco cuya longitud es igual al radio. La longitud de una circunferencia será de 2
radianes.
Esta relación es la que se utilizó para los primeros cálculos del radio de la esfera terrestre obteniendo una
magnitud aproximada de la dimensión de la tierra. La longitud del arco AB puede ser medida sobre el
16
terreno y z puede determinarse mediante observaciones astronómicas .
Distancia entre dos puntos situados en el mismo meridiano: d = R. 0en radianes. La distancia
viene dada por la diferencia de latitudes. El radio del meridiano es el mismo de la esfera.
Distancia entre dos puntos situados en el mismo paralelo: d = R. cos.  , en radianes.
Esta distancia se calcula de forma semejante a la anterior sustituyendo la diferencia de latitudes por la
diferencia de longitudes y el radio de la esfera por el radio del paralelo.
Radio de un paralelo: R0 = R cos
Todo plano que divida a la tierra en dos partes iguales, que intersecte a la esfera pasando por su centro,
definirá en la superficie de esta un círculo de radio R (denominado círculo máximo). La distancia entre dos
puntos cualesquiera se define sobre el arco de círculo máximo que los une. Recibe el nombre de distancia
ortodrómica.
1
Distancia ortodrómica: d = R. cos- sen .sen + cos .cos .cos]
1 2 1 2
ACIMUT (Z): ángulo esférico formado por el meridiano y un círculo máximo en un punto.
Elaborado a partir de Alves, fig. pág. 62.

Dados dos puntos A y B el acimut para llegar a medido desde B es el definido por el círculo máximo que
pasa por A y B y el meridiano de B. El acimut para llegar al punto B medido desde A es el definido por el
círculo máximo que pasa por A y B y el meridiano de A. Se corresponde con el camino más corto entre esos
dos puntos, de aplicación tradicional en navegación marítima y aérea. Cuando los dos puntos se sitúan en
las antípodas (están separados por un arco de 180º), existe una infinidad de rutas ortodrómicas. Por eso,
partiendo del punto A y manteniendo cualquier derrota, se llega al punto antípoda B. En Cartografía se
define convencionalmente la dirección entre dos puntos midiendo el ángulo según el sentido de las agujas
del reloj y utilizando, únicamente, el Polo Norte como referencia. Si el meridiano es un meridiano
astronómico lo denominamos acimut geográfico Zg. Si es el magnético, acimut magnético Zm. La
diferencia entre ambos valores será la declinación (δ).
Zg – Zm = δ
Para el cálculo del acimut empleamos la fórmula:
Ocot Z = cos a.tanb.cosec - sena.cot
La convergencia de meridianos es función de la latitud y puede calcularse como convergencia angular ()
16
Las primeras mediciones de arcos terrestres se realizaban sobre un mismo meridiano o en puntos muy próximos pues
en ese caso la determinación de Z se simplifica. Con el perfeccionamiento de la medida del tiempo, se hizo posible la
determinación precisa de la longitud. El procedimiento no quedó condicionado por la dirección norte-sur. A principios del
siglo XVII (1615) el holandés Willebrord Snell (Snellius) introdujo la medición del arco mediante la observación de los
ángulos de triángulos consecutivos (lo que se conoce como triangulación) que proporcionó mayor precisión en el cálculo
de distancias y en el del radio de la tierra.
o como convergencia lineal (c). Para la definición de c, además de la latitud, habrá que conocer la
distancia entre los meridianos.
 = d.tan /R
c = d.l.tan/R
El rumbo de una distancia ortodrómica no suele ser constante. En los puntos en los que se alcanza la latitud
máxima o mínima el rumbo es igual a 90º o a 270º. Son los denominados vértices de la ortodromia.
Sus coordenadas (V y V) se calculan según:
Ocos V = senR cos1
Ocos (V -1) = cot V. tan 1
El rumbo inicial de una distancia ortodrómica viene dado por:
–1
Orumbo inicial = sen (cos2 sen . sen d)
La distancia descrita por los recorridos que mantienen un rumbo constante se llama distancia loxodrómica
que en el ecuador y sobre los meridianos coincide con la ortodromia. Salvo esos dos casos, la distancia
entre dos puntos a lo largo de la línea loxodrómica es mayor que la distancia ortodrómica.

Elección del radio de la tierra esférica.

En esta aproximación, el tamaño de la tierra vendrá dado por el radio de la esfera (R), la línea recta que une
el centro con cada uno de los puntos de la superficie. Como la tierra no es una esfera, la determinación del
radio dependerá de los objetivos prácticos que se persigan. No se plantean de igual manera los cálculos
necesarios en una proyección o los realizados en navegación, por ejemplo. Maling, utilizando el Elipsoide
17
Internacional como figura de ajuste, proporciona un cuadro en donde se obtienen quince mediciones
distintas del radio que suponen una variación entre el valor máximo y el valor mínimo que alcanza los
34.751 metros. Como parámetros obtenidos del GRS80 la IUGG adopta los que indicamos en color rojo al
lado de los de Maling. Los tres primeros procedimientos se basan en la elipse.
MALING (unidades, metros) IUGG-GRS80 (unidades, metros)
1. Semieje mayor, a: 6.378.388. 6 378 137.0
La esfera suele situarse en este caso siendo tangencial al elipsoide en el ecuador. Es denominada, por
algunos, esfera de Jacobi:
2. Media aritmética de a y b: 6.367.650.

3. Media geométrica de a y b: 6.367.641.
El cuarto y el quinto procedimiento consideran el elipsoide.
18
4. Media aritmética de tres semiejes (2a+b) /3: 6.371.229. 6 371 008.7714
1/3
5. Media geométrica de tres semiejes (2a.b) : 6.371.221.
Podemos utilizar como formas de cálculo del radio terrestre la media de los radios calculados para unas
latitudes de 15º, 45º o 60º.
6. Media geométrica para  = 15º: 6.359.778.
A continuación, utilizamos los cálculos del radio de curvatura del meridiano o del radio de curvatura
transverso a 45º o a 60º de latitud:
9.  para  = 45º: 6.367.586.
10.  para  = 45º: 6.389.135.
11.  para  = 60º: 6.383.727.
12.  para  = 60º: 6.394.529.
Si calculamos el volumen del elipsoide y hallamos el radio de una esfera que tenga un volumen equivalente:
13. Esfera de igual volumen: 6.371.221 6 371 000.7900
Si, en lugar del volumen hallamos el radio de una esfera para que su superficie coincida con la superficie de
referencia, estamos ante la que se conoce como esfera autálica.
14. Esfera autálica: 6.371.228 6 371 007.1810
Una esfera rectificada es aquella cuyos meridianos tienen la misma longitud que los meridianos del
elipsoide.
15. Esfera rectificada: 6.367.655
16. Radio de curvatura polar (a 2 / b): 6 399 593.6259
6 399 593,6258 (WGS84)
17
Sobre el Elipsode Internacional y otros elipsoides anteriores al GRS80 ver apéndice.
18
Radio promedio en la manera definida por la IUGG.
Coordenadas geodésicas.
ELIPSOIDE: figura matemática generada por la rotación de una elipse sobre uno de sus ejes. El elipsoide
D
usado como esferoide tiene como eje de rotación el eje menor de la elipse. Los dos puntos definidos en la
superficie del elipsoide por el eje menor se denominan polos.
Desde que se comprobó el aplanamiento debido a la rotación y que la figura de la tierra, más que una
esfera es un cuerpo próximo a ella (un esferoide), el elipsoide de revolución es utilizado como una aproxi-
mación matemática a la curvatura terrestre que mejora la simple utilización de la esfera, aunque la tierra no
es tampoco un elipsoide exacto. El elipsoide de revolución sigue siendo después de las observaciones
proporcionadas por los satélites artificiales la figura que permite generalizar geométricamente una forma
regular de la tierra. Existen irregularidades, descritas como ondulaciones, pero no son de una gran amplitud.
Además, están las irregularidades del relieve geográfico y la propia asimetría de la tierra.
Para los cálculos geométricos a escalas grandes o medias, en los denominados mapas geográficos, es
necesario el uso de un elipsoide. En las representaciones a pequeña escala no llegan a percibirse por el ojo
humano las diferencias entre una tierra esférica y una elipsoidal siendo suficiente el uso de la primera.
Incluso la topografía clásica adoptaba un plano como figura de la tierra. De todas formas, en cartografía,
siempre se emplean coordenadas referidas sobre el elipsoide (coordenadas geodésicas) aunque se utilice
posteriormente una esfera para los cálculos geométricos o para las operaciones de proyección.
SUPERFICIE DE REFERENCIA GEODÉSICA es la utilizada en las coordenadas que muestran la posición
relativa de los puntos de la tierra. Generalmente se usa un elipsoide.
SUPERFICIE DE REFERENCIA CARTOGRÁFICA es la utilizada en la construcción de proyecciones
usándose el elipsoide, la esfera e, incluso, el plano.
SUPERFICIE DE CÁLCULO. En la navegación los cálculos suelen obtenerse sobre una superficie esférica y
los antiguos topógrafos suelen utilizar un plano, aunque las coordenadas estén referidas a un esferoide.
SUPERFICIE DEL TERRENO y SUPERFICIE OCEANOGRÁFICA: la superficie real o natural de la tierra.
Geometría de la elipse.
La descripción matemática de una elipse viene dada por la ecuación:
2 2 2 2
O x /a + y /b = 1
O
Si OF1 es la distancia a los focos, a es la longitud del semieje mayor de la elipse y b la del semieje menor,
el tamaño o dimensión máxima del esferoide vendrá dado por a. Su forma geométrica la indicaremos por la
desviación respecto a una circunferencia a la que denominamos aplanamiento o elipticidad (f):
elipticidad f = (a-b) /a
O
O que, generalmente, suele referirse mediante su recíproca 1/f

Otra forma de indicar la desviación respecto a una circunferencia es mediante la relación entre la distancia
de los focos y el semieje mayor que es denominada excentricidad (e):
0e = OF1/a,
O que también se define mediante:
2 2 2 2 2 2
e = (a -b ) /a o, sustituyendo, e = 2f - f
La denominada segunda excentricidad (e´), referida al semieje b, viene dada según:
2 2 2 2
e´ = (a - b ) /b
En la circunferencia, el aplanamiento y la excentricidad serán igual a cero (a = b).

b
F1
1
a
O
La distancia al centro de los puntos situados en la superficie del esferoide no es constante como en la
esfera, sino que varía entre un máximo = a y un mínimo = b. Los planos que cortan al esferoide pasando
por su centro definen elipses a excepción del plano que comprende el semieje a, que define una
circunferencia. También los planos paralelos a este último, que cortan al esferoide sin pasar por su centro,
definen circunferencias, pero en el resto de los casos nos definirán elipses.
A diferencia de la esfera, en donde para un mismo valor de z obteníamos valores iguales de la longitud del
arco, independientemente de la posición en la que nos encontráramos, en el esferoide una misma longitud
del arco entre dos puntos de su superficie se corresponderá con una distancia angular que varía en las
diferentes posiciones en las que se establezcan dichos puntos. La curvatura del esferoide varía en cada
punto de su superficie siendo mínima en la circunferencia de radio = a y siendo máxima en los polos.
2
El radio de curvatura polar (c) viene dado por la expresión: c = a /b
La medida del radio de curvatura en cada uno de los puntos del esferoide es necesaria para el cálculo de
distancias. Analizaremos esta cuestión más adelante.
Referencias geográficas sobre el elipsoide.
La longitud en el elipsoide se define de igual manera que en la esfera. En cambio, en la definición de la

latitud se nos presenta una opción. Si en un elipsoide medimos la latitud como el ángulo referido al centro
del mismo, la denominamos latitud geocéntrica representándola por . Si la medimos en la perpendicular
del punto la denominamos latitud geodésica (). En la mayor parte de los casos las coordenadas
19
geográficas suelen indicar la latitud geodésica .
LONGITUD GEODÉSICA (): ángulo formado por el plano meridiano origen con el plano meridiano de un
punto dado, tomado como positivo hacia el este
19
El IGN define,
 Latitud geográfica (φ): ángulo medido sobre el plano meridiano que contiene al punto entre el plano ecuatorial y
la normal al elipsoide en P.
 Longitud geográfica (λ): ángulo medido sobre el plano ecuatorial entre el meridiano origen y el plano meridiano
que pasa por P.
Nosotros, siguiendo criterios aceptados en otros países, preferimos los nombres de latitud geodésica y longitud
geodésica. Las acompañamos de la latitud geocéntrica (). No se induce a posibles errores en los cálculos de
distancias y direcciones geodésicas.
Elaborado a partir de la figura de la pág. 66 Maling.
LATITUD GEODÉSICA (): ángulo formado por la dirección de la perpendicular al elipsoide que pasa por un
punto dado con el plano ecuatorial, tomando como positivo hacia el norte
LATITUD GEOCÉNTRICA (): ángulo formado por el plano ecuatorial y la línea que une un punto dado con
el centro del elipsoide, tomando como positivo hacia el norte
Cálculo de distancias y direcciones en la tierra elipsoidal.

La obtención de distancias y de ángulos sobre el elipsoide terrestre solía evitarse por la complejidad de los
cálculos y dado que las diferencias con los resultados obtenidos sobre la tierra esférica no tenían incidencia
práctica en la mayoría de los casos. Incluso los geodestas utilizan en ocasiones en lugar del elipsoide una
esfera tangente al área del cálculo que suele denominarse esfera local. La informática ha facilitado estos
cálculos.
Distancia d en el mismo paralelo (de radio R0) y con una diferencia de longitud de :
d =  . R0
Distancia a lo largo de un meridiano entre dos puntos de latitudes 1y 2:
2 4 2 4
d = a {(1– 0,25e –0,046875 e ) - (0,375e +0,09375 e ) (sen 2.2 - sen2.1) + 0,05859375
4
e (sen 4.2 - sen 4.1)}
El tercer tipo de distancia, la distancia más corta entre dos puntos cualesquiera P 1 (1, 2) y P2 (2, 2) es la
denominada distancia geodésica (dg):
Distancia dg entre dos puntos cualesquiera P1 (1, 2) y P2 (2, 2)
dg =ʃR m+R pcos  (d/d )  d
2 2 2 2 1/2

Siendo Rm y Rp los radios de curvatura del meridiano y el radio de curvatura transverso

Evolución de las determinaciones de las superficies de referencia20.

PRINCIPALES DETERMINACIONES DE LA FIGURA DE LA TIERRA
Longitud de los semiejes

------------------- -------------------------
Fecha Nombre a b Elipticidad Utilización
_____ ________ _________ _______ _____________ ______________________________
1687 Newton l/231 -

1738 Maupertuis 6397300 l/191 -
1749 Bouguer l/266 -
1799 Commission des Poids
et Mésures 6375738.7 l/334.29 -
1810 Delambrea 6376428 l/311.5 anteriormente en Bélgica
Delambreb 6375653 6356564 1/334.0
1817 Plessis 6376523 6355863 l/308.65 anteriormente en Francia
1819 Walbeck 6376896 6355833 1/302.8
1830 Everest 6377276.345 6356075.413 1/300.8017 India, Birmania, Ceilán, Malasia
(parcialmente)
1841 Bessel 6377397.155 6356078.963 1/299.1528128 la mayor parte de Europa Central, Chile,
China, Asia Sudeste
1849 Airy 6377563.396 6356256.9 1/299.3249646 Gran Bretaña
1858 Clarke 6378294 6356618 l/294.261 anteriormente en Australia
1860 Struve 6378298,3 1/294.73 España
1866 Clarke 6378206.4 6356585 l/298.9786982 Norteamérica (NAD 27)
1876 Andrae 6377104.43 6356762 l/300.0 Dinamarca y anteriormente en Islandia
1880 Clarke 6378249.145 6356514.9 1/293.465 Francia y la mayor parte de África
continental
1893 Zhdanov 6377717 6356433 1/299.6
1907 Helmert 6378200 6356818 l/298.30 anteriormente en Egipto
1909 Hayford o
Internacional (1924) 6378388 6356911.946 1/297.0 Todo el mundo excepto Norteamérica,
África y otras áreas pequeñas
1940 Krasovsky 6378245 6356863.019 1/298.3 antigua URSS y algunos de sus países
satélites
1956 Army Map Service 6 378 270.0 6 356 794.343 l/297.0
1960 Fischer (Mercury Datum) 6 378 166.0 6 356 784.284 l/298.3
1960 World Geodetic System
(WGS60) 6378165 6356783 1/298.3 anteriormente en Australia
1961 Kaula 6378163 l/298.24
1965 Internacional Astrono-
mical Union (IAU65) 6378160 6356774.719 1/298.25 Australia
1965 Applied Physics Labora-
tory (APL 4-5) 6378137 l/298.25
1965 Naval Weapons Labora-
tory (NWL9D) 6378145 6 356 759.8 l/298.25
(WGS66) 6378145 6356 760 1/298.25
1967 Geodetic Reference
System (GRS67) o 6378160.00 6356774.161 1/298.2472
International Astrono-
mical Union (IAU68)
1968 Fischer (Modified
Mercury Datum) 6378150 6 356 768.955 l/298.3
1972 World Geodetic
System (WGS72) 6378135 6 356 750.52 l/298.26
1979 Lerch el al. 6378139 l/298.257
1980 Geodetic Reference
System (GRS80) o 6378137.00 6356752.3141 1/298.257222101 Norteamérica (NAD 83) y de creciente
Nuevo Internacional implantación en los países de la OTAN y
en utilizaciones con GPS
(WGS84) EPSG:4326 6 378 137.00 l/298.257223563 De creciente implantación en los países de
la OTAN y en utilizaciones con GPS
1985 Engelis 6 378 136.05 l/298.2566
a: según Maling
b: según Maslov
Las principales figuras usadas en cartografía básica se refieren en negrita.
20
“A lo largo de la historia diversos elipsoides se han utilizado para definir el sistema de referencia de cada país, de tal
forma que se define aquel que mejor se ajuste al geoide.” según dice el IGN, aunque debemos considerar que no siempre se
utilizaron criterios científicos. El GRS80 se ha adoptado internacionalmente como la figura de referencia y los numerosos
sistemas estatales o nacionales van cayendo en desuso. La ICAO (International Civil Aviation Orgamization) los denomina
datums y sistemas locales cuyo origen y orientación de los ejes son arbitrarios.
Elipsoide equipotencial GRS80. Determinación de elipsoides de ajuste.

Mediante la observación de las órbitas de los satélites artificiales y aquellas mediciones astronómicas que
hemos enunciado, se ha llegado a la determinación del elipsoide conocido con las siglas GRS80 (Geodetic
21
Reference System 1980) que es el elipsoide geocéntrico de la Unión Internacional de Geodesia y
Geofísica (IUGG) para la georrefenciación adoptada internacionalmente. Hasta el sistema GRS80, los
elipsoides se definían por sus parámetros geométricos. El radio ecuatorial a (con el que se establecía su
tamaño), el aplanamiento f, el radio polar b, la relación b/a, eran parámetros geométricos primarios o
22
derivados. A partir del GRS80 la IUGG adoptó magnitudes físicas y así se eligieron inicialmente:
0a, GM, J2 y ω.
La elipticidad y el resto de magnitudes geométricas pasaron a ser magnitudes derivadas. También se
derivaron otras magnitudes físicas. De esta manera, además de un único parámetro geométrico (a), el
sistema de referencia atiende a las propiedades físicas de la tierra como son su masa total geocéntrica,
incluida la atmósfera (M), el factor de la forma dinámica (J2), parámetro de relativa complejidad, y la
velocidad angular de la rotación sobre su eje menor (ω). Son estos el parámetro geométrico y los
parámetros físicos básicos o primarios. Pero la constante gravitacional (GM) es un parámetro conocido con
mayor precisión que la constante gravitatoria (G) o la masa de la tierra (M). Se especifica, también, que el
sistema es geocéntrico coincidiendo con el centro de masas. La orientación y el meridiano de referencia
siguen los criterios ITRS que ya hemos visto.
ELIPSOIDE GEODÉSICO: elipsoide de rotación, normalmente elegido para que se ajuste a la forma de la
tierra tanto como sea posible, ya sea global o localmente. Por ello, distinguimos:
ELIPSOIDE GLOBAL es el que adapta el Geoide al conjunto de la tierra. Actualmente se emplea el
GRS80. Es el convenido internacionalmente y el que utilizaremos en el cálculo de coordenadas.
ELIPSOIDE REGIONAL o LOCAL es el que ajusta la forma de la tierra sólo en un ámbito terrestre.
La tierra tiene una superficie asimétrica, irregular y ondulada. En cambio, el elipsoide es una figura
geométrica regular. Por ello, el elipsoide que mejor ajusta en una región no es necesariamente el que mejor
ajusta en otra. El primer elipsoide internacionalmente aceptado, fue el que, por ello, es conocido como
Elipsoide Internacional de 1924. Cuando utilizábamos un elipsoide regional, el ajuste podía superar diferen-
cias de 300 m.
Muchos servicios cartográficos estatales continúan utilizando elipsoides que difieren en forma y tamaño,
orientación y posición relativa con el GRS80 (Nuevo Internacional) generando innumerables problemas de
transformación de coordenadas. Las ventajas de un único sistema internacional son muy claras.
De todas formas, a grandes escalas como en la ingeniería y los catastros, deben utilizarse elipsoides regio-
23
nales para que no se introduzcan grandes distorsiones en las mediciones .
ELIPSOIDE GEOCÉNTRICO es aquel cuyo centro geométrico se corresponde con el centro de masas de la
tierra. Actualmente se usa el GRS80 como elipsoide geocéntrico.
24
DATUM (plural data): superficie con valores de coordenadas constantes. Las coordenadas del GRS80
son las que hacen coincidir su centro con el centro de la tierra. Su eje menor también coincide con el eje de
rotación. En esta situación, es el datum convenido actualmente.
21
El WGS84 Development Committee adoptó el GRS80. Después de la publicación inicial, introdujo ligeras
modificaciones que son de importancia cuando queremos una elevada precisión en el cálculo de las órbitas de los satél-
ites pero que no intervienen significativamente en cartografía geográfica.
22
La resolución nº 7 de la XVII Asamblea General de la IUGG adopta un elipsoide equipotencial geocentrico que tiene
como constantes:
Radio ecuatorial de la tierra a = 6.378.137 m, constante gravitacional geocéntrica (incluyendo la atmósfera) GM =
8 3 -2
3.986.005 x 10 m s , factor de la forma dinámica de la tierra, excluyendo la deformación permanente de las mareas: J2
= 108 263 x 10 , velocidad angular de la tierra:  = 7292 115 x 10 rad s .
-8 -11 -1
Además de la normalización en el uso de fórmulas de cálculo, se adopta que el eje menor del elipsoide de referencia
será paralelo a la dirección establecida por el Conventional International Origin y que el primer meridiano será paralelo
0
al meridiano de longitud cero del BIH.

23
Las0normas internacionales hablan de0un Datum de Ingenieria. El Catastro no utiliza en España0un sistema propio.
24
En Ordnance Survey (2002) se dice, contundentemente, en pág.15: “Datum is the most familiar term amongst
surveyors (…). TRS is a more modern term for the same thing”.
25
La definición en INSPIRE no es tan directa pues está afectada por el deseo de interoperabilidad.
DATUM (INSPIRE): un parámetro o conjunto de parámetros que definen la posición del origen, la escala y
la orientación de un sistema de coordenadas.
Al igual que en ISO 19111, también se concreta el Datum Geodésico.
DATUM GEODÉSICO (INSPIRE): un datum que describe la relación de un sistema de coordenadas con la
tierra.
Para los elipsoides regionales:
DATUM REGIONAL: conjunto mínimo de parámetros necesarios para la definición de la localización y la
orientación de un sistema regional o local respecto a un sistema y un marco globales.
De forma más restrictiva,
DATUM HORIZONTAL: usualmente, elipsoide determinado con una localización terrestre utilizado en la
determinación de las coordenadas  y . Utilizaremos el GRS80, determinado como elipsoide geocéntrico
horizontal (en realidad es tridimensional).
26
DATUM VERTICAL es la superficie utilizada con elevación o altitud cero. Conjunto de datos de elevación
cuya precisión permite utilizarlos como referencia de otros datos altimétricos.
En la actualidad, la determinación del datum se construye mediante observaciones que definan su posición,
su orientación y la deflección de la vertical o diferencia entre las normales en el geoide y el elipsoide, de tal
manera que sea igual a cero en un punto origen. Tradicionalmente, el datum quedaba fijado con las
coordenadas del origen (posición), la orientación, y los parámetros del elipsoide. El datum vertical se
obtenía mediante observación de las variaciones de las mareas en un período de tiempo. La observación
podía ser directa o utilizando el nivel medio del mar de uno o más puntos en el citado período de tiempo.
Sistemas ETRS89 y WGS84.

“A partir de las series temporales de resultados del IERS, se ha puesto de manifiesto que la Placa
Continental Europea mantiene un movimiento bastante uniforme, de unos 3 cm por año, con relación al
ITRS, con excepción del extremo sureste de Europa (Grecia, Turquía). Por esta razón, con el fin de
mantener unas coordenadas razonablemente estables para Europa, la Subcomisión EUREF decidió definir
un Sistema ligado a la placa europea. Este sistema (datum) se denomina ETRS, o ETRS89, ya que fue
idéntico al ITRF en el año 1989,0. Desde entonces, las coordenadas ETRS89 ajustadas con relación a la
Placa Europea, han modificado sus valores con respecto a los expresados en ITRF. Sin embargo, esta
modificación es bien conocida, controlada por IERS y EUREF, y son posibles las transformaciones entre
unas y otras con exactitud de 1 cm para la mayor parte”. Así pues, el control en Europa está ligado a la
parte estable de la placa europea nombrándose como ETRS89 (European Terrestrial Reference System
1989).
27
EPN (EUREF GNSS PERMANENT NETWORK) red de estaciones con coordenadas de elevada precisión
que se utilizan para la actualización permanente del sistema ETRS89. A cargo de EUREF, subcomisión de
IAG.
25
(REGLAMENTO (UE) Nº 1089/2010 DE LA COMISIÓN). en conformidad con EN ISO 19111. También en INSPIRE
D2. 8.I.1_v3.2 apartado 2.3.0con referencia más detallada0a EN ISO 19111:2007, Geographic information — Spatial
referencing by coordinates.
26
Además, en ISO 19111 se dice: “un datum para la ingeniería es el que no es geodésico ni vertical”. De momento,
:0
podemos no considerarlo en geografía. Las definiciones de datum en ISO son las del cuadro:
Datum parameter or set of parameters that define the position of the origin, the scale, ISO 19111:2007(E)
and the orientation of a coordinate system 19130:2010(E)
(Es la definición que coincide con INSPIRE) 19136:2007(E)
Datum parameter or set of parameters that serve as a reference or basis for the ISO 19116:2004(E)
calculation of other parameters
27
El EPN Central Bureau está en el Royal Observatory of Belgium, http://www.epncb.oma.be.”
“A six-parameter transformation between ETRS89 (European TRS 1989, the Europe-fixed version of WGS84) and the
various versions of the ITRS (International Terrestrial Reference System) is published by IERS (International Earth
Rotation Service) via their Internet site. (…) These transformations are not suitable for ITRFs before ITRF94 (these are
no longer in common use)” Ordnance Survey (2000), pág. 33.
ETRS89 emplea el elipsoide GRS80 cuyo centro coincide con el centro de masas y es el origen del sistema
de coordenadas X, Y, Z. Es el elipsoide utilizado también por el sistema americano WGS84 (World Geodetic
System 1984), adherido al anterior con ligeras modificaciones.
28
Sus elipsoides (datums de superficie) tienen los siguientes parámetros geométricos
Oa 1/f Obb
Elipsoide
GRS 80 6 378 137,0 m 298,257 222 100 882 711 243 ≈ 6 356 752,314 140 m
WGS 84 6 378 137,0 m 298,257 223 563 ≈ 6 356 752,314 245 m
con una pequeñísima diferencia en f, por problemas de cálculo, que da lugar a una ligera modificación en el
O O
semieje polar. También en el WGS 84 se utiliza el IERS Reference Meridian y se calculan, como constantes
O O
geométricas derivadas:
O f = 0, 003352810681183637418;
O relación entre ejes b / a = 0, 996647189318816362;

2 2 5
O primera excentricidad (a − b ) = 521.854,0097 m.;
2 2
segunda excentricidad a − b / a = 0,0818191910435;
Os egu ntricidad
cuadrante de un meridiano = 10.001.965,7293 m.

Os egu
Además, se obtienen varias medidas del radio de la tierra esférica, como ya hemos visto.
29
En la actualidad, ha mejorado el valor del semieje a pero se ha preferido mantener el valor anterior para
no tener que recalcular las observaciones, La elipticidad pasa a ser uno de los parámetros de la definición
30
del elipsoide que fue modificada en WGS84 con el cálculo de nueve decimales .
28
Puede verse, por ejemplo, en www.ngs.noaa.gov (XIONG LI & HANS-JÜRGEN GÖTZE (2001)).
29
Las nuevas estimaciones aplican “satellite radar altimeter data, satellite laser ranging data, and Doppler satellite data.
Results from these determinations (…), range from 6378134.9 to 6378137.0 meters”, en NGA, STND.0036 (2014), pag.
2-4. En la práctica el cálculo de coordenadas mediante GNSS no requiere un elipsoide mejorado.
30
Como dicen de sí mismos: “The WGS 84 represents the best global geodetic reference system for the Earth available
at this time for practical applications of mapping, charting, geopositioning, and navigation”, NGA.STND.0036 (2014),
pág.1-1.
2. GEORREFERENCIACIÓN Y REGISTRO DE LA ALTIMETRÍA.
FIGURAS DE REFERENCIA. GEOIDE.

Ya hemos hablado de algunas superficies geométricas de referencia terrestre entre las que destacamos el
elipsoide. Es la superficie convencional en la que se definen las coordenadas latitud y longitud ( En ella
se calculan operaciones geométricas como son las distancias y las direcciones entre puntos o como es el
cálculo de algunas superficies. Sirve para indicar la orientación (está orientada al N) y para las
transformaciones geométricas de traslación, rotación y cambio de escala. En estas operaciones son
necesarias las coordenadas sobre un elipsoide que se utiliza como una figura regular de ajuste. Con menor
precisión, también se usan la esfera o, incluso, los topógrafos han venido usando el plano (el plano
horizontal o plano tangente).
La tierra no es una figura geométrica regular (no es plana ni tampoco es lisa). En la descripción de su forma
y tamaño se utilizaron superficies físicas del potencial de gravedad (W) debido a las insuficiencias de los
elipsoides geométricos. Desde finales del siglo XVIII ya se había pensado que el nivel medio de los
océanos, en las áreas no cubiertas por el hielo, eliminando los efectos de las corrientes marinas, del oleaje
o de la acción del viento, es parte de una superficie utilizable como inicial figura de la tierra. En los océanos,
evaluados en un 70%, se observaba un equilibrio en las acciones de la fuerza centrífuga, debida a la
rotación, y de la fuerza de la gravedad, debida a la masa terrestre. Dicho de otra forma, se suponía que los
océanos eran parte de una superficie con potencial de gravedad constante. A esta superficie equipotencial
pero geométricamente irregular la denominábamos geoide utilizando un término introducido por J. B. Listing
en 1878 y la denotamos con W 0 pues se empleaba de referencia inicial. Las direcciones de la plomada,
también llamadas líneas verticales (aunque son curvadas), cruzaban perpendicularmente a esta superficie
en toda su extensión. Después, la prolongábamos idealmente bajo el interior de los continentes, también
con la referencia W 0, aunque entonces su descripción era más difícil. Teníamos que efectuar mediciones del
potencial de gravedad en cada punto y reducirlo al que tendría al nivel del mar. Se suponía, falsamente, el
conocimiento de la densidad de la corteza para esa reducción.
Comparándolo con el elipsoide, aparecía con ondulaciones que señalaban la distinta distribución de masas
en los distintos puntos de la tierra. El agua de los océanos describe un trazado que no es plano pues es
continuo, cerrado y convexo, lo podemos ajustar mediante un elipsoide, pero describiendo las ondulaciones
(sin aristas) que señalan sus diferentes formas. Lo mismo en su prolongación en los continentes.
Precisamente, la representación del geoide se expresaba por su separación con respecto a un elipsoide de
referencia que es la que se conoce como ondulación geoidal N.
GEOIDE: superficie física equipotencial del campo de gravedad terrestre que ajusta el nivel promedio del
31
mar y que se expresa mediante ondulaciones sobre el elipsoide . Suele utilizarse como superficie de
referencia (W 0, nivel = 0) en la medición de altitudes.
“Debido a que la figura del geoide depende de la distribución de la masa en el interior de la tierra, desco-
32
nocida por nosotros; rigurosamente hablando es indeterminable” . Las irregularidades geológicas hacen
difícil su descripción en términos matemáticos o físicos. Además, varía temporalmente por la modificación
de las densidades internas de la tierra y por la redistribución de masas en su superficie (particularmente las
redistribuciones hidrológicas). Teóricamente, existe un solo geoide, aunque en su determinación existan dis-
tintas aproximaciones y versiones. Es una figura de referencia dinámica, en evolución, a diferencia de las
superficies geométricas regulares a las que suponemos con estabilidad en sus parámetros.
Pese a sus innegables imprecisiones, se conserva el nombre por su significado físico teórico y por su histo-
ria geodésica. Por inercia, en cartografía podemos continuar llamando geoide a los datos sobre la gravedad
que definen un valor equipotencial W 0 y que pueden utilizarse de referencia inicial para el cálculo de
altitudes, aunque no estén obligatoriamente calculados en el interior de las masas continentales. A veces,
este conjunto de datos se indica en la bibliografía como “geoide” (con entrecomillado). En todo caso,
siempre indicaremos las ondulaciones sobre el elipsoide, como decimos. El potencial del geoide lo
31
La definición de la NGA es más reducida. Se accede en: http://earth-info.nga.mil/GandG/coordsys/definitions. html. La
definición de NOAA en http://www.ngs.noaa.gov/corbin/class_description/Geodetic_Tidal_Datums 0811. shtml.
32
ZAKATOV, P. (1997), pág. 14.
Georreferenciación y registro de la altimetría. Figuras de referencia. 23
2 2 2 2
expresamos en m /s . Por ello, dos superficies equipotenciales diferentes en 1 m /s están separadas en
2 2 -2 -2
torno a 0.1 m dado que la gravedad es de unos 9.8 m/s (0.1 m = 1.0 m s / 9.8 m s ).
Pese a todo esto, el estudio de las densidades interiores tiene otras aplicaciones que van más allá de lo car-
tográfico y no deben abandonarse. En los estudios geofísicos, los datos de gravedad se utilizan en la dis-
tribución de masas en el interior de la tierra, así como su cambio con el tiempo. En los estudios geo-
dinámicos se centran en fenómenos como la isostasia, los movimientos tectónicos, las variaciones en la
rotación de la tierra y, especialmente en los últimos años en nuestro país, en las deformaciones volcánicas.
Entre las aplicaciones metrológicas, se encontraban también la determinación del patrón de masa, el kilo-
gramo, o la determinación de la constante de gravitación universal (G).
ONDULACIÓN O ALTURA DEL GEOIDE distancia + N o -N que separa al geoide del elipsoide de ajuste
(GRS80) en cada uno de sus puntos. Se mide en la perpendicular geométrica al elipsoide.
Ondulación del geoide
Ondulación del geoide
Perpendicular al elipsoide
Línea de plomada, perpendicular al
geoide; curva de máxima pendiente
0
Deflección de la vertical
Globalmente, la diferencia en la elevación respecto al elipsoide geocéntrico puede situarse en + 100 m. Esta
diferencia responde a la llamada ondulación del geoide que denotamos con + N. En España, los valores son
positivos siempre.
“Molodensky en 1941 fue capaz de demostrar que la superficie física de la tierra puede determinarse a partir
de medidas geodésicas solas, sin usar la densidad de la corteza terrestre. Esto requiere que el [antiguo]
concepto de geoide sea abandonado… Ahora las anomalías de la gravedad y las desviaciones de la vertical
se refieren al terreno y no al nivel del mar; las anomalías de la altitud al nivel del terreno toman el lugar de
33
las ondulaciones del geoide” .
GRAVEDAD (g): combinación de la fuerza centrífuga (debida a la rotación) y de la gravitación (causada por
2 2
la masa). Su valor se incrementa desde ~9.78 m/s en el ecuador hasta ~9.83 m/s en los polos. En cada
punto, esos valores se modifican según la densidad de las masas interiores y según la altitud topográfica.
Es la fuerza que interviene en los cálculos de las superficies equipotenciales y de las altitudes. Se mide con
gravímetros.
ELIPSOIDE EQUIPOTENCIAL MEDIO TERRESTRE O ELIPSOIDE DE NIVEL.

Cuando en las observaciones utilizamos el elipsoide, no se puede alcanzar una gran precisión altitudinal
pese a sus ventajas como figura geométrica que permite otros cálculos. Se debe a las diferencias
elipsoidales en el ajuste a la tierra y a que la normal geométrica no coincide con la perpendicular física
sobre el geoide, con la línea de plomada o curva de máximas pendientes a las superficies de nivel. A esta
diferencia la llamamos deflección de la vertical. “Cuando nombramos una superficie de nivel estamos
33
HEISKANEN, W y MORITZ, H. (1985), págs. 287-288. En apéndice tratamos algo más extensamente la teoría de
Molodensky, verdadera “revolución copernicana” en geodesia.
hablando de una superficie que es perpendicular a la dirección de la gravedad en todos sus puntos. Se
supone que, generalmente, la dirección de la gravedad es hacia el centro de la tierra, pero con variaciones
complejas en dirección y magnitud según los emplazamientos, aunque estemos a escala local o a muy
reducida escala. Estas variaciones, muy pequeñas y que exigen un equipo de medida especializado, se
deben a la irregular distribución de masas sobre la superficie terrestre (valles y colinas) y también a las
variaciones de densidad interna. Por lo tanto, todas las superficies de nivel son en realidad irregulares
34
(bumpy) y complejas” . Las diferentes líneas de plomada no son líneas paralelas como un teórico de la
geometría nos diría. Las líneas de plomada de todos los puntos de la tierra, aunque no se cortan, confluyen
en el centro de masas y tampoco son líneas rectas pues deben sus trayectorias a la distribución de las
densidades en el interior. Es decir, en las mediciones habrá que tener en cuenta la distribución de masas
del planeta y su campo gravitatorio.
DEFLECCIÓN DE LA VERTICAL (ξ, η): ángulo formado en cada punto por las direcciones de la gravedad
(dirección de la plomada) y la perpendicular geométrica al elipsoide. Tiene dos componentes: uno norte-sur
(ξ) y otro este-oeste (η).
35
ELIPSOIDE EQUIPOTENCIAL o ELIPSOIDE DE NIVEL : elipsoide definido para ser equipotencial en su
superficie.
Si utilizamos el elipsoide geométrico como referencia para las coordenadas geográficas  y , le asignamos,
además, una superficie equipotencial (deberá contar con iguales valores del potencial gravitacional para
todos sus puntos). Dado el elipsoide de ajuste GRS80, le adjudicamos una superficie equipotencial, U =
constante. La función potencial U es llamada gravedad potencial normal (teóricamente). Está determinada,
únicamente, mediante la superficie del elipsoide (semiejes a y b), la masa de su cuerpo sólido M y la veloci-
dad angular ω. Para que sirva de forma sencilla de comparación, no intervienen las diferentes densidades
36
internas .
El campo gravífico de la tierra lo descomponemos en un campo normal, matemáticamente simple, y un
campo perturbador. Para el primero utilizamos, arbitrariamente, el elipsoide GRS80 sabiendo que las des-
viaciones de su gravedad asignada respecto a los valores reales (las perturbaciones) son pequeñas. “El
campo de gravedad de un elipsoide es de importancia práctica fundamental porque es fácil de manejar
matemáticamente y las desviaciones del campo gravífico real respecto al campo elipsoidal normal son tan
pequeñas que pueden considerarse lineales”. De otra forma, el problema geodésico sería de difícil solución;
el elipsoide va a ser la forma normal del geoide. Todavía se dice que si el elipsoide de nivel tiene una super-
37
ficie U = constante, “corresponde exáctamente al geoide, definido como una superficie W = constante” .
ELIPSOIDE MEDIO TERRESTRE: elipsoide equipotencial en el que la masa M, la diferencia en los
momentos de inercia y la velocidad angular w son las de la tierra. La gravedad potencial U=W 0.
Desde hace unos 250 años en geodesia se utiliza la teoría del elipsoide equipotencial, geocéntrico, con la
misma masa que la tierra, rotando con la misma velocidad angular y con su eje mayor coincidente. La tierra
se ajusta a una figura regular en la que no existirían océanos ni montañas, no existirían diferencias en las
densidades de las rocas ni tampoco diversas composiciones en su interior (no existirían la corteza, el manto
o el núcleo, por ejemplo). Ese elipsoide, llamado elipsoide normal, está asociado a un campo de gravedad
teórico conocido también como campo de gravedad normal. Su superficie es equipotencial y puede usarse
como superficie de referencia en los cálculos de la magnitud y dirección de las coordenadas de la posición y
de la gravedad.
De esta manera, el elipsoide equipotencial GRS80 se define incluyendo la masa total de la tierra y su
atmósfera y se utiliza en dos ocasiones: 1) como superficie de referencia para las coordenadas cuyos cálcu-
los no son geométricos sino físicos y 2) para la obtención de medidas gravimétricas y de geodesia satelital.
Es decir, también se utiliza en el cálculo de la latitud y la longitud para los puntos de la superficie terrestre y
34
Traducimos de ORDNANCE SURVEY (2002), pág. 5.
35
“La teoría del elipsoide equipotencial se debe, inicialmente, a Pizzeti en 1894. Fue elaborada posteriormente por
Somigliana en 1929. Esta teoría sirvió, además, como base de la International Gravity Formula adoptada en la
Asamblea General de Estocolmo en 1930”, traducimos de MORITZ, H. Geodetic Reference System 1980, pág. 129.
36
Para todo ello, también se encuentra un resumen de referencia en MORITZ, H. op. cit. Traducimos de la misma
página 129, “de acuerdo con el teorema de Stokes-Poincaré, completamente independiente de la distribución de la
densidad interna”. Sobre estos mismos asuntos, ICAO (2ª ed. 2002), apartado 3.1.
37
HEISKANEN, W. y MORITZ, H. (1979), op. cit. pág.64. El subrayado es mío. Además, he utilizado la página 81.
en la geodesia que utiliza satélites artificiales. Se define un campo de gravedad normal (teórico o
convencional) en la superficie de la tierra e interviene en los cálculos sin contemplar tampoco las
38
variaciones en la densidad atmosférica . El GRS80 es el que mejor ajusta a la tierra y el que se utiliza en
las mediciones con GPS, aunque desde su adopción como datum, los movimientos de las placas han
cambiado la superficie de la tierra (como ya habíamos visto).
En resúmen, los parámetros del GRS80 modificados según el WGS84 son los del cuadro:
39
Parámetros notación magnitud
________________________________________________________________ _______________________________________________________ _________________________
Semieje mayor a = 6 378 137, 00 metros

Elipticidad 1/f = 298. 257 223 563 adimensional
a estos parámetros geométricos se les asocian los parámetros del potencial de gravedad del geoide:
0
Velocidad angular  = 7 292 115,010 -11

rad/s
y el parámetro que ya habíamos indicado en la determinación ITRF
0
8 3 2
Constante gravitacional incluida GM = 3 986 004,41810 m /s
la masa de la atmósfera terrestre
0
+07 2 2
Potencial de gravedad normal U0= 6,26368517146 × 10 m /s
CÁLCULO DE LA GRAVEDAD NORMAL.

Definiendo de esa manera las características físicas del elipsoide de ajuste, al gradiente de la función
potencial normal (U) lo denominamos gravedad normal (γ) cuya magnitud sobre la superficie (cuando h=0)
tiene un cálculo sencillo (fórmula de Somigliana):
0donde
2 2 -5
Las unidades en m/s que en miligales (abreviadamente mgal) las obtenemos con: 1 m/s = 10 mgal.
Según esa fórmula, el valor de la gravedad en un punto de la superficie del elipsoide depende de su forma
y tamaño (representados por a, b y e en las fórmulas anteriores) y de unos valores de la gravedad polar y
ecuatorial (aceptados por convenio). Además, varía en función de la latitud.
Parámetros de la gravedad en el elipsoide GRS80 aceptados internacionalmente:
e= gravedad normal en el ecuador = 9,7803253359 m/s2,
p= gravedad normal en el polo = 9.8321849379 m/s2.
Constante k = 1.931852652458 × 10-03
U0 = 6,26368517146 × 10+07m2/s2
 = gravedad promedio = 9,7976432223 m/s2
Con estos parámetros, hemos elaborado en apéndice una hoja de cálculo con la que obtenemos la
gravedad normal para cada punto definido por la latitud.
Si en lugar de mediciones geométricas, utilizamos la gravedad:
LATITUD ASTRONÓMICA (): ángulo formado por la dirección de la gravedad en un punto dado con el
plano ecuatorial, tomado como positivo hacia el norte.
LONGITUD ASTRONÓMICA (): ángulo formado por el plano meridiano origen con el plano meridiano
celeste de un punto dado, tomado como positivo hacia el este. El plano meridiano celeste es el plano que
contiene las direcciones del eje polar y de la gravedad en un punto dado.
Cuando la deflección de la vertical es cero, la latitud y la longitud astronómicas son iguales respectivamente
a la latitud y la longitud geodésicas.
38
“Si se deben tener en cuenta los efectos atmosféricos, han de aplicarse correcciones a los valores medidos”,
traducido de NGA: STND.0036_1.0.0_WGS84, pág. 3-2.
39
De esos parámetros se han derivado estos otros: Semieje menor b: 6.356.752,3142 m, Primera excentricidad e:
-2 2 -3
8,1819190842622*10 , Cuadrado de la primera excentricidad e : 6.69437999014*10 , Coeficiente armónico zonal de
-3
segundo grado del potencial gravitacional (normalizado) C: 2-0. 48416677498510 . Sobre estos y otros parámetros
ver: MORITZ, H. “Geodetic Reference System 1980”, citado en apéndice.
ALTITUD GEOGRÁFICA Y DATUM VERTICAL.
La posición de un punto quedará descrita con las coordenadas tridimensionales longitud y latitud (,) y con
el valor de la altitud para el que empleamos los símbolos h (altura geométrica) y H (con distintos procedi-
mientos para el cálculo de las alturas gravimétricas). La latitud y la longitud hay que completarlas con la
distribución de las altitudes que es otro componente fundamental en los análisis geográficos, aunque no sea
una coordenada sino un valor distinto. La obtención de un modelo topográfico, de una referencia de las
elevaciones, es una de las primeras tareas del cartógrafo. Sólo teniendo en cuenta los valores de la altitud
podremos cartografiar acertadamente la información geográfica. De esa manera, la latitud y la longitud, los
parámetros de la figura de ajuste y la altitud son los valores necesarios para una completa
georreferenciación. Hablamos, engañosamente, de coordenadas planas o bidimensionales y de coordena-
das tridimensionales. Aunque son dos coordenadas angulares, la latitud y la longitud no son bidimen-
40
sionales, como habíamos visto .
Puesto que la altitud es necesaria para fijar la posición terrestre de un punto, llega a decirse que es la ter-
cera coordenada geográfica, aunque esto sea impropio pues es un valor numérico diferente. Tan sólo es un
dato que expresa el vector de la elevación. En el sistema geográfico, dada la posición de un punto () po-
demos conocer la latitud y la longitud de un segundo punto situado a cualquier distancia (n), (n) o
(nn). Pero dada la altitud h de un primer punto y una distancia no podemos calcular la altitud del se-
gúndo punto.
Si bien la altitud no es la tercera coordenada, la georreferenciación ha de completarse con la descripción de
la superficie topográfica. Los valores de la elevación del terreno son la base de aquellos fenómenos geográ-
ficos en los que no basta la simple referencia a un elipsoide. Únicamente con las coordenadas geográficas y
sus correspondientes coordenadas proyectadas, la tierra se reduciría a una figura regular de superficie nive-
lada o lisa. Para su descripción, lo usual en cartografía es el registro del valor de la altitud asignándolo a las
coordenadas utilizadas en la localización.
De esta manera, a las coordenadas geográficas, latitud y longitud, asociamos posteriormente los valores del
relieve terrestre con el cálculo de la altitud para que no sea una figura geométrica regular. Aunque para los
cálculos de distancias y direcciones son necesarias coordenadas sobre un elipsoide, una esfera o, incluso,
sobre un plano, la distribución geográfica de las altitudes es otro componente fundamental en los análisis
geográficos. A partir de la altitud calcularemos la orientación, las pendientes y otros valores del relieve.
Es indudable que el terreno es la superficie sobre la que se describen los fenómenos físicos y humanos del
área emergida. En ella, además de la altitud y de sus magnitudes derivadas (pendiente y su variabilidad,
orientación, escorrentía), realizaremos el análisis hidrológico, geológico o geomorfológico. En geografía es
fundamental la referencia al relieve, a las altitudes, a la topografía. Tanto en las actuaciones de gestión o de
ingenieria territorial como en el entendimiento científico o cultural de lo geográfico. No es casual que los
nombrados como mapas básicos se conozcan como mapas topográficos indicando que las elevaciones son
41
el fundamento, la base, de los restantes elementos representados .
SUPERFICIE DEL TERRENO y SUPERFICIE OCEANOGRÁFICA: en conjunto, son dos superficies que
forman la superficie real o natural de la tierra. Están establecidas mediante puntos cuyas coordenadas geo-
gráficas se miden respecto al elipsoide ( ) y para cuya altitud y batimetría utilizamos una superficie de
referencia llamada datum vertical.
40
http://inspire.ec.europa.eu/documents/Data_Specifications/INSPIRE_DataSpecification_RS_v3.2.pdf
A diferencia de lo que suele considerarse bidimensional bajo criterios espaciales, en ISO 19107 se define: Dimension.
The attribute "dimension" is the length of coordinate sequence (the number of entries). This is determined by the
reference system.” Con criterios informáticos se llega a decir que las dos entradas en un registro de coordenadas
angulares se definen como bidimensionales. Se olvidan del elipsoide y de sus dimensiones. Tanto en ISO como en
INSPIRE se acrecienta el equívoco cuando se llama coordenadas horizontales a las dos coordenadas necesarias en
una posición no proyectada. Todo el mundo sabe que la tierra no es plana, aunque la podamos proyectar.
41
Sobre estas cuestiones ESPIAGO, J. La gestión de la información geoográfica, en BB
42
DATUM VERTICAL en teoría es la superficie utilizada con elevación o altitud cero. En la práctica es un
conjunto de datos de elevación cuya precisión permite utilizarlos como referencia de otros datos altimétricos
o batimétricos, aunque no describan una superficie equipotencial. Es decir, al igual que en los sistemas de
coordenadas que habíamos visto, en el datum vertical intervienen dos situaciones: una definición o
referencia teórica y un marco para su realización o materialización.
En la descripción del relieve conviene la utilización de un datum diferente al datum de la posición pues la
altitud no debe calcularse como una magnitud de la geometría regular. Como siguiente operación del
sistema de coordenadas entre las que nombrábamos al principio señalábamos la “determinación de un
datum geodésico para el cálculo de la altitud, llamado datum vertical. Puede ser distinto al datum utilizado
para las coordenadas (,) llamado, desacertadamente, datum horizontal”. Según ISO, hablamos de
sistemas de coordenadas simples o compuestos cuando utilizamos un solo datum o cuando utilizamos
43
también un segundo datum como datum vertical .
SISTEMA DE COORDENADAS COMPUESTO: Sistema de referencia de coordenadas que, en la descrip-

ción de la posición, utiliza dos sistemas independientes de coordenadas de referencia. Uno para el compo-
nente horizontal y otro para el componente vertical.
42
Además, en ISO 19111 se dice: “un datum para la ingeniería es el que no es geodésico ni vertical”. De momento,
podemos no considerarlo en geografía. Contraponiendo el calificativo, al elipsoide de ajuste lo llaman datum horizontal
con una nefasta visión geométrica plana. Las coordenadas geográficas no son planimétricas.
43
Traducimos la definición de [EN ISO 19111:2007, Geographic information-Spatial referencing by coordinates.”,
utilizada por INSPIRE D2.8. I.1-v3.2 apartado 2.3.
DETERMINACIÓN SATELITAL DEL CAMPO GRAVITATORIO. MODELOS DE GRAVE-

DAD EN LA SUPERFICIE TERRESTRE, EGM2008. ANOMALÍAS DE LA GRAVEDAD.
El antiguo geoide era de compleja determinación por el desconocimiento de las densidades internas. En
cambio, desde la década de 1960 se cuenta con la ayuda de satélites artificiales. En sus órbitas se
describen otras superficies con potencial de gravedad constante que disminuyen su valor a medida que nos
alejamos del centro. Un sistema de referencia satelital orbitando la tierra siempre lo hace respecto al centro
de gravedad terrestre. También sus superficies físicas equipotenciales son convexas y no se cortan nunca,
pudiendo utilizarse como superficies de nivel.
Con los satélites se persigue una ajustada determinación tridimensional y la medición del campo gravitatorio
que se compara con la gravedad normal del GRS80. Destacan los sistemas GNSS y la misión de los
satélites GRACE (Gravity Recovery And Climate Experiment), especialmente dedicada al campo de
44
gravedad. Han permitido que la precisión sea próxima al metro en el cálculo de la dimensión máxima de
W 0. También el resto de parámetros ha ido obteniendo una precisión mayor (a excepción de G).
Elaborado a partir de la fig. pág. 49 HEISKANEN y MORITZ.
SUPERFICIE DE NIVEL o SUPERFICIE EQUIPOTENCIAL (W): superficie que tiene un potencial constante
y que en todos sus puntos es perpendicular a la dirección de la gravedad. Expresa un nivel físico pues en
esas superficies el agua no se movería. Existe un número infinito y no se intersectan. La distancia entre dos
45
superficies es menor en los polos que en el ecuador .
Las superficies W, aunque llegan a ser analíticas, no son superficies físicas que se describan mediante una
ecuación regular o un modelo geométrico regular. Presentan desigualdades debidas a la diferente distri-
bución de masas en el interior terrestre (si la tierra fuera un cuerpo homogéneo y regular, las superficies W
se corresponderían con elipsoides de revolución equipotenciales), al achatamiento de los polos, a la
atmósfera, a la acción del sol, de los planetas y de la luna y a los efectos electromagnéticos. Pese a ello,
son superficies físicas básicas en la descripción terrestre que ayudan al entendimiento de la forma real de la
tierra sin incluir su topografía. También serán la principal descripción del tamaño y la orientación de la tierra
global sin que existan otras superficies no satelitales que puedan sustituirlas. Suele decirse que son esferoi-
des pues tampoco son elipsoides regulares. El análisis especializado de los datos orbitales expresa la gra-
vedad exterior de la tierra mediante series de armónicos esféricos en las que se eliminan los factores
distintos a los específicamente terrestres. Los coeficientes de los términos armónicos son función de las
distintas perturbaciones orbitales. Solamente unos pocos dependen de la forma de la tierra.
Su cálculo corresponde en la actualidad a unos escasos organismos que emplean las técnicas espaciales
que habíamos nombrado anteriormente entre las que destacan los satélites GRACE. Al igual que las nor-
44
“The first decade of the new millennium has been called “The Decade of Geopotentials” and has seen the launch of
three dedicated gravity field mapping missions: CHAMP [Reigber et al., 1996] launched in July 2000, GRACE [GRACE,
1998] launched in March 2002, and GOCE [European Space Agency, 1999] launched in March 2009.” PAVLIS, et al.
(2012), págs 1 y 2
45
Traducido de https://confluence.qps.nl/display/KBE/Geoid+Model+Glossary. Se añade la siguiente nota: “The mean
sea level is not an equipotential surface because oceans have a permanent topography that is caused by temperature,
salinity, currents, etc.”
mas internacionales ITRS cuyas determinaciones ITRF hemos adoptado en la georreferenciación, queda
fuera del alcance de la mayoría de instituciones públicas (autonómicas o estatales) y de los profesionales
tradicionales, aunque todos pueden acceder a sus resultados en direcciones electrónicas. Pese a todo, las
observaciones satelitales han de complementarse con los cálculos de la gravedad terrestre obtenidos con
gravímetros y con el estudio de la topografía o la altimetría obtenidas en los institutos tradicionales.
Anteriormente, el sistema americano WGS84 había elaborado el todavía denominado modelo de geoide
global llamado también WGS84 o, más específicamente, Earth Gravitational Model 1984 (EGM84) que
presentaba 200 km de resolución. Después, proporcionó en 1996 otro modelo conocido con las siglas
46
EGM96 con 100 km de resolución que fué actualizado en 2004 pero manteniendo su antiguo nombre. Se
destacaban anomalías gravimétricas en los océanos que tampoco se comportaban como superficies
enteramente normales (aunque fueran equipotenciales). Las ondulaciones se encontraban en el rango entre
−105 m hasta + 85 m.
47
Actualmente, el modelo gravimétrico global es el EGM2008, recomendado por el IERS , en el que se
alcanza nuevamente una precisión mayor mediante el uso de satélites específicos, medidas de las
anomalías de la gravedad en tierra y cálculos de altimetría satelital. Con una resolución próxima a 10 km,
48
llega a alcanzar 4,6 millones de términos en la expansión esférica . Los dos anteriores 32.757 y 130.317,
respectivamente. Se le calcula una precisión del orden ± 0.05 m (océanos) y ± 0.5 m (áreas continentales).
EGM2008 Geoid Undulations on a Global 2.5’ x 2.5’ Grid, with respect to the WGS 84 Ellipsoid. Unit is Meter.
Esta precisión aumenta en las regiones planas disminuyendo en las áreas de relieve montañoso. La
49
deflección de la vertical alcanzada en los USA y en Australia es de 1,1 a 1,3 arco segundos .
Además, se accede a una síntesis de esos armónicos esféricos con valores gravimétricos en mallas de 1 x
50
1 minuto o de 2,5 x 2,5 minutos de resolucion . En su descripción se utiliza una figura mediante armónicos
46
“spherical harmonic coefficients complete to degree and order 360”
47
“The official Earth Gravitational Model EGM2008 has been publicly released by the U.S. National Geospatial-
Intelligence Agency (NGA) EGM Development Team. This gravitational model is complete to spherical harmonic degree
and order 2159, and contains additional coefficients extending to degree 2190 and order 2159”. Los datos son “a global
set of area-mean free-air gravity anomalies… formed by merging terrestrial, altimetry-derived, and airborne gravity
data… supplemented with gravitational information implied by the topography… A single global model may satisfy the
requirements of a very wide range of applications.” PAVLIS et al. (2012) pág. 1. Se permite “full access to the model's
coefficients and other descriptive files with additional details”. La URL del EGM08 es http://earth-info.nga.mil/GandG/
wgs84/gravitymod/egm2008. Para los formatos SIG: EGM2008 GIS Data Page
48
“For many applications the complete series is unnecessarily complex and is truncated after a few (perhaps several
dozen) terms”.
49
Debido a esas anomalías de la gravedad en los océanos “their use was restricted to certain coastal areas, and to
areas where significant ocean surface variability makes altimetry derived gravity anomalies less reliable, while marine
data of verifiable quality exist”. PAVLIS et al. (2012), pág. 3.
esféricos pero esos datos se transforman también a un siste ma raster o en cuadrícula, como el de la
0
imagen anterior. Es decir, los datos se ofrecen interpolados en una rejilla de puntos equidistantes. En la ilus-
tración, vemos que este modelo global EGM2008 presenta ondulaciones en el rango -102 m (en el Océano
Índico) y +78 m (en el Pacífico, al norte de Australia). Su cartografía hipsométrica y batimétrica refleja la
distribución de la densidad en el manto. En las zonas de subducción, el geoide presenta una depresión de
15 a 20 m de profundidad en donde el trazado de sus ondulaciones se corresponde con el de las dorsales
oceánicas.
En ficheros cada 45’ y ondulaciones en una malla de 2,5’ se descargan formatos grid activando las
51
cuadrículas siguientes:
EGM2008 GIS Data: Global 2.5 Minute Geoid Undulations. Activando sobre la malla (45º X 45º) se descarga un fichero zip,
http://earth-info.nga.mil/GandG/wgs84/gravitymod/egm2008/egm08_gis.html
La precisión y la utilidad del EGM2008 pueden distinguirse en dos ámbitos distintos, según se trate de las
superficies oceánicas o las del área continental. En las oceánicas, se usa para el cálculo de la Dinamic
Ocean Topography (DOT). Se llega en ellas a un modelo gravitacional global de gran precisión y elevada
52
resolución . En las áreas continentales tiene una precisión de los datos que es menor. De todas formas, se
utiliza en la determinación de coordenadas astronómicas y en aquéllas altitudes que, utilizando GPS, no exi-
gen, teóricamente, tareas de nivelación.
GEOIDE GLOBAL: está formado mediante ondas del campo de gravedad de gran y de medio radio.
53
GEOIDE LOCAL : formado mediante ondas de gran, medio y corto radio.
50
Escrito en FORTRAN se conoce como hsynth_WGS84. f. Además, se usan dos programas de interpolación llamados
interp_1min.f e interp_2p5min.f, respectivamente. Existe un ejecutable windows hsynth_WGS84.exe y otras ayudas al
usuario. “Note that the harmonic synthesis software (provided)… applies a constant, zero-degree term of -41 cm to all
geoid undulations computed using EGM2008 with the height_anomaly-to-geoid_undulation correction model (also
provided).
Similarly, all pre-computed geoid undulations incorporate this constant zero-degree term. This term converts geoid
undulations that are intrinsically referenced to an ideal mean-earth ellipsoid into undulations that are referenced to WGS
84. The value of -41 cm derives from a mean-earth ellipsoid for which the estimated parameters in the Tide Free system
are: a=6378136.58 m and 1/f=298.257686”. Recordemos que los valores iniciales eran: a=6 378 137 m y 1/f= 298.257
223 563. Se publicó en PAVLIS, N. y otros, (2012) op. cit.
51
Es lo que llaman “gis data format” suponiendo esos autores que en esta herramienta solamente se consideran
formatos raster. Como dice el IGN, “Para la obtención de los modelos de geoide válidos para España se obtuvieron dos
ficheros de ondulación EGM2008 (uno para Península, Baleares, Ceuta y Melilla y otro para Canarias), con un paso de
malla de 1' x 1' con los siguientes límites: • Península, Baleares, Ceuta y Melilla: Longitud: 9º 30' W - 4º 30' E Latitud:
35º N - 44 º N • Canarias: Longitud: 18º 30' W - 13º W Latitud: 27º 30' N - 29º 30' N”. INSTITUTO GEOGRÁFICO NA-
CIONAL. SUBDIRECCIÓN GENERAL DE ASTRONOMÍA, GEODESIA Y GEOFÍSICA. CENTRO DE OBSERVACIO-
NES GEODÉSICAS: s. a.): El nuevo modelo de geoide para España EGM08 – REDNAP, pág. 12. Se descarga un .pdf
en ftp://ftp.geodesia.ign.es/.
52
PAVLIS, N. y otros (2012). The development and… op. cit. apéndice, págs. 1 y 2.
53
“For local engineering applications and cadastre-surveying, global geoids are not as accurate as needed. For such
applications, local geoid models are calculated. These can only be developed using local field measurements. They offer
centimetre accuracy over several hundred kilometres, with a high resolution. Local geoids are not suitable for height
comparison over large distances since they are based on different origins and reference heights (different equipotential
levels)”. ZAKATOV, P. (1997), pág. 14.
La forma usual para mostrar el campo gravitatorio del EGM08 es mediante ondulaciones sobre el elipsoide
por lo que se sigue llamando geoide. Con los mismos datos pero comparándolos con los valores de la
gravedad normal obtenemos la diferencia entre la gravedad convencional y la realmente obtenida. Reciben
el nombre de anomalías de la gravedad en la superficie terrestre (no en el interior de la corteza) y también
54
free air gravity anomalies .
55
En rojo, áreas en donde los movimientos de las placas se asocian con terremotos y volcanes
Las ondulaciones se registran en metros, las anomalías en el intervalo + 100 mGals. Es destacable la cone-
xión entre la topografía, especialmente la marina, y las anomalías en la distribución regular de la gravedad.
Como dice el pie de la imagen anterior, se describen numerosos fenómenos geofísicos como las fosas y
dorsales oceánicas, las zonas de subducción y de fracturación o las cadenas submarinas. Se destacan las
elevadas anomalías que recorren buena parte del litoral del Pacífico y del Caribe. También en la costa del
Mediterráneo y del Índico asiático. Por ello, en cartografía, estos últimos datos debieran utilizarse en las tra-
dicionales altitudes al nivel del mar. En general, las anomalías han de tenerse en cuenta cuando indicamos
la simplista distribución equipotencial normal. Destacamos el carácter lineal de las áreas de tensión o máxi-
mo esfuerzo, de color rojo en la imagen inferior. Es inmediata su conexión con las anomalías de valores
máximos.
54
http://earth-info.nga.mil/GandG/wgs84/gravitymod/egm2008/anomalies_dov.html. Se acompañan de los datos sobre
deflección de la vertical (DOV).
55
El subtítulo original de la imagen es “Áreas Geologic, geodetic, and earthquake data help determine the zones on
Earth where the crustal motion diverges from rigid plate motion. Areas of high strain, in red, experience increased
earthquakes and volcanoes. SOURCE: Kreemer et al., 2000.”
TIPOLOGÍA DE ALTITUDES. ANOMALÍAS DE LA ALTITUD. OTRAS FIGURAS DE LA

TIERRA (TELUROIDE. CUASIGEOIDE).
Dependiendo de la superficie de referencia que adoptemos y del procedimiento utilizado, tendremos diferen-
tes tipos de altitud.
ALTURA ELIPSOIDAL (h): es la distancia, medida de forma geométrica y ortogonal, entre un punto P y el
elipsoide que es utilizado como datum vertical. Dicho de otra forma, es la elevación sobre el elipsoide de
56
referencia. A veces aparece nombrada como ALTURA GEOMÉTRICA o ALTURA GEODÉSICA .
En esta aproximación, la altura del punto P de la superficie terrestre se mediría geométricamente a través
de la perpendicular al elipsoide en el que se determina el punto P’. Este punto P’ tiene altitud = 0 pues está
situado sobre el elipsoide que podemos utilizar como referencia geométrica teórica. Los puntos P y P’ tienen
las mismas coordenadas  y  pero distinta altura. Tienen la misma posición geográfica aunque distinta
posición en la vertical.
Es un término equívoco pues esa distancia al elipsoide no indica, necesariamente, la altitud de la superficie
terrestre dado que el elipsoide no se ajusta de igual manera en todos sus puntos ni es una superficie de
nivel. No tiene en cuenta el campo de gravedad terrestre y, por lo tanto, esta altitud no define la circulación
57
del agua lo que ha de considerarse en morfología . Esa distancia h es una necesidad geométrica derivada
de la relación de las coordenadas tridimensionales X, Y, Z con las coordenadas elipsoidales   que
necesitan un valorh para la descripción de los puntos que no se encuentran sobre el elipsoide. En la
actualidad, es un cálculo geométrico a partir de las mediciones GNSS, no es un cálculo físico. Para la
obtención de valores alltitudinales de índole física son insuficientes las figuras geométricas.
ALTURA DE LA SUPERFICIE DEL MAR (Sea Surface Height, SSH): Es la altura elipsoidal en
la superficie oceánica. Puede observarse a partir de altimetría radar de satélites (Topex/
Poseidon, ERS-1, Jason, etc.).
ALTURA ORTOMÉTRICA (H): es la distancia de un punto P y el datum vertical medida a lo largo de la línea
de fuerza que es perpendicular a las diferentes superficies equipotenciales. Si conocemos las ondulaciones
del geoide, llamamos altura ortométrica (H) a la diferencia entre la altura elipsoidal (h) y la ondulación (N).
H = h - N.
Esa distancia de P a unos datos o a una figura de referencia a la que seguimos llamando geoide la
indicamos con H y es la que también acostumbramos llamar altitud sobre el nivel del mar.
Es una definición teórica. En la práctica, para un punto P, la altura elipsoidal ( hp) puede obtenerse a partir
0
de observaciones GNSS, la altura ortométrica (Hp) es observable mediante operaciones de nivelación,

corregidas por observaciones de gravimetría, y la altura del geoide u ondulación (Np) se calcula en los
modelos de gravedad terrestre (EGM). Esa definición teórica es una aproximación ya que no considera
diferencias entre las perpendiculares geométricas o físicas, ni consideran los diferentes factores de escala
según los tipos de observaciones. De todas formas, para propósitos cartográficos (geográficos o no) el error
producido por esta aproximación puede ser normalmente ignorado pues es de pocos milímetros.
TOPOGRAFÍA DE LA SUPERFICIE DEL MAR (Sea Surface Topography, SST): Separación
entre la superficie del mar y el geoide. En los océanos es equivalente a la altura ortométrica.
SST = SSH - N.
56
En ISO 19111:2007, 4.19 “ellipsoidal height: geodetic height h. Distance of a point from the ellipsoid (4.1.12)
measured along the perpendicular from the ellipsoid to this point, positive if upwards or outside of the ellipsoid. Note 1 to
entry: Only used as part of a three-dimensional ellipsoidal coordinate system (4.1.13) and never on its own”.
57
“Gravity defines the vertical (up/down) dimension in the real (non-ellipsoidal) world” en http://www. hydrometronics.
com/downloads/What's Up.pdf.
NÚMERO GEOPOTENCIAL (C): Para una localización dada, es la diferencia entre su superficie
equipotencial (W i) y la superficie equipotencial utilizada como referencia (W 0) medida a lo largo de la línea
de la plomada.
DIFERENCIA EN EL NÚMERO GEOPOTENCIAL (ΔCij): Para dos localizaciones distintas en la superficie
de la tierra (i, j), es la diferencia entre sus superficies equipotenciales W j − W i. Puede calcularse por
nivelación.
ALTURA NORMAL (H*): es la diferencia del potencial de un punto respecto del potencial al nivel del mar.
Tiene en cuenta la gravedad normal medida a lo largo de la línea de plomada.
La diferencia entre la altura normal y la altura ortométrica suele ser significativa en altas elevaciones.
58
ANOMALÍAS DE LA ALTITUD (: diferencia entre las alturas normales y la altura sobre el elipsoide .
Permiten la conversión de las alturas elipsoidales (h) en alturas normales (H*):
h = H* + ,
siendola anomalía de la altitud .
59
TELUROIDE. Las distancias al elipsoide medidas como altitudes normales H* y no a través de cálculos geo-
métricos describen una figura llamada teluroide cuya ecuación es, precisamente : h = H* +  .
60
A diferencia del geoide, el teluroide es una superficie “observable” aunque no es una superficie equipo-
tencial. Las posiciones de sus puntos sobre el terreno se miden usando mediciones de la gravedad, nivela-
ciones y observaciones astronómicas. Los resultados se comparan con la altitud elipsoidal h para cualquier
61
latitud proporcionando las llamadas Molodensky free-air gravity anomalies.
Sobre los océanos se pensaba que siempre = N pero en las áreas continentales “la anomalía de la altitud
(es, en general, mayor que la correspondiente ondulación del geoide N”.
Se supone que sobre los océanos:  - N = 0
Sobre los continentes:  - N = H - H*
La altura ortométrica (H) que resultaba de la ondulación o altura del geoide (+ N, -N) que era la separación
entre el geoide y el elipsoide, se calcula de forma separada a las anteriores.
“Observemos… que también podemos dibujar las anomalías de la altitud  por encima del elipsoide .De
esta forma, obtenemos una superficie que es idéntica al geoide sobre los océanos porque allí = N y es muy
58
“The distance between points P and Q, measured along the ellipsoidal normal is called height anomaly, and
corresponds to the geoid undulation that is used in the conventional approach”. “Over the ocean, height anomalies are
virtually identical to geoid undulations; over land their difference is a function of the Bouguer anomaly and the elevation “,
HEISKANEN, W y MORITZ, H, (1967), pág. 292 y sección 8–13].
59
“En Europa existen tres tipos de altitudes (…) normales, ortométricas y normal-ortométricas. Ejemplos de ortométricas
serían las de Bélgica, Dinamarca, Finlandia, Italia, Suiza o la Península Ibérica. Hoy en día, las altitudes normales son
usadas en Francia, Alemania, Suecia y la mayoría de países de Europa Oriental. En Noruega, Austria y en los países de
la antigua Yugoslavia se usan la normal-ortométricas”. Además, hemos de considerar las alturas elipsoidales y otras
magnitudes que describimos en apéndice.
60
En Geoid Model Glossary (Jos Verheijen, 2015):”The telluroid was proposed by Molodenskii to avoid the complex
determination of the topographical density and vertical gradient of gravity, which are necessary components in geoid
modelling”. Se determina por nivelación, combinándola con medidas de la gravedad. HEISKANEN, W. y MORITZ, H.
(1985), op. cit. pág. 292. Otra definición es “Is a surface whose normal potential U at every point Q is equal to the actual
gravity potential W at the corresponding surface point P, with the points P and Q being situated on the same line that is
normal to the ellipsoidal [HEISKANEN y MORITZ, 1967, p. 292].
61
HEISKANEN, W y MORITZ, H., (1967), sección 8–3 citados por PAVLIS et al (2012) pág. 4.
próxima al geoide por las otras partes. Esta superficie ha sido llamada cuasigeoide por Molodensky”…”el
62
cuasigeoide no es una superficie de nivel, no tiene significado físico alguno” ,”pero proporciona una
63
visualización conveniente” y “refleja la topografía” .
Heyskanen y Moritz (1985), fig pág. 328.
CUASIGEOIDE: superficie paralela al teluroide obtenida sumando las anomalías de la gravedad al

elipsoide. Es una función de la altura normal y no tiene un significado físico.
En topografías suaves o discretas, la diferencia entre el geoide y el cuasigeoide es muy pequeña. En áreas
montañosas llega a varios decímetros o, incluso, al metro. El geoide es el datum de la altitud ortométrica y
el cuasigeoide de la altura normal.
Considerando esas dos superficies (teluroide y cuasigeoide), volvemos a nuevas definiciones de los mismos
términos anteriores:
ANOMALÍA DE LA ALTITUD (ζ): es la separación entre el teluroide y la superficie del terreno o
la superficie oceanográfica.
ANOMALÍA DE LA ALTITUD (ζ): es la separación entre el elipsoide y el cuasigeoide.
ALTURA NORMAL (H*): es la elevación del teluroide sobre el elipsoide.
ALTURA NORMAL (H*): es la elevación sobre el cuasigeoide.
62
HEISKANEN, W y MORITZ, H, (1985), págs. 293-294.
63
Los entrecomillados en HEISKANEN, W. y MORITZ, H. (1979), págs. 328, 294 y 327, respectivamente.
MODIFICACIONES REGIONALES DEL EGM2008. DATUM VERTICAL EUROPEO

(EVRS, European Vertical Reference System). EGM08- REDNAP
Para las coordenadas de la posición geodésica () adoptamos el elipsoide pero no existe un acuerdo
internacional sobre una figura o unos datos que sirvan de referencia global para el cálculo de la altitud. No
64
existe un único datum vertical reconocido internacionalmente . Por ello, los valores de altitud no son direc-
tamente comparables. Se utilizan las distintas superficies que representan la forma y tamaño de la tierra de
cuya descripción hemos hablado. Se controlan con acuerdos regionales que, para nosotros, serán los de la
Unión Europea y, más concretamente, los requerimientos de España. En ambos casos se parte del modelo
gravimétrico global elaborado en USA conocido como EGM2008 pero modificándolo con los datos obtenidos
para una red de puntos calculados con gravímetros de precisión, con nivelaciones y con los cálculos
65
GNSS . Ya hemos dicho, que, también en la práctica, se llega a la definición de un conjunto de datos de
referencia para los que se emplea el nombre geoide W 0 y que utilizamos como datum vertical regional o
66
local . Además, el antiguo nombre de geoide sigue siendo útil para la comprensión de otras figuras que se
han usado en el cálculo altimétrico (teluroide y cuasigeoide).
Hasta la geodesia satelital se utilizaban redes de nivelación geométrica estatales a partir de la observación
de uno o de varios puntos establecidos, preferentemente, en la parte marina. La superficie equipotencial
usada como datum vertical se observaba con las variaciones de las mareas en un período de tiempo. La
observación podía ser directa o utilizando el nivel medio0del mar. Posteriormente, en esa nivelación se
añadían algunos cálculos astronómicos y medidas de la gravedad también en el interior continental. En las
áreas continentales se decía que eran puntos del geoide (alturas sobre el nivel del mar). Como decimos,
hasta fechas recientes eran infraestructuras aplicables en el ámbito de los estados que servían para el
control de la posición y de la altitud. La determinación del atributo de la altitud para cada punto , no es
67
una tarea sencilla ni individual. Para esas situaciones locales se definía el datum de mareas como datum
vertical.
DATUM DE MAREAS (tidal datum): Niveles de marea que se utilizan en la determinación de superficies de
referencia de la batimetría marina y0como datos base en la determinación de la altitud terrestre. En el
68
ámbito emergido, eran la base de las tareas de nivelación. Llamados también planos de marea .
Su cálculo exige mediciones continuadas (que suelen estar normalizadas) en estaciones de marea espe-
cíficas. Tradicionalmente, se emplea el nivel medio del mar (NMM o MSL por sus siglas en inglés) pero
existen otras mediciones de las fases de marea (Mean Low Water, MLW; Mean High Water, MHW).
NIVEL MEDIO DEL MAR: (NMM, mean sea level, MSL) altura promedio de la superficie del mar en una
estación mareográfica durante un periodo de 19 años. Determinada, generalmente, mediante lecturas
69
horarias a partir de un nivel de referencia fijo predeterminado . Se intenta la obtención de datos promedio
en las variaciones de las mareas motivadas por las fuerzas gravitatorias del sol y la luna. De todos modos,
ese promedio se presenta con picos y valles respecto al geoide aunque se filtren los efectos de marea de
corto plazo.
64
La herramienta EPSG, para la transformación de coordenadas, (ver apéndice) incluye unos 130 datums verticales
usados en el mundo. Algunos países tienen más de uno.
65
En apéndice estudiamos las distintas redes geodésicas.
66
La definición ISO es la del cuadro:
ISO19111:2007(E)
vertical datum datum describing the relation of gravity-related heights or depths to the Earth
ISO19116:2004(E)
67
“Tidal datums are used as references to measure local water levels and should not be extended into areas having
differing oceanographic characteristics without substantiating measurements”. También, “tidal datums are also the basis
for establishing privately owned land, state owned land, territorial sea, exclusive economic zone, and high seas
boundaries”. Lógicamente: “In order that they may be recovered when needed, such datums are referenced to fixed
points known as bench marks”. Todo en Center for Operational Oceanographic Products and Services https://tidesand
currents.noaa. gov/datum_options.html. En http://www.ign.es/ign /hwide/redMareografica.do se insiste en el carácter
local del datum de mareas: “Are local references and should not be extended into areas with different oceanographic
characteristics” y “Are referenced to fixed and stable points on land (Bench Marks)”.
68
USIGS GLOSSARY (NUG) (1997). Traducimos “Specific tide levels which are used as surfaces of reference for depth
measurements in the sea and as a base for the determination of elevation on land. Many different datums have been
used, particularly for leveling operations. Also called tidal datum plane”.
69
SP-32 (OHI en su 5ª edición, 1994).
Como decimos, el NMM ha sido un dato utilizado como referencia principal para las elevaciones altimétricas
70
antes del uso de las tecnologías satelitales y del empleo de gravímetros de precisión . Se llegaba a la
definición de la superficie potencial en un punto que por acuerdo institucional se adoptaba para la definición
de todos los demás puntos de un territorio dado a los que se asignaba el valor W 0. Pero, ahora se sabe, de
esta manera no se permite la descripción de una superficie equipotencial sobre los océanos (no coincide
exactamente con una superficie de nivel). La superficie del océano presenta pequeñas elevaciones o valles
71
en el rango de -2,0 m a + 2,0 m respecto al geoide . Aunque no se permita la descripción de una superficie
real, sigue utilizándose como una superficie teórica de fácil imaginación que, convencionalmente, será al
cero de elevación (o de nivel), ya que la superficie del mar está disponible desde casi todas partes
DATUM VERTICAL EUROPEO – EVRS (European Vertical Reference System): sistema vertical teórico,
relacionado con un nivel cero de las mareas que emplea las diferencias de gravedad para las altitudes.
Ya desde 1994 la EUREF perseguía la unificación de las redes de nivelación estatales europeas adoptando
un datum vertical único de resolución decimétrica. Se creó una red europea llamada UELN (United
European Levelling Network) en cuyos nodos se calculan los números geopotenciales y las alturas
normales. Para la conversión de valores geopotenciales en valores normales se utiliza el campo de
72
gravedad del GRS80 y las coordenadas geodésicas se expresan en valores ETRS89 .
En el año 2000 se concretó o materializó el Marco Europeo de Referencia Vertical (EVRF2000) formado por
GPS estacionarios, usados para la obtención de alturas elipsoidales y utilizando también esos GPS para la
nivelación, añadiéndoles observaciones del nivel del mar con mareógrafos y midiendo las variaciones de la
73
gravedad. El sistema mareal cero, de acuerdo con las resoluciones del IAG , era el llamado Normaal
Amsterdams Peil (NAP). Se estableció la relación entre ese datum NAP y los pertenecientes a cada país
europeo. Posteriormente, se añadieron las redes de nivelación de nuevos países y se calculó una nueva
realización que es la que se conoce como EVRF2007 que incluye la velocidad de desplazamiento vertical.
Utiliza un total de trece puntos de referencia principal con mareógrafos en el Báltico, Mar del Norte,
Mediterráneo, Mar Negro y Océano Atlántico con diferencias de varios decímetros.
Resumiéndolo, el sistema europeo cumple con las cuatro convenciones siguientes:
 El datum vertical es una superficie equipotencial en la que el campo de gravedad es constante:
W 0= W 0E= constante
Esa superficie equipotencial pasa por el mareógrafo de Amsterdam.
 La unidad de longitud es el metro. La de tiempo, el segundo. Se obtienen mediante el correspondiente
modelo de la relatividad.
 Los componentes de la altura son las diferencias ΔW P entre el potencial W P del campo de gravedad de
la tierra para el punto P y el potencial W 0E del EVRS en el nivel convencional cero. La diferencia potencial
−ΔW P se nombra también como el número geopotencial CP. Así: −ΔW P = CP =W 0E −W P
 Los números geopotenciales son equivalentes a las altitudes normales, con la condición de que el
campo de gravedad quede especificado.
 El EVRS es un sistema mareal cero.
70
“Conventionally, the mean sea level has been used as the reference surface for topographic elevations for
generations. Regional, national and continental vertical datums were realized using geodetic levelling observations that
were constrained to mean sea level as determined by tide gauge measurements”, http://www.nrcan.gc.ca.
71
“The mean ocean surface has slight hills and valleys similar to land topography, but much smoother. Globally, these
hills and valleys range from -2.0 m to +2.0 m with respect to the geoid… This separation between the geoid and the sea
surface is the Sea Surface Topography (SST) or Dynamic Ocean Topography (DOT).” Llegan a la definición
(normalizada) de Sea Surface Topography (SST): It is the separation between the geoid and ocean surface. SST can
be determined from Sea Surface Height (SSH) measured from satellite radar altimetry and a geoid height (N): SST =
SSH - N. SST is positive if the ocean surface is above the geoid. It is equivalent to the orthometric height of the ocean
surface. (Unit: m), http://www.nrcan.gc.ca. Las normas también hablan de Sea Surface Height (SSH).
72
Es, por lo tanto, un sistema de coordenadas compuesto conocido como ETRS89-EVRS. “These systems [ETRS89 y
EVRS] are the basis for geo-referencing in Europe both on a national as on an international level and they have been
recommended for adoption by the European Union, EuroControl and EuroGeographics, the consortium of the European
National Mapping and Cadastre Agencies”, pág, 1 de BRUYNINX, C y ROOSBEEK, F. (2006).
73
AIS-AIMS G/6-SN/5, Attachment A. Las siglas de autor corresponden a Aeronautical Information Services (AIS) y
Aeronautical Information Management (AIM). Posteriormente, la red de nivelación de Holanda no consideró el NAP y no
se puede relacionar un ajuste a este punto. Para algunos “esto ofrece una oportunidad única de pasar del concepto de
un único punto para la definición del datum a una especificación del mismo multi-estación”, en www.upv.es/unigeo
74
EVRF2007. Marco del datum vertical para Europa aprobado desde 2008 por la IAG-EUREF . Es un
75
sistema de referencia vertical cinemático y gravimétrico .
El European Terrestrial Reference System 89 (ETRS89) es la referencia geocéntrica con precisión
centimétrica para las coordenadas geodésicas en Europa. El European Vertical Reference System (EVRS)
cumple la misma función para la altitud.
EGM08- REDNAP
En España, “el Sistema de Referencia Altimétrico tomará como referencia el nivel medio del mar
Mediterráneo en Alicante para la península y las respectivas referencias mareográficas locales para el resto
de los territorios. En todos los casos, el sistema altimétrico queda materializado por la nueva Red Española
76
de Nivelación de Alta Precisión (REDNAP) con altitudes ortométricas Helmert (H)" . Corresponde al IGN “la
planificación y gestión de las redes geodésicas terrestres de ámbito nacional, de la red de nivelación de alta
precisión y de la red de mareógrafos, así como el desarrollo de aplicaciones sobre sistemas de navegación
y posicionamiento”. Las medidas de la gravedad absoluta exigen un instrumental al alcance de un número
77
limitado de instituciones que conforman una red internacional de observatorios puntuales . Se densifica esa
red mediante la utilización de otro instrumental que mide las diferencias de gravedad entre puntos
(mediciones relativas).
Ese Decreto fue anterior a 2008 aunque, para el territorio del Instituto Geográfico Nacional (IGN), “la calidad
del modelo gravimétrico EGM2008 era tal que las diferencias en la mayor parte de los casos… no
78
superaban el decímetro” . De todos modos, se corrigió mediante su combinación con los datos de la
nivelación geométrica (REDNAP), con observaciones gravimétricas y con datos GNSS. El modelo
resultante, conocido como EGM08-REDNAP tiene diferencias centrimétricas. Sobre los datos utilizados, se
79
publica el siguiente cartograma :
74
Es una determinación práctica, un marco (EVRF) del EVRS (S de System). Recomendado en Europa, “to express
gravity-related heights” (Commission Regulation (EU) No 1089/2010). Cuando se utliliza el ETRF89 como sistema
elipsoidal estamos ante un “compound coordinate reference system ETRS89/EVRS” como dice INSPIRE, Specification
on Coordinate reference systems, págs. 8 y 3. La página EVRS en https://evrs.bkg.bund.de/Subsites/EVRS/EN/Home/
home.html. “The European Vertical Reference System (EVRS) is an example of a system that has been built to reflect
the globalisation of Geographic Information System (GIS) applications and the need for continental-wide, consistent
height information… is realised in the European Vertical Reference Frame (EVRF) by the geopotential numbers and
Normal Heights of nodal points of the United European Levelling Network … in relation to the Normaal Amsterdams Peils
(NAP). The geopotential number at NAP is zero”. Vertical Reference… pág. 3. “Guarantee of a 1 cm accuracy level for
datum and network realization” en SACHER, M. (2015), pág. 13. Sobre el ámbito de aplización: for the vertical
component on land, within its geographical scope (continental Europe). Outside the geographical scope of EVRS (e.g.
overseas territories), other vertical reference systems related to the Earth gravity field shall be used, which shall be
defined and documented according to EN ISO 19111 and ISO 19127 by the Member State concerned. The Earth
Gravitational Model (EGM Version 2008) is recommended in that cases” INSPIRE D2.8.II.1_v3.0, pág VI.
75
“Classical geodetic positioning deals with stationary points (objects). In recent times, however, geodetic positioning
has found its role also in positioning moving objects, such as ships, aircraft and cars. This application became known as
kinematic positioning, and it is understood as being the real-time positioning part of navigation” VANÍCEK, Petr: An
Online Tutorial in Geodesy, http://einstein.gge.unb.ca.
76
Según el Real Decreto 1071/2007 (Artículo 4). “En la Península, las altitudes, referidas al nivel medio del Mar
Mediterráneo en Alicante, tienen el código EPSG 5782. En Baleares, las altitudes están referidas al nivel medio local del
Mar Mediterráneo en las islas. Código EPSG: 5782. En Ceuta y Melilla, las altitudes están referidas al nivel medio del
Mar Mediterráneo. Código EPSG: 5782. En Canarias, las altitudes están referidas al nivel medio del Océano Atlántico en
las islas”, (procedente de los metadatos del MDT 200 del IGN).
77
“The stablishment of these global gravity networks has been coordinated by the IUGG and the current network is
known as the International Gravity Standardisation Network 1971 (ISGN 1971)”, Vanicek y Krakiwsky (1986), p. 535. En
España, como dice la página del IGN, “el servicio de Gravimetría de la Subdirección de Astronomía, Geofísica y
Aplicaciones Espaciales del Instituto Geográfico Nacional tiene como misión el establecimiento de la infraestructura
básica de Gravimetría, y entre otras, el establecimiento de las llamadas redes gravimétricas absolutas de orden cero y
primer orden.
78
IGN (2010): El nuevo modelo de geoide para España EGM08–REDNAP, pág.9.
79
En ftp://ftp.geodesia.ign.es/documentos/EsquemaRGeodesicas_201501.jpg se descarga un mapa sobre el asunto.
Existe un visor sobre las redes para el posicionamiento en ftp://ftp.geodesia.ign.es/documentos/EsquemaRGeodesicas_
201501.jpg

Las mediciones son independientes para cada una de las 25.000 señales en una red que se extiende sobre
unos 18.000 km. Cada punto se concreta en una marca material (monumento) y los valores se gestionan en
una base de datos de acceso informático. Actualmente, las coordenadas geodésicas se han transformado a
ETRS89 desde los sistemas antiguos. El modelo raster resultante simplifica el modelo EGM2008 original en
el que se registraban algunos detalles geológicos de interés (las unidades son metros).
MODELO FINAL EGM08- REDNAP EGM08

Espiago 2017 Cartografía para La Geografía - Georreferenciación

Cargado por

Copyright:

Formatos disponibles

Espiago 2017 Cartografía para La Geografía - Georreferenciación

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Espiago 2017 Cartografía para La Geografía - Georreferenciación

Cargado por

Copyright:

Formatos disponibles

CARTOGRAFÍA PARA LA GEOGRAFÍA

Sistema Internacional de Referencia Terrestre (ITRS - ITRF).

2. GEORREFERENCIACIÓN Y REGISTRO DE LA ALTIMETRÍA.

Altitud geográfica y datum vertical.

Tipología de altitudes. anomalías de la altitud.

Modificaciones regionales del EGM2008.

1. GEORREFERENCIACIÓN Y SISTEMAS DE COORDENADAS.

Coordenadas cartesianas según ISO 19111

2) su materialización (determinación práctica) para un momento temporal y

Georreferenciación y sistemas de coordenadas. 3

Las últimas fases del proceso son:

Meridiano cero Meridiano local

Sistema Internacional de Referencia Terrestre (ITRS - ITRF).

B) Determinación práctica del sistema (ITRF).

El movimiento de los polos terrestres presenta un componente anual de trayectoria elíptica y un

C) Control del sistema ITRF.

Georreferenciación y sistemas de coordenadas. 10

Elección de una figura de ajuste.

Fig. pág. 50 Maling

Georreferenciación y sistemas de coordenadas. 11

Cálculo de distancias y direcciones en la tierra esférica.

Elaborado a partir de Alves, fig. pág. 62.

Georreferenciación y sistemas de coordenadas. 14

Elección del radio de la tierra esférica.

2. Media aritmética de a y b: 6.367.650.

O que, generalmente, suele referirse mediante su recíproca 1/f

Georreferenciación y sistemas de coordenadas. 16

Referencias geográficas sobre el elipsoide.

La longitud en el elipsoide se define de igual manera que en la esfera. En cambio, en la definición de la

Elaborado a partir de la figura de la pág. 66 Maling.

Cálculo de distancias y direcciones en la tierra elipsoidal.

Siendo Rm y Rp los radios de curvatura del meridiano y el radio de curvatura transverso

Georreferenciación y sistemas de coordenadas. 18

Evolución de las determinaciones de las superficies de referencia20.

Longitud de los semiejes

1687 Newton l/231 -

Elipsoide equipotencial GRS80. Determinación de elipsoides de ajuste.

al meridiano de longitud cero del BIH.

Sistemas ETRS89 y WGS84.

O relación entre ejes b / a = 0, 996647189318816362;

cuadrante de un meridiano = 10.001.965,7293 m.

2. GEORREFERENCIACIÓN Y REGISTRO DE LA ALTIMETRÍA.

FIGURAS DE REFERENCIA. GEOIDE.

Ondulación del geoide

Ondulación del geoide

ELIPSOIDE EQUIPOTENCIAL MEDIO TERRESTRE O ELIPSOIDE DE NIVEL.

Semieje mayor a = 6 378 137, 00 metros

Velocidad angular  = 7 292 115,010 -11

CÁLCULO DE LA GRAVEDAD NORMAL.

ALTITUD GEOGRÁFICA Y DATUM VERTICAL.

SISTEMA DE COORDENADAS COMPUESTO: Sistema de referencia de coordenadas que, en la descrip-

DETERMINACIÓN SATELITAL DEL CAMPO GRAVITATORIO. MODELOS DE GRAVE-

Elaborado a partir de la fig. pág. 49 HEISKANEN y MORITZ.

TIPOLOGÍA DE ALTITUDES. ANOMALÍAS DE LA ALTITUD. OTRAS FIGURAS DE LA

de observaciones GNSS, la altura ortométrica (Hp) es observable mediante operaciones de nivelación,

Heyskanen y Moritz (1985), fig pág. 328.

CUASIGEOIDE: superficie paralela al teluroide obtenida sumando las anomalías de la gravedad al

MODIFICACIONES REGIONALES DEL EGM2008. DATUM VERTICAL EUROPEO

Georreferenciación y registro de la altimetría. Figuras de referencia. 38