读取CloudSAT HDF Swath数据,绘图分上下两部分,上面是时间和高度维的Radar Reflectivity Factor二维图,下面是卫星轨迹图。

示例程序:

# Add file
f = addfile('D:/Temp/hdf/2010128055614_21420_CS_2B-GEOPROF_GRANULE_P_R04_E03.hdf') # Read data
vname = 'Radar_Reflectivity'
v_data = f[vname]
data = v_data[:,:]
v_height = f['Height']
height = v_height[0,:]
time = f['Profile_time'][:]
lon = f['Longitude'][:]
lat = f['Latitude'][:] # Read attributes
long_name = v_data.attrvalue('long_name')[0]
scale_factor = v_data.attrvalue('factor')[0]
valid_min = v_data.attrvalue('valid_range')[0]
valid_max = v_data.attrvalue('valid_range')[1]
units = v_data.attrvalue('units')[0]
units_h = v_height.attrvalue('units')[0] # Apply scale factor
valid_max = valid_max / scale_factor
valid_min = valid_min / scale_factor
data = data / scale_factor
data[data>valid_max] = nan
data[data<valid_min] = nan
data = transpose(data)
data = data[::-1,:] # Make a split window plot
subplot(2, 1, 1) # Contour the data
levs = arange(-38, 50, 2)
layer = imshow(time, height[::-1], data, levs)
colorbar(layer)
title('Radar Reflectivity Factor')
xlabel('Seconds since the start of the granule. (seconds)')
ylabel('Height (m)') # The 2nd plot is the trajectory
subplot(2, 1, 2)
axesm()
lworld = shaperead('D:/Temp/map/country1.shp')
geoshow(lworld, edgecolor='k')
plotm(lon, lat, '-b', linewidth=4)
#scatterm(lon, lat, lon, size=4, edge=False, facecolor='b')
scatterm(lon[0], lat[0], size=6, facecolor='r')
xlim(-180, 180)
ylim(-90, 90)
title('Trajectory of Flight Path (starting point in red)')

最新文章

  1. checkbox的readonly效果
  2. 20145235 《Java程序设计》第9周学习总结
  3. 覆盖率测试工具gcov的前端工具_LCOV
  4. Sprint第二个冲刺(第四天)
  5. 【软件多国语言】一个demo
  6. AIDL与stub
  7. centos7和windows7双系统安装
  8. Unity3d Realtime Dynamic Volume Clouds Rendering
  9. 把war包放到Tomcat安装文件夹下,不能直接訪问的解决方式
  10. gets()函数
  11. php将中文插入数据库出现乱码
  12. 《高性能Javascript》读书笔记-2
  13. 醒醒吧!互联网的真正未来不是AI,更不是VR,AR,而是区块链
  14. 【jQuery】复选框的全选、反选,推断哪些复选框被选中
  15. bzoj 1855: [Scoi2010]股票交易
  16. 20171201 - macOS High Sierra 神级 bug
  17. 前端使用node.js的http-server开启一个本地服务器
  18. Page Visibility(网页可见性) API与登录同步引导页实例页面
  19. ezmorph将一种对象转换成另外一种对象
  20. LPDMvvmKit 系列之 UITableView 的改造

热门文章

  1. Ajax跨域解决方案大全
  2. [Java数据结构]HashSet,LinkedHashSet,TreeeSet
  3. Linux:crond(crontab)定时任务
  4. python 3 for嵌套
  5. openstack核心组件——horizon Web管理界面(10)
  6. Java程序员必备基础:JDK 5-15都有哪些经典新特性
  7. springboot之启动端口指定
  8. java8--排序
  9. 垃圾回收器(GC)
  10. vscode链接git拉取项目