处理NDVI数据SetNull函数使用

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")
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值