2013-04-22 406 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