用matplotlib画等高线图详解


Posted in Python onDecember 14, 2017

等高线图是在地理课中讲述山峰山谷时绘制的图形,在机器学习中也会被用在绘制梯度下降算法的图形中。

因为等高线的图有三个信息:x,y以及x,y所对应的高度值。

这个高度值的计算我们用一个函数来表述:

计算x,y坐标对应的高度值
def f(x, y):
 return (1-x/2+x**5+y**3) * np.exp(-x**2-y**2)

这个函数看起来挺复杂的,但我们这里只是为了能够获得一个高度值,因此其中函数代表什么意义不用关心,只要知道输入一个x,y,输出一个高度值就可以了。

要画出等高线,核心函数是plt.contourf(),但在这个函数中输入的参数是x,y对应的网格数据以及此网格对应的高度值,因此还需要调用np.meshgrid(x,y)把x,y值转换成网格数据才行,这样完整的代码如下:

画等高线的代码如下:

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

# 计算x,y坐标对应的高度值
def f(x, y):
 return (1-x/2+x**5+y**3) * np.exp(-x**2-y**2)

# 生成x,y的数据
n = 256
x = np.linspace(-3, 3, n)
y = np.linspace(-3, 3, n)

# 把x,y数据生成mesh网格状的数据,因为等高线的显示是在网格的基础上添加上高度值
X, Y = np.meshgrid(x, y)

# 填充等高线
plt.contourf(X, Y, f(X, Y))
# 显示图表
plt.show()

上述代码显示的图形为:

用matplotlib画等高线图详解

这颜色有点太冷了,我们想显示热力图,那只要在plt.contourf()函数中添加属性cmap=plt.cm.hot就能显示热力图,其中cmap代表为color map,我们把color map映射成hot(热力图),此处关键代码为:

# 填充等高线
plt.contourf(X, Y, f(X, Y), cmap=plt.cm.hot)

显示的图为:

用matplotlib画等高线图详解

是否显示得挺热的。:)

上面是用plt.contourf()填充了等高线,但还有一种方式是可以直接显示等高线,而不是填充的方式,例如:

C = plt.contour(X, Y, f(X, Y), 20)

这里20代表的是显示等高线的密集程度,数值越大,画的等高线数就越多。

这样显示的图形为:

用matplotlib画等高线图详解

当然,如果我们不调用前面的plt.contourf()函数,则就会直接显示等高线。

最后我们想在等高线中添加上标注值:

plt.clabel(C, inline=True, fontsize=12)

显示的图为:

用matplotlib画等高线图详解

完整的代码为:

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

# 计算x,y坐标对应的高度值
def f(x, y):
 return (1-x/2+x**5+y**3) * np.exp(-x**2-y**2)

# 生成x,y的数据
n = 256
x = np.linspace(-3, 3, n)
y = np.linspace(-3, 3, n)

# 把x,y数据生成mesh网格状的数据,因为等高线的显示是在网格的基础上添加上高度值
X, Y = np.meshgrid(x, y)

# 填充等高线
plt.contourf(X, Y, f(X, Y), 20, cmap=plt.cm.hot)
# 添加等高线
C = plt.contour(X, Y, f(X, Y), 20)
plt.clabel(C, inline=True, fontsize=12)
# 显示图表
plt.show()

总结

以上就是本文关于用matplotlib画等高线图详解的全部内容,希望对大家有所帮助。感兴趣的朋友可以继续参阅本站:

如有不足之处,欢迎留言指出。感谢朋友们对本站的支持!

Python 相关文章推荐
python的类方法和静态方法
Dec 13 Python
Python中在for循环中嵌套使用if和else语句的技巧
Jun 20 Python
Python开发的实用计算器完整实例
May 10 Python
Python 2.x如何设置命令执行的超时时间实例
Oct 19 Python
python+opencv识别图片中的圆形
Mar 25 Python
Python实现 版本号对比功能的实例代码
Apr 18 Python
Python动态参数/命名空间/函数嵌套/global和nonlocal
May 29 Python
Python字典中的值为列表或字典的构造实例
Dec 16 Python
浅谈tensorflow之内存暴涨问题
Feb 05 Python
python通过文本在一个图中画多条线的实例
Feb 21 Python
翻转数列python实现,求前n项和,并能输出整个数列的案例
May 03 Python
python实现扫雷游戏的示例
Oct 20 Python
Python数据结构与算法之图的广度优先与深度优先搜索算法示例
Dec 14 #Python
python实现发送邮件功能代码
Dec 14 #Python
python正则实现计算器功能
Dec 14 #Python
Python数据结构与算法之使用队列解决小猫钓鱼问题
Dec 14 #Python
python实现BackPropagation算法
Dec 14 #Python
python实现随机梯度下降(SGD)
Mar 24 #Python
Python实现将一个正整数分解质因数的方法分析
Dec 14 #Python
You might like
PHP统计目录下的文件总数及代码行数(去除注释及空行)
2011/01/17 PHP
php程序总是提示验证码输入有误解决方案
2015/01/07 PHP
小试JQuery的AutoComplete插件
2011/05/04 Javascript
Javascript中的delete操作符详细介绍
2014/06/06 Javascript
jquery JSON的解析方式示例介绍
2014/07/27 Javascript
jQuery实现的多级下拉菜单效果代码
2015/08/24 Javascript
canvas绘制环形进度条
2017/02/23 Javascript
微信小程序 密码输入(源码下载)
2017/06/27 Javascript
JS实现的计数排序与基数排序算法示例
2017/12/04 Javascript
Angular实现的敏感文字自动过滤与提示功能示例
2017/12/29 Javascript
如何使用electron-builder及electron-updater给项目配置自动更新
2018/12/24 Javascript
如何优雅地在vue中添加权限控制示例详解
2019/03/07 Javascript
JS script脚本中async和defer区别详解
2020/06/24 Javascript
JS实现京东商品分类侧边栏
2020/12/11 Javascript
JavaScript canvas实现雨滴特效
2021/01/10 Javascript
Python解释执行原理分析
2014/08/22 Python
python实现识别相似图片小结
2016/02/22 Python
Python二叉搜索树与双向链表转换实现方法
2016/04/29 Python
利用Python查看目录中的文件示例详解
2017/08/28 Python
python写入已存在的excel数据实例
2018/05/03 Python
浅谈Python中的全局锁(GIL)问题
2019/01/11 Python
Python tempfile模块生成临时文件和临时目录
2020/09/30 Python
Pytorch - TORCH.NN.INIT 参数初始化的操作
2021/02/27 Python
CSS3教程:background-clip和background-origin
2008/10/17 HTML / CSS
做一个能自适应高度的textarea的示例代码
2019/09/06 HTML / CSS
雅诗兰黛美国官网:Estee Lauder美国
2016/07/21 全球购物
特步官方商城:Xtep
2017/03/21 全球购物
研究生自我鉴定范文
2013/10/30 职场文书
奶茶专卖店创业计划书
2014/01/18 职场文书
论文评语大全
2014/04/29 职场文书
旅游节目策划方案
2014/05/26 职场文书
生物科学专业毕业生求职信
2014/06/02 职场文书
四风专项整治工作情况汇报
2014/10/28 职场文书
市级三好学生评语
2014/12/29 职场文书
浅析MySQL如何实现事务隔离
2021/06/26 MySQL
MySQL批量更新不同表中的数据
2022/05/11 MySQL