python批量对遥感影像进行归一化与数据清洗

掌握遥感数据清洗与归一化

目录

  1. 引言:为什么要进行数据清洗与归一化?
  2. [Step 1:遥感影像的读取](#Step 1:遥感影像的读取)
  3. [Step 2:处理无效数据(NoData 值)](#Step 2:处理无效数据(NoData 值))
  4. [Step 3:数据归一化处理](#Step 3:数据归一化处理)
  5. [Step 4:保存处理后的影像](#Step 4:保存处理后的影像)
  6. [Step 5:整合与执行](#Step 5:整合与执行)
  7. 完整代码
  8. 总结

引言:为什么要进行数据清洗与归一化?

遥感数据广泛应用于地球观测、环境监测、农业分析等领域,为我们提供了宝贵的地球表面信息。然而,处理遥感数据并不像直接打开照片那么简单。

遥感数据的复杂性使得我们无法直接对数据进行分析,原因包括:

  1. 噪声与空值 :遥感图像的某些像素由于传感器故障或云层遮挡可能无法获得有效数据,这些无效像素通常被标记为 NoData。如果不处理,可能会导致后续分析出现偏差。
  2. 数据尺度问题:遥感影像的数值范围取决于传感器及其使用的测量标准。例如,有些数据的反射率值可能在 0 到 1 之间,而有些则是 0 到 10000。这会导致在处理多波段图像或融合多源数据时出现不一致的问题。

因此,我们需要对遥感数据进行清洗和归一化处理,以保证数据的一致性和可靠性。接下来,我们将通过具体代码演示,如何使用 Python 对遥感数据进行有效的预处理。

环境配置

首先,确保你已经安装了 GDALnumpy 库,这两个库是本教程中遥感数据处理的核心工具。

安装命令:

bash 复制代码
pip install gdal numpy

接下来,开始我们的代码讲解。

Step 1:遥感影像的读取

遥感数据通常以 GeoTIFF 格式存储,它是一种包含地理信息的多波段图像文件格式。首先,我们需要使用 GDAL 库读取 .tif 文件并提取图像数据。

代码示例:

python 复制代码
from osgeo import gdal
import numpy as np

def get_img(dataset_img):
    width = dataset_img.RasterXSize  # 获取图像宽度(列数)
    height = dataset_img.RasterYSize  # 获取图像高度(行数)
    bands = dataset_img.RasterCount  # 获取波段数

    print(f"图像大小:{width} x {height},波段数:{bands}")
    
    proj = dataset_img.GetProjection()  # 获取投影信息
    geotrans = dataset_img.GetGeoTransform()  # 获取仿射变换矩阵
    
    img = dataset_img.ReadAsArray(0, 0, width, height)  # 读取图像数据
    return img, proj, geotrans

在这段代码中,我们使用 GDAL 来打开遥感影像并读取其行数、列数和波段数,同时提取图像的投影信息和仿射矩阵。接下来我们还需要处理每个波段的数据,确保无效像素被正确处理。

Step 2:处理无效数据(NoData 值)

在遥感影像的每个波段中,可能存在传感器无法获取的无效数据。这些数据通常被标记为 NoData,为了避免影响分析,我们需要将这些无效数据值设为 0 或其他适当的替代值。

代码

python 复制代码
def clean_nodata(image, dataset_img):
    for i in range(image.shape[0]):  # 遍历每个波段
        band = dataset_img.GetRasterBand(i + 1)
        nodata_value = band.GetNoDataValue()  # 获取该波段的 NoData 值
...

完整文章请查看该链接:

https://mp.weixin.qq.com/s/eTuDlztP1UkTFtSikqBSvg

相关推荐
起名字真南5 分钟前
【OJ题解】C++实现字符串大数相乘:无BigInteger库的字符串乘积解决方案
开发语言·c++·leetcode
tyler_download16 分钟前
golang 实现比特币内核:实现基于椭圆曲线的数字签名和验证
开发语言·数据库·golang
小小小~16 分钟前
qt5将程序打包并使用
开发语言·qt
hlsd#17 分钟前
go mod 依赖管理
开发语言·后端·golang
哇咔咔哇咔17 分钟前
【科普】conda、virtualenv, venv分别是什么?它们之间有什么区别?
python·conda·virtualenv
小春学渗透18 分钟前
Day107:代码审计-PHP模型开发篇&MVC层&RCE执行&文件对比法&1day分析&0day验证
开发语言·安全·web安全·php·mvc
杜杜的man21 分钟前
【go从零单排】迭代器(Iterators)
开发语言·算法·golang
亦世凡华、21 分钟前
【启程Golang之旅】从零开始构建可扩展的微服务架构
开发语言·经验分享·后端·golang
测试界的酸菜鱼36 分钟前
C# NUnit 框架:高效使用指南
开发语言·c#·log4j
GDAL36 分钟前
lua入门教程 :模块和包
开发语言·junit·lua