如何将来自不同陆地卫星平台的2个特定波段导入python?

2022-10-02 02:26:37 发布

您现在位置:Python中文网/ 问答频道 /正文

我有一个文件夹的解压缩文件从陆地卫星5号,7号和8号。我想导入红色和近红外波段运行NDVI。但是,这意味着我需要为陆地卫星8号引入波段4和5,为陆地卫星5号和7号引入波段3和4。我很难编写一个代码来导入这些波段。我对python完全陌生,所以这可能有点离谱,但我有:

import os
import arcpy
import re

from arcpy import env
from arcpy.sa import *

mydir = r"E:\Thesis\005005-006004\005005\US_Landsat_4-8_ARD"

存储索引的目录

rasters = r"E:\Thesis\Processing\005005"

图像列表;对应的xml文件

P_2Band = []

for root,dirs,files in os.walk(mydir):
    for name in files:
        if name.startswith("LC08")
            if name.endswith("4.tif") or if name.endswith("5.tif"):
                mypath = root+"\\"+name
                P_2Band.append(mypath)
                print(name)

        elif name.startswith("LE07") or if name.startswith("LT05")
            if name.endswith("3.tif" or "4.tif"):
                mypath = root+"\\"+name
                P_Metadata.append(mypath)
                print(name)   

谢谢你的帮助


Tags: ornamefromimportifos波段rootarcpytifmydirthesismypathstartswithendswith