Python+Matplotlib制作超高分辨率动态气象/海洋要素可视化动图

Animated version of Gulf Stream sea surface temperatures done with

@oceanographer from yesterday - inspired by the moving media #artinstallations @refikanadol! Shading and texture comes from simulated zonal currents #SciArt @neuro_rish@matplotlib #FESOM2

FESOM2 SST shaded by U component of velocity

Example of scalar field visualization shaded by corresponding vector field.Data used are from FESOM2 modelsimulations with 1km resolutionin Gulf Stream area. Data areinterpolated to regular 1/100 degree grid.

数据空间上进行了插值，大概就是1km左右的格点分辨率了。

链接: (数据+代码+测试图)

游客，如果您要查看本帖隐藏内容请回复

代码:

[Python] 纯文本查看 复制代码

importxarray asxr
importnumpy asnp
importmatplotlib.pylab asplt
importmatplotlib.cm ascm
###Download data:
#wget [url]https://zenodo.org/record/4759091/files/GF_FESOM2_testdata.tar.gz[/url]
###Unpack data:
#gunzip GF_FESOM2_testdata.tar.gz
#tar -xvf GF_FESOM2_testdata.tar
###Open netCDF files
sst = xr.open_dataset('./GF_FESOM2_testdata/sst.nc')
u = xr.open_dataset('./GF_FESOM2_testdata/u_surf.nc')
v = xr.open_dataset('./GF_FESOM2_testdata/v_surf.nc')
###Plot the figure with U shading
fig, ax = plt.subplots(
1,
1,
constrained_layout=True,
figsize=(10,10),
)
ax.imshow(np.flipud(sst.sst.values), cmap=cm.RdBu_r,  vmin=12, vmax=24)
ax.imshow(np.flipud(u.u_surf.values), alpha=0.3, cmap=cm.gray, vmin=-.3, vmax=0.3)
ax.axis('off');
### plt.savefig('./out_figure.png', dpi=300)
###Plot the figure with V shading
fig, ax = plt.subplots(
1,
1,
constrained_layout=True,
figsize=(10,10),
)
ax.imshow(np.flipud(sst.sst.values), cmap=cm.RdBu_r,  vmin=12, vmax=24)
ax.imshow(np.flipud(v.v_surf.values), alpha=0.3, cmap=cm.gray, vmin=-.3, vmax=0.3)
ax.axis('off');
###Rotate "light source" for U by 45 degrees
direct = np.rad2deg(np.arctan2(v.v_surf.values, u.u_surf.values))
speed_rot = np.hypot(u.u_surf.values, v.v_surf.values)
myangle = 45
U = speed_rot * np.cos(np.deg2rad(myangle - direct))
# V = speed_rot * np.sin(np.deg2rad(myangle - direct))
fig, ax = plt.subplots(
1,
1,
constrained_layout=True,
figsize=(10,10),
)
ax.imshow(np.flipud(sst.sst.values), cmap=cm.RdBu_r,  vmin=12, vmax=24)
ax.imshow(np.flipud(U), alpha=0.3, cmap=cm.gray, vmin=-.3, vmax=0.3)
ax.axis('off');

测试图:

把高分辨率的图拼起来，帧数至少得30 FPS看起来视频就很丝滑了！

                               
登录/注册后可看大图

% r( M& M' j* J2 Y/ |* L                               
登录/注册后可看大图

                               
登录/注册后可看大图