MATLAB Answers

Automatic Signal segmentation for feature extraction

112 views (last 30 days)
Radons
Radons on 21 Jan 2018
Answered: Hristo Zhivomirov on 13 Oct 2019
Hi,
Does anyone know how to do signal segmentation on the raw signal? I need to segment the raw signal into 8 different segments so that i can do feature extraction on individual segments.
  2 Comments

Sign in to comment.

Accepted Answer

Image Analyst
Image Analyst on 21 Jan 2018
Attach some signals so we can play with them.
Assuming the "tall" parts of the signal are not uniformly spaced (in which case you could extract them trivially with indexing at fixed indexes), then you need to find the center of the quiet parts. I'd start by thresholding the absolute value of the signal to find the low/quiet parts. You should have 9 of them, unless your tall signal parts hit the edge of your range. If you don't have 9, then I'd call imdilate repeatedly until you get exactly 9 quiet regions. Then I'd call regionprops() to get the centroid (middle index) of the quiet parts, and then I'd loop over the regions extracting the tall regions into cells. They need to be in cells because each region might have a different number of elements. Something like (untested)
quietParts = abs(signal) < threshold;
[~, numRegions] = bwlabel(quietParts);
while numRegions > 9
quietParts = imdilate(quietParts, [1,1,1]);
[~, numRegions] = bwlabel(quietParts);
end
% Now we should have the quiet parts. Find the centroids.
props = regionprops(quietParts, 'Centroid');
% Extract each tall signal part into a cell
for k = 1 : length(props)-1
index1 = round(props(k).Centroid;
index2 = round(props(k+1).Centroid;
individualSignals{k} = signal(index1:index2);
end
Note, the code above requires the Image Processing Toolbox.
  11 Comments
Image Analyst
Image Analyst on 25 Jan 2018
You can try sgolayfilt with a shorter window to smooth it, then perhaps use findchangepts() to identify where it starts climbing.

Sign in to comment.

More Answers (2)

Jehad Ali
Jehad Ali on 30 Mar 2019
Can any one help me how to write code segmentation for arterial blood pressure pulse in matlab ?
  1 Comment
Image Analyst
Image Analyst on 30 Mar 2019
Yours is a new question, not an "Answer" to Radons question. Start a new thread.

Sign in to comment.


Hristo Zhivomirov
Hristo Zhivomirov on 13 Oct 2019
Hi, Radons!
I think that this the most suitable, easy-to-use and straightforward way.
Also, you can use the example.m file as reference.
All best,
Hristo

Community Treasure Hunt

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

Start Hunting!

Translated by