1998~2014 年工企所处的东西部属性及其与胡焕庸线的距离

之前给大家讲解过如何使用 R 语言判断工企所处的南北方属性及计算其与秦岭淮河线的距离，并且分享了计算好的数据：

如何判断工企业在南方还是北方 & 计算每个工企业距离秦岭淮河线的距离？：https://rstata.duanshu.com/#/course/8a4e1e28c1e54c30a7d6996823648242

最近有培训班的小伙伴问是否也可以计算下工企所处的东西部属性及其与胡焕庸线的距离，这不就来了！

数据概览

通常我们使用胡焕庸线划分东西部，胡焕庸线也即黑河-腾冲线，是我国地理学家胡焕庸在 1935 年提出的划分我国人口密度的对比线，例如之前使用 R 语言绘制的这幅图人口密度图：

该图来源于课程：https://rstata.duanshu.com/#/brief/course/8f921af271a04b6db79c52ce48a04d49

数据格式方面，我提供了 dta 格式的工企的东西部及其与胡焕庸线的距离，时间范围是 1998～2014 年，经过计算之后一共得到了 471.8 万条结果：

包含如下变量：

• gqid
• 经度
• 纬度
• 东西部
• 与胡焕庸线的距离_km

下图以 2014 年为例，展示了中国工业企业与胡焕庸线的距离：

处理方法

通常我们使用胡焕庸线划分东西部，不过胡焕庸线是线条元素，不能用于经纬度的判别。所以我们需要生成东西部区域的多边形。一种思路就是先生成中国范围的方框，然后延长胡焕庸线把这个方框切分成两个区域：

使用胡焕庸线划分东西之后，就可以对工企进行东西部分类了。首先把工企地理位置数据转换成 sf 对象，再计算工企距离胡焕庸线的距离，就可以得到附件中的工企的东西部及其与胡焕庸线的距离.dta了。

对处理过程感兴趣的小伙伴可以学习下面的课程：

1. 使用 R 语言计算胡焕庸线两侧的省市区县（含结果）：https://rstata.duanshu.com/#/course/657a959f57ce4f9a82336994691a5eb3
2. Stata：如何根据工企经纬度判断所处的南北方、东西部：https://rstata.duanshu.com/#/course/0fbc213d6d2e405bbce9d7e7a9557801

附件中也提供了详细的代码以供参考～

点击这里跳转到 RStata 短书平台获取附件：1998~2014 年工企所处的东西部属性及其与胡焕庸线的距离