import fnmatch
import os
import arcpy
FCIout_path = '**********************'
FCIout_path02 = '*****************************'
for i in range(1, 86):
print(i)
_B12 = fnmatch.filter(os.listdir(FCIout_path), str(i) + '_B12_*.tif')
_B11 = fnmatch.filter(os.listdir(FCIout_path), str(i) + '_B11_*.tif')
_B8A = fnmatch.filter(os.listdir(FCIout_path), str(i) + '_B8A_*.tif')
B12 = FCIout_path + _B12[0]
B11 = FCIout_path + _B11[0]
B8A = FCIout_path + _B8A[0]
s = [B12, B11, B8A]
v1_FCI_tif = FCIout_path02 + str(i) + '_FCI.tif'
arcpy.CompositeBands_management(s, v1_FCI_tif)
Arcpy-波段合成
最新推荐文章于 2023-08-20 13:37:43 发布