matlab读取气象数据,基于MATLAB实现3种气象数据的读取和绘图

第 29卷第 4期 2007年 11月 广 东 气 象 GuangdongMeteorology Vol. 29　　No. 4 November　2007 基于 MATLAB实现 3种气象数据的读取和绘图 刘 　峰 , 刘 　娟 , 李 　萍 , 林 　智 (民航中南空中交通管理局气象中心 ,广东广州 　510406) 　　摘 　要 :针对 M ICAPS、 GR IB和 NetCDF三种气象数据的特点 ,利用 MATLAB软件强大的数值计算和图形制作功能 ,介绍实现数据读取以及在地图上绘制等值线、填色等值线、风矢量和流线 4种图形的功能和方法。 关键词 :计算机应用 ;MATLAB绘图 ;气象数据 ;背景地图 中图分类号 : TP317　　　　文献标识码 :B　　　文章编号 : 1007 - 6190 (2007) 04 - 0045 - 03 收稿日期: 2007 - 08 - 09 作者简介:刘 　峰 (1972年生 ) ,男 ,工程师 ,硕士 ,主要从事航空气象预报和研究。 　　随着气象业务的不断发展 ,庞大的气象数据要求标准和高效的数据存储和交换格式。M ICAPS、GR IB 和NetCDF等已经成为通用的数据格式。使用这些数据前 , 需要借助专门软件或者自行编制程序完成数据读取和处理 [1 - 4 ] ,这对使用者提出了较高的要求。可以考虑使用具备数值计算和可视化能力的 MATLAB[5 ]软件实现数据提取以及可视化输出。 利用 MATLAB软件强大的数值计算和图形制作功能可实现对 M ICAPS、GR IB 和 NetCDF三种气象数据的读取 ,以及在地图上绘制等值线、填色等值线、风矢量和流线 4种图形的功能。本文介绍利用 MATLAB在 W indows 环境下 ,如何实现 M ICAPS、GR IB和 NetCDF三种气象数据的读取 ,以及物理量图绘制的方法 ,为 MATLAB在气象领域的应用提供一些参考。 1　M ICAPS数据 M ICAPS的数据文件为文本文件 [6 ]。M ICAPS数据有 2种 ,一种是格点值 ;另一种是离散点资料。由于 MAT2LAB 具备对离散数据插值成格点数据的功能 ,因此下面介绍对第 2种资料的处理。 选用的数据是第 2类数据格式 ,即高空全要素填图。数据包括 2部分 : 第 1部分是文件头 ,第 2部分是数据 ,包括 :区站号、经度、纬度、拔海高度、级别、高度、温度、温度露点差、风向、风速。 数据的预处理分 2步。首先去除文件头 ,即将第 1和第 2 行删去。然后将缺测数据删去 , 即将要素值是 “9999”的行删去。经过预处理的数据是 10列实数。 1. 1　读取资料和插值 使用 MATLAB的函数 fopen()和 fscanf()读取预处理后的 M ICAPS数据。代码为 : fid = fopen(‘06050608. 000’ , ‘r’) ; data_all = fscanf(fid,‘% g % g % g % g % g % g % g % g % g % g’ , [10 inf]) ; %读取 10列数据lon_1 = data_all(2, : ) ; 　%第 2列是经度lat_1 = data_all(3, : ) ; 　%第 3列是纬度height_1 = data_all(6, : ) ; 　%第 6列是位势高度 wd = data_all(9, : ) ; 　%第 9列是风向ws = data_all(10, : ) ; 　%第 10列是风速使用插值函数 gr