自己写的一个TIF波段自动融合的练习代码

@自己写的一个波段自动融合的代码 ，python的
作用是把landsat波段全融合并重命名。放在文件夹里直接双击就好。

import arcpy
import os
from arcpy import env
inpath=os.path.abspath('.')
env.workspace = inpath
fcs = arcpy.ListFiles()
bands=[]
for fc in fcs:
 f=arcpy.Describe(fc)
 if f.extension == "TIF": bands.append(f.name)
arcpy.CompositeBands_management(bands, bands[-1].rstrip('TIF.B0123456789')+".tif")
#print bands[-1].rstrip('TIF.B0123456789')
print "done"