利用R语言从高德地图生成AOI(Area Of Interests)

高德地图开放了POI的WEB API接口（也就是搜索那个API接口）。然而很多时候我们并不满足于将地物抽象成一个点。在城市这个尺度分析的时候，大型的住宅小区抽象成点往往就难以满足需求了，更何况还有更大面积的景区、学校之类。这个时候，兴趣区（Area Of Interests）就是我们需要的东西了。遗憾的是，高德并没有开放这个API接口，百度貌似也没有。但还是可以通过高德的地物详情查询接口获取兴趣区的边界坐标信息，很多人已经摸索出了解决方案，基本都是基于Python的，我这里狗尾续貂，用R语言来解决一下。最后生成的类型是sf。

---

p_load("httr",'tidyverse','jsonlite','sf','rgdal','rlist')
library('Rgctc2')

需要用的到的包环境，Rgctc2是我为了在高德坐标系和wgs84坐标系之间转换搞的一个小包，可以在这里下载。

获取AOI ID

前人似乎都是通过查看页面源代码的方式来获取ID的，似乎高德提供了一个很简便的方法。