MeteoInfoLab脚本示例：GPM Swath数据

MeteoInfo

GPM (Global Precipitation Measurement) 卫星swath数据读取及绘图的例子：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = '1C.F19.SSMIS.XCAL2015-P.20160105-S214106-E232259.009078.V03A.HDF5'
   
4. fn = folder + fns
   
5. f = addfile(fn)
   
6. lon = f['Longitude'][:,:]
   
7. lat = f['Latitude'][:,:]
   
8. vname = 'Tc'
   
9. v_data = f[vname]
   
10. data = v_data[:,:,0]
    
11. data[data<=-9999.9] = nan
    
12. long_name = v_data.attrvalue('LongName')[0]
    
13. units = v_data.attrvalue('Units')[0]
    
14. #Plot
    
15. axesm()
    
16. mlayer = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(mlayer, edgecolor='k')
    
18. levs = arange(40, 90, 1)
    
19. layer = pcolorm(lon, lat, data, 20)
    
20. colorbar(layer, orientation='horizontal', aspect=40, label=units)
    
21. title('{0}{1}'.format(fns, long_name+' (nchannel1=0)'))
    
22. axism()

#Add data file<br /> folder = 'D:/Temp/hdf/'<br /> fns = '1C.F19.SSMIS.XCAL2015-P.20160105-S214106-E232259.009078.V03A.HDF5'<br /> fn = folder + fns<br /> f = addfile(fn)<br /> lon = f['Longitude'][:,:]<br /> lat = f['Latitude'][:,:]<br /> vname = 'Tc'<br /> v_data = f[vname]<br /> data = v_data[:,:,0]<br /> data[data<=-9999.9] = nan<br /> long_name = v_data.attrvalue('LongName')[0]<br /> units = v_data.attrvalue('Units')[0]<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer, edgecolor='k')<br /> levs = arange(40, 90, 1)<br /> layer = pcolorm(lon, lat, data, 20)<br /> colorbar(layer, orientation='horizontal', aspect=40, label=units)<br /> title('{0}{1}'.format(fns, long_name+' (nchannel1=0)'))<br /> axism()

freekiller

强大                     

风之精灵

太厉害了   

kmdqlb

王老师威武

hxl110110

哇三，高大山！

aas314

好厉害的图  学习了  

刘迪

我是兰州大学的学生，最近在学GPM卫星数据，遇到一些问题，想和您交流学习一下，就是最简单的经纬度的维数我还没有看懂，不知道是怎么对应起来的，方便留个联系方式么？谢谢您了！

MeteoInfo

刘迪 发表于 2018-1-4 13:36
我是兰州大学的学生，最近在学GPM卫星数据，遇到一些问题，想和您交流学习一下，就是最简单的经纬度的维数 ...

论坛不就是交流的地方吗？你也可以加入MeteoInfo QQ群交流。

刘迪

MeteoInfo 发表于 2018-1-4 20:09
论坛不就是交流的地方吗？你也可以加入MeteoInfo QQ群交流。

王老师，我现在用GPM卫星的Ku波段的二级数据在挑选深对流个例，之前在文献上看到说是 ：雷达反射率因子达到20dBZ的高度大于15Km 就可以确定为深对流个例，我现在在编写程序，有一些疑问，我看了GPM的数据说明，但还是不太确定表示雷达反射率因子以及高度的量是哪一个，请老师指点一下，谢谢您！

MeteoInfo

刘迪 发表于 2018-1-7 11:43
王老师，我现在用GPM卫星的Ku波段的二级数据在挑选深对流个例，之前在文献上看到说是 ：雷达反射率因子达 ...

脚本只是示例读取绘制相关数据的能力，数据本身并没有研究过，相信数据的说明文档应该能说清楚吧。