NDVI数据剔除无效值/赋为空值

# -*- coding: utf-8 -*-
import arcpy
from arcpy import env
from arcpy.sa import *
 
arcpy.CheckOutExtension("ImageAnalyst") #检查许可
arcpy.CheckOutExtension("spatial") #检查许可
env.workspace = r"F:\EVI_TEST" #数据存放目录
 
outfile=r'F:\output' #输出数据存放路径
 
whereClause = "VALUE < 0.1" #设空函数条件
 
for file1 in arcpy.ListFiles("*.tif"): #循环列表下的tif文件
    s=file1 #原始文件名
    print s
    C_outfile=outfile+'\\'+ s  #输出的文件名，完整路径
    print C_outfile
    inRaster= Raster(s)
 
    outSetNull = SetNull(inRaster, inRaster, whereClause)   # value<0.1  的值设为空值. 
 
    outSetNull.save(C_outfile) #保存数据
 
    print 'over'
    

转载：https://blog.csdn.net/u012131430/article/details/82779231