Python+matplotlib+scipy站点数据绘制气象分布图
Posted thisyan
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Python+matplotlib+scipy站点数据绘制气象分布图相关的知识,希望对你有一定的参考价值。
站点数据绘制气象分布图
数据集任意,本随笔只记录方法,不多请确保数据集中包含以下几个属性(站点编号no且假设有n个,站点经度longitude,站点纬度latitude,站点观测数据data)
1. 取站点经纬度:lon,lat
2. 取站点的观测数据集:data (这个数据维度与站点数量同,即1*n)
3.插值准备:
(1)准备两个列表用于构造网格矩阵:
olon = np.linspace(108,115,97)
olat = np.linspace(24,31,97)
(2)构造网格矩阵
olon,olat = np.meshgrid(olon,olat)
4. 做插值:(需要到入Rbf函数:from scipy.interpolate import Rbf)
func = Rbf(lon,lat,data,function=‘linear‘) rain_data_new = func(olon,olat)
5. 经纬度与地图比例做映射:
xx,yy = m(olon,olat)
6. 上色:
c = m.contourf(xx,yy,rain_data_new)
7. 添加色标:
m.colorbar(c)
另,效果图展示:(尚未做白化处理)
以上是关于Python+matplotlib+scipy站点数据绘制气象分布图的主要内容,如果未能解决你的问题,请参考以下文章
Python中安装numpy matplotlib scipy
在 Python 中输出绘图时迭代线性回归(SciPy 和 MatPlotLib)
Windows10+Python3下安装NumPy+SciPy+Matplotlib
Python地理数据处理环境的搭建和准备(numpy,scipy,matplotlib)