使用大型光栅时出现奇怪的条纹

Wal*_*tiz 6 gis r raster geospatial r-raster

我有一个来自 ASTER 数据库的大高程栅格,使用raster::mosaic(). 每个瓦片代表地球的 1 度 x 1 度部分,估计栅格分辨率为 30m^2。我在下面包含了我用来镶嵌它们的代码,但我不认为这是问题所在。

Mosaic <- do.call(mosaic,c(list of rasters, tolerance = 1, fun=mean)
Run Code Online (Sandbox Code Playgroud)

运行后,我会得到以下海拔 tif,您可以在此处下载

高程栅格图像

然后我使用 测量地形坚固性raster::terrain(),使用地形坚固性指数作为我的方法:

TRI <- terrain(Mosaic, opt="TRI")
Run Code Online (Sandbox Code Playgroud)

这产生:

三光栅

如您所见,地图中有我无法解释的微弱条纹。由于我的研究范围很广,我需要将此栅格聚合为较粗的分辨率。因此,我也在 raster 包中使用了聚合函数。

TRI_Agg <- aggregate(TRI,fact = 255, fun = mean)
Run Code Online (Sandbox Code Playgroud)

这将产生以下栅格。 TRI聚合

在这里,您可以看到条纹在整个地图中更加明显且相当一致。我尝试在其他研究领域这样做,但问题仍然存在(即世界各地的光栅图块都会发生这种情况)。我无法解释为什么会出现这些条纹。我尝试了不同的方法来解决这个问题,包括进行邻域分析(如果有的话,填充缺失值)和在将它们拼接在一起之前处理瓷砖。这些都没有奏效。

当我之前聚合栅格时,粗糙度计算,条纹不会出现,这让我相信使用aggregate(). 但是,对于我的特定研究问题,我无法使用聚合栅格计算 TRI。