Python基于Arcpy对多时间尺度NDVI进行最大值合成

Python基于Arcpy对多时间尺度NDVI进行最大值合成

1. 前期数据准备

2. 处理代码

前言：在现实研究中很容易拿到日尺度、月尺度的NDVI栅格影像，繁多的栅格数据需要按照研究的目的进行重新聚合，手动在Arcgis中处理不仅耗时而且容易出错，因此使用Arcpy对同规格数据聚合十分必要。

**前期数据准备：**文件夹中准备所有需要进行最大值合成的TIFF格式数据，合成过程中对文件夹中所有的栅格进行合成。

代码如下：
import arcpy
from arcpy import env
import glob
import os

path = r"C:\NDVI"

os.chdir(path)
env.workspace = path
rasters = glob.glob("*.tif")
name = str(rasters).replace("’, '", “;”)[2:-2]
out_name = “NDVI_Composite.tif”
arcpy.MosaicToNewRaster_management(input_rasters=name,
output_location=path,
raster_dataset_name_with_extension=out_name,
pixel_type=“32_BIT_FLOAT”,
number_of_bands=1, mosaic_method=“MAXIMUM”)

代码运行需注意：
1.修改path变量为自己NDVI所在的文件夹；
2.输出最大值合成的tiff影像在path路径中；
3.输出的最大值合成影像命名由out_name变量决定， 默认为"NDVI_Composite.tif"。

以上只是根据个人需要进行编码，如有不妥之处，请大家批评指正。

如果大家觉得有帮助，请收藏和点赞，读者的支持是我们创作的源源动力。