麋鹿分布图制作(二)—— 用Python和R在地图上打点

栏目: Python · 发布时间: 5年前

内容简介:上期用Python调用百度API得到地点的经纬度, 今天就将这些点标注在地图上。 首先我们用R语言来标注。 文中用到的文件已上传到画出的地图如下:接下来用python的basemap画一次。

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

library(rgdal)
library(ggplot2)
library(mapproj)
library(xlsx)

# 设置一下工作路径
setwd('D:/class/zonglv')
# 读取文件
milu = read.xlsx('milu.xlsx', 1, encoding = 'UTF-8')
# 读取中国地图
china_map <- readOGR('./mapdata/bou2_4p.shp')
# 转化为数据框
china_map <- fortify(china_map)

p = ggplot() +
     geom_polygon(data=china_map, aes(x=long, y=lat, group = group), colour="grey60", size = 0.25) +
     geom_point(data = milu, aes(x = lon,y = lat, color = time, shape = species), size = 1.5) + 
    coord_map("polyconic")
p

画出的地图如下: 麋鹿分布图制作(二)—— 用Python和R在地图上打点

接下来用python的basemap画一次。

import time
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon
import pandas as pd
import numpy as np

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

ml = pd.read_excel('milu.xlsx')
# 我们还是先把数据读入进来
plt.figure(figsize = (16,8))
m = Basemap(llcrnrlon=77,llcrnrlat=14,urcrnrlon=140,urcrnrlat=51,projection='lcc',lat_1=33,lat_2=45,lon_0=100)
m.drawcoastlines()
m.drawcountries(linewidth=1.5)
m.readshapefile('shapefile/gadm36_CHN_1', 'states', drawbounds=True)
lats = list(ml['lat'])
lons = list(ml['lon'])
x, y = m(lons, lats)
m.scattter(x,y,3,marker = 'o', color = 'red')
plt.show()

麋鹿分布图制作(二)—— 用Python和R在地图上打点

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

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

版权声明: 作者保留权利。文章为作者独立观点,不代表数据人网立场。严禁修改,转载请注明原文链接:http://shujuren.org/article/775.html

数据人网: 数据人学习,交流和分享的平台,诚邀您创造和分享数据知识,共建和共享数据智库。


以上所述就是小编给大家介绍的《麋鹿分布图制作(二)—— 用Python和R在地图上打点》,希望对大家有所帮助,如果大家有任何疑问请给我留言,小编会及时回复大家的。在此也非常感谢大家对 码农网 的支持!

查看所有标签

本站部分资源来源于网络,本站转载出于传递更多信息之目的,版权归原作者或者来源机构所有,如转载稿涉及版权问题,请联系我们

创新公司

创新公司

[美]艾德·卡特姆、埃米·华莱士 / 靳婷婷 / 中信出版社 / 2015-2 / 49.00元

●《玩具总动员》《海底总动员》《机器人瓦力》《飞屋环游记》等14部脍炙人口的动画长片, 近30次奥斯卡奖, 7部奥斯卡最佳动画长片,7次金球奖; ●几乎每一部电影一上映都位居票房榜首,所有电影都曾进入影史票房总榜前50,每一部电影都是商业与艺术的双赢。 ●即便新兴动画公司不断涌现,皮克斯始终保持动画界的王者之位,这一切背后的秘密就在于:不断推动创新的创意管理方式。 你可以从本书......一起来看看 《创新公司》 这本书的介绍吧!

RGB转16进制工具
RGB转16进制工具

RGB HEX 互转工具

SHA 加密
SHA 加密

SHA 加密工具

HEX CMYK 转换工具
HEX CMYK 转换工具

HEX CMYK 互转工具