栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 软件开发 > 后端开发 > Python

Python:arcpy批量重分类实现

Python 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

Python:arcpy批量重分类实现

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:UsersNDVI"
    outputPath = r"C:Users\reclass"
    reclass(inputPath, outputPath)

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

转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/844425.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号