arcpy 实现DEM综合 对DEM进行不同分辨率重采样

import arcpy
from arcpy import env
from arcpy.sa import *  # 取出空间分析工具中的所有工具
from glob import glob  # 方便遍历shp文件
import sys
import os
arcpy.env.overwriteOutput =True
##-----------DEM综合----
print '------------------开始对DEM进行重采样--------------------'
# Local variables:
path = r'G:\0pythonprogress\public\demresample'
arcpy.env.workspace = path
inputDem = r"\90MDEM.tif"
outputDem = r"G:\0pythonprogress\public\demresample\CUBIC"
def DemResample (frequency,interval):
    mard = 90
    for i in range(0,frequency):
        mard =mard + interval
        print '当前的DEM分辨率为'+str(mard)+'米'
        output = outputDem+r'\{Mard}MDem.tif'.format(Mard=mard)
        arcpy.Resample_management(inputDem, output, "{Mrad} {Mrad}".format(Mrad = mard), "CUBIC")
        print("{Mrad}".format(Mrad = mard)+'米'+'已完成')
DemResample(100,30)