Parametros Sistema de Georreferenciacion CMT 2010
Parametros Sistema de Georreferenciacion CMT 2010
Parametros Sistema de Georreferenciacion CMT 2010
SATELITAL GPS
FOTOGRAFÍA
AÉREA
DIGITAL
ESTUDIOS
TOPOGRÁFICOS
w w w . c e l i s - s t . c o m
SISTEMA DE GEOREFERENCIACIÓN
CMT - 2010
PARAMETROS DE CÁLCULO
RECOMENDACIONES PARA
SANTIAGO DE CALI
SEPTIEMBRE DE 2010
1
PARAMETROS DE CALCULO Y
TRANSFORMACION DEL SISTEMA CMT
2
Geoide
Es la superficie teórica de la tierra que une todos los puntos que tienen igual
gravedad, esta superficie no es uniforme sino que presenta una serie de
irregularidades, causadas por la distinta composición mineral del interior de la
tierra y de sus distintas densidades, lo que implica que para cada punto de la
superficie terrestre exista una distancia distinta desde el centro de la tierra al punto
del geoide.
Figura 1.
Fuente de las coordenadas geográficas y la proyección UTM – Prof. Ignacio Alonzo Fernández-
Coppel – Universidad de Valladolid.
3
Elipsoide
Figura 2.
Fuente de las coordenadas geográficas y la proyección UTM – Prof. Ignacio Alonzo Fernández-
Coppel – Universidad de Valladolid
Figura 3.
4
Figura 4.
Fuente Tesis “Evaluación de la Alteración Superficial por Influencia de la Altura y Localización en
una Proyección Universal Transversal de Mercator UTM”
Marjorie Tapie Prieto
Juan Luis Ramírez Lara
Universidad Tecnológica Metropolitana - Chile
Figura 5.
Fuente de las coordenadas geográficas y la proyección UTM – Prof. Ignacio Alonzo Fernández-
Coppel – Universidad de Valladolid
5
Coordenadas Elipsoidales
“La longitud () es el ángulo, medido sobre el plano ecuatorial, entre el meridiano
de referencia (normalmente Greenwich) y el meridiano del punto de interés, es
positiva al este de Greenwich y negativa hacia el oeste. Su rango se define
mediante -180º + 180º o 180º W 180º E, lo que también
equivale a 0º 360º.”
Figura 6.
Fuente de las coordenadas geográficas y la proyección UTM – Prof. Ignacio Alonzo Fernández-
Coppel – Universidad de Valladolid.
6
Coordenadas Planas.
7
Proyección Cartográfica Gauss-Krüger
Esta es una representación conforme del elipsoide sobre un plano, es decir, que el
ángulo formado entre dos líneas sobre la superficie terrestre se mantiene al ser
éstas proyectadas sobre el plano. Los meridianos y paralelos se interceptan
perpendicularmente, pero no son líneas rectas, sino curvas complejas, excepto el
meridiano central (de tangencia) y el paralelo de referencia. La escala de
representación permanece constante sobre el meridiano central, pero ésta varía al
alejarse aquel, introduciendo deformaciones en función de la longitud ().
Latitud 04°35’46.32150” N
Longitud 77°04’39.02850” W
Falso Norte 1’000.000
Falso Este 1’000.000
Factor de Escala 1.00000
Figura 7.
Sistema de proyección cartográfica Gauss-Krüger
Fuente IGAC
8
Proyección Cartesiana
Latitud 03°26’30.77999” N
Longitud 76°31’14.02499” W
Falso Norte 872.364,63
Falso Este 1’061.900,18
Plano de Proyección 1.000,00
Factor de Escala 1.00000
Figura 8.
Fuente IGAC (Instituto Geográfico Agustín Codazzi)
9
Certificado Placa IGAC
10
Generalidades del Sistema CMT
El sistema Catastro Municipal Topografía (CMT) fue creado inicialmente como una
base de puntos de control para el ajuste del recientemente plano digital de Cali del
año 1993, este proyecto fue auspiciado por EMCALI, Planeación Municipal y
Catastro Municipal.
11
Con base en estas conclusiones, un grupo de profesionales conformado por:
12
Certificado Placa IGAC
13
Certificado Placa IGAC
14
AMARRE PROYECTO TERRAZAS A PUNTOS IGAC
15
Estos puntos se calculan y se ajustan obteniéndose los siguientes resultados.
Con estas coordenadas elipsoidales se realiza una nueva red hacia 08 puntos
CMT existentes y distribuidos en la ciudad, se escogen estos puntos por permitir
un adecuado rastreo de señales GPS, conllevando menores saltos de ciclo e
indicadores bajos de ruido.
16
AMARRE PROYECTO TERRAZAS A PUNTOS CMT
17
Los resultados obtenidos fueron los siguientes:
18
Ya obteniendo las coordenadas elipsoidales ajustadas a sistema Magna-Sirgas
IGAC, se inicia a desarrollar los cálculos para determinar el origen probable del
sistema CMT-2010, evaluándose el comportamiento del sistema en planas
cartesianas y planas Gauss Kruger para determinar los parámetros a seguir y
determinar los componentes del sistema nuevo.
Se determinaba cada punto CMT como posible origen del sistema plano
cartesiano y se evaluaban los desplazamientos a las coordenadas teóricas del
sistema CMT-1993 original, estas coordenadas planas están descritas en las
cartillas planeación y son las que siempre se determinaron inicialmente como las
coordenadas del punto CMT estudiado.
19
2- Análisis de resultados con origen probable en CMT 098.
20
5- Análisis de resultados con origen probable en CMT 062.
21
8- Análisis de resultados con origen probable en CMT 129.
22
Estos resultados nos permiten llegar a las siguientes conclusiones:
El punto que mejor se comporta como origen cartesiano del sistema CMT-
2010 es el punto denominado CMT 101, ya que sus desviaciones sin las de
mejor aproximación a las teóricas del sistema CMT original de 1993.
Todos los cálculos se han hecho con base en los parámetros de cálculo
como sistema plano cartesiano, se realizó un nuevo cálculo para deducir si
podría pertenecer a un sistema plano Gauss Krüger, y el resultado nos
permitió deducir que por sus desviaciones altas, el sistema original CMT-
1993 siempre fue concebido como sistema plano cartesiano, a continuación
se presenta el cuadro con los resultados Gauss Krüger obtenidos.
23
10- Análisis de resultados con origen probable en CMT 101, para sistema plano
Gauss Krüger.
11- Análisis de resultados con origen probable en CMT 101, con base en la
desviación media del error y altura del sistema según cota geométrica
teórica.
12- Análisis de resultados con origen probable en CMT 101, con base en la
desviación media del error y el promedio de altura elipsoidal.
24
13- Análisis de resultados con origen probable en CMT 101, con base en la
desviación según determinación de promedio determinado por métodos
robustos.
14- Análisis de resultados con origen probable en CMT 101, con base en la
desviación según determinación de promedio determinado por métodos
robustos.
25
Con base en todos estos análisis la conclusión final para determinar el nuevo
origen del sistema CMT-2010, y el cual permitirá desarrollar trabajos con GPS en
la ciudad de Cali permitiendo cálculos y traslaciones a los sistemas Magna-Sirgas
existentes, es el origen deducido en el cuadro No.11, es decir se respetan las
coordenadas elipsoidales calculadas inicialmente a partir de los puntos IGAC
certificados, pero se desplazan las planas cartesianas CMT-1993 para ajustar
mejor los desplazamientos a los puntos CMT evaluados en las campañas GPS,
los parámetros de cálculo definitivos se presentan a continuación;
Constante gravitacional
Constante definida GM = 3 896 005 x 108 m3 s‐2
geocéntrica (GM)
Factor de aplanamiento
Constante definida J2 = 108 263 x 10‐8
dinámico (J2)
Velocidad angular de rotación
Constante definida = 7 292 115 x 1011 rad s‐1
()
Cálculo iterativo a partir
Primera excentricidad (e2) e2 = 6,694 380 022 90 x 10‐3
de ɑ, GM , J2,
26
Plano de Proyección = 966.560
Términos Auxiliares
= o
m = (
NO = a / √1‐e2 sin2
N = a / √1‐e2 sin2p
Oa(1-e) / (1-esin2 )3/2
ma(1-e) / (1-esin2 )3/2
27
Coordenada Norte
rad((tan radcos(
[ m
= O rad((tan radcos(
[ mx[ m
Coordenada Este
radcos
[
= radcos[
28
Modelo Matemático para conversión de coordenadas planas Cartesianas
Sistema CMT-2010 a Elipsoidales (WGS84)
Datos para Cálculo
Latitud del punto de cálculo p
Longitud del punto de cálculo p
Términos Auxiliares
= N ‐ No
= E ‐ Eo
O = a / √1‐e2 sin2
= a / √1‐e2 sin2p
Oa(1-e) / (1-esin2 )3/2
Latitud del punto buscado p
a(1-e2))]
tan / (2 M N)
a)
= a(1-e2))]tan / (2 M N) x a)
p = o +
Longitud del punto buscado p
= / ( cos 1 + /a])
p = o +
29
Con el origen ya determinado, se proceden a calcular las coordenadas tanto
planas cartesianas sistema CMT-2010 y planas cartesianas Magna-Sirgas origen
Cali de los puntos componentes del proyecto base.
30
Todo este análisis previo se utilizó para determinar las coordenadas del proyecto
PAVICOL, en el cual se materializaron y posicionaron 20 puntos en campaña GPS
realizada el mes de agosto del presente año.
Los resultados de los cálculos arrojaron los siguientes resultados.
17- Cuadro de coordenadas elipsoidales puntos PAVICOL.
31
18- Cuadro resumen de coordenadas planas sistema CMT-2010 proyecto
PAVICOL.
32
Aproximaciones al Cálculo de Cota
Con base en todos los cálculos realizados a partir de alturas elipsoidales, se
analizó el comportamiento en altura vertical, esto con el fin de evaluar el valor la
ondulación que nos permitiera trabajar con altura amarrada al sistema CMT, los
detalles del análisis se resumen en el siguiente cuadro.
33
Del cuadro anterior se puede deducir que:
Debido a las discrepancias en cota del sistema CMT 1993, en las cuales en un
mismo circuito pueden encontrarse diferencias hasta de 0.20 centímetros en un
radio menor a 300 metros, estas se ven reflejadas en los valores altos de
desviaciones respecto a una media inicial calculada.
Se puede optar por un cálculo más riguroso como el método de estimador robusto
de hubert, pero ello no nos garantiza la estabilidad en la cota, ya que parte del
sistema CMT 1993, el cuál como hemos reiterado en el presente documento, tuvo
y tiene grandes problemas en precisión, tanto en coordenadas x,y como en cota.
34
20- Cuadro resumen de análisis previo a ajuste.
35
21- Matriz de ecuaciones de condición.
36
22- Cuadro resumen de circuitos ajustados.
37
Recomendaciones Finales
38
Se agradece a las siguientes firmas por haber participado en el proyecto.
39