关于三维插值及四维绘图，请指导！！！

ling19880131

VV.mat (89.58 KB, 下载次数: 1)

2015-4-13 19:25 上传

点击文件名下载附件 下载了一个关于海流V分量的nc文件，基本信息如下：
NetCDF-3 Classic E:\vcur.2012.nc {
dimensions:
time = UNLIMITED ; (12 currently)
level = 40 ;
lat = 418 ;
lon = 360 ;
variables:
   single level(level), shape = [40]
   level:long_name = "depth below sea level" ;
   level:axis = "Z" ;
single lon(lon), shape = [360]
   lon:units = "degrees_east" ;
   lon:actual_range = 1.000000 360.000000 f;
   lon:standard_name = "longitude" ;
   lon:axis = "X" ;
single lat(lat), shape = [418]
   lat:units = "degrees_north" ;
   lat:actual_range = -74.000000 64.999001 f;
   lat:standard_name = "latitude" ;
   lat:axis = "Y" ;
double time(time), shape = [12]
   time:units = "days since 1800-01-01 00:00:0.0" ;
   time:axis = "T" ;
int16 vcur(time,level,lat,lon), shape = [12 40 418 360]
   vcur:missing_value = 32767 s;
   vcur:var_desc = "v of current" ;
   vcur:statistic = "Monthly Mean" ;
    vcur:long_name = "v-component of current" ;
   vcur:units = "m/s" ;

使用interp3将vcur数据插值到0.2个精度*0.1个纬度，深度插值到5米，所用程序如下：
clc;
clear all;
close all;
% nc_dump('E:\vcur.2012.nc')
lon=nc_varget('E:\vcur.2012.nc','lon',124,11)';%取精度范围125-135E
lat=nc_varget('E:\vcur.2012.nc','lat',222,76)';%取纬度范围0-25N
time=nc_varget('E:\vcur.2012.nc','time')';
level=nc_varget('E:\vcur.2012.nc','level',0,33)';%取深度5-1500米
vcur=nc_varget('E:\vcur.2012.nc','vcur',[0,0,222,124],[1,33,76,11]);%对应取第1个月的vcur数据
save('E:\VV.mat','lon','lat','time','level','vcur');%变量全部保存在VV。mat中，见附件
xi=min(lon):0.2:max(lon);%插值到0.2个经度
yi=min(lat):0.1:max(lat);%插值到0.1个纬度
zi=min(level):5:max(level);%插值到5米
[yy,zz,xx]=meshgrid(yi,zi,xi);
V=interp3(lat,level,lon,vcur,yy,zz,xx);

到这里都没有问题，但是想画关于海流的四维海图，如下图所示，x，y，z分别表示经度、纬度和深度轴，而
用颜色表示海流的值




我使用的是surf函数，但总是报错，可能是函数使用方法不正确，报错信息如下：



请问各位，我改如何画出海流四维图呢？？？

又是那隻貓

数据的维数应该出问题了

ling19880131

又是那隻貓 发表于 2015-4-13 20:41
数据的维数应该出问题了

知道是维数的问题，但是怎么修改呢

李帅

楼主，搞定这个问题了没，告诉我一声