Estoy usando QGIS 2.0.1. Quiero calcular las áreas (en km^2) de todos los municipios de México. Entiendo que necesito convertir mi mapa de WGS84 a UTM para obtener unidades en km, pero no estoy seguro de cómo explicar el hecho de que México abarca seis zonas UTM (de la 11N a la 16N), o qué CRS debo elegir de los cientos disponibles.
¿Qué CRS debería usar para esta tarea?
Yo evitaría utilizando la proyección UTM para comparar las áreas de seis zonas diferentes debido a las diferentes distorsiones de área. En cambio, el Proyección acimutal de igual área de Lambert (LAEA) es la mejor opción, porque preserva el área. Por ejemplo, puede utilizar el "Atlas Nacional de EE.UU. Área Equal" ( EPSG:2163 ) o puede definir uno personalizado para México simplemente cambiando las coordenadas del origen natural utilizado en el SIR.
¿Utiliza metros o grados decimales para las unidades? Obtengo exactamente las mismas áreas cuando uso EPSG:2163 o WGS 84.
Tengo datos sobre las áreas (en km2) de un subconjunto de municipios, que puedo utilizar para comprobar mis resultados. Por ejemplo, la superficie del municipio 1002 es de 547,22 km2. Cuando lo calculo en QGIS, obtengo 0,048057 (utilizando EPSG:2163 o WGS84). Así que creo que QGIS me está dando áreas en grados decimales, porque el factor de conversión es 106,71^2 (=547,22/0,048057). Pero este factor de conversión varía según los municipios. ¿Cómo puedo hacer que QGIS muestre la superficie en km2?
Mi Project Properties | General | Measurements son los ajustes: Elipsoide: Ninguno/ Planimétrico; Distancia: metros; Área: km cuadrados. Pero sigo obteniendo lo mismo que antes.
Si quieres la superficie de todos los municipios de todo el país no puedes utilizar UTM. Esto se debe a que México se encuentra en las zonas UTM 11 a 16. Si quieres calcular áreas en km2 debes usar la proyección cónica conforme de Lambert del INEGI. En QGis 2.14 puedes encontrar que las EPSG 6362 (ITRF92) y 6372 (ITRF2008)tienen los parámetros específicos para esta proyección que son:
+proj=lcc +lat_1=17,5 +lat_2=29,5 +lat_0=12 +lon_0=-102 +x_0=2500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0 +units=m +no_defs
Como puede ver, estos parámetros definen las unidades de medida de las coordenadas en metros. Si quieres tener los mismos parámetros en otro datum debes crear un CRS personalizado (settings/custom CRS). Por ejemplo:
Cónica Conforme de Lambert del INEGI, datum WGS84:
+proj=lcc +lat_1=17,5 +lat_2=29,5 +lat_0=12 +lon_0=-102 +x_0=2500000 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0 +units=m +no_defs
Cónica Conforme de Lambert del INEGI, datum NAD27:
+proj=lcc +lat_1=17,5 +lat_2=29,5 +lat_0=12 +lon_0=-102 +x_0=2500000 +y_0=0 +ellps=clrk66 +towgs84=-12,130,190,0,0,0,0 +units=m +no_defs
Espero que sea de ayuda
Gerardo
He actualizado a 2.14.13, pero no encuentro EPSG:6263 o EPSG:6272 . Intentaré hacer una proyección personalizada.
Tengo la 2.14.12 y tengo los dos SIR que menciono. Mi 2.14.12 es una instalación stand allone en win 64. También tengo la instalación 2.18.5 Osgeo4w. Puedo encontrar los CRS no sólo por el número, si escribo México en el filtro puedo ver todos los CRS que tienen un México en el nombre (incluso Nuevo México). Los códigos correctos son 6362 y 6372 no 6263 o 6272....
Lo siento, escribí mal los códigos. Pero tanto en 2.14.13 como en 2.18.5, 6362 y 6372 no aparecen. Estoy en Linux Mint 17.3.
Para áreas pequeñas, como ya has dicho, lo mejor sería elegir la proyección de la zona UTM correspondiente.
En cambio, si se trata de todo el país, creo que lo mejor sería elegir la zona UTM de la capital. Para este caso concreto, yo utilizaría la zona UTM 14 Norte ( EPSG:32614 ), que parece ser el que corresponde a la Ciudad de México.
¿Cuál es la opción de PolyGeo? Utilizando el elipsoide también parece que el área de salida en grados decimales.
I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.
0 votos
Así es como calculo las áreas: primero, establezco la capa CRS en WGS84; luego, uso
save aspara crear una nueva capa con el CRS deseado. NO utiliceset layer CRS. A continuación, utilice la calculadora de campos (tabla de atributos, ctrl-i) con$area.0 votos
Gerardo sugiere utilizar la cónica conformada de lambert. afalciano recomienda la azimutal de lambert de igual área. Asistente de proyección dice que hay que usar la cónica de igual área de Albers. ¿Cómo debo elegir entre ellas?