matlab读取grb、grb2数据并导出、绘图

grb数据是grib码文件, GRIB  码是与计算机无关的压缩的二进制编码,主要用来表示数值天气预报的产品资料。现行的 GRIB  码版本有 GRIB 1 和 GRIB 2 两种格式。本文讲解nctoolbox的基本使用,并读取grb文件数据输出到excel文件中,并绘制相关图件。
本文使用nctoolbox工具进行操作并扩展其功能(安装包见csdn: http://download.csdn.net/download/thesunrize/9757036;Github: https://github.com/Miaor/nctoolbox-extended/ )。  
解压后文件目录如下:
matlab读取grb、grb2数据并导出、绘图_第1张图片
其中readfile.m文件完成数据的导入、提取、输出至excel、属性打印以及绘图功能,代码如下:
function  readfile
setup_nctoolbox %初始化nctoolbox,操作前请检查电脑是否安装java1.6及以上
[filename, pathname] = uigetfile( '*.GRB' 'choose a GRB file' ); %选择grb文件,支持grb、grb2以及nc文件
if  isequal(filename,0)    
    msgbox( 'you choose nothing' );
else
    pathfile=fullfile(pathname, filename); %获取文件路径
ds= ncdataset(pathfile);%读取数据集
    ds.netcdf %display file details
 
%绘图
GPMData = ds.data(ds.variables{2});
GPMData = squeeze(GPMData);
Temp = ds.time( 'time' , ds.data( 'time' , 1));
lat = ds.data( 'lat' );
lon = ds.data( 'lon' );
  
surf(lon,lat,GPMData);
xlabel( 'lontitude' )
ylabel( 'latitude' )
zlabel( 'GPM' )
shading interp;
title({datestr(Temp), '3d' }); 
colorbar;
 
%数据导出到excel
  name = { 'lontitude' , 'latitude' , 'GPM' };
  xlswrite(strcat(filename, '.xlsx' ), name,1, 'A1' )
  xlswrite(strcat(filename, '.xlsx' ), lon,1, 'A2' )
  xlswrite(strcat(filename, '.xlsx' ), lat,1, 'B2' )
  xlswrite(strcat(filename, '.xlsx' ), GPMData,1, 'C2' )
      
msgbox( 'finished!' );
end



将readfile.m文件拖入控制台即可运行,效果如下:

选择需要导入的grb文件
matlab读取grb、grb2数据并导出、绘图_第2张图片


计算完成后,在nctoolbox目录下会生成以grb源文件名称相同的excel文件,同时将数据可视化,并在控制台输出grb文件相关信息:

matlab读取grb、grb2数据并导出、绘图_第3张图片
matlab读取grb、grb2数据并导出、绘图_第4张图片


补充:nctoolbox目录下demo文件夹中提供了很多案例,可使用runtests.m文件演示所有案例,根据需要对案例进行修改。
matlab读取grb、grb2数据并导出、绘图_第5张图片




你可能感兴趣的:(MATLAB)