1.格点数据有时需要插值到另一个范围和分辨率的网格上,可以用 linint2 函数来插值。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=91352
2.格点插值到站点
格点数据插值到站点主要有两种方法:双线性插值和最近距离,算法都很简单,MeteoInfoLab中插值到站点有几种方法:(a)利用DimDataFile的tostation方法,(b)利用DimArray的tostation方法,(c)利用interp2d插值函数。推荐使用interp2d方法,该方法中的kind参数缺省为’linear’双线性插值,也可以设置为kind='neareast’最近距离插值(其实就是找离站点最近的格点将其值赋给站点)。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=54054
3.站点数据读取插值
用MeteoInfoLab脚本读取站点数据并匹配经纬度,绘制散点图,插值为格点数据并绘制等值线图的过程。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=43091
4.多站点数据平均
日常用到的站点数据(比如micaps 1类地面全要素观测)每个时次包含的站点都不一样,要同时计算多个站点数据的平均值就比较困难。MeteoInfoLab提供了Series类可以对一维数组设置index,方便此类操作。这里给出一个计算micaps 1类数据多个站点日平均温度的例子供参考。需要首先确定计算哪些站点(在文件cma_stations_base.csv中),然后在提取micaps文件中站点和温度数组,组合成为Series数据类型,再从中提取出设定站点数据(如果某个设定站点在micaps数据中不存在则为nan),放入一个数据列表中,最后对数据列表进行平均即可。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=51290