Cubic spline interpolation in interp1 doesn't remove NaN
10 views (last 30 days)
Show older comments
I have functional near-infrared spectroscopy timeseries data which sometimes contain NaNs and/or infinite values in some channels.
Baed on literature, spline interpolation seems to be the method of choice to replace these.
I use interp1 as follows:
Y.hbo(isnan(Y.hbo)) = interp1(find(~isnan(Y.hbo)), Y.hbo(~isnan(Y.hbo)), find(isnan(Y.hbo)), 'spline');
Y ist the structure containing the different hemoglobine signals, with rows indicating sample points and columns containing channels'data. (I attached the mat file.)
NaNs are replaced, e.g., with method 'nearest', 'linear', or 'pchip'. Only 'spline' doesn't replace any NaN.
Any idea what is going wrong?
2 Comments
Mathieu NOE
on 10 Jun 2021
hello
which matlab release are you using
I do remember that in older versions, interp1 with spline method was not able to cope with NaN values
that may have been changed as my doc of R2020b does not speak about this limitation;
Accepted Answer
Walter Roberson
on 10 Jun 2021
Y.hbo is 3781 by 24. You are using logical indexing on it, which is going to give you a column vector result in which all of the columns are mixed together. You then ask to do interpolation on all of that at the same time. This is not going to give you the kinds of answers that you want. Consider that if you have a NaN first thing in a column, you would be interpolating from a different column.
You should be interpolating channel by channel. And because NaN could occur at edges, you should specify extrapolation.
More Answers (1)
Steven Lord
on 10 Jun 2021
The fillmissing function offers several methods to use in filling in the missing data (NaN in the case of your numeric data.) One of those options is 'spline'.
0 Comments
See Also
Categories
Find more on Spline Construction in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!