爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 22028|回复: 33

MeteoInfoLab脚本示例:CALIPSO Backscatter data

[复制链接]

新浪微博达人勋

发表于 2016-2-19 12:42:42 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 MeteoInfo 于 2016-2-19 14:12 编辑

读取和绘制CALIPSO卫星遥感数据中的Total Attenuated Backscatter变量。

  1. # Add file
  2. fn = 'D:/Temp/hdf/CAL_LID_L1-ValStage1-V3-01.2007-06-12T03-42-18ZN.hdf'
  3. f = addfile(fn)

  4. # Read data
  5. vname = 'Total_Attenuated_Backscatter_532'
  6. var = f[vname]
  7. data = var[:1000,:]
  8. data = rot90(data)
  9. lats = f['latitude'][:1000,0]
  10. latstrs = []
  11. for lat in lats:
  12.     latstrs.append('%.1f' % lat)
  13. height = f['metadata'][None]
  14. m = height.findMember('Lidar_Data_Altitudes')
  15. h = array(height.getArray(0, m))
  16. h = h[::-1]
  17. hstrs = []
  18. for hh in h:
  19.     hstrs.append('%.1f' % hh)

  20. # Plot
  21. levs = [0.0001,0.0002,0.0003,0.0004,0.0005,0.0006,0.0007,0.0008,0.0009,\
  22.     0.001,0.0015,0.002,0.0025,0.003,0.0035,0.004,0.0045,0.005,0.0055,0.006,\
  23.     0.0065,0.007,0.0075,0.008,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,01]
  24. layer = imshow(data, levs, cmap='MPL_terrain')
  25. xticks(data.dimvalue(1), latstrs)
  26. yticks(data.dimvalue(0), hstrs)
  27. xlabel('Latitude')
  28. ylabel('Altitude (km)')
  29. colorbar(layer, label=r'$\rm{km}^{-1}$ \rm{sr}$^{-1}')
  30. basename = os.path.basename(fn)
  31. title('{0}\n{1}'.format(basename, vname))


calipso_l1.png

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-2-19 15:46:38 | 显示全部楼层
本帖最后由 MeteoInfo 于 2019-10-30 15:47 编辑

基于CAPLISO官网上MatLab例子中的kathys_lidar_colors.m程序,写了一个生成和官网色标一致的色标生成脚本。本例中使用的数据文件来自此网页:http://ccplot.org/api/ ,利用新色标生成的图形如下:
calipso_l1.png

脚本程序文件:   CALIPSO_L1.py (1.46 KB, 下载次数: 34)

CALIPSO_colors.py

5.04 KB, 下载次数: 21, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-3-15 08:01:36 | 显示全部楼层
谢谢王老师分享,漂亮哈
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-17 16:42:47 | 显示全部楼层
请问楼主这句代码“from mipylib.midata import array”如何实现啊,我运行出来总是提示 ImportError: No module named midata
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-17 18:17:48 | 显示全部楼层
2020Doraemon 发表于 2016-4-17 16:42
请问楼主这句代码“from mipylib.midata import array”如何实现啊,我运行出来总是提示 ImportError: No m ...

midata改为minum
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-18 10:24:37 | 显示全部楼层
请问楼主,L1的后向散射数据,质量控制参数如何设置呀?我才开始学,还请老师不吝赐教,谢谢谢谢~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-18 13:47:49 | 显示全部楼层
直来直往 发表于 2016-4-18 10:24
请问楼主,L1的后向散射数据,质量控制参数如何设置呀?我才开始学,还请老师不吝赐教,谢谢谢谢~ ...

我也没有真正用过此类数据,帮不了你
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-18 20:44:16 | 显示全部楼层

谢谢楼主!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-21 16:09:56 | 显示全部楼层
MeteoInfo 发表于 2016-4-18 13:47
我也没有真正用过此类数据,帮不了你

好吧,还是谢谢啦~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-28 10:53:24 | 显示全部楼层
老师 ,这个可以用于2级数据吗?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表