ArcPy最大值法由NDVI月度数据转化为年度数据

最新推荐文章于 2023-08-01 18:57:03 发布

冬_冬_

最新推荐文章于 2023-08-01 18:57:03 发布

阅读量3.1k

点赞数 2

分类专栏：数据处理笔记文章标签： arcgis NDVI python

本文链接：https://blog.csdn.net/weixin_46629224/article/details/121802039

版权

数据处理笔记专栏收录该内容

21 篇文章 9 订阅

订阅专栏

下载中国5km分辨率逐月NDVI数据集（1982-2020）

配置python环境

编写代码

下载地址：国家地球系统科学数据中心数据详细信息http://www.geodata.cn/data/datadetails.html?dataguid=239118756960240&docid=474

配置python环境

配置方法可以参考我的另一篇文章

Pycharm中安装arcpy_冬_冬_的博客-CSDN博客_arcpy安装

编写代码

#!/usr/bin/env python
# -*- coding: UTF-8 -*-
'''
@Project ：arcpy_WSI 
@File    ：NDVI_sum.py
@Author  ：Ryo
@Date    ：2021/12/8 19:55 
'''
'''This program is for calculate annual NDVI dataset by monthly dataset '''
import arcpy
arcpy.env.overwriteOutput = True
from arcpy import env
from arcpy.sa import *
import os
import re
env.workspace = r'E:\BaiduNetdiskDownload\new WSI\NDVI.mdb'
luccpath=r'E:\BaiduNetdiskDownload\new WSI\Catchment_LUCC'
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput = 1
arcpy.env.overwriteOutput = 1

# Convert tiff dataset to rasters
for i in range(1,468):
    year=1982+i//12
    month=i%12+1
    print(year,month)
    if year<2000:#Extract data after 2000
        continue
    ##Create raster layer from single raster dataset with clipping feature
    arcpy.MakeRasterLayer_management(r"E:\FireFoxDownload\CDRMVCCN.tiff", "rdlayer", "#", "", str(i))
    arcpy.CopyRaster_management("rdlayer",str(year)+str(month))

for year in range(2000,2021):# turn monthly data to annualdata by max
    # Check out the ArcGIS Spatial Analyst extension license
    arcpy.CheckOutExtension("Spatial")

    # Execute CellStatistics
    outCellStatistics = CellStatistics([Raster(str(year)+str(month)) for month in range(1,13)], "MAXIMUM", "DATA")

    # Save the output
    outCellStatistics.save("NDVI"+str(year))
    print(year,'finished')