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

Question

我有一个来自 Landsat 5、Landsat 7 和 Landsat 8 的解压缩文件的文件夹。我想导入红色和 NIR 波段来运行 NDVI。但是，这意味着我需要为 Landsat 8 引入波段 4 和 5，为 Landsat 的 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)   

谢谢你的帮助！