Charger la base de données MIT-BIH Arythmie ECG sur MATLAB

Vous avez besoin du programme rddata.m (script MATLab) de ce site internet. Le programme peut être trouvé ici. rddata.m est probablement le seul programme que vous aurez besoin de lire les signaux ecg. Je me souviens d'avoir utilisé ce programme et de la base de données me suis pas trop longtemps.

Il y a un tutoriel pour l'utilisation de matlab pour lire les données. tutoriel pour matlab utilisateur

1. installer "La WFDB Toolbox pour Matlab" à partir du lien ci-dessus. Ajouter le dossier de la boîte à outils pour le chemin d'accès dans matlab.

2. Télécharger le signal ECG. Assurez-vous de télécharger '.atr', '.dat' and '.hea' ensemble pour le signal vous sont traiter avec.

3. De commande de matlab est comme suit :
   
[tm,signal,Fs]=rdsamp( filename , 1 ) ;
[ann,type]=rdann( filename , 'atr' ) ;

   Remarque: pour le signal '101', son nom est '101'. Et vous pouvez consulter le détail des informations sur rdsamp et rdann du tutoriel.

Utilisation ATM pour extraire .tapis de comme décrit par Kamtal (maintenant connu Rashid).
Toutefois, notez que pour voir la .info fichier dans certains cas, vous devez cliquer sur la flèche

Après que j'ai poussé en avant pour les développeurs icinous avons obtenu des améliorations dans la documentation ici dans la Section 4.

Si ils sont tous les entiers de l'intervalle [-2^N, 2^N-1 ] ou [ 0, 2^N ], ils sont probablement numérique. Comparer les valeurs pour voir si elles sont dans la gamme physiologique du signal analysé. Par exemple, si l'en-tête stipule que le signal est un ECG stockées dans milivolts, qui a généralement une amplitude d'environ 2 mv, un signal d'entiers allant de -32000 à 32000 n'est probablement pas de vous donner le physique de l'ECG dans milivolts...

Si elles ne sont pas entiers puis ils sont physiques. Une fois de plus, vous pouvez rapidement comparer les valeurs pour voir si elles sont dans la gamme physiologique du signal analysé.

0-9-10 bdcea - unités physiques

Nous dire que les signaux soient en unités physiques " lorsque les valeurs sont utilisées pour représenter les vrai valeurs de la vie d'aussi près que possible, bien que, évidemment, tout sur l'ordinateur est numérique et discret plutôt que de l'analogique et du continu. Cela inclut nos précieux 64 bits en virgule flottante double précision des valeurs, mais c'est aussi proche que nous pouvons obtenir et déjà très proche de la physique réelle des valeurs, donc nous référer à eux comme "physique".

-

Par exemple, si un 15 bits du signal est recueilli par l'intermédiaire d'un dispositif de capture de l', Physionet sera probablement stocker en tant que 16 bits du signal. Chacun des 16 bits du bloc stocke une valeur entière comprise entre -2^15 et 2^15-1, et à l'aide du gain et de l'offset indiqué dans l'en-tête pour chaque canal, l'origine physique du signal peut être tracée pour le traitement.

Les unités par défaut sont maintenant unités physiques où la base et de gain doit être ajoutée indiqué dans l'en-tête pour chaque canal, de sorte que le signal physique peut être tracée pour le traitement.

% rawUnits
%       A 1x1 integer (default: 0). Returns tm and signal as vectors
%       according to the following values:
%               rawUnits=0 - Uses Java Native Interface to directly fetch  data, returning signal in physical units with double precision.
%               rawUnits=1 -returns tm ( millisecond precision only! ) and signal in physical units with 64 bit (double) floating point precision
%               rawUnits=2 -returns tm ( millisecond precision only! ) and signal in physical units with 32 bit (single) floating point  precision
%               rawUnits=3 -returns both tm and signal as 16 bit integers (short). Use Fs to convert tm to seconds.
%               rawUnits=4 -returns both tm and signal as 64 bit integers (long). Use Fs to convert tm to seconds.

rawUnits=1rawUnits=2 aussi utiliser des unités physiques.
rawUnits=3rawUnits=4 utilisation puis de nouveau analogique/numérique où vous avez besoin de supprimer de la base et du gain.
Si vous utilisez rawUnits=1 ou rawUnits=2, vous avez besoin d'ajuster de base et d'acquérir où base = 1024 et gain = 200

# Kamtal's method in considering base and gain
load('201m.mat');
val = (val - 1024)/200;    % you have to remove "base" and "gain"
ECGsignal = val(1,16:950); % select the lead (Lead I)

Voir le .info fichier ci-dessous où vous pouvez obtenir le base et gain. Il est également l'unité de mV qui suggère que les valeurs doivent être près de 2 après la base de gain d'opérations.

<0-9-9 bdcea - analogique/numérique unités de base et d'acquérir par défaut; maintenant seulement rawUnits=3,4 pour les unités analogiques

Après la sélection de l'ATM, vous devriez être en mesure de voir la liste dans laquelle vous pouvez sélectionner .info fichier après l'exportation, comme décrit dans Kamtal de réponse.
Le .info fichier indique à supprimer les base et gain à partir des données avant de les utiliser

Source: record mitdb/201  Start: [00:02:10.000]
val has 2 rows (signals) and 3600 columns (samples/signal)
Duration:     0:10
Sampling frequency: 360 Hz  Sampling interval: 0.002777777778 sec
Row     Signal  Gain    Base    Units
1       MLII    200     1024    mV
2       V1      200     1024    mV

To convert from raw units to the physical units shown
above, subtract 'base' and divide by 'gain'.

La comparaison de mauvaises réponses! [Obsolète]

Kamtal (maintenant appelé Rashid) réponse est sur le vieux bdcea système utilisé unités numériques sans enlèvement de base et d'acquérir

# Kamtal's method in considering base and gain
load('201m.mat');
val = (val - 1024)/200;    % you have to remove "base" and "gain"
ECGsignal = val(1,16:950); % select the lead (Lead I)

# Method without considering base and gain
load('201m.mat');
ECGsignal2 = val(1,16:950); 

# http://www.mathworks.com/matlabcentral/fileexchange/10502-image-overlay
imshow(imoverlay(ECGsignal, ECGsignal2, uint8([255,0,0])))

et vous obtenez la différence entre ma méthode et de sa méthode