2013-04-22 64 views

回答

7

可用GRIB2在NCTOOLBOX for Matlab,你可以打开一个GRIB2文件就像本地的NetCDF文件或远程OPeNDAP数据集:

% download data 
! wget ftp://polar.ncep.noaa.gov/pub/history/waves/multi_1.at_4m.dp.200607.grb2 

% create ncgeodataset object 
nc=ncgeodataset('multi_1.at_4m.dp.200607.grb2'); 

% list variables 
nc.variables 

% create geovariable object 
dirvar=nc.geovariable('Primary_wave_direction_degree_true_surface'); 

% get data at 1st time step 
dir=dirvar.data(1,:,:); 

% get grid at 1st time step 
g=dirvar.grid_interop(1,:,:); 

% plot 
pcolorjw(g.lon,g.lat,dir); 
title(datestr(g.time)) 

enter image description here

+0

非常感谢你丰富。你的回答非常有帮助。我还想知道如何在给定坐标处绘制时间序列波形数据(即波高)?目前WWIII模型给出了每个月的波浪数据,如何将每月的时间序列波浪图组合成一个全年的时间序列?例如,我想在一个给定的位置绘制一个整年的时间序列的波高,例如1998年。 – user2309188 2013-04-23 22:27:38

+0

非常感谢您能否提供帮助。问候,伊恩 – user2309188 2013-04-23 22:28:18

+0

伊恩,你应该问这是一个单独的问题 – 2013-04-24 11:24:36