在使用Python GDAL将MODIS的HDF数据转换为tif之后,可使用Python中的ArcGIS Mosaic to New Raster工具来拼接多幅影像。或者也可继续使用GDAL 来拼接。数据文件名如:MOD11A2.A2000049.h23v04.006.2015058135050.tif
# -*- coding: utf-8 -*-
"""
Created on Thu Aug 8 22:05:41 2019
@author: zoe
"""
import arcpy
import glob
import os
#指定影像所在目录
inws =r'E:\MOD11A2_LandSurfaceTem\2000_0208\try' #输入路径
outws = r'E:\MOD11A2_LandSurfaceTem\2000_0208\out' #输出路径
mask = r"E:\map\A_TLM_2\unionCopy2.prj" #投影文件,虽然在转换HDF为tif时已经重新投影过
#判断日期函数,由于MODIS数据文件是以年内天数命名,为了将其转为日期,并给镶嵌后的文件重新命名,如2001001变为2001_01_01
def out_date(year,day):
fir_day = datetime.datetime(year,1,1)
zone = datetime.timedelta(days=day-1)
return datetime.datetime.strftime(fir_day + zone, "%Y_%m_%d")
#利用glob包,将inws下的所有tif文件读存放到rasters中
rasters = glob.glob(os.path.join(inws, "*.tif"