- 积分
- 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');
|
|