Mathematical, Algorithmic and Software Support for Phonocardiographic Signal Processing to Detect Mitral Insufficiency of Human Heart Valves

Abstract

Mathematical support for phonocardiographic signal processing has been developed based on a mathematical model in the form of a periodically correlated stochastic process and a component processing method. On the basis of mathematical support, algorithmic and software was developed in the Matlab environment for automated systems of computer diagnostics of the functional state of the valves of the human heart when mitral insufficiency is detected. Mitral valve insufficiency of the human heart, phonocardiographic signal, mathematical support, periodically correlated stochastic process, component method, algorithmic support, Proceedings ITTAP'2023: 3rd International Workshop on Information Technologies: Theoretical and Applied Problems, November 22-24, Proceedings

1. Introduction

According to the WHO and the Medical Association of Cardiologists, it has been established that the development trend of human heart valve lesions occupies a dominant place among all cardiovascular diseases in Europe [ 1-19 ]. The primary causes of such lesions are congenital (arising during fetal development) and acquired (complications of rheumatism, infectious endocarditis or after surgical treatment of mitral stenosis) defects of the human heart.

In medical practice, such basic non-invasive instrumental research methods as cardiography [ 1,2 ], auscultation [17,18], echocardiography [ 4,6,13,14,15 ] and phonocardiography [17,19] are used to diagnose mitral insufficiency. Phonocardiography along with other methods makes it possible to objectify unclear sound phenomena of human heart valves (degree of weakening and presence of heart sounds, intensity, duration and form of systolic murmur) (Sacks, Roberts and Evans) according to registered phonocardiographic signals (PCG signals) (Figure 1). with the noise of systole (indicated by the pointing arrow)

2020 Copyright for this paper by its authors. CEUR

ceur-ws.org

Of particular value is the possibility of long-term (over a year or more) objective dynamic monitoring of the PCG signal during the formation of a valve defect. To implement the phonocardiography method in medical practice, hardware devices in the form of computer phonocardiographs are used, such as Cardio+ (Ukraine, Metekol LLC), PCG-02 (India, Mumbai), Techbook Scientech 2356 (India, Scientech Technologies Pvt . Ltd.) and Audio-Technica ATR288W (Japan). Well-known phonocardiographs based on the principle of obtaining diagnostic information about the condition of human heart valves are built according to the scheme "bioobject (human heart valves)-mathematical support (mathematical model-processing method)-algorithmic support-softwarediagnosis result" [21,24,26,28,29,38,39,40]. In phonocardiographs, software and algorithmic support, which form the result of diagnostics, are implemented on the correlation, spectral and synphase methods of FCS processing, the core of which is a mathematical model in the form of a deterministic harmonic function (Dodge and Cabot), periodic functions of different frequencies (Manheimer), almost periodic functions (Kasyrskyi G.I.), a stationary stochastic process (Metin Akay (Houston)), a mixture of a stationary stochastic process and a deterministic function (Metin Akay (Houston)), periodically correlated stochastic process (Osukhivska H.M., Dragan Ya.P. [22,23,35,36], Palianytsia Y.B. [31,32]), relaxation multipulsator (Osukhivska H.M.) [37]. Such mathematical support (models and methods developed on their basis) of processing by its structure and properties does not provide tracking the temporal dynamics of changes in the amplitude and phase indicators of the functioning of the valves of the human heart, which is important in detecting the manifestations of mitral insufficiency over time.

Therefore, the development of mathematical support (a processing method based on its adequate mathematical model) and algorithmic-software processing of PCG signal for detecting mitral valve insufficiency of the human heart is an urgent scientific task. 2. Mathematical support for phonocardiographic signal processing

The experimentally recorded implementation of the PCG signal in the case of mitral insufficiency is shown in Figure 2.

The PCG signal in mitral insufficiency as a stochastic and cyclical process by the nature of the correlation function satisfies the conditions of cyclicity of the average characteristics and their variations: mt  T   mt  , rt  T , s  T   rt, s. (1) (2)

In this case, T corresponds to the duration of one PCG signal cycle in mitral insufficiency. Expressions (1-2) indicate the cyclic nature of the characteristics, namely that the values of the average characteristics and their variations due to the cardiac cycle are repeated, but with changes caused by chance.

PCG signal in the case of mitral class insufficiency provides the study of its type of harmonicity as a process of stochastic and cyclic nature in the form of PСSP from the standpoint of the energy approach (developer Ya.P. Dragan [20,22,23,30,31,32,33,34]).

Since the PCG signal in mitral insufficiency is presented in a stochastic process with cyclic characters of a finite nature, it can be represented through the expression of translational components [20]: (3) (6)  t     H t  pT  k t   qt  ,

pZ kZ where is a H  single function that indicates the place of generation of the PCG signal tone in mitral insufficiency; T - is the cycle period of the PCG signal in case of mitral insufficiency;  k t PCG signal tones in case of mitral insufficiency within the cardiac cycle; k - cardiac cycle number.

But based on the equivalence of representations through translational components (3) and stationary components, the expression of the PCG signal model in mitral insufficiency has the form [2023,28,31,32]: ik 2 t

T , t R  t    k t e (4)

kZ where are  k t the steady-state components of the PCG signal in mitral insufficiency.

Mitral insufficiency by nature has the appearance of a noise similar in shape to white noise. Such similarity should be reflected in model (4) additively in the form of an expression: ik 2 t  t    k t e T  nt  , t R (5)

kZ where nt is white noise as a sign of mitral insufficiency in the realization of the PCG signal (according to the nature of its formation).

Expression (5) ensures the use of the component method for the processing of the PCG signal in the case of mitral insufficiency as a stochastic process with cyclic characters of a finite nature.

The component method of PCG signal processing based on model (5) ensures the calculation of correlation components according to the expression [20,24,27,28]:  Bk u  1 T 0 0

 t  u  kT  t  kT e T 0 ik 2 t

T dt , t R 0 0 where is  t    t   m t  a centered PCG signal with mitral valve insufficiency; m t  - the average PCG signal in the case of mitral valve insufficiency. 3. Algorithmic support for phonocardiographic signal processing Having considered the expression (6), the main operations of its implementation are highlighted: 1. Finding the T period of the FСG signal in the case of mitral valve insufficiency [25]. 2. Calculation of the estimate of the mathematical expectation m of the PCG signal in the case of mitral valve insufficiency. insufficiency. 5. Counting the components Bˆk u  of the PCG signal in case of mitral valve insufficiency.

In Figure 3, we show all the listed operations in the form of an algorithm, which displays all the operations of implementing the component method of processing. 4. Software and results of phonocardiographic signal processing

Software with a graphical user interface was developed in the MATLAB environment to automate the process of processing the PCG signal when mitral insufficiency is detected.

The results of processing PCG signals in the normal state (Patient A) and in the presence of mitral insufficiency (Patient B) are shown in Figure 4-7.

A comparison of the magnified view of the components from above is shown in Figure 8.

(Patient A - norm)

According to the results of the analysis of the correlation components, a slight change in the power intensity of the components is visible, namely, with mitral insufficiency, the power is more scattered than in the normal state. The components allow only visual detection of changes in the operation of human heart valves without quantitative evaluation.

For the quantitative evaluation of correlation (spectral) components, the evaluation of their averaging by time shifts was used (Figure 9-10).

In Figure 9-10 it is visually and quantitatively visible that the maximum values of the averaged components for a patient in a physiologically normal state (patient A) and a patient with mitral valve insufficiency (patient B) are concentrated at the same frequencies (identical in structure), but the amplitude values of the estimates slightly differ from each other, namely, for patients A and B, the maximum amplitudes are concentrated on the 1st and 7th components. However, an increased power level of the 2nd, 3rd, 4th, and 5th components, as well as an increase in the maximums of the 1st and 7th components, were noticed in patient B. Such changes quantitatively reflect the manifestations of mitral valve insufficiency of the human heart.

Therefore, the calculation of correlation (spectral) components is a quantitative parameter for detecting early changes in the functioning of the valves of the human heart, in particular valve insufficiency.

5. Conclusions

The work substantiates the mathematical model of the PCG signal as a mixture of a periodically correlated stochastic process and noise for solving the problem of detecting the insufficiency of the mitral valves of the human heart. A method and algorithm for processing the PCG signal based on the component method was developed, which provided a procedure for detecting the insufficiency of the mitral valves of the human heart based on the averaged estimates of the correlation components. Software was developed using the MATLAB tool for processing the empirical data of the PCG signal, which provided the process of calculating correlation components as diagnostic signs that quantitatively reflect indicators of mitral valve insufficiency of the human heart. 6. References

References

[1] Baumgartner H , Falk V , Bax JJ , De Bonis M , Hamm C , Holm JP , et al. ESC/EACTS Guidelines for the management of valvular heart disease . Eur Heart J . 2017 ; 38 ( 36 ): 2739 - 2791 . Available from: doi: 10.1093/eurheartj/ehx391.

[2] Dal-Bianco JP , Aikawa E , Bischoff J , et al. Myocardial infarction alters adaptation of the tethered mitral valve . J Am Coll Cardiol . 2016 ; 67 ( 3 ): 275 - 287 . doi: 10 .1016/j.jacc. 2015 . 10 .092.

[3] Deferm S , Bertrand PB , Verbrugge FH , et al. Atrial functional mitral regurgitation: JACC review topic of the week . J Am Coll Cardiol . 2019 ; 73 ( 19 ): 2465 - 2476 . doi: 10 .1016/j.jacc. 2019 . 02 .061.

[4] Deferm S , Bertrand PB , Verhaert D , et al. Mitral annular dynamics in AF versus sinus rhythm: Novel insights into the mechanism of AFMR . JACC Cardiovasc Imaging . 2021 ; S1936 - 878X ( 21 ) 00440 - X . doi: 10 .1016/j. jcmg. 2021 . 05 .019.

[5] Donal , E ; Panis, V ( October 2021 ). "Interaction between mitral valve apparatus and left ventricle. Functional mitral regurgitation: A brief state-of-the-art overview" . Advances in Clinical and Experimental Medicine . 30 ( 10 ): 991 - 997 . doi: 10 .17219/acem/143324. PMID 34714608. S2CID 240154628.

[6] Gertz ZM , Raina A , Saghy L , et al. Evidence of atrial functional mitral regurgitation due to atrial fibrillation: Reversal with arrhythmia control . J Am Coll Cardiol . 2011 ; 58 ( 14 ): 1474 - 1481 . doi: 10 .1016/j.jacc. 2011 . 06 .032.

[7] Grayburn PA , Sannino A , Packer M. Proportionate and disproportionate functional mitral regurgitation: A new conceptual framework that reconciles the results of the MITRA-FR and COAPT trials . JACC Cardiovasc Imaging . 2019 ; 12 ( 2 ): 353 - 362 . doi: 10 .1016/j.jcmg. 2018 . 11 .006.

[8] Eng MH , Wang DD . Transseptal transcatheter mitral valve replacement for post- surgical mitral failures . Intervent Cardiol Rev . 2018 ; 13 ( 2 ): 77 - 80 . Available from: doi: 10.15420/icr. 2017 : 16 : 3 .

[9] Iung B , Baron G , Butchartv EG , Delahaye F , Gohlke-Baerwolf CW , Leavang OW , et al. A prospective survey of patients with valvular heart disease in Europe: The Euro Heart Survey on Valvular Heart Disease . Eur Heart J . 2003 ; 24 : 1231 - 1243 . Available from: doi: 10.1016/ 50195 - 668X ( 03 ) 00201 - X .