几个要点
1、原始数据是站点数据,故要用插值法转成格点数据
2、涉及到指定地图范围的色斑图,必须用到地图白化,借鉴原文:哔哩哔哩 搜索摸鱼气象& Pthon 白化绘图基础,这里就不粘贴了,涉嫌广告。
3、绘制等值线用函数Axes.contour(*args, data=None, **kwargs)
等值线填图用函数Axes.contourf(*args, data=None, **kwargs)
先贴一张成品图:强降水次数分布图
代码涉及到的数据有一个降水的excel文档,和一个本市的分县shp文档
# -*- coding: utf-8 -*-
import numpy as np
import pandas as pd
from scipy.interpolate import Rbf # 径向基函数 : 将站点信息插到格点上 用于绘制等值线
import matplotlib.pyplot as plt
import matplotlib.colors as colors
import matplotlib as mpl
import cartopy.crs as ccrs # 默认投影
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
import geopandas as gpd
from matplotlib.path import Path
fr