基于MATLAB实现3种气象数据读取和绘图.pdf

第29卷第4期 广 东 气 象 v01．29No．4 2007年11月 GuangdongMeteorology Novenlber2007 基于MATLAB实现3种气象数据的读取和绘图 刘峰。刘娟．李萍．林智 (民航中南空中交通管理局气象中心．广东广州S10406) 摘要：针对MICAPS、GRIB和NetCDF三种气象数据的特点。利用MATLAB软件强大的数值计算 和图形制作功能，介绍实现数据读取以及在地图上绘制等值线、填色等值线、风矢量和流线4种图形的 功能和方法。 关键词：计算机应用；MATLAB绘图；气象数据；背景地图 中图分类号：TF317 文献标识码：B 文章编号：1007—6190(2007J04—0045—03 随着气象业务的不断发展，庞大的气象数据要求标 准和高效的数据存储和交换格式。MICAPS、GRIB和g％g％g’[10iIlf])；％读取10列数据 lon NetCDF等已经成为通用的数据格式。使用这些数据前． l=data_all(2，：)；％第2列是经度 需要借助专门软件或者自行编制程序完成数据读取和处 lat_l=data_all(3．：)；％第3列是纬度 all(6，：)；％第6列足位势高度 理”-．】．这对使用者提出了较高的要求。可以考虑使用具 height一1=data 备数值计算和可视化能力的MATLAB”。软件实现数据提 wd=data_all(9，：)；％第9列是风向 取以及可视化输出。 ws=data_all(10．：)；％第10列是风速 利用MATLAB软件强大的数值计算和图形制作功能 使用插值函数griddata()将离散点值插值为格点值。 代码为： 可实现对MICAPS、GRIB和NetCDF三种气象数据的读 取，以及在地图上绘制等值线、填色等值线、风矢量和漉 确定格点经纬度范围 线4种图形的功能。本文介绍利用MATLAB在Windows I，lon，lat，‘cubic’)； zl=扣ddata(Ion．1，lat_l。height 环境下，如何实现MICAP$、GRIB和NetCDF三种气象数 1．2绘制气象圈 据的读取，以及物理量图绘制的方法．为MATLAB在气象 领域的应用提供一些参考。 陈金辉”o研究了使用MATLAB绘制气象图，但没有 给出如何使用地图为底图。哥伦比亚大学开发了M—Map 1 MICAPs数据 软件包。主要功能是支持18种地图投影方式、10分辨率 MICAP$的数据文件为文本文件”’。MICAPS数据有 的全球地形数据集和1／4。分辨率的海岸线数据集。图1 2种。一种是格点值；另一种是离散点资料。由于MAT- 给出了基于软件包绘制的等值线图和风矢量图。 LAB具备对离散数据插值成格点数据的功能，因此下面 制作方法如下： 介绍对第2种资料的处理。 设置地图投影方式、范围和海岸线使用函数m一州 选用的数据是第2类数据格式，即高空全要素填图。 ()和re_coast()。代码为：