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

Zoe_左

于 2018-09-19 22:09:11 发布

阅读量8.2k

点赞数 5

分类专栏： python 文章标签： python MOD11A2 去除无效值和缺省值

本文链接：https://blog.csdn.net/u012131430/article/details/82779231

版权

本文介绍如何使用Python处理MOD11A2地表温度数据，通过MRT工具进行数据拼接和投影转换后，进一步去除数据中的缺省值和无效值，确保输出数据的值范围在摄氏度℃之间。

摘要由CSDN通过智能技术生成

输入数据为使用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'

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