使用ENVI读取Earth Engine下载的图像波段名的方法

最新推荐文章于 2023-03-13 14:08:41 发布

desertsTsung

最新推荐文章于 2023-03-13 14:08:41 发布

阅读量1.6k

点赞数 2

分类专栏： Earth Engine Python

本文链接：https://blog.csdn.net/desertstsung/article/details/101613405

版权

引言

Google旗下Earth Engine平台的用户越来越多，高效云计算赢得了众多学术科研人员的青睐。熟悉这个平台的读者应该都知道，如果使用代码var ndvi = image.normalizedDifference(['B5', 'B4']).rename('NDVI')得到一个图像，并将其下载到本地，使用常用的遥感软件ENVI是无法读取到rename修改后的波段名的。

简介

基于以上原因，我使用Python3+GDAL写了一段简单的代码，可以通过输入tif图像，输出其对应的hdr头文件，这样就可以使用ENVI读取出它的波段名了。

代码

包括两个函数，getInfo实现读取波段数和波段名，wHDR实现写入相应的hdr头文件。

# -*- coding:utf-8 -*-
import os.path as path
import sys
from osgeo import gdal
from gdalconst import *

#获取波段数和波段名的函数
#参数ifn：输入TIFF文件名
#返回包括波段数和波段名的list
def getInfo(ifn):
    info = []
    ds = gdal.Open(ifn, GA_ReadOnly)
    
    if ds is None:
        print(