上期用Python调用百度API得到地点的经纬度， 今天就将这些点标注在地图上。
首先我们用R语言来标注。
文中用到的文件已上传到Github。

1. library(rgdal)

2. library(ggplot2)

3. library(mapproj)

4. library(xlsx)

6. # 设置一下工作路径

7. setwd('D:/class/zonglv')

8. # 读取文件

9. milu = read.xlsx('milu.xlsx', 1, encoding = 'UTF-8')

10. # 读取中国地图

11. china_map <- readOGR('./mapdata/bou2_4p.shp')

12. # 转化为数据框

13. china_map <- fortify(china_map)

15. p = ggplot() +

16.     geom_polygon(data=china_map, aes(x=long, y=lat, group = group), colour="grey60", size = 0.25) +

17.     geom_point(data = milu, aes(x = lon,y = lat, color = time, shape = species), size = 1.5) +

18.    coord_map("polyconic")

19. p

画出的地图如下：

接下来用python的basemap画一次。

1. import time

2. from mpl_toolkits.basemap import Basemap

3. import matplotlib.pyplot as plt

4. from matplotlib.patches import Polygon

5. import pandas as pd

6. import numpy as np

在这些调用的模块中，basemap的安装可能会难一点，这里请参考Basemap手册，另外，Basemap tutorial也不错，参考这两个链接应该就可以画很多图了，但是我在跑Basemap的时候很多关键的代码都会出现Can’t set attribute的报错，目前还没有找到原因，勉强跑了一个图出来，如果有大神看到这个帖，麻烦告知一下解决方法。

1. ml = pd.read_excel('milu.xlsx')

2. # 我们还是先把数据读入进来

3. plt.figure(figsize = (16,8))

4. m = Basemap(llcrnrlon=77,llcrnrlat=14,urcrnrlon=140,urcrnrlat=51,projection='lcc',lat_1=33,lat_2=45,lon_0=100)

5. m.drawcoastlines()

6. m.drawcountries(linewidth=1.5)

7. m.readshapefile('shapefile/gadm36_CHN_1', 'states', drawbounds=True)

8. lats = list(ml['lat'])

9. lons = list(ml['lon'])

10. x, y = m(lons, lats)

11. m.scattter(x,y,3,marker = 'o', color = 'red')

12. plt.show()

在这里需要注意的是，这个代码无论是在Jupyter notebook还是在sublime text 3中运行都是会报错的，但是在jupyter notebook中是可以跑出地图来的，如上图。

在两个程序之中跑的代码，我们可以看到，有些点的位置是不准确的，我不知道这是由于百度地图API给的经纬度不够准确还是我下载的shp文件和百度地图不太兼容，如果我以后找到原因会更新本文。