请教关于NC格式的雷达数据读取问题

hbl1021

我下载了NC格式的雷达图，但是不会可视化处理，请各位大神帮忙看一下怎样在编程软件中看到NC文件中的具体数据。下图为下载数据时的雷达图，附件为NC格式文件。
HBR_Z9010_20240614001800.zip (55.88 KB, 下载次数: 2)

2024-6-14 09:24 上传

点击文件名下载附件

MeteoInfo

供参考

1. fn = r'D:\Temp\binary\radar\HBR_Z9010_20240614001800.nc'
   
2. f = addfile(fn)
   
3. v = f['hbr']
   
4. data = v[:]
   
5. lat = v.attrvalue('rlat')[0] / 10000.
   
6. lon = v.attrvalue('rlon')[0] / 10000.
   
7. height = v.attrvalue('rhei')[0]
   
8. scale = v.attrvalue('scale')[0]
   
9. res = v.attrvalue('res')[0]
   
10. 
    
11. data = data / scale
    
12. data[data<-50] = nan
    
13. 
    
14. #Set azimuth equidistant projection
    
15. proj = projinfo(proj='aeqd', lon_0=lon, lat_0=lat)
    
16. yn, xn = data.shape
    
17. xmin = -(xn / 2 * res)
    
18. ymin = -(yn / 2 * res)
    
19. x = arange1(xmin, xn, res)
    
20. y = arange1(ymin, yn, res)
    
21. 
    
22. #Plot        
    
23. ax = axesm(facecolor='b', projection=proj)
    
24. geoshow('cn_province', edgecolor=None, facecolor=[230,230,230])
    
25. geoshow('cn_province', edgecolor=[80,80,80])
    
26. city = geoshow('cn_cities', facecolor='r', size=8)
    
27. city.addlabels('NAME', fontname=u'黑体', fontsize=16, yoffset=18)
    
28. levs = [5,10,15,20,25,30,35,40,45,50,55,60,65,70]
    
29. cols = [(255,255,255,0),(102,255,255),(0,162,232),(86,225,250),(3,207,14),
    
30.     (26,152,7),(255,242,0),(217,172,113),(255,147,74),(255,0,0),
    
31.     (204,0,0),(155,0,0),(236,21,236),(130,11,130),(184,108,208)]
    
32. imshow(x, y, data, levs, colors=cols, proj=proj, zorder=1)
    
33. colorbar(shrink=0.8, label='dBZ', labelloc='top')
    
34. #Plot circles
    
35. cx, cy = project(lon, lat, toproj=proj)
    
36. rr = array([50, 100, 150, 200, 250, 300]) * 1000
    
37. for rd in rr:
    
38.     ax.add_circle((cx, cy), rd, edgecolor='r')
    
39. plot([cy,cy], [cx-rd,cx+rd], color='r', proj=proj)
    
40. plot([cy+rd,cy-rd], [cx,cx], color='r', proj=proj)
    
41. xlim(cx - rd, cx + rd)
    
42. ylim(cy - rd, cy + rd)

fn = r'D:\Temp\binary\radar\HBR_Z9010_20240614001800.nc'<br /> f = addfile(fn)<br /> v = f['hbr']<br /> data = v[:]<br /> lat = v.attrvalue('rlat')[0] / 10000.<br /> lon = v.attrvalue('rlon')[0] / 10000.<br /> height = v.attrvalue('rhei')[0]<br /> scale = v.attrvalue('scale')[0]<br /> res = v.attrvalue('res')[0]<br /> <br /> data = data / scale<br /> data[data<-50] = nan<br /> <br /> #Set azimuth equidistant projection<br /> proj = projinfo(proj='aeqd', lon_0=lon, lat_0=lat)<br /> yn, xn = data.shape<br /> xmin = -(xn / 2 * res)<br /> ymin = -(yn / 2 * res)<br /> x = arange1(xmin, xn, res)<br /> y = arange1(ymin, yn, res)<br /> <br /> #Plot        <br /> ax = axesm(facecolor='b', projection=proj)<br /> geoshow('cn_province', edgecolor=None, facecolor=[230,230,230])<br /> geoshow('cn_province', edgecolor=[80,80,80])<br /> city = geoshow('cn_cities', facecolor='r', size=8)<br /> city.addlabels('NAME', fontname=u'黑体', fontsize=16, yoffset=18)<br /> levs = [5,10,15,20,25,30,35,40,45,50,55,60,65,70]<br /> cols = [(255,255,255,0),(102,255,255),(0,162,232),(86,225,250),(3,207,14),<br />     (26,152,7),(255,242,0),(217,172,113),(255,147,74),(255,0,0),<br />     (204,0,0),(155,0,0),(236,21,236),(130,11,130),(184,108,208)]<br /> imshow(x, y, data, levs, colors=cols, proj=proj, zorder=1)<br /> colorbar(shrink=0.8, label='dBZ', labelloc='top')<br /> #Plot circles<br /> cx, cy = project(lon, lat, toproj=proj)<br /> rr = array([50, 100, 150, 200, 250, 300]) * 1000<br /> for rd in rr:<br />     ax.add_circle((cx, cy), rd, edgecolor='r')<br /> plot([cy,cy], [cx-rd,cx+rd], color='r', proj=proj)<br /> plot([cy+rd,cy-rd], [cx,cx], color='r', proj=proj)<br /> xlim(cx - rd, cx + rd)<br /> ylim(cy - rd, cy + rd)

TANG-堂

还有这个，您也可以看一下 不知道是否对您有帮助。https://pycwr.readthedocs.io/_/downloads/en/latest/pdf/