NDVI批量處理排除值-3000

代碼如下：

import arcpy
from arcpy import env
from arcpy.sa import *
import sys
reload(sys)
sys.setdefaultencoding('utf8')

arcpy.CheckOutExtension("ImageAnalyst")  # 檢查許可
arcpy.CheckOutExtension("spatial")  # 檢查許可
env.workspace = r"D:\NDVI\NDVI樣地\MOSAIC2016"  # 數據存放目錄

outfile = r'D:\NDVI\排除值\paichu2016'  # 輸出數據存放路徑

whereClause = "VALUE = -3000"  # 設空函數條件

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) * 0.0001  # 去除無效值並乘以0.0001

    outSetNull.save(C_outfile)  # 保存數據

    print ('over')

 

 

這樣就可以批量處理MODIS下載下來的-3000的無效值，同時可以乘以0.0001

具體前一步參考青燈常伴古佛的博客：

https://www.cnblogs.com/9587cgq/p/12671881.html