Cut a signal into heart cycle pieces

Question

Jawad Jawadi 2018-12-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437530-cut-a-signal-into-heart-cycle-pieces

评论： Chris Turnes 2019-1-7

Hi guys,

I have a signal from a photophlethysmographic (PPT) sensor showing contours of the heart pulse.

My signal needs to be cut/subdivided into signal pieces for a single heart beat. Each single heart beats needs to be saved in the workspace.

So when the slope of the signal is getting values different than 0, the data should be cut and saved to workspace.

Is there any possibilities I could do so ? Does anybody have suggestions ?

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

Jawad Jawadi 2019-1-3

Hi Naman,

first of all , thank you for your reply.

If this is my signal:

Then I would like to cut the signal into following segments:

In order to cut the signal into single undulations:

The goal is cut the signal according to the above mentioned undulations and save it to the workspace.

Later I want to normalize the signals and stack them to calculate and disply a density map.

Here is my code so far:

________________________________________________________

dataset = xlsread('data.xlsx','EA0165','A1:G1578');

Time=dataset(:,1);

A1=dataset(:,2);

V2=dataset(:,3);

A3=dataset(:,4);

V4=dataset(:,5);

Timesystol=dataset(:,6);

A = fillmissing(V2,'linear');

B = movmedian(A,3);

detrended_B=detrend(B);

trend = B - detrended_B;

C = movmedian(detrended_B,3);

D =movmedian(C,[5 3]);

figure(4)

hold on

plot(D,'g');

____________________________________________________________

Can you help me ?

John D'Errico 2019-1-3

编辑：John D'Errico 2019-1-3

This is not an answer, since I don't have your data signal, nor even the signal processing toolbox. But since it looks like you want to use the local minima to break the signal up, just use findpeaks.

Find the locations of each local minimum. help findpeaks
Use those locations to break your signal into pices, probably as independent cells orf a cell array.
Subtract off the local average value for each section. This will detrend the problem, helping them to now overlap.

You do NOT want to create the segments as separate variables in your workspace. That is a particular place in programming hell where you do not want to live. Instead, learn to use arrays, and here, cell arrays or even structures.

请先登录，再进行评论。

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

Answer 1

Chris Turnes 2019-1-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437530-cut-a-signal-into-heart-cycle-pieces#answer_354833

在 MATLAB Online 中打开

example.mat

I don't have your data, but using something that looks roughly similar (attached), here's an approach. Furthering John's comment, if you have access to R2017b or beyond, you can use islocalmin to find the minima.

>> tf = islocalmin(x, 'FlatSelection', 'center');
>> t = 1:length(x);
>> plot(t, x, t(tf), x(tf), 'rx')
>> legend('Input', 'Local minima')

Once you have them, you can do a cute little trick to segment the signal: you can call cumsum to get unique labels for each heartbeat:

>> glabels = 1 + cumsum(tf);

With these labels, you can uniquely identify any individual heartbeat you've detected without having to explicitly construct each one. For example, to plot the 4th heartbeat you've segmented:

>> idxFourth = glabels == 4;
>> plot(t, x, t(idxFourth), x(idxFourth));
>> legend('All', 'Fourth')

You could also mean-center each heartbeat as John suggests by using the grouptransform function of R2018b (though you'll need to pack your data into a table):

>> Tx = table(x', glabels', 'VariableNames', { 'Data', 'Heartbeat' });
>> Txmc = grouptransform(Tx, 'Heartbeat', 'meancenter');
>> plot(t, Txmc.Data)

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

Jawad Jawadi 2019-1-7

Dear Chris,

thanks for your quick reply.

The function grouptransform seems not to be accepted by MATLAB.

I get the reply that the function or variable 'grouptransform' is undefined.

Do I have to previously define grouptransform ?

I am on Matlab 2018a.

Chris Turnes 2019-1-7

在 MATLAB Online 中打开

grouptransform was introduced in R2018b, unfortunately, so it seems you won't have access to that function. An alternative approach is to simply loop through the number of groups you have and manually re-normalize using indexing. For example, something like

for k = 1:max(glabels)
    thisGroup = D(glabels == k);
    D(glabels == k) = thisGroup - mean(thisGroup);
end

请先登录，再进行评论。

Cut a signal into heart cycle pieces

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

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Cut a signal into heart cycle pieces

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

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

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

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