Interpolating missing values

Question

Jenny 2011-11-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20185-interpolating-missing-values

回答： Jennifer Rebbin 2023-9-20

I have a matrix with NaNs in some of the cells. I need to find all cells of the matrix with NaN values, then interpolate the missing values using the first value before and after the NaNs, and put the interpolated values in the place of the NaN values within the matrix. NaNs may be only one at a time or several in a row, the first and last value around NaN sections may span more than one row.

I have looked at find and interp1 in Matlab help, but can't make it work for my situation. For example: how do you make find locate NaN values? And with interp1, I am not sure what all the variables should be in this case. Any tips?

2 个评论
显示无隐藏无

Walter Roberson 2011-11-3

find(isnan(X))

Jenny 2011-11-6

thanks!

请先登录，再进行评论。

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

Answer 1

Jennifer Rebbin 2023-9-20

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20185-interpolating-missing-values#answer_1314647

Starting in R2023b, you can fill missing entries in 2-D data using the fillmissing2 function.

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

请先登录，再进行评论。

Answer 2

Fangjun Jiang 2011-11-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20185-interpolating-missing-values#answer_26662

在 MATLAB Online 中打开

Interpolation on a 1-D vector is not that hard. Maybe it could be expanded to 2-D matrix too. See if interp1() is all you need.

%%original data
x=1:0.1:10;
y=sin(x);
ind=round(ceil(90*rand(10,1)));
y(ind)=nan;
figure(1);plot(x,y);
%processing
ind=~isnan(y);
NewY=interp1(x(ind),y(ind),x);
figure(2);plot(x,NewY);

For 2-D

%%original good data
[x,y,z]=peaks;
figure(1); surf(x,y,z);
%missing some z values
[m,n]=size(z);
ind_m=ceil((m-1)*rand(10,1));
ind_n=ceil((n-1)*rand(10,1));
Incomplete_z=z;
Incomplete_z(ind_m,ind_n)=nan;
figure(2);surf(x,y,Incomplete_z);
%interpolate to recover z
index=isnan(Incomplete_z);
row_index=any(index,2);
col_index=any(index,1);
NewX=x(~row_index,~col_index);
NewY=y(~row_index,~col_index);
NewZ=Incomplete_z(~row_index,~col_index);
Interpolate_z=interp2(NewX,NewY,NewZ,x,y);
figure(3);surf(x,y,Interpolate_z)

2 个评论
显示无隐藏无

Jenny 2011-11-6

I understand what you've done here, but I am not sure how to make it fit if I have a m by n matrix with randomly dispersed nan values in it. I don't know what to set my x,Y, and xi to in the equation yi = interp1(x,Y,xi). I can use find=(isnan(X')) to find all the indeces of the nan's by row, but I don't know where that fits in the equation above.