資源簡(jiǎn)介
本資源中對(duì)數(shù)據(jù)做了求平均值，標(biāo)準(zhǔn)差，累積離差，極差處理最終得到的結(jié)果是極差與標(biāo)準(zhǔn)差的比值，要想得到Hurst指數(shù)，還需要進(jìn)行最小二乘法代碼。
代碼片段和文件信息
#?-*-?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?“-------------------標(biāo)準(zhǔn)差-------------------------“
????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?“-----------------評(píng)論
共有 條評(píng)論