compression.m

ECG压缩中的直接域压缩CORTES算法

clc;clear all;
%------ SPECIFY DATA ------------------------------------------------------
PATH= 'E:\lipai\ECG\MIT-BIH Arrhythmia Database'; % path, where data are saved
% RecordName = 'cu01';         % record name
RecordName = '100';         % record name
SAMPLES2READ=1500/2;
%SAMPLES2READ=127232/2;
%SAMPLES2READ=3000;         % number of samples to be read
                            % in case of more than one signal:
                            % 2*SAMPLES2READ samples are read
HEADERFILE=[RecordName, '.hea'];      % header-file in text format
ATRFILE=[RecordName, '.atr'];         % attributes-file in binary format
DATAFILE=[RecordName, '.dat'];        % data-file

                            
%------ LOAD DATA --------------------------------------------------
InfoForRead.path = PATH;
InfoForRead.headerfile = HEADERFILE;
InfoForRead.atrfile = ATRFILE;
InfoForRead.datafile = DATAFILE;
InfoForRead.samples2read = SAMPLES2READ;

record = rdrecord_hailang(InfoForRead);
%[sfreq, NoSig, Signal, Time, AnnTime, AnnType] = rdrecord_hailang(PATH, HEADERFILE, ATRFILE, DATAFILE, SAMPLES2READ);

sfreq = record.sfreq;
NoSig = record.nosig;
Signal = record.signal;
Time = record.time;
AnnTime = record.anntime;
AnnType = record.anntype;
%-------------COMPRESSION--------------------------------------------------
for i=1:2
    for k=1:1499/2
        delta((k+1),i)=Signal((k+1),i)-Signal(k,i);
    end
    delta(1,i)=Signal(1,i);
end
%----------------------DECODE----------------------------------------------
for a=1:2
    V(1,a)=delta(1,a);
    for b=2:1500/2
        V(b,a)=delta(b,a)+V((b-1),a);
    end
end
%-------------------DISPLAY------------------------------------------------
figure(1); clf, box on, hold on
plot(Time, Signal(:,1),'r',Time, Signal(:,2),'b');
figure;
plot(Time, V(:,1),'r',Time, V(:,2),'b');
figure;
plot(Time, delta(:,1),'r',Time, delta(:,2),'b');