MeteoInfoLab脚本示例:CloudSAT Swath HDF数据
读取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)')
转载于:https://www.cnblogs.com/yaqiang/p/4881096.html
MeteoInfoLab脚本示例:CloudSAT Swath HDF数据相关推荐
- MeteoInfoLab脚本示例:AIRS Swath HDF数据
例子中的AIRS Swath HDF数据在Polar Stereographic(南极)投影中接近矩形,需要先从数据中读出经纬度及相关数据数组,利用surfacem函数绘制Swath数据(散点),在s ...
- MeteoInfoLab脚本示例:OMI Swath HDF数据
这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_ ...
- MeteoInfoLab脚本示例:TRMM 3B43 HDF数据
TRMM 3B43是卫星观测月平均降水量产品,是HDF的格点数据.需要注意的是数据中降水变量维的顺序里经度维在前纬度维在后,这与通常的设置(纬度维在前经度维在后)相反,需要对获取的二维数组进行转置,使 ...
- MeteoInfoLab脚本示例:MODIS Sinusoidal投影HDF数据
MODIS卫星很多陆面数据都是Sinusoidal投影,数据被分为一个个10*10度(赤道地区)的瓦片(http://modis-land.gsfc.nasa.gov/MODLAND_grid.htm ...
- MeteoInfoLab脚本示例:OMI Grid HDF数据
OMI卫星格点数据的例子,全球臭氧柱总量分布. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_L3-OMTO3e_2005m ...
- MeteoInfoLab脚本示例:获取气团轨迹每个节点的气象数据
读取HYSPLIT输出的轨迹数据文件和相应时间的气象数据文件,生成轨迹图层,循环每条轨迹的节点,读出该节点的经度.纬度.气压.时间,通过对气象数据插值获得该节点的气象数据. 脚本程序: #------ ...
- MeteoInfoLab脚本示例:多Y轴图
数据范围相差比较大的数据序列进行对比的时候多Y轴图就很重要了.MeteoInfoLab中提供了一个twinx函数来根据已有的坐标系(Axes)生成一个新的Axes,这个命令会使得已有的Axes不绘制右 ...
- MeteoInfoLab脚本示例:站点填图
打开包含站点填图的站点数据文件(比如micaps 1)之后,用文件对象的smodeldata函数获取StationModel数据对象,然后用stationmodel函数绘制站点填图图层. 脚本程序: ...
- MeteoInfoLab脚本示例:Streamline流线图
绘制Stramline流线图的函数是streamline,需要两个变量(U/V分量或者风向/风速). 脚本程序: f = addfile('D:/Temp/GrADS/model.ctl') u = ...
最新文章
- invalid use of incomplete type ‘class B‘
- IDC:2015年第四季度全球WLAN市场增幅最高
- python绘制3d图-python3利用Axes3D库画3D模型图
- 高退出低留存:六年百万数据透析,想颠覆传统教育的MOOC怎么了?
- [BLE--Link Layer]物理信道
- [ActionScript 3.0] 对数组中的元素进行排序Array.sort()的方法
- Web框架之Django_01初识(三大主流web框架、Django安装、Django项目创建方式及其相关配置、Django基础三件套:HttpResponse、render、redirect)
- VHDL基本门电路和数值比较器的设计
- 软件工程师 VS 硬件工程师
- pgsql数据库默认配置事务类型_postgreSql最佳配置详解(connection 申请、回收策略)...
- object-c中对文件和url操作
- DevExpress v15.1:CodeRush控件功能增强
- 给P40让路!华为Mate 30 5G降至这个价,还贵吗?
- 制造业数字化转型的启明星——低代码开发平台
- 教大家一种迅雷下载百度云文件方法
- 阿酷三合一版_阿酷插件下载|3dmax阿酷插件下载 v3.2 开源版 - 比克尔下载
- 2021苏州大学计算机考研分数,2021苏州大学考研分数线已公布
- 基于网络安全的Docker逃逸
- Codeforces_(A).Array
- 陈省身文集51——闭黎曼流形高斯-博内公式的一个简单的内蕴证明
热门文章
- java web 获取 mac 地址_java web根据访问的ip地址获取MAC地址
- win8仅仅是个试验品
- 织梦 dedecms 大批量文章sitemap生成方法
- linux firefox提示“firefox is already running”的解决方法
- Echarts笔记:全国主要城市空气质量地图
- (转载)typed用法总结
- Springboot入门之集成MybatisPlus
- 手机上如何修复老照片?老照片高清修复的方法
- 安装tpcc mysql_MySQL:安装tpcc 标准测试包
- 软考高级信息系统项目管理师——单代号网络图,双代号网络图,双代号时标网络图,CPM关键路径法计算