Python 批量拼接Landsat遥感影像

pycharm 调用arcpy进行遥感影像批量拼接

#!/usr/local/bin/python
# -*- coding:utf-8 -*-
import arcpy,os
def mosaic(path='.',name=''):
    #指定工作目录，即存放影像的目录
    outFolder = r'H:\wangyu\pcos\2000useful\07'
    i=0
    for item in os.listdir(path):  #读出path中所有的目录为item
        item_path=os.path.join(path,item)  #加深文件夹
        if os.path.isdir(item_path):  #判断加深之后的目录是否是文件夹
            mosaic(item_path,name)#如果是在继续寻找更深层，目前为止都是对读出的第一个目录进行操作
        elif os.path.isfile(item_path):  # 判断加深之后的目录中是否包含文件
            if name in item:  # 判断name（特征）是否在item中
                a,b=os.path.split(item_path)
                bb = item[15:21]
                sanjiangping=os.path.join(outFolder,bb+'sanjiangping.tif')
                if os.path.exists(sanjiangping):
                    i=i+1
                    print(bb+'sanjiangping.tif'+'yicunzai')