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

第 29卷第 4 期 广 东 气 象 Vol. 29　　No. 4 2007年 11月 Guangdong M eteorology Novem ber　2007 基于 MATLAB 实现 3种气象数据的读取和绘图 刘 　峰 , 刘 　娟 , 李 　萍 , 林 　智 (民航中南空中交通管理局气象中心 ,广东广州 　510406) 　　摘 　要 :针对 M ICA PS、GR IB 和 N etCD F三种气象数据的特点 ,利用 MA TLAB 软件强大的数值计算 和图形制作功能 ,介绍实现数据读取以及在地图上绘制等值线 、填色等值线 、风矢量和流线 4 种图形的 功能和方法 。 关键词 :计算机应用 ;MA TLAB 绘图;气象数据 ;背景地图 中图分类号 : TP3 17　　　　文献标识码 : B 　　　文章编号 : 1007 - 6 190 (2007) 04 - 0045 - 03 　　随着气象业务的不断发展 ,庞大的气象数据要求标 data_all = fscanf ( fid, ‘% g % g % g % g % g % g % g % 准和高效 的数据存储和交换格式 。M ICA PS、GR IB 和 g % g % g ’, [ 10 inf ] ) ; %读取 10 列数据 N etCD F等已经成为通用的数据格式 。使用这些数据前 , lon_ 1 = data_all (2 , : ) ; 　%第 2 列是经度 需要借助专门软件或者 自行编制程序完成数据读取和处 lat_ 1 = data_all ( 3, : ) ; 　%第 3列是纬度 理 [ 1 - 4 ] ,这对使用者提出了较高的要求 。可以考虑使用具 heigh t_ 1 = data_all ( 6, : ) ; 　%第 6列是位势高度 备数值计算和可视化能力的 MA TLAB [ 5 ] 软件实现数据提 wd = data_all ( 9, : ) ; 　%第 9列是风向 取以及可视化输出。 w s = data_all ( 10, : ) ; 　%第 10 列是风速 利用 MA TLAB 软件强大的数值计算和图形制作功能 使用插值函数 griddata ( )将离散点值插值为格点值 。 可实现对 M ICA PS、GR IB 和 N etCD F 三种气象数据 的读 代码为 : 取 , 以及在地图上绘制等值线 、填色等值线 、风矢量和流 [ lon, lat ] = m eshgrid ( [ 90: 1: 140 ] , [ 15: 1: 50 ] ) ; % 线 4种图形的功能 。本文介绍利用 MA TLAB 在 W indow s 确定格点经纬度范围 ( ) 环境下 ,如何实现 M ICA PS、GR IB 和 N etCD F 三种气象数 z1 = griddata lon_ 1, lat_ 1, height_ 1, lon, lat, ‘cub ic ’; 据的读取 , 以及物理量图绘制的方法 ,为 MA TLAB 在气象 1.