Automatic Signal segmentation for feature extraction

14 次查看（过去 30 天）

显示更早的评论

Radons 2018-1-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/378207-automatic-signal-segmentation-for-feature-extraction

评论： Image Analyst 2022-11-4

采纳的回答： Image Analyst

Signal.mat

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 个评论
显示无隐藏无

Star Strider 2018-1-21

What is the corresponding time vector that corresponds to those samples, or the sampling frequency?

Radons 2018-1-21

It is sampled at 1Khz

请先登录，再进行评论。

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

采纳的回答

Image Analyst 2018-1-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/378207-automatic-signal-segmentation-for-feature-extraction#answer_301083

在 MATLAB Online 中打开

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.

19 个评论
显示 17更早的评论隐藏 17更早的评论

Cey 2022-11-4

thank you for your help @Walter Roberson

Image Analyst 2022-11-4

@Cey by now I imagine you've found one. You can just use trial and error until you get the right one. If you have any more questions, start your own discussion thread so we're not bugging @Radons on his 4 years old question with emails about activity on it.

请先登录，再进行评论。