清空筛选条件
清空筛选条件

Regarding Interpolation( vertical)

4 次查看（过去 30 天）
Uday
Uday 2011-7-18
I have two data sets in which one is 4D data sets and another one is 2D.
1) first data: 120( longitude) x 90(latitude) x 34 (vertical level) x 8 (Time) 2) second data: 15( vertical level) x 469( number of observations)
I would like to interpolate 34 vertical levels to 15 vertical levels. so somebody knows how to fix this problem?

请先登录，再回答此问题。

回答（2 个）

Sean de Wolski
Sean de Wolski 2011-7-18
A = rand(5,5,34,5);
[ii jj kknew LL] = ndgrid(1:5,1:5,linspace(1,34,15),1:5);
C = interpn(A,ii,jj,kknew,LL);
You could just use interpn directly, thought this will use information from other dimensions.
  4 个评论
显示 2更早的评论
Uday
Uday 2011-7-19
Actually I have two different data-sets 1) 4D data (120(longitude) X 90(latitude) x 34( pressure levels) x 8 (time)), 2) 2D 15( number of vertical levels) x 469 ( number of observations those varies always). My main objective to interpolate vertical levels of two data-sets ( which is 34 & 15). so for this first of all I need to get 15 levels from 4D data. I tried with 4D data-sets the method which you have mentioned earlier, but I got error "Undefined function or method 'colon' for input arguments of type 'double' and attributes 'full 4d real".
Sean de Wolski
Sean de Wolski 2011-7-19
Show us _exactly_ what you wrote.

请先登录，再进行评论。

Andrei Bobrov
Andrei Bobrov 2011-7-18
In = DataInput; % Array 4D
[a b oldc d] = size(In);
newc = 15;
Out = zeros([a b newc d]);
for j4 = 1:d
d1 = reshape(permute(In(:,:,:,j4),[3 2 1]),oldc,[]);
Out(:,:,:,j4) = permute(reshape(cell2mat(arrayfun(@(i1)interp1((1:oldc)',d1(:,i1),linspace(1,oldc,newc)'),1:a*b,'un',0)),newc,b,a),[3 2 1]);
end

类别

在 Help Center 和 File Exchange 中查找有关 Resizing and Reshaping Matrices 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by