代码分享 Python |批量删除栅格影像异常值

任务举例：下载的NPP影响最大值为63355，通过查看栅格数据直方图要删除异常值。

数据准备：文件夹(工作空间)-含有需要批量处理的tif影像

工具准备：Python arcpy 环境
操作：替换代码中的路径即可

# -*- coding: UTF-8 -*-
import arcpy
import os
import glob
import arcpy
from arcpy.sa import *
​
arcpy.CheckOutExtension("ImageAnalyst")  # 检查许可
arcpy.CheckOutExtension("spatial")
​
arcpy.env.workspace = "E:/1/"  #工作空间
rasters = arcpy.ListRasters("*", "tif") #遍历工作空间中的tif格式数据
​
outPath = 'E:/2'
​
whereClause = "VALUE > 3550 "  # 去除异常值
​
# 循环rasters中的所有影像，进行去异常值操作
for ras in rasters:
    outrasters = outPath + str(ras)#更改输出栅格名字
    outSetNull = SetNull(ras, ras, whereClause)  # 去除异常值
    outSetNull.save(outrasters)
    print(str(ras))
print("All project is OK！")

END 更多精彩尽在GeoLab 219