CC BY
9Journal of Siberian Federal University. Engineering & Technologies 6 (2014 7) 648-654
УДК 621.396.44
Method and Device for Compression of the Non Stationary Signals
Sergey P. Panko*, Andrey V. Mishurov and Valentina G. Andyuseva
Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041, Russia
Received 17.06.2014, received in revised form 21.07.2014, accepted 20.08.2014
In the report a new method and the operation of the device ofstatistical compression ofnon-stationary pulse signals are described. The pulse signals are characterized by, first of all, weak irregular period of signal, e.g. electro cardio signals (ECS). It is shown that construction and processing the ECS histogram provide high value of the compression.
Keywords: signal compression, ECS, electro cardio signals.
Метод и устройство
для сжатия нестационарных сигналов
С.П. Панько, А.В. Мишуров, В.Г. Андюсева
Сибирский федеральный университет Россия, 660041, Красноярск, пр., Свободный, 79
Описан новый способ и представлена работа устройства статистического сжатия нестационарных импульсных сигналов. Импульсные сигналы характеризуются в первую очередь слабым нерегулярным периодом сигнала, например электрокардиосигналы (ЭКС). Показано, что создание и обработки гистограммы ЭКС обеспечивают высокое значение сжатия.
Ключевые слова: сжатия сигналов, ЭКС, электрокардиосигналы.
Research motivation for the signal compression
The strategic task for signal compression has a long history and it has become the most popular in the era of mass distribution of digital transfer information because it provides, the greatest reduced traffic.
© Siberian Federal University. All rights reserved
* Corresponding author E-mail address: [email protected]
Compression of digitized almost periodic non-stationary ECS signals is a very important area of telemedicine services, especially in remote Holter monitoring associated with the transmission over loaded communication networks.
Non-stationary ECS results fro m the influence of respiratory and other human activities on the heart, violations of ehythm, whica are called "arrhythmoa". Arrhythmias are expressed in irregular occurrence of cardiac contractions ("extrasystoles") or in a sharp and short-term increase in the duration of RR interval. Individual features of the patient influence greatly the unsteadiness of ECS.
Research of ECG signal compression is motivated, first of all, by decrease of traffic and its costs.
Effectiveness af the compressio n procedures
The effectiveness of the compression procedure is defined by the following:
• compression Ratio (CR) is defined as the ratio of the number of bits required to represent the original (Org) and compressed signals;
• the differences befween restored (Rec) signal and original (Org):
PRD = 100% pjyo-^eoP
where N is the length of window in the number of samples in the signal.
It is cfear thet using rompression does ttot lead to the lors of information PRD
A brief overview of the most common compression techniques
Theoretical and applied researches in this field have been studied by many experts, the results of the researcW are published in mano works. (One of the esrliest eompression techniques of digitize ECS is the AZTEC (the Amplitude Z one Epoch Time Coding) [1]. The ECS is considered to consist of flat aoeas and "tilts".
AZTEC provides a high cempression ratio (CR = 10:1 af a frequency of samples 500 Hz, resolution ADC 12 t)it/sample), but the quality oe the nertored signal is very low.
In [2] they reduce the average value of RR intervals from a fixed realization of the RR intervals. Residua- daea is encoded u sing first differe nces.
In [3] tee curva-ure of tire digitized signal at the points of each signal "turn" (breaks) is evaluated in the distinguished point Curvature at the pnirt (X, Y) is an arbitrary curve which is non-linear combinaiion og the first and reeond derivatives. In two--dimensional space it is determined by the focmute:
(1)
=0%.
Cur if-^---(2)
Subject of the transmission, is not all the digitized signal sampling, but only specific points. Reconstruction of the ECS is mede by interpolation oC transferred distinguished points with a straight line, polynomial, cubic spline or nnother functions.
Passing parameters of the approximation of a Bezier curve also provides a reduction in traffic, —ut, of foutse, increas es the amount of computational pro c e dure s in the proce ss of restoration of the initial signal.
Great attention of researchers is paid to the use of wavelet-transformation properties to reduce traffic, and particularly the application of the SPIHT algorithm to encode the coefficients of the wavelet transform of the ECG signal [5]. SPIHT Algorithm (Set Partitioning in Hierarchical Trees) is incorporated method of coding and is the further development of the EZW algorithm (The embedded zero tree wavelet).
In [5] studied a two-step procedure compression, the first of which is decorrelation of the compressed signal, and the second is traditionally entropy encoding. Table 1 summarizes the results of compression using the various techniques described above.
Analysis of the results and development of ECS compression technologies suggests that one stage compression procedure reached its compression ratio limit. Therefore, further efforts should be focused on search of multi-stage compression.
One of the possible ways is following two-step compression. First look to Table 2, which contains, for example.
A brief sampling of the MIT-BIH basis is presented in Table 2 [6]. Left column lists the values of heart rate range HRR (58...72) beat/minutes. Digits of 100...104 indicate the number of record in MIT-BIH database. The number of occurrences of RR intervals specified in the heart rate during the measurement of the order of 1 minute is indicated on the intersections of HRR and record number.
As clearly seen, the number of identical duration of periods of RR changes greatly. If we assume that the 368 indicates the number of the most often intervals, it is clear that the ECS No. 102 has a much more concentrated form than, e.g. record No.101.
Operation of the device
The sensor (S) converts low electric potential from patient's skin into ECS by filtering frequency bandwidth of up to 200 Hz and gain factor around 1000 in BFA block.
BFA Output is connected to the ADC - analog-digital Converter - and Fmax block that generates a short pulse when reaching the ECS maximum. The pulse indicates the maximum of QRS complex.
There is a sequence of short pulses with a period equal to the duration of the RR intervals in BFA output. Each RR-interval in the sequence has a unique number corresponding to the ordinal number that is stored in the counter.
ADC starts with the frequency near 1 kHz. Digital codes of samples of the input signal from output of ADC, enter the filter FDBL, which eliminates drift baseline (patent of Russian Federation № 2417050).Estimation of the deviation of the EСS middle line in the area fee of QRS complex, predicts a contour shift within the nearest QRS complex and to compensate the drift with fairly high accuracy.
Table 1
Algorithm PRD CR
Wavelet 18.2 (RMS) 21.4
SPHIT 6.49 20
JPEG2000 3.26 20
The combined algorithm [4] 1.17 18.27
Table 2
HRR Record number from MIT- BIH
100 101 102 103 104
58 1 135 0 0 0
59 3 187 0 1 0
60 7 230 0 1 0
61 4 188 0 5 0
62 4 218 2 10 0
63 4 133 0 32 0
64 3 127 3 79 4
65 0 89 2 143 9
66 3 70 17 179 8
67 3 46 8 231 22
68 5 33 38 213 33
69 20 40 143 269 79
70 56 36 219 200 103
71 137 24 292 199 146
72 229 28 368 188 298
DC
S BFA Fmax —> MP
ADC
FDBL
BC SH
Meml-
4
Huffman coder or other
■ Mem 2
FS
a b
e
Start
Fig. 1. The Structural scheme of the device for ECS compression by statistical processing
Digital equivalents of all samples ECS obtained in the process of measuring and purified from the drift contour are stored in memory node (Mem1). Frequency synthesizer FS is triggered by a signal Start and generates a sequence of all the control pulses, marked by the letters a...j. The duration of each RRi interval is measured in block MP (Meter period) by counting the number of quantizing pulses along the "e" line.
The SH block produces a histogram in coordinates: quantile of histogram (bin, bunker) and the number of RRi hits in bin. To plot the histogram we need to have the duration of the bin proportional to 1 beat/min. Whole range of possible values of RR-intervals is 50...200 beats/min. At obtaining RRi, not a multiple of an integer, rounding by the usual rules. RRi intervals are accumulated in the bins corresponding to their numbers. Frequency distribution of RRi intervals in the form of a histogram is shown in Fig. 2. Initial data for histogram are taken from the database of the MIT-BIH [6]. As seen from Fig. 2, there is r significant eonrentration of s RRi interval duration. There is no need to send samples oU all IfRi intervals, bur transmit only a sample. The recovery of the original ECS in the receiver essures a balance oh intehvals in accordance with tSieir orninal oumber.
However, the same period duration of RR intervals do not mean identical goes inside the interval, so you mi ust mapping RR members within the bin on the criteria of maximum correlation coefficient e.
For it the control unit BC reads rooms RR tntervals from each bin. In accordance with tne chosen number interval from a block of memory Mem 1 read audio samples within this interval, and the correlator Cor calculated the correlation coefficient between adjacent intervals.
TOe exceeding of the predefined threshold e0 decision on the pnoximity of bhe RR-intehvals, which allows counts one RR interval and rooms close io this value other RR internals within the bin.
Fig. 3 shows the dependence of quantity of RR-intervals, close by the criterion threshold e0 within one bin.
In the block Dec decided on the proximity of RR intervals, relntivcly to which thh correlation coeffi cient is calculated.
Fig . 2. Frequency distribution of RR intervals, where: M - the number of RR intervals of equal duration; HR - RR interval duration, beats / min.
Fig. 3. Number of RR interval s in the aisle s of the hopper after comparison with a correlation coefficient s0 = 0,955, where: M - the number of RR intervals of equal duration; HR - RR interval duration, beats / min.
In the result of the simulation algorithm described above compression on the ECS records from the database of the MIT-BIH [6] was obtained CR=6,269 by a factor of s<,= 0,95 and PRD=0.05 and CR=20.746) for factor s0= 0,90 and PRD=0.1. It should br noted that tlie efficiency of the considered algorithm dtpends heavily from statistics of the sigaal you want to compress.
Conclusion
The rrHability erf temotr me dical monitoring systems depends on the success of data transmission, in particular, the compression efficiency of non-stationury signals. For that purpose the most promising is to usr of multistage procedures.
Work is performed with financial support of the Ministry of Education and Science of the Russian Federation in the Stierian federaluniversity (Contract No. 02.G25.31.0041).
References
[1] Cox J.R. et al. AZTEC: A Preprocessing Program for Real-Time ECG Rhythm Analysis, IEEE Trans. Biomed. Eng., 15, 128-129,1968.
[2] Hamilton P. S, Tompkins W.J. Compression of the ambulatory ECG by average beat subtraction and residual differencing. IEEE Trans Biomed ling. 199)1 Mart; 38(3): 253-9.
[3] Sweeney; R.J. Signal compressing based on curvature parameters. US Patent Application 20090t92395. July 30, 2009 ,
[4] BritoM., Henriques J., Carvalho P., Ribeiro B. and AntunesM. An ECG compression approach base on a segment dictionary and Besier approximatios. 15th European Signal Processing Conferrncr (EUSIPCO 2007), Poznan, Polfnd, September n7, 2S07.
[4] Ridha Iskandar I, Wayan Simri W. Compression of ECG Signal Using Neural Network Predictor and Huffman Coding . Universitas Gunadarma (Indonesia), Proceeding Seminar Ilmiah Nasional KOMMIT, 24-11-2010.
[5] PooyanM., Taheri A., Moazami-GoudarziM., Saboori I. Wavelet Compression of ECG Signals Using SPIHT Algorithm. International Journal of Information and Communication Engineering 1:4, 2005.
[6] Portoles L.B. Lossles Compression of ECG signals. Performance Analysis in a Wireless Network. 2009-09-15. LiTH-IMT-ERASMUS-R-09/34- SE.
