Open Access

Recognition of lung volume condition based on phase space mapping using electrical impedance tomography


Cite

Introduction

A great number of people die annually due to various kinds of lung diseases, which can cause impariment and respiratory disorders that are mostly caused by air pollution, tobacco consumption, genetic tribulation, etc. Lungs are the organs of the body which expand and contract hundreds of times durung the day separating carbon dioxide from oxygen. Their disorders can cause serious problems in the whole respiratory system [1]. These pulmonary diseases are namely: Asthma, chronic obstructive pulmonary diseases (COPD), pneumonia, tuberculosis, lung cancer, etc. Asthma is an illness which can be detected by cough and origins from the inflammation of the respiratory tract. It is mostly caused by external factors such as dust, the breath of animals or some specific plants. This illness appears by the reaction of the body’s immune system [2,3].

Chronic obstructive pulmonary diseases (COPD) consists of bronchitis and emphysema. Chronic bronchitis create sputum, however, emphysema may cause other kinds of lung disorders. The majority of patients usually suffer from both [4]. Chronic bronchitis is created by a mucus overproduction of the lung cells, and it declines the lung function, which increases the risk of airway obstructions, raises the risk of respiratory infection and could even cause death [5].

Pneumonia is an infectious diseases among adults in which the bacteria enters to the body from the individual’s respiratory tracts [6].

Tuberculosis is another bacterial and also contagious disease which spreads through blood or lymph node [7].

Furthermore, lung cancer is an illness which is contained from myriad forms. The most prevalent of those takes place in the center of the lung lobes [8].

Electrical Impedance Tomography (EIT) is a non-invasive imaging modality which possesess a lot of potential in industry and the medical field. This method applies electrical stimulation and related measurements at the body surface to extract the electrical characteristics of internal tissues.

The distribution of electrical conductivity within a tissue could be captured by applying an electrical current to the tissue and measuring the surface’s voltage divisions [9].

Furthermore, EIT can be applied bedside as a continuous monitoring technique, as it only requires small devices and does not expose patients or users to any ionized radiation. The EIT technique was introduced in the early 1980s by Barber and Brown. Soon after, a broad spectrum of possible applications in medicine was suggested, ranging from gastric emptying to brain function monitoring and from breast imaging to lung function assessment. This method was tested on acute respiratory distress syndrome by measuring surface voltages frequently using electrodes [10, 11].

In electrical impedance tomography, current is injected from electrodes placed arround the chest, then the resulting surface voltages are measured [12]. In 1985, clinical diagnostic software was provided using EIT, which is a convenient method for lung imaging [13]. In 1999, by using the EIT method it was concluded that the electrical properties for each section of a tissue is different from the others [14]. According to research by Victorino and his colleagues, measured imbalance in air ventilation in each different region of the lungs was due to pulmonary damages. By comparing EIT and CT methods, the EIT was introduced as a great procedure for diagnosis of lungs disorders [15]. EIT images were identified as well as CT images and as a result, EIT had the advantages of both being non-invasive and free of radiation. [17, 18]. In 2016, research was done on a number of anesthetized rats and EIT images of cortical neural activities were obtained [19]. Additionlly, some studies were assessed in the fields of the various patterns of the EIT measurement to maximize the heart rate signal in parallel with lungs by Tobias Menden and his colleague in 2017. They found that the heart rate signal was completely different from the respiration signal [20]. In 2017 Chitturi and his colleague focused on important techniques to upgrade EIT's spatial system [21]. Hanan and her colleagues have performed another research and found out an essential tool for understanding neural circuits by using EIT [22].

We used EIT signal as a non-invasive method for evaluating lung volumes and comparing the differences between normal and abnormal lungs with mathematic techniques such as Phase Space methods.

Materials and methods
Data base

Data used in this paper includes 82 registered signals recorded from 52 patients with lung disorders collected by Dr. Hashemian from the Masih Daneshvari Hospital of Tehran. The sampling frequency is 20 Hz and data are mostely recorded from patients who have had a surgery and were using artificial respiration. For recording the data, a belt with 16 electrodes attached to the device was utilized. This belt have got the ability of adjustment according to the size of an adult’s body and closes below the patient’s breast, right around the lung. The outputs consist of electrical impedance signals of four regions of lungs according to fig. 1. These regions are numbered from 1 to 4. A sample of registrated signals from each region using the EIT system, is depicted in fig. 2. The total of four region’s signals (ROI1 to ROI4) make up the Global signal.

Fig. 1

Lungs Classification to four regions.

Fig. 2

Electrical Impedance Tomography of lung’s four regions.

Suggested method

In this paper, the state of lung volume after surgery has been considered before non-linear signal processing. Comparing the results can make a significant contribution for doctors and aiding them to a better detection of lungs function. Looking further, since EIT signals with time and space accuracy cannot be detected by lung’s volume which was obtained before and after surgery, a better recognition of lung volume condition named Phase Space Method was used to introduce an indicator for comparing the states of the lung’s volume throughout and after the surgery.

Phase space

Recent studies about medical signal processing illustrate that due to the nature of the non-linearity of the EIT signal, linear methods do not succeed in extracting useful information from the signal and cannot classify with high accuracy between normal signals and patient’s signals. Thus, it can be observed that today, in various studies the phase space method is used to diagnose among normal signals and patient signals [23, 24, 25].

Assuming that the input signal of s(n) = s1, s2, …, sn has been defined for making phase space, the s(n) signal was divided into two signals called x1 and x2. In this formula τ is the delay parameter.

x1:s2.s3.snx2:s1.s2.snτ$$\begin{align} & {{x}_{1}}:{{s}_{2}}.\,{{s}_{3}}.\,\,\ldots \,\,{{s}_{n}} \\ & {{x}_{2}}:{{s}_{1}}.\,{{s}_{2}}.\,\,\ldots \,\,{{s}_{n-\tau }} \\ \end{align}$$

Next, for creating the phase space, x1, x2 and τ = 1 are drawn in accordance with each other. τ = 1 means that every sample is possessed in accordance with the next sample of itself. For making a distinction, different phase spaces must be introduced to the features of this mapping. In this paper, deviation criteria for short-term and long-term changes are namely: SD1 (Standard Deviation) and SD2 which are the features of this mapping. SD1 and SD2 are dispersions of points perpendicular to the axis of line-of-identity. Formulas(2) show how to calculate SD1 and SD2:

SD1=22std(x1x2)SD2=2std(x1)20.5std(x1x2)2$$\begin{matrix} S{{D}_{1}}=\frac{\sqrt{2}}{2}std\left( {{x}_{1}}-{{x}_{2}} \right) \\ S{{D}_{2}}=\sqrt{2std{{\left( {{x}_{1}} \right)}^{2}}-0.5std{{\left( {{x}_{1}}-{{x}_{2}} \right)}^{2}}} \\ \end{matrix}$$

Figure 3 shows the phase space for one EIT signal sample and parameters introduced in this mapping.

Fig. 3

Phase space mapping.

Informed consent

Informed consent has been obtained from all individuals included in this study.

Ethical approval

The research related to human use has been complied with all relevant national regulations, institutional policies and in accordance with the tenets of the Helsinki Declaration, and has been approved by the authors’ institutional review board or equivalent committee.

Results

As mentioned before, the temporal data of volume signals cannot be well differentiated by the arrival of lung volume to the desired state. In this paper, the idea of applying the nonlinear method and converting the lung volume signal into phase space maps was presented. Figures 4 and 5 delineate the phase space mapping of a suitable and an unsuitable Global Signal respectively. As it can be observed, no information can be extracted to understand the condition of lungs’ volume by visual comparison of Figs. 3 and 4.

Fig. 4

A phase space mapping for suitable lung volume.

Fig. 5

A phase space mapping for unsuitable lung volume.

Furthermore, the calculation of SD1 and SD2 parameters of two modes (suitable and unsuitable lung volume) have also showed similar results. Table 1 presents the average results and the standard deviation of the SD1 and SD2 features in the case which lung volume is present in two appropriate and inappropriate situations. The statistical analysis were used to distinguish between two groups (before and after surgery). And as a result, the value of destination between the two groups with P - value was considered less than 0.05. All the analyses were obtained using MATLAB 2018.

The average and the standard deviation of the SD1 and SD2 parameters with τ = 1.

Lung Volume ConditionsSD1SD2
Mean ± Standard deviationMean ± Standard deviation
Unacceptable89.23 ± 10.451914.71 ± 203.61
Suitable94.43 ± 11.7611673.52± 174.11
Implementation of the delay parameter in phase space mapping

The Phase Space Mapping introduced in the previous section is plotted with the delay parameter value (τ), equal to 1. For a better discrimination between two groups, Phase Space Mapping was reconstructed by applying different delay parameters. Values in Table 2 indicate the mean and standard deviations of the SD1 and SD2 criterion for the delay parameter of 2 to 6. The statistical comparison of the results denotes the meaningful distinction of the SD1 property with increasment of τ. So that when the volume is abnormal and the amount of τ is considered 6, the average value of SD1 is twice that of the normal case (P-value < 0.05). This means that the EIT signal in the case that a person's lung is in a suitable condition shows a higher uniformity in the four regions.

The average and the standard deviation of the SD1 and SD2 parameters with τ = 2.

Lung Volume ConditionsSD1SD2
Mean ± Standard deviationMean ± Standard deviation
Unacceptable96.51 ± 12.102104.23 ± 223.25
Suitable112.67 ± 24.091716.06± 213.43

The average and the standard deviation of the SD1 and SD2 parameters with τ = 3.

Lung Volume ConditionsSD1SD2
Mean ± Standard deviationMean ± Standard deviation
Unacceptable118.19 ± 15.752217.25 ± 228.51
Suitable158.09 ± 34.121835.12± 216.41

The average and the standard deviation of the SD1 and SD2 parameters with τ = 4.

Lung Volume ConditionsSD1SD2
Mean ± Standard deviationMean ± Standard deviation
Unacceptable132.54 ± 21/142375.78 ± 234.01
Suitable216.37 ± 41/621932.64± 220.17

The average and the standard deviation of the SD1 and SD2 parameters with τ = 5.

Lung Volume ConditionsSD1SD2
Mean ± Standard deviationMean ± Standard deviation
Unacceptable148.20 ± 24.172418.23 ± 239.67
Suitable297.41 ± 49.252098.27 ± 226.15

The average and the standard deviation of the SD1 and SD2 parameters with τ = 6.

Lung Volume ConditionsSD1SD2
Mean ± Standard deviationMean ± Standard deviation
Unacceptable156.71 ± 26.012457.98 ± 245.54
Suitable342.32 ± 57.752113.18 ± 230.67

Figures 6 and 7 illustrate feature changes of the SD1 and SD2 with a leap in τ. The comparison of the obtained results shows that by increasing the delay parameter (τ), the value of the SD1 parameter rises, and when the amount of τ remained stable at 6, there is a significant shift between appropriate (342.57 ± 32.75) and inappropriate (156.71 ± 26.01) lung volume status. It is to be noted that the SD2 parameter does not show significant changes in the two desirable and undesirable lung conditions as τ rises.

Fig 6

Comparison of the SD1 parameter for different value of the τ.

Fig. 7

Comparison of the SD2 parameter for different value of the τ.

Discussion

In this paper, we use the method of Phase Space Mapping to compare the lung volume situations in different patients. With the help of the diagnostic computer system presented in this paper, a better understanding of the condition of the lung volume after surgery can be obtained. The results of the proposed method can provide valuable information about the function and the amount of patients’ lung volumes for the information of doctors. In Gonzales’s research in 2016, the use of 3D television systems was significantly more effective in the 3D EIT imaging than the Gaussian smoothing rule, which can be used in imaging of electrical tomography and optical diffusion tomography [9]. In a study in 2015, EIT, which controls lung’s condition over time, is a significant step towards detection of lung diseases and can demonstrate patient conditions better than any other imaging method [10]. In 2012, Smith decided that in diseases like ARDS, patients are in serious need of artificial respiration with specific conditions (ECMO), for this very reason, it is highly significant for such patients’ lungs’ condition to be exhibited continuously in order to inform doctors about the changes that need to be applied in the ventilator device settings such as (HFOV) [11]. According to the data collected from the pulmonary investigation in 2000, the smallest changes in the adsorbed oxygen are visible in detail [13]. In a paper from 1999, which was authored by Cheney and his colleagues, the mathematical model relations are defined for the EIT to be more accurate in some existing theories referred to by the adaptive current tomography (ACT), and after studies, an algorithm was selected that was the most effective one for making images [14].

In 2004, according to comparisons between dynamic computerized tomography in a heterogeneous population of critically ill patients under mechanical ventilation in different moments, Victorino and his colleagues concluded that, among variables derived from computerized tomography, changes in absolute air content was best explained by the integral of impedance changes inside the regions of interest [15]. Pelosi and his colleagues also carried out a study in which each lung was divided into 10 X-rays, and the hydrostatic pressure (SPL) was calculated separately, and the result was that the hydrostatic pressure in the lungs of ARDS patients is higher than the hydrostatic pressure in the lungs of healthy individuals [16]. In 2000, Newman and his colleagues tested some EIT experiments on 6 pigs to find similarities between them and their ARDS patients. It was concluded that in order to prevent lung collapse during artificial respiration, PEEP should be higher than 20 cm of water in less than 0.6 seconds while a PEEP with maximum of 15 cm of water in 1.4 seconds is necessary for these patients. Therefore, due to the differences in time, it cannot be stated that there are any similarities between patients and the pigs [17]. In another paper published in 2016, different tests have been done on brain systems of two mice in which, only one of them being injected in their brain to investigate the differences, and as a result, no changes were observed in any of them [18].

In 2017, Tobias and his colleagues compared the signals received from hearts and lungs, and the result was that the frequency spectrum of the heart had various zig-zag patterns, at lower measurement frequencies there were smaller patterns and vice versa [20]. In 2017, focusing on the previous research, the result was found that in a breast cancer case peak-to-peak distances represent disturbances of the breast, and Array Processing technique can be a step towards the development of spatial resolution for future EIT applications [21]. In 2018, Hanan and his colleagues have found two frequency responses in two impedance signals associated with epilepsy in mouse brain cortex, which is due to the ability of imaging the neural activity by this method of imaging as well as the metabolic changes [22].

Therefore, by applying the proposed method in this paper, we used informations about the individuals’ lung volumes with chronic respiratory diseases and accurate assessmentss of the different methods to treat respiratory system disorders. This method can reduce the patient’s desired treatment duration in addition to saving various therapeutic costs, which might be caused by using poor treatment methods that are also likely to give adverse consequences.

Conflict of Interest

The authors whose names are listed in this manuscript has no affiliation with or involvement in any organization or entity with any financial interest or nonfinancial interest in the subject matter or materials discussed in this paper. Authors have also declared that no conflict of interest exists.