Python读取nc数据和绘图

lorraineni

{:eb303:}{:eb303:}{:eb303:}感谢楼主无私分享

Orient2016

谢谢 下载下来好好学习

Orient2016

Orient2016

{:eb303:}{:eb303:}

wency

画填色加等值线图的程序里，为什么要写lats=lats[::-1]这一句呢

平流层的萝卜

wency 发表于 2018-1-24 14:59
画填色加等值线图的程序里，为什么要写lats=lats[::-1]这一句呢

南北倒置，目的让 lats里高纬在前，低纬在后

mockings

学习了，感谢楼主分享

mockings

请教下楼主，我刚开始学习Python，读取nc文件中obj.variables中的variables这个属性的定义是在netcdf4里面默认的吗？我该怎么查找这个属性的定义。多谢楼主

schliezer

你好，我在叠图的时候遇到了点问题，打扰请教一下。第一个是不知道怎么把经纬度网格线去掉（看起来碍眼睛，原谅我的强迫症），mdrawparallels和mdrawmeridians函数似乎都没有提供可以直接操作的方法。第二个不知道怎么叠地图，虽然我单独可以正确画出来地图。网上也没有找到明确的示例程序。我放一下自己的程序和出的图：

1. import numpy as np
   
2. import  matplotlib.pyplot as plt
   
3. from scipy.io import loadmat
   
4. from matplotlib.patches import Polygon
   
5. from scipy.interpolate import griddata
   
6. from mpl_toolkits.basemap import Basemap
   
7. lonlat=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\lonlat.mat')
   
8. AnnualRain_Ratio=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\AnnualRain_Ratio.mat')
   
9. Rain_Ratio=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\Rain_Ratio.mat')
   
10. StnName=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\StnName.mat')
    
11. sample=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\sample.mat')
    
12. Token=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\Token.mat')
    
13. New_Rain_Ratio=loadmat('C:\Users\CQQXT\Documents\MATLAB\BP_forecast_onebyone\New_Rain_Ratio.mat')
    
14. grid_lon,grid_lat=np.meshgrid(np.arange(105,110.5,0.1),np.arange(28,32.5,0.1))
    
15. lon_lat=np.transpose(np.array([lonlat['lonlat'][:,0],lonlat['lonlat'][:,1]]))
    
16. def PlotCqMap():
    
17.    fig=plt.figure()
    
18.    ax1=fig.add_axes([0.1,0.1,0.8,0.8])
    
19.    bmap=Basemap(llcrnrlon=105,llcrnrlat=28,urcrnrlon=111,urcrnrlat=33)
    
20.    shp_info=bmap.readshapefile('F:\\2016\\Youth_Foundation\\CHN_adm_shp\\CHN_adm3','states',drawbounds=False)
    
21.    for info, shp in zip(bmap.states_info,bmap.states):
    
22.        proid=info['NAME_1']
    
23.        if proid=='Chongqing':
    
24.            poly=Polygon(shp,facecolor='w',edgecolor='k',lw=0.5)
    
25.            ax1.add_patch(poly)
    
26.    bmap.drawparallels(np.arange(28,33,2),color='k',labels=[1,0,0,0])
    
27.    bmap.drawmeridians(np.arange(105,110,2),color='k',labels=[0,0,0,1])
    
28.    return
    
29. for i in range(6):
    
30.     plt.subplot(2,3,i+1)
    
31.     agency=griddata(lon_lat,AnnualRain_Ratio['AnnualRain_Ratio'][:,i],(grid_lon,grid_lat),method='cubic')
    
32.     plt.imshow(agency)
    
33.     PlotCqMap
    
34.     plt.title('cubic interpolation')
    

复制代码

hui86613296

真心不错 很好啊