使用 QGIS修复缺失数据的栅格
处理栅格数据时,有时可能需要处理数据间隙。这些可能是传感器故障、处理错误或数据损坏的结果。以下是航拍图像中数据间隙(即无数据值)的示例。
图片来源:© 土地委员会(COLA);桑给巴尔革命政府 (RGoZ),从 OpenAerialMap 下载。(注意:数据间隙是使用python脚本模拟的,不是原始数据集的一部分)
如果数据间隙很小,则可以通过插入来自相邻像素的值来有效解决。我将概述解决此问题的 2 种方法。第一个使用 QGIS,另一个使用纯 Python。
此处显示的方法使用该gdal_fillnodata
工具应用反距离加权插值和平滑。正如文档中所指出的,这适用于填充连续栅格数据(例如高程)中的缺失区域。它也适用于不同数据(例如航拍图像)中的非常小的差距。如果您希望插入点数据以创建栅格,则应改用该gdal_grid
工具。这在 QGIS 中的处理 → 工具箱 → 栅格分析 → 网格...
修复 QGIS 中的数据缺口
GDAL 带有一个工具 gdal_fillnodata,可以从 QGIS 的处理工具箱中使用。
如果源栅格设置了无数据值并且与缺失数据值相同,则可以跳过此步骤。否则,第一步是将栅格的无数据值设置为数据间隙的像素值。从 Processing → ToolBox,搜索并找到Translate(转换格式)工具
在我们的示例中,无数据像素值为 0。为将指定的无数据值分配给输出波段选项设置值 0,并为转换后的栅格输入文件名。
现在我们准备好从 Processing Toolbox运行Fill nodata工具
此工具一次适用于 1 个频段。选择波段 1(红色)。设置最大距离以搜索要插值的值到1,因为我们只有 1 个像素间隙。将输出另存为01_red.tif并单击运行。保存带有01_等前缀的文件很重要,因为下一步将按文件名的字母顺序合并波段。
对波段 2(绿色)和波段 2(蓝色)重复该过程,为它们选择合适的文件名。您应该有 3 个没有填充数据值的单独栅格。现在我们可以将它们合并到一个文件中。从处理工具箱中搜索并找到合并工具。
在合并工具中,选择所有 3 个单独的栅格。选中将每个输入文件放入单独的带框。输入输出的文件名,然后单击运行。
生成的合并栅格将具有 3 个波段,无数据间隙将填充来自相邻像素的内插值。
这是显示版本之前和之后的动画。
文章来源: blog.csdn.net,作者:此星光明2021年博客之星云计算Top3,版权归原作者所有,如需转载,请联系作者。
原文链接:blog.csdn.net/qq_31988139/article/details/121442638
- 点赞
- 收藏
- 关注作者
评论(0)