我想在一个循环中做以下事情,这样效率会更高:
data_Clim1 = griddata(lonClim, latClim, ice_Clim(:,:,1), X, Y);
data_Clim2 = griddata(lonClim, latClim, ice_Clim(:,:,2), X, Y);
data_Clim3 = griddata(lonClim, latClim, ice_Clim(:,:,3), X, Y);
data_Clim4 = griddata(lonClim, latClim, ice_Clim(:,:,4), X, Y);
data_Clim5 = griddata(lonClim, latClim, ice_Clim(:,:,5), X, Y);
data_Clim6 = griddata(lonClim, latClim, ice_Clim(:,:,6), X, Y);
data_Clim7 = griddata(lonClim, latClim, ice_Clim(:,:,7), X, Y);
data_Clim8 = griddata(lonClim, latClim, ice_Clim(:,:,8), X, Y);
data_Clim9 = griddata(lonClim, latClim, ice_Clim(:,:,9), X, Y);
data_Clim10 = griddata(lonClim, latClim, ice_Clim(:,:,10), X, Y);
data_Clim11 = griddata(lonClim, latClim, ice_Clim(:,:,11), X, Y);
data_Clim12 = griddata(lonClim, latClim, ice_Clim(:,:,12), X, Y);
% Add together matrices
data_Clim = [data_Clim1; data_Clim2; data_Clim3; data_Clim4; ...
data_Clim5; data_Clim6; data_Clim7; data_Clim8; data_Clim9; ...
data_Clim10; data_Clim11; data_Clim12];
% Reshape to 3D matrix
data_Clim = reshape(data_Clim, 1437, 159, 12); 基本上,我有一个360x161x12的矩阵,我需要将它插值到1437x159x12。然而,现在,使用网格数据,我一次只能插值z维中的一个。然后我组合所有独立的矩阵(有12个,因为它是每个月一个),并将其重塑为我需要的维度。
我假设while循环或for循环会不那么麻烦。
发布于 2013-06-25 01:43:17
尝试如下所示:
dataClim = zeros(1437,159,12); // pre-allocation for speed
for i=1:12
// Reshape each sub-matrix
dataClim(:,:,i) = reshape(griddata(lonClim, latClim, ice_Clim(:,:,i), X, Y),1437,159);
endhttps://stackoverflow.com/questions/17281433
复制相似问题