查看原文
其他

【GIS基础】城市功能混合程度计算

十六便士 数读城事 2022-04-25
絮絮叨叨:

不想干正事的时候就想起了欠的债,虽然这个真的没有技术含量,但是还是欢迎大家找我讨论互相学习!

——投稿人的垃圾话

数读菌插句话:你这匿名谁知道你是谁==怎么找你讨论互相学习???


 
01混合度计算◐ 

1.1 概念与计算

熵:信息论中度量随机事件在某项实验中的不确定程度的概念。

计算公式:

H(X)表示随机变量X的熵; Pi为X取Xi的概率 。显然 , 熵值越大 , 不肯定性越大;熵值越小 , 不肯定性越小 。


◐ 

1.2 参考案例

我们可以依据熵来计算地块中城市用地的混乱程度,即混合程度。

例如:城市中有两块地域:地域A内 55%是居住用地,25%是商业用地,20%是工业用地;地域B内35%是商业用地,65%是居住用地。

可以用信息熵(X)来比较两块地的混合程度:

  H(A)= -0.55lg0.55 -0.25lg0.25 -0.20lg0.20 ≈ 0.433

  H(B)  = -0.35log0.35 -0.65log0.65 ≈ 0.281

  H(A)> H(B)

熵值越高,混合度越高 — 根据熵值显示 , 地域I混合程度高。

缺少用地数据的时候,我们可以用poi的个数来代替。

 

02操作过程


以下以上海市杨浦区为例,演示操作过程。

◐ 

2.1 加载数据

导入上海市杨浦区的行政边界、路网数据,及poi点数据(投影后)。为方便计算,这里只计算六类公共服务设施:办公、金融、休闲、住宿、餐饮、购物(没有为什么,随便选的,嘻嘻)。 

图1基础数据


◐ 

2.2 分割地块

编辑杨浦行政区的图层,选中所有道路,点击“分割面”将杨浦区分割成地块(可以先单独导出一个文件命名为地块方便后续操作)。

  

图2|分割地块操作示意


◐ 

2.3 计算poi个数

先在地块要素中添加六个字段(可批量),再分别空间连接每类设施的poi至地块要素中,计算每个地块中各类设施poi的个数。(这里我一直想找到批量连接并计算的方法,未果,求大佬教教我)

图3|添加字段操作示意

图4|空间连接计算各地块poi个数


◐ 

2.4 计算poi比例

批量添加字段后,先计算每个地块中所有poi个数的总和,再分别计算各地块中,每类poi的占比。(有很多总和=0的地块是道路,计算比例的时候我们选择总和>0的进行计算即可)

图5|空间连接计算各地块poi个数

便于后续计算,令比例=0的数都计算为一个极小值0.000001。 

图6|计算字段


◐ 

2.5 计算混合度

最后依据公式计算熵值,即可得到各地块的混合程度分布。结果在此就不分析了,设施种类不齐全,也未进行数据清洗、整理,同时还欠缺考虑规模等因素,因此不保证结果的科学性和准确性。

图7|计算混合度

图8|混合度分布示意

另外,在缺少路网数据的时候,或者想做的更精细的话,可以用渔网作为底图,方法如上。

9混合分布示意基于网格最后,还可以进一步进行空间相关性分析,如,聚类和异常值分析。图10|聚类和异常值分析


作者|十六便士

排版|Moon

校阅|数读菌



那今天就到这里结束啦,欢迎留言讨论。文中的图片文字未经许可不要随便“引用”。

如果可以的话,希望能够转发分享,点个在看并且点个,给个赞赏~~也欢迎规范转载~

也希望大家和我多留言互动啊!(据说这样可以增加我的推送在你的订阅号里出现的概率)

历史文章推荐:

【瞎七八扯】如何确定“真实”人口数量、分布以及流动?

【瞎七八扯】灯塔国全球军事基地的分布

【数据分享】中国生态服务功能重要性与中国生态系统敏感性

【数据整理】2019年行政区划调整的数据更新与分享

【数据分享】如何获取北京旅游景区AOI数据?

【数据分享】如何获取全国矢量地铁站点与线路数据?

【瞎七八扯】全国同名的行政区划:市县同名(修订增补)

【瞎七八扯】仅仅是“巧合”的同名行政区划

需要你的“分享”和“在看”


END>

如需全文转载文章、投稿或者合作

可添加微信

(回复超慢!!!)

(不要添加我问各种问题,我大概率不会的==)

(入群请一定要备注入群)

(添加后会在晚上非工作时间通过,请稍安勿躁)


公众号


微博

▼ 点击阅读原文,使用关键词搜索历史文章


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存