利用Python实现ARCGIS中栅格计算器的简单功能

1.加法操作

import arcpy
from arcpy import env

env.workspace = "E:\\work\\cppservice\\RasterCalc" #工作空间可以设置也可不设置
inRaster1 = "E:\\work\\cppservice\\finalResult\\propertyGower.tif"
inRaster2 = "E:\\work\\cppservice\\finalResult\\propertyGowerCsharp.tif"
arcpy.CheckOutExtension("3D")
arcpy.Plus_3d(inRaster1,inRaster2,"E:\\work\\cppservice\\RasterCalc\\solim-csharp.tif")

2.减法操作

import arcpy
from arcpy import env

env.workspace = "E:\\work\\cppservice\\RasterCalc" #工作空间可以设置也可不设置
inRaster1 = "E:\\work\\cppservice\\finalResult\\propertyGower.tif"
inRaster2 = "E:\\work\\cppservice\\finalResult\\propertyGowerCsharp.tif"
arcpy.CheckOutExtension("3D")
arcpy.Minus_3d(inRaster1,inRaster2,"E:\\work\\cppservice\\RasterCalc\\solim-csharp.tif")

3.乘法
修改一下最后一行的函数即可

arcpy.Times_3d()

4.当一个数除以零时，输出结果为 NoData。
输入的数据类型将决定输出的数据类型：
如果两个输入均为整数，则将执行整数除法，并且输出结果为整数。例如，如果 3 除以 2，则输出为 1。
如果两个输入中的任何一个为浮点型，则将执行浮点型除法，并且输出结果也为浮点型值。例如，如果 3 除以 2.0，则输出为 1.5。

arcpy.Divide_3d()

参考自：http://wenku.baidu.com/link?url=Ei4IxTmqbshZ5nBYa962sjhuMVKfa1UY3o9ADGPXWXwjaLXZ-Wz0auRsHkNr2FaZ5CbfYuOCz1FmvRaK4EACAG-fCkIaqqdbl4VkW5-87TW