Remove evenly distributed spikes from curve

2 次查看（过去 30 天）
Jan Heinsoth
Jan Heinsoth 2021-4-27
评论： Star Strider 2021-4-29
Hi, I got some data from the lab giving the height data of holes in a sample I put under a microscope, so I could get a profile of the sample. Since I expect regularly occurring holes in the sample i tried to find the places with the steepest slope, by maxima and minima of the gradient of the data. I found that the gradient curve has many evenly distributed spikes on top of a function that is otherwise really good looking. I tried a smoothing spline to get rid of the spikes, but it also replicates them. I'll give you my code and the graph below. Please help me to get rid of those spikes. Thanks in advance.
% reading in the data from the file
calib5 = importdata(['C:\Users\jan\OneDrive - Carl von Ossietzky Univers',...
'ität Oldenburg\UniShit\Engineering Physics\Semester III - winter 20',...
'20-2021\Laboratory Project I\Data from the Lab\5_100x_zStapelBericht.xlsx']);
% calculating the gradient
c5_1_dDx = gradient(calib5.data.Profil1(:,2))./gradient(calib5.data.Profil1(:,1));
hold on
% creation of the plot above
plot(calib5.data.Profil1(:,1),sum([[diff(calib5.data.Profil1(:,2));0],...
[0;diff(calib5.data.Profil1(:,2))]],2)) % plot of the gradient replicated using diff()
plot(calib5.data.Profil1(:,1),c5_1_dDx) % plot of the gradient
hold off
PS: I really need to find the prominent local maxima and minima of this and not just the spikes.
  4 个评论
显示 2更早的评论
Jonas
Jonas 2021-4-27
编辑：Jonas 2021-4-27
maybe a moving median is also a good alternative, looks like the peaks you want to remove are only one sample long, so you can try a moving median with a window length of 3
Rik
Rik 2021-4-29
Comment posted as flag by @Jan Heinsoth:
The proposed lowpass worked as well.

请先登录，再进行评论。

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

采纳的回答

Star Strider
Star Strider 2021-4-27
See if the filloutliers function will do what you want.
  2 个评论
Jan Heinsoth
Jan Heinsoth 2021-4-29
编辑：Jan Heinsoth 2021-4-29
Thank you, It did not perfectly smooth the signal, but now it is smooth enough for my purposes.
But the lowpass filter proposed by @Jonas did also work for me.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Smoothing and Denoising 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by