- 积分
 - 286
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2014-4-14
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
 
 楼主 |
发表于 2014-7-28 21:42:21
|
显示全部楼层
 
 
 
 
画的是时间-高度剖面图,画出来锯齿感太强,我刚学matleb,用interp1插值了半天一直不对,麻烦帮我看下应该怎么写插值程序。 
我写的程序是这样的:ncinf=netcdf.open('f:\uwnd.mon.mean.nc','NOWRITE'); 
lon=ncread('f:\uwnd.mon.mean.nc','lon'); 
lat=ncread('f:\uwnd.mon.mean.nc','lat'); 
lev=ncread('f:\uwnd.mon.mean.nc','level'); 
time=ncread('f:\uwnd.mon.mean.nc','time'); 
uwnd=ncread('f:\uwnd.mon.mean.nc','uwnd',[1 1 1 1],[1 1 8 5],[1 1 1 1]); 
a=permute(uwnd,[3 4 1 2]); 
[time,level]=meshgrid(double(time(1:5)),double(lev(1:8))); 
contourf(time,level,a); 
set(gca, 'YDir', 'reverse'); 
 |   
 
 
 
 |