Python解决ENVI无法打开landsat8 level2级产品

Vorüarn0362

已于 2023-03-14 21:44:27 修改

阅读量1k

点赞数 1

分类专栏：遥感CV 文章标签： python 遥感 Landsat8 Powered by 金山文档

于 2023-03-14 21:40:50 首次发布

本文链接：https://blog.csdn.net/Sgappti/article/details/129540960

版权

遥感CV 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

原理，将MTL.txt后缀的文件的第一行改为GROUP = L1_METADATA_FILE，倒数第二行改为END_GROUP = L1_METADATA_FILE，即可。手动修改太累了，我们写一个Python程序自动化实现上述操作，省时省力。

路径结构如下图所示。

import os


ROOT='H:/2023_Files/Dataset/Landsat/NewAdded_0314_2023'

print('start correct!')
LC8_L1TP_FOLDERS=[(ROOT+'/'+i) for i in os.listdir(ROOT) if i.find('.')==-1]

MTL_FILES=[]
for LC8_FOLDER in LC8_L1TP_FOLDERS:

    curMTL_NAME=LC8_FOLDER+'/'+LC8_FOLDER.split('/')[-1]+'_MTL.txt'
    MTL_FILES.append(curMTL_NAME)

for fpath in MTL_FILES:
    print('fpath:',fpath)
    with open(fpath,'r') as file:
        conf_list=file.readlines()
    file.close()   
    # index=[x for x in range(len(conf_list))if 'LEVEL1'in conf_list[x]]
    fo = open(fpath, "w")
    
    conf_list[0] = 'GROUP = L1_METADATA_FILE\n'
    conf_list[-2]='END_GROUP = L1_METADATA_FILE\n'
    for i in range(len(conf_list)):
        fo.writelines(conf_list[i])         
    fo.close()

print('done!!!')