Interpolation and missing values

Question

Chad Greene 2014-2-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/117174-interpolation-and-missing-values

评论： Chad Greene 2014-3-11

I have a series of data given by an x array and a corresponding y array. The arrays are equal in size. The x array is approximately equally spaced, but it has a few small gaps, and some very large gaps. I would like to interpolate to get numeric values in the small gaps, but I'd like to replace the big gaps with NaNs.

Here is an example of some small gaps, and one big gap. I'd like to replace the obviously incorrect interpolation with NaNs. How can I do this?

% create some data:
x = 0:.02:15; 
y = sin(x); 
% create some holes in the data
x([3 25 32:33 200:280 410:415]) = []; 
y([3 25 32:33 200:280 410:415]) = []; 
plot(x,y,'ko'); hold on
% the interpolation I'd like to improve: 
xi = 0:.015:15;
yi = interp1(x,y,xi,'cubic'); 
plot(xi,yi,'b.')

3 个评论
显示 1更早的评论隐藏 1更早的评论

Chad Greene 2014-2-21

在 MATLAB Online 中打开

Here is a loop that almost solves the problem. I'd rather avoid loops, but the biggest problem with this solution is the little tails that persist:

x = 0:.02:15; 
y = sin(x); 
x([3 25 32:33 200:280 410:415]) = []; 
y([3 25 32:33 200:280 410:415]) = []; 
plot(x,y,'ko'); hold on
xi = 0:.015:15;
yi = interp1(x,y,xi,'cubic'); 
for n = 1:length(xi)
    if min(abs(xi(n)-x))>.1
        yi(n) = NaN; 
    end
end
plot(xi,yi,'b.')

Above, a threshold is set to set to change any values more than 0.1 x units away from an x data point to NaN.

Chad Greene 2014-2-21

I do have the Image Processing Toolbox, and I'm reading the documentation right now. Application of bwareaopen and regionprops is not immediately apparent.

请先登录，再进行评论。

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

Answer 1

Paul 2014-2-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/117174-interpolation-and-missing-values#answer_125329

在 MATLAB Online 中打开

x_gap = x(2:end)-x(1:end-1);
ind=find(x_gap>0.025);
ind_int=[];
for j=1:numel(ind)
    ind_int = [ind_int,find((xi>x(ind(j)) & xi<x(ind(j)+1)))];
end
yi(ind_int)=NaN;

4 个评论
显示 2更早的评论隐藏 2更早的评论

Image Analyst 2014-2-22

在 MATLAB Online 中打开

Well I couldn't get it to work but you did so that's all that counts:

clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 18;
% create some data:
x = 0:.02:15; 
y = sin(x); 
% Create some holes in the data
x([3 25 32:33 200:280 410:415]) = []; 
y([3 25 32:33 200:280 410:415]) = []; 
% Plot data
subplot(3,1,1);
plot(x, y, 'ro-');
grid on;
title('Original Signal', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
xi = 0:.015:15;
x_gap = x(2:end)-x(1:end-1);
subplot(3,1,2);
bar(x_gap);
title('x_gap', 'FontSize', fontSize, 'Interpreter', 'None');
ind=find(x_gap>0.025)
ind_int=[];
for j=1:numel(ind)
    ind_int = [ind_int,find((xi>x(ind(j)) & xi<x(ind(j)+1)))];
end
yi(ind_int)=NaN;
% Plot interpolated arrays.
subplot(3,1,3);
title('Interpolated Signal', 'FontSize', fontSize);
plot(xi, yi, 'bs-');

And I won't bother with using regionprops() since you already found a working solution.

Paul 2014-2-23

在 MATLAB Online 中打开

add this:

yi = interp1(x,y,xi,'cubic');

after xi = 0:.015:15; and then it works.

Chad Greene 2014-3-11

I turned Paul's solution into a function: http://www.mathworks.com/matlabcentral/fileexchange/45842. Thank you both for your help.

请先登录，再进行评论。

Interpolation and missing values

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Interpolation and missing values

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

3 个评论
显示 1更早的评论隐藏 1更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论