之前给大家分享过 2007~2016 年税调企业的地理位置数据：
2007~2016 年税调企业地理信息数据（含经纬度及其所处的省市区县）：https://rstata.duanshu.com/#/brief/course/76d38022cd004b09b2aa09647936beb0
最近有小伙伴想要计算税调企业与所在区县税务局的距离，这不就来了！
税务局经纬度
各区县税务局的经纬度是根据“各区县的名字 + 税务局”使用百度地图地理编码解析得到。
税务局经纬度.dta
下图展示了各区县税务局地理分布：
使用百度地图进行地理编码可以学习这两个课程：
使用 Stata 进行地理编码：地址解析经纬度、坐标转换 & 根据经纬度判断所处的省市区县：https://rstata.duanshu.com/#/brief/course/537300af1a9947edb758789785c600f3 使用 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 培训班获取全部课程和以会员价获取数据资料（10元/份）详情可阅读这篇推文：数据处理、图表绘制、效率分析与计量经济学如何学习～
详情可点击阅读原文进入 RStata 学院了解（从首页的会员卡专区即可查看和购买会员卡）。
更多关于 RStata 培训班的信息可添加微信号 r_stata 咨询：
附件下载（点击文末的阅读原文即可跳转）：
https://rstata.duanshu.com/#/brief/course/8719f12a1c5e47119e5d9e1fb1aa9220
