python绘制气象等值线图_利用Python插值绘制等值线图

最近需要根据有限的站位点绘制插值等值线图，在网上中文搜索一通，只发现了这货Matplot Basemap 画湖北地图、插值、等值线，要么就是对这货的转载，这货不提供数据的形式，但是基本的代码思路还是不错的，于是继续转向google英文，搜到了如下的回答，我对代码做了注释，已备别人查询，关于文中提到的数据是txt格式的，我也直接将数据贴在下面了。总结一下：在地图上绘制等值线：

确定基本的绘图框架；

获取采集数据，与地图做映射，并根据映射数据插值；scipy.interpolate.griddata包插值比较快，常用的三种插值方法为liner(基于三角形的线性插补法)，cubic(基于三角形的三次插补法),nearest( 最近邻居插补法)，这些方法定义了匹配数据点的曲面类型，'cubic' 方法生成平滑曲面，而 'linear' 和 'nearest' 分别具有一阶导数和零阶导数不连续。

根据栅格插值数据绘图

# -*- coding: utf-8 -*-

@author: Adwiy Wang

import numpy as np

import pandas as pd

from matplotlib.mlab import griddata

from mpl_toolkits.basemap import Basemap

import matplotlib.pyplot as plt

from matplotlib.colors import Normalize

from scipy.interpolate import gr