Error in Interp1 function when searching for a corresponding matrix

Question

Yun Inn 2012-7-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/44365-error-in-interp1-function-when-searching-for-a-corresponding-matrix

I have eight 1024X1024 matrices corresponding to 8 thicknesses (0,5,10,15,20,25,30,35cm). My code below aims to retrieve a correct matrix for a given thickness, e.g. Th12...

I am getting this error:

??? Error using ==> interp1 at 259 The values of X should be distinct.

Error in ==> Predict_EPIDM at 122 CorrFactor = interp1(col,A,Th);

What is wrong with the code? __________________________________________________________________________

% There are 8 PT 0,5,10,15,20,25,30,35cm
PT=PhysicalThickness; 
% A is matrices 1024X1024 for the 8PTs. Size(A) returns 1024X1024X8. 
A = cat(3, A{:}); 
% define the column number nearest to the input Th
Th=12;
colindex=min(find(PT>=Th));
% return the 'in between' A matrix
col=A(:,colindex)+((A(:,colindex)-A(:,colindex-1))*(Th-PT(colindex))/(PT(colindex)-PT(colindex-1)));
% Using the RDepid value to find MUx
CorrFactor = interp1(col,A,Th);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Andrei Bobrov 2012-7-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/44365-error-in-interp1-function-when-searching-for-a-corresponding-matrix#answer_54389

编辑：Andrei Bobrov 2012-7-24

在 MATLAB Online 中打开

Let A - cell array your matrices.

A = arrayfun(@(ii)randi(10,5),1:8,'un',0);
Da = cellfun(@(x)x(:),A,'un',0);
Dn = cat(2,Da{:});
Th = [0,5,10,15,20,25,30,35];
PP = arrayfun(@(x)interp1(Th,Dn(x,:),'linear','pp'),(1:size(Dn,1))','un',0);
outfun = @(th)reshape(cellfun(@(x)ppval(x,th),PP),size(A{1}));
outfun(12);

or

t = 12;
idx = find(Th > t,1,'first')-[1 0];
out = diff(cat(3,A{idx}),[],3)/diff(Th(idx))*(t - Th(idx(1)))+A{idx(1)};

2 个评论
显示无隐藏无

Yun Inn 2012-7-24

Thank you very much! I used the second option and it works!

Yun Inn 2012-8-22

编辑：Yun Inn 2012-8-22

Andrei,

I chose option 2 of the solution that u gave. If Th remains as [0,5,10,15,20,25,30,35], but 't' is now a matrix, how do I make each pixel search for the right value? Also, some pixels in 't' matrix contain value smaller than 0 and bigger than 35.

请先登录，再进行评论。