基于arcpy实现去除异常值和缩小值

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

---

前言

提示：这里可以添加本文要记录的大概内容：
下载的数据有效值是1-10000，且需要乘以0.0001恢复原值

---

实现代码

代码如下（示例）：

# -*- coding: utf-8 -*-
import arcpy
from arcpy.sa import *

arcpy.CheckOutExtension("ImageAnalyst")  # 检查许可
arcpy.CheckOutExtension("spatial")  # 检查许可
s="E:\\DATA\\CCI_AP_SM\\monthmean/"
arcpy.env.workspace = s # 工作空间
rasters = arcpy.ListRasters("*", "tif")  # 遍历工作空间中的tif格式数据
a="E:\\DATA\\CCI_AP_SM\\monthmean0.0001\\"
whereClause = "VALUE < 1"   # 去除异常值

# 循环rasters中的所有影像，进行去异常值操作
for ras in rasters:
    outrasters = a + str(ras)  # 更改输出栅格名字
    outSetNull = SetNull(ras, ras, whereClause)# 去除异常值
    ra_DY = Times(outSetNull, 0.0001)
    ra_DY.save(outrasters)
    print(str(ras))
print("All null project is OK！")

总结

自认为比较简单快捷的处理方式，以上只是根据个人需要进行编码，仅供参考，自留保存用。