python FY2E卫星HDF5格式的TBB数据画图

chiqu296

1、数据为国家卫星中心下载的HDF5是全球标称数据，2288*2288的格点2、每个格点对应一个一个经纬度。对应的经纬度数据文件可在卫星中心下载
3、读取hdf5数据使用h5py库
4、cmap为使用ncl的色标，论坛里有，自行查找

1. import h5py
   
2. import cmaps
   
3. import numpy as np
   
4. import matplotlib.pyplot as plt
   
5. from mpl_toolkits.basemap import Basemap
   
6. 
   
7. 
   
8. 
   
9. # read hdf5 file
   
10. f = h5py.File('data/FY2E_TBB_IR1_NOM_20170912_1030.hdf','r')
    
11. # 查看文件信息
    
12. # print(list(f.keys()))
    
13. data = f.get('FY2E TBB Hourly Product')
    
14. # print(list(data.attrs.items()))
    
15. data = np.array(data)
    
16. '''
    
17. 处理缺值，数据里有效数据的最大值为340，最小为160
    
18. print(list(data.attrs.items()))即可查看
    
19. '''
    
20. data = np.ma.masked_outside(data,160,340)
    
21. info = f.get('NomFileInfo')[0]
    
22. '''
    
23. 卫星的中心经度 info[4]
    
24. 中心纬度为0
    
25. '''
    
26. lonCenter = info[4]
    
27. # read lonlat msg
    
28. loc = np.fromfile('data/NOM_ITG_2288_2288(0E0N)_LE.dat',dtype=np.float32,count=-1)
    
29. latlon = loc.reshape(2,2288,2288)
    
30. lon = np.array(latlon[0]+lonCenter)
    
31. lat = np.array(latlon[1])
    
32. '''
    
33. 经纬度信息里300为缺值
    
34. '''
    
35. lon = np.ma.masked_where(lon>300,lon)
    
36. lat = np.ma.masked_values(lat,300)
    
37. # 画图
    
38. fig = plt.figure(figsize=(18.5, 10.5))
    
39. m=Basemap(projection='cyl',llcrnrlat=10,llcrnrlon=90,urcrnrlat=35,urcrnrlon=140)
    
40. m.drawcoastlines(color='blue')
    
41. m.drawstates(color='blue')
    
42. m.drawcountries(color='blue')
    
43. x,y = m(lon, lat) # 将lats / lons转换为地图投影坐标
    
44. cf = m.contourf(x,y,data,levels=np.linspace(np.min(data),np.max(data),400),cmap=cmaps.MPL_RdGy )
    
45. cbar = m.colorbar(cf ,location='right',size='3%' ,pad='2%')
    
46. cbar.set_label('Brightness Temperature ( K )', fontsize=12, fontweight='bold')
    
47. m.drawmeridians(np.arange(90, 140, 5),labels=[0, 0, 0, 1])
    
48. m.drawparallels(np.arange(10, 35, 5),labels=[1, 0, 0, 0])
    
49. plt.title('Brightness Temperature(20170912_1030) ',fontweight='bold', fontsize=14)
    
50. plt.savefig('test.png')
    
51. plt.show()

复制代码

tulalang

想问一下，出的图的格式是什么，可以直接附在论文上么

chiqu296

tulalang 发表于 2017-9-29 20:12
想问一下，出的图的格式是什么，可以直接附在论文上么

上面出图为png  matplotlib可以出位图和矢量图，矢量图能满足出版的的要求

wade

地图不规范，有解决方案吗

chiqu296

wade 发表于 2017-9-30 07:16
地图不规范，有解决方案吗

可以加载自己的shapefile文件

tulalang

chiqu296 发表于 2017-9-30 07:14
上面出图为png  matplotlib可以出位图和矢量图，矢量图能满足出版的的要求

那确实很方便。。。赶紧学习python~谢谢楼主的帖子和回复了~

小其其格

赞！就是python自带的地图有点粗

平流层的萝卜

小其其格 发表于 2017-10-1 16:17
赞！就是python自带的地图有点粗

是的，而且python自带的中国地图国界是不标准的。

小其其格

平流层的萝卜 发表于 2017-10-2 15:43
是的，而且python自带的中国地图国界是不标准的。

cmaps我在ubuntu下怎么安装不成功呀！

小肥妞

谢谢楼主分享