MeteoInfoLab脚本示例:TRMM 2A12 HDF数据
2024-08-30 06:45:06
TRMM 2A12 HDF数据是卫星观测的SWATH数据(轨道数据),比格点数据处理起来要麻烦一些。数据的经纬度保存在geolocation变量中,需要先将经纬度数据读出来(均为2维数组),然后读取云水含量数据(cldWater)。虽然都是2维数组,但并不是格点数据,相当于2维的散点数据,点数相当多(3019*208)如果用散点图来绘制的话会非常的慢,需要将其插值为格点数据。插值方法建议选择nearest,该方法速度最快。插值后的格点数据用imshowm函数显示为图像就很快了。
脚本程序:
#Add data file
folder = 'D:/Temp/hdf/'
fns = '2A12_090711_66394_6.HDF'
fn = folder + fns
f = addfile(fn)
#Get data variable
geo = f['geolocation']
lat = geo[:,:,0]
lon = geo[:,:,1]
water = f['cldWater'][:,:,9]
water = water / 1000
#Interpolate to grid data
gwater, x_g, y_g = griddata((lon, lat), water, method='neareast')
#Plot
axesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
levs = arange(0, 10, 0.5)
layer = imshowm(x_g, y_g, gwater, levs)
colorbar(layer, orientation='horizontal')
title(fns + ' (cldWater level=9)')
axism([-180, 180, -90, 90])
最新文章
- ASP.NET MVC5+EF6+EasyUI 后台管理系统(40)-精准在线人数统计实现-【过滤器+Cache】
- RabbitMQ 通过记日志来看routingkey
- ng-init
- 常用的MIME类型(资源的媒体类型)
- pyqt显示指定范围的数字
- uva 12100 Printer Queue
- 转载【ViewPager+Fragment】ViewPager中切换界面Fragment被销毁的问题分析
- 对AppStore中的项目进行评分(转载)
- ActiveMQ 503错误
- ActiveMQ_Windows版本的安装部署
- BZOJ_3969_[WF2013]Low Power_二分答案
- Centos7创建支持ssh服务的docker镜像
- GitHub如何删除一个仓库
- ES - es为什么要移除type?
- Arukas.io云主机安装CentOS
- LeetCode 33 Search in Rotated Sorted Array(循环有序数组中进行查找操作)
- 如何Python写一个安卓APP
- CSS中的各种width(宽度)
- JSTL 、 OGNL 与 EL
- CMMI5级——原因分析及解决方案(Causal Analysis and Resolution)