打开包含站点填图的站点数据文件(比如micaps 1)之后,用文件对象的smodeldata函数获取StationModel数据对象,然后用stationmodel函数绘制站点填图图层.脚本程序: #Set data folders basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo' datadir = os.path.join(basedir, 'sample/MICAPS') mapdir = os.path.join…

读取风场U/V变量数据,可以从U/V计算出风速:speed = sqrt(u*u+v*v).quiverm函数用来绘制风场矢量图,参数中包括U/V变量,如果要绘制彩色风场还需要第三个变量,这里是风速speed,10是颜色的分级数,size=8指定了风场箭头的基准大小. 脚本程序: f = addfile('D:/Temp/GrADS/model.ctl') u = f['U'][0,[500],[10,60],[60,140]] v = f['V'][0,[500],[10,60],[60,14…

数据范围相差比较大的数据序列进行对比的时候多Y轴图就很重要了.MeteoInfoLab中提供了一个twinx函数来根据已有的坐标系(Axes)生成一个新的Axes,这个命令会使得已有的Axes不绘制右边的Y轴,而新生成的Axes只绘制右边的Y轴.yaxis函数可以对某个Axes的Y轴进行相应的设置,第一个参数是某个Axes的对象,可以设置Y轴的颜色(color)和偏移(shift),对于超过2个Y轴的图形来说第3个Y轴就必须设置偏移以避免压盖.双Y轴图: 多Y轴图: 脚本程序: ax1 = ax…

最近的一篇文章中需要绘制多Y轴图形,Excel只能做双Y轴图,又尝试了Origin,SigmaPlot,Igor等软件,手动做起来相当繁琐,批量做更是觉得费劲,干脆尝试在MeteoInfoLab里实现这样的功能.目前已经初具成效:…

在脚本中使用中文需要指明是unicode编码,即在含有中文的字符串前加u,比如:u'中文'.还需要将字体指定为一种中文字体.详见下面的例子.脚本程序: x = [1,2,3,4] y = [1,4,9,16] plot(x, y, 'r--*') ylabel(u'Y轴', fontname=u'宋体', fontsize=16) xlabel(u'X轴', fontname=u'宋体', fontsize=16) axis([0,5,0,20]) title(u'示 例 图 形', fontn…

1.已有数据绘图如下,其中网络流量的单位是M Bytes/s,与另外两组数据的单位(时间)不同,现在要为其添加右侧Y轴. 2.首先选中该图像,找到工具条,点击第三个按钮“Add Right-Y Layer”,这样就为其添加了副纵坐标轴,如图: 3.然后在图形上点击右键,选中Plot Setup菜单,弹出如下的窗口: 上图中可以看出,新添加的RightY下无任何数据,下面将要放置到右侧Y轴上绘制的数据添加到该坐标轴下,首先选择所在的worsheet,然后指定X\Y数据列,选中RightY1行,再点…

有时候我们会遇到一种图表,就是X轴一致,可是Y轴的数据相差很大.如下图中,年龄和收入就不是一个数量级,在图表中显示的时候,“年龄”的曲线根本看不到(表中数据仅供举例): 解决的方法就是使用双Y轴显示,把年龄的数据可以显示到另一边.步骤如下: 1.在图表工具条中选中数据值较小的那个系列.本例中的"年龄". 2.在图中的系列(下图中的四个小黑方框上面)上,单击右键.选择“数据系列格式”. 3.选择“数据系列格式”后,在如下图页面,选择“次坐标轴”.o(∩_∩)o…确定,看到结果了吗? 4.…

AIRS Grid HDF数据是HDF4 EOS格式,数据地理坐标信息可以被MeteoInfo自动识别,脚本程序更为简单.需要注意的是读取数据时Y轴是反向的(卫星数据通常如此).脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'AIRS.2002.08.01.L3.RetStd_H031.v4.0.21.0.G06104133732.hdf' fn = folder + fns f = addfile(fn) vname = 'Tempera…

FY3A卫星有HDF格式的AOD产品数据,全球范围,分辨率为0.05度.读取数据文件变量后要重新设定X/Y维,数据是Y轴反向的,且需要除以10000得到AOD值. 脚本程序: #Add data file fn = 'D:/Temp/hdf/FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140331_AOAM_5000M_MS.HDF' f = addfile(fn) #Get data variable v = f['Aerosol_Optical_Thickness_of…

尝试编写MeteoInfoLab脚本计算垂直螺旋度,结果未经验证. 脚本程序: print 'Open data files...' f_uwnd = addfile('D:/Temp/nc/uwnd.2011.nc') f_vwnd = addfile('D:/Temp/nc/vwnd.2011.nc') f_omega = addfile('D:/Temp/nc/omega.2011.nc') print 'Calculate vertical helicity...' tidx = 173…