北京市街道边界矢量数据获取方法

小猿猴GISer 1周前

以下文章来源于GIS点滴 ，作者James

GIS点滴

每天进步一点点，让收获的满足感，冲淡压力生活下的浮躁和焦虑，获得真正的专业技术能力！

Hello,大家好，我是James.

我们知道：数据是GIS比较重要的组成部分，没有基础数据的支撑，其他工作无从展开。可能你曾为了一份行政边界数据在群里或者论坛各种跪求，但都石沉大海，膝盖也已长出了老茧，今天我以一个近期遇到的数据爬取为例，手把手教你抓取地图矢量数据。

也是偶然看到天地图北京的官网有行政区、街道办的边界数据

惯性使然打开F12，惊奇地看一串串熟悉的坐标

曾经搜遍论坛网站各个角落，今天在官方却能轻而易举的找到，真是太激动

于是先研究了url请求结构，最终得到了北京16个区和下属街道办的边界信息

接下来跟大家分享抓取过程，

1、首先进入首页（http://beijing.tianditu.gov.cn/bjtdt-mapcenter/）

2、点击【行政区划】，F12我们可以看到在http://beijing.tianditu.gov.cn/tdt-web-sysmanager/services/catalog/special/getLoginOutCatalogSpecials/38ed87a3fe8e4bb1a01949c47baf5244请求结果中包含了北京下属的16个区！

3、展开其中一个后，看到更为详细的信息

4、点击东城区，我发现在http://beijing.tianditu.gov.cn/tdt-web-sdks/biz/httpproxy/httpproxy.jsp?url=http://www.bjmap.gov.cn:8080/dfc/services/sgssfs/2220?request=getfeature&filter=+GBCODE%3D%27110101%27请求结果中有不得了的坐标串，明显是行政区划的边界节点坐标串嘛

5、继续点击建国门街道，同样在http://beijing.tianditu.gov.cn/tdt-web-sdks/biz/httpproxy/httpproxy.jsp?url=http://www.bjmap.gov.cn:8080/dfc/services/sgssfs/2220?request=getfeature&filter=+GBCODE%3D%27110101001%27看到街道办的坐标串，完美！

 最后得到结果如下：

北京市辖区边界

北京市辖区街道办边界

接下来上python代码

1、首先惯例，我使用的是python2,首先声明编码，导入request、json解析模块

2、接下来需要写4个函数方法，

主要是：

（1）解析、构造url地址

（2）组织爬取的信息

3、然后是调取以上函数方法，写入至txt文件

4、最后分别得到行政区和街道办边界节点的txt文本