¿Cuál es la forma correcta de reproyectar una trama de población?

Question

Estoy trabajando con un raster de población que está en el sistema de coordenadas WGS84, y quiero llevarlo a un sistema de coordenadas proyectadas (Albers Equal Area), utilizando ArcGIS 10.1.

Según los metadatos, las celdas de la cuadrícula tienen un tamaño de "30 segundos de arco (1 km)" (en un lado presumiblemente) y sus valores representan el recuento de la población en esa zona de 1 km^2.

Cuando utilizo el Proyecto Raster, la herramienta rellena automáticamente el parámetro "Tamaño de celda" con un valor en torno a los 830 metros. Cuando comparo ventanas equivalentes de 3 por 3 de celdas de datos entre las versiones WGS84 y Albers, me encuentro con que los valores han cambiado.

Aquí está WGS84 (ignora la cuarta columna de celdas por ahora):

Aquí está el equivalente de Albers de la ventana de 3 por 3 de arriba:

Obsérvese que, cuando el de WGS84 se proyecta a Albers, los valores de 3ª columna son "desplazados" por los valores del 4ª columna. Supongo que está haciendo algún tipo de remuestreo del vecino más cercano debido al tamaño alterado de las celdas. Los rangos en los 2 rasters también difieren en ~2,000 en el extremo superior.

Otros puntos de confusión:

1. ¿Cómo pueden los metadatos afirmar que 30 segundos de arco = 1 km cuando la distancia de un segundo de arco varía en función de la longitud? ¿Se refiere a la conversión en el ecuador? Tenía la impresión de que en el ecuador, 1 segundo de arco = 30 m, por lo que 30 segundos de arco = 900 m. ¿Los metadatos afirman simultáneamente y de forma implícita "en el ecuador" y redondean?
2. En la herramienta Project Raster de ArcGIS, ¿de dónde procede el tamaño de celda por defecto de 830,4934540? ¿Significa esto que cada celda debe tener 830 m de lado, lo que equivaldría a unos 0,69 km^2? ¿Debo modificarlo para que se ajuste más al valor de los metadatos de 1 km^2 por celda?

¿Cuál es la forma correcta de realizar esta reproyección sin dejar de respetar los recuentos de población originales? ¿Por qué se alteran claramente mis valores para que no haya una correspondencia uno a uno y cuál es la mejor práctica para tener en cuenta esto?

Answer 1

Teniendo en cuenta sus ejemplos anteriores, debe estar utilizando la opción "Vecino más cercano" para realizar el remuestreo en la herramienta Raster del proyecto. Por eso algunos valores "desaparecen". Qué más puede hacer, toma el valor más cercano (medido por la posición central de la celda) y lo asigna a la nueva celda en el nuevo sistema de coordenadas. El histograma de los valores de las celdas entre la antigua y la nueva nunca puede ser el mismo al proyectar. Intente utilizar las opciones bilineal o cúbica para obtener un "promedio" de las entradas para el valor de la nueva celda. Si todavía quiere que sean números enteros, entonces después puede hacer un "int" en la salida. En cuanto al tamaño de las celdas, se calcula por defecto como se ha mencionado anteriormente. Puede establecer el tamaño de la celda y el punto de origen en la herramienta o utilizando los ajustes del entorno.