MATLAB Answers


How to interpolate over long section of bad data?

Katie Colfer 님이 질문을 제출함. 16 Oct 2019 20:20
최근 활동 Katie Colfer 님이 댓글을 추가함. 18 Oct 2019 12:33
Hello all,
I'm very new to MATLAB (and coding in general) so forgive me if I don't explain the problem well. Basically I have a large dataset from an experiment during which the analyzer equipment malfunctioned. I want to interpolate over a long section of data for 2 variables (VO2 and VCO2). Here's a sample of the data (VO2; starting with 3 good data points, 3 empty data points when the equipment malfunctioned, and then 3 data points following the empty cells that I don't believe are accurate):
Screen Shot 2019-10-16 at 4.16.07 PM.png
I've tried this code so far:
x = 1:4076;
yVO2 = VO2(1:4076);
yVCO2 = VCO2(1:4076);
xq = x;
splineVO2 = spline(x,yVO2);
splineVCO2 = spline(x,yVCO2);
interpVO2 = ppval(splineVO2,xq);
interpVCO2 = ppval(splineVCO2,xq);
And the returned interpolated values match the original values. I've tried this same code with only cells I want to interpolate over (2407:2822) and I get the same result. I've also tried this same code but with xq as an input argument in the spline function ( splineVO2 = spline(x,yVO2,xq); ) and have gotten the same result.
I'm guessing I'm doing something wrong with defining xq but I'm not sure.
Thanks in advance!

  댓글 수: 0

로그인 to comment.

답변 수: 1

John D'Errico 님의 답변 16 Oct 2019 23:25

I had to laugh. No, your problem is not how you defined xq, but more in how you defined x! At least you were thinking in sort of the right direction.
The point is, you don't want to define the set of all values as x and y for the spline. Only use the ones where you have valid data for x and y.
x = 1:4076;
yVO2 = VO2(1:4076);
gooddatalocs = (yVo2 ~= 0); % you don't really need the parens here.
baddatalocs = ~gooddatalocs;
% use only the good data to create an interpolant
splineVO2 = spline(x(gooddatalocs),yVO2(gooddatalocs));
% now interpolate, but at the bad points
yVO2(baddatalocs) = ppval(splineVO2,x(baddatalocs));

  댓글 수: 3

Thank you for this! This is super helpful and makes a lot of sense. However, my 'bad' data points over which I want to interpolate aren't just the 0 points, it's a string of data beginning at my 0 points. Basically, I want to interpolate 2407:2822. I've been trying to modify your code to define gooddatalocs as basically yVO2(1:2406,1) and yVO2(2823:4076) with no luck. Any thoughts?
So your bad data is comprised of the set of indices
baddatalocs = 2807:2822;
Then your good data locations are everything else. HINT: you could use setdiff to find the good stuff. Or you could just recognize the good data is [1:2406,2823:4076] and build it directly.
Got it. I figured it out. Thanks so much!

로그인 to comment.

Translated by