python_MODIS HDF数据转为tif并拼接图像

首先是转换格式，原始hdf数据是这样

# -*- coding:utf-8 -*-
import os
import arcpy
from arcpy import env

sourceDir=u'E:' #中文路径 输入
targetDir=u'E:' #英文路径 输出

arcpy.CheckOutExtension("Spatial")
env.workspace = sourceDir
arcpy.env.scratchWorkspace = sourceDir
hdfList = arcpy.ListRasters('*','hdf')
for hdf in hdfList:
    print hdf
    eviName=os.path.basename(hdf).replace('hdf','tif')
    outname=targetDir+'\\'+eviName
    print outname
    data1=arcpy.ExtractSubDataset_management(hdf,outname, "Eight_Day_Snow_Cover") #Eight_Day_Snow_Cover是子数据集名称
##	
print 'all done'

接下来是获取每个HDF/TIF 年内天数数据，生成列表或文件，以方便后期判断合成日期是否一致

import os
import glob
import pandas as