最近有小伙伴需要计算各城市市辖区的空间形态紧凑度指标，该指标的计算原理是这样的（方颖、白秀叶，2022）：
除了市辖区的，我还分别计算了 2001～2021 年各省份、各城市和各区县的空间形态紧凑度指标。这里用到的中国土地覆盖数据是之前分享过的：
2001~2021 年各省市区县不同土地覆被类型的区域面积：https://rstata.duanshu.com/#/course/9f6a61dfdc0948b1b30b1e272de7e766
其中，土地覆被类型 13 代表城市和建成区。
数据概览
附件中包含如下文件：
2001~2021年各城市空间形态紧凑性指标面板数据.dta 2001~2021年各城市市辖区空间形态紧凑性指标面板数据.dta 2001~2021年各区县空间形态紧凑性指标面板数据.dta 2001~2021年各省份空间形态紧凑性指标面板数据.dta
除了省市区县、年份信息，各文件都包含了 avg_distance_km(建成区内任意两点之间距离的均值)、adk_ec(建成区等面积圆内任意两点之间距离均值)、spatial_compactness(城市空间形态紧凑度)变量。以 2001~2021年各区县空间形态紧凑性指标面板数据 为例，该数据包含了 58568 条观测值，预览如下：
绘图展示
为了更直观地展示空间形态紧凑度指标，我还绘制了下面这些图。
各城市空间形态紧凑度
下图展示了 2021 年中国各城市空间形态的紧凑度：
各区县空间形态紧凑度
下图展示了 2021 年中国各区县空间形态的紧凑度：
连续性检查
2001~2021 年各城市年空间形态紧凑度变化展示如下：
连续性还不错！
不同来源土地数据计算结果对比
之前我还提供过一份根据中国土地利用数据（1980-2015）：
1980~2015 年各省市区县不同土地利用类型的区域面积：https://rstata.duanshu.com/#/course/0b37661545e1491e8e3887466163a1a7
计算得到的空间形态紧凑度结果：
1980~2015 年各省市区县空间形态紧凑度指标面板数据：https://rstata.duanshu.com/#/course/9e6bacd816324d5ea15004cd0d4c15e3
本文采用了另一种来源的土地覆盖数据，两种不同来源的土地数据计算得到的各城市空间紧凑度结果对比如下：
处理方法
对空间形态紧凑度感兴趣的小伙伴可以学习下面这个课程：
名师讲堂｜使用 R 语言计算各城市空间形态的紧凑度：https://rstata.duanshu.com/#/course/5ebfe03c91774471b483104e9477e70f
获取数据
是不是感觉很硬核！欢迎报名 RStata 培训班获取全部课程和以会员价获取数据资料（10元/份）详情可阅读这篇推文：数据处理、图表绘制、效率分析与计量经济学如何学习～
详情可点击阅读原文进入 RStata 学院了解（从首页的会员卡专区即可查看和购买会员卡）。
更多关于 RStata 培训班的信息可添加微信号 r_stata 咨询：
附件下载（点击文末的阅读原文即可跳转）：
https://rstata.duanshu.com/#/brief/course/fd68e2b4b73749c8862dcf399ac63f84
