- 积分
- 55948
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 MeteoInfo 于 2016-2-19 14:12 编辑
读取和绘制CALIPSO卫星遥感数据中的Total Attenuated Backscatter变量。
- # Add file
- fn = 'D:/Temp/hdf/CAL_LID_L1-ValStage1-V3-01.2007-06-12T03-42-18ZN.hdf'
- f = addfile(fn)
- # Read data
- vname = 'Total_Attenuated_Backscatter_532'
- var = f[vname]
- data = var[:1000,:]
- data = rot90(data)
- lats = f['latitude'][:1000,0]
- latstrs = []
- for lat in lats:
- latstrs.append('%.1f' % lat)
- height = f['metadata'][None]
- m = height.findMember('Lidar_Data_Altitudes')
- h = array(height.getArray(0, m))
- h = h[::-1]
- hstrs = []
- for hh in h:
- hstrs.append('%.1f' % hh)
- # Plot
- levs = [0.0001,0.0002,0.0003,0.0004,0.0005,0.0006,0.0007,0.0008,0.0009,\
- 0.001,0.0015,0.002,0.0025,0.003,0.0035,0.004,0.0045,0.005,0.0055,0.006,\
- 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]
- layer = imshow(data, levs, cmap='MPL_terrain')
- xticks(data.dimvalue(1), latstrs)
- yticks(data.dimvalue(0), hstrs)
- xlabel('Latitude')
- ylabel('Altitude (km)')
- colorbar(layer, label=r'$\rm{km}^{-1}$ \rm{sr}$^{-1}')
- basename = os.path.basename(fn)
- title('{0}\n{1}'.format(basename, vname))
|
|