画图时阴影图被地图遮盖

schliezer

地图遮住了阴影，但我明明是先画的地图，另外，有没有办法可以只保留第一列子图的坐标？所有坐标挤在一起很难看。代码和图：

1. 2# -*- coding: utf-8 -*-
   
2. """
   
3. Created on Thu Jan 25 09:18:00 2018
   
4. 
   
5. """
   
6. 
   
7. import numpy as np
   
8. import  matplotlib.pyplot as plt
   
9. from scipy.io import loadmat
   
10. from matplotlib.patches import Polygon
    
11. from scipy.interpolate import griddata
    
12. from mpl_toolkits.basemap import Basemap
    
13. from matplotlib.font_manager import FontProperties
    
14. font=FontProperties(fname=r'c:\windows\fonts\simsun.ttc',size=14)
    
15. lonlat=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\lonlat.mat')
    
16. AnnualRain_Ratio=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\AnnualRain_Ratio.mat')
    
17. Rain_Ratio=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\Rain_Ratio.mat')
    
18. StnName=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\StnName.mat')
    
19. sample=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\sample.mat')
    
20. Token=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\Token.mat')
    
21. New_Rain_Ratio=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\New_Rain_Ratio.mat')
    
22. grid_lon,grid_lat=np.meshgrid(np.arange(105,110.5,0.1),np.arange(28,32.5,0.1))
    
23. lon_lat=np.transpose(np.array([lonlat['lonlat'][:,0],lonlat['lonlat'][:,1]]))
    
24. fig=plt.figure()
    
25. for i in range(6):
    
26.     names=locals()
    
27.     names['ax%s'%(i+1)]=fig.add_subplot(2,3,i+1)
    
28. #    ax0=fig.add_axes([0.85,0.10,0.04,0.80])
    
29.     map=Basemap(projection='cyl',resolution='l',llcrnrlon=105,llcrnrlat=28,urcrnrlon=111,urcrnrlat=33,ax=names['ax%s'%(i+1)])
    
30.     shapeinfo=map.readshapefile('F:\\2016\\Youth_Foundation\\CHN_adm_shp\\CHN_adm3','states',drawbounds=False)
    
31.     for info, shp in zip(map.states_info,map.states):
    
32.         proid=info['NAME_1']
    
33.         if proid=='Chongqing':
    
34.            poly=Polygon(shp,facecolor='w',edgecolor='k',lw=0.5)
    
35.            names['ax%s'%(i+1)].add_patch(poly)
    
36.    map.drawparallels(np.arange(28,33,2),color='k',labels=[1,0,0,0])
    
37.    map.drawmeridians(np.arange(105,110,2),color='k',labels=[0,0,0,1])
    
38.    agency=griddata(lon_lat,AnnualRain_Ratio['AnnualRain_Ratio'][:,i],(grid_lon,grid_lat),method='cubic')
    
39.    cf=map.contourf(grid_lon,grid_lat,agency)
    
40. #    cb=map.colorbar(cf,'right',size='5%')
    
41. #    plt.title(u'降水频率')

复制代码