meteoinfo处理最新的calipso数据没有Feature_Classification_Flags变量

chongzika

http://mailman.ucar.edu/pipermai ... January/011273.html
http://mailman.ucar.edu/pipermai ... January/011348.html
同样ncl也有如此的问题；貌似是数据已经更新了？https://catalog.data.gov/dataset/cal-lid-l2-vfm-standard-v4-10
--------------------------------------------------------------
脚本是画[size=13.3333px]CAL_LID_L2_VFM-Standard-V4-10的脚本： http://hdfeos.org/forums/showthread.php?t=710&page=3
Feature_Classification_Flags is not a variable name
Traceback (most recent call last):
  File "vfm_meteoinfo1.py", line 8, in <module>
    var = f[vname]
  File "/data/c03n02/cliu/software/meteinfo/MeteoInfo-1.4.8R7/MeteoInfo/pylib/mipylib/dataset/dimdatafile.py", line 45, in __getitem__
    raise ValueError()
ValueError

独孤酒见

王老师很厉害，太服了。谢谢王老师。我最近一直在鼓捣calipso的数据，前几天一直在琢磨您在今年2月1号的报告里的读后散的脚本还有就这个帖子里的脚本。好像隐约琢磨出点味道了。
第26行    alt<i> = -0.5 + i*0.03有个小问题，把<>改为（）

我现在琢磨两个问题，请王老师指导。
1、横轴坐标能把经度加上吗？
2、在那个后散的脚本里面，能否在图上叠加上另外的一个变量，比如Surface_Wind_Speeds等这种变量

请王老师在得空的时候指导，谢谢王老师。

MeteoInfo

你附上的数据里是有 Feature_Classification_Flags 变量的

独孤酒见

楼主搞定没？我现在也在搞，一头雾水，读不出来啊。

MeteoInfo

独孤酒见 发表于 2018-4-24 13:50
楼主搞定没？我现在也在搞，一头雾水，读不出来啊。

1. # Add file
   
2. fn = 'D:/Temp/hdf/CAL_LID_L2_VFM-Standard-V4-10.2017-06-09T18-19-56ZN.hdf'
   
3. f = addfile(fn)
   
4. 
   
5. # Read data
   
6. vname = 'Feature_Classification_Flags'
   
7. var = f[vname]
   
8. data = var[:,:]
   
9. lat = f['Latitude'][:,0]
   
10. 
    
11. # Extract Feature Type only through bitmask.
    
12. data = data & 7
    
13. 
    
14. # Subset latitude values for the region of interest.
    
15. lat = lat[1000:2000]
    
16. size = lat.shape[0]
    
17. 
    
18. data2d = data[1000:2000, 1165:]  # -0.5km to  8.2km
    
19. data3d = reshape(data2d, (size, 15, 290))
    
20. data = data3d[:,0,:]
    
21. 
    
22. # Generate altitude data according to file specification [1].
    
23. alt = zeros(290)
    
24. # -0.5km to 8.2km
    
25. for i in range (0, 290):
    
26.     alt = -0.5 + i*0.03
    
27. 
    
28. # Plot
    
29. levs = arange(8)
    
30. cols = [(255,255,255),(0,0,255),(51,255,255),(255,153,0),(255,255,0),(0,255,0),(127,127,127),(0,0,0)]
    
31. ls = makesymbolspec('image', levels=levs, colors=cols)
    
32. layer = imshow(rot90(data, 1)[:,::-1], symbolspec=ls, extent=[lat[-1],lat[0],alt[0],alt[-1]])
    
33. colorbar(layer, ticklabels=['Invalid', 'Clear Air', 'Cloud', 'Aerosol', 'Strato Feature', 'Surface', 'Subsurface', 'No Signal'])
    
34. basename = os.path.basename(fn)
    
35. title([basename, 'Feature Type (Bits 1-3) in Feature Classification Flag'])
    
36. xlabel('Latitude (degrees north)')
    
37. ylabel('Altitude (km)')
    
38. xaxis(tickin=False)
    
39. yaxis(tickin=False)

# Add file<br /> fn = 'D:/Temp/hdf/CAL_LID_L2_VFM-Standard-V4-10.2017-06-09T18-19-56ZN.hdf'<br /> f = addfile(fn)<br /> <br /> # Read data<br /> vname = 'Feature_Classification_Flags'<br /> var = f[vname]<br /> data = var[:,:]<br /> lat = f['Latitude'][:,0]<br /> <br /> # Extract Feature Type only through bitmask.<br /> data = data & 7<br /> <br /> # Subset latitude values for the region of interest.<br /> lat = lat[1000:2000]<br /> size = lat.shape[0]<br /> <br /> data2d = data[1000:2000, 1165:]  # -0.5km to  8.2km<br /> data3d = reshape(data2d, (size, 15, 290))<br /> data = data3d[:,0,:]<br /> <br /> # Generate altitude data according to file specification [1].<br /> alt = zeros(290)<br /> # -0.5km to 8.2km<br /> for i in range (0, 290):<br />     alt<i> = -0.5 + i*0.03<br /> <br /> # Plot<br /> levs = arange(8)<br /> cols = [(255,255,255),(0,0,255),(51,255,255),(255,153,0),(255,255,0),(0,255,0),(127,127,127),(0,0,0)]<br /> ls = makesymbolspec('image', levels=levs, colors=cols)<br /> layer = imshow(rot90(data, 1)[:,::-1], symbolspec=ls, extent=[lat[-1],lat[0],alt[0],alt[-1]])<br /> colorbar(layer, ticklabels=['Invalid', 'Clear Air', 'Cloud', 'Aerosol', 'Strato Feature', 'Surface', 'Subsurface', 'No Signal'])<br /> basename = os.path.basename(fn)<br /> title([basename, 'Feature Type (Bits 1-3) in Feature Classification Flag'])<br /> xlabel('Latitude (degrees north)')<br /> ylabel('Altitude (km)')<br /> xaxis(tickin=False)<br /> yaxis(tickin=False)

独孤酒见

MeteoInfo 发表于 2018-4-24 18:02

王老师很厉害，太服了。谢谢王老师。我最近一直在鼓捣calipso的数据，前几天一直在琢磨您在今年2月1号的报告里的读后散的脚本还有就这个帖子里的脚本。好像隐约琢磨出点味道了。
第26行    alt<i> = -0.5 + i*0.03有个小问题，把<>改为（）

我现在琢磨两个问题，请王老师指导。
1、横轴坐标能把经度加上吗？
2、在那个后散的脚本里面，能否在图上叠加上另外的一个变量，比如Surface_Wind_Speeds等这种变量

请王老师在得空的时候指导，谢谢王老师。

lwqqqqq

想问下第20行代码data = data3d[:,0,:]，为什么是0，不是0-14中的其他值？

天天paper

独孤酒见 发表于 2018-4-26 10:27
王老师很厉害，太服了。谢谢王老师。我最近一直在鼓捣calipso的数据，前几天一直在琢磨您在今年2月1号的报 ...

老师，您好，您的横坐标上把经度加上了吗？我也遇到了这样的问题

晴天

天天paper 发表于 2021-2-16 12:33
老师，您好，您的横坐标上把经度加上了吗？我也遇到了这样的问题

参考王老师后向散射的那个程序可以把经度加上