8. 4波段遥感影像获取前3波段和计算添加ndvi

最新推荐文章于 2025-11-09 22:28:08 发布

原创

最新推荐文章于 2025-11-09 22:28:08 发布 · 1.5k 阅读

15 ·

CC 4.0 BY-SA版权

文章标签：

#python #深度学习 #矩阵

本文介绍了如何使用GDAL库配合numpy对遥感数据进行读取，包括4波段数据的前3波段提取，以及添加NDVI波段并以tif格式保存的操作步骤。通过实例演示了如何进行数据处理和坐标保存，适合初学者理解遥感数据处理流程。

gdal读取了遥感数据后就变成了numpy数组了，之后就是numpy的操作了，numpy操作了之后再加坐标保存就行了，实际上很简单读取和写入看下面的连接就够了
(4条消息) 5. gdal实现对遥感影像的读写，信息统计和波段选择_xiaotiig的博客-优快云博客
https://blog.youkuaiyun.com/xiaotiig/article/details/118721427

1 提取4波段的前3波段

# 将4波段的遥感影像提取出前3波段

import os
import gdal
import numpy as np

#  读取tif数据集
def readTif(fileName):
    dataset = gdal.Open(fileName)
    if dataset == None:
        print(fileName + "文件无法打开")
    return dataset

#  保存tif文件函数
def writeTiff(im_data, im_geotrans, im_proj, path):
    if 'int8' in im_data.dtype.name:
        datatype = gdal.GDT_Byte
    elif 'int16' in im_data.dtype.name:
        datatype = gdal.GDT_UInt16
    else:
        datatype = gdal.GDT_Float32
    if len(im_data.shape) == 3:
        im_bands, im_height, im_width = im_data