旧版本｜2007～2016 年税调企业与所在区县税务局的距离

新版本在这里：https://rstata.duanshu.com/#/brief/course/c619040c1b164f8da5d9b58c70e29a5c

之前给大家分享过 2007~2016 年税调企业的地理位置数据：

2007~2016 年税调企业地理信息数据（含经纬度及其所处的省市区县）：https://rstata.duanshu.com/#/brief/course/76d38022cd004b09b2aa09647936beb0

最近有小伙伴想要计算税调企业与所在区县税务局的距离，这不就来了！

税务局经纬度

各区县税务局的经纬度是根据“各区县的名字 + 税务局”使用百度地图地理编码解析得到。

• 税务局经纬度.dta

下图展示了各区县税务局地理分布：

使用百度地图进行地理编码可以学习这两个课程：

1. 使用 Stata 进行地理编码：地址解析经纬度、坐标转换 & 根据经纬度判断所处的省市区县：https://rstata.duanshu.com/#/brief/course/537300af1a9947edb758789785c600f3
2. 使用 R 语言进行地理编码：地址解析经纬度、坐标转换 & 根据经纬度判断所处的省市区县：https://rstata.duanshu.com/#/brief/course/6f8633b486ed49e398c486d0d9f0f59a

税调企业与所在区县税务局距离

解析得到税务局经纬度之后，就可以计算 2007~2016 年税调企业与所在区县税务局的距离了：

• 税调企业与所在区县税务局的距离_km.dta

数据包含如下变量：

sdid、与税务局的大圆距离_km、税务局经度、经度、纬度、省、省代码、市、市代码、县、县代码、税务局纬度、法人代码、企业名称

下图展示了 2007~2016 年各年税调企业与税务局的平均距离变化：

距离的计算在 Stata 中可以使用 geodist 命令，在 R 语言中可以使用 sf 包，对地理计算感兴趣的小伙伴可以学习平台上的这个课程：

如何计算区县质心的地理距离：R 语言和 Stata 的对比使用：https://rstata.duanshu.com/#/brief/course/99d0753092194983833382aca15942d5

点击这里跳转到 RStata 短书平台获取附件：旧版本｜2007～2016 年税调企业与所在区县税务局的距离