使用GRASS平滑DEM？

我有一个SRTM DEM，我想从中创建阴影。我在GRASS中创建了阴影浮雕，效果非常好，但是有点粗糙，因为该区域接近平坦并且DEM的分辨率为90m。

我想要使​​DEM更平滑以生成平滑的阴影浮雕。是否有算法或插值方法可以做到这一点？

这是获得想法的阴影浮雕，我想弄平这些小凸起：

r.denoise来自GRASS AddOns Wiki的 John Stevenson's怎么样：

使用Sun Xianfang 的降噪算法，r.denoise降噪（平滑/去斑）地形数据，特别是从雷达数据（包括SRTM）导出的DEM 。它旨在保留锐利的边缘并在对原始数据进行最小更改的情况下进行去噪。

我从这个网站上进一步读到（我也对上面的动画表示赞赏），一种更通用的方法是使用Esri ASCII Grid文件。mdenoise（可从Sun的网站下载）的位置必须位于PATH变量中（例如Windows用户：将OSGeo4w或FWTools安装MDenoise.exe到该bin文件夹中）。然后，例如，您可以使用以下shell命令来处理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

降噪已获得GNU许可，请参见此处

美国国家公园管理局首席制图师汤姆·帕特森（Tom Patterson）提供了一些出色的教程，介绍如何使用DEM数据制作漂亮的阴影浮雕。他的部分工作流程涉及使用Natural Scene Designer和Adobe Photoshop。

对于我自己的工作流程，在渲染阴影浮雕之前，我喜欢使用GDAL重新采样DEM的大小。这通常有助于减少细节和噪音，更不用说文件大小了。我在Github上有一个教程，演示了如何做到这一点。

基本上，该过程是为输出DEM指定小于原始文件大小的新宽度和/或高度（以像素为单位）。例如，执行以下操作：

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

将DEM的宽度从原来的3800像素减小到3000像素。将高度设置为0将使GDAL根据原始纵横比确定新文件的最佳高度。

我使用了GIMP插件小波降噪工具，并获得了良好的效果和快速的效果：

GIMP插件小波降噪

您可以在DEM上运行一个简单的过滤器-例如5x5移动窗口平均过滤器。

在ArcMap中，我喜欢使用Spatial Analyst>焦点统计工具在制作轮廓之前使激光雷达平滑。这与我相信的上述功能类似。