Cómo calcular el área en metros cuadrados - solución de problemas de proyección en un conjunto de datos global en QGIS

Question

Estoy utilizando la versión 2.12.3 de QGIS.

Tengo 3 shapefiles de polígonos: A, B y la intersección de A y B, a la que podemos llamar C.

Lo único que quiero hacer es añadir un campo a la tabla de atributos de cada shapefile que indique el área (en metros cuadrados o kilómetros cuadrados) de cada característica en esa capa. A, B y C tienen cada uno más de 10,000 características individuales (polígonos de áreas protegidas), distribuidas en todos los continentes.

He utilizado una proyección global (WGS 84/NSIDC EASE Grid-Global EPSG:3975, también he probado con 3410), que figura en 'Sistemas de Coordenadas Proyectadas', por lo que asumo que está proyectada y no es geográfica, y la he establecido como el sistema de coordenadas para el proyecto y cada una de las capas.

Cuando utilizo la herramienta de cálculo de campos y la función $area, los valores resultantes de campo definitivamente no están en metros cuadrados o kilómetros cuadrados y no se parecen en nada al área tal como la entiendo.

Lo he intentado varias veces después de trastear con configuraciones que estoy bastante seguro de que no entiendo, y he terminado con valores realmente pequeños (por ejemplo, 0.0000002345) o conjuntos aleatorios de valores que incluyen números negativos.

¿Alguien puede identificar dónde me equivoqué y sugerir cómo puedo completar esta tarea bastante simple?

Answer 1

WGS 84/NSIDC EASE Grid-Global parece ser una proyección estereográfica, la cual no conserva áreas. Esta es la razón por la que no puedes obtener resultados significativos.

La solución más sencilla es reproyectar tu capa a una proyección UTM adecuada para tu área de interés.

Answer 2

Todos los archivos shp utilizados tienen el EPSG 31468.

Ejecuto el modelo, agrego un campo y calculo el área, la suma de todos los registros será de 10275 m²

Cuando hago lo mismo con la calculadora de campos desde la interfaz gráfica, la suma de todos los registros será de 10311 m². Cuando uso otro GIS, siempre obtengo este resultado