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

相关推荐
艾莉丝努力练剑10 分钟前
【LeetCode&数据结构】单链表的应用——反转链表问题、链表的中间节点问题详解
c语言·开发语言·数据结构·学习·算法·leetcode·链表
橡晟4 小时前
深度学习入门:让神经网络变得“深不可测“⚡(二)
人工智能·python·深度学习·机器学习·计算机视觉
墨尘游子4 小时前
神经网络的层与块
人工智能·python·深度学习·机器学习
倔强青铜34 小时前
苦练Python第18天:Python异常处理锦囊
开发语言·python
u_topian5 小时前
【个人笔记】Qt使用的一些易错问题
开发语言·笔记·qt
企鹅与蟒蛇5 小时前
Ubuntu-25.04 Wayland桌面环境安装Anaconda3之后无法启动anaconda-navigator问题解决
linux·运维·python·ubuntu·anaconda
autobaba5 小时前
编写bat文件自动打开chrome浏览器,并通过selenium抓取浏览器操作chrome
chrome·python·selenium·rpa
珊瑚里的鱼5 小时前
LeetCode 692题解 | 前K个高频单词
开发语言·c++·算法·leetcode·职场和发展·学习方法
AI+程序员在路上5 小时前
QTextCodec的功能及其在Qt5及Qt6中的演变
开发语言·c++·qt
xingshanchang5 小时前
Matlab的命令行窗口内容的记录-利用diary记录日志/保存命令窗口输出
开发语言·matlab