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

相关推荐
武当豆豆34 分钟前
C++编程学习阶段性总结
开发语言·c++
学不动CV了1 小时前
C语言32个关键字
c语言·开发语言·arm开发·单片机·算法
你怎么知道我是队长2 小时前
python-enumrate函数
开发语言·chrome·python
小屁孩大帅-杨一凡2 小时前
如何解决ThreadLocal内存泄漏问题?
java·开发语言·jvm·算法
大熋2 小时前
Playwright Python 教程:网页自动化
开发语言·python·自动化
赟赟、嵌入式2 小时前
imx6ul Qt运行qml报错This plugin does not support createPlatformOpenGLContext!
开发语言·qt
A7bert7772 小时前
【YOLOv8-obb部署至RK3588】模型训练→转换RKNN→开发板部署
linux·c++·人工智能·python·yolo
cdg==吃蛋糕3 小时前
selenium 使用方法
开发语言·python
Y1nhl3 小时前
力扣_二叉树的BFS_python版本
python·算法·leetcode·职场和发展·宽度优先
爱掉发的小李3 小时前
前端开发中的输出问题
开发语言·前端·javascript