逐像元趋势分析+F值检验

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

@author: CSY
“““

import arcpy
arcpy.CheckOutExtension（“spatial“）
from arcpy.sa import *
import math
arcpy.CheckOutExtension（“spatial“）
workspace1=“K:/Data3/growth“
arcpy.env.workspace=unicode（workspace1“utf8“）
ys=arcpy.ListRasters（“2*““tif“）

#生成时间序列图层
i=0 #初始化i
for y in ys:
    #生成时间序列图层
    i=i+1
    out=“K:/Data3/growth/“
    if i<10:
        out1=unicode（out“utf8“）+“30“+str（i）+“.tif“
    else:
        out1=unicode（out“utf8“）+“3“+str（i）+“.tif“
    outx=Con（y>-10iy）
    outx.save（out1