ARCGIS栅格计算器SetNull Python脚本应用

JH_Kong

已于 2022-04-16 18:02:00 修改

阅读量2.2k

点赞数 3

文章标签： python

于 2022-04-09 22:21:00 首次发布

本文链接：https://blog.youkuaiyun.com/qq_33524950/article/details/124069483

版权

这篇博客介绍了如何使用Python的arcpy模块批量处理栅格数据，将所有值为-999的栅格单元设置为NoData。通过提供一段代码示例，展示了如何运用RasterCalculatorSa函数实现这一操作，适用于需要对大量栅格文件进行相同修改的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

论文需要把栅格数据某些不符合条件的栅格设为NoData，但是要批量操作，去网上搜了一顿终于解决了这个问题，下面附上代码（我的需求是把-999的栅格设为Nodata）

import arcpy
from arcpy.sa import *
arcpy.env.workspace = "你的栅格文件路径"

#搜索路径中的栅格文件
rasterlist = arcpy.ListRasters("*","img/tiff")
output_path = "栅格文件输出路径"

arcpy.CheckOutExtension("Spatial")

for raster in rasterlist:
    raster1 = str(raster)
    inraster = "你的栅格文件路径" + raster1
    expression = 'SetNull(' + "'"  + inraster + "'" + '==-999,' + "'" + inraster + "'" + ')'
    arcpy.gp.RasterCalculator_sa(expression, output_path + raster1[:-4] + ".tif")

需要在python 2.x中运行arcpy。

如果有用就请点个赞吧！！