MeteoInfoLab脚本示例：AIRS Swath HDF数据

例子中的AIRS Swath HDF数据在Polar Stereographic（南极）投影中接近矩形，需要先从数据中读出经纬度及相关数据数组，利用surfacem函数绘制Swath数据（散点），在surfacem函数中经纬度数据会被投影到目的投影坐标系（在axesm函数中定义），并从2维散点数据生成surface数据（也就是格点数据）用所谓的surface插值方法。2维散点数据可以组成格网（非矩形格网），遍历目标格点数据的每个格点，该格点落在格网中哪个网格中，就将此网格的散点值赋给该格点。有些拗口，其实算法不是很复杂，和MatLab中的surfacem函数的功能类似。

脚本程序：

#Add data file
folder = 'D:/Temp/hdf/'
fns = 'AIRS.2002.12.31.001.L2.CC_H.v4.0.21.0.G06100185050.hdf'
fn = folder + fns
f = addfile(fn)
lon_v = f['Longitude']
lat_v = f['Latitude']
lon = lon_v[:,:]
lat = lat_v[:,:]
vname = 'radiances'
rad_v = f[vname]
rad = rad_v[:,:,567]
#Plot
plot,proj = axesm(proj='stere', lat_0=-90, gridline=True, griddx=30, griddy=30)
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
levs = arange(40, 90, 1)
#slayer = scatterm(lon, lat, rad, levs, edge=True)
layer = surfacem(lon, lat, rad, levs)
colorbar(layer, orientation='horizontal')
title('{0}\n {1}'.format(fns, vname))
axism()

转载于:https://www.cnblogs.com/yaqiang/p/4625449.html

MeteoInfoLab脚本示例：AIRS Swath HDF数据相关推荐

MeteoInfoLab脚本示例：OMI Swath HDF数据
这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_ ...
MeteoInfoLab脚本示例：CloudSAT Swath HDF数据
读取CloudSAT HDF Swath数据,绘图分上下两部分,上面是时间和高度维的Radar Reflectivity Factor二维图,下面是卫星轨迹图. 示例程序: # Add file f ...
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 = ...

MeteoInfoLab脚本示例：AIRS Swath HDF数据

MeteoInfoLab脚本示例：AIRS Swath HDF数据相关推荐

最新文章

热门文章