Título Principal/ Main Heading
SISTEMAS DE COORDENADAS ASOCIADOS A MICROGEODESIA Y
LA GENERACIÓN DE ECUACION DE UN PROYECTO: UN ANALISIS
COMPARATIVO ENTRE METODOS PROYECTIVOS Y ESPACIALES.
Nombre de autor / Name
Rodrigo Severino Arce
Institución / Institution
Zenit Ingenieros Geomensores Ltda.
Dirección / Address
[email protected] Avenida José Miguel Carrera 5161 of. 305 San Miguel
Resumen / Abstract
Cuando se busca algo, la pregunta es categórica ¿Donde está?, y su respuesta se
compone de información precisa, asociada a la ubicación Sin duda, es la pregunta
que pueden hacer las coordenadas de un proyecto vial a las de una cartografía
general, o bien a las de un plano de riesgos naturales. Al conversar estas fuentes de
información entre si se encuentran con una torre de Babel, distintas lenguas para
decir “donde esta”. Sin duda nos estamos adentrando ya específicamente en el
profundo mundo de los sistemas de coordenadas y como se traducen unos a otros,
la georreferenciacion de las fuentes de información. Este trabajo relaciona dos
conocidas herramientas: La transformación de Coordenadas (UTM a PTL) y por otro
lado la transformación de Datum (Bursa-Wolf). El primero con aplicaciones más
directas en proyectos de ingeniería y el segundo con aplicaciones en Geodesia y
finalmente en Cartografía general. Esta es una versión de la traducción que
permite de un modo simple y hasta popular hacer conversar fluidamente a estas
grandes fuentes de información haciendo mínima la expresión de burocracia digital
en la que a menudo se puede caer para su adecuada traducción o se podría decir
también que corresponde a, “Sistemas de Coordenadas Asociados a Microgeodesia
y La Generación de Ecuación de un Proyecto: Un Análisis Comparativo Entre
Métodos Proyectivos Y Espaciales” Se denomina métodos proyectivos a la
Transformación de Coordenadas provenientes de la Proyección Cartográficas
Universal Transversa de Mercator (UTM), al Sistema de Coordenadas en un Plano
Topográfico Local, ampliamente utilizada en Chile. La Importancia de esta relación
se descubre en la aplicación de trabajos Geodésicos, Cartográficos y Topográficos
que requieran de un sistema único el cual a su vez este vinculado a la red
Geodésica Nacional, en palabras simples, crear un sistema de transformación de
coordenadas que haga Compatible proyectos similares. Muchos proyectos
Mineros o de concesiones viales se ven beneficiados de esta nueva manera de
enfrentar los sistemas de Coordenadas. Por otra parte se entiende por
Microgeodesia, a todas aquellas Herramientas Geodésicas que satisfacen fines
Topográficos. Desde este punto de vista la Geodesia Geométrica es quien aporta
la mayoría de los criterios empleados en Proyectos de Grandes Obras Civiles. Sin
embargo, con el correr de los últimos años vemos en la Geodesia Espacial una
inagotable fuente Teórica que contribuye a través de técnicas como los sistemas de
posicionamiento Global (GPS). Finalmente se entregará una herramienta que
permite construir polinomios de transformación o “la ecuación de un proyecto”
simplificando de este modo el traspaso de información georreferenciada en
diferentes sistemas.
Antecedentes Generales
Una ecuación de proyecto, es un polinomio simple, que permite transformar un set de
coordenadas de un sistema a otro. La importancia radica en la simpleza del cálculo y la dificultad
radica en la definición apropiada y debidamente validada de los factores que gobiernan el polinomio
o la denominada “ecuación de proyecto”. Para comenzar tomaremos Herramientas de la
Microgeodesia, entendiendo que la microgeodesia involucra los elementos Geodésicos que
satisfacen fines Topográficos. Desde este punto de vista la Geodesia Geométrica es quien aporta
la mayoría de los criterios empleados en Proyectos de Grandes Obras Civiles. Sin embargo, con el
correr de los últimos años vemos en la Geodesia espacial una inagotable fuente Teórica que
contribuye a través de técnicas como los sistemas de posicionamiento Global (GPS), al desarrollo
de aplicaciones efectivas y eficientes. En este contexto el método de Bursa Wolf para
transformación de Datum, y tratamiento espacial de los sistemas de coordenadas Cartesianas nos
entrega a través de modo simplificado la forma que tendrá la ecuación de proyecto.
Por otra parte se debe recalcar, que en Proyectos de Ingeniería, sea en la Explotación Minera, o
bien en las concesiones Viales, hoy en día se encuentran nuevas dificultades, que las normas
vigentes no alcanzan a actualizar o bien los estándares cambian según lo especifico de la
herramienta o la tarea. Sin embargo, en la búsqueda del "aseguramiento de la calidad" es donde los
sistemas están implementando nuevos métodos de control y nuevas técnicas. El objetivo del
presente documento es relacionar dos conocidas herramientas que permiten la transformación de
sistemas de coordenadas para aplicaciones en proyectos de Ingeniería. Para llevar a cabo este fin,
hemos denominado como Método Proyectivo a la Transformación de Coordenadas provenientes de
la Proyección Cartográficas Universal Transversa de Mercator (UTM), ampliamente utilizada en
Chile. Estas coordenadas son Transformadas en Coordenadas en Planos Topográficos Locales
(PTL). El análisis de esta transformación al ser comparada estadísticamente con “la ecuación de
proyecto”, será la hipótesis que nos permitirá validar los parámetros del polinomio en cuestión, para
cada aplicación o zona de trabajo.
El método de transformación de Datum conocido por su autor Bursa Wolf, nos permite realizar la
transformación entre dos sistemas de coordenadas considerando sus aspectos Geométricos y
aplicando los criterios de Traslación, rotación y factor de escala. Si extrapolamos su uso a la
transformación de sistemas planos como UTM a PTL, nos permite encontrar una expresión para
cada proyecto, pudiendo realizar un análisis estadístico para la calidad de la transformación y la
definición de los parámetros de dicha expresión. Es aplicable a sistemas UTM que provienen de
distinto Datum Geodésico, compatibilizando de modo empírico y real los sistemas de Coordenadas.
Finalmente se puede establecer un sistema único para cada proyecto, siendo de vital importancia
pues evita la burocracia digital, esto es, facilita el traspaso de información entre diversas
aplicaciones. A continuación presentaremos los aspectos teóricos que nos permiten construir la
ecuación del proyecto y su análisis respectivamente.
Transformación de Datum. Bursa Wolf
Se considera un sistema geodésico con coordenadas rectangulares (x, y, z.) en donde el origen
esta en el centro del elipsoide. Los parámetros de este elipsoide serán el radio ecuatorial a y el
achatamiento f. Conociendo x, y, z, a, y f podemos calcular las coordenadas geodésicas latitud,
longitud y altura de los puntos en el sistema viejo , luego podemos considerar un nuevo sistema de
coordenadas rectangulares x, y, z. Cuyo origen puede estar en el centro de masa de la tierra
conociendo los parámetros a y f de este nuevo elipsoide de referencia , podemos calcular las
coordenadas geodésicas de los puntos en el nuevo sistema esto sistemas o pueden definirse de la
siguiente forma.
Debemos tener presente que se deben transformar las coordenadas geodésicas en coordenadas
cartesianas:
X = (N+h) * Cos(ϕ) * Cos(λ)
Y = (N+h) * Cos(ϕ) * Sen(λ)
Z = (N*(1-e)+h) * Cos(ϕ)
Luego X,Y,Z del datum origen se transformaran a X2,Y2,Z2 en el Datums
destino
X2 = Dx Cos(X,X2) Cos(Y,X2) Cos(Z,X2) X
Y2 = Dy + Cos(X,Y2) Cos(Y,Y2) Cos(Z,Y2) * Y
Z2 = Dz Cos(X,Z2) Cos(Y,Z2) Cos(Z,Z2) Z
Primero Consideramos las rotaciones considerando los tres ángulos de Euler. ω, ξ,
ψ.
Los efectos de rotación involucran: a los ejes del sistema viejo girando o
inclinándose con respecto a los ejes del sistema nuevo
Cos ω Sen ω 0
R3(ω)= -Sen ω Cos ω 0
0 0 1
1 0 0
R3(ω)= 0 Cos ξ Sen ξ
0 -Sen ξ Cos ξ
Cos ψ 0 -Sen ψ
R3(ω)= 0 1 0
Sen ψ 0 Cos ψ
Estas matrices se multiplican para obtener
R2(ψ) * R1(ξ) * R3(ω) =Cos ωCos ψ - Sen ω Sen ξ Sen ψ Sen ω Cos ψ + Cos
ω Sen ξ Sen ψ - Cos ξ Sen ψ-Sen ω Cos ξ Cos ω Cos ξ Sen ξ Cos ωSen ψ - Sen
ω Sen ξ Cos ψ Sen ω Sen ψ - Cos ω Sen ξ Cos ψ - Cos ξ Cos ψ
Escribiendo las ecuaciones individuales
X2 = Dx + X + ωY - ψ-z + X∆L
Y2 = Dy + Y + ωX - ξ-z + X∆L
Z2 = Dz + Z + ψZ - ξ-z + X∆L
∆L = Factor de escala
∆L = (S2 - S1) / S
DX, DY, DZ Son desvíos de traslación entre el sistema nuevo y el sistema viejo
teniendo como origen el sistema nuevo.
Como obtener los ángulos de rotación:
Se considerarán dos puntos i y j.
Donde Sik = distancia de cuerda entre las unidades i y k
Ahora se deberán calcular los cosenos direccionales correspondientes para la
misma línea basándonos en las coordenadas X,Y,Z estos valores se despejan como
aik, bik, cik. De donde mediante una conveniente aproximación se obtienen las
expresiones finales para obtener las coordenadas en el nuevo sistema de
coordenadas:
X ∆x 1 ω −ψ x
Y = ∆y + − ω 1 ε * y
Z ∆z ψ − ε 1 z
Obteniendo las relaciones individuales:
X = ∆x + x + ωy −ψz + x∆L
Y = ∆y + y − ωx + εz + y∆L
Z = ∆z + z + ψx − εy + z ∆L
Este sistema de relaciones individuales Finalmente nos proporciona una ecuación de la forma:
X = T + (+dL) Rx
En la cual podemos operar a partir de coordenadas Nuevas para obtener coordenadas antiguas o
bien emplear el arreglo de obtener los parámetros a partir de set de coordenadas conocidas (nuevas
y antiguas) a través de un ajuste estadístico con datos empíricos.
Esta forma nos permitirá construir una ecuación que, ejemplo puede tener la forma de:
NORTE UTM NUEVO = NORTEUTMANTIGUO*FACTOR1 +ESTEUTMANTIGUO*FACTOR2 +FACTOR3
ESTE UTM NUEVO = ESTEUTMANTIGUO*FACTOR4 +NORTEUTMANTIGUO*FACTOR5 +FACTOR5
Donde los factores serán validos para una de terminada Zona, esto puede ser comprobado a partir
de la transformación de Coordenadas UTM a Coordenadas en un plano Topográfico Local.
Transformación UTM a Topográficas PTL
La transformación de poligonales UTM (medidas con GPS generalmente) en poligonales
topográficas es un proceso donde cada distancia plana definida en el sistema UTM se transforma a
través de sucesivas operaciones en una distancia en el plano topográfico. Este plano topográfico
siempre será definido por su altitud. En la aplicación UTM esta altitud se llama altitud del plano de
referencia o altura de proyecto. Los ángulos formados por dos lados consecutivos de una poligonal
se ven afectados por las reducciones esféricas en la ocasión del cálculo de la proyección UTM. De
esta manera, en la transformación para el plano topográfico, estas reducciones afectan la dirección
de los lados de la poligonal, es decir, recaen sobre los azimut de cada lado de la poligonal. Los
azimut y las distancias que así se obtienen forman la poligonal topográfica, que mantiene una
proporcionalidad con la poligonal UTM. Para iniciar el proceso de transformación se deben elegir
de antemano dos puntos consecutivos: uno para origen de las coordenadas transformadas y otro
para servir de referencia (Azimutal) de las transformaciones. La referencia debe ser un punto
inmediatamente anterior al origen, considerando el sentido del recorrido de la poligonal.
E2
E1
E3
E0 E4
TOPOGRÁFICO
E2
E1
E3
E0 E4
UTM
Figura 1.0.- Plano UTM y Plano Topográfico Local PTL
El origen para la transformación no es necesario que sea la estación de partida de la poligonal
UTM. Se puede elegir cualquier estación de la poligonal, entretanto se recomienda observar los
siguientes criterios:
• Las Poligonales son poligonales teóricas con recorridos teóricos, basta con conocer las
coordenadas de los puntos, posteriormente se relacionan y se les asigna un recorrido, el cual puede
estar asociado con la realidad geoespacial del proyecto.
• No elegir la estación de partida o de llegada de la poligonal UTM para servir de origen de
las transformaciones, a menos que estas estaciones formen parte del recorrido normal de la
poligonal. Generalmente las referencias se ignoran en las transformaciones, como veremos más
adelante.
• Si el area de interés de la poligonal topográfica se restringe a determinado tramo, elija para
origen de las transformaciones la estación que esté más próxima al tramo. Las eventuales
distorsiones serán mínimas en las áreas próximas al origen.
Con las coordenadas de origen definidas, las distancias y azimutes son sucesivamente
transformados, generando las nuevas coordenadas hasta que la transformación llegue al punto de
origen.
A continuación se presentan distintos casos que se pueden presentar en poligonales.
POLIGONAL CERRADA
Para este tipo de poligonal no hay modificación de su conformación:
ORIGEN
E1
E0
REFERENCIA
E4 E3
TOPOGRÁFICO
E1
E0
REFERENCIA
E4 E3
UTM
Figura 2.0.- Poligonal Cerrada Teórica para Transformación.
POLIGONAL VINCULADA
Para este tipo se creará analíticamente un lado extra para la poligonal transformada, uniendo
siempre la última estación a la estación de partida.
VISUAL CALCULADA
E4
E3 E2
TOPOGRÁFICA
E0
LLEGADA
E4 PARTIDA
REFERENCIA
E3 E2
UTM
Figura 3.0.- Poligonal Vinculada Teórica para Transformación.
POLIGONAL ABIERTA
En este caso, la última estación de la poligonal será unida a la estación de partida por medios
analíticos. La poligonal transformada siempre tendrá la forma de una figura cerrada.
ORIGEN
VISUAL CALCULADA
REFERENCIA E0
E4
E3 E2
TOPOGRÁFICA
E0
LLEGADA
PARTIDA
E4
E3 E2
UTM
Figura 4.0.- Poligonal Abierta Teórica para Transformación.
Transformación por Tramos
Como hemos dicho en otro punto de este manual, la transformación se puede hacer en un único
plano de referencia o por tramos involucrando a varios planos escalonados.
Cada tramo tendrá definida una altitud de referencia propia y todos los planos se referirán a un
único origen de coordenadas.
Las figuras siguientes ayudan a aclarar el tema abordado.
Poligonales
A
E
B
D
C
PLANO 700.000
A
E
B
D
C
PLANO 600.000
A
E
B
D
C
PLANO 500.000
A
E
B
D
C
PLANO UTM
Figura 5.0.- Poligonales por tramos Teóricas para la Transformación.
La figura representa la poligonal UTM y la transformación para topográfica ejecutada en varios
planos, cada plano esta definido por la altura de proyecto como se aprecia en la figura.
La estación final de un tramo (Estaciones C y E en la figura) sirve de unión entre dos planos y
sus coordenadas valen para ambos planos. Como cada plano constituye por sí mismo, una
referencia, las coordenadas de un determinado plano no deben intervenir con las del otro plano,
excepto en el caso de que la diferencia entre los planos (diferencia de altitud) sea pequeña y que las
distorsiones ocasionadas del pasaje de un plano a otro sean despreciables.
Las distancias y ángulos que forman parte de un tramo, se transforman según el plano definido
para el tramo. Por lo tanto, las coordenadas transportadas desde el origen y que pasen por varios
planos, sufrirán distorsiones impuestas por la diferencia de altitud entre los planos, y por
consecuencia, al volver al punto de origen las coordenadas finales se verán afectadas por la suma
de las distorsiones.
Para estos casos, en el momento en que no haya vínculo entre las coordenadas, (es decir, que no
se estén transportando), no habrá elemento de unión entre los planos. Las coordenadas finales de
cada tramo se calculan también para el tramo siguiente, formando conjuntos distintos de
tramo/coordenadas para cada plano definido.
F G A
E
C
B
PLANO 700.000 A
F G
E
C
B
PLANO 600.000
G A
E F
C
B
PLANO 500.000
F G A REFERENC
E
ORIGEN
C
B
PLANO UTM
Figura 6.0.- Radiación Teórica para la Transformación.
Cada punto transportado queda vinculado al origen y a la referencia de la transformación, y se
puede intervenir libremente con estos puntos y también con los demás puntos que pertenezcan al
mismo plano. La interacción entre coordenadas de planos diferentes se debe restringir a los puntos
situados en planos que presenten diferencias de altitudes poco significativas.
Proceso de cálculos de UTM a Topográficas
El proceso de cálculo es conocido y se encuentra a vuestra disposición en
www.severino.s5.com/UTM-topo.html
Resultados
Se ha tomado una muestra particular para verificar la precisión de la “ecuación
de proyecto”, la cual comprende Latitud 25º y Longitud 69º en una extensión de
50 Kilómetros en Sentido Este Oeste y 50 Kilómetros en Sentido Norte Sur
Coordenadas UTM
Vértice Norte Este Cota
Pívot 7235015.017 500000.000 87.455
Calaje 7285015.017 500000.000 145.090
A3 7235015.017 525000.000 111.079
A4 7285015.017 525000.000 87.455
A5 7285015.017 475000.000 145.090
A6 7235015.017 475000.000 111.079
Tabla 1.0 Coordenadas UTM Muestra
Se adopto el criterio de Radiación en un Plano, para la generación del cálculo
UTM a PTL, obteniéndose los siguientes resultados:
Coordenadas Topográficas
Vértice Norte Este
Pívot 7235015.017 500000.000
Calaje 7285035.928 500000.000
A3 7235015.017 525010.391
A4 7285035.671 525010.649
A5 7285035.671 474989.351
A6 7235015.017 474989.609
Tabla 2.0 Coordenadas PTL calculadas
A partir de este calculo teórico, se realizo un ajuste empleando Criterios Bursa
Wolf, para determinar las deformaciones debido a Traslación, rotación y factor de
escala, obteniéndose los siguientes parámetros:
Parámetros Calculados
A 0.99958298070642100
B -0.00000000000086918
K1 208.50965309984400000
K2 3017.20106665319000000
Tabla 3.0 Factores de Conversión “Ecuación de Proyecto”
Si recalculamos las Coordenadas UTM considerando los factores de Conversión,
empleando una simple planilla de cálculo, se tienen los siguientes resultados:
Coordenadas UTM
Vertice Norte Este Norte Calculado Este Calculado Residuos Norte R e s i d u o s E s t e
Pivot 7235015.017 500000.000 7235015.078 500000.000 0.061 0.000
Calaje 7285015.017 500000.000 7285015.129 500000.000 0.112 0.000
A3 7235015.017 525000.000 7235015.078 524999.961 0.061 -0.039
A4 7285015.017 525000.000 7285014.871 525000.219 -0.146 0.219
A5 7285015.017 475000.000 7285014.871 474999.781 -0.146 -0.219
A6 7235015.017 475000.000 7235015.078 475000.039 0.061 0.039
Tabla 4.0 Comparación entre muestra de coordenadas UTM y Coordenadas
UTM obtenidas de la ecuación de Proyecto.
De acuerdo a este calculo se tiene una precisión de 1 : 250 000 y la dispersión de
los datos se puede apreciar en gráficos, según su variación en la coordenada Norte
y en la Coordenada este respectivamente.
Se realizo esta comprobación en diferentes Latitudes de Chile obteniéndose
similares resultados, estas se encuentran en Internet en www.severino.s5.com/utm-
topo.htm.
Conclusión
Se evidencia una precisión de 1 : 250 000 del calculo propuesto, siendo aplicable
por una simple planilla de calculo, abarcando una grilla de 50Km. Las proyecciones
de este tipo de cálculo permiten realizar transformación de coordenadas en
diferentes sistemas, se puede ampliar el rango de acción de este cálculo a
proyectos de mayor extensión pudiendo integrar fuentes de información.
Bibliografía
Richard H. Rapp, 1981. Geodesia Geométrica Departamento de Ciencias
Geodésicas, Universidad Estatal de Ohio.
Santiago y Cintra, 1987. Manual de Usuario Software Topograph, Brasil
Zakatov, 1970. Curso Geodesia Superior. Editorial MIR.
