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/ ）。

解压后文件目录如下：

其中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文件 
   
    计算完成后，在nctoolbox目录下会生成以grb源文件名称相同的excel文件，同时将数据可视化，并在控制台输出grb文件相关信息： 
   
    补充：nctoolbox目录下demo文件夹中提供了很多案例，可使用runtests.m文件演示所有案例，根据需要对案例进行修改。

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

猜你喜欢