生态遥感指数的mk分析的准备工作

在获得逐年逐月的裁剪后的rsei图层后，首先需要将需要分析的波段进行波段合成
（这里不得不吐槽我的正版arcgis10.8，无论做几个波段的波段合成，结果出来后在arcgis文件夹里的图层都是band1，而且是一个东西，值全是一样的，但使用信息查看工具看到的数值又是正确的，根本就无法再提取原始波段；但在另一台电脑上就是正常显示，这就是正版的妖力吗？？？？）

安装mk分析python包，网址为https://www.cnpython.com/pypi/pymannkendall，pip即可
接下来开始做分析的准备工作；我不知道原因，但被指导要这么做👇

首先制作用来裁剪栅格图像的渔网
使用 创建渔网 工具，创建一个栅格图像范围的渔网，大小为栅格像元的大小，生成储存渔网的gdb
然后按属性分割，即可获取裁剪后的多个面要素
使用

选择中的按位置选择，将与栅格图像相交的面提取出来，就制作好了需要用来裁剪栅格图像的矢量面块块
然后就把栅格图像裁成一块一块

import arcpy
from arcpy import env
from arcpy.sa import *
import os, random, arcpy
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput = 1

arcpy.env.snapRaster = r"D:\out\cmp_rsei.tif"
for i in range(1, 77):

    normalname = r"D:\cmp_rsei.tif"
    outResult2 = r'H:\split_rsei\No' + str(i).zfill(2) + '.tif'
    mask = 'H:/split_kuai.gdb/T' + str(i)

    print (mask)
    outExtractByMask2 = ExtractByMask(normalname, mask)
    outExtractByMask2.save(outResult2)
    print(str(i).zfill(3), ' is finished.')

感觉有点难表达清楚，就直接放arcpy代码吧哈哈哈

将裁剪好的栅格图像进行批量mk分析即可，最后将生成的斜率小块图进行镶嵌至新栅格，就完成啦

看客也许读的云里雾里，但写这个主要是为了自己巩固下😂