python_MOD11A2 地表温度数据 剔除缺省值和无效值

   输入数据为使用MRT 拼接并转换投影后的tif文件

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

arcpy.CheckOutExtension("ImageAnalyst") #检查许可
arcpy.CheckOutExtension("spatial") #检查许可
env.workspace = r"G:\MOD11A2_LandSurfaceTem\data\ExtractData\2002\200201" #数据存放目录

outfile=r'G:\MOD11A2_LandSurfaceTem\data\delete_NullData\LSTDay_Nodata\2002\200201' #输出数据存放路径

whereClause = "VALUE < 7500 OR VALUE > 65535" #设空函数条件

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.02-273.15   # value>=65535 or <=7500  的值设为空值. T(℃)=K-273.15  0.02为缩放因子
 
    outSetNull.save(C_outfile) #保存数据

    print 'over'
    

  输出数据值范围，单位为摄氏度 ℃