¿Suavizar DEM usando GRASS?

¿Qué tal el r.denoise de John Stevenson, de la wiki de complementos de GRASS :

  

r.denoise denoises (alisar / despojar) datos topográficos , especial   DEM derivados de los datos de radar (incluido SRTM), que utilizan Xianfang Sun's   algoritmo de eliminación de ruido . Está diseñado para preservar los bordes afilados y para   denoise con cambios mínimos en los datos originales.

Leímásde este sitio web (también doy crédito por la animación anterior) ) que un método más genérico sería usar un archivo de Esri ASCII Grid. La ubicación de mdenoise (descargada de sitio web de Sun ) debe estar en su PATH variable (por ejemplo, usuarios de Windows: suelte MDenoise.exe en la carpeta bin con su instalación de OSGeo4w o FWTools). Luego, por ejemplo, puede usar el siguiente comando de shell para procesar el archivo de cuadrícula ASCII:

# gdal_translate -of AAIGrid my_dem.tif my_dem.asc      # convert to .asc
mdenoise -i my_dem.asc -n 5 -t 0.99 -o my_dem_DN.asc    # denoise
# gdal_translate -of GTiff my_dem_DN.asc my_dem_DN.tif  # convert back to .tif

Denoise está bajo la licencia GNU, consulte aquí

Tom Patterson, el principal cartógrafo del Servicio de Parques Nacionales de EE. UU. tiene excelentes tutoriales sobre cómo trabajar con datos de DEM para crear hermosos relieves sombreados . Parte de su flujo de trabajo consiste en utilizar Natural Scene Designer y Adobe Photoshop.

Para mi propio flujo de trabajo, me gusta usar GDAL para volver a muestrear el tamaño del DEM antes de representar un relieve sombreado. Esto a menudo ayuda a reducir la cantidad de detalles y el ruido, sin mencionar el tamaño del archivo. Tengo un tutorial sobre Github que demuestra cómo hacerlo.

Básicamente, el proceso consiste en especificar un nuevo ancho y / o alto (en píxeles) para un DEM de salida que es más pequeño que el tamaño del archivo original. Por ejemplo haciendo:

gdalwarp -ts 3000 0 -r bilinear kings_canyon_2228.tif kings_canyon_2228_rs.tif

Reducirá el ancho del DEM a 3000 píxeles del ancho original de 3800 píxeles. Establecer la altura en 0 permitirá a GDAL determinar la mejor altura del nuevo archivo según la relación de aspecto del original.

Puede ejecutar un filtro simple en el DEM, por ejemplo, un filtro de promedio de ventana móvil de 5x5.

Usé la herramienta de reducción de ruido wavelet del complemento GIMP y obtengo buenos resultados y rápido:

Reducción de ruido de Wavelet del complemento GIMP

En ArcMap, me gusta usar la herramienta de estadísticas focales Spatial Analyst > para suavizar el lidar antes de hacer contornos ... esto también podría ayudar en su caso de sombreado. Es una función similar a la de la publicación anterior, creo.