画图报错问题

lynpatty

想用计算出来的温度指数画北美地区的图，用这个代码画前几年都没问题，单画到这个年份就报错“bboxes can not be empty”,但是确认过这个区域的数据不会是0，所以想问问这个是什么问题呢？报错：2.6718372428853816Traceback (most recent call last):

  File "/Users/untitled2.py", line 134, in <module>
    cf = ax.contourf(lon,lat,md_values_plot,extend='both')

  File "/Users/opt/anaconda3/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py", line 1404, in contourf
    for col in result.collections

  File "/Users/opt/anaconda3/lib/python3.7/site-packages/matplotlib/transforms.py", line 701, in union
    raise ValueError("'bboxes' cannot be empty")

ValueError: 'bboxes' cannot be empty

1. year_se = 1991
   
2. mon_se = 7
   
3. md_values_plot=np.nansum(md_values[(year==year_se)&(mon==mon_se),:,:],axis = 0)
   
4. md_values_plot = np.where(md_values_plot==0,np.nan,md_values_plot)
   
5. md_index=np.nanmax(md_values_plot)
   
6. print(md_index)
   
7. 
   
8. proj = ccrs.PlateCarree()
   
9. fig = plt.figure(figsize=(4, 4), dpi=550)  # 创建画布
   
10. extent=[-130,-70,25,50]
    
11. ax = fig.subplots(1, 1, subplot_kw={'projection': proj})  # 创建子图
    
12. ax.set_extent(extent)
    
13. ax.add_feature(cfeature.LAND)####添加陆地######
    
14. ax.add_feature(cfeature.COASTLINE,lw=0.3)#####添加海岸线#########
    
15. ax.add_feature(cfeature.RIVERS,lw=0.25)#####添加河流######
    
16. ax.add_feature(cfeature.LAKES)######添加湖泊#####
    
17. ax.add_feature(cfeature.BORDERS, linestyle='-',lw=0.25)####不推荐，我国丢失了藏南、台湾等领土############
    
18. ax.add_feature(cfeature.OCEAN)######添加海洋########
    
19. ax.add_feature(cfeature.OCEAN.with_scale('10m'))
    
20. ax.add_feature(cfeature.LAND.with_scale('10m'))
    
21. ax.add_feature(cfeature.RIVERS.with_scale('10m'),lw=0.6)
    
22. ax.add_feature(cfeature.LAKES.with_scale('10m'))
    
23. ax.add_feature(cfeature.BORDERS.with_scale('50m'), linestyle='-',lw=0.6)
    
24. ax.add_feature(cfeature.COASTLINE.with_scale('10m'),lw=0.5)
    
25. gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, linewidth=0.2, color='k', alpha=0.5, linestyle='--')
    
26. gl.xlabels_top = False ##关闭上侧坐标显示
    
27. gl.ylabels_right = False ##关闭右侧坐标显示
    
28. gl.xlabel_style={'size':5}
    
29. gl.ylabel_style={'size':5}
    
30. gl.xformatter = LONGITUDE_FORMATTER ##坐标刻度转换为经纬度样式
    
31. gl.yformatter = LATITUDE_FORMATTER
    
32. gl.xlocator = mticker.FixedLocator(np.arange(extent[0], extent[1]+0.1, 5))
    
33. gl.ylocator = mticker.FixedLocator(np.arange(extent[2], extent[3]+0.1, 5))
    
34. #plot
    
35. 
    
36. 
    
37. Colors = ('#FFFFFF','#C0C0C0','#696969','#FFD700','#FF6347','#FF0000','#A0522D','#8B0000')
    
38. cf = ax.contourf(lon,lat,md_values_plot,extend='both')
    
39. cbar=plt.colorbar(cf,ax=ax,orientation='horizontal',pad=0.15,shrink=1)
    
40. cbar.ax.tick_params(labelsize=8)
    
41. cbar.set_label('HWMId',fontdict={'size':12})
    
42. 
    
43. cbar.ax.set_title('West U.S 1991 July')
    
44. plt.show()

复制代码

解脱

绘图的时候没有色阶参数。