MeteoInfo 发表于 2011-7-13 13:00
micaps第2类数据是高空全要素填图，应该是站点格式数据。micaps第11类格式是专门的格点矢量数据。不明白楼主 ...
谢谢MeteoInfo的指点
从Micaps当初定义的数据格式来看，第二类要素的确是高空全要素填图，也当然是站点格式的数据，因为当初定义数据格式是这样的：
文件头： diamond 2 数据说明（字符串） 年 月 日 时次 层次 总站点数（均为整数） 注：此类数据用于规范的高空填图 数据： 区站号（长整数） 经度 纬度 拔海高度（均为浮点数） 站点级别（整数） 高度 温度 温度露点差 风向 风速（均为浮点数） =================
所以我这样做是不通用的，当初写这个是因为群里有人就是想转换的EC的850hpa的风场预报，我也很奇怪为什么要把这个数据存为diamond2格式，因为这个EC的数据其实是2.5*2.5的格点数据，下面是数据片段
diamond 2 11070720_024时效EC_850hPa风场
11 7 7 12 850 5328
1 0.0 90.0 1 1 9999 9999 9999 74 11
2 2.5 90.0 2 1 9999 9999 9999 76 11
3 5.0 90.0 3 1 9999 9999 9999 79 11
4 7.5 90.0 4 1 9999 9999 9999 81 11
5 10.0 90.0 5 1 9999 9999 9999 84 11
6 12.5 90.0 6 1 9999 9999 9999 86 11
7 15.0 90.0 7 1 9999 9999 9999 89 11
8 17.5 90.0 8 1 9999 9999 9999 91 11
9 20.0 90.0 9 1 9999 9999 9999 94 11
10 22.5 90.0 10 1 9999 9999 9999 96 11
11 25.0 90.0 11 1 9999 9999 9999 99 11
12 27.5 90.0 12 1 9999 9999 9999 101 11
13 30.0 90.0 13 1 9999 9999 9999 104 11
14 32.5 90.0 14 1 9999 9999 9999 106 11
15 35.0 90.0 15 1 9999 9999 9999 109 11
16 37.5 90.0 16 1 9999 9999 9999 111 11
17 40.0 90.0 17 1 9999 9999 9999 114 11
================================
可以看出，所谓的站号全部被赋值为类似格点ID（lon:144 lat:37），数据是从纬度最高处，经度最小处开始按照先经度后纬度的顺序开始排列，因此直接当做格点数据就可以写入grads进行画图，当然，如果是数据格式中定义的那种当然不能这么读取啦。
所以这样读取只适合类似EC850风场预报的这种文件。
我在其他地方也看过类似的讨论，说也可以按站点读入，然后再grads中进行插值，我还没有验证两者有什么区别。
我上传了一个示例数据，射月楼主有兴趣的话可以试一下，Meteoinfo好像还没能直接读取diamond2格式，不知道我理解的对不对呢
11070720.024
(291.45 KB, 下载次数: 44)
| 
[复制链接]
窥视卡
雷达卡
楼主
显身卡
发表于 2011-7-13 14:55:58
顶一个，收藏啦！