Fusión de datos ráster con diferentes resoluciones

Question

Necesito fusionar diferentes conjuntos de datos ráster que tienen una capa y contienen valores discretos a una pila de capas. Las resoluciones de los datos pueden tener cualquier valor, por lo que no sólo son múltiplos de las resoluciones de los otros rásters de entrada, por lo que hay problemas de alineación. Además, las coordenadas de las esquinas no son necesariamente idénticas.

¿Alguien tiene una idea creativa de cuál es la mejor manera de hacer la fusión sin perder información ya que una pila de capas sólo puede tener capas con la misma resolución? He utilizado gdal_merge.py y, por supuesto, hay bastante pérdida de información a través del remuestreo del vecino más cercano (aunque debo utilizar el vecino más cercano ya que quiero mantener los valores discretos).

Una cosa que se me ocurrió fue poligonizar cada trama y luego unirlas y rasterizar el polígono resultante con una resolución adecuada. Esto mantendría la mayor parte de la información, pero me temo que no funcionará correctamente con conjuntos de datos grandes. ¿Quizás alguien tenga experiencia con esto?

Answer 1

No es necesario convertir la trama en polígonos.

Querrás hacerlo en varias pasadas. En primer lugar, usted quiere volver a muestrear el conjunto de datos de menor resolución a tan cerca como usted puede conseguir y un múltiplo del conjunto de datos de mayor resolución utilizando la interpolación bilineal que debe obtener un mejor aspecto de remuestreo. Después, haz la fusión de forma normal usando tu vecino más cercano (que creo que gdal_merge.py sólo puede hacer de todas formas).

Para el primer paso necesitarías usar gdalwarp con el -r bilinear o -r cubic opción. Esto le dará un remuestreo decente para continuar con su fusión. gdalwarp también puede hacer el vecino más cercano si usted quería usar eso en su segunda pasada.