[特殊字符] 第十四讲 | 空间异质性检验与地统计局部指标(LISA)应用

📌 关键词:空间异质性、LISA、局部Moran's I、空间聚集、冷热点分析、GeoDa、R语言


🧠 导语:空间现象为何"不一样"?

在地理学与农学研究中,我们经常遇到"某地污染严重,而邻近区域却很轻微"的现象。这种空间格局的不均匀性,正是**空间异质性(Spatial Heterogeneity)**的体现。而揭示这些"异质"的核心利器,就是------LISA(Local Indicators of Spatial Association)局部空间自相关分析


🔍 一、空间异质性是什么?

空间异质性指的是某种地理变量(如土壤肥力、污染浓度等)在地理空间中呈现不规则分布,具体体现在:

  • 位置差异性(Spatial Discontinuity)

  • 局部聚集性(Local Clustering)

  • 不稳定的空间关系

我们需要借助LISA进行定量刻画!


🧪 二、LISA(局部空间自相关)指标详解

LISA 指标衡量的是某个点与其邻居的空间关系是否显著偏离"空间随机性",主要包括:

指标 含义
局部Moran's I 衡量某位置值与邻居值的相似性
局部Getis-Ord G* 用于发现局部"高值"或"低值"集中区域(热点分析)
局部Geary's C 类似Moran's I,强调差异性

🧰 三、R语言实战:LISA分析步骤

📦 所需R包

复制代码
library(spdep)
library(sf)
library(tmap)

🧮 示例代码流程

复制代码
# 读取空间数据
shp <- st_read("your_shapefile.shp")

# 构建邻接矩阵
nb <- poly2nb(shp)
lw <- nb2listw(nb, style="W")

# 计算局部Moran's I
local_moran <- localmoran(shp$target_variable, lw)

# 将结果加入数据框
shp$Ii <- local_moran[,1]
shp$p_value <- local_moran[,5]

# 可视化:显著性图
tmap_mode("view")
tm_shape(shp) +
  tm_fill("Ii", style="quantile", palette="-RdBu", title="Local Moran's I") +
  tm_borders()

📌 四、结果解读:四象限含义

LISA 结果常被划分为 4 个象限:

象限 类型 含义
I象限 高-高(HH) 高值区域被高值包围(热点)
II象限 低-低(LL) 低值被低值包围(冷点)
III象限 高-低(HL) 异常值,高值被低值包围(离群)
IV象限 低-高(LH) 异常值,低值被高值包围(离群)

🧑‍🔬 五、应用场景推荐

  • 🧫 土壤污染高风险区识别

  • 🌾 农作物高产区空间聚类分析

  • 🌍 区域经济发展空间不均衡性检验

  • 🏙 城市热岛效应的热点识别


✍️ 学术论文写作句式参考

"利用局部Moran's I指标对土壤有机碳含量进行空间聚集性分析,结果显示研究区呈现明显的高-高聚集特征,主要集中于中部农田密集区。"


📚 总结

✅ LISA 是空间分析中极具辨识力的工具

✅ 通过空间聚集和离群识别,助力精准干预决策

✅ 搭配GeoDa软件或R语言工具可轻松实现全流程分析


📌 下一讲预告 :空间面板数据建模与时间-空间交互分析

📎 如果你喜欢这个系列,欢迎收藏 + 点赞 + 评论互动~

如需本节课的R脚本 + 示例数据 + 图形PDF,可留言获取。

相关推荐
杨超越luckly3 小时前
HTML应用指南:利用POST请求获取全国中国工商银行网点位置信息
前端·arcgis·html·银行·数据可视化·网点
大翻哥哥15 小时前
Python地理空间数据分析:从地图绘制到智能城市应用
开发语言·python·数据分析
Lucky_Turtle18 小时前
【electron】一、安装,打包配置
javascript·arcgis·electron
用户Taobaoapi201421 小时前
淘宝店铺所有商品详情API技术文档
大数据·数据挖掘·数据分析
GIS小小研究僧1 天前
同一个栅格数据,为何在QGIS和ArcGIS Pro中打开后显示的数值范围不同?
arcgis·gis·qgis·rs
南飞测绘视界1 天前
【实操教学】ArcGIS 如何进行定义坐标系
arcgis·坐标系·坐标系定义
wanzhong23331 天前
ArcGIS学习-12 实战-综合案例
学习·arcgis
医工交叉实验工坊1 天前
R 语言 ComplexUpset 包实战:替代 Venn 图的高级集合可视化方案
算法·r语言
青春不败 177-3266-05201 天前
R语言贝叶斯方法在生态环境领域中的高阶技术应用
r语言·贝叶斯·生态学·科研绘图·结构方程·环境科学·混合效应
星空1 天前
网络安全监控中心
安全·web安全·arcgis·地图·前段