arcGIS python两个栅格数据,栅格计算器进行栅格代数计算。
下载链接:多个栅格数据进行代数运算实例demo
信息熵(区位熵)计算公式:
#-*- coding:UTF-8 -*-
import arcpy
from arcpy import env
from arcpy.sa import *env.workspace = r"C:\Users\HP\Desktop\20210323\tweets.gdb"#CheckOut liscen SpatialAnalysis
arcpy.CheckOutExtension("ImageAnalyst")
arcpy.CheckOutExtension("Spatial")
env.overwriteOutput = 1outDACA = r"C:\Users\HP\Desktop\20210323\tweets.gdb\temp_DACA_KDE"outMudslide = r"C:\Users\HP\Desktop\20210323\tweets.gdb\temp_mudslide_KDE"
meanValue =0.5
z=0.00000001resultMap = Log2(((DACA_Raster + z)/(Mudslide_Raster + z))* 1/meanValue)
result_KDE = r"C:\Users\HP\Desktop\20210323\tweets.gdb\result_KDE"
resultMap.save(result_KDE)
print("---output result KDE Successful-----")