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)
就这样吧!!!!!!!!!!!!!!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)