matplotlib-base_matplotlib3d柱状图

matplotlib-base_matplotlib3d柱状图作者:Vamei出处:http://www.cnblogs.com/vamei欢迎转载,也请保留这段声明。谢谢!在数据可视化过程中,我们常常需要将数据根据其采集的地理位置在地图上显示出来。比如说我

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

作者:Vamei 出处:http://www.cnblogs.com/vamei 欢迎转载,也请保留这段声明。谢谢!

 

在数据可视化过程中,我们常常需要将数据根据其采集的地理位置在地图上显示出来。比如说我们会想要在地图上画出城市,飞机的航线,乃至于军事基地等等。通常来说,一个地理信息系统都会带有这样的功能。今天我们讨论如何在Python上实现,并且使用免费的工具包。

 

matplotlib是Python常用的数据绘制包。它基于numpy的数组运算功能。matplotlib绘图功能强大,可以轻易的画出各种统计图形,比如散点图,条行图,饼图等。matplotlib常与numpy和scipy相配合,用于许多研究领域。他们是免费工具,但其功能足可以与科研界的大佬Matlab竞争。

Basemap是Matplotlib的一个子包,负责地图绘制。在数据可视化过程中,我们常需要将数据在地图上画出来。比如说我们在地图上画出城市人口,飞机航线,军事基地,矿藏分布等等。这样的地理绘图有助于读者理解空间相关的信息。

 

我们下面用Basemap画出亚洲主要城市的人口。如下图,人口的数量用圆圈的大小表示:

matplotlib-base_matplotlib3d柱状图

数据如下(我从Wikipedia上整理的,你可以随意使用)。将数据保存在文件major_city:

Shanghai 23019148  31.23N  121.47E  China
Mumbai   12478447  18.96N  72.82E   India
Karachi  13050000  24.86N  67.01E   Pakistan
Delhi    16314838  28.67N  77.21E   India
Manila   11855975  14.62N  120.97E  Philippines
Seoul    23616000  37.56N  126.99E  Korea(South)
Jakarta  28019545   6.18S  106.83E  Indonesia
Tokyo    35682460  35.67N  139.77E  Japan
Peking   19612368  39.91N  116.39E  China

第一列是城市名,第二列是人口,第三第四列为纬度和经度,最后一列为所在国家。

 

下面是我的Python代码,用以绘制上面的地图:

matplotlib-base_matplotlib3d柱状图matplotlib-base_matplotlib3d柱状图View Code

# Written by Vamei, http://www.cnblogs.com/vamei/ # Feel free to use or modify this script.

from mpl_toolkits.basemap import Basemap import matplotlib.pyplot as plt import numpy as np #============================================# read data
names = [] pops = [] lats = [] lons = [] countries = [] for line in file("../data/major_city"): info = line.split() names.append(info[0]) pops.append(float(info[1])) lat = float(info[2][:-1]) if info[2][-1] == 'S': lat = -lat lats.append(lat) lon = float(info[3][:-1]) if info[3][-1] == 'W': lon = -lon + 360.0 lons.append(lon) country = info[4] countries.append(country) #============================================ # set up map projection with # use low resolution coastlines.
map = Basemap(projection='ortho',lat_0=35,lon_0=120,resolution='l') # draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=0.25) map.drawcountries(linewidth=0.25) # draw the edge of the map projection region (the projection limb)
map.drawmapboundary(fill_color='#689CD2') # draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0,360,30)) map.drawparallels(np.arange(-90,90,30)) # Fill continent wit a different color
map.fillcontinents(color='#BF9E30',lake_color='#689CD2',zorder=0) # compute native map projection coordinates of lat/lon grid.
x, y = map(lons, lats) max_pop = max(pops) # Plot each city in a loop. # Set some parameters
size_factor = 80.0 y_offset = 15.0 rotation = 30
for i,j,k,name in zip(x,y,pops,names): size = size_factor*k/max_pop cs = map.scatter(i,j,s=size,marker='o',color='#FF5600') plt.text(i,j+y_offset,name,rotation=rotation,fontsize=10) plt.title('Major Cities in Asia & Population') plt.show()

程序分为两个部分,第一部分为从文件读取数据并处理。第二部分才是真正用basemap绘图。

 

地图的大小、投影方法等重要信息,是在Basemap()的调用中实现的:

map = Basemap(projection=ortho,lat_0=35,lon_0=120,resolution=l)

projection参数规定了投影方法。改变投影方法,绘图结果也将非常不同。

 

城市所在位置是经纬度。我们想要把经纬度对应图像的像素点,需要转换:

x, y = map(lons, lats)

这个语句转换为图像上的位置。

 

最后,调用绘制散点图的方法scatter():

cs = map.scatter(i,j,s=size,marker=o,color=#FF5600)

在地图上画出数据。

 

总结

matplotlib中的Basemap是很好用的,具有专业标准的地图绘制工具。它可以与matplotlib的一般绘图功能结合,在地图上绘制数据。

 

文中需要的软件包:

numpy, matplotlib, mpl_toolkits

可以使用pip来下载安装各个包。

 

在Ubuntu的repository中,你可以找到python-matplotlib包。

 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/167708.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • intel处理器_bios cppc

    intel处理器_bios cppcC++string类的find()函数怎么用

    2022年10月23日
  • busybox 安装mysql_busybox 安装问题解决「建议收藏」

    busybox 安装mysql_busybox 安装问题解决「建议收藏」直接编译错误1。loginutils/passwd.c:93:16:error:storagesizeof‘rlimit_fsize’isn’tknown解决方法:在busybox根目录下查找到文件:find-namelibbb.h在libbb.h中添加#include1。MakefileMakefile:431:***mixedimplicitandnormalr…

  • Java集合汇总篇「建议收藏」

    Java集合汇总篇「建议收藏」一.集合框架Java集合框架一些列的接口和类来实现很多常见的数据结构和算法,例如LinkedList就是集合框架提供的实现了双向链表的数据结构,关于这一篇文章建议大家收藏,我会不断地完善和扩充它的内容,例如最下面的系列文章我以后也会对它进行不断的更新集合框架的接口集合框架提供了很多接口,这些接口都包含了特定的方法来实现对集合上的特定操作)我们将要学习这些接口以及子接口和它们的各种实现类,在开始之前我们先简单学习一下这些广泛运用的接口,可以看到整个集合框架,总共有三个顶级接口Collecti

  • python中for循环语句例子_python怎么循环1到8不要4

    python中for循环语句例子_python怎么循环1到8不要4这篇文章主要介绍了python中关于for循环使用过程中的碎碎念,需要的朋友可以参考下为什么要挑战自己在代码里不写forloop?因为这样可以迫使你去使用比较高级、地道的语法或库。文中以python为例子,讲了不少大家其实在别人的代码里都见过、但自己很少用的语法。这是一个挑战。我要你避免在任何情况下写for循环。同样的,我也要你找到一种场景——除了用for循环以外,用其他方法写都太难。请分享你的…

  • matlab中axis的使用

    matlab中axis的使用%提示disp('该功能练习axis功能');%初始化快捷式数组x=-2*pi:pi/20:2*pi;y=sin(x);plot(x,y);title('sin(x)图形

  • initramfs 命令_linux initramfs

    initramfs 命令_linux initramfs原文地址:http://www.cnblogs.com/pied/archive/2013/01/29/2880718.html这个是翻译来的,原文地址:http://www.landley.net/writing/rootfs-howto.html怎样使用initramfs 工作过程简述在2.6kernel启动时,它把rootfs作为它的第一个文件系统挂载(注意:这里的rootfs是真名!!!不…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号