【已解决】python极射赤面投影出现白条

Iree_Aurora

想画一下某一天100hPa纬向风场（等值线）和滤波后风场（填色）的图，但是填色图有空白的这部分，请问各位大佬，下面是哪里设置错了啊……

file = nc.Dataset('uwnd.1979.nc')
uwnd = np.array(file.variables['uwnd'][:])
lat = file.variables['lat'][:]
lon = file.variables['lon'][:]
Uwnd = uwnd[:, 11, :, :]

#======绘图=======
fig = plt.figure(figsize=(12, 8))
leftlon, rightlon, lowerlat, upperlat = (0, 360, 0, 90)

proj = ccrs.Orthographic(central_longitude=-90, central_latitude=90)
data_proj = ccrs.PlateCarree(central_longitude=0)
ax = fig.add_subplot(projection=proj)
ax.set_global()

ax.add_feature(cfeature.LAKES, alpha=0.5) #添加湖泊
ax.coastlines(linewidth=0.5) #添加海岸线
# ax.gridlines(draw_labels=True, x_inline=False, y_inline=False)

lon_formatter = LongitudeFormatter(zero_direction_label=True)
lat_formatter = LatitudeFormatter()
ax.xaxis.set_major_formatter(lon_formatter)
ax.yaxis.set_major_formatter(lat_formatter)
ax.gridlines(color='gray', linestyle='--', xlocs=np.arange(-180, 200, 20))
ax.gridlines(draw_labels=True)
cf= ax.contourf(lon, lat, Uwnd[160, :, :], transform=data_proj)
plt.clabel(cl, fmt='%i')
"""添加色标"""
plt.colorbar(cf, pad=0.042, orientation='horizontal', shrink=0.87, aspect=35)
plt.tick_params(labelsize=13)  
plt.show()

1099221723

这是数据的空缺，例如数据经度是0-359，那么359-360是不会画的，因为没有这个区间的数据，解决也很简单，自己再359后面补一个360，数据都是0 的，就好了

Iree_Aurora

1099221723 发表于 2021-11-30 12:26
这是数据的空缺，例如数据经度是0-359，那么359-360是不会画的，因为没有这个区间的数据，解决也很简单，自 ...

谢谢您的回答！但是我刚接触python不久，还是不太理解您的意思，可以具体写一下吗，麻烦您了

edwardli

7. 零度0°经线数据的处理
from cartopy.util import add_cyclic_point
lon=np.arange(0.125,360.125,0.25)
data, lon = add_cyclic_point(data, coord=lon) #插值0°经线也有数据
lon=np.arange(0,360.25,0.25)

Iree_Aurora

edwardli 发表于 2021-11-30 15:22
7. 零度0°经线数据的处理
from cartopy.util import add_cyclic_point
lon=np.arange(0.125,360.125,0.2 ...

谢谢！学习了学习了！！

易烊启迪

Iree_Aurora 发表于 2021-11-30 18:17
谢谢！学习了学习了！！

楼主您好，如果只想画出经度为0-360，纬度是40-90的地图应该怎么修改代码。

来自苹果APP客户端

yangzi88552

易烊启迪 发表于 2022-8-22 14:40
楼主您好，如果只想画出经度为0-360，纬度是40-90的地图应该怎么修改代码。

楼主程序中的leftlon, rightlon, lowerlat, upperlat = (0, 360, 0, 90) ，改为你想要的纬度就可以了，比如你的改为leftlon, rightlon, lowerlat, upperlat = (0, 360, 40, 90)