How to find the first edge (not minimum value) after the peak?

Question

Vanessa Phung 2017-7-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/350356-how-to-find-the-first-edge-not-minimum-value-after-the-peak

回答： Star Strider 2017-7-26

The subplot 1 shows a line profile across a density column. I want to find the radius of this column by getting the coordinate of peak(max) and the first edges. The problem is that the edge is not at minimum value, therefore I can't use min().

The subplot 2 is obtained from findpeaks().

figure; imagesc(FT_phase_shift), colorbar;title('Phase shift'); aline=FT_phase_shift(1:600,700); amax=max(aline);amin=min(aline) pmax =find(aline==amax(:,1));pmin=find(aline==amin(:,1));x_max = aline(pmax);

% find peaks and location [pks, locs]= findpeaks(aline); figure;subplot (211); plot(aline); subplot(212); plot(locs, pks);

Please comment.Thank you.

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2017-7-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/350356-how-to-find-the-first-edge-not-minimum-value-after-the-peak#answer_275617

在 MATLAB Online 中打开

I would use findpeaks on the negative of your signal:

[negpks, neglocs]= findpeaks( -aline );

Then take the first ‘neglocs’ value less than the maximum peak ‘locs’ value, and the first ‘neglocs’ value greater than the the maximum peak ‘locs’ value. Use the corresponding ‘negpks’ values (as ‘-negpks’) to find the values you want.

Example —

x = linspace(0, 10*pi);                                 % Create Data
aline = abs(sin(x-5*pi)./(x-5*pi))*7;                   % Create Data
[pks, locs]= findpeaks(aline, 'MinPeakHeight',5);
[negpks, neglocs]= findpeaks( -aline );
minidx = find(neglocs < locs, 1, 'last');
maxidx = find(neglocs > locs, 1, 'first');
minloc = neglocs(minidx);
maxloc = neglocs(maxidx);
figure(1)
plot(x,aline,'-b', x(locs),pks,'^r')
hold on
plot(x([minloc maxloc]), aline([minloc maxloc]), 'vr')
hold off