[关闭]
@Yano 2018-02-06T08:27:44.000000Z 字数 5853 阅读 2304

爬取北京摩拜单车信息(附分析过程和详细代码)

爬虫


目的

统计北京摩拜单车的数量、分布、类型等信息。

先贴出结果:

分析过程

Fiddler 抓包

去年11月份可以在微信中抓取摩拜的小程序,但是现在不行了。当时微信小程序的API很简陋,利用代理可以直接抓取。但是现在试下挂上代理小程序都打不开了。

后来使用摩拜的APP抓取成功,附图。

Fiddler 截图1

获取附近单车信息的接口是:http://app.mobike.com/api/nearby/v3/nearbyBikeInfo

同时 response 是个 json,解析出来是这样的:

Fiddler 返回值1

其中包括 bike 和 mpl 数组,bike 即为最终要获取的单车详细信息,结构如下:

Fiddler 返回值2

而 mpl 应该是摩拜推出的新功能,表示「停车地点」,在这些地点停车有机会获得 NIKE 新春贺礼 85 折优惠券,结构如下:

Fiddler 返回值3

经纬度查询

在 Fiddler 中抓取 request 发现纬度和经度分别为:39.977524,116.305863,通过经纬度查询的网站发现定位极其精准:

定位1

查询 response 中 mpl 信息的一个地址:39.975823, 116.305252

发现确实是地图上的位置:

Postman 验证

把请求参数放到 postman 中,发现能够直接给出结果:

通过调试发现 post 的几个重要参数作用如下:

  1. scope:爬取的范围,下面会分析;
  2. longitude:经度,可以精确到小数点后7位;
  3. latitude:纬度;
  4. bikenum:获取单车信息的最大数量,这个值设置大一点比较好,爬取更快;
  5. biketype:0应该是爬取所有类型。

经纬度与距离的换算

经纬度在不同地区,每度距离差是不同的,任意两点距离计算公式为

d=111.12cos{1/[sinΦAsinΦB十cosΦAcosΦBcos(λB—λA)]}

其中A点经度,纬度分别为λA和ΦA,B点的经度、纬度分别为λB和ΦB,d为距离.

所以可知每度大概为111千米

比如北京 B = 40、cosB = 0.766,经度变化1度,则东西方向距离变化 85.567km。

确定爬取范围

决定爬取北京五环内的摩拜单车数量,确定经纬度:

北四环西北 39.9777778168,116.2875366211
北四环东南 39.8401771328,116.4787673950

确定步长、单次爬取单车数量

这两个参数是相互影响的,单次步长越小,结果越精确,但是爬取时间就会更长;增大步长,增加单次爬取单车数量可以缩短时间,但是单位面积的单车数量并不是无限增大的。通过测试发现 scope 参数也不是可以无限增大,貌似到500就无效了。

我确定这两个参数的方法:在地图上找2个点,固定经度调节纬度,看两者爬取的单车是否有重复。

测试了4个地点:

十字路口1:
longitude=116.30581114680909&latitude=39.98106569393767

十字路口2:
longitude=116.30613220306786&latitude=39.978998558907335

维亚大厦:
longitude=116.30608720311963&latitude=39.97747017430741

苏州街十字路口:
longitude=116.30627577120181&latitude=39.975675035347415

发现纬度差0.003时,两次爬取结果会有33辆单车的重叠(每个点大约有180辆单车的信息)。确定经纬度移动的步长是0.003,而单次爬取单车数量为500(99%的地点都到不了500个单车)。

header 参数

发现 header 有很多参数,甚至还有 time 和 utctime 两个时间参数。摩拜会存取每个时刻单车的静态信息吗?显然不会,这样的存储成本太高了,通过程序验证确实是这样。对于其它参数,也可以自行验证。

代码

主体代码很简洁,就是 post 发送请求,每次在地图上移动0.003的经度或纬度,直至遍历完地图上的正方形。

  1. @Test
  2. public void testGetMoBike() throws IOException {
  3. OkHttpClient client = new OkHttpClient();
  4. RequestBody body = RequestBody.create(GSON, "{}");
  5. double longitudeMin = 116.2875;
  6. double longitudeMax = 116.4787;
  7. double latitudeMin = 39.8402;
  8. double latitudeMax = 39.9778;
  9. double step = 0.003;
  10. for (double longitude = longitudeMin; longitude <= longitudeMax; longitude += step) {
  11. for (double latitude = latitudeMin; latitude <= latitudeMax; latitude += step) {
  12. // 获取对应经纬度的结果
  13. List<BikeDo> bikeDos = getBikeDos(client, body, longitude, latitude);
  14. // 将每辆单车信息都存入redis
  15. restoreBikeInfo(bikeDos);
  16. // 随机休眠,防止被封
  17. sleepRandomSeconds(2);
  18. System.out.println("经度:" + longitude + " 纬度:" + latitude + " 结束,爬取总个数:" + jedis.hlen(MAP_KEY));
  19. }
  20. }
  21. }

使用 redis 的 map 存储结果并去重,并打印出每一步的信息。这样可以在断网等异常情况时,从异常处重新开始爬取。redis 在本机可以达到 10万次/s 的连接量,就没必要在 Java 中维护这些代码了。

其余代码如下(删除了个人相关的信息,所以代码不完整):

  1. private static Jedis jedis = RedisUtil.getJedis();
  2. private static final String MAP_KEY = "mobike";
  3. private static final Random RANDOM = new Random();
  4. private void restoreBikeInfo(List<BikeDo> bikeDos) {
  5. for (BikeDo bikeDo : bikeDos) {
  6. // 将每一个 bike 存入 redis 中
  7. if (jedis.hexists(MAP_KEY, bikeDo.getDistId())) {
  8. continue;
  9. }
  10. jedis.hset(MAP_KEY, bikeDo.getDistId(), JSON.toJSONString(bikeDo));
  11. }
  12. }
  13. private void sleepRandomSeconds(int seconds) {
  14. try {
  15. Thread.sleep(RANDOM.nextInt(seconds) * 1000);
  16. } catch (InterruptedException e) {
  17. e.printStackTrace();
  18. }
  19. }
  20. private List<BikeDo> getBikeDos(OkHttpClient client, RequestBody body, double longitude, double latitude) throws IOException {
  21. Request request = createRequest(body, longitude, latitude);
  22. Response response = client.newCall(request).execute();
  23. String responseJson = response.body().string();
  24. JSONArray parse = (JSONArray) JSONObject.parseObject(responseJson).get("bike");
  25. return JSON.parseArray(parse.toJSONString(), BikeDo.class);
  26. }
  27. private Request createRequest(RequestBody body, double longitude, double latitude) {
  28. return new Request.Builder()
  29. .url("http://app.mobike.com/api/nearby/v3/nearbyBikeInfo?scope=700&sign=&client_id=android&biketype=0&longitude=" + longitude + "&latitude=" + latitude + "&userid=&bikenum=600")
  30. .post(body)
  31. .addHeader("version", "6.7.1")
  32. .addHeader("versioncode", "1649")
  33. .addHeader("platform", "1")
  34. .addHeader("mainsource", "4002")
  35. .addHeader("subsource", "3")
  36. .addHeader("os", "25")
  37. .addHeader("lang", "zh")
  38. .addHeader("time", DateTime.now().getMillis() + "")
  39. .addHeader("country", "0")
  40. .addHeader("eption", "685f8")
  41. .addHeader("deviceresolution", "1080X1920")
  42. .addHeader("utctime", DateTime.now().getMillis() / 1000 + "")
  43. .addHeader("mobileno", "")
  44. .addHeader("accesstoken", "")
  45. .addHeader("uuid", "")
  46. .addHeader("longitude", longitude + "")
  47. .addHeader("latitude", latitude + "")
  48. .addHeader("content-type", "application/x-www-form-urlencoded")
  49. .addHeader("content-length", "203")
  50. .addHeader("host", "app.mobike.com")
  51. .addHeader("accept-encoding", "gzip")
  52. .addHeader("user-agent", "okhttp/3.9.1")
  53. .addHeader("connection", "Keep-Alive")
  54. .addHeader("cache-control", "no-cache")
  55. .addHeader("postman-token", "fc17ae5a-32e9-d2bb-eb1f-43dd6760aea6")
  56. .build();
  57. }

bike 类:

  1. public class BikeDo {
  2. private String distId;
  3. private double distX;
  4. private double distY;
  5. private int distNum;
  6. private int distance;
  7. private int biketype;
  8. private int type;

结果

大概1个半小时,爬取到了345679个摩拜单车。

摩拜单车型号[1 x 248627, 2 x 96840, 3 x 2, 999 x 210]:

分布热力图

使用百度地图的API。

经纬度模糊

因为爬取到单车的地理位置太精确了,在地图上画热力分布精确到小数点后5位即可。所以需要统计在每个经纬度下的单车数量,并将点打印出符合百度地图API的格式(统计时将经纬度合并成一个String,否则复杂度太高)。

  1. @Test
  2. public void testMobikeHeat() {
  3. String REDIS_KEY = "mobike";
  4. Map<String, String> bikesMap = jedis.hgetAll(REDIS_KEY);
  5. Multiset<String> multiset = HashMultiset.create();
  6. for (String key : bikesMap.keySet()) {
  7. BikeDo bike = JSON.parseObject(bikesMap.get(key), BikeDo.class);
  8. String mixKey = String.format("%.6f", bike.getDistX()) + String.format("%.6f", bike.getDistY());
  9. multiset.add(mixKey);
  10. }
  11. File file = new File("mobike.txt");
  12. BufferedWriter out = null;
  13. try {
  14. out = new BufferedWriter(new OutputStreamWriter(new FileOutputStream(file, true)));
  15. for (String ele : multiset.elementSet()) {
  16. out.write("{\"lng\":"+ele.substring(0, 10)+",\"lat\":"+ele.substring(10)+",\"count\":"+multiset.count(ele)+"}," + "\r\n");
  17. }
  18. } catch (Exception e) {
  19. e.printStackTrace();
  20. } finally {
  21. try {
  22. out.close();
  23. } catch (IOException e) {
  24. e.printStackTrace();
  25. }
  26. }
  27. }

结果

结果如下,可以看出中间位置的故宫、五环西北角的圆明园、颐和园均没有单车,热力图大体是正确的。

影响精度的因素

  1. 五环位置不准确,爬取完才发现五环右边不是方形,爬取结果少了一块……
  2. 爬取过程中,正在被骑行的车没有统计到
  3. 在爬取过程中,车恰好从未爬取位置移动到了已爬取位置
  4. 报修的车无法统计
  5. 在爬取的过程中,断网了1个小时……时间间隔过长,精确度受影响
  6. 爬取是按照0.003*0.003的经纬度方形区域,没有精确计算是否能够覆盖所有区域,可能四个方形因为爬取单车数量的因素,导致区域之间有死角。

如何获取更精确的数据

  1. 多次爬取,去重,改变步长。这样可以爬取到因为骑行而未被爬取到的车。
  2. 多线程爬取,减少时间。

总结

爬取到北京五环内有345679辆单车(爬取区域少了一点),估计整个北京有45~50万单车吧。

添加新批注
在作者公开此批注前,只有你和作者可见。
回复批注