Fundamentos Georref v1
Fundamentos Georref v1
Fundamentos Georref v1
es
Julio - 2017
Introducción.
Sistemas de coordenadas utilizadas en la georreferenciación
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 formas, 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:
Coordenadas cartesianas:
ISO 19111
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 localización de
los satélites artificiales utilizados en geografía. Es decir, las coordenadas cartesianas tridimensionales 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 navegació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
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.
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 coordenadas
geográficas angulares (φ,λ), todo punto puede referirse con dos sistemas diferentes. Uno lo sitúa en el espacio
el otro, respecto a la tierra1.
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).
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.
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 se utiliza el
conocido como International Terrestrial Reference System (ITRS) 2. Es un sistema ideal o convencional con el
que adoptamos unas propiedades teóricas, matemáticas y físicas.
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.
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.
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.
La concreción práctica (la realización o materialización) de este sistema teórico recibe el nombre de International
Terrestrial Reference Frame (en siglas ITRF). Ha conocido distintas versiones 3 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 determinación
práctica correspondía al Bureau International de l´Heure (BIH)4. 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 Systems
Service (IERS)5 con precisión de 0.005“ sobre el CTP. Este polo es conocido como IERS Reference Pole (IRP)
siendo el utilizado en la georreferenciación6.
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 el origen, la escala y la orientació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ón del
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).
6 http://www.iers.org/
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 componentes
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.
4. Como punto de la superficie instituido convencionalmente para la definición del plano meridiano origen (0° 0'
0")7 se elige el Observatorio de Greenwich8.
Se observan permanentemente las coordenadas de estaciones terrestres 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
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), 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.
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 rígidas (high-quality
sites) o sobre áreas de deformación10. 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 estructura profunda de la tierra y a la aplicación de
modelos atmosféricos.
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.
GGOS: An observing system of layered infraestructure (based on the Fig. 9.2. and Fig.9.9. in chapter 9
(Rothacher et al., 2009) of PLAG and PEARLDMAN (2009)(c) Springer-Verlag Berlin Heidelberg 2009)
11 http://cddis.gsfc.nasa.gov/Techniques/GNSS/GNSS_Overview.html
12Para 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, R (2010).
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 glaciarismo,
los embalses de agua, el agua atmosférica, el movimiento terrestre y las variaciones del campo de gravedad.
Si a las técnicas espaciales añadimos las antiguas técnicas “terrestres” como la nivelación, la fotogrametría, la
gravimetría o las estaciones que miden el nivel del mar, se estudia entonces el sistema en sus implicaciones
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 como la migración de los polos y el estudio de la rotación terrestre y la gravedad son los
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).
La más sencilla entre las aproximaciones cartográficas a la figura de la tierra es la de una esfera que presenta
irregularidades del terreno, aunque, a determinadas escalas, estas irregularidades no son de gran magnitud.
Con un radio terrestre de aproximadamente 6.37114 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 es 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, nombramos polos a los dos puntos en su intersección con la superficie de la esfera. El plano
perpendicular al eje y que pasa por el centro define un 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.
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 de λ, 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 o rejilla 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 (δ).
Coordenadas geográficas
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
siendo z en radianes
15La 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.
Distancia entre dos puntos situados en el mismo meridiano: d = R ∙ Δφ, siendo Δφ en 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φ ∙ Δλ , siendo Δλ 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.
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.
ACIMUT (Z): ángulo esférico formado por el meridiano y un círculo máximo en un punto.
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.
La convergencia de meridianos es función de la latitud y puede calcularse como convergencia angular (θ) 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.
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.
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 Internacional 17
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 el cuadro de la página siguiente al lado
de los de Maling. Los tres primeros procedimientos se basan en la elipse.
17 Sobre el Elipsode Internacional y otros elipsoides anteriores al GRS80 ver la tabla de principales determinaciones de la figura de la
tierra en este mismo capítulo.
Podemos utilizar para el cálculo del radio terrestre la media de los radios calculados para las latitudes de 15º, 45º o 60º:
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:
Si calculamos el volumen del elipsoide y hallamos el radio de una esfera que tenga un volumen equivalente:
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:
Una esfera rectificada es aquella cuyos meridianos tienen la misma longitud que los meridianos del elipsoide:
ELIPSOIDE: figura matemática generada por la rotación de una elipse sobre uno de sus ejes. El elipsoide 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 aproximació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 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.
Geometría de la elipse
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
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):
e = OF1/a,
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.
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.
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
LONGITUD GEODÉSICA (λ): ángulo formado por el plano meridiano origen con el plano meridiano de un punto
dado, tomado como positivo hacia el este
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
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 utilizaron una esfera tangente al área del cálculo que
suele denominarse esfera local. La informática ha facilitado estos cálculos.
Distancia d entre dos puntos situados en el mismo paralelo (de radio R0) y con una diferencia de longitud de Δλ:
d = Δλ ∙ R0
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):
dg =ʃ[R2m+R2pcos2φ (dλ/d λ)2]1/2d φ
Siendo Rm y Rp los radios de curvatura del meridiano y el radio de curvatura transverso.
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.
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 se han utilizado diversos elipsoides 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.
a, 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 (J 2), 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 diferencias de 300 m.
21El 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élites 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 = 3.986.005
x 108 m3 s-2, factor de la forma dinámica de la tierra, excluyendo la deformación permanente de las mareas: J2 = 108 263 x
10-8, velocidad angular de la tierra: ω = 7292 115 x 10-11 rad s-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 al meridiano de longitud cero
del BIH.
De todas formas, a grandes escalas como en la ingeniería y los catastros, deberían utilizarse elipsoides
regionales para que no se introduzcan grandes distorsiones en las mediciones 23.
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.
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.
La definición en INSPIRE25 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.
DATUM GEODÉSICO (INSPIRE): un datum que describe la relación de un sistema de coordenadas con la
tierra.
DATUM HORIZONTAL: usualmente, elipsoide determinado con una localización terrestre utilizado en la
determinación de las coordenadas φ y λ.
23 Las normas internacionales hablan de un Datum de Ingenieria. El Catastro no utiliza en España un sistema propio.
24En 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(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. con referencia más detallada a EN ISO 19111:2007, Geographic information — Spatial referencing by coordinates.
DATUM VERTICAL26 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.
“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ó 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).
EPN (EUREF GNSS PERMANENT NETWORK)27 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.
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.
26Ademá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. 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, and the ISO 19111:2007(E)
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 calculation of other ISO 19116:2004(E)
parameters
27El 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.
Elipsoide a 1/f b
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
semieje polar. También en el WGS 84 se utiliza el IERS Reference Meridian y se calculan, como constantes
geométricas derivadas:
f = 0, 003352810681183637418;
relación entre ejes b / a = 0, 996647189318816362;
primera excentricidad (a21/f− b2)5 = 521.854,0097 m.;
segunda excentricidad a2 − b2 / a = 0,0818191910435;
b
cuadrante de un meridiano = 10.001.965,7293 m
Además, se obtienen varias medidas del radio de la tierra esférica, como ya hemos visto.
En la actualidad, ha mejorado el valor del semieje a29 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 del elipsoide que fue
modificada en WGS84 con el cálculo de nueve decimales30.
Mallas geográficas
MALLA, RETÍCULA O REJILLA GEOGRÁFICA. Red de meridianos y paralelos representada mediante líneas
con valores seleccionados a intervalos regulares
28 Puede verse, por ejemplo, en www.ngs.noaa.gov (XIONG LI & HANS-JÜRGEN GÖTZE (2001)).
29Las 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.
30Como 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.
Los arcos de 1º de paralelo disminuyen con el aumento de la latitud pues los círculos tienen un radio que es
menor desde el ecuador a los polos. En la península el arco de 36º N mide 90.163,69 metros y el de 44º N se
reduce a 80.206,19 metros. Al contrario, los meridianos, debido a la elipticidad, aumentan la longitud de sus
arcos desde el ecuador a los polos. El arco de 1º entre 43º N y 44º N mide 111.102,49 metros mientras que el
definido entre 36º N y 37º N se reduce a 110.949,77 metros.
Consecuencia de estas reducciones y aumentos son los valores del perímetro de los trapecios sobre el
elipsoide. Serán iguales los dos lados definidos por los meridianos en cada franja, pero el lado del paralelo
superior es mayor que el arco del inferior. Los valores extremos son 401.111,97 metros en la franja entre 36º N y
37º N y 383.952,13 metros entre 43º N y 44º N31.
Estas son algunas de las características principales de la malla o retícula obtenidas en el GRS_80 que es la
figura curva, cerrada y contínua que ajusta geométricamente la tierra y es la que ha de proyectarse. Son
características geográficas que señalan las diferencias en una malla de líneas curvadas, que no es enteramente
regular, se entiende dentro de la simetría elipsoidal y con un espaciado que es constante entre los paralelos y
cuyos meridianos tienen su mayor separación en el ecuador y sus puntos de convergencia en los polos. Aunque
tenga una extensa descripción es de fácil estudio. Sirve de referencia para la evaluación de las distorsiones en
las proyecciones sobre un plano. Hemos visto en la península que hay una diversidad de situaciones en la malla
elipsoidal que atienden a las diferencias naturales entre las franjas. Se acrecientan en las celdas proyectadas.
GEOREF (también llamado WGRS, World Geographic Reference System). Malla geográfica (no proyectada) de
meridianos y paralelos cada 15º cuyas celdas tienen una denominación.
En GEOREF, las celdas delimitadas por los 24 meridianos se nombran con las letras A- Z del vocabulario inglés,
exceptuando la I y la O. Se comienza en el meridiano -180º y se finaliza en 180º. Las celdas delimitadas por los
12 paralelos, separados por el ecuador32, se nombran con A- M, también con las mismas excepciones. Se
comienza ahora en -90º y se finaliza en 90º.
31Para los cálculos utilizamos una rejilla de 1º en ETRS_89 dentro de lo que ArcMap señala como WKID 4258 (EPSG) con meridiano
central 0º y con parámetros indicados para el GRS_80: a= 6.378.137,0 m, b= 6.356.752,314140356 m y aplanamiento inverso =
298,257222101. El valor de radianes por unidad es 0,017453292519943295. No proporciona otros parámetros pero estos coinciden
con los de MORITZ, H. Geodetic… No hemos averiguado los procedimientos de cálculo, pero no hemos encontrado grandes
diferencias con las formulaciones de una hoja de cálculo incluida en un fichero adjunto. No hay ninguna diferencia en los valores del
paralelo (son arcos de circunferencias) ni en los arcos de meridiano (afectados por la elipticidad). Existen en los cálculos de áreas
(aunque obtenemos una precisión de nueve decimales en lugar de tres). Aceptamos los cálculos de ArcMap en esta y en ocasiones
posteriores pues es una herramienta que es de utilidad en otros aspectos del análisis geográfico y las representaciones
cartográficas..
Las proyecciones de esta rejilla, que mostraremos más adelante (mallas UTM, LCC y LAEA), no siguen el procedimiento de la
proyección ortográfica de Tissot.
32Llamadas, a veces, zonas longitudinales (24) y bandas de latitud (12). Surge “for air defense and strategic air operations”
(http://earth-info.nga.mil/GandG/coordsys/grids/georef.pdf).
A B C D E F G H J K L M N P Q R S T U V W X Y Z
CC Javier Espiago - SCUAM
Cada uno de esos cuadrángulos tiene su propia malla de 1º cuyas celdas vuelven a nombrase de la A- Q en
ambos sentidos.
A B C D E F G H J K L M N P Q
A
B
C
D
E
F
G
H
MJ
J
K
L
M
N
P
Q
CC Javier Espiago - SCUAM
De esa forma, cada celda de 1º se identifica con cuatro letras. Las dos primeras del cuadrángulo de 15º y las
dos segundas de los cuadrángulos de 1º.
La malla continúa con la división en 60 minutos de latitud y de longitud y con denominaciones numéricas de las
celdas ahora comenzando de E-O y de N-S.
La posición con resolución de un minuto se expresa por lo tanto mediante cuatro letras y dos números. Como
son celdas geográficas, esa resolución en metros es desigual. Crece según nos acercamos a los polos.
MAIDENHEAD (también QTH) 34: malla geográfica de 20º de longitud (18 meridianos) y 10º de latitud (18
paralelos) con 324 celdas (18 x 18) con una denominación alfanumérica.
La denominación MAIDENHEAD se indica mediante dos letras desde AA a RR (llamadas el Campo). La primera
expresa la longitud y la segunda, la latitud. El origen de la malla es 180º O y 90º S con alfabetización y
numeración correlativas de O- E y S-N. En la Península Ibérica los dos primeros caracteres son IN, JN, IM e IL
(como se ve, se utiliza la I y también la O).
-180 -140 -100 -60 -20 0 20 60 100 140 180
90 R
80 QP
70 O
60 N
50
M
40
L
30
K
20
10
J
I
0
H
-10
-20
G
-30 F
-40 E
-50 DC
-60 B
-70 A
-80
-90
A B C D E F G H I J K L M N O P Q R
CC Javier Espiago - SCUAM
Cada cuadrángulo se subdivide mediante meridianos de 2º y paralelos de 1º en 10 x 10= 100 celdas a las que
se asignan dos números de 00 a 99. A las dos letras del Campo se les añaden estos dos números del que ahora
llaman el Cuadrado. Por ejemplo, el primer cuadrángulo AA (que en la figura proyectada aparece como un
rectángulo) adquiere la numeración:
9
8
7
6
5
4
3
2
1
0
0 1 2 3 4 5 6 7 8 9
CC Javier Espiago - SCUAM
33“A disadvantage of any standard system of location is that the enemy, if he intercepts one of our messages using the system, can
interpret the message and find our location. It is possible and can be eliminated by using an authorized low-level numerical code to
express locations. Army Regulation 380-40 outlines the procedures for obtaining authorized codes”, http://www.map-reading.com/ch4-
4.php.
34Utilizada en el sistema radiofónico de localizacion Q conocido también como Sistema de Localización QTH procedente del QRA
ahora en desuso.
X
60,0’ W
57,5’
V
55,0’
52,5’
U
50,0’ T
47,5’ S
45,0’ R
42,5’ Q
40,0’ P
37,5’ O
35,0’ N
32,5’ M
30,0’ L
27,5’ K
25,0’ J
22,5’ I
20,0’
G
17,5’
15,0’
H
12,5’ F
10,0’ E
7,5’ D
5,0’ C
2,5’ B
0,0’ A
A B C D E F G H I J K L M N O P Q R S T U V W X
Si se desea una referencia de mayor precisión, se procede a una nueva subdivisión en 10 x 10 = 100 celdas que
se numeran del 0 al 99 como habíamos hecho anteriormente. También pueden repetirse otra subdivisión
alfabética y una nueva rejilla numérica36.
GARS (Global Area Reference System) pese al nombre es un sistema de localización de los USA, aunque
abarca todo el globo37 definiendo una malla geográfica de 30’ x 30’.
En el caso de GARS, la identificación es con cinco caracteres (tres numéricos y dos alfabéticos). Los tres
primeros desde 000 al 720 nombran los husos de los meridianos cada 30’ desde 180º E. Con los dos caracteres
indicamos las bandas de los paralelos cada 30’ comenzando en el polo sur desde AA hasta QZ.
Estas celdas de 30’ se dividen de nuevo en cuatro cuadrángulos de 15’ x 15’ que se nombran secuencialmente
de O – E y de N – S (1, 2, 3 y 4). Con el mismo procedimiento añadimos una malla de nueve celdas de 5’ x 5’
numeradas del 1 al 9. El tedioso resultado “resembles a telephone keypad”. Se incluye el siguiente gráfico
35Son localizadores para intercambios radiofónicos de 10-GHz. “This level of locator detail is what people generally mean when
referring to grid squares” en http://www.vklogger.com/grid_squares_info.php. También sirve para rastreo de telefonía celular, software
SIG y GPS. Existen numerosos convertidores con precisión. También sirven esas herramientas informáticas. En cambio, no
contamos con buenos textos de explicación y suficientes gráficos o imágenes. Aunque se centra en Australia, esa página en red
puede ser de utilidad.
36Esbozado el detalladísimo procedimiento teórico, conviene acudir a herramientas informatizadas. No es infrecuente que haya
radioaficionados que encuentren sus localizaciones empleando elipsoides variados que dificultan la comparación de las trabajosas
coordenadas.
37“For use across the United States Department of Defense” “to facilitate air-to-ground coordination, deconfliction, integration, and
synchronization”. http://earth-info.nga.mil/GandG/coordsys/grids/gars.html también para lo que sigue. Tuvo como antecedente el
CGRS Common Geographic Reference System dividido en numerosos esquemas que no eran fácilmente compatibles.
Con estas tres mallas (GEOREF, MAIDENHEAD, GARS) se muestra una localización de celdas que también
recibe el nombre de coordenadas, aunque no se permitan los cálculos propios de un sistema enteramente
numérico sobre puntos. En cartografía temática las superficies de las celdas encuentran frecuentes aplicaciones.
No siempre se necesita que vayan acompañadas por esas denominaciones siempre engorrosas.
Mediante un convertidor38 ilustramos con los casos de referencia 40°32'43.7284", -003°41'58.7942" del Servicio
de Cartografía-SCUAM y 40°28′55″N y 3°21′51″O en Alcalá de Henares:
38Utilizamos www.earthpoint.us/Convert.aspx herramienta vinculada a Google Earth cuyos resultados se cartografían con sus mapas
e imágenes.Con menor precisión, http://www.egloff.eu/googlemap_v3/carto.php que no exige una instalación previa.
Para las divisiones geométricas del ámbito terrestre de España, además de la malla geográfica y de las mallas
proyectadas del sistema UTM, se han decretado otras que atienden a la división en hojas en las diferentes
escalas 39 . Son mallas vectoriales de arcos regulares de meridianos y paralelos que se utilizaban para la
impresión de los mapas oficiales. Actualmente, también sirven como ámbitos para el registro o almacenamiento
de datos. Sigue siendo necesaria la división territorial para la compilación y como unidades territoriales de
proyección.
Son las hojas del Instituto Geográfico Nacional utilizadas como referencia en otras cartografías siendo acertada
la división en coordenadas elipsoidales que facilita el almacenamiento en bases de datos. Los geógrafos
acostumbrados a trabajar con hojas 1: 50.000 las pueden utilizar como referencia que se manifiesta sobre el
terreno sin la completa abstracción de meridianos y paralelos.
Hemos calculado las longitudes de esos arcos geográficos y las medidas de las rectas una vez proyectados en
el sistema UTM. Vuelve a ser un útil cálculo numérico de la distorsión.
Para la escala 1: 50.000 (MTN_50) se usa una malla geográfica analítica de 20’ en longitud y 10’ en latitud,
siguiendo los meridianos y paralelos. Su origen se sitúa en un punto del mar a -9º 51' 15'' y 44º 00' 00'',
respectivamente. Para la designación de cada una de las celdas (hojas) se utilizan procedimientos de
numeración correlativa del Instituto Geográfico Nacional y de numeración columna-fila del Centro Geográfico del
Ejército de Tierra. No se ha decretado la designación de cada hoja mediante el nombre de la entidad de
población con mayor población, pero también se usa como referencia.
La malla a escala 1:25.000 tiene un tamaño de hojas de 5' x 10' dividiendo la cuadrícula MTN50 en cuatro
cuartos. A escala 1:10.000 el tamaño es de 2' 30'' x 5' (división en dieciséis partes) y a escala 1:5.000, de 1' 15''
x 2' 30'' (sesenta y cuatro partes)40. Como ejemplo, para la hoja 559 el RD 1071/2007 proporciona el siguiente
gráfico:
En el paralelo 36º 0’ 0,00’’ de latitud miden 30.054,56 m de Este a Oeste y en la latitud máxima de 43º
50’ 0,00’’ se reducen hasta 26.735,40 m. Esta última son hojas de una fila teórica enteramente
marítima. Las tres hojas continentales de la segunda fila (43º 50’) miden 26.810,12 m.
La extensión de los arcos de meridiano límite de las hojas varía desde 18.493,42 m hasta 18.517,89 m.
Esta es la división en hojas para la compilación y proyección oficial que hoy emplea UTM ETRS_89 en los tres
husos de España peninsular y Baleares. No tendría que utilizarse otro sistema 41.
Javier Espiago
Servicio de Cartografía de la Universidad Autónoma de Madrid.
Julio - 2017