python如何读取Micaps第四类数据

Eegle

请问大家，python该如何读取Micaps中的第四类数据呢？
搜了一下，河南新乡气象局的朱枫同学写过一个Diamond4的库
但，感觉说的不是很明白，也不够好用的，很多没看明白的地方。
在此，求助大神们。
感谢！
祝您学习愉快，一切安康！

Masterpiece

1. import numpy as np
   
2. from dateutil.parser import *
   
3. import scipy.ndimage as sn
   
4. import matplotlib.pyplot as plt
   
5. from mpl_toolkits.basemap import Basemap
   
6. 
   
7. file_name='17033020.000'
   
8. raw=np.loadtxt(file_name,dtype=np.str,delimiter='/n')
   
9. 
   
10. raw=' '.join(raw.tolist())
    
11. raw=raw.split()
    
12. 
    
13. title=raw[:3]
    
14. 
    
15. time_str='20'+'-'.join(raw[3:3+4])
    
16. file_time=parse(time_str)
    
17. duration,level=int(raw[7]),int(raw[8])
    
18. xint,yint=float(raw[9]),float(raw[10])
    
19. start_lon,end_lon=float(raw[11]),float(raw[12])
    
20. start_lat,end_lat=float(raw[13]),float(raw[14])
    
21. xsize,ysize=int(raw[15]),int(raw[16])
    
22. lineint,start_level,end_level=int(raw[17]),int(raw[18]),int(raw[19])
    
23. smooth,boldlevel=int(raw[20]),raw[21]
    
24. 
    
25. data=np.array(raw[22:],dtype=float).reshape(ysize,xsize)
    
26. 
    
27. fig, ax = plt.subplots(figsize=(14,9))
    
28. 
    
29. m = Basemap(projection='cyl',
    
30.             llcrnrlat=end_lat, urcrnrlat=start_lat,
    
31.             llcrnrlon=start_lon, urcrnrlon=end_lon,
    
32.             resolution='l')
    
33. m.drawcoastlines(linewidth=0.3, color='black')
    
34. 
    
35. lon=np.linspace(start_lon,end_lon,xsize)
    
36. lat=np.linspace(end_lat,start_lat,ysize)[::-1]
    
37. 
    
38. if smooth:
    
39.     factor=6
    
40. else:
    
41.     factor=1
    
42. 
    
43. lons, lats = np.meshgrid(sn.zoom(lon, factor), sn.zoom(lat, factor))
    
44. x, y = m(lons, lats)
    
45. c1=m.contour(x,y,sn.zoom(data, factor),
    
46.              levels=np.arange(start_level,end_level+lineint,lineint),colors='k',zorder=1)
    
47. ax.clabel(c1, fmt='%d', inline=True, fontsize=12, inline_spacing=8)
    
48. 
    
49. m.readshapefile('D:/pyproj/shp/CHN_2019/province', 'province', color='r',linewidth=0.6,zorder=0)
    
50. 
    
51. plt.title('%s      Level:%d'%(file_time.strftime('%Y-%m-%d %H:%M'),level),size=14, loc='left')
    
52. plt.savefig(file_name+'.png',dpi=300,bbox_inches='tight')
    
53. 
    
54. plt.show()

复制代码

Eegle

Masterpiece 发表于 2020-5-27 23:28

好的 感谢。我自己试试看

天边的那边

都是大神呀，佩服

Rainch

  Basemap库安装不上

Rainch

Rainch 发表于 2020-11-20 15:00
Basemap库安装不上

windows系统下安装Basemap库 需要查看当前python支持的版本后，下载相应的安装文件，否则 会出现 is not a supported wheel on this platform 错误提示。
只向下兼容。如果python是3.8版本的，你下载 （basemap-1.2.2-cp39-cp39-win_amd64.whl ）是徒劳的。
只能下载basemap-1.2.2-cp38-cp38-win_amd64.whl
查看命令：（适用python3.X）
pip debug --verbose

hanqianli

if smooth:
    factor=6
else:
    factor=1


请问楼主   这句时什么意思