利用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

转载自:https://blog.csdn.net/hnyzwtf/article/details/44850981

You may also like...