Segmenter un signal ECG

Version imprimable

3 pièce(s) jointe(s)

Bonjour / bonsoir à tous

Je suis actuellement en train de traiter un signal comportant 10 "battements cardiaques"
Je souhaiterais diviser mon vecteur en dix autres entités de manière "automatique".
Dans un premier temps j'ai appliqué un seuil d'amplitude pour supprimer le bruit lorsqu'il n'y a pas de perturbation importante en amplitude correspondant à mes signaux cardiaques.

Voici le vecteur initial, suivi d'un autre vecteur sans le bruit quand le signal est au "repos".
Pièce jointe 146519Pièce jointe 146520

Voici le visuel des deux signaux
Pièce jointe 146523

Le code pour générer tout ça
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 v1 = pas_mike_150k_01.Y(1,1).Data; j=1; k=1; m =0; Nv1 = length(v1); for i = 1 : Nv1 Tv1(1,i) = j/10000; j = j+1; end for j =1:Nv1 if abs(v1(1,j)) >= 0.025 V1(1,k) = v1(1,j); k = k+1; m = 1; elseif abs(v1(1,j)) < 0.025 && m <=10 m = m+1; V1(1,k) = v1(1,j); k = k+1; elseif abs(v1(1,j)) < 0.025 && m >10 V1(1,k) = 0; k = k+1; m = m+1; end end figure(4) plot(Tv1,v1,'b','LineWidth',3) grid on hold on plot(Tv1,V1,'r','LineWidth',3)

28/05/2014, 01h35
Paul-Ceria

Citation:

Je souhaiterais diviser mon vecteur en dix autres entités de manière "automatique".

Concrètement qu'est ce que tu souhaiterais comme résultat final ?
28/05/2014, 16h30
Mizukage Tobi

Je voudrais récupérer dans 10 vecteurs différents mais de taille égale les dix "battements cardiaques" du signal initial.

J’espère que ça te convient, par contre mon code est brouillon :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
clear;clc;
load('pas_mike_150k_01.mat')
v1 = pas_mike_150k_01.Y(1,1).Data;
j=1; k=1; m =0;
Nv1 = length(v1);
for i = 1 : Nv1
    Tv1(1,i) = j/10000; j = j+1;
end
 
for j =1:Nv1
    if abs(v1(1,j)) >= 0.025
        V1(1,k) = v1(1,j);
        k = k+1; m = 1;
    elseif abs(v1(1,j)) < 0.025 && m <=10
        m = m+1;
        V1(1,k) = v1(1,j);
        k = k+1;
    elseif abs(v1(1,j)) < 0.025 && m >10
        V1(1,k) = 0;
        k = k+1;
        m = m+1;
    end
end
 
figure(4)
% plot(Tv1,v1,'b','LineWidth',3)
grid on
hold on
plot(Tv1,V1,'r','LineWidth',3)
 
%% Battements dans vecteurs
vecteur=struct( 'vect', {[]});
tempbattement=1800;
nvecteur=0;
j=0;
 
nbbat=ceil(size(find(V1(1:size(Tv1,2)))~=0,2)/tempbattement)+2;
fprintf('Nombres de battements : %d \n', nbbat)
k=1;
 
for i=1:size(Tv1,2)
    j=j+i;
    if k<nbbat+1
        if V1(j)~=0
            nvecteur=nvecteur+1;
            vecteur(nvecteur).vect(1,:)=Tv1(j-i+1:j-i+1+tempbattement);
            vecteur(nvecteur).vect(2,:)=V1(j-i+1:j-i+1+tempbattement);
 
            quiver(Tv1(j-i+1),0,0,1,'b');
            quiver(Tv1(j-i+1+tempbattement),0,0,1,'g');
 
            j=j+tempbattement;k=k+1;
        end
    end
end
 
for i=1:nbbat;
    figure
    plot(vecteur(i).vect(1,:),vecteur(i).vect(2,:))
end