做区域填色图时出错

mockings

因为刚接触python，想在青藏高原地区绘制站点温度的填色图，在参考了气象家园和官网上basemap的教程后，自己编写了程序，但是一直报错

1. load_sta = 'e:/data/data.mat'
   
2. station = sio.loadmat(load_sta)
   
3. #print (type(station))
   
4. #print station.keys()
   
5. #print (station['data'])
   
6. data=station['data']   
   
7. #print (data[1,1])
   
8. lat=data[:,1]
   
9. lon=data[:,2]
   
10. #print (lon)
    
11. lon,lat=np.meshgrid(lon,lat)  #网格数据
    
12. 
    
13. 
    
14. 
    
15. # 画图
    
16. fig = plt.figure(figsize=(16,16))
    
17. # 设置全局字体等
    
18. plt.rc('font',size=15,weight='bold')
    
19. plt.subplot(121)
    
20. m = Basemap(projection='cyl',llcrnrlat=25,llcrnrlon=70,urcrnrlat=45,urcrnrlon=110)
    
21. m = Basemap(projection='cyl',llcrnrlat=25,llcrnrlon=70,urcrnrlat=45,urcrnrlon=110)
    
22. # 加载中国地图文件
    
23. m.readshapefile('e:\\python\\CHN_adm_shp\\CHN_adm0','CHN_adm0.shp',linewidth=1,color ='black')
    
24. # 加载青藏高原边界文件
    
25. m.readshapefile('e:\\python\\DBATP/DBATP_Line','DBATP_Line.shp',linewidth = 2,color='red')
    
26. x,y=m(lon,lat)
    
27. 
    
28. levels = np.linspace(3000,np.max(juping),50)
    
29. cf =m.contourf(x,y,juping,levels=levels)   ##contourf(x,y,数据，参数，)如果默认true则x，y成经纬度
    
30. cbar = m.colorbar(cf,location='bottom',ticks=np.linspace(3000,np.max(juping),10),label='m')
    
31. 
    
32. 
    
33. lon_num = np.arange(70,111,5)
    
34. lon_name = ['70','75','80','85','90','95','100','105','110']
    
35. lat_num =  np.arange(25,46,5)
    
36. lat_name = ['25','30','35','40','45']
    
37. plt.yticks(lat_num,lat_name)
    
38. plt.xticks(lon_num,lon_name)

复制代码

data是站点的经纬度资料，程序我判断（因为改为批注后就没错了）是在cf=m.contourf  那一句出错，报错信息是：UnicodeDecodeError: 'ascii' codec can't decode byte 0xe6 in position 58: ordinal not in range(128)。     我百度后说是因为数据文件里面有中文字符，但是我juping是一个自建的纯数字数组，然后尝试了在程序开头添加import sys
reload(sys)
sys.setdefaultencoding('utf-8')，仍然不管用。
请求各位大神帮忙，非常感谢！！

Lancelot

那就print juping看看

mockings

Lancelot 发表于 2018-4-2 21:18
那就print juping看看

1. ave=np.zeros((71,12))               #建立（71,12）初始元素为0数组
   
2. juping=np.zeros(((71,12,55)))
   
3. for i in range (0,70):               
   
4.    for j in range (0,11):
   
5.       ave[i,j]=np.mean(c[i,j,:])  #求站点每月的年平均
   
6.       for k in range (0,54):
   
7.           juping[i,j,k]=c[i,j,k]-ave[i,j]

复制代码

感谢感谢，确实是juping数据出了问题。这是我求juping的代码。当时检查的时候ave前几个数据都没有出错，今天检查时发现ave后几个数据均为0，如果单独np.mean(c[i,j,:]）则可以正确求出平均，不知道是哪里出了问题。

Lancelot

mockings 发表于 2018-4-3 09:37
感谢感谢，确实是juping数据出了问题。这是我求juping的代码。当时检查的时候ave前几个数据都没有出错 ...

对缺测的处理是什么 mask吗？
自己检查一下缺测的处理方法和对应函数
还有 我不懂你这个ave的后几个数据指什么 请自己确认一下