最近在做保护生物学的论文，查论文的过程中发现之前的麋鹿分布图是手绘的，看起来真的是。。。，因此我决定重新画一下这个图。

这个地图的绘制我分为两个部分：（一）调用百度API查找地点的经纬度，（二）用python的basemap模块绘制分布图。

我手录了麋鹿数据（点击下载）。

在使用百度地图之前，我们需要先注册、申请一个百度地图的API，具体操作见教程。

1. import json

2. from urllib.request import urlopen, quote

3. import requests

4. import pandas as pd

5. pd.set_option('display.max_columns', None)

6. # 显示所有的行

8. def getlnglat(address):

9.    url = 'http://api.map.baidu.com/geocoder/v2/'

10.    output = 'json'

11.    ak = '99flG4xrb8OrWR2fOGBfLxqB9DdkhC90'

12.    # 这个是申请的百度地图API的密码

13.    add = quote(address)

14.    #由于本文城市变量为中文，为防止乱码，先用quote进行编码

15.    uri = url + '?' + 'address=' + add  + '&output=' + output + '&ak=' + ak

16.    req = urlopen(uri)

17.    res = req.read().decode()

18.    #将其他编码的字符串解码成unicode

19.    temp = json.loads(res)

20.    #对json数据进行解析

21.    return temp

23. milu = pd.read_excel('miludistrbute.xlsx')

24. milu['lat'] = None

25. milu['lon'] = None

26. # 添加经纬度的两列，在R里面可以不用提前做俩空的列的，直接添加就行了，但是python里还是要加。

28. for i in range(0,len(milu.index)):

30.    if 'result' in getlnglat(milu.iloc[i][3]):

31.    # 这个是防止有些地址查询不到，就会被迫中断循环

32.        milu['lon'][i] = getlnglat(milu.iloc[i][3])['result']['location']['lng']

33.        milu['lat'][i] = getlnglat(milu.iloc[i][3])['result']['location']['lat']

34.        print(milu.iloc[i])

35.    else:

36.        continue

38. milu.to_excel('milu1.xlsx')

39. # 最后保存一下这个数据