Loading [MathJax]/jax/output/CommonHTML/config.js
前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >专栏 >使用 QGIS修复缺失数据的栅格

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

作者头像
此星光明
发布于 2024-02-02 03:51:58
发布于 2024-02-02 03:51:58
6360
举报

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

图片来源:© 土地委员会(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 个波段,无数据间隙将填充来自相邻像素的内插值。

这是显示版本之前和之后的动画。

本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2024-02-01,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
暂无评论
推荐阅读
编辑精选文章
换一批
PostGIS导入导出栅格数据
上一篇博文PostGIS导入导出ESRI Shapefile数据介绍了如何导入空间矢量数据到PostgreSQL中,紧接上一篇,本文将介绍如何使用PostGIS导入导出空间栅格数据。 实验环境和上一篇的博文一样: Ubuntu 16.04 LTS 。 数据为全球影像(ESRI ArcGIS提供的示例数据),下载连接:全球影像百度网盘下载 可以使用GDAL的gdalinfo命令查看其详细信息:
卡尔曼和玻尔兹曼谁曼
2019/01/22
4.9K0
PostGIS导入导出栅格数据
QGIS软件项目的开发框架
QGIS(Quantum GIS)是一个领先的开源桌面地理信息系统(GIS)软件,其开发框架主要基于以下核心技术和设计理念。
数字孪生开发者
2025/05/20
1480
QGIS软件项目的开发框架
手把手教你使用QGIS制作地图
QGIS是一款开源免费的地理信息系统软件,虽然比不上商业的ArcGIS软件,但是QGIS免费而且跨平台,值得学习!
卡尔曼和玻尔兹曼谁曼
2020/04/09
8.7K0
手把手教你使用QGIS制作地图
Python | GDAL处理影像
注意读取数据的数组下标不要越界!GDAL并不会自动帮你处理下标越界的问题,它只会报错。因此特别当你想用部分读取的方式处理一个很大的文件时,对边界的处理需要你特别的注意,必须正好读完不能越界也不能少读。
GIS与遥感开发平台
2022/04/29
4.7K0
Python | GDAL处理影像
Google earth engine——如何导入栅格数据?
开始上传后,“资源摄取”任务会出现在代码编辑器右侧的“任务”选项卡上。将鼠标悬停在任务管理器中的任务上会显示 ? 可用于检查上传状态的图标。要取消上传,请单击任务旁边的旋转图标。摄取完成后,资产将出现在您的用户文件夹中,并带有image 图标。
此星光明
2024/02/02
3530
Google earth engine——如何导入栅格数据?
QGIS小妙招-创建TAC边界
之前分享过 使用 mapinfo 创建 TAC 边界的方法,现如今,mapinfo 已经被禁用了,开源软件 QGIS 在通信网络优化工作中开始大量应用,今天分享一下如何使用 QGIS 创建 TAC 边界。
披头
2022/08/22
1.9K0
QGIS小妙招-创建TAC边界
Python提取大量栅格文件各波段的时间序列与数值变化
  本文介绍基于Python语言,读取文件夹下大量栅格遥感影像文件,并基于给定的一个像元,提取该像元对应的全部遥感影像文件中,指定多个波段的数值;修改其中不在给定范围内的异常值,并计算像元数值在每一景遥感影像中变化的差值;最终将这些数据保存为一个新的Excel表格文件的方法。
疯狂学习GIS
2024/07/22
2820
Python提取大量栅格文件各波段的时间序列与数值变化
GIS数据读取与转换
前面给大家分享了GIS格式总结,今天讲一下GIS数据的读取和转换,主要基于ArcGIS和QGIS这两个用的最广泛的GIS软件平台来讲。
郭好奇同学
2021/05/28
2.9K0
GIS数据读取与转换
遥感数据机器学习的准备工作:python将栅格数据提取至EXCEL
大部分我们处理的降水、气温等栅格数据的格式是nc形式,需要我们将他转换成栅格数据并导入至Arcgis中,进行下一步操作。
Twcat_tree
2023/02/28
1.6K0
遥感数据机器学习的准备工作:python将栅格数据提取至EXCEL
使用Rasterio读取栅格数据
有没有觉得用GDAL的Python绑定书写的代码很不Pythonic,强迫症的你可能有些忍受不了。不过,没关系,MapBox旗下的开源库Rasterio帮我们解决了这个痛点。
卡尔曼和玻尔兹曼谁曼
2019/01/22
2.2K0
PostGIS批量导入栅格数据
如果我们有一批以文件存储的影像数据如何利用PostGIS批量的导出到PostgreSQL数据库中进行管理呢? (单个数据的导入参见我上篇博文:PostGIS导入导出栅格数据)
卡尔曼和玻尔兹曼谁曼
2019/01/22
2.4K0
PostGIS批量导入栅格数据
GDAL简介
Geospatial Data Abstraction Library (GDAL)是使用C/C++语言编写的用于读写空间数据的一套跨平台开源库。现有的大部分GIS或者遥感平台,不论是商业软件ArcGIS,ENVI还是开源软件GRASS,QGIS,都使用了GDAL作为底层构建库。
卡尔曼和玻尔兹曼谁曼
2019/01/22
2.5K0
Python中gdal栅格影像读取计算与写入及质量评估QA波段筛选掩膜
  本文介绍基于Python语言gdal等模块对遥感影像加以处理的详细代码与操作。
疯狂学习GIS
2023/07/09
4000
Python中gdal栅格影像读取计算与写入及质量评估QA波段筛选掩膜
geotrellis使用(四)geotrellis数据处理部分细节
       前面写了几篇博客介绍了Geotrellis的简单使用,具体链接在文后,今天我主要介绍一下Geotrellis在数据处理的过程中需要注意的细节,或者一些简单的经验技巧以供参考。 一、直接操作本地Geotiff   如果不想将tiff数据切割成瓦片存放到集群中,也可以直接使用Geotrellis操作本地geotiff文件,可以直接使用SinglebandGeoTiff读取单波段的tiff,使用MultibandGeoTiff读取多波段tiff。 val geotiff = SinglebandG
魏守峰
2018/04/28
1.2K0
Python中gdal实现多幅栅格影像批量绘制直方图
  现需要对多幅栅格数据文件进行直方图绘制,具体绘制内容即各栅格图像像素数值的分布情况;所有栅格数据都保存在同一目标路径下,且均为.tif格式;而目标路径下具有其它非.tif格式的文件,以及不需要进行直方图绘制的.tif格式文件,因此需要在绘制前对目标路径下的文件列表加以筛选,只保留需要绘制直方图的栅格文件。
疯狂学习GIS
2023/07/09
3430
Python中gdal实现多幅栅格影像批量绘制直方图
大栅格数据如何更快运算
这两周我在使用python进行大量的栅格数据的运算,在运算过程中遇到了数据量超级大但算力不足的问题。通过这两周的探索,也慢慢找到了一些加快栅格数据计算的方法,和读者分享。
自学气象人
2023/06/21
4510
大栅格数据如何更快运算
python与地理空间分析(一)
在气象数据分析中,地理空间要素是一个必须考虑的关键特征项,也是重要的影响因素。例如气温会随着海拔的升高而降低,地形的坡向朝向也会影响风速的分布,此外,典型的地形会形成特定的气候条件,也是数据挖掘中可以利用的区域划分标准。数据分析中,地理空间分析往往能提供有效的信息,辅助进行决策。随着航空遥感行业的发展,积累的卫星数据也成为了数据挖掘的重要数据来源。 地理空间分析有好多软件可以支持,包括Arcgis,QGIS等软件平台,本系列文章将会着重分享python在地理空间分析的应用。主要包括地理空间数据的介绍,常用的python包,对矢量数据的处理,对栅格数据的处理,以及常用的算法和示例。 地理空间数据包括几十种文件格式和数据库结构,而且还在不断更新和迭代,无法一一列举。本文将讨论一些常用的地理空间数据,对地理空间分析的对象做一个大概的了解。 地理空间数据最重要的组成部分:
MeteoAI
2019/07/24
8.4K0
基于C++ GDAL依据像元个数批量裁剪栅格图像
  本文介绍基于C++语言的GDAL模块,按照给定的像元行数与列数,批量裁剪大量多波段栅格遥感影像文件,并将所得到的裁剪后新的多波段遥感影像文件保存在指定路径中的方法。
疯狂学习GIS
2024/08/01
2050
基于C++ GDAL依据像元个数批量裁剪栅格图像
QGIS软件系统性能优化
QGIS作为一个强大的桌面GIS软件,在处理大量数据或进行复杂分析时,性能优化至关重要。以下是一些针对QGIS软件系统性能优化的策略和技巧。
数字孪生开发者
2025/05/20
1530
QGIS软件系统性能优化
geotrellis使用(二十七)栅格数据色彩渲染
目录 前言 复杂原因及思路分析 实现过程 总结 一、前言        今天我们来研究一下这个看似简单的问题,在地理信息系统中颜色渲染应当是最基本的操作和功能,比如我们将一幅Landsat数据拖拽到Arcgis或者QGis等软件中,软件会自动为我们呈现出漂亮的图案,一切看似来的那么容易,但是在分布式海量空间数据的情况下实现色彩渲染操作实在也是要了命的。今天我们就接着上一篇文章中的数据处理(权且将色彩渲染归结到数据处理中)来介绍一下如何在Geotrellis中为栅格数据渲染漂亮的色彩。 二、复杂原因及思路分析
魏守峰
2018/04/28
1.3K0
相关推荐
PostGIS导入导出栅格数据
更多 >
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档