pycharm进行NDVI最大合成月尺度数据（超详细！！！）

一、NDVI数据来源

NDVI的数据是北大的GIMMS 标准化差异植被指数 1982-2022，我选择因为它比较长。网址：https://zenodo.org/records/8253971。

图是盗的（不会画）

二、在pycharm中导入arcpy

注意： 开始前，保证电脑已安装 Arc gis 。找到Arc Gis 自带的python程序的位置。（我电脑的位置  C:\Python27\ArcGIS10.2\python.exe）

1、在pycharm中的右下角找到Python X.X 版本，左击--->

2、Interpreter Settings---->

3、Add Interpreter---> 

4、Add Local Interpreter ... ---->

5、Existing ---> ··· ---->

6、选择路径 -----> OK

三、最大合成法合成月尺度NDVI数据代码

代码在下面，批量循环的代码写好了上传资源，在下面放链接。需要的兄弟们，收藏一下，等我消息。

# encoding: utf-8
# 上述代码解决 SyntaxError: Non-ASCII character '\xe6' in file问题

import arcpy
arcpy.CheckOutExtension("spatial")  # 检查模块许可
from arcpy.sa import *

# Set environment settings
workspace = r"D:/Users/shuow/Desktop/NDVI/NDVI-10KM/PKU_GIMMS_NDVI_AVHRR_MODIS_consolidated_1982_1990" # NDVI数据存放目录
arcpy.env.workspace = workspace
#####
su12 = ["PKU_GIMMS_NDVI_V1.2_19820101.tif","PKU_GIMMS_NDVI_V1.2_19820102.tif"]  #1982年1月NDVI数据1-15天和16-30天
outCellStatistics = CellStatistics(su12, "MAXIMUM", "DATA") # 合成2002年12个月的NDVI数据，方法为 MAXIMUM ， 无数据的地方以有数据的最大值为准DATA
outpath = 'D:/Users/shuow/Desktop/test/ndvi/' + str("PKU_GIMMS_NDVI_V1.2_19820101.tif") #输出路径
outCellStatistics.save(outpath) #保存合成后的NDVI数据

放一个成功结算画面和输出的结果图。

---

批量处理数据代码链接：https://download.csdn.net/download/PETTER_WEN/88870389?spm=1001.2014.3001.5501

---