已解决-python读取风云2号卫星TBB数据出图报错

蔡小亦

已解决，经纬度数据和亮温数据对应有问题。


参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=56315获取TBB数据，用cartopy对TBB数据可视化。
FY2E、FY2G数据均能作等值线图，见附件。但作填色图就报错：ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer
搜了很久还是没法解决，代码如下，请大神指教。

1. import h5py
   
2. import numpy as np
   
3. import matplotlib.pyplot as plt
   
4. import cartopy.crs as ccrs
   
5. import pandas as pd
   
6. 
   
7. f = h5py.File('./data/TBB/FY2G_TBB_IR1_NOM_20230507_0400.hdf', 'r')
   
8. for key in f.keys():
   
9.     print(key)
   
10. 
    
11. a = f.get('FY2G TBB Hourly Product')
    
12. data = np.array(a)
    
13. data = np.ma.masked_outside(data,160,340)
    
14. info = f.get('NomFileInfo')[0]
    
15. print(info)
    
16. 
    
17. lonCenter = info[4]
    
18. loc = np.fromfile('./data/TBB/NOM_ITG_2288_2288(0E0N)_LE.dat',dtype=np.float32,count=-1)
    
19. latlon = loc.reshape(2,2288,2288)
    
20. lon = np.array(latlon[0]+lonCenter)
    
21. lat = np.array(latlon[1])
    
22. lon = np.ma.masked_where(lon>300,lon)
    
23. lat = np.ma.masked_values(lat,300)
    
24. 
    
25. lon_df = pd.DataFrame(lon)
    
26. lat_df = pd.DataFrame(lat)
    
27. data_df = pd.DataFrame(data)
    
28. 
    
29. # 绘图
    
30. fig=plt.figure(figsize=(12,8))
    
31. proj = ccrs.PlateCarree()
    
32. ax = fig.add_subplot(1,1,1,projection = proj)
    
33. ax.set_extent([105,120,18,28],crs=ccrs.PlateCarree())
    
34. ax.coastlines()
    
35. #levels = np.arange(160, 340 ,10),cmap = 'YlGnBu_r'
    
36. 
    
37. ax.contour(lon_df,lat_df,data_df)
    
38. 
    
39. # 填色图报错？？？
    
40. # ax.contourf(lon_df,lat_df,data_df,levels=np.linspace(np.min(data),np.max(data),400) )
    
41. 
    
42. plt.savefig('tbb.png',format='png',dpi=200)
    
43. print('done')

复制代码

cxf2524

可以问一下如何解决的吗？

蔡小亦

cxf2524 发表于 2023-10-26 21:08
可以问一下如何解决的吗？

有的经纬度点是没有亮温值的（可能是因为在边缘），所以要把这些值剔除。