arcGIS python两个栅格数据，栅格计算器进行栅格代数计算。信息熵（区位熵）计算

原创已于 2022-11-02 21:59:47 修改 · 4k 阅读

13 ·

CC 4.0 BY-SA版权

文章标签：

#gdb #arcgis #数据库

于 2021-03-23 16:01:19 首次发布

python 同时被 3 个专栏收录

55 篇文章

订阅专栏

arcpy

37 篇文章

订阅专栏

arcgis

30 篇文章

订阅专栏

本文介绍如何使用Python与ArcGIS进行栅格数据的代数运算，如计算信息熵。通过实例演示了如何加载数据，应用Log2函数并保存结果。适合GIS和Python开发者学习高级分析技巧。

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 = 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-----")