Aparicio RV
Aparicio RV
Aparicio RV
TESIS
ELABORADO POR
ASESOR
LIMA- PERÚ
2018
© 2018, Universidad Nacional de Ingeniería. Todos los derechos reservados
“El autor autoriza a la UNI a reproducir de la Tesis en su totalidad o en parte, con fines
estrictamente académicos.”
Aparicio Ramírez, Vittorio Francesco
[email protected]
922767983
UNIVERSIDAD NACIONAL DE INGENIERÍA
DEDICATORIA
FACULTAD DE INGENIERÍA CIVIL
Dedicatoria
ÍNDICE
ABSTRACT __________________________________________________________ 3
PRÓLOGO __________________________________________________________ 4
RESUMEN:__________________________________________________________ 5
ÍNDICE DE FIGURAS __________________________________________________ 6
ÍNDICE DE TABLAS ___________________________________________________ 8
INTRODUCCIÓN _____________________________________________________ 9
CAPÍTULO I: GENERALIDADES_________________________________________ 10
1.1 ANTECEDENTES __________________________________________________ 10
1.2 JUSTIFICACIÓN ___________________________________________________ 12
1.3 PLANTEAMIENTO DEL PROBLEMA___________________________________ 12
1.4 DEFINICIÓN DE LOS OBJETIVOS _____________________________________ 12
1.4.1 Objetivo principal _____________________________________________________ 12
1.4.2 Objetivos específicos ___________________________________________________ 12
CAPÍTULO II: MARCO TEÓRICO ________________________________________ 13
2.1 CONCEPTOS FUNDAMENTALES _____________________________________ 13
2.1.1 Superficie topográfica __________________________________________________ 13
2.1.2 Geoide ______________________________________________________________ 13
2.1.3 Elipsoide de revolución _________________________________________________ 14
2.1.4 Altura Ortométrica (H) _________________________________________________ 15
2.1.5 Altura Elipsoidal (h) ____________________________________________________ 15
2.1.6 Ondulación Geoidal (N) _________________________________________________ 15
2.1.7 Sistema de referencia __________________________________________________ 16
2.1.8 Sistema astronómico global _____________________________________________ 16
2.1.9 Sistema elipsoidal global________________________________________________ 17
2.1.10 Proyecciones cartográficas ______________________________________________ 18
2.2 PROPIEDADES DE LAS PROYECCIONES CARTOGRÁFICAS ________________ 19
2.2.1 Proyección equidistante ________________________________________________ 19
2.2.2 Proyección conforme __________________________________________________ 19
2.2.3 Proyección equivalente_________________________________________________ 19
2.3 PROYECCIONES CILÍNDRICAS _______________________________________ 20
2.3.1 Proyección de MERCATOR ______________________________________________ 20
2.3.2 Proyección transversal de MERCATOR ____________________________________ 20
2.3.3 Proyección Universal Transversal de MERCATOR (UTM) ______________________ 21
2.3.4 FACTOR DE ESCALA (K ESCALA) __________________________________________ 21
2.4 ANÁLISIS DEL FACTOR DE ESCALA ___________________________________ 26
2.4.1 Radio de curvatura del meridiano en el punto “P” (ρ): Z GEODÉSICO = 0° ___________ 26
2.4.2 Radio de curvatura de la primera vertical en el punto “P” (GRAN NORMAL N): Z
GEODÉSICO=90° _______________________________________________________________ 26
2.4.3 Radio medio de curvatura (r) ____________________________________________ 27
2.4.4 Análisis del factor de escala en un punto (K escala) __________________________ 27
ABSTRACT
The evolution of geodesy, from it’s origins with Pythagoras and Aristotle, it has
sought the geometric shape that is closest to the surface of the Earth. Currently is
accepted that the shape closest to the earth is called the geoid and is defined as
the geometric place which is in equilibrium with the forces of attraction of the Earth,
other stars and the centripetal force generated by the rotation of the Earth.
Coordinates, angles and directions UTM, are projections from the geodetic surface
of reference surface and therefore does not accurately represent exactly the
reality, for this reason to reframe these coordinates it is necessary to go through a
transformation before that can represent the topographic coordinates starting from
UTM.
PRÓLOGO
Para este fin se optó por presentar las bases matemáticas y términos a utilizar,
posteriormente exponer la manera en que estas distintas variables se relacionan
y generan, luego de una transformación, nuevas coordenadas que podrán ser
utilizadas en campo; mostrando además todo el proceso de transformación y el
comportamiento de las variables durante dicho proceso.
Finalmente es válido mencionar que la presente Tesis tiene como marco teórico
los distintos términos y fórmulas matemáticas que han sido desarrollados durante
los últimos años al estudio de la Geodesia, tal es el caso de los aportes de
Eratóstenes de Alejandría (276-195 a.C), Jean Richter (1630-1696), Friedrich
Robert Helmert (1843-1917) y John Fillmore Hayford (1869-1925), en este sentido
es grato poder utilizar este conocimiento y aplicarlo al desarrollo del avance en el
campo de la geodesia en el Perú.
RESUMEN:
ÍNDICE DE FIGURAS
ÍNDICE DE TABLAS
INTRODUCCIÓN
Debido a que el método para la obtención de las coordenadas de los hitos, puntos
de control, se realiza mediante la utilización de un GPS diferencial de alta
precisión, estos puntos de control serán entregados en coordenadas UTM más no
en coordenadas topográficas, que son las utilizadas en campo.
CAPÍTULO I: GENERALIDADES
Personas de todo el mundo buscaron durante mucho tiempo el elipsoide que más
se aproximase a la forma de la Tierra, siendo esto imposible por las irregularidades
que ésta presenta. Pero gracias a los aportes del renombrado geodesta y
matemático alemán Friedrich Robert Helmert (1843 -1917) y su presentación de
los fundamentos para desarrollar un modelo de determinación del geoide se
produce la transición al concepto actual de la figura de la Tierra.
En 1909 el eminente geodesta John Fillmore Hayford (1869 -1925) publicó los
resultados para un elipsoide con DATUM en EEUU. En el año 1924 la asamblea
General de la Asociación Internacional de Geodesia, celebrada en Madrid, adopto
este elipsoide pasando a formar parte de las “Constantes Internacionales de
1927”.
El obtener diferentes resultados bajo los mismos datos iniciales, genera replanteos
erróneos que en muchos casos afectan predios vecinos al proyecto.
Pues un mal replanteo genera además de una traslación de los puntos del
proyecto, una rotación de los mismos.
2.1.2 Geoi de
Figura 2: Geoide.
Fuente: http://definicion.de/geoide/
Figura 3: Elipse.
Fuente: Jorge Mendoza Dueñas (2017). Topografías Técnicas Modernas, Página 568, Lima Perú, Editorial
Maraucano E.I.R.L
𝒂𝒂 − 𝒃𝒃 𝒂𝒂 − 𝒃𝒃
Achatamiento 𝒇𝒇 = 𝒇𝒇 = ( 1)
𝒂𝒂 𝒂𝒂
Fuente: Jorge del Río San José (2010) -Introducción al tratamiento de datos espaciales en hidrología.
Disponible en https://books.google.com.pe/books?id=87BeTf_cWS8C
Fuente: Jorge del Río San José (2010) -Introducción al tratamiento de datos espaciales en hidrología, figura
mejorada para el presente trabajo. Disponible en https://books.google.com.pe/books?id=87BeTf_cWS8C
Fuente: Jorge Mendoza Dueñas (2017). Topografías Técnicas Modernas, Página 592, Lima Perú, Editorial
Maraucano E.I.R.L
Fuente: Jorge Mendoza Dueñas (2017). Topografías Técnicas Modernas, Página 594, Lima Perú, Editorial
Maraucano E.I.R.L
Tiene la cualidad de mantener la distancia real entre dos puntos situados sobre la
superficie del Elipsoide.
Tiene la cualidad de conservar los ángulos formados por dos líneas, tanto en el
elipsoide como en el plano cartográfico; sin embargo es importante puntualizar
que no existe ninguna proyección conforme que mantenga dicha propiedad en
todo el elipsoide. Este tipo de proyecciones conserva la forma de las figuras, pero
no el tamaño de éstas.
Fuente: Jorge Mendoza Dueñas (2017). Topografías Técnicas Modernas, Página 603, Lima Perú, Editorial
Maraucano E.I.R.L
Fuente: Sánchez Menéndez, F. (2009). Geodesia Y Cartografía: Los Conceptos y su Aplicación Práctica.
Disponible en https://books.google.com.pe/books?id=dlUurxOj4A8C
Fuente: Sánchez Menéndez, F. (2009). Geodesia Y Cartografía: Los Conceptos y su Aplicación Práctica.
Disponible en https://books.google.com.pe/books?id=dlUurxOj4A8C
Es aquel valor que permite proyectar la longitud medida entre dos puntos en el
elipsoide de referencia sobre el plano cartográfico.
Lp = KEscalaxLo (4)
Fuente: Fernández Coppel, I. (2013). Las Coordenadas Geográficas y las Proyecciones UTM, Universidad de
Valladolid, España.
Esta proyección tiene valides entre las latitudes 84° Norte y 80° Sur. En áreas
polares no es conveniente el uso de estas coordenadas, debido a que a medida
que los paralelos se van separando las latitudes crecen.
Fuente: Daza Vásquez, L. (2004). Informe de Suficiencia: Aplicaciones de Navegadores GPS en la Ingeniería,
Universidad Nacional de Ingeniería, Perú.
Fuente: Daza Vásquez, L. (2004). Informe de Suficiencia: Aplicaciones de Navegadores GPS en la Ingeniería,
Universidad Nacional de Ingeniería, Perú.
Así pues, la latitud geodésica Φ, es el ángulo limitado por la normal ρ con el plano
ecuatorial.
a(1 − e2 )
ρ=
(1 − e2 . sen2 ɸ)1.5 (5)
a
N=
(1 − e2 . sen2 ɸ)0.5 (6)
Fuente: Mendoza Dueñas, J. (2013). Propiedades de las Proyecciones Cartográficas, Curso de Geodesia
Satelital, Universidad Nacional de Ingeniería, Perú.
ds
k=
ds´ (8)
ds´ = �dx 2 + dy 2
(9)
Dónde:
• ρ: Radio de curvatura del meridiano en el punto A.
• N: Radio de la gran normal en “A”
• ɸ: Latitud geodésica en A
Dónde:
q = 0.000001 X
(12)
1 + e´2 . cos 2 ɸ
P=� �. 1012
2. N 2 . k 0 2 (14)
Figura 17: Cálculo del factor de escala con el programa Geo uní 2015.
Fuente: Elaboración Propia
La Figura
2.5 17 muestra
TRAN SFORla interfaz
MAC del IÓN
programa
D GeoUni2015,
E C OORD el cual fue elaborado
ENADA S y diseñado por Vittorio Aparicio Ramírez con la
asesoría del Ing. Jorge Mendoza Dueñas, con fin académico para la presente Tesis y posteriores. Fue calibrado con tablas en Excel
Para el presente proyecto se usará el Datum WGS84 como referencia para todos
los cálculos.
La Figura 18, muestra las principales relaciones a emplearse durante los cálculos
y la Tabla 2 los principales parámetros del Elipsoide WGS 84 así como las tablas
de la 3 a la 5 muestran los parámetros a utilizar en la transformación de
coordenadas UTM a Geodésicas.
Fuente: Jorge del Río San José (2010) -Introducción al tratamiento de datos espaciales en hidrología.
Disponible en https://books.google.com.pe/books?id=87BeTf_cWS8C
Fuente: Mendoza Dueñas, J. (2013). Propiedades de las Proyecciones Cartográficas, Curso de Geodesia
Satelital, Universidad Nacional de Ingeniería, Perú.
𝐚𝐚
𝐍𝐍𝟏𝟏 = 𝐂𝐂𝟏𝟏 = 𝐞𝐞´𝟐𝟐 . 𝐜𝐜𝐜𝐜𝐜𝐜 𝟐𝟐 𝛟𝛟𝟏𝟏
(15) (16)
�𝟏𝟏 − 𝐞𝐞𝟐𝟐 . 𝐬𝐬𝐬𝐬𝐬𝐬𝟐𝟐 𝛟𝛟𝟏𝟏
𝐚𝐚 Este − 500 000
𝐑𝐑 𝟏𝟏 = (17) D1 = (18)
�𝟏𝟏 − 𝐞𝐞𝟐𝟐 . 𝐬𝐬𝐬𝐬𝐬𝐬𝟐𝟐 𝛟𝛟𝟏𝟏 0.9996
𝟏𝟏 − (𝟏𝟏 − 𝐞𝐞𝟐𝟐 )𝟎𝟎.𝟓𝟓 Norte − 10 000 000
𝐞𝐞𝟏𝟏 = (19) M1 = (20)
𝟏𝟏 + (𝟏𝟏 − 𝐞𝐞𝟐𝟐 )𝟎𝟎.𝟓𝟓 0.9996
M1
µ= (22)
e2 3. e4 5. e6
𝐓𝐓𝟏𝟏 = 𝐭𝐭𝐭𝐭 𝟐𝟐 𝛟𝛟𝟏𝟏 (21) a(1 − −
4 64
−
256
)
donde 𝑒𝑒2 es el cuadrado de la primera excentricidad
Fuente: Mendoza Dueñas, J. (2013). Transformación de Coordenadas UTM a Geodésicas, Curso de Geodesia
Satelital, Universidad Nacional de Ingeniería, Perú.
𝐃𝐃𝟒𝟒 N1 + tg ϕ1
𝐐𝐐 = −�𝟓𝟓 + 𝟑𝟑. 𝐓𝐓𝟏𝟏 + 𝟏𝟏𝟏𝟏. 𝐂𝐂𝟏𝟏 − 𝟒𝟒. 𝐂𝐂𝟏𝟏𝟐𝟐 − 𝟗𝟗. 𝐞𝐞´𝟐𝟐 �. (26) Δϕ = − (P + Q + S) (27)
𝟐𝟐𝟐𝟐 R1
Fuente: Mendoza Dueñas, J. (2013). Transformación de Coordenadas UTM a Geodésicas, Curso de Geodesia
Satelital, Universidad Nacional de Ingeniería, Perú.
Cálculo de la Longitud λ
Cálculo de Δ λ (Radianes)
𝐃𝐃𝟑𝟑 𝐃𝐃𝟓𝟓
𝐉𝐉𝐉𝐉 = 𝐃𝐃𝟏𝟏 − (𝟏𝟏 + 𝟐𝟐 𝐓𝐓𝟏𝟏 + 𝐂𝐂𝟏𝟏 ). (30) 𝐗𝐗𝐗𝐗 = �𝟓𝟓 − 𝟐𝟐. 𝐂𝐂𝟏𝟏 + 𝟐𝟐𝟐𝟐.𝐓𝐓𝟏𝟏 − 𝟑𝟑. 𝐂𝐂𝟏𝟏 𝟐𝟐 + 𝟖𝟖. 𝐞𝐞´𝟐𝟐 + 𝟐𝟐𝟐𝟐. 𝐓𝐓𝟏𝟏 �. (31)
𝟔𝟔 𝟏𝟏𝟏𝟏𝟏𝟏
𝐉𝐉𝐉𝐉 + 𝐗𝐗𝐗𝐗
𝚫𝚫𝚫𝚫 = � � (32)
𝐜𝐜𝐜𝐜𝐜𝐜 ɸ𝟏𝟏
Fuente: Mendoza Dueñas, J. (2013). Transformación de Coordenadas UTM a Geodésicas, Curso de Geodesia
Satelital, Universidad Nacional de Ingeniería, Perú.
𝛥𝛥λ. 180°
λ = λ0 + � � (33)
𝜋𝜋
Ejemplo:
Cálculo de la Zona
Sea: P=Zona.
λ
P = Truncar � + 31� (34)
6
Cálculo de Δ λ
Δ λ = λ − λ0 (36)
Cálculo de E (Este)
Tabla 6: Fórmulas para el cálculo del Este a partir de las coordenadas geodésicas.
𝐚𝐚
𝐭𝐭 = 𝐭𝐭𝐭𝐭ɸ (37) 𝐍𝐍 = = 𝐑𝐑𝐑𝐑𝐑𝐑𝐑𝐑𝐑𝐑 𝐝𝐝𝐝𝐝 𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂𝐂 𝐍𝐍𝐍𝐍𝐍𝐍𝐍𝐍𝐍𝐍𝐍𝐍 (38)
�𝟏𝟏 − 𝐞𝐞𝟐𝟐 . 𝐬𝐬𝐬𝐬𝐬𝐬𝟐𝟐 ɸ
(Δ λ. cos ɸ)3 . N(1 − t 2 + n2 ) (Δ λ. cos ɸ)5 .N(5 − 18. t 2 + t 4)
𝐧𝐧𝟐𝟐 = 𝐞𝐞´𝟐𝟐 . 𝐜𝐜𝐜𝐜𝐜𝐜 𝟐𝟐 ɸ (39) E´ = ( Δ λ.cos ɸ)N + + (40)
6 120
𝐄𝐄 = 𝟓𝟓𝟓𝟓𝟓𝟓 𝟎𝟎𝟎𝟎𝟎𝟎 + 𝟎𝟎. 𝟗𝟗𝟗𝟗𝟗𝟗 𝟔𝟔 . 𝐄𝐄´ (41)
Fuente: Mendoza Dueñas, J. (2013). Transformación de Coordenadas Geodésicas a UTM, Curso de Geodesia
Satelital, Universidad Nacional de Ingeniería, Perú.
Tabla 7: Fórmulas para el cálculo del Norte a partir de las coordenadas geodésicas.
𝐀𝐀𝐀𝐀 = 𝐚𝐚. (𝐀𝐀 𝟎𝟎. ɸ − 𝐀𝐀 𝟐𝟐 .𝐬𝐬𝐬𝐬𝐬𝐬 𝟐𝟐ɸ + 𝐀𝐀 𝟒𝟒. 𝐬𝐬𝐬𝐬𝐬𝐬 𝟒𝟒ɸ − 𝐀𝐀 𝟔𝟔 . 𝐬𝐬𝐬𝐬𝐬𝐬 𝟔𝟔ɸ) (46)
(𝚫𝚫 𝛌𝛌. 𝐜𝐜𝐜𝐜𝐜𝐜 ɸ)𝟐𝟐 . 𝐍𝐍. 𝐭𝐭 (𝚫𝚫 𝛌𝛌. 𝐜𝐜𝐜𝐜𝐜𝐜 ɸ)𝟒𝟒 . 𝐍𝐍. 𝐭𝐭. (𝟓𝟓 − 𝐭𝐭 𝟐𝟐 + 𝟗𝟗𝐧𝐧𝟐𝟐 + 𝟒𝟒𝐧𝐧𝟒𝟒) (𝚫𝚫 𝛌𝛌. 𝐜𝐜𝐜𝐜𝐜𝐜ɸ)𝟔𝟔 . 𝐍𝐍. 𝐭𝐭. (𝟔𝟔𝟔𝟔 − 𝟓𝟓𝟓𝟓𝐭𝐭 𝟐𝟐 + 𝐭𝐭 𝟒𝟒 )
𝐍𝐍´ = 𝐀𝐀𝐀𝐀 + + + (47)
𝟐𝟐 𝟐𝟐𝟐𝟐 𝟕𝟕𝟕𝟕𝟕𝟕
Ejemplo:
Ejemplo:
Dado que dicha longitud se desarrolla en un plano; su cálculo está gobernado por
la fórmula aplicada al plano cartesiano y – x, como lo muestra la Figura 23.
Dónde:
Para:
ka + kb
Kescala = (60)
2
Lc
L0 = (61)
Kescala
Es aquel valor que permite proyectar la longitud medida entre dos puntos en el
terreno (distancia reducida al horizonte) sobre el elipsoide de referencia.
Dónde:
Lo = Kelevacion x Lt (62)
R−M
Kelevacion = (63)
R+h
Ejemplo:
La dirección de una línea AB, está determinada por el ángulo horizontal (Ɵ) que
forma respecto a un sistema de coordenadas establecidas convencionalmente.
Comúnmente la dirección de una línea de referencia se determina mediante el
Azimut o Rumbo, donde el azimut es el Angulo horizontal formado por el Norte y
la línea de referencia, la Figura 29 ilustra el azimut entre dos puntos.
El M.G. de un punto del elipsoide de referencia, es la elipse que pasa por dicho
punto y por los polos norte y sur de dicho elipsoide, la Figura 30 muestra dicho
meridiano geográfico para un punto A cualquiera.
t = tgɸ (68)
Dónde:
(T − t)AB = −ΔN. (2X1 + X 2 ). P. 6.8755. 10−8 (70)
ΔN = NB − NA (71)
La Figura 36 muestra el azimut geográfico de una línea AB, el cual se calcula del
siguiente modo:
𝑍𝑍𝑍𝑍 = 𝑇𝑇 + ɣ (75)
Ejemplo:
Figura 37: Cálculo del Azimut de Cuadricula, Azimut Proyectado y Azimut Verdadero.
Ejemplo:
Para los siguientes datos halle las coordenadas topográficas de B tomando como
coordenadas topográficas de A las coordenadas UTM de A.
Tabla 10: Coordenadas (UTM ZONA 18 SUR – WGS84) para el control del Circuito 1.
separación al vértice
Descripción Punto Norte Este h
de control inicial
Vértice Azimut Inicial AZ66 8523718,957 596920,182 4182,078 1033,327
Vértice control Inicial PB66 8523648,917 597951,133 4221,378 0
Vértice control Final AZ65 8521325,372 599720,167 4266,709 2920,333
Vértice Azimut Final PB65 8521193,196 600568,487 4221,822 3589,026
Tabla 11: Coordenadas (UTM ZONA 18 SUR – WGS84) para el control del Circuito 2.
separación al
Descripción Punto Norte Este h vértice
de control inicial
Vértice Azimut Inicial AZ65 8521325,372 599720,167 4266,709 858,555
Vértice control Inicial PB65 8521193,196 600568,487 4221,822 0
Vértice control Final AZ64 8519390,272 601183,306 4224,746 1904,872
Vértice Azimut Final PB64 8518539,939 601578,364 4201,037 2838,948
Tabla 12: Coordenadas (UTM ZONA 18 SUR – WGS84) para el control de Circuito 3.
separación al
Descripción Punto Norte Este h vértice
de control inicial
Vértice Azimut AZ61 8516036,123 608235,972 4155,914
456,094
Inicial
Vértice control PB61 8516015,177 608691,585 4174,951
0
Inicial
Vértice control Final AZ60 8516254,893 611399,967 4136,368 2718,97
Vértice Azimut Final PB60 8516977,283 611928,908 4123,136 3377,263
Para este caso el circuito fue elaborado mediante trabajo de campo dentro del
campus de la Universidad Nacional de Ingeniería, se utilizaron los siguientes
puntos de control para el replanteo de las coordenadas Topográficas a UTM.
Se aclara que éstas son coordenadas UTM, la Figura 42 muestra las coordenadas
utilizadas en el trabajo, las mismas que se encuentran en la Tabla 13.
Tabla 13: Coordenadas (UTM ZONA 18 SUR – WGS84) para el control de Circuito 4 .
separación
al vértice
Descripción Punto Este Norte h
de control
inicial
Vértice Azimut Inicial P 277,083.711 8,670,035.367 104.739 45,989
Vértice control Inicial A 277,047.761 8,670,006.686 114.478 0
Vértice control Final H 276,852.569 8,670,664.165 105.493 685,841
Vértice Azimut Final Q 276,852.887 8,670,856.731 102.211 872,097
Los pasos para el replanteo son relativamente fáciles de seguir y pueden ser
programados en cualquier lenguaje computacional.
Una vez obtenidos los valores latitud, longitud y radio de gran escala en “A”, punto
en estudio, podemos calcular el factor de escala para la coordenada de inicio “A”,
posteriormente con los mismos pasos calculamos el factor de escala para el punto
“B”, coordenada final. Estos dos factores se promedian para tener un factor de
escala final “K de escala” el cual relaciona la longitud geodésica con la longitud de
cuadricula y posteriormente servirá para calcular la longitud topográfica.
En este sentido, el factor “K” será cercano a la unidad y por ende podemos
considerar que las coordenadas topográficas y UTM serán cercanas en su valor,
no obstante, esto no quiere decir que sean las mismas, por este motivo podemos
considerar este hecho como punto de partida para iniciar un proceso iterativo con
el fin de encontrar las coordenadas UTM que darán origen a las coordenadas
topográficas.
Se comparan las coordenadas del punto “B´” con las coordenadas que le dieron
origen “B”, la diferencia generada es restada a las coordenadas UTM del punto de
origen “B” obteniendo así un nuevo punto de origen “B´´” con sus respectivas
coordenadas UTM, se verifica que la diferencia entre B´ y B sea menor a 0.001,
caso contrario a partir de B´´ se generan nuevos valores de B´ y con este un nuevo
valor para B´´, restando la diferencia entre B´ y B a las coordenadas B´´, se vuelve
a verificar que la diferencia entre estas coordenadas B´ y B sea menor a 0.001,
de esta manera se genera un ciclo iterativo que termina cuando las coordenadas
del punto “B´´” den como resultado las coordenadas del punto “B”, las cuáles son
las coordenadas topográficas dato, coordenadas iniciales.
En tal sentido, para la primera iteración será necesario restar al valor de “a.r” el
valor “a”, valor de la separación inicial entre B´´ y UTM B, obteniendo como
resultado una separación entre UTM B y B´´2 de “a. (r-1)”, luego se llamará “f” al
incremento de la variable “r” entre los puntos de la siguiente iteración, en tal razón
ahora tendremos una distancia de separación para el punto topográfico de B y el
nuevo B´3 de “a.(r-1).(r+f)”, este valor será restado de las coordenadas UTM que
le dieron origen B´´2, obteniendo así nuevas coordenadas UTM B´´3, donde la
separación de estas coordenadas y las UTM B será de “a.(r-1)2+a.(r-1).f”, luego
debido a que en la última iteración las coordenadas UTM B´´3 generadas distan
de A menos que las UTM B distan del mismo punto, B´´3 se encuentra a la
izquierda de UTM B, la distorsión habrá disminuido con respecto de la obtenida
entre UTM B y topográfica de B, llamaremos a esta disminución en la distorsión
“g”, en este sentido la separación entre puntos topográfico B y B´4 será de
"(𝐚𝐚. (𝐫𝐫 − 𝟏𝟏)𝟐𝟐 + 𝐚𝐚(𝐫𝐫 − 𝟏𝟏). 𝐟𝐟).(𝐫𝐫 − 𝐠𝐠)", si continuamos con el análisis encontraremos
que la nueva separación entre las coordenadas UTM será de "𝐚𝐚. (𝐫𝐫 − 𝟏𝟏)𝟑𝟑 +
𝐚𝐚. (𝐫𝐫 − 𝟏𝟏)𝟐𝟐 . (𝐟𝐟 − 𝐠𝐠) + 𝐚𝐚. (𝒓𝒓 − 𝟏𝟏). 𝐟𝐟. 𝐠𝐠".
Figura 46: Comportamiento del proceso iterativo utilizado en el programa informático analizando las diferencias entre las coordenadas UTM y
topográficas del punto B con un K < 1.
Figura 48: Comportamiento del proceso iterativo utilizado en el programa informático analizando las diferencias entre las coordenadas
UTM y topográficas del punto B con un K > 1.
Por este motivo, en el presente trabajo se utilizarán tres métodos para la obtención
de las coordenadas UTM, el primer método es el descrito anteriormente en donde
mediante iteraciones repetidas se logra obtener la coordenada UTM de origen, el
segundo método consiste en calcular el “K” promedio generado entre los puntos
de control de la poligonal abierta utilizada para la corrección de los mismos, los
cuales son necesarios en todo trabajo de topografía, y posteriormente mediante
las distancias topográficas y este “K promedio” apoyándonos en que en puntos
cercanos, distancias de separación menores a 3 km, los valores del factor
combinado que transforma las longitudes topográficas a longitudes de cuadricula
serán similares para dichos puntos y que la corrección por curvatura es mínima,
el tercer método consiste en Calcular las coordenadas UTM mediante un proceso
de iteraciones repetidas y posteriormente llevar estas coordenadas a un sistema
de corrección mediante poligonal abierta.
Por otro lado, vale la pena aclarar que aunque el programa creado durante la
presente tesis tiene la opción de considerar en los cálculos la convergencia de
meridianos como variable dentro de sus líneas de código, para el presente trabajo
se consideró no utilizarla, dado que esta variable representa la diferencia entre el
norte de cuadricula y el norte geográfico, y debido a que el presente trabajó de
investigación se orientó con respecto a las coordenadas UTM, creando una
convergencia de meridianos nula, no es necesaria su utilización; es más, al optar
por considerarla estaremos girando el sistema de coordenadas un valor igual a
esta variable; para poder comprobar que el programa utilizado genera valores
coherentes incluso considerando la convergencia de meridianos, sería necesario
correr el programa con puntos orientados tomando como referencia al norte
geográfico, lo cual deja abierta la posibilidad de esta comprobación para futuros
trabajos de investigación.
Por otro lado, el error angular es pequeño, esto debido a que se está trabajando
sobre una proyección conforme, por lo que los ángulos topográficas y UTM serán
muy similares, mas no serán iguales pues al trabajar con las coordenadas UTM
como si fueran topográficas estamos obviando parámetros de corrección como la
corrección por curvatura y la convergencia de meridianos que generan un ligero,
pero importante cambio en los cálculos finales.
Vertice Angulo Horizontal Lado Distancia Dis. Acumulad Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
1,582 -2,127
Admisibles
Error Angular -0,00195246 0,00734931 OK Pto Este Norte Linea Azimut
Error Lineal 2,65062604 AZ66 596920,182 8523718,957 AZ66--PB66 93,887
Error Relativo 1/1223,128 10000 Falla PB66 597951,133 8523648,917 AZ65--PB65 278,856
AZ65 599720,167 8521325,372
PB65 600568,487 8521193,196
Figura 50: Ajuste de la Poligonal abierta del Circuito 1 sin una transformación previa.
Vertice Angulo Horizontal Lado Distancia Dis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
-0,080 0,089
Admisibles
Error Angular -0,001533727 0,00734931 OK Pto Este Norte Linea Azimut
Error Lineal 0,119814844 AZ66 596918,9583 8523715,26 AZ66--PB66 93,678
Error Relativo 1/27058,876 10000 OK PB66 597951,133 8523648,917 AZ65--PB65 278,647
AZ65 599730,3077 8521329,634
PB65 600579,8967 8521200,438
Figura 51: Ajuste de la Poligonal abierta del Circuito 1 con una transformación previa.
Circuito 2
Circuito 3
Vertice Angulo Horizontal Lado Distancia Dis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
0,058 0,032
Admisibles
Error Angular -0,00061275 0,0087841 OK Pto Este Norte Linea Azimut
Error Lineal 0,065898078 AZ65 599718,8745 8521322,312 AZ65--PB65 98,641
Error Relativo 1/40913,94 10000 OK PB65 600568,487 8521193,196 AZ64--PB64 334,866
AZ64 601190,6435 8519390,895
PB64 601589,2551 8518541,256
Figura 52: Ajuste de la Poligonal abierta del Circuito 2 con una transformación previa.
Fuente: Elaboración Propia
Vertice Angulo Horizontal Lado Distancia Dis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
0,026 -0,053
Admisibles
Error Angular -0,00101244 0,00680414 OK Pto Este Norte Linea Azimut
Error Lineal 0,058701887 AZ64 601179,7076 8519389,562 AZ64--PB64 154,863
Error Relativo 1/51408,381 10000 OK PB64 601578,364 8518539,939 AZ63--PB63 12,444
AZ63 604024,0049 8517367,113
PB63 603895,0403 8516782,684
Figura 53: Ajuste de la Poligonal abierta del Circuito 3 con una transformación previa.
Para este fin se optó por realizar un levantamiento topográfico con estación total
en una poligonal de 10 puntos (P, A, B, C, D, E, F, G, H, Q), la Tabla 17 muestra
las coordenadas obtenidas mediante un GPS diferencial de alta precisión mientras
la Tabla 18 muestra los puntos de control para este circuito, además de un
levantamiento de estos puntos con un GPS de alta precisión obteniendo así para
cada punto las coordenadas UTM provenientes de visaciones satelitales, este
trabajo se desarrolló dentro del campus de la Universidad Nacional de Ingeniería.
Tabla 17: Puntos GPS utilizados en el trabajo (UTM Zona 18 SUR – WGS84).
Tabla 21: Coordenadas UTM generadas a partir de las coordenadas topográficas dato.
K Promedio = 1.000198676
Vertice Angulo Horizontal Lado Distancia Distancia UTMDis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
0.019 -0.005
Admisibles
Error Angular 3.25'' 28.28'' OK Pto Este Norte Linea Azimut
Error Lineal 0.01995103 Az1 277083.711 8670035.37 Az1--PB1 231.417
Error Relativo 1/43901 10000 OK PB1 277047.761 8670006.69 PB8--Az2 180.095
PB8 276852.569 8670664.17
Az2 276852.887 8670856.73
Figura 55: Ajuste de la Poligonal abierta del Circuito 4 con una transformación previa de las longitudes.
Vale decir que aunque este método a diferencia del anterior es factible para la
gran mayoría de personas con algún conocimiento técnico del tema, no podemos
asegurar la precisión que tendrá este método en grandes longitudes con respecto
del primer método, pero aun así en distancias menores a 1 km se mantiene dentro
de los parámetros permitidos en trabajos topográficos; observamos además que
luego de corregidas las coordenadas obtenemos errores, en el mejor de los casos,
con respecto de las coordenadas dato del orden de 1/93 046, lo que quiere decir
que la aproximación una vez corregida la poligonal es muy alta.
Por otro lado, si hacemos una corrección de poligonal sin una transformación
previa obtendremos un error de 1/5928, el cual supera el límite máximo permitido
de 1/10 000 como lo podemos apreciar en la Figura 56.
Vertice Angulo Horizontal Lado Distancia Distancia UTMDis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
0.058 -0.136
Admisibles
Error Angular 3.25" 28.28" OK Pto Este Norte Linea Azimut
Error Lineal 0.14771362 Az1 277083.711 8670035.37 Az1--PB1 231.417
Error Relativo 1/5928.411 10000 Falla PB1 277047.761 8670006.69 PB8--Az2 180.095
PB8 276852.569 8670664.17
Az2 276852.887 8670856.73
Figura 56: Ajuste de la Poligonal abierta del Circuito 4 sin una transformación previa.
Tabla 27: Corrección de errores para datos calculados mediante iteraciones sucesivas.
Distancia Error
Norte Este Distancias
∆N ∆T Acumulada relativo
Punto Corregido Corregido Norte (m) Este (m) ∆E (m) (m)
(m) (D.Ac/∆T)
(m) (m)
(m) (m)
BIBLIOGRAFÍA
ANEXOS:
-0,080 0,089
Admisibles
Error Angular -0,001533727 0,00734931 OK Pto Este Norte Linea Azimut
Error Lineal 0,119814844 AZ66 596918,9583 8523715,26 AZ66--PB66 93,678
Error Relativo 1/27058,876 10000 OK PB66 597951,133 8523648,917 AZ65--PB65 278,647
AZ65 599730,3077 8521329,634
PB65 600579,8967 8521200,438
0,058 0,032
Admisibles
Error Angular -0,00061275 0,0087841 OK Pto Este Norte Linea Azimut
Error Lineal 0,065898078 AZ65 599718,8745 8521322,312 AZ65--PB65 98,641
Error Relativo 1/40913,94 10000 OK PB65 600568,487 8521193,196 AZ64--PB64 334,866
AZ64 601190,6435 8519390,895
PB64 601589,2551 8518541,256
Vertice Angulo Horizontal Lado Distancia Dis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
0,026 -0,053
Admisibles
Error Angular -0,00101244 0,00680414 OK Pto Este Norte Linea Azimut
Error Lineal 0,058701887 AZ64 601179,7076 8519389,562 AZ64--PB64 154,863
Error Relativo 1/51408,381 10000 OK PB64 601578,364 8518539,939 AZ63--PB63 12,444
AZ63 604024,0049 8517367,113
PB63 603895,0403 8516782,684
Vertice Angulo Horizontal Lado Distancia Distancia UTM Dis. Acumulada Azimut Azimut Corr. Este Norte Y corr. X corr. Este Corr. Norte Corr.
Grados Minutos Segundos
0.019 -0.005
Admisibles
Error Angular 3.25" 28.28" OK Pto Este Norte Linea Azimut
Error Lineal 0.0199509 P 277083.711 8670035.367 P--A 231.417
Error Relativo 1/43901 10000 OK A 277047.761 8670006.686 H--Q 180.095
H 276852.569 8670664.165
Q 276852.887 8670856.731