我有一个文件夹的解压缩文件从陆地卫星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"
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)
谢谢你的帮助
目前没有回答
相关问题 更多 >
编程相关推荐