用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 filesfor file in files:if os.path.splitext(file)[1]=='.tif':print fileout_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 fileout_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 fileout_raster=arcpy.sa.Power(file,2)out_raster.save("E:\\suanli\\0509\\2020.0510" + os.sep +"2-2\\"+file)

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

用Python实现栅格计算器中的Con和Power相关推荐

python实现栅格计算器功能
对于Python不能直接调用arcgis中的栅格计算器工具.但是可以通过使用Raster函数实现. 如下面代码实现了每年12个月做一个累加,然后输出一个年值的功能: for year in range ...
栅格计算器函数之Con
Con函数是condition(条件)的缩写,其作用语法是: Con(条件,条件为真执行语句,条件为假执行语句[可选]) 或Con(输入栅格,条件为真执行语句,条件为假执行语句[可选],逻辑表达式) ...
Python实现ARCGIS栅格计算器con函数功能
目的现有某地区土地利用类型图,图中不同的数值代表不同的土地利用类型,如像元值为20的,我们重新赋值为0,像元值为123的,我们重新赋值为15等等. 图1 原始图像图2 结果图像代码如下 from ...
arcgis利用栅格计算器修改DEM高程数据
有时候高程数据有误,需要修正,我们可以使用栅格计算器来进行修改前在ArcScene中如下图 [img]http://dl2.iteye.com/upload/attachment/0098/3382 ...
ArcGIS中栅格计算器常用函数的使用
先看下目录一. ArcGlS中栅格计算器简介栅格计算器位置及界面布局二. ArcGIS棚格计算器中Con()函数的使用 (一) Con()函数的用法 (二) ArcGIS 案例操作三. Arc ...
arcgis栅格计算器python教程_arcgis 栅格计算器（Spatial Analyst/Raster Calculator）
栅格计算器中用得到$$相关函数 $$NROWS: the number of rows in the analysis window (行数) $$NCOLS: the number of colum ...
arcgis 栅格计算器（Spatial Analyst/Raster Calculator）
原文地址:点击打开链接栅格计算器中用得到$$相关函数 $$NROWS: the number of rows in the analysis window (行数) $$NCOLS: the num ...
arcgis 栅格数据邻域计算_ArcGIS栅格计算器小知识
1.栅格计算器中,表达式中的图层名要用[]括起来,运算符前后都要加空格,不然会出错. 2．Con命令命令使用方法有两种: 一是使用栅格计算工具(Spatial Analyst/Raster Calc ...
ArcGIS 栅格计算器con（）函数用法
函数的基本语句:con(条件(语句),真(语句),假(语句)) 若"真(语句)"或"假(语句)"存在缺失,则对应栅格单元的值赋值为空(NoData). 三处语句 ...

用Python实现栅格计算器中的Con和Power

用Python实现栅格计算器中的Con和Power相关推荐

最新文章

热门文章