#?-*-?coding:?utf-8?-*-?????
import?arcpy
from?arcpy?import?env
from?arcpy.sa?import?*

#讀取數(shù)據(jù)

workspace1=“K:/Data3/growth“
arcpy.env.workspace=unicode（workspace1“utf8“）
YS=arcpy.ListRasters（“2*““tif“）#選擇數(shù)據(jù)開頭為2的數(shù)據(jù)
XS=arcpy.ListRasters（“3*““tif“）#選擇數(shù)據(jù)開頭為3的數(shù)據(jù)
ZS=arcpy.ListRasters（“1*““tif“）#選擇數(shù)據(jù)開頭為1的數(shù)據(jù)
#?Check?out?the?ArcGIS?Spatial?Analyst?extension?license
arcpy.CheckOutExtension（“Spatial“）
‘‘‘
#生成時間序列圖層
i=0?#初始化i
for?y?in?YS:
????#生成時間序列圖層
????i=i+1
????out=“C:/Users/Admin/Desktop/v3/LST/“
????if?i<10:
????????out1=unicode（out“utf8“）+“30“+str（i）+“.tif“
????else:
????????out1=unicode（out“utf8“）+“3“+str（i）+“.tif“
????outx=Con（y>-3iy）
????outx.save（out1）
#讀取時間序列圖層????
arcpy.env.workspace=unicode（workspace1“utf8“）
XS=arcpy.ListRasters（“R*““tif“）
‘‘‘
#計