分享一个meteoinfo文本格式站点数据格式的批处理程序

平流层的萝卜

     最近有了任务，要用micaps数据画若干时次的24h降水量图，并叠加上850hPa高度场和风场。一开始用的软件是grads，结果grads要想画填色图或者等值线图必须得经过插值过程，但问题恰恰出在插值函数oacres上面，无论如何修改插值半径，总会有高值区被平滑掉。正巧接触到了meteoinfo这个的软件（看了王亚强老师的图文详解文本格式站点数据绘制等值线图），由于用meteoinfo仅仅做等值线需要插值，而填色图不需要插值。所以可以直接用填色表示降水量而避免了平滑的损失。（meteoinfo很给力，而且是一个中国人研发的国产软件，顿时提升了民族自豪感~~~~，支持国货！！ ）。
      下面直入正题：图文详解文本格式站点数据绘制等值线图提出，要作图，首先要把micaps第一类数据转化成所需要的csv格式（第一行文件名，第二行以后的第一~四列分别为站点号、经度、纬度和一种物理量），不过micaps第一类数据格式复杂，单个处理成csv格式还好，手动即可，但是很多个micaps文件的处理则很麻烦。所以，我自己写了一个批处理的fortran文件，用于将多个micaps第一类文件批处理成csv格式（见附件）。想必会方便一些。程序做了尽可能详尽的注释，用于方便使用者阅读并根据自己的情况修改。


      作为一个完完全全的meteoinfo的新手，还需要很多学习的东西，尤其是meteoinfo的批处理脚本用python语言怎么写，以后还要和大家多多学习。欢迎大家批评指正。

MeteoInfo

谢谢分享！不过MeteoInfo直接支持MICAPS第一类数据格式，无需再转为CSV文件。

MeteoInfo

另外，你说的填色图不需要插值不是很明白，站点数据作等值线及等值线填色图都需要先插值为格点数据。

kongfeng0824

请问站点资料使用插值到格点，如果站点太少，不知道会出什么问题吧？

平流层的萝卜

MeteoInfo 发表于 2013-3-22 23:06

                               
登录/注册后可看大图

另外，你说的填色图不需要插值不是很明白，站点数据作等值线及等值线填色图都需要先插值为格点数据。

谢谢。是我弄错了，果然都是要插值的，还需要再学习一下{:eb308:}。嗯，了解可以直接支持micaps第一类数据，看了您的帖子，就想试试其他的方法，结果是一样的。正在学习meteoinfo插值，还有很多不明白的地方要学习！

平流层的萝卜

kongfeng0824 发表于 2013-3-23 03:10

                               
登录/注册后可看大图

请问站点资料使用插值到格点，如果站点太少，不知道会出什么问题吧？

只试过grads里面，如果站点太稀疏，会产生虚大值，不知道meteoinfo会不会

張顧.煒Zane

最近开始看看meteoinfo  小马哥记得经常交流 啊

平流层的萝卜

張顧.煒Zane 发表于 2013-3-23 21:06

                               
登录/注册后可看大图

最近开始看看meteoinfo  小马哥记得经常交流 啊

ok，常交流啊

kongfeng0824

平流层的萝卜 发表于 2013-3-23 21:05

                               
登录/注册后可看大图

只试过grads里面，如果站点太稀疏，会产生虚大值，不知道meteoinfo会不会

我一般差值会用arcgis或者surfer

平流层的萝卜

kongfeng0824 发表于 2013-3-23 23:53

                               
登录/注册后可看大图

我一般差值会用arcgis或者surfer

那表示不懂