import arcpy
import math
import os
from arcpy import env
from arcpy.sa import *
## 检验模块是否激活
arcpy.CheckOutExtension("Spatial")
arcpy.env.workspace = "F:/data/hhy/hhy_NDVI/hhy/extract_1982-2015"
## 输入文件路径
in_path = "F:/data/hhy/hhy_NDVI/hhy/extract_1982-2015/"
ext_name = "tif"
## 遍历输入文件夹下的所有文件
for in_file in os.listdir(in_path):
if in_file.endswith(ext_name):
f_name = os.path.splitext(in_file)
print(f_name)
in_raster = in_path + f_name[0] + ".tif"
print(in_raster)
out_path_name = "F:/data/hhy/hhy_NDVI/hhy/SetNull_out/" + str(f_name[0]) + ".tif"
print(out_path_name)
outSetNull = SetNull(in_raster, in_raster, "VALUE < 0") * 0.0001
"""
返回值:条件为真,返回NoData。条件为假,返回第二个参数。
第一个参数:输入文件
第二个参数:条件为假,则输出第二个参数
第三个参数:逻辑表达式
"""
outSetNull.save(out_path_name)
print("finish")
print("All finish")
处理NDVI数据SetNull函数使用
最新推荐文章于 2024-06-16 14:00:00 发布