用Python实现栅格计算器中的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