Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 1 de 17

Miniature Inertial Sensor Attitude Heading Reference System with GPS, AHRS, IMU, Tilt sensors Yellowbrick 3 GPS Tracker A truly global tracker & messenger. Tracking and messaging over Irid GPS Tracker + Software Want to Build Tracking Service? Easy & Save 50% GPRS Cost. www.bofa

:: Conversin de datum con el modelo de 7 parmetros Burs

Hace unos meses publiqu en esta pgina un artculo detallando el proceso de transformacin de coordenadas geodsicas trataba el tema del traspaso de coordenadas que estuvieran referidas a datums diferentes. Es decir, tratbamos el desarr UTM referidas a un elipsoide determinado (el que fuera), y cuyo destino eran sus coordenadas geodsicas equivalentes so proceso inverso. No se abordaba, sin embargo, la posibilidad de que las coordenadas origen estuvieran dadas sobre un elips diferente. Este proceso, conocido como cambio de datum o traspaso de datum, est enunciado en muchos libros de texto y pginas we mtodos que existen para su realizacin con sus correspondientes ecuaciones; yo no he llegado a encontrar un texto con proceso de principio a fin. Quiz, debido a la extensin del proceso, que es larga, la mayor parte de los textos encontrados las ecuaciones, pero no entran en la totalidad del proceso. De esta forma, se pasan por alto numerosos detalles que son imp ejecucin. S que hay bastantes dudas sobre este tema, as que voy a intentar escribir un texto que describa TODO el proceso de p detallado, y adjuntado al final del artculo una hoja de clculo en formato Microsoft Excel con la implementacin dije en el artculo de conversin de coordenadas, el hecho de adjuntar una hoja de clculo no pretende que sta sirva como sirva para comprender el funcionamiento de las ecuaciones y que facilite a los programadores el desarrollo de los conversor de las ecuaciones es ms fcil de entender que poner cdigo fuente directamente.
Si vas a hacer un uso comercial del contenido, puedes realizar una donacin voluntaria al proyecto en http://r contenido del texto que vas a ver te resulta til y te ayuda a aprender, te pido que cites esta pgina web en tu trabajo.

En cuanto a la estructuracin del artculo, realizaremos en primer lugar un ligero encuadre terico de qu es del datum y lo q lugar, enunciaremos los distintos mtodos que existen para realizar traspasos entre datums y, en tercer lugar, realizaremos mtodos: el modelo de ecuaciones de 7 parmetros de Bursa-Wolf, que es uno de los procedimientos ms precisos para esta

1.CONCEPTOS PREVIOS: DATUM Y CAMBIO DE DATUM.

Es necesario realizar un encuadre terico de lo que significa y lo que implica el traspaso de datum. Para ello, es nece superficies de referencia que existen en cartografa.

Superficies de Referencia en Cartografa.

El proceso de formacin de la cartografa implica una abstraccin de la realidad. Necesariamente, tenemos que simplif introducirla y grafiarla sobre un mapa. Indudablemente, este proceso de abstraccin implica muchos niveles del proceso ca ms profundos desde el punto de vista terico, consiste en simplificar la propia superficie de referencia sobre la que trabajam Como es sabido, los mapas necesitan de un modelo matemtico (sistema de proyeccin) que permita el paso de una superfi a un medio plano como es el papel. Para poder aplicar cualquiera de los modelos matemticos que permiten tal paso (o complejidad de la superficie de la tierra a una superficie de aproximacin que sea modelizable matemticamente. Sera impo la tierra, que est plagada de accidentes geogrficos como sierras, valles, acantilados, etc. en una sucesin de complejidad in A partir de esta necesidad de encontrar una primera simplificacin de la superficie de referencia de la tierra, los geodestas concepto de geoide. El geoide podemos imaginarlo como la superficie que observaramos si el mar estuviera en calm imaginariamente por debajo de los ocanos. Se trata de una superficie equipotencial, es decir, en la cual el potencial de la

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, ...

29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 2 de 17

aunque no existe fsicamente (los mares nunca estn en calma total, sin mareas y tampoco atraviesan los ocanos por mediciones empricas.

Esta medicin del geoide no es fcil de hacer, e histricamente ha sido uno de los principales problemas de la geodesia imaginario de cmo podramos precisar esa superficie. Imaginemos que tenemos un mar en calma y sin mareas e imagine ese mar. Sostenemos una cuerda con un peso (una plomada); esa plomada sealar la vertical del lugar (o vertical prolongarla hasta tocar en un punto a la superficie del mar. Justo por ese punto de toque pasara la superficie del geoide, geoide la superficie de ste es perpendicular a la direccin de la vertical del lugar:

Y si realizramos infinitas medidas como la anterior alrededor de todo el planeta? cuantas ms medidas realizramos en superficie del geoide. En cada punto del planeta, la vertical astronmica sealara en una direccin diferente:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 3 de 17

Obviamente, este procedimiento imaginario no es viable desde el punto de vista prctico... Sin embargo, se han desa evolucionado sobre todo con el advenimiento de la geodesia por satlite. Por ejemplo, esta es una imagen del geoide depresiones y las zonas rojas como elevaciones:

Ampliar imagen

Montaas y valles en el mar? s, algo parecido... En la configuracin del geoide intervienen muchos factores: la diferente de que rige los complejos movimientos de los materiales fundidos del magma y ncleo del planeta, la fuerza centrfuga de la ro la gravedad se presente con un patrn irregular. Por eso es tan difcil delimitar el paso del geoide, y por eso es tan impo geodesia fsica. Ah est el uno de los puntos de conexin entre la fsica y la cartografa. Al observar la imagen de arriba, vemos que el geoide tampoco es una superficie muy regular. Ms bien tendra la forma de esa configuracin, no debe ser fcil modelizar su forma matemticamente. Otra rama de la geodesia, la geodesia geomtrica, se ha encargado de conseguir una superficie de aproximacin ms elipsoide. El elipsoide es una figura tridimensional que intenta reproducir la geometra del geoide, pero con la particularidad de ecuaciones. Tiene su origen en una elipse a la que imaginariamente hacemos rotar entorno a su eje ms corto:

A partir del giro de la elipse sobre su eje menor, se genera el elipsoide de revolucin que ya es una figura tridimensional. E longitud de sus dos semiejes. A partir de ellos se sacan otros valores de su geometra muy utilizados en los problemas ge segunda excentricidad, etc. En realidad, comprender la geometra de la elipse y del elipsoide abre las puertas a comprender Sobre el elipsoide se miden las coordenadas geodsicas, expresadas en grados, minutos y segundos sexagesimales. Par geodsica, que es la vertical perpendicular al elipsoide. A partir de ella se fija nada menos que la latitud geodsica (repres smbolos son en realidad la misma letra griega cuya lectura se pronuncia "fi").

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 4 de 17

Recapitulando, tenemos tres superficies a considerar: la superficie que vemos de la tierra, el geoide y el elipsoide. Estas t variando sus posiciones relativas constantemente. A la distancia existente entre un punto del elipsoide y su punto homnim geoide, desviacin geoidal, o altura geoidal:

Y esta diferencia es muy importante para muchas operaciones cartogrficas y geodsicas. Por poner uno de los ejemplos vemos en nuestros mapas (curvas de nivel, puntos acotados, etc.) se establece normalmente en el geoide. En concreto, se en Espaa el origen de las altitudes cartogrficas es el nivel medio del mar en Alicante (NMMA). A estas altitudes medidas ortomtricas. Por otro lado, existen las altitudes medidas a partir del elipsoide y que se vienen usando en los ltimos aos y cada vez m GPS. Estas altitudes presentan desfases notables con respecto a las ortomtricas, puesto que como hemos dicho, l constantemente a lo largo de toda la tierra, y eso sucede con cualquier elipsoide considerado. La diferencia -ondulacin del g metros. Al punto tomado como referencia para medir las altitudes cartogrficas o cotas Z se le denomina datum vertical. Vemos q datum vertical el que ocupa el presente artculo. El datum que nos interesa en este texto es el datum horizontal.

Qu es el datum cartogrfico o datum horizontal?

Bueno, ya hemos dicho que la superficie del elipsoide trata de ajustarse a la superficie del geoide en la medida de lo posible. versiones distintas del elipsoide (cientos de elipsoides, cada uno con sus parmetros de semieje mayor, semieje menor y ubicar tridimensionalmente el elipsoide que consideremos con respecto al geoide? Cmo situar uno con respecto al otro? pu Eso es precisamente lo que hace el datum: fija el elipsoide a utilizar con sus parmetros geomtricos correspondientes elipsoide asociado a un datum), y establece el punto en el que la vertical geodsica coincide con la vertical astronmica. Es d relacin entre geoide y elipsoide, que normalmente se corresponde con un punto de tangencia entre ambas superficies. relacin se le denomina punto astronmico fundamental. En el caso de Espaa, hasta fechas recientes el datum que s datum europeo de 1950 (ED50), que toma como elipsoide de referencia el elipsoide de Hayford (o elipsoide internacional de y que establece el punto astronmico fundamental en la ciudad alemana de Potsdam; en concreto, en la Torre de Helmert campus universitario de esa ciudad, y que se puede ver en la foto de la derecha. No es casualidad que sea un observator lugar en que se pone en relacin el geoide con el elipsoide. Lgicamente, es necesario conocer con mucha precisin la ver para orientar con precisin el elipsoide, pues es el punto inicial de materializacin del sistema geodsico de referenci implicacin an ms profunda: si la vertical astronmica y la vertical geodsica coinciden, significa que tambin las coordena las coordenadas astronmicas coincidirn. De ah que los puntos fundamentales de los datums estn frecuentemente astronmicos y de ah procede la estrecha relacin entre la astronoma y la geodesia, que se plasma en la rama de la geodesia astronmica. La siguiente idea importante es que un mismo punto de la superficie de la tierra, tiene diferentes coordenadas geodsica datum elegido:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 5 de 17

As, las coordenadas tomadas con GPS (referidas al datum WGS84), pueden diferir cientos de metros de su posicin rea tomadas. Puede que los valores de las coordenadas sean de buena calidad, pero si el usuario equivoca el datum de referen metros:

Por tanto, se deduce que conocer el datum de referencia de las coordenadas es un parmetro fundamental para mante capacidad de traducir datos cartogrficos de un datum a otro es una habilidad que todo gegrafo ha de conocer.

2. TIPOS DE TRANSFORMACIONES DE DATUM.

Existen diversos mtodos para realizar el traspaso de unas coordenadas en un datum a otro de destino. Unas operan con co centro geomtrico del elipsoide de referencia considerado (no hay que confundir estas coordenadas geocntricas X Y Z con c UTM); otras, operan directamente sobre las coordenadas geodsicas expresadas en grados, minutos y segundos; por ltim proyectadas, si es necesario. Las que operan con las coordenadas geodsicas directamente son las Ecuaciones de Molodensky (o Molodenski, porq ecuaciones existe una versin simplificada denominada comnmente Abriged-Molodensky, aunque a veces tambin se secas, lo cual introduce una cierta confusin. La versin reducida, viene a ser comparable a lo que veremos a continuacin decir, vlida para precisiones bajas. La transformacin de Molodensky completa es la que ofrece mayores garantas de precis En cuanto a las transformaciones que operan con coordenadas proyectadas directamente (si es necesario), se encuent emplean expresiones de distinto grado de polinomio para hacer el mejor ajuste entre dos sistemas de coordenadas que no

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 6 de 17

La versin ms sencilla de estas transformaciones sera una transformacin afn, que en esencia es la aplicacin de transformacin de este tipo poco puede hacer por nosotros ante un problema tan complicado como un cambio de datum... ms complejas s que pueden resolver un complejo cambio de datum. De hecho, en Espaa hace algunas dcadas se tuvo q fue migrar toda la cartografa del datum de Madrid 1870 (referido al meridiano de Madrid, concretamente con punto fundam Astronmico de esta ciudad y que empleaba el elipsode de Struve), al datum europeo unificado de 1950 con punto fundam Hayford. El paso de datum se resolvi usando trasformaciones polinmicas desarrolladas por el Servicio Geogrfico del Ejrc para el territorio penisular espaol y estn parcialmente documentadas aqu. En aquella poca, adems, coincidi esta prob de representacin de la proyeccin polidrica que se usaba hasta entonces a la adopcin de la proyeccin UTM. Slo vamos a entrar en profundidad en las transformaciones que utilizan coordenadas cartesianas rectangulares. Veamos qu Las transformaciones que operan con coordenadas cartesianas, requieren que estas coordenadas estn referidas al centro g en trminos geogrficos se denominan coordenadas ECEF (Earth Centered-Earth Fixed), que en castellano podramos traduci se llaman tambin coordenadas cartesianas o coordenadas rectangulares). Podemos imaginar el marco de referencia de esta de ejes con centro en el elipsoide de revolucin:

As, cuando queremos hacer el traspaso entre coordenadas referidas a dos datums diferentes, lo primero que tenemos qu lugar, pasarlas a geocntricas; en tercer lugar, hacer el cambio de datum mediante la compensacin de las diferencia tridimensionales geocntricos:

Luego veremos en detalle cmo se ejecuta todo este proceso, pero conceptualmente podemos ver que supone tener que que se descomponen en un incremento de distancia en X, un incremento en Y y otro en Z. Tambin supone tener que com

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 7 de 17

homnimo en el otro sistema. Si los sistemas de ejes fueran paralelos, estos ngulos seran inexistentes, pero si no lo son que presentan entre s. Y por ltimo, la diferencia de tamao entre las unidades de un sistema y otro, que se denomina fac diferencias de la siguiente manera:

En el grfico est sin simbolizar la diferencia en escala entre ambos sistemas. Imaginemos que el sistema de ejes colorea embargo, el otro sistema tiene una unidad equivalente a 0.997 metros; el factor de escala para pasar del sistema verde al multiplicar las coordenadas verdes por 1.003 En consecuencia, con la aproximacin aqu expuesta son 7 los parmetros que tenemos que considerar para hacer el incrementos de distancia en X, en Y y en Z, rotaciones en los ejes X, Y, Z, y por ltimo el factor de escala. La mayor parte de la diferencia entre ambos sistemas se compensa teniendo en cuenta la diferencia de los centros (los desp Z). De hecho, una transformacin sencilla de datums consiste en pasar las coordenadas geodsicas a geocntricas y una ve distancia al origen. Este tipo de transformacin, suele denominarse de 3 parmetros y es slo recomendable para aplicacione Por ejemplo, el conversor de coordenadas en lnea que aparece a la izquierda de esta pgina utiliza una transformacin de 3 datum. En ese conversor de coordenadas, la conversin en el mismo datum es totalmente precisa, pero como digo si trasl Por qu entonces se usa una transformacin de 3 parmetros si no es precisa? Porque utiliza unos parmetros que son existen tablas que proporcionan las diferencias de origen entre datums: en este ejemplo tenemos las diferencias de varios da Para aplicar una transformacin de 7 parmetros, sin embargo, no valen unos parmetros genricos de rotaciones de ejes unos concretos de la zona sobre la que estamos trabajando. Estos parmetros lgicamente tienen que hacer referencia a la vamos a traducir y de su calidad dependen directamente los resultados.

3. TRANSFORMACIONES DE DATUM CON 7 PARMETROS.

Dentro de las transformaciones de 7 parmetros usadas en geodesia, existen dos alternativas de formulacin y una cierta casi idnticas, son las ecuaciones de Bursa-Wolf y la transformacin de Helmert. Son prcticamente iguales y slo presenta de rotaciones y sus signos:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 8 de 17

En cuanto a los nombres, existen denominaciones alternativas que crean una cierta confusin. As, hay fuentes que a transformacin de Bursa-Wolf y el nombre alternativo Position Vector a la de Helmert (ese criterio asumo yo en este artcul estos nombres cruzados e incluso se cruzan tambin las denominaciones originales con respecto a las fomulaciones de las datos de transformaciones del EPSG (European Petroleum Survey Group, ahora llamado OGP), existe confusin y disp importante. Lo importante es saber que dentro de la formulacin Bursa-Wolf existe otra variante, y que la forma de contar l o contrario a las agujas del reloj) es importante para el funcionamiento de las ecuaciones. Tambin, hay que considerar que de transformacin: consideran un datum origen y un datum destino, y por lo tanto si queremos revertir la conversin hay qu Todas estas consideraciones obligan a una cierta cautela cuando recibimos un conjunto de parmetros para realizar una t saber a qu tipo de transformacin pertenece un conjunto de parmetros slo con verlos, y tambin hay que tener cui metadatos, aunque existan, no garantizan esta parte del trabajo, pues como he dicho anteriormente existen diferencias a n ser b. Se hace necesario, pues, comprobar con las dos variantes de formulacin anteriormente expuestas, y hacer pruebas v coordenadas conocidas en los dos sistemas, podemos ver si el funcionamiento es correcto o no, y cul es el error que nos es Por ltimo, y antes de cerrar este punto sobre las transformaciones de 7 parmetros hay que advertir tambin de sus transformaciones de siete parmetros que hemos visto realizan una rotacin tomando como punto de rotacin el centro del e los parmetros de desplazamiento y giro a partir de una serie de puntos de control con coordenadas en ambos sistemas, ex rotaciones. Para eliminar esta dependencia entre los parmetros, se utiliza una tercera variacin de las ecuaciones de Badekas, que con la adopcin de 3 parmetros adicionales pasa a usar un total de 10 en el proceso. Estos tres nuevos parmetros de la variacin Molodensky-Badekas definen el punto de rotacin (Xp, Yp, Zp), que puede se de la zona en la que se estn traduciendo coordenadas (centroide del polgono de actuacin). Este punto se puede sacar media de las X de todos los puntos de control, la media de las Y y la media de las Z, naturalmente con las coordenadas en el La variante Molodensky-Badekas queda como sigue:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 9 de 17

Como es lgico, para usar esta transformacin tenemos que tener parmetros deducidos al efecto y no vale aplicar parmet La incorporacin de un punto de rotacin particularizado requiere de parmetros tambin con este condicionante extraidos co

4. LAS ECUACIONES DE BURSA-WOLF EN FUNCIONAMIENTO: UN EJEMPLO DE

Una vez presentado todo el encuadre terico anterior, estamos en condiciones de iniciar un ejemplo de conversin de dat hemos entrado. Utilizaremos la transformacin de Bursa-Wolf, propuesta por Bursa en 1962 ("The theory of the Determination of the Non Initial Astronomical and Geodetic meridians from observation of Artificial Satellitea", en Studia Geophysica et Geodetica, ("Geometric connection and re-orientation of three-dimensional triangulation nets" en Bulletin Godsique vol. 68, pp. 65-16 da por los paquetes GIS y aunque arriba hemos hecho matizaciones sobre su uso, ofrece un nivel de calidad ms que suficie El ejemplo lo haremos traduciendo las coordenadas de un vrtice geodsico del centro de Espaa: el vrtice de Carbonera. en dos datum distintos, a partir de clculos realizados por el Instituto Geogrfico Nacional. Se trata de uno de los vrtices d en ED50 y coordenadas en ETRS89 (~WGS84):
Vrtice de Carbonera Datum ED50 Datum ETRS89 X UTM 448611.14 448500.79 Y UTM Altitud -3 35' 5 -3 35' 5 4377788.61 718.50 (ortomtrica, sobre NMMA) 4377580.93 771.46 (elipsoidal, sobre GRS80)

Como se ve, tenemos coordenadas en UTM y en geodsicas, y en los dos datum con los que vamos a trabajar. Hay que dec unificado que se usa hoy da en Europa, es en la prctica equivalente a WGS84, por lo que el presente ejercicio es asimilab datum WGS84 a ED50. Se trata en todos los casos de coordenadas seguras, en el sentido de que han sido calculadas con una precisin centimt escogido hacer el ejercicio con un vrtice geodsico situado en el centro de la Pennsula Ibrica, donde se supone que la re adems los parmetros de transformacin (tambin facilitados por el IGN) pueden funcionar mejor. Dichos parmetros calculados por el IGN se distribuyen en tres conjuntos, cada uno aplicable a una zona de la Pennsula Ib

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 10 de 17

Los desplazamientos se dan en metros, las rotaciones en segundos sexagesimales (que luego tendremos que transformar a milln (que luego tendremos que dividir por 1 milln, o lo que es lo mismo, multiplicar por 10-6). Estos parmetros sirven error puede estar entorno al metro o metro y medio. Los parmetros pueden ser de mayor precisin cuanto ms reducido es
Cmo y dnde conseguir los parmetros de transformacin? Normalmente los usuarios han de solicitarlos para cada zona concreta del globo a su institucin nacional o regional encarg de Espaa, existe una aplicacin especfica para servir parmetros en la web del Centro Nacional de Informacin Geog parmetros seguros, pero tambin podemos deducirlos a partir de los modelos de ecuaciones anteriormente presentados. coordenadas en los dos sistemas y pasaramos a despejar las incgnitas que seran los trminos de rotacin, traslacin y abordemos prcticamente este tema.

En este artculo asumiremos los parmetros denominados "Pennsula" de la tabla arriba expuesta, que como hemos dicho sirv El ejercicio que vamos a hacer procede de coordenadas UTM en ETRS89 (~WGS84) y tiene destino en coordenadas UTM p Para hacer este proceso, necesitamos realizar las siguiente secuencia general:

Paso de UTM en datum origen a Geodsicas en el mismo datum origen.

En un artculo de esta web se explicaba cmo realizar el paso de UTM a Geodsicas en el mismo datum, por lo que no descr dicho procedimiento, puede acceder al artculo aqu. Tan slo baste decir que las coordenadas geodsicas correspondientes al vrtice de Carbonera son, segn el mtodo descrito

Vrtice de Carbonera Datum ETRS89

X UTM 448500.79

Y UTM 4377580.93

-3 35' 57''7336

39 32' 46

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 11 de 17

Existe una mnima diferencia en las geodsicas calculadas por este procedimiento y las geodsicas anunciadas por el I como digo es mnima y totalmente compatible con la aplicacin, se debe al diferente procedimiento de clculo.

El valor de las altitudes en el proceso y su clculo.

Antes de seguir adelante con el proceso, es necesario realizar algunas consideraciones sobre las altitudes. Para operar conocer las altitudes del punto pero sobre el elipsoide. Antes hemos dicho que el elipsoide presenta una diferencia con Tambin hemos dicho, que las altitudes que nos muestran los mapas suelen ser altitudes medidas con origen en la superfic que los GPS nos suelen dar las altitudes sobre el elipsoide (altitudes elipsoidales). En la transformacin nos hacen falta altitudes elipsoidales, que podemos calcular si conocemos la ondulacin del geoide en e

Dado que:

donde h es la altitud elipsoidal, H es la altitud ortomtrica y N es la ondulacin del geoide. Podemos obtener una aproximacin a las alturas elipsoidales utilizando un modelo de geoide en lnea. Para ello, recomien Geospatial-Intelligence Agency. En este modelo, a partir de las coordenadas geodsicas en WGS84 (~ETRS89), podemo elipsoide WGS84 (~GRS80) en ese punto:

Ver detalle del uso del modelo EGM96

Con este modelo, obtenemos la altura geoidal que es de 53.26 m, luego la altura elipsoidal ser:

El nmero 718.50 sale de la altitud ortomtrica dada por el IGN o consultable en los mapas, y que en el caso de Espaa Alicante (NMMA). El valor h de altitud elipsoidal calculado por este procedimiento no difiere mucho con el calculado por e aceptado como vlido para el tipo de aplicacin de precisin media que suponemos en este artculo. Tambin podramos realizadas con GPS de alta precisin sobre el punto.

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 12 de 17

Por lo tanto, ya tenemos tambin altitud elipsoidal Z para el punto. Pero... cul es el valor realmente de las altitudes en pueden convertir coordenadas de un datum a otro sin conocer las altitudes elipsoidales? Bien, yo creo que el valor de las altitudes elipsoidales en el proceso de conversin y si estas son estrictamente necesarias o que nos ocupa en el presente artculo en el que manejamos parmetros vlidos para una superficie muy grande de territorio transformacin con una precisin del orden de un metro o metro y medio, las altitudes elipsoidales no son totalmente nec tiene un rango de mil y algo metros, es mnimo comparado con los aportes de los ejes X e Y, que tienen rangos mucho mayo He plasmado en la siguiente tabla los resultados comparados de las transformaciones para varios puntos de la Pennsula Ibr datum utilizando altitud elipsoidal y sin usar altitud (osea, suponiendo h = 0). Los errores se sacan de comparar las coorde seguras ED50 del IGN:

Como se puede ver, los resultados casi no varan si usamos altitudes o suponemos estas igual a 0. Las disminuciones o coordenadas ED50 facilitadas por el IGN) son mnimas.
Insisto en que esta afirmacin es slo vlida para el caso que nos ocupa. No estoy afirmando que las altitudes elipsoid modelo Bursa-Wolf. En otras situaciones diferentes al escenario planteado en este artculo las altitudes elipsoidales s cuando operamos con parmetros para reas de territorio ms reducidas.

El recurso de suponer las altitudes elipsoidales iguales a 0 est siempre ah, pero depende del usuario evaluar la incidencia d perfectamente vlido operar con una altitud = 0 y otras veces no. En este artculo s utilizaremos la altitud elipsoidal de 77 resto del proceso.

Paso de geodsicas a geocntricas en datum origen.

Anteriormente ya hemos explicado lo que supone pasar las coordenadas geodsicas a geocntricas. Veamos a continuacin cartesianas geocntricas (ECEF):

Donde:

Para una mayor claridad expositiva, pasar a escribir en negro las ecuaciones originales y en color azul las aplicacione concretos del ejemplo.

Conocemos tambin los datos correspondientes al elipsoide en el que estn las coordenadas de origen. Como hemos dic GRS80, que es prcticamente igual al elipsoide de referencia del sistema WGS84. Por lo tanto, tomaremos los parmetr GRS80:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 13 de 17

Comenzamos con el clculo de N:

Aplicando al ejemplo que nos ocupa:

Es importante destacar que tanto la latitud como la longitud han sido pasadas a grados sexagesimales en notacin dec la forma: grados decimales = grados + minutos/60 + segundos/60/60. Como se puede ver, los datos angulares as obtenido tambin modificados para introducirlos en radianes. Por eso los multiplicamos por el nmero (pi) y dividimos por 180. D pasar todas las medidas angulares a grados decimales y luego a radianes. Ahora, ya podemos obtener las coordenadas geocntricas cartesianas, segn las formulaciones expuestas arriba:

Cambio de datum con la aplicacin del modelo Bursa-Wolf.

Ahora que ya disponemos de las coordenadas geocntricas cartesianas ECEF podemos aplicar el modelo de Bursa-Wolf directamente sobre coordenadas proyectadas UTM X,Y,Z, por lo que es imprescindible pasar por las geocntricas. Ya hemos dicho antes que la formulacin del modelo es como sigue:

De lo que se deduce que las soluciones de X' Y' y Z' seran:

Hay que tener en cuenta que los parmetros de transformacin nos proporcionan los giros Rx, Ry, Rz en segundos sexag grados decimales y luego a radianes (grados decimales = segundos/60/60; radianes = grados decimales * / 180). Tam escala nos viene en partes por milln, por lo que tendremos que dividir la cantidad que nos facilitan por un milln para pode esto, transformamos los parmetros originales:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 14 de 17

Con los parmetros ya transformados, procedemos a resolver el cambio de datum:

Estas coordenadas siguen siendo cartesianas geocntricas ECEF, pero estn ya en el datum de destino ED50. Hemos rea queda es el proceso inverso de conversin de geocntricas a geodsicas y posteriormente a UTM, con lo que terminaremos e

Paso de geocntricas a geodsicas en datum destino.

Para el paso de coordenadas geocntricas a geodsicas, utilizaremos las siguientes frmulas:

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 15 de 17

Como vemos, se trata de sacar tres parmetros: latitud geodsica, longitud geodsica y altitud elipsoidal. Para ello, no ob parmetros que son necesarios. Como ya hemos realizado el traspaso de datum, estamos operando en el datum de de parmetros del semieje mayor y semieje menor del elipsoide asociado, que es el elipsoide de Hayford, tambin llamado In parmetros son:

Aplicando las frmulas expuestas para el clculo de la excentricidad al cuadrado y la segunda excentricidad al cuadarado, ten

Tambin podemos calcular los dos valores auxiliares p y :

Con estos valores, ya estamos en disposicin de calcular la latitud geodsica (), la longitud geodsica (), segn las formula

Una vez que tenemos latitud y longitud geodsicas, podemos calcular tambin el valor de N:

Conociendo N, calculamos tambin la altitud geodsica (h):

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 16 de 17

Como cabe suponer, la latitud y longitud geodsicas obtenidas estn expresadas en radianes, y para obtener su valor en obtenidas por el cociente de /180:

Si pasamos los grados sexagesimales expresados en notacin decimal a expresarlos en grados, minutos y segundos, obtenem

Paso de Geodsicas a UTM en el mismo datum destino.

En un artculo de esta web se explicaba cmo realizar el paso de Geodsicas a UTM con origen y destino en un mismo datum lector precisa conocer dicho procedimiento, puede acceder al artculo aqu. Tan slo baste decir que las coordenadas UTM resultantes de la conversin de las geodsicas anteriormente expuestas result
Vrtice de Carbonera Datum ED50 (convertidas con transformacin de 7 parmetros) X UTM 448610.60 Y UTM 4377788.16 -3 35' 53''072

Con estas coordenadas UTM ya convertidas en datum destino hemos concluido el ejemplo de cambio de datum.

5. CONVERSOR DE COORDENADAS EN HOJA DE CLCULO EXCEL.

He creado una hoja de clculo con la implementacin de las frmulas comentadas. Esta hoja, realizada en formato Microsoft Excel, permite ver cmo se van calculando todos los parmetros en cadena cuando se introducen unas determinadas coordenadas a convertir. En ella se incluyen todos los clculos descritos (incluido el paso de UTM a geodsicas y de geodsicas a UTM). Incorporan tres juegos distintos de parmetros adecuados todos ellos para la Espaa peninsular, si bien el usuario puede introducir sus propios parmetros de transformacin para su regin del globo. Tambin se pueden seleccionar diferentes elipsoides a partir de una lista predefinida y con posibilidad de incluir otros nuevos por parte del usuario. Descarga la Hoja de Clculo con la Conversin de Datum.
Si vas a hacer un uso comercial del contenido, realiza una donacin voluntaria al proyecto en http://recursos.gabrielortiz.com/donaciones.htm. Si el contenido del texto te resulta til y te ayuda a aprender, te pido que cites esta pgina web en tu trabajo.

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012

Conversin de datum con el modelo de 7 parmetros Bursa-Wolf

Pgina 17 de 17

6. CONCLUSIONES.
Es fundamental para el gegrafo conocer el datum en el que se encuentra cada conjunto de datos, as como tener la ha datum a otro, conociendo las implicaciones, riesgos y condicionantes involucrados en el proceso. Existen mltiples procedimientos para realizar los cambios de datum. El modelo Bursa-Wolf con transformacin de 7 parme para este procedimiento ya que tiene capacidad suficiente para realizar transformaciones de media y alta precisin. Si comparamos el resultado de las UTM obtenidas en el ejercicio realizado con las coordenadas calculadas por el IGN de fo que los errores cometidos son de 0.54 m en el valor de X y de 0.45 m en el valor de Y. Es decir, el punto calculado con nu punto calculado por procedimientos ms precisos por el IGN. Esta magnitud de error es un excelente resultado si tenemos e en el proceso, que como se ha dicho sirven para la mayor parte del mbito de la Pennsula Ibrica y tienen una precisin anu Los resultados son tan buenos porque estamos operando con un vrtice geodsico correspondiente a la parte central de la mantiene una calidad mucho mayor que muchos de los vrtices geodsicos de las zonas perifricas costeras de Espaa, d ptima la antigua red de triangulacin por la presencia del mar impone unas mayores tensiones en la red. As, hemos realiza geodsica espaola situados en estas condiciones perifricas y los errores en la transformacin de datum con la metodolog mayores, incluso por encima de la precisin estimada para los parmetros utilizados: Garita (en Cantabria) con un error de error de 2.26 m, Miramundo (en Huelva) con un error de 1.90 m. Esta situacin contrasta con las precisiones de los vrtices comportamiento de los parmetros de transformacin es sustancialmente superior. Ya hemos hablado del vrtice de Carbon Cerrillos (en Madrid) con un error de 0.35 m. Esto demuestra la importancia de realizar la transformacin de datum con unos parmetros adecuados a nuestra zona y a metodologa expuesta anteriormente si no disponemos de los parmetros adecuados y con el nivel de precisin apropiado pa Por ltimo, termino como empec: pidindote que si el contenido de este artculo te ha resultado til y te ha ay pgina web en tu trabajo. Si el contenido te sirve para una aplicacin comercial re http://recursos.gabrielortiz.com/donaciones.htm. Para cualquier comentario, puedes encontrarme en gis@gabrielortiz.com. Gabriel Ortiz.

Esta informacin te ha sido til? Aydanos a mantener la pgina con una donacin:

GabrielOrtiz.com

mhtml:file://E:\manual\GIS - SIG ArcGIS, ArcView, MapInfo, Geomedia, MapServer, A... 29-02-2012