Arcgis 栅格数据批量镶嵌 python代码

# -*- coding: utf-8 -*- # Name: SpatialJoin_Example2.py # Description: Join attributes of cities to states based on spatial relationships. # Requirements: os module # Import system modules # coding:utf-8 import sys # import sys # reload(sys) import time import os import arcpy import glob import itertools from calendar import monthrange import datetime from datetime import timedelta from dateutil.relativedelta import relativedelta from os.path import basename print(sys.getdefaultencoding()) def encode(s): return s.decode("utf-8", "ignore") # Set local variables # gridWorkspace = encode(r"D:\宣老师实验\process.gdb") # print(gridWorkspace) # grid3KmFeatureClass = os.path.join(gridWorkspace,"grid_3KM") # clipping_geometry_featureclass=os.path.join(gridWorkspace,"grid_10KM") # clipping_cjs_geometry_featureclass = os.path.join(gridWorkspace,"csj_xzqh") # 处理数据的根目录 processWksFolder = encode(r"F:\大三下\小学期\印度\AOD") # processWksFolder = encode(r"D:\GWRexperiment\DataProcessing_zhy\TempData") # 原始数据 # src_hdf_folder = os.path.join(processWksFolder,"0src_hdf"); # encode(r"D:\Study\Learning\GWR\Data\AOD\MOD04_2L\2015\1_geoloc_wgs84_tif\*") # tifs = geoloc_wgs84_tif_folder+"\*" # daysList = glob.glob(tifs) # print(daysList) # 输出结果 # daily_mosaic_proj_folder =os.path.join(processWksFolder,"temp_daily") #投影数据 mosaic_proj_monthly_tif_folder = os.path.join(processWksFolder, "temp_monthly") # m每月数据 for i in range(6, 9): # tif数据 geoloc_wgs84_tif_folder = os.path.join(processWksFolder, str(i) + "_10k") daysfiles = [] for dayfile in os.listdir(geoloc_wgs84_tif_folder): daysfiles.append(os.path.join(geoloc_wgs84_tif_folder, dayfile)) print(daysfiles) monthly_mosaic_file = os.path.join(mosaic_proj_monthly_tif_folder, str(i).zfill(3) + "_10k.tif") # if(os.path.exists(monthly_mosaic_file)): # continue print("镶嵌多天数据至一个月") arcpy.CreateRasterDataset_management(mosaic_proj_monthly_tif_folder, str(i).zfill(3) + "_10k.tif", "", "16_BIT_SIGNED", "", 1, "", "", "", "", "") arcpy.Mosaic_management(daysfiles, monthly_mosaic_file, "MEAN", "MATCH", -9999, -9999, "NONE", 0, "NONE") # arcpy.WorkspaceToRasterDataset_management(geoloc_wgs84_tif_folder, monthly_mosaic_file,"INCLUDE_SUBDIRECTORIES", "MEAN","MATCH", -28672, -28672)