Antialiasing raster

Question

Estoy buscando una manera de suavizado entero DEM datos (hgt archivo) para flotar GeoTIFF. De esa manera quiero evitar gdaldem hillshade artefacto líneas que usted puede ver en esta imagen:

No quiero desenfoque DEM datos para mantener el resultado agudo. De hecho, estoy buscando una herramienta para interpolar píxeles en el siguiente (o mejor). Para simplificar la ilustración utiliza sólo 1D, 2D no:

Edit: Añadir imagen para la elevación cambiando su dirección. La entrada es de color azul (pen), la producción es discontinua gris (lápiz):

Apuesto a que alguien se ha inventado ya dicho algoritmo ;-). Yo soy "sólo" buscando una herramienta que implementa.

Respecto a la 2D, el algoritmo se aplica en ambas direcciones y el resultado sería el promedio de la misma.

Answer 1

He creado una herramienta con el algoritmo descrito.

El resultado (animando datos DEM originales / con alias):

Tenga en cuenta que los datos de DEM ruidosos (p. Ej., NASA SRTM 1 arc-second Global) no tienen este problema. Son "solo" ruidosos ;-).

Answer 2

Parece que usted está en necesidad de focal estadísticas como la media. Asumiendo que usted está utilizando GDAL Ejecución de una Focal Media explica que usted tendrá que usar gdal_calc.py para escribir su propia focal de la función. Para el objetivo, creo que el promedio debe hacer muy bien.