arcGIS python两个栅格数据,栅格计算器进行栅格代数计算。
信息熵(区位熵)计算公式:
#-*- 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 = 1
outDACA = 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.00000001
resultMap = 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-----")