使用 QGIS修复缺失数据的栅格

举报
此星光明 发表于 2022/04/14 00:23:22 2022/04/14
【摘要】 处理栅格数据时,有时可能需要处理数据间隙。这些可能是传感器故障、处理错误或数据损坏的结果。以下是航拍图像中数据间隙(即无数据值)的示例。  图片来源:© 土地委员会(COLA);桑给巴尔革命政府 (RGoZ),从 OpenAerialMap 下载。(注意:数据间隙是使用python脚本模拟的,不是原始数据集的一部分) ...

处理栅格数据时,有时可能需要处理数据间隙。这些可能是传感器故障、处理错误或数据损坏的结果。以下是航拍图像中数据间隙(即无数据值)的示例。

 图片来源:© 土地委员会(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

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。