python 如何用xarray和cartopy使流线连续？

fdbelqdn 于 2023-03-16 发布在 Python

关注(0)|答案(1)|浏览(140)

我现在使用xarray数据集在全球Map上用cartopy绘制流线，但是流线在穿越dateline的时候不连续，有没有人知道怎么解决这个问题，代码（只有streamplot部分）和输出图如下：

wnd = xr.merge([u[::2,::2],v[::2,::2]])
proj = ccrs.PlateCarree(central_longitude = 180)
fig, ax = plt.subplots(1,1,
                     subplot_kw={'projection': proj},
                     figsize=(12,7))
wind_plt = wnd.plot.streamplot(ax=ax,
                               transform=ccrs.PlateCarree(),
                               x='longitude', y='latitude',
                               u='u', v='v',
                               arrowsize=2,arrowstyle='->',
                               color="black",
                               )

ax.coastlines(color = 'dimgray')
ax.set_extent([0,360,-90,90],crs=ccrs.PlateCarree())
ax.set_xticks(np.arange(0, 361,90), crs=ccrs.PlateCarree())
ax.set_yticks(np.arange(-90,91,45), crs=ccrs.PlateCarree())  
ax.set_xlabel(' ')
ax.set_ylabel(' ')
ax.set_title(' ')
ax.text(70, 30, '03-12', fontsize = 14, transform=ccrs.PlateCarree())

ax.get_xaxis().set_visible(False)
ax.get_yaxis().set_visible(False)
plt.axis('off')

plt.savefig('0312.png', bbox_inches = 'tight', pad_inches = 0)
plt.close()

谢谢你，欢迎所有的意见！

python

来源：https://stackoverflow.com/questions/75729486/how-to-make-streamline-with-xarray-and-cartopy-continuous

1条答案

按热度按时间

ux6nzvsh1#

我怀疑您的流线（流线图数据：X，Y，U，V）未创建为填充绘图区所需的完整范围。或者，未正确创建。
这是我的演示代码，它显示了来自一组良好数据的完整流图，供您尝试。
请注意，我的数据流图是完整的，因此，将没有差距周围的日期变更线。

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import numpy as np

noproj = projection=ccrs.PlateCarree()
myproj = projection=ccrs.PlateCarree(central_longitude=180)

# Creating dataset
w = 180
h = 90
Y, X = np.mgrid[-h:h:100j, -w:w:100j]
U = -1 - X**2 + Y
V = 1 + X - Y**2
#speed = np.sqrt(U**2 + V**2)

fig = plt.figure(figsize=[5,6.7])

# Plot on ordinary projection
ax1 = fig.add_subplot(211, projection = noproj)
ax1.streamplot(X, Y, U, V, density=[0.4, 0.8], transform=noproj)
ax1.set_extent([-180, 180, -90, 90], crs=noproj)
ax1.coastlines(resolution='110m', alpha=0.3)
ax1.gridlines(draw_labels=True)
ax1.set_title("Zero-meridian-centered")

# Plot on projection with shifted center
ax2 = fig.add_subplot(212, projection = myproj)
ax2.streamplot(X, Y, U, V, density=[0.4, 0.8], transform=noproj)
ax2.set_extent([-180, 180, -90, 90], crs=noproj)
ax2.coastlines(resolution='110m', alpha=0.3)
ax2.gridlines(draw_labels=True)
ax2.set_title("Dateline-centered")
plt.show()

赞(0）回复(0）举报 2023-03-16

我来回答

python 如何用xarray和cartopy使流线连续？

1条答案

相关问题

热门标签

最新问答