python_MOD13A3 NDVI数据剔除无效值

最新推荐文章于 2024-08-15 11:17:44 发布

Zoe_左

最新推荐文章于 2024-08-15 11:17:44 发布

阅读量5k

点赞数 3

分类专栏： python 文章标签： python 去除无效值栅格数据

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

版权

python 专栏收录该内容

16 篇文章 4 订阅

订阅专栏

本文介绍了一种利用ArcPy中的SetNull函数去除栅格数据中无效值(-3000)的方法，并通过乘以0.0001进行数值调整。此过程包括检查许可、指定输入栅格文件及无效值条件、执行SetNull操作并保存处理后的栅格。

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

无效值为 -3000，使用setNull函数来执行

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

arcpy.CheckOutExtension("ImageAnalyst") #检查许可
arcpy.CheckOutExtension("spatial") #检查许可
whereClause = "VALUE = -3000" #无效值
inRaster=r'G:\MOD13A3_NDVI\MaxMonth_NDVI\2002_MaxNDVI.tif' #输入栅格数据

outSetNull = SetNull(inRaster, inRaster, whereClause) * 0.0001 #去除无效值并乘以0.0001
outname=r"G:\MOD13A3_NDVI\去除无效值\try1.tif" #输出路径
outSetNull.save(outname) #保存数据
print 'over'