Sistema Geodésico Nacional UTM

1 Sistema Geodésico Nacional UTMDirección General de Geog...
Author: Bayardo Esquivel
0 downloads 2 Views

1 Sistema Geodésico Nacional UTMDirección General de Geografía y Medio Ambiente Subdirección Estatal de Geografía Departamento de Geodesia Sistema Geodésico Nacional UTM Junio 2014

2 NTG-001 Sistema Geodésico NacionalContenido Introducción Conceptos Básicos Proyecciones Cartográficas UTM Practica con navegador NTG-001 Sistema Geodésico Nacional DGGMA

3 Introducción

4 Introducción Es de vital importancia, cuando se combinan datos de diferentes fuentes en un SIG, o bien cuando se van a capturar nuevos datos, tener presente la referencia geodésica de cada información, así como la proyección cartográfica utilizada en el caso de datos proyectados, como pueden ser los de un mapa en papel. Las Estaciones Geodésicas en sus vertiente horizontal, vertical y gravimétrica forman parte del Grupo de Datos del Marco de Referencia Geodésico del Subsistema Nacional de Información Geográfica y del Medio Ambiente; DGGMA

5 Introducción La referencia geodésica “elipsoide y datum” determina la relación entre las coordenadas de los datos y su localización en la superficie terrestre. Así, el desconocimiento o uso incorrecto de la referenciación geodésica hace que datos de diferentes fuentes no coincidan entre sí. Es típica una desviación de unos 300 m. Tal sería el caso de intentar combinar datos de GPS con mapas catastrales o mapas topográficos sin tener en cuenta las distintas referencias. Las Estaciones Geodésicas en sus vertiente horizontal, vertical y gravimétrica forman parte del Grupo de Datos del Marco de Referencia Geodésico del Subsistema Nacional de Información Geográfica y del Medio Ambiente;

6 Introducción Incluso si los datos tienen un sistema de referencia común, el desconocimiento de éste puede acarrear problemas, si, por ejemplo, se requiere convertir coordenadas geográficas a U.T.M o viceversa. Hay que destacar que las herramientas que los SIG proporcionan al usuario para tratar estos temas son desiguales según los casos. Las Estaciones Geodésicas en sus vertiente horizontal, vertical y gravimétrica forman parte del Grupo de Datos del Marco de Referencia Geodésico del Subsistema Nacional de Información Geográfica y del Medio Ambiente;

7 Conceptos Básicos

8 Conceptos Básicos GeoideEl geoide es una superficie equipotencial que representa el nivel medio del mar y sirve como referencia para determinar alturas ortométricas.

9 Conceptos Básicos La Geodesia tiene un doble propósito: Como ciencia, la determinación de la forma y dimensiones del globo terrestre y de su campo gravimétrico externo Como aplicación práctica, la determinación de la posición geográfica de puntos cerca o sobre la superficie terrestre, en apoyo principalmente a la generación de cartografía. Tiene también otras aplicaciones, como exploraciòn de recursos energéticos y mineros, y en cualquier aplicación que requiera ser referenciada geográficamente.

10 Conceptos Básicos Uno de los propósitos fundamentales de la Geodesia, es la ubicación espacial precisa de los objetos que se encuentran en, sobre o cerca de la superficie de la Tierra, por lo que se convierte en elemento básico, que garantizará la referencia de los datos y productos estadísticos y geográficos generados por las Unidades del Estado, que integran el SNIEG

11 Conceptos Básicos Incluso la superficie del mar es irregular.La Forma de la Tierra.- Nuestro planeta en realidad tiene una forma mucho más complicada de lo que se pudiese lograr con una representación matemática de alguna figura que se le asemeje. Incluso la superficie del mar es irregular.

12 Conceptos Básicos La elipse es el lugar geométrico de los puntos del plano tales que la suma de las distancias a dos puntos fijos llamados focos es una constante positiva. Una elipse es la curva cerrada que resulta al cortar la superficie de un cono por un plano oblicuo al eje de simetría –con ángulo mayor que el de la generatriz respecto del eje de revolución. Una elipse que gira alrededor de su eje menor genera un esferoide achatado, mientras que una elipse que gira alrededor de su eje principal genera un esferoide alargado.

13 Conceptos Básicos Elipsoide.-Es un sólido geométrico generado por la rotación de una elipse sobre uno de sus ejes. Rotación del eje “b” Eje Polar Rotación del eje “a” Eje Ecuatorial Elipsoide de Referencia.- El elipsoide empleado como la mejor aproximación local o global de la forma de la Tierra;

14 Conceptos Básicos La Forma de la Tierra.-El elipsoide es la mejor aproximación global o local de la forma de la Tierra, con un radio levemente más grande en su Ecuador (eje “a”) que en los Polos (eje “b”). Elipsoide Esferoide Entonces es posible representar la posición de cualquier punto sobre el planeta, con valores o magnitudes (coordenadas) asociados a una figura de revolución.

15 Conceptos Básicos Las Tres Superficies.-El Terreno o Superficie Terrestre es una superficie irregular y es el lugar físico donde se realizan diversas mediciones (estación, punto, vértice, marca, etcétera). El Geoide es una superficie equipotencial del campo gravitatorio terrestre, que se aproxima al nivel medio del mar en los océanos. El Elipsoide es la representación matemática de la superficie terrestre (global o local).

16 Conceptos Básicos Dátum.-Cualquier cantidad o conjunto de ellas que sirve como referencia para calcular otras. NAD27 Dátum Norteamericano … NAD83 Dátum Norteamericano … WGS84 Sistema Geodésico Mundial … ITRFxx Marco Internacional de Referencia Terrestre … (0, 0, 0) NAVD29 Dátum Vertical Norteamericano … NAVD88 Dátum Vertical Norteamericano … International Gravity Standardization Net IGSN71 Red Internacional de Estandarización de la Gravedad …

17 Conceptos Básicos La Tierra, El Dátum y El Elipsoide.- Dátum ElipsoideRadio Ecuatorial (a) Radio Polar (b) Radio Medio Terrestre 6378 Km 6357 Km 6371 Km (a-b) = 21 Km Dátum Elipsoide Semieje Mayor (a) Semieje Menor (b) NAD27 CLARKE 1866 m m NAD83 GRS80 m m WGS84 m ITRFxx

18 Conceptos Básicos Sistema de Referencia.-Consiste de una serie de prescripciones y convenciones junto con un modelo requerido para definir en cualquier momento un sistema de ejes coordenados.

19 Conceptos Básicos Marco Coordenado de Referencia.-Es la construcción física de cualquier sistema de referencia mediante el establecimiento de puntos sobre la superficie terrestre, directamente accesibles por ocupación u observación.

20 Conceptos Básicos Coordenadas.-Cualquier punto en la Tierra puede ser expresado o representado por magnitudes únicas como: Coordenadas Geodésicas (φ, λ, h), Coordenadas Cartesianas (X, Y, Z) o Coordenadas de Proyección (X-UTM, Y-UTM y H- NMM = N,E,O). Coordenadas Cartesianas Coordenadas Geodésicas Coordenadas UTM

21 Proyecciones Cartográficas

22 Proyecciones CartográficasExpresión matemática para producir todo o una parte de un cuerpo esférico como la Tierra sobre una superficie plana o una figura geométrica.

23 Proyecciones CartográficasTodo mapa está en un determinado sistema de proyección, que responde a la necesidad de representar en una forma sistemática la superficie terrestre, con sus detalles, sobre la superficie del mapa.

24 Proyecciones CartográficasGeométricamente podemos ver la deformación al proyectar la elipse en un plano. En la imagen es claro que b > a.

25 Proyecciones CartográficasIntentando reducir la deformación (esto es, a b), hacemos tangente el plano al elipsoide:

26 Proyecciones CartográficasOtra opción es hacer secante el plano al elipsoide:

27 Proyecciones CartográficasRequisitos de las proyecciones En términos generales, se requiere de una proyección que se satisfagan los siguientes requisitos: a) Mantenimiento de la escala (equidistancia), b) Preservación de las áreas (equivalencia), c) Conservación de las formas (ortomorfismo) d) Exactitud en las direcciones. Estas cualidades no pueden ser satisfechas simultáneamente, y así, toda proyección operativa es una solución de compromiso entre éstas. Su escogencia depende del propósito de uso de la proyección, según el tipo de mapa.

28 Proyecciones CartográficasCLASIFICACIÓN, POR TIPO DE SUPERFICIE Cónicas: Traslada la información de la esfera terrestre a un cono, tomando como punto focal uno de los polos. Cilíndricas: Usa un cilindro tangente o secante a la esfera terrestre, en su forma normal se coloca el cilindro de tal manera que el paralelo de contacto es el ecuador. Transversa Oblicua Normal Paralelos Estándar Meridiano Central La malla de meridianos y paralelos se dibuja proyectándolos sobre la superficie especificada suponiendo un foco de luz que se encuentra en el centro del globo. Cónica: Hay una distorsión asimétrica que también afecta a las zonas polares, pero ofrece mayor precisión en el hemisferio que corresponde al polo que se haya tomado como foco.

29 Proyecciones CartográficasLas proyecciones utilizadas en el INEGI son: Cónica Conforme de Lambert. Universal Transversa de Mercator. Transversa Modificada Ejidal. Proyección Cónica Proyección Cilíndrica.

30 Proyección Uniersal Transversa de Mercator

31 UTM Hacia finales de los cuarenta el servicio cartográfico de la armada de los Estados Unidos adopto el uso de la proyección Transversa de Mercator, con la variante de que se generaron franjas de 6 en longitud, 3 a cada lado del meridiano central. Este sistema cartográfico recibió el poco modesto nombre de Proyección Universal Transversa de Mercator o UTM. El sistema UTM esta compuesto por sesenta zonas numeradas, usando los modelos matemáticos de la proyección Gauss Kr¨uger. En la actualidad, el 80% de las áreas continentales esta cubierta con cartografía basada en los modelos elipsoidales de la UTM.

32 UTM La proyección UTM esta definida de la siguiente manera:Esta basada en la proyección Gauss–Kr¨uger (proyección conforme) El mundo esta cubierto por franjas de 6 en longitud, con un meridiano central a cada 3. Cada meridiano central tiene como Falso Este metros. En el meridiano central el factor de escala es constante e igual a 0,9996

33 UTM Sistema de representación de la superficie de la Tierra sobre un plano basado en una superficie cilíndrica que es secante y en una dirección perpendicular al eje de rotación terrestre. Divide a la tierra en 60 zonas de 6 grados de ancho cada una. Es una proyección cilíndrica transversa secante. La proyección de referencia es la Gauss-Krüger. Las zonas se numeran de la uno a la sesenta, cubriendo la zona uno la longitud 180 Oeste a 174 Oeste, incrementándose el numero de zona hacia el Oeste. La zona sesenta abarca de la longitud 174 Este a 180 Este. El Sistema cubre de la latitud 84 Norte a 80 Sur. 120 W 84 W

34 UTM Es una versión especializada de la Transversa Mercator.El globo está dividido en 60 fajas al Norte y al Sur de 6º de longitud cada una. Cada faja tiene su propio meridiano central Los límites de cada zona son 84º al N y 80º al S, ocurriendo la división entre N y S en el ecuador. El origen de cada zona es el meridiano central y el ecuador.

35 UTM Método de Proyección: cilindro transversoLíneas Estándar: dos líneas paralelas y a aproximadamente 180 km a ambos lados del meridiano central de la zona UTM.

36 UTM PROPIEDADES: FORMA: Conforme. Mantienen las pequeñas formas. Mínima distorsión de las formas grandes dentro de la zona. AREA: Mínima distorsión dentro de cada zona UTM. DIRECCION: Los ángulos locales son verdaderos DISTANCIA: La escala es constante a lo largo del meridiano central, pero a un factor de escala de 0,9996 se reduce la distorsión lateral dentro de cada zona.

37 UTM Esta proyección fue definida hacia finales de los cuarenta por el servicio cartográfico de la armada de los Estados Unidos, la cual tiene las siguientes especificaciones: Cada meridiano central tiene como Falso Este 500,000 metros. En el meridiano central el factor de escala es constante e igual a El origen de las ordenadas es el Ecuador con valor 0 metros para latitudes norte y 10’000,000 metros la latitudes sur. ZONA-14 500,000 200,000 800,000 102° W 96° W 99° W 10’000,000 m Incrementando hacia el Norte 0 m Disminuyendo hacia el Sur Ecuador REMARCAR que los datos en la proyección UTM van acompañados de la zona en la cual son generados. En el meridiano central, la convergencia es constante e igual a 0

38 UTM Representación de la Republica Mexicana.Ejemplo de tres vértices diferentes pero con valor en Norte y Este iguales. ZONA11 ZONA12 ZONA13 ZONA14 ZONA15 ZONA16 117 W 111 W 105 W 99 W 93 W 87 W Note que los tres vértices son diferentes en coordenadas geodésicas sin embargo son iguales al ser proyectados con la UTM pero ubicados en diferente zona.

39 UTM Representación de la Republica Mexicana considerando que se encuentra en una zona determinada. Este Norte Zona 16 Zona 11

40 Taller Ejecutivo del Sistema IRIS 4.010/04/2017 UTM La República Mexicana se ubica desde la zona 11 hasta la 16. Por cuestiones de organización INEGI emplea fajas horizontales que abarcan de la letra “D” a la “I” para el mosaico nacional. Cada cuadrante mide 6° por 4°. G13 Sistemas de Información Geográfica

41 UTM La cartografía topográfica oficial 1: 50,000 y 1: se edita en el sistema de proyección cartográfica UTM (Universal Transversa Mercator), siendo una proyección conforme, es decir que mantiene los ángulos y la semejanza de figuras superficiales comprendidas dentro de elementos infinitesimales de la superficie terrestre; esta propiedad hace de la proyección UTM una de las más convenientes para la resolución sobre el plano de los problemas topográficos.

42 UTM Uno de los fines de las transformaciones (conversiones) cartográficas es simplificar los cálculos, afectando a las figuras del elipsoide de las oportunas correcciones para poder sustituirlas por sus transformadas en el plano, y de esta forma poder resolver los problemas por trigonometría plana rectilínea.

43 Lalt. =22°45´00” Long. =102° 30´00” X= Y=

44 3° 6° ZONA-14 200,000 800,000 500,000 Disminuyendo hacia el Sur 102° WIncrementando hacia el Norte 0 m Disminuyendo hacia el Sur Ecuador

45 Pánfilo Natera Ojocaliente Villa Gonzalez Noria de Ángeles Luis Moya Loreto Villa Garcia

46

47

48 Geografícas a UTM

49 UTM Dos son los problemas prácticos que normalmente se le presentan al usuario de la cartografía, para trabajos topográficos, que son la determinación de direcciones y distancias. La dirección de la recta determinada por dos puntos A y B de la proyección (fig. 1) se define por el ángulo que esta recta forma con la parte positiva del eje de las Y. Este ángulo se llama orientación y no se debe confundir con el acimut geográfico, que como sabemos, es el ángulo formado por dirección AB con la meridiana geográfica. Esta orientación puede medirse directamente en el plano, con un transportador, si en el plano está dibujada la cuadricula; pero, si se desea mayor precisión debe calcularse analíticamente. Para ello de la figura 1, tenemos: tg OAB = XB - XA YB - YA

50 UTM La distancia que separa ambos puntos en la proyección se puede calcular por cualquiera de las fórmulas conocidas: DAB = √(XB - XA)2 + (YB -YA)2 DAB = XB - XA sen OAB DAB = YB - YA cos OAB

51 UTM Las fórmulas (for.6) nos dan el valor de la distancia entre los puntos A y B en la proyección, pero, en alguna ocasión se nos puede presentar la necesidad de conocer la distancia que realmente existe entre ellos en el elipsoide, que defiere de aquella debido al módulo de anamorfosis lineal K, conocido como factor de escala de la proyección UTM, y además, por la corrección de altitud, conocido como el factor de altitud para trabajos topográficos. El valor de K viene dado por la expresión: K = Δ l ´ = Proyección Δ l terreno

52 UTM Por lo tanto la distancia entre los puntos en el terreno será:Δ l = Δ l ´ K El valor de K puede obtenerse en las Tablas de la proyección o bien, puede calcularse con suficiente aproximación para cualquier punto, por la expresión: K = ( * q2 ) Siendo: q = (Xm – 500,000) * 10-6 En la que, Xm, es la abscisa del punto medio de la distancia expresado en metros.

53 UTM Ejemplo: Calcular la distancia existente entre los puntos A y B sobre el elipsoide, sabiendo que las coordenadas de dichos puntos son: A: X = , Y = B: X = , Y = La distancia en la proyección UTM es, según la primera fórmula de las XB - XA = YB - YA = D = √ = √ = m El valor de K lo obtenemos por la expresión para ello calculamos los valores de: Xm = = 2 q = ( – ) * 10-6 =

54 UTM K = * ( *q2) = * = Y por lo tanto la distancia en el elipsoide será De = = metros

55 UTM La corrección por altitud la podemos obtener del estudio de la figura, de la que se deduce: Dt = R + H De R Y de aquí: Dt = De * R + H R En la que H es la altitud media de los puntos A y B, y R el radio de la Tierra, ambos valores expresados en metros, tomándose para R el valor de metros, radio promedio

56 UTM Esta corrección es pequeña en general, pero no debe despreciarse ya que para altitudes de 700 metros, alcanza valores de cierta consideración, especialmente dignos de tenerse en cuenta cuando se realizan medidas con aparatos de medición electrónica de distancias, capaces de medir distancias de kilómetros con muy pocos centímetros de error. EJEMPLO: Obtener la distancia que separa en el terreno a los puntos A y B del ejemplo anterior; sabiendo que altitud media es de metros. De la fórmula se tiene: Dt = = mts

57 Proyección UTM El calculo inicia con las formulas que se indican a continuación. Para determinar la zona UTM correspondiente: Y el meridiano central La diferencia de longitudes se obtiene a partir de:

58 Proyección UTM Para la delta Este: Y

59 Proyección UTM Para la delta Norte: Y

60 Proyección UTM Para el calculo de la convergencia de cuadricula:

61 Proyección UTM y para el factor de escala:

62 Proyección UTM Finalmente el factor de escala y las coordenadas se obtienen con:

63 Proyección UTM Donde:

64 Proyección UTM La distancia meridional se obtiene a partir de:

65 Proyección UTM

66 Proyección UTM Para realizar la transformación de coordenadas se realizan los siguientes pasos: Se calcula el valor del meridiano central con la expresión (11). 2. Calculo de la diferencia de longitudes, mediante la expresión (13). 3. Calculo del radio de curvatura del primer vertical (N), con la expresión (15) 4. Calculo de las deltas Norte y Este con las formulas (16) y (14). 5. Obtener el valor de la distancia meridional mediante la expresión (23) 6. Calcular las coordenadas en la proyección UTM, mediante las expresiones (21)

67 Dichas coordenadas son las siguientes:CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO). Partimos en primer lugar de las coordenadas geográficas-geodésicas del vértice con el que haremos el ejemplo, que como he dicho antes es el vértice de Llatías. Los datos de este vértice están en principio en geodésicas sobre el elipsoide de Hayford (también llamado Internacional de 1909 o Internacional de 1924). Dichas coordenadas son las siguientes: También vamos a necesitar los datos básicos de la geometría del elipsoide de Hayford. Cuando digo datos básicos me refiero al semieje mayor (a) y al semieje menor (b). A partir de estos datos, aprenderemos a deducir otros parámetros de la geometría del elipsoide que nos harán falta en el proceso de conversión de coordenadas. Así, los datos referentes a los semiejes del elipsoide Hayford son:

68 CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO).Calculamos la excentricidad, la segunda excentricidad, el radio polar de curvatura y el aplanamiento: Aprovechamos para calcular también el cuadrado de la segunda excentricidad, pues nos hará falta en muchos pasos posteriores:

69 Seguimos con el radio polar de curvatura y el aplanamiento:CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO). Seguimos con el radio polar de curvatura y el aplanamiento: Lo primero que hacemos es convertir los grados sexagesimales (grados, minutos y segundos) a grados sexagesimales expresados en notación decimal (lo que se suele denominar normalmente "grados decimales"). Para ello operamos de la siguiente forma:

70 CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO).Una vez que tenemos la longitud y la latitud en grados decimales, procedemos a su paso a radianes, pues la mayor parte de los pasos posteriores se realizarán con entrada de datos en radianes. Operamos para ello de la forma: El siguiente paso es calcular el signo de la longitud. Para ello el proceso lógico es muy sencillo:

71 CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO).Una vez tenemos preparados los datos de longitud y latitud, podemos calcular el huso o zona UTM (UTM Zone) donde caen las coordenadas a convertir, con operaciones muy sencillas:

72 CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO).Con el huso ya conocido, el siguiente paso es obtener el meridiano central de dicho huso. El meridiano central es la línea de tangencia del cilindro transverso. Pero antes de seguir con los cálculos e introducir más conceptos, vamos a repasar algunos de los elementos principales de la proyección UTM. Así, conviene recordar que en la proyección UTM el cilindro transverso que se usa como superficie desarrollable, se va girando virtualmente para definir los diferentes husos (60) que rodean la tierra. Se empiezan a contar los husos por el antimeridiano de Greenwich y por eso la parte central de España cae en el huso 30, por estar en el lado opuesto del inicio de la numeración de husos, que queda al otro lado de la tierra. El meridiano central del huso es muy importante porque es el origen de las coordenadas X. Como el meridiano central dejaría la parte del huso situada a su izquierda con coordenadas X negativas, para evitar eso, se suma a todas las coordenadas X la cantidad de Esto hace que no existan valores negativos para las coordenadas X, puesto que se ha realizado un retranqueo del eje X de 500 km. Y algo semejante se hace para los valores de Y, cuyo origen es el ecuador. Como el ecuador está normalmente más lejos que el meridiano central del huso, las coordenadas Y suelen tener un guarismo más (en el caso de España, las Y son mayores que 4 millones). Si el ecuador es el origen de las Y, toda la parte situada al sur del mismo tendría coordenadas negativas. Para evitar eso, se suma el valor a los valores de Y, pero sólo en el caso de que se trate de coordenadas pertenecientes al hemisferio sur; si las coordenadas pertenecen al hemisferio norte, no se tocan los valores Y. Volviendo con el meridiano central del huso, éste también tiene la particularidad de que es automecoico. En teoría, para cualquier latitud que caiga dentro del rango de operación de la proyección UTM (intervalo entre los 84° N y los 80° S), el punto de

73 CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO).Expuestos estos conceptos, para saber mínimamente lo que estamos calculando, vamos a retomar los cálculos donde los habíamos dejado. Habíamos dicho que el siguiente paso es obtener el meridiano central del huso en el que caen las coordenadas geodésicas sobre las que operamos. La operación es muy sencilla: Ahora calculamos la distancia angular que existe entre la longitud del punto con el que operamos y el meridiano central del huso (véase la figura anterior). Es muy importante señalar que ambos datos tienen que ser introducidos en radianes. La longitud ya la habíamos traducido a radianes antes, pero no así el valor del meridiano central que acabamos de calcular. Para convertirlo a radianes multiplicamos por Pi y dividimos por 180:

74 CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO).A continuación debemos calcular una serie de parámetros que van encadenados unos a otros y que son el núcleo de las ecuaciones de Coticchia-Surace. Son muchas operaciones pero vereis que el proceso es muy rutinario y fácilmente programable:

75 Cálculo Final de Coordenadas:CONVERSIÓN DE GEOGRÁFICAS A UTM (PROBLEMA DIRECTO). Cálculo Final de Coordenadas: Una vez disponemos de todos los parámetros anteriores calculados, procedemos a la solución de las coordenadas UTM finales, de la forma: Para el caso de la solución de Y es muy importante recordar que si la latitud de las coordenadas geodésicas con las que operamos pertenece al hemisferio sur deberemos sumar el valor  al resultado obtenido. Como en el caso del ejemplo estamos operando con latitudes al norte del Ecuador, no realizamos tal operación:

76 Practica

77 Lista de puntos Geograficas Lista de puntosID Norte (Metro) Este (Metro) Altura (Metro) Código de característica 0001 N22°45' " O102°34' " INEG N21°51' " O102°17' " IZAC N22°46' " O102°36' " Lista de puntos ID Norte (Metro) Este (Metro) Elevación (Metro) Código de característica 0001 INEG IZAC

78