资源简介
本资源适用于对arcpy有所了解，针对需要对遥感栅格数据做相关分析和偏相关的高校学生
代码片段和文件信息
# -*- coding: utf-8 -*-
import arcpy
from arcpy import env
from arcpy.sa import *
#读取数据
workspace1=“K:/Data3/growth“
arcpy.env.workspace=unicode（workspace1“utf8“）
YS=arcpy.ListRasters（“2*““tif“）#选择数据开头为2的数据
XS=arcpy.ListRasters（“3*““tif“）#选择数据开头为3的数据
ZS=arcpy.ListRasters（“1*““tif“）#选择数据开头为1的数据
# 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“）
‘‘‘
#计
评论
共有 条评论