Python:arcpy批量重分类实现

Python:arcpy批量重分类实现,第1张

Python:arcpy批量重分类实现

今天没有写新的东西,把之前写的一个小代码上传一下吧,免得之后自己也为忘记。利用arcpy分别对20年的NDVI进行重分类,代码很简单。。。。

重分类何需这么复杂,大可不必写那么多,虾面我直接贴代码吧!

# -*- coding: UTF-8 -*-
import os
import arcpy
from arcpy.sa import *

'''
我需要对NDVI,**,**,**,**进行重分类,
本代码对20年的NDVI进行重分类,我的数据最大不超0.9,在我的文件夹下共有20张tif,hhh,2021.12.13写!
'''
def reclass(inputPath, outputPath):
    arcpy.env.workspace = inputPath
    ras_List = arcpy.ListRasters("*", "tif")
    reclassField = "VALUE"
    '''
    RemapRange([[分类下限,分类上限,新赋的值],[以此类推],...])
    '''
    remap_NDVI = RemapRange([[-1, 0, 1], [0, 0.1, 2], [0.1, 0.2, 3], [0.2, 0.3, 4], [0.3, 0.4, 5], [0.4, 0.5, 6],
                             [0.5, 0.6, 7], [0.6, 0.7, 8], [0.7, 0.8, 9], [0.8, 0.9, 10]])

    for ras in ras_List:
        ra = os.path.join(inputPath, ras)
        outRecalssify = Reclassify(ra, reclassField, remap_NDVI, "NODATA")
        out_ra = os.path.join(outputPath, ras[0:9] + "reclass.tif")
        outRecalssify.save(out_ra)
        #print(out_ra)

if __name__ == '__main__':

    inputPath = r"C:\Users\NDVI"
    outputPath = r"C:\Users\reclass"
    reclass(inputPath, outputPath)

就这样吧!!!!!!!!!!!!!!

欢迎分享,转载请注明来源:内存溢出

原文地址: https://outofmemory.cn/langs/789174.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-05-05
下一篇 2022-05-05

发表评论

登录后才能评论

评论列表(0条)

保存