用Python实现网格计算器中的con和power,栅格,Con,Power


一般我们经常用到栅格计算器中的Con函数,如将插值中的某些误差数据的修改,下面这个是在栅格计算器中的条件公式:

Con("Raster"<0,0,"Raster")

于是,我们就可以直接用Python调用Arcpy来使用这个Con函数
直接使用IDLE运行代码

import os
import arcpy
from arcpy.sa import *
# 导入扩展模块ArcPy Spatial Analyst
arcpy.CheckOutExtension("spatial")
# 检查并授权许可给模块
arcpy.env.overwriteOutput=1
arcpy.env.workspace="E:\\suanli\\0509"
files = arcpy.ListRasters()
print files

for file in files:
    if os.path.splitext(file)[1]=='.tif':
        print file
        out_raster=arcpy.sa.Con(file<0,0,file)
        out_raster.save("E:\\suanli\\0509" + os.sep + "2-1\\"+file)

运行结束后发现栅格并没有按照Con条件公式变化,修改公式

Con("Raster",0,"Raster","value <= 0 ")
# Con(条件栅格,正确时的值,错误时的值,条件表达式)
for file in files:
    if os.path.splitext(file)[1]=='.tif':
        print file
        out_raster=arcpy.sa.Con(file, 0, file, "value <= 0 ")
        out_raster.save("E:\\suanli\\0509" + os.sep + "2-1\\"+file)

成功运行!
同样Power函数也可用同样的方法运行

for file in files:
    if os.path.splitext(file)[1]=='.tif':  
        print file
        out_raster=arcpy.sa.Power(file,2)
        out_raster.save("E:\\suanli\\0509\\2020.0510" + os.sep +"2-2\\"+file)

调用过程中容易出现错误 999998 ,文件存储出错, 修改路径中的中文,恢复默认路径 ,如果还是出现错误,就 直接复制代码到ArcGIS中的Arcpy运行即可