本资源基于arcpy的hurst指数算法实现，对需要做地理方面有所帮助

# -*- coding: utf-8 -*-
“““
Created on Thu Jun  1 13:06:45 2017

@author: CSY
“““

import arcpy
arcpy.CheckOutExtension（“spatial“）
from arcpy.sa import *
outpath=“I:/hents分析/S_2000-20017“
arcpy.CheckOutExtension（“spatial“）
arcpy.env.workspace=outpath
rasters=arcpy.ListRasters（“*““tif“）

for i in range（1len（rasters）+1）:
    
    #print “-------------------平均值-------------------------“
    Sum=0
    for ras in rasters[0:i]:     
        Sum=Sum+Raster（ras）
    mean=Sum/i
    
    #print “-------------------标准差-------------------------“
    s=0
    for fs in rasters[0:i]:
        s=s+（fs-mean）**2

    S= SquareRoot（s/i）
    outpath1=“/S/“
    sty1=“S“+str（i）+“.tif“
    out1=outpath+outpath1+sty1
    S.save（out1） 
    print sty1
    #print “-----------------