雷达基数据三维格点插值及体绘制

MeteoInfo

生成 x, y, z 格点数据坐标，利用 get_grid_3d_data 函数将多层扫描的雷达基数据插值为三维格点数据，然后利用volumplot函数进行三维体绘制。

1. import time
   
2. 
   
3. st = time.clock()
   
4. fn = r'D:\Temp\binary\radar\Z_RADR_I_Z9898_20190828181529_O_DOR_SAD_CAP_FMT.bin'
   
5. f = addfile(fn)
   
6. slat = f.attrvalue('StationLatitude')[0]
   
7. slon = f.attrvalue('StationLongitude')[0]
   
8. h = f.attrvalue('AntennaHeight')[0]
   
9. x = arange(-200000, 200001, 1000)
   
10. y = arange(-200000, 200001, 1000)
    
11. z = arange(1000, 20001, 500)
    
12. 
    
13. #Get 3d grid data
    
14. data = f.get_grid_3d_data('dBZ', x, y, z, h)
    
15. 
    
16. #Set azimuth equidistant projection
    
17. proj = projinfo(proj='aeqd', lon_0=slon, lat_0=slat)
    
18. 
    
19. # Plot test
    
20. axes3d(projection=proj, orthographic=False, aspect='xy_equal', facecolor='k',
    
21.     fgcolor='w', axis=True, clip_plane=False)
    
22. geoshow('world_topo.jpg')
    
23. geoshow('cn_province')
    
24. volumeplot(x, y, z, data, vmin=0, vmax=60, cmap='NCV_bright')
    
25. colorbar(shrink=0.8, tickcolor='w')
    
26. xlim(x[0], x[-1])
    
27. ylim(y[0], y[-1])
    
28. zlim(0, z[-1])
    
29. xlabel('Longitude')
    
30. ylabel('Latitude')
    
31. zlabel('Height (m)')
    
32. 
    
33. et = time.clock()
    
34. print('Running time: {} Seconds'.format(et-st))

import time<br /> <br /> st = time.clock()<br /> fn = r'D:\Temp\binary\radar\Z_RADR_I_Z9898_20190828181529_O_DOR_SAD_CAP_FMT.bin'<br /> f = addfile(fn)<br /> slat = f.attrvalue('StationLatitude')[0]<br /> slon = f.attrvalue('StationLongitude')[0]<br /> h = f.attrvalue('AntennaHeight')[0]<br /> x = arange(-200000, 200001, 1000)<br /> y = arange(-200000, 200001, 1000)<br /> z = arange(1000, 20001, 500)<br /> <br /> #Get 3d grid data<br /> data = f.get_grid_3d_data('dBZ', x, y, z, h)<br /> <br /> #Set azimuth equidistant projection<br /> proj = projinfo(proj='aeqd', lon_0=slon, lat_0=slat)<br /> <br /> # Plot test<br /> axes3d(projection=proj, orthographic=False, aspect='xy_equal', facecolor='k',<br />     fgcolor='w', axis=True, clip_plane=False)<br /> geoshow('world_topo.jpg')<br /> geoshow('cn_province')<br /> volumeplot(x, y, z, data, vmin=0, vmax=60, cmap='NCV_bright')<br /> colorbar(shrink=0.8, tickcolor='w')<br /> xlim(x[0], x[-1])<br /> ylim(y[0], y[-1])<br /> zlim(0, z[-1])<br /> xlabel('Longitude')<br /> ylabel('Latitude')<br /> zlabel('Height (m)')<br /> <br /> et = time.clock()<br /> print('Running time: {} Seconds'.format(et-st))

威武霸气的潘潘

感谢楼主分享！

ayzqs

功能越来越强大了

司机是猪

beiruoyiqiao

感谢分享！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！

我是sswwkk

请问这是什么编程软件啊

我爱气象111

我运行提示没有times这个函数库，应该怎么安装呢

MeteoInfo

我是sswwkk 发表于 2024-9-18 16:27
请问这是什么编程软件啊

参考这里：http://www.meteothink.org/docs/meteoinfolab/milab_cn/index.html

MeteoInfo

我爱气象111 发表于 2024-9-19 08:20
我运行提示没有times这个函数库，应该怎么安装呢

是 time  ，不是  times

我爱气象111

MeteoInfo 发表于 2024-9-19 11:10
是 time  ，不是  times

难怪，哈哈哈哈，谢谢大佬