MeteoInfo

摘要：站点数据绘制等值线需要首先将站点数据插值为格点数据，MeteoInfo中提供了反距离权法（IDW）和cressman两个方法，其中IDW方法可以有插值半径的选项。这里示例读取一个MICAPS第一类数据（地面全要素观测），获取6小时累积降水数据（Precipitation6h），然后用站点数据的gri... 阅读全文

摘要：添加南海脚图实际上就是增加一个地图坐标系（Axesm），在用axesm函数创建地图坐标系的时候可以指定position参数（即位置和大小），第1和第2个参数是坐标系的左下角x, y位置，后两个参数是坐标系的宽度和高度，单位都是Figure的归一化坐标（0 - 1）。由于地图坐标系在Figure缩放时... 阅读全文

摘要：在用axesm函数创建地图坐标系的时候可以指定地图投影（设置projinfo参数），地图投影可以通过projinfo函数来创建，里面的参数依据proj4投影字符串，可以参考此网页：http://remotesensing.org/geotiff/proj_list。投影后的地图坐标系显示范围的设置可... 阅读全文

摘要：应用最广泛的的地图数据应该是shape格式，网络上有很多免费下载资源。MeteoInfoLab中读取shape文件的函数是shaperead，参数即文件名，返回数据包含图形和属性信息的图层对象。矢量图层按照图元类型通常分点、线、面三种，都可以通过geoshow函数加入到地图坐标系（Axesm）中。在... 阅读全文

摘要：这里演示从micaps第一类数据（地面全要素观测）中读取一个变量（用DimDataFile类的stationdata方法），然后maskout掉中国区域之外的数据，利用scatterm函数绘制散点图。脚本程序：f = addfile_micaps('D:/Temp/micaps/10101414.0... 阅读全文

摘要：inpollygon函数是用来判断带坐标（x/y）的数据是否在某个或者一组多边形（Polygon）中，返回的结果中如果做多边形内则值为1，否则值为-1。下面一个例子演示了利用一个shape文件和inpolygon函数生成这种0、1数据。需要下载最新的MeteoInfo版本（1.2.6R1）：http... 阅读全文

摘要：增加了对AWX卫星数据格式的支持（和C#版本的功能相当）。在MeteoInfo中打开AWX数据：在MeteoInfoLab中打开AWX数据： 阅读全文

摘要：MeteoInfoLab提供一个线性拟合函数linregress，参数是参与拟合的两个数据序列，返回拟合的斜率、截距和相关系数。有了上述拟合参数可以用polyval函数生成拟合数据（直线）。然后可以将数据、拟合线、公式等绘图。脚本程序：fn = os.path.join('D:/KeyData/PM... 阅读全文

摘要：数据范围相差比较大的数据序列进行对比的时候多Y轴图就很重要了。MeteoInfoLab中提供了一个twinx函数来根据已有的坐标系（Axes）生成一个新的Axes，这个命令会使得已有的Axes不绘制右边的Y轴，而新生成的Axes只绘制右边的Y轴。yaxis函数可以对某个Axes的Y轴进行相应的设置，... 阅读全文

摘要：LaTeX是排版常用的语法，科学计算软件中也常用它来写数学公式（比如MatLab, Matplotlib等），MeteoInfo通过调用JMathLaTeX库也可以实现这样的功能。LaTeX的语法介绍可以参考此网页：http://matplotlib.org/users/mathtext.html下... 阅读全文

摘要：subplot命令可以自动等间距分配多个坐标系（Axes），命令中有三个参数，前两个定义了行数和列数，第三个指定了当前的坐标系（Axes），绘图命令是作用在当前坐标系中的。脚本程序：def f(t): return exp(-t) * cos(2*pi*t)t1 = arange(0., 5.... 阅读全文

摘要：绘图的时候首先要有坐标系（Axes），可以用axes命令创建，如果没有创建在绘图时会自动创建一个。参数里的position是用来置顶坐标系的图形（figure）中的位置的，通过位置置顶，可以将多个坐标系组合在一个figure中。位置是以figure为单位坐标（0 - 1之间），前两个参数是左下角的x... 阅读全文

摘要：气象数据基本为多维数据（通常是4维，空间3维加时间维），只让数据中一维可变，其它维均固定即可提取一维数据。比如此例中固定了时间维、高度维、纬度维，只保留经度维可变：hgt = f['hgt'][0,[500],[4],[180,360]]。固定某一维可以用序号（从0开始），比如此例中的时间维为0，既... 阅读全文

摘要：绘制格点数据的散点图，用scaterm函数。脚本程序：f = addfile('D:/Temp/GrADS/model.ctl')ps = f['PS'][0,(10,60),(60,140)]axesm()mlayer = shaperead('D:/Temp/map/country1.shp')... 阅读全文

摘要：读取风场U/V变量数据，可以从U/V计算出风速：speed = sqrt(u*u+v*v)。quiverm函数用来绘制风场矢量图，参数中包括U/V变量，如果要绘制彩色风场还需要第三个变量，这里是风速speed，10是颜色的分级数，size=8指定了风场箭头的基准大小。脚本程序：f = addfile... 阅读全文

摘要：Maskout通常有两种类型：Maskout图形和Maskout数据。这里是Maskout图形的示例。需要用shaperead读取地图数据形成图层作为Maskout图层（这里是中国的行政区域china.shp），然后用masklayer函数对数据图层进行屏蔽，第一个参数是用来Maskout的图层，第... 阅读全文

摘要：FLEXPART是一个类似HYSPLIT的扩散模式，它输出的netcdf文件参照了WRF，可惜全局属性没有写全，比如只有一个投影名称（例如Lambert），没有相关的投影参数：中央经度，标准纬度等等。必须查阅WRF的头文件才能重建投影（为什么不照猫画虎把属性写全呢？）。数据的经纬度坐标是有的，但在L... 阅读全文

摘要：此例中的降水文本文件下载自“中国气象科学数据共享服务网”（http://cdc.nmic.cn/sksj.do?method=ssrjscp），其实是ESRI的文本格点数据格式。对于文本数据，MeteoInfoLab中有asciiread函数来读取数据，第一个参数是文件名，然后可以有一些可选参数，h... 阅读全文

摘要：MODIS的气溶胶光学厚度（AOD）产品应用很广，数据可以在Giovanni上下载：http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html。有HDF, netCDF和文本格式，这里示例打开netCDF格式数据，提取AOD数据并绘图。脚本程... 阅读全文

摘要：用ncep数据计算水汽通量散度的脚本。需要air, uwnd, vwnd和rhum变量。数据是4维数据，需要固定时间维和高度维，数据中纬度维的数据是反向的，因此读取时需要特殊的设置（::-1）。脚本中用到几个内置的函数：cdiff, hdivg和magnitude，和GrADS中同名函数的作用一样。... 阅读全文

摘要：提升了MeteoInfoLab脚本数据处理能力，比如双Y轴图、多Y轴图、数组计算、坐标投影计算等。这里给出几个示例图，以后有空了会将示例脚本程序整理放在网上。坐标投影计算：双Y轴图：多Y轴图：多Y轴图（实际数据）：数据线性拟合：水汽通量散度计算： 阅读全文

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

摘要：在JLaTeXMath库（http://forge.scilab.org/index.php/p/jlatexmath/）的支持下，实现了利用LaTeX语法显示特殊符号和数学公式的功能。需要在字符串首尾加“$”字符。在MeteoInfo桌面软件中：在MeteoInfoLab中：相应的脚本程序：def... 阅读全文