前言
2022.9学习绘图
利用python进行气象绘图,本文为学习绘制期间记录笔记,分为上、下两部分:处理数据和图像绘制。处理数据流程:读入olr资料,截取夏季/冬季数据,进行10-30dButterWorth带通滤波,计算方差,输出nc资料。
绘制用到资料orl.day.mean.nc
参考文献
[1] Qian Y , Hsu P C , Kazuyoshi K . New real-time indices for the quasi-biweekly oscillation over the Asian summer monsoon region[J]. Climate Dynamics, 2019.
->图1 Variance of 10–30-day BP-filtered OLR (shading; units:W2m-4) and propagation vectors of 10–30-day-filtered OLR for a the annual-mean state, b boreal summer (May– Oct.) and c boreal winter (Nov.– Apr.) during 1980–2012.
代码
滤波原文用的是30年数据,考虑到内存和计算时长本文仅用了1999年1年的数据计算,大概分布也是一致的,如果感兴趣也可以用30年算算,会更趋近于原文。
import time
time_start = time.process_time()
import numpy as np
import xarray as xr
import datetime as dt #处理datetime64格式数据使用,可以解析时间格式
from scipy import signal
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.mpl.ticker as cticker
#######################################################
## 读取olr数据
#######################################################
f = xr.open_dataset('E:/LACS/learn/python/replot/data/olr.day.mean.nc')
# print(f) #xarray.Dataset 对象
twStrt, twLast <