Matplotlib Toolkits:地图绘制工具
Matplotlib Toolkits:地图绘制工具
有没有一种可以直接在详细地图(如谷歌地图)上绘制上百万坐标点的工具???谷歌地图坐标点多了也不能绘制了。
Basemap
(Not distributed with matplotlib,要手动安装)
Plots data on map projections, with continental and political boundaries, see basemap docs.
安装
先下载anacondaPython 3版本
conda install -c scitools basemap
[Installing¶]
使用
分辨率设置
if resolution == 'c': area_thresh = 10000. elif resolution == 'l': area_thresh = 1000. elif resolution == 'i': area_thresh = 100. elif resolution == 'h': area_thresh = 10. elif resolution == 'f': area_thresh = 1.
[basemap]
[地理信息可视化——Python之matplotlib Basemap简介]
[So You’d Like To Make a Map Using Python]
静态地图示例
# coding=utf8 import os import sys import matplotlib.pyplot as plt import pandas as pd CWD ] os.makedirs(os.path.join(CWD, 'middlewares'), exist_ok=True) sys.path.append(os.path.join(CWD, '../../..')) df = pd.read_pickle(os.path.join(CWD, 'middlewares/df.pkl')) # ca datasets ll = df[['longitude', 'latitude']].values print(len(ll)) x, y ], ll[] from mpl_toolkits.basemap import Basemap # 创建一个地图用于绘制。我们使用的是墨卡托投影,并显示整个世界。 m , , , , , resolution='c') # 绘制海岸线,以及地图的边缘 m.drawcoastlines() m.drawmapboundary() m.drawcountries() m.drawstates() m.drawcounties() x, y = m(y, x) m.scatter(x, y, , marker='.', color='r') plt.show()
Debug
m.drawcounties() UnicodeDecodeError: 'utf-8' codec can't decode byte 0xf1 in position 2: invalid continuation byte
解决:shapefile.py源码有错:修改编码为latin-1就可以了,也是醉了。。。
return v.decode('latin-1') # return v.decode('utf-8')
[[译]比一比:Python的七个数据可视化工具 静态地图]
Cartopy
评价:安装麻烦,map功能也不能使用了,估计开发者也抛弃了。。。真是跪了。。。lz已弃用
(Not distributed with matplotlib)
Cartopy is a Python package designed to make drawing maps for data analysis and visualisation as easy as possible.
Cartopy makes use of the powerful PROJ.4, numpy and shapely libraries and has a simple and intuitivedrawing interface to matplotlib for creating publication quality maps.
Some of the key features of cartopy are:
- object oriented projection definitions
- point, line, vector, polygon and image transformations between projections
- integration to expose advanced mapping in matplotlib with a simple and intuitive interface
- powerful vector data handling by integrating shapefile reading with Shapely capabilities
An alternative mapping library written for matplotlib v1.2 and beyond.Cartopy builds on top ofmatplotlib to provide object oriented map projection definitions and closeintegration with Shapely for powerful yet easy-to-use vector data processingtools.
Note: shapely,是geos的python封装,而geos是jts的c++移植版本。[python地理数据处理库geopy]
安装
先下载anacondaPython 3版本
conda install -c scitools cartopy
[Installing Cartopy¶]
下载配角:地图数据
下载地址:http://www.naturalearthdata.com/downloads/
里面有三种分辨率的shape地图数据可选,方便起见,分别下载三种分辨率中的physical数据中的Coastline和Land数据,每种数据下载后都是一个压缩包,如下载了1:10分辨率的physical中的coastline数据压缩包:ne_10m_coastline.zip,解压缩后有6个文件,其中“ne_10m_coastline.README”和“ne_10m_coastline.VERSION”可直接删除,剩余4个,进行改名操作,扩展名前面的文件名,如“ne_10m_coastline”,修改为“10m_coastline”,即去掉“ne_”,4个文件分别这样更改。再下载1:50和1:110的文件分别进行此操作。所有地图文件下载、解压、更名完毕后,拷贝到一个文件夹下。我的文件夹列表如下图,把这些文件全选(注意进入文件夹目录,全选文件,不带文件夹),复制粘贴到D:\Program Files\WinPython-32bit-2.7.9.3\settings\.local\share\cartopy\shapefiles\natural_earth\physical 目录下(该目录根据自己所装的python而定,运行(1)中的程序后,程序会自动创建physical文件夹,具体该文件夹在哪,搜索电脑文件找找看),我安装的是winpython2.7.9.3,physical目录就位于上面这个目录中,所以我把所有shape地图文件拷贝到了该physical目录下。
或者下载详细地图:http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_admin_1_states_provinces.zip解压到:c:\Users\《用户》.local\share\cartopy\shapefiles\natural_earth\cultural\
使用
最近不能使用map服务了:
As of July 11, 2016, direct tile access has been discontinued.Please visit our blog post for more information:http://goo.gl/xB0xXt
[Cartopy]
[examples: Cartopy gallery]
谷歌地图python接口gmplot
10w以上数据坐标点绘制js就基本崩了不能用了。
Plotting data on Google Maps, the easy way. A matplotlib-likeinterface to generate the HTML and javascript to render all thedata you'd like on top of Google Maps. Several plotting methodsmake creating exploratory map views effortless.
# coding=utf8 import os import sys import pandas as pd CWD ] df = pd.read_pickle(os.path.join(CWD, 'middlewares/df.pkl')) # ca datasets ll = df[['longitude', 'latitude']].values print(len(ll)) x, y ], ll[] import gmplot gmap ) # gmap.plot(latitudes, longitudes, 'cornflowerblue', edge_width=10) gmap.scatter(x, y, ) # gmap.heatmap(heat_lats, heat_lngs) gmap.draw("map.html")
[gmplot]
高德地图amap python接口
1)读取文件信息,这里用pandas库以dataframe形式.
2)高德地图添加点有50数量限制,所以lz就直接弃用了。。。
3)将查询经纬度后的结果存为csv
百度地图python接口
其它工具
[使用地图工具,给出经纬度,画出受影响的中国省份,(比如台风经过的path、各省经济问题热点图。。。)]
from: Matplotlib Toolkits:地图绘制工具
ref: [Mapping Toolkits]
Matplotlib Toolkits:地图绘制工具的更多相关文章
- arcgis api for javascript 学习(五) 实现地图绘制工具
1.本文实现的功能为鼠标与地图之间的交互,能够在地图上绘制不同形状的图形 2.代码部分主要讲到的为Graphic函数的相关功能 <!DOCTYPE html> <html> & ...
- ArcGisJS实现地图常用工具条、距离测量和面积测量(非官方实例)
常用地图工具包括:平移.拉框缩小.拉框放大.全图.距离测量.面积测量.清除标记,距离测量.面积测量没有使用官方自带的组件代码. 1.距离测量 2.面积测量 3.源代码 <!DOCTYPE htm ...
- 矢量图绘制工具Svg-edit调整画布的大小
矢量图绘制工具Svg-edit调整画布的大小 ------------------------------ ------------------------
- C#小程序飞行棋地图绘制
1. 初始化地图,在绘制时可先将地图进行初始化,用数组来存储关卡的位置,然后利用循环给地图中 关卡所在处赋予代表关卡的值. 关键代码如下 /// <summary> /// 初始化游戏地图 ...
- opengl典型例程立方体投影与地图绘制
立方体投影 http://www.cnblogs.com/opengl/p/3790450.html 地图绘制 http://www.cnblogs.com/opengl/p/3790354.html
- 让IE8支持HTML5及canvas功能!chart.js图表绘制工具库IE8上兼容方案
第一步,我们加上对html5的支持. <!--[if IE]> <script src="/public/html5.js" type="text/ja ...
- Graphviz--图形绘制工具
可以试试ProcessOn.com, 一个专业在线画流程图的工具,使用Chrome&Firefox浏览器,支持快捷键,非常方便.快捷键如下:Ctrl+A 全选,当移动整张图时非常方便Ctrl+ ...
- matplotlib之scatter绘制散点
# 使用matplotlib.pyplot.scatter绘制散点 import matplotlib.pyplot as plt from pylab import mpl # 设置默认字体,解决中 ...
- 「工具」三分钟了解一款在线流程绘制工具:Whimsical
Whimsical 是一款在线流程绘制工具,只需要一个浏览器就随时随地绘制精美的流程图.除了流程图(Flowcharts)功能,官方还推出了线框图(Wireframes).便利贴(Sticky Not ...
随机推荐
- 学习react系列(八)—— mixins迁移
先来介绍一下mixins(混入) 先来看一段代码: const mixin = function(obj, mixins) { const newObj = obj; newObj.prototype ...
- 技巧:如何提升Oracle用户密码的安全性
环境:Oracle 11.2.0.4 客户需求:主要背景是数据库中有很多业务用户名,且由于部分用户缺乏安全意识,甚至直接将自己的密码设置为和用户名一样,目前客户期望密码设置不要过于简单,最起码别和用户 ...
- window10下安装linux虚拟机
一.准备工具 虚拟机.镜像.putty 1.安装虚拟机 VMware Workstation Pro 安装成功之后需要输入密钥,请点击以下链接 http://www.360doc.com/conten ...
- [LeetCode] Bold Words in String 字符串中的加粗单词
Given a set of keywords words and a string S, make all appearances of all keywords in S bold. Any le ...
- 实验吧_Guess Next Session&Once More(代码审计)
Guess Next Session 看题目提示,是一道代码审计: <?php session_start(); if (isset ($_GET['password'])) { if ($_G ...
- [Vijos 2024]无向图最短路径
Description 无向图最短路径问题,是图论中最经典也是最基础的问题之一.本题我们考虑一个有 $n$ 个结点的无向图 $G$.$G$ 是简单完全图,也就是说 $G$ 中没有自环,也没有重边,但任 ...
- ●Joyoi Normal
题链: http://www.joyoi.cn/problem/tyvj-1953题解: 定义d(u,v)这个函数,满足: d(u,v)=1,当且仅当在点分树中,u是v的祖先 d(u,v)=0,其它情 ...
- ●洛谷P3687 [ZJOI2017]仙人掌
题链: https://www.luogu.org/problemnew/show/P3687题解: 计数DP,树形DP. (首先对于这个图来说,如果初始就不是仙人掌,那么就直接输出0) 然后由于本来 ...
- ●CodeForces 280D k-Maximum Subsequence Sum
题链: http://codeforces.com/problemset/problem/280/D 题解: 神题,巨恶心.(把原来的那个dp题升级为:序列带修 + 多次询问区间[l,r]内取不超过k ...
- hdu 5012(bfs)
题意:给你2个 骰子,让你通过翻转使第一个变成第二个,求最少翻转数 思路:bfs #include<cstdio> #include<iostream> #include< ...