请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9178|回复: 0

[求助] 画图时阴影图被地图遮盖

[复制链接]

新浪微博达人勋

发表于 2018-2-5 16:49:50 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 schliezer 于 2018-2-5 16:54 编辑

地图遮住了阴影,但我明明是先画的地图,另外,有没有办法可以只保留第一列子图的坐标?所有坐标挤在一起很难看。代码和图:
  1. 2# -*- coding: utf-8 -*-
  2. """
  3. Created on Thu Jan 25 09:18:00 2018

  4. """

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




Figure_1.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表