Open AccessArticle

A Prototype Framework Design for Assisting the Detection of Atrial Fibrillation Using a Generic Low-Cost Biomedical Sensor

Jesús Pérez-Valero

¹,

Antonio-Javier Garcia-Sanchez

^1,*

Manuel Ruiz Marín

and

Joan Garcia-Haro

Department of Information and Communication Technologies, Universidad Politécnica de Cartagena (UPCT), Campus Muralla del Mar, E-30202 Cartagena, Spain

Department of Quantitative Methods, Law and Modern Languages, Universidad Politécnica de Cartagena (UPCT), E-30202 Cartagena, Spain

Author to whom correspondence should be addressed.

Sensors 2020, 20(3), 896; https://doi.org/10.3390/s20030896

Submission received: 10 December 2019 / Revised: 24 January 2020 / Accepted: 4 February 2020 / Published: 7 February 2020

(This article belongs to the Special Issue Advances in Wireless Sensor Networks Applications for Health and Environment)

Download

Browse Figures

Figure 1
SRPs of a 102-length RR interval time series for an embedding dimension m = 3 in a normal sinus (a) and atrial fibrillation (b) patient. "> Figure 2
Block diagram showing the main modules that shape the prototype application developed. After the entire procedure, the objective is to determine if a patient is with normal sinus rhythm (NS) or with atrial fibrillation (AF). "> Figure 3
Block diagram showing the steps of the digitalization procedure. "> Figure 4
Horizontal scanning (green line) vs vertical scanning (red line). Arrows represent the direction of the scanning process. "> Figure 5
Transformation process of the signal from PDF or JPEG formats (a), then a binary mask image is obtained (b) and after that, some pixels are removed (c). Finally, a shrink operation is applied to capture only useful information (d), and then the output digitalized signal (e) is shown. "> Figure 6
Transformation of the signal before (a) and after (b) applying the filter function "> Figure 7
Kardia Mobile device by AliveCorc (a) and MAC800 by General Electric Healthcare (b); both employed for the calibration procedure. "> Figure 8
Electrocardiograph (ECG) placement configuration in both devices MAC800 (a) and Kardia Mobile (b) "> Figure 9
Different ECG deviations from (a–d) of the Kardia Mobile (blue) with respect to the MAC800 (red) equipment. "> Figure 10
RR interval data for 20 samples from 20 different patients measured with Kardia Mobile (blue), MAC800 (red) and corrected Kardia Mobile (green). "> Figure A1
Functional diagram of the application implemented in MATLAB (a). Remote access to the developed application in the cloud (b). ">

Versions Notes

Abstract

Cardiovascular diseases are the leading cause of death around the world. As a result, low-cost biomedical sensors have been gaining importance in business and research over the last few decades. Their main benefits include their small size, light weight, portability and low power consumption. Despite these advantages, they are not generally used for clinical monitoring mainly because of their low accuracy in data acquisition. In this emerging technological context, this paper contributes by discussing a methodology to help practitioners build a prototype framework based on a low-cost commercial sensor. The resulting application consists of four modules; namely, a digitalization module whose input is an electrocardiograph signal in portable document format (PDF) or joint photographic expert group format (JPEG), a module to further process and filter the digitalized signal, a selectable data calibration module and, finally, a module implementing a classification algorithm to distinguish between individuals with normal sinus rhythms and those with atrial fibrillation. This last module employs our recently published symbolic recurrence quantification analysis (SRQA) algorithm on a time series of RR intervals. Moreover, we show that the algorithm applies to any biomedical low-cost sensor, achieving good results without requiring any calibration of the raw data acquired. In addition, it has been validated with a well-accepted public electrocardiograph (ECG) data base, obtaining 87.65%, 91.84%, and 91.31% in terms of sensitivity, specificity and accuracy, respectively.

Keywords:

biomedical sensor; RR interval; R peak, ECG; logistic regression

1. Introduction

Over the last few decades, life expectancy has increased considerably due to better quality in living standards and constant progress in medicine, among other causes [1]. The consequential increase in population has significant social impact in terms of healthcare requirements. In this regard, healthcare cost containment is a crucial factor which implies strict control of medicines, medical equipment, and hospital care expenditures [2]. There is, thus, a need to develop and implement groundbreaking devices able to provide medical services at an affordable cost. This is particularly true for populations with limited access to basic healthcare services.

Cardiovascular diseases (CVDs) are the leading cause of death in the world [3]. In 2016, over 17.9 million people died of CVDs, representing 31% of all global deaths by disease. Of these deaths, 85% are due to heart attacks and strokes, and about three quarters take place in countries with low to middle average incomes [4].

To address these troubling statistics, a possible approach consists of thoroughly studying the functioning of patients’ hearts using portable devices. They can gather an enormous amount of data through recurrent heart tests which can be carried out at home, while patients sleep or go about their daily activities. So, monitoring patients with portable sensors is a way to improve life quality; especially for elder individuals and people with mobility problems, or to extend access to a health protection system. This type of monitoring provides an efficient and cost-effective alternative to on-site clinical monitoring [5]. Systems equipped with low-cost portable sensors are useful and feasible diagnostic tools for the real-time monitoring of important physiological indicators in such domestic/external facilities [6]. In the same context, wearable sensors are gaining popularity and attracting the attention of many researchers and technological companies [7]. Consequently, many leading corporations such as Apple and Samsung are spending large amounts of money on developing wearable devices focused on heart fitness. The set of features integrated in a single wrist-worn device will provide relevant information about the individual in almost real time. The main drawback to this approach is its low accuracy in ECG data acquisition. Interestingly, our methodology does not require ECG data calibration to provide good results.

The problem of data quality in ECG sensors has also been tackled in [8], where authors developed a single-chip based wearable ECG hardware system for monitoring patients and offer a high-quality ECG signal by applying three types of filters. Nevertheless, they did not focus on a heart disease classification approach. The authors of Reference [9] proposed a multi-sensor data fusion scenario in order to improve the quality of heart disease detection. They used a kernel random forest ensemble with several time and frequency domain features for their classification method. Even though they attained 98% accuracy working on a continuous stream of data, they incur in network bandwidth and battery consumption. Contrastingly, our design is single-device oriented, uses only RR intervals (time between two consecutive heartbeats) as input, and requires no network configuration.

Parallel to the development of these portable sensors for health monitoring, a large set of algorithms has been implemented to understand, detect, and classify different diseases. For instance, during the last decade, the design of different algorithms for the analysis of electrocardiogram signals has attracted the interest of many researchers. The detection and classification of different heart arrhythmias, in particular atrial fibrillation (AF), has played a central role since AF is a leading preventable cause of recurrent stroke for which early detection and treatment are critical.

Some relevant algorithms for the detection of AF are the ones based on the analysis of the P-wave [10,11], or those that utilize a large set of features obtained from ECGs in an artificial neural network [12,13,14] or in a deep learning approach [15,16]. Furthermore, it is worth mentioning that several recent clinical trials on large populations focusing on AF detection have been carried out [17,18,19]. In those clinical studies ECG data were obtained with portable sensors achieving outstanding results.

However, some limitations of these algorithms are: the high dependence on the robustness of the training data, the requirement for a filter preprocessing phase, and the fact that they are built upon a large set of features that are usually computationally complex to analyze, among others.

In order to overcome some of these limitations, several algorithms for the automatic detection of AF, based only on RR intervals (time between two consecutive heart beats), have been introduced [20,21,22].

Our work is aimed at developing a generic application and its associated methodology to be able to apply a predictive algorithm over the RR intervals acquired from any low cost biomedical sensor. In principle, this will not be easy to carry out due to the expected low precision of most of the portable sensors on the market. Portable sensors or wearables for this purpose, such as Apple Watch Series 4 [23], MySignals [24] or Qardio [25], provide electrocardiogram data in PDF format, which could be used as the input to our application. Regarding the low accuracy of these devices, it would be necessary to figure out a calibration procedure. In order to rapidly create a fully functional prototype, and taking low-cost and publicly available data into consideration, we use a commercially available portable electrocardiograph sensor called Kardia Mobile from AliveCor [26]. Albeit, Kardia Mobile has some limitations recently published [27,28] in comparison to the data accuracy of a 12-lead gold standard ECG recording. The MAC800 Healthcare [29] is employed as the 12-lead gold ground truth standard ECG equipment. For reproducibility purposes, the source code of the application developed, as well as its corresponding updates, can be found in Reference [30]. Furthermore, a possible MATLAB functional diagram is shown in Appendix A.

In addition, in this paper, we use our recently proposed predictive algorithm to detect atrial fibrillation, which is based on Symbolic Recurrence Quantification Analysis and can be found in Reference [31]. In order to maintain a self-contained paper, its more important features will also be summarized here. The use of Symbolic Recurrence Quantification Analysis (SRQA) [32] will provide us with a framework to analyze dynamic changes in a time series. This time series will be extracted from electrocardiogram (ECG) signals. More concretely, RR intervals are the measurement processed by our algorithm. The use of RR intervals as the unique input for our algorithm reduces the complexity of the approach and allows us to create a robust, simple and reliable solution for short and long term ECG monitoring. Specifically, we constructed a logistic regression algorithm that employs SRQA over RR interval time series to detect AF. The resulting approach can differentiate between normal sinus individuals and patients with AF. Furthermore, in order to show the power of our algorithm, we performed a logistic regression model that provides outstanding results in terms of power accuracy.

Our main contributions are: (i) the integration in an application of our computing time efficient algorithm; (ii) the fact that our algorithm only employs RR intervals as input so it can virtually operate with any low-cost sensor; and (iii) predictive power accuracy not being affected by a data calibration procedure.

The rest of this paper is organized as follows. Section 2 outlines the symbolic dynamics approach as well as the logistic model for the classification algorithm. Section 3 elaborates the methodology for applying the algorithm to generic low-cost devices. This section also presents the digitalization phases, the determination of RR intervals, the selectable calibration procedure, and the devices employed for such purpose. Section 4 describes the data employed for training the calibration method and the data used for the validation of the predictive algorithm to detect AF. Finally, Section 5 concludes.

2. Normal Sinus and Atrial Fibrillation Classification Scheme

Recently, we have consistently applied symbolic recurrence quantification analysis (SQRA) for the detection of atrial fibrillation [31]. SRQA provides a powerful framework to study the dynamic behavior of a system. In particular, we developed a novel algorithm based on a logistic model that we called ReAD-AF (Recurrence Analysis to Detect Atrial Fibrilation). The algorithm is intended to determine whether a patient is in a normal state (normal sinus rhythm NS) or with AF. The solution is based on the RR interval of an ECG signal as the unique input, requiring minimum pre-processing computational cost. Specifically, we constructed a logistic regression algorithm that employs SRQA over RR interval time series. Here, we summarize its operation and main features. Additionally, the pseudocode implementing the algorithm is included in Appendix B. The proposed approach uses a symbolization procedure based on neighboring permutations. That is, given an RR interval time series

{x_{t}}_{t = 1}^{T}

, its space representation can be built by means of Taken’s time delay method [33], so the space vector will be

{\bar{x}}_{t} = (x_{t}, x_{t + 1}, \dots, x_{t + m - 1})

for an embedding dimension m. A symmetric group (

S_{m}

) of order m! is defined as a group that contains all the permutations of length m. In turn, each permutation is defined as a symbol

π

. It should be mentioned that each vectorial time series from

{\bar{x}}_{t}

is called m-history. Finally, a symbolization map

S ({\bar{x}}_{t})

transforms the m-history of the vectorial time series into a sequence of symbols

π

representing ordinal patterns of length m, so that

{\bar{x}}_{t} \mapsto S ({\bar{x}}_{t}) = π .

(1)

As an example, the transformation of the RR interval sample set (2) into a sequence of symbols (3) is shown as follows.

x_{t} = (0.55, 0.65, 0.60, 0.70, 0.65, 0.50) .

(2)

Taking into account that if

m = 3

, the symmetric group is formed of

m! = 3! = 6

symbols, so

S_{3} = \{(2, 1, 0), (2, 0, 1), (1, 2, 0), (1, 0, 2), (0, 1, 2), (0, 2, 1)\} .

(3)

Each m-history of length 3 corresponds to a symbol (and an assigned color) of the symmetric group

S_{3}

. Thus, for

t = 1

{\bar{x}}_{1} = (0.55, 0.65, 0.60)

, it can be observed that

x_{t + 0} = 0.55 < x_{t + 2} = 0.60 < x_{t + 1} = 0.65

which means that

{\bar{x}}_{1}

corresponds to a symbol of type (0,2,1). Similarly, for

t = 2

{\bar{x}}_{2} = (0.65, 0.60, 0.70)

, it is observed that

x_{t + 1} = 0.60 < x_{t + 0} = 0.65 < x_{t + 2} = 0.70

, so

{\bar{x}}_{2}

corresponds to a symbol of type (1,0,2).

Now, two symbols

S ({\bar{x}}_{t})

and

S ({\bar{x}}_{s})

at different instants in time (t and s) are recurrent if and only if

S ({\bar{x}}_{t})

S ({\bar{x}}_{s})

. So, an indicator function is defined as

S R P_{t s} = \{\begin{matrix} k if S ({\bar{x}}_{t}) = S ({\bar{x}}_{s}) = π_{k}, \\ 0 otherwise . \end{matrix}

(4)

Therefore, the Symbolic Recurrence Plot (SRP) is defined as a matrix that contains the recurrences of all the symbols provided by the indicator function. The SRP becomes a powerful tool for analyzing the dynamic behavior of the time series and can be represented in a graph where the axes represent the time indexes of the time series. Each colored dot at coordinates (t,s) in the graph means that the m-histories

{\bar{x}}_{t}

and

{\bar{x}}_{s}

are recurrent to the corresponding symbol. White dots mean that the m-histories are not recurrent. This enables the display of the portion of the phase space that is being visited by the system. If a dynamic change occurs, then the colored distribution of the SRP will also change.

Figure 1 illustrates an example of two SRPs plotted from a 102-lengh RR interval time series from a patient with normal sinus rhythm and another with atrial fibrillation. We can observe that the colored distribution is different in both graphics, which means that the dynamic behavior changes in both time series. The normal sinus rhythm SRP presents a much more structured distribution, with a remarkable predominant color (black) of symbolic recurrence. In contrast, the atrial fibrillation SRP does not exhibit a structured pattern. The symbols appear in a random way without following any particular distribution. These random patterns characterize the atrial fibrillation patient due to the non-periodic state of the RR intervals. However, the normal sinus rhythm patient presents a characteristic set of symbols that repeat periodically.

To figure out the probability of a patient being categorized as being in a NS or with AF, a logistic model is created. To estimate the model, it is necessary to define a set of covariates. These covariates are divided into two groups. The first group contains those based on the diagonal and vertical lines of the SRP, as well as the recurrence rates of each symbol. The second group is based on the distribution of the RR intervals and comprises the mean, median, Pearson coefficient of variation, and the coefficient of variation of the median. These covariates allow the discrimination between normal sinus rhythm and atrial fibrillation patients. The logistic model is generated through a receiver operating characteristic (ROC) curve analysis in order to compute a probability threshold. So, a patient with an estimated probability above this threshold is classified as having AF.

In order to build this curve, three metrics are assessed. Namely, sensitivity, specificity, and accuracy. To this end, the number of true positive (TP), true negative (TN), false negative (FN), false positive (FP), sensitivity (Se) or the true positive rate

S e = \frac{T P}{T P + F N}

, and the false positive rate

F P R

given by

F P R = \frac{F P}{T N + F P}

have to be defined. In addition, specificity (Sp) is also defined as

S p = 1 - F P R

. Finally, accuracy (ACC) is determined as

A C C = \frac{T P + T N}{T P + F N + F P + T N}

. Taking all these metrics into account, the ROC curve is constructed using the

(F P R, S e)

points, and the optimal threshold

τ

parameter is the one that minimizes the distance between the points (FPR, Se) and (0,1),

τ = \underset{τ}{arg min} {F P R_{τ}^{2} + {(1 - S e_{τ)})}^{2}} .

(5)

In the original paper [31] it was shown that sensitivity, specificity and accuracy increased with the window size selected. Sensitivity always reached values above

96 %

and specificity slightly lower with values above

94.8 %

for the smallest window sizes. Finally, accuracy achieved values of

95.4 %

presenting good predictive power.

To validate the algorithm, a k-fold cross-validation procedure was performed. The logistic model was fitted with a training set and further evaluated with a test set. The entire data set is publicly available and provided by Physiobank [34].

The ReAD-AF algorithm was shown to be robust and capable of discriminating AF patients with high precision. Indeed, the ReAD-AF algorithm developed can work with only the RR interval feature as input, which implies a reduction of the computational cost of the model, enabling portability to mobile platforms such as smartphones.

In addition to this classification algorithm based only on RR intervals as input, we want to further contribute with a methodology to help practitioners develop an application based on algorithms that can be applied to the growing low-cost biomedical sensor field. Thus, the next section illustrates a framework for such purpose.

3. Prototype Application

In this section we look deeper into the application modules developed. A general block diagram of this application is depicted in Figure 2. A detailed description of these main modules as well as how they work is discussed below.

The first module uses an ECG signal (acquired from a low-cost device) as input in PDF or JPEG format and digitalizes it. Satisfying the well-known sampling theorem, the module takes as its main parameter double the bandwidth of the device under consideration. Selecting double the bandwidth is enough to be able to obtain all the signal information. The background color of the electrocardiogram and the duration of the ECG are also important parameters to take into account. Since direct access to data captured by a particular sensor usually follows proprietary protocols, this digitalization procedure enables the implementation of a generic application, regardless of the low-cost commercial device finally selected to provide the input ECG.

The next module is responsible for detecting the inter-beat (RR) intervals also known as heart rate variability (HRV) in an ECG. This metric measures the specific changes in time (or variability) between successive heartbeats. This time variation (usually expressed in milliseconds or seconds) is controlled by a primitive part of the nervous system called the autonomic nervous system (ANS). RR intervals are one of the indicators of a person’s state of health, fitness, and recovery. There are many algorithms available in the open literature for detecting heart rate variability (HRV). For instance, the well-known Pan and Tompkins algorithm [35] detects QRS waves, or QRS complex, based on the slope, amplitude, and width of the ECG signal. Our algorithm is divided into two stages; the pre-processing and the decision-making process. In the pre-processing phase, the signal is prepared by removing noise, smoothing the signal, and amplifying the QRS slope. Then, in the decision stage only signal peaks are obtained. It is also possible to apply a set of well-proven tools developed by MIT [36,37], enabling the detection of QRS complex and RR intervals. In our prototype, the implementation of the RR interval method is carried out with a native MATLAB function (findpeaks).

Once the RR intervals have been acquired, our application could take them directly to further feed a third module for calibrating these data obtained by a low-cost portable device against professional 12-lead gold standard electrocardiograph equipment in order to achieve more precise data and later diagnosis. This calibration is usually based on a linear regression model aimed at improving the accuracy of the data acquired by the low-cost device finally employed. The calibration will allow us to analyze the ECG signal with the highest degree of accuracy. Interestingly, we prove that measures obtained from SRQA are invariant under monotonic transformations (see Appendix C), which means that the symbolic recurrence measures employed as covariates in our algorithm are not affected in those cases in which the calibration processes are based on linear regression. That is, raw uncalibrated data can be used to feed our algorithm directly. However, for practitioners who might use different approaches, a possible calibration method will also be introduced below.

The last module, implementing the ReAD-AF scheme, is in charge of assisting the AF classification. Its power of classification is addressed in Section 4.2.

All of these modules are intended to illustrate a general methodology since, as commented above, commercial sensor applications are mainly based on proprietary protocols. That is, the data acquired is not directly accessible, and only an ECG representation, commonly in PDF or JPEG format, is available.

3.1. Module 1: ECG Digitization

To analyze the ECG signals comprehensively, we first require a technique that enables the digitalization of a graph in PDF or JPEG formats, which are the most common output formats of commercially available sensors. This will provide independence to the data acquisition platform and simplicity in the generalization of our methodology. For this purpose, we have customized a procedure implementing this functionality [38]. The sequence of actions for the digitization process of the ECG signal in PDF or JPEG format are shown in Figure 3, in addition to a detailed description of each one.

The first phase of our proposal is the scanning of the ECG signal. To do this, a MATLAB-scan-function (imread) is applied at a resolution of 600 dpi which takes a PDF file containing the ECG as an input parameter.

To perform the scan a portion of the PDF document must first be selected. The complete PDF document contains additional information, such as the patient’s name, or his/her date of birth, among other things. Therefore, it is necessary to prune the signal to be digitalized.

In the image binarization phase, RGB color images are converted into binary images. This process reduces the number of colors at the binary level, resulting in a clear reduction of the memory needed to store the signal. This function also simplifies the processing of the binarized image in comparison to the original image. However, ECG charts are always represented graphically in a grid, which interferes with the binarization process of the image. One way to eliminate or erase this grid in the background of the image is to set an RGB threshold. In detail, all interfering colors are enclosed in the threshold and transformed into white before starting the binarization mechanism. The end result is a white background with the ECG signal in a selectable (e.g., blue) color.

Some noise might appear, characterized by isolated black and white pixels in the binarized image. This may impair the subsequent analysis of the data. Therefore, a filtering function has also been implemented to eliminate these error pixels. To extract the ECG signal graph from the two-dimensional binarized image, horizontal scanning is usually employed. However, our module implements vertical scanning to identify the ECG signal pixels. The vertical scanning is selected to achieve faster and more accurate digitalization. While with horizontal scanning pixel location identification requires an iterative process, with vertical scanning, each pixel directly represents a (x, y) position coordinate of the graph, which results in a one-dimensional vector. Figure 4, below, illustrates both types of scanning.

In order to offer greater data accuracy, vectors double the number of columns of pixels analyzed, because some peaks are plotted as successive vertical black lines in scenarios with low resolution images. Each component of these vectors is a complex number. The two components of the same column have an identical real part, which is the column index of the pixel matrix. The two consecutive imaginary parts (of the same real part) quantify the upper and lower limits of the vertical black line, that is, of the ECG signal under study. For example, the 2D-vector

T = [\dots; 12 + 27 i; 12 + 28 i; \dots]

means that the vertical black line spans from line 27 to line 28 in column 12.

The digitalization algorithm searches the data from the first pixel (bottom left) to the last one (top right). Once the original axes have been deleted from the ECG graph, the value of the bottom vertical and the first column positions are used as a reference for the ECG scanning.

To convert the above mentioned 2D vector into a one-dimensional vector, the algorithm computes the modal distance (

μ

) between the imaginary part of two successive components; that is, the number of vertical pixels comprising the signal in a specific column. If the difference between the imaginary components of the two coordinates is within the limit

μ

, then only the imaginary part of the second element is included in the 1D version of the vector. This value is assumed as a reference (

λ

) to the R peak identification in the next column analysis. If this difference is greater than

μ

, the module of the difference between the

λ

and each of the two components is calculated. The stored 1D value is the greatest obtained value. One can often find portions of the graph where the drawing is no longer continuous. In order to provide a consistent 1D vector, the components in these positions are estimated through linear interpolation. Finally, the digital signal is plotted taking into account this final one dimensional vector and the sampling frequency.

For a better understanding of the digitalization process, the pseudocode of the employed methods is also shown below in Algorithm 1, as well as an explanation of the sequence of actions. The input is an ECG in PDF or JPEG format. Then, step 2 allows the reading of the input file. Step 3 applies an interactive image pruning operation to the image previously read. After this, the original image is transformed into a binary one. This is denoted in step 4, during which a binary mask and a composite image are generated. The composite image shows the original RGB pixel values under the mask. Step 5 indicates that portions with fewer than P pixels, in our case 1260 pixels, are removed from the binary image, resulting in a clearer signal where all the pixels but the ECG signal curve are eliminated. Then, a shrink function is applied in step 6 for a better union between points. Finally, vertical scanning is implemented in step 7 for the remaining steps to convert the binarized image into a 1D vector, and the resulting output signal is then plotted in MATLAB format (.fig).

Algorithm 1 PDF to MATLAB

1:: Input: PDF file (.pdf)
2:: $r \leftarrow$ read(PDF)
3:: $p \leftarrow$ prune(r)
4:: $b \leftarrow$ createBinaryMask(p)
5:: $b w \leftarrow$ bwAreaOpen( $b, 1260$ )
6:: $m \leftarrow$ bwMorph( $b w, s h r i n k, 2.5$ )
7:: $v e c t o r \leftarrow$ pixelToVector(m)
8:: Output: MATLAB signal (.fig)

Figure 5 shows the digitalization procedure step by step. At each step, the image is carried out by a transformation to be able to obtain the final digitalized signal.

3.2. Module 2: Signal Processing and Filtering

Once the digital output signal is obtained, the processing phase to attain the RR intervals starts. First, in order to smooth the data obtained from Module 1 in Section 3.1 we apply the MATLAB function filter(b,a,x) that uses the rational transfer function

Y (z) = \frac{b (1) + b (2) z^{- 1} + \dots + b (n_{b} + 1) z^{- n_{b}}}{a (1) + a (2) z^{- 1} + \dots + a (n_{a} + 1) z^{- n_{a}}} X (z),

(6)

with

a = 1

and

b = [0.2 0.2 0.2 0.2 0.2]

providing a filter of fourth order.

As an example, the functionality of the filter is illustrated in Figure 6.

The next stage of the analysis of the ECG signal is peak detection. To do this, the findpeaks function of MATLAB is employed, adjusting some parameters of this function according to the type of signal under consideration. This function returns a vector with the local maxima (peaks) of the input signal vector (data), together with the position index in which the peak is located.

The resulting vector allows the computation of RR intervals in ECG signals (see Figure 6, red signal); which, as previously mentioned, are essential values to provide as input to our classification algorithm.

3.3. Module 3: Calibration Procedure

In order to exemplify the calibration procedure, we employ a low-cost, portable, commercial device called Kardia Mobile (AliveCor Inc., San Francisco, CA, USA). It enables one-lead ECG recording. Note that any other low-cost alternative device could be used, provided that it offers an output ECG in PDF or JPEG format, as is commonly the case. Kardia by AliveCor is one of a family of mobile, clinical-quality electrocardiogram (ECG) recorders.

Furthermore, we employ professional 12-lead gold standard ECG ground truth equipment, generally used in many hospitals to monitor patients. In particular, the MAC800 by General Electric Healthcare is the one selected. The MAC 800 system provides the ability to store, analyze and display ECGs in a huge variety of formats with utmost accuracy. Both devices are shown in Figure 7.

Kardia Mobile and MAC800 are capable of monitoring patients and representing their electrocardiograms. However, the 12-lead MAC800 is a much more accurate device than the single lead Kardia Mobile. Therefore, to be comparable, data gathered by Kardia Mobile have to be calibrated, taking as a reference the data acquired by the MAC800 gold standard in order to attain the same quality results.

To achieve comparable results, we have to adjust the data acquired by the low-cost device with respect to the data measured by the more powerful ECG equipment (gold standard). The calibration procedure is carried out by comparing the morphological signals in both devices, which is the signal contained in the Lead I of each ECG. The lead placement configuration in both devices is shown in Figure 8. Due to the fact that Kardia Mobile only has a Lead I configuration, which goes from the RA (right arm) to the LA (left arm), the lead placement of the MAC800 needs to have a Lead I as well. In our case, the Mason-Likar placement [39] is employed in MAC800 because this configuration contains the Lead I signal (other leads are not used), necessary to be compared with the Lead I signal of the Kardia Mobile. The wires and adhesive electrodes of the MAC800 are much less vulnerable to artifacts and variability, providing much cleaner and more stable ECG recording with respect to the Kardia Mobile ECG (see Figure 8). In fact, we have observed in our study that the PQ and QT intervals of the Kardia Mobile had different lengths compared to the same intervals of the MAC800 gold standard (see Figure 9). This means that the morphology of the Kardia Mobile ECG waveform slightly differs from the MAC800 signal. Consequently, it could be expected that the deviation in the waveform will affect the precision of Kardia Mobile’s RR interval data measurements, and therefore, the atrial fibrillation detection procedure.

For a deeper insight into the accuracy of the Kardia Mobile in comparison to the MAC800 equipment, a set of ECGs is taken for the same patient using both devices (Kardia Mobile and MAC800) simultaneously.

Although the absolute error of the Kardia Mobile acquired data is relatively low compared to the MAC800 recordings, it is not negligible, and the corresponding deviation should be taken into consideration. Figure 9 shows four examples of each Lead-I ECG representing the difference between the Kardia Mobile signal in blue and the MAC800 Healthcare signal in red. It should be noted that the scale in the y-axis is notably different between both devices, due to the higher gain of the MAC800 with respect to the Kardia Mobile. Note, however, that this particularly significant gap does not affect the calibration procedure since the RR interval is the metric of interest.

In order to provide a better estimation in the calibration procedure for the Kardia Mobile measurements of RR intervals -namely

R R_{K M}

- we propose a common modification of

R R_{K M}

that reduces the deviations of the measurements given by MAC800 (

R R_{M}

) by means of a linear regression:

R R_{M} = a + b R R_{K M} .

(7)

The estimated values for the linear regression parameters in Equation (7) are

a = - 0.011794

and

b = 1.0679

, with a coefficient of determination of

R^{2} = 0.91

, providing an excellent goodness of fit measure. Therefore, according to these estimated values, to calibrate the RR interval samples given by Kardia Mobile, we have to increase each measurement

R R_{K M}

by 6.79% and then use a global correction of all the samples of

- 0.011794

Figure 10 illustrates how the proposed adjustment reduces the deviation of the Kardia Mobile RR intervals with respect to the RR intervals of the MAC800. Although the model has been trained with data sets of 20 patients (400 RR interval data samples), the calibration is intended to be applied to each patient individually. Thus, if we measure the RR intervals of an unseen patient with the Kardia Mobile, it will be transformed into other more precise RR intervals after applying the calibration procedure.

It is important to emphasize that we have described the calibration procedure (and its validation below) for illustrative purposes for future practitioners. The procedure is intended to improve the accuracy of the RR intervals extracted from any low-cost biomedical sensor. Nevertheless, our specific mathematical framework (ReAD-AF) is slightly affected by linear transformations and the difference before and after the calibration can be calculated (see Appendix C). However, this difference only affects the predictive power (accuracy) of the algorithm in the fourth decimal at most, as can be seen in Table 1.

4. Results

4.1. Validation of the Calibration Procedure

In order to validate the calibration procedure, we have measured the RR intervals of 20 different patients (60% men, mean age

25.3 \pm 6.2

years) with both the Kardia Mobile and the MAC800 devices, gathering 20 samples from each patient, with a total sample size of 400. Each patient was asked to provide formal consent for the ECG data acquisition. All subjects agreed to be involved in the study. The clinical state of the subjects was stable at the time of ECG acquisition. The center in which the tests were taken is the Technical University of Cartagena, whose informed consent form was approved by the local bioethics committee. After analyzing the ECGs, all subjects were with normal sinus rhythms at the time of Kardia Mobile and MAC800 ECG data acquisition. There were no patients with sustained arrhythmia. Thus, we trained a linear regression model with the RR intervals acquired from the MAC800 gold standard to correct the RR intervals provided by Kardia Mobile.

As an estimation of the error after calibration, the mean squared error (MSE) metric is employed. Specifically, given a sample

{(X_{i}, Y_{i})}_{i = 1}^{n}

of size n of the bidimensional variable

(X, Y)

, the MSE measures the mean squared differences between the two variables; that is:

M S E (X, Y) = \frac{1}{n} \sum_{i = 1}^{n} {(X_{i} - Y_{i})}^{2}

(8)

The MSE metric ranges from zero to infinite, especially penalizing those data points that are far from each other. In our case, the RR intervals of both the Kardia Mobile and MAC800 devices are compared by using the MSE. First, the mean squared error of the RR interval samples from 20 patients between the Kardia Mobile (without calibration,

R R_{K M}

) and the MAC800 (

R R_{M}

) device is calculated, achieving the following result (9):

M S E (R R_{K M}, R R_{M}) = \frac{1}{n} \sum_{i = 1}^{n} {(R R_{K M i} - R R_{M i})}^{2} = 0.001971 .

(9)

Second, the mean squared error of the RR interval samples from the 20 patients between the calibrated Kardia Mobile (

R R_{K M}^{*}

) and the MAC800 device is also computed, obtaining the following result (10):

M S E (R R_{K M}^{*}, R R_{M}) = \frac{1}{n} \sum_{i = 1}^{n} {(R R_{K M i}^{*} - R R_{M i})}^{2} = 0.000355 .

(10)

Observing these results, it is concluded that the MSE of the calibrated device is much lower. This means that when the calibration process is applied, the resulting values are closer to the values of the RR interval samples taken by the MAC800 gold standard professional equipment.

4.2. Validation of the Classification Model

For the atrial fibrillation detection procedure, our ReAD-AF algorithm was trained with publicly available ECGs taken from Physionet Challenge 2017 [40]. As a description of the challenge database, the training set contains single lead ECG recordings lasting from 9 s to just over 60 s, and the test set contains 3658 ECG recordings of similar duration. There are four different classes of ECG waveforms—normal rhythm, AF rhythm, other rhythm, and noisy recordings. We have selected the ECGs classified as AF and NS since the other signals are not of interest in this study. It should be noted that these reference ECG data were actually transferred from the Kardia Mobile device company to the Physionet Challenge 2017, which, without any loss of generality, was one of the reasons for selecting Kardia Mobile as the low-cost device to use in this study. Furthermore, the whole set of ECGs were calibrated through the calibration procedure mentioned in Section 3.3. Thus, we tested our algorithm with both sets of ECGs (raw and calibrated). To provide validation of the classification model, three metrics are employed. These metrics are sensitivity, specificity, and accuracy, as defined in Section 2. Sensitivity is a measure of the proportion of patients with AF correctly diagnosed or identified. Specificity refers to the proportion of healthy patients that are appropriately identified. Finally, accuracy measures the share of data taken from NS or AF patients that is correctly identified. The final performance attained by our algorithm is presented in Table 1.

As discussed in the previous section, and due to the fact that SRQA is invariant under monotonic transformations, Table 1 shows that the sensitivity, specificity, and accuracy attained by our logistic model in both, raw and calibrated data sets are quite similar, which means that our approach has good power performance even with non-calibrated data. This invariance property of the SRQA approach enables the obtaining of robust results with commercial low-cost sensors, because their inherent low-accuracy does not affect the final NS or AF classification.

In order to compare the results attained by our algorithm with the existing methods in literature, Table 2 shows a comparison of the sensitivity, specificity, and accuracy of some relevant algorithms.

Table 2 illustrates that our algorithm employs a relatively small number of features, which translates into a more efficient computational cost of the method. Note, however, that even using a smaller number of features, it attains good performance power in terms of sensitivity, sensibility, and accuracy; even better (in some cases) than other algorithms. More concretely, the ReAD-AF algorithm detects AF in 87.65% of the cases, discerns healthy patients in 91.86% of the cases, and its predictive power (accuracy) is 91.33%.

5. Conclusions

In this paper, we have described the methodology for designing a prototype, but fully functional application, which meets the following requirements: (i) it is robust against low-quality data, therefore enabling the use of multiple low-cost wrist-worn sensors generally characterized by low accuracy in comparison to professional electrocardiographs; and (ii) it is effective for analyzing dynamic changes based on a time series of RR intervals.

Symbolic recurrence quantification analysis (SRQA) is applied to implement a predictive scheme for the detection of AF by using, in principle, any commercially available sensor. Our approach only employs the RR intervals taken from an ECG signal in PDF or JPEG format, which further reduces the complexity of the pre-processing phase, enabling an efficient algorithm. Despite its good power performance for AF detection, it does not take into account the medical history of the patient and is restricted to the detection of AF. As a future research line we plan to extend the logistic model to a multinomial scenario where several types of arrhythmia (the most prevalent ones) will be classified and explained with new covariates extracted from the medical history of the patient.

In addition, it should be mentioned that our first preliminary tests show the computational efficiency of the entire application here introduced, which is a fundamental step towards its portability to mobile devices and, thus, its wide potential to be provisioned in those communities without sufficient resources to access medical care.

Finally, our study is intended to offer practitioners a foundation for further development of apps in the emerging digital health field based on low-cost sensors.

Author Contributions

J.P.-V. and A.-J.G.-S. conceived the idea; A.-J.G.-S. and J.G.-H. proposed the scientific methodology; J.P.-V. conceived the application and the calibration procedure; M.R.M. contributed to the derivation of the analytical procedure; A.-J.G.-S. and J.G.-H. improved technical issues. All authors participated in the elaboration of the manuscript and its revision. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This research has been supported by the projects AIM, ref. TEC2016-76465-C2-1-R (AEI/FEDER, UE), e-DIVITA, ref.20509/PDC/18 (Proof of Concept, 2018) and ATENTO, ref. 20889/PI/18 (Fundación Seneca, Región de Murcia). This study is part of the collaborative activities carried out under the program Groups of Excellence of the region of Murcia, the Fundación Séneca, Science and Technology Agency of the region of Murcia (Spain) project 19884/GERM/15.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Functional diagram of the application implemented in MATLAB (a). Remote access to the developed application in the cloud (b).

Appendix B

Algorithm A1 ReAD-AF

1: Input: RR interval data set of a patient
2: $r a w_d a t a \leftarrow$ load(RR)
3: $d a t a \leftarrow$ splitRR(raw_data, w) ▹ Split RR intervals in windows of size w
4: $m e a n R R \leftarrow$ mean(data) ▹ Mean
5: $m e d i a n R R \leftarrow$ median(data) ▹ Median
6: $V R R \leftarrow$ std(data)/mean(data) ▹ Pearson coefficient of variation
7: $V m e R R \leftarrow$ sum(abs(data-medianRR))/medianRR; ▹ Mean absolute dispersion
8: SRR: symbolic recurrence rate
9: SRP: symbolic recurrence plots
10: Vinc(1,3): Shannon entropy prob. distrib. vertical line lengths; SRP for increasing symbol {i}
11: Vinc(1,4): Mean length of vertical lines; SRP for increasing symbol {i}
12: Vdec(1,3): Shannon entropy prob. distrib. vertical line lengths; SRP for decreasing symbol {i}
13: Vdec(1,4): Mean length of vertical lines; SRP for decreasing symbol {i}
14: DETt: The percentage of recurrence points which form diagonal lines; SRP for all symbols
15: ENTRt: Shannon entropy of the probability distribution of the diagonal; SRP for all symbols
16: $M a t r i x \leftarrow$ func_SRP(data, m) ▹ Matrix SRR for all symbols
17: $[V i n c (1, 3), V i n c (1, 4)] \leftarrow$ recu_SRQA(Matrix {i})
18: $[V d e c (1, 3), V d e c (1, 4)] \leftarrow$ recu_SRQA(Matrix {i})
19: $[D E T t, E N T R t] \leftarrow$ recu_SRQA(Matrix)
20: measure: Storage of all metrics in a matrix
21: $m e a s u r e \leftarrow$ [meanRR,medianRR,VRR,VmeRR,DETt, ENTRt,Vinc(1,3),Vdec(1,3),Vinc(1,4),Vdec(1,4)]
22: pos_af: The boolean vector identifying the AF windows
23: $m o d \leftarrow$ fitglm(measure,’Distribution’,’binomial’,pos_af) ▹ Compute the logistic model
24: $p r o b \leftarrow$ mod.Fitted.Probability ▹ Compute the estimated probability of AF
25: $T l i s t \leftarrow$ [0.001:0.001:1]
26: for j=1 to j=lenght(Tlist) do
27: $t h s \leftarrow$ Tlist(j) ▹ Compute the threshold
28: $T P \leftarrow$ length(find(prob(find(pos_af>0))>ths)) ▹ True Positive
29: $T N \leftarrow$ length(find(prob(find(pos_af==0))<ths)) ▹ True Negative
30: $F N \leftarrow$ length(find(prob(find(pos_af>0))<ths)) ▹ False Negative
31: $F P \leftarrow$ length(find(prob(find(pos_af==0))>ths)) ▹ False Positive
32: $T P R \leftarrow$ TP/(TP+FN) ▹ Sensitivity or True Positive Rate
33: $F P R \leftarrow$ FP/(FP+TN) ▹ False Positive Rate
34: $S P C \leftarrow$ 1-FPR ▹ Specificity
35: $A C C \leftarrow$ (TP+TN)/(TP + FN + FP + TN) ▹ Accuracy
36: end for
37: $τ$ : Optimal threshold parameter
38: $τ \leftarrow m i n {F P R_{τ}^{2} + {(1 - S e_{τ})}^{2}}$ ▹ Compute the optimal threshold for AF classification
39: if $p r o b \geq τ$ do
40: patient classified as AF
41: else
42: patient classified as NS
43: end if
44: Output: Classified patient

Appendix C

Let

{x_{t}}_{t \in I}

be a real valued time series and

f : R \to R

a monotonic real map transforming

x_{t} \to f (x_{t}) = y_{t}

. Since f is a monotonic map, it preserves the order relations

″ < ″

and

″ > ″

, that is if

x_{t} < x_{t + k}

(respectively

x_{t} > x_{t + k}

) then

f (x_{t}) = y_{t} < f (x_{t + k}) = y_{t + k}

(respectively

f (x_{t}) = y_{t} > f (x_{t + k}) = y_{t + k}

). Therefore,

S ({\bar{x}}_{t}) = S ({\bar{y}}_{t})

for all

t \in I

and thus, the symbolic recurrence measures associated to the time series

{x_{t}}_{t \in I}

and

{f (x_{t}) = y_{t}}_{t \in I}

are the same.

On the other hand, the covariates used in the logistic model for AF classification based on

R R

interval measures given in [31] are divided into two groups. The first group is formed of covariates coming from symbolic recurrence measures. The second group is formed of the arithmetic mean (

\bar{R R}

), median (

M e_{R R}

), Pearson coefficient of variation (

V_{R R} = \frac{s_{R R}}{\bar{R R}}

where

s_{R R}

is the standard deviation of

R R

) and coeficient of variation of the median (

V_{M e} = \frac{D_{M e_{R R}}}{M e}

where

D_{M e_{R R}}

is the mean absolute deviation of

R R

with respect to

M e_{R R}

Let the log-odds of the logistic model be

l o g (\frac{p}{1 - p}) = α_{0} + α_{1} \bar{R R} + α_{2} M e_{R R} + α_{3} V + α_{4} V_{M e} + \sum_{i} β_{i} s r m_{i}

(A1)

where

s r m_{i}

are symbolic recurrence measures. Denote by

p_{K M}

and

p_{M}

the probabilities of AF for the Kardia Mobile measurements of RR interval and its adjustment given by a linear transformation, namely

R R_{K M}

and

R R_{M} = a + b R R_{K M}

, respectively. It is straightforward to check that

{\bar{R R}}_{M} = a + b {\bar{R R}}_{K M}

M e_{R R_{M}} = a + b M e_{R R_{K M}}

s_{R R_{M}} = | b | s_{R R_{K M}}

, and

D_{M e_{R R_{M}}} = b D_{M e_{R R_{K M}}}

. Then, since the symbolic recurrence measures are invariant under monotonic transformations, it follows that the difference of the log-odds for

R R_{M}

and

R R_{K M}

is given by:

\begin{matrix} l o g (\frac{p_{M}}{1 - p_{M}}) - l o g (\frac{p_{K M}}{1 - p_{K M}}) = & a (α_{1} + α_{2}) + (b - 1) α_{1} {\bar{R R}}_{K M} + (b - 1) α_{2} M e_{K M} \\ - \frac{α_{3} a}{a + b {\bar{R R}}_{K M}} V_{R R_{K M}} - \frac{α_{4} a}{a + b M e_{K M}} V_{M e_{R R_{K M}}} . \end{matrix}

(A2)

Therefore, when a is close to 0 and b is close to 1, as is the case in our adjustment (

a = - 0.011794

and

b = 1.0679

), the difference of the log-odds computed in (A2) is very small, which means that the predicted power of the logistic model is not affected by the proposed linear adjustment.

References

Centers for Disease Control and Prevention. The State of Aging and Health in America 2013. Atlanta, GA: Centers for Disease Control and Prevention; US Dept of Health and Human Services: Washington, DC, USA, 2013. [Google Scholar]
Stadhouders, N.; Kruse, F.; Tanke, M.; Koolman, X.; Jeurissen, P. Effective healthcare cost-containment policies: A systematic review. Health Policy 2019, 123, 71–79. [Google Scholar] [CrossRef] [PubMed]
Biglu, B.H.; Mohammad-Hossein, G.M.; Biglu, S. Cardiovascular diseases in the mirror of science. J. Cardiovasc. Thorac. Res. 2016, 8, 158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cardiovascular Diseases (CVDs). World Health Organization. Available online: https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases- (accessed on 21 June 2019).
Deen, M.J. Information and Communications Technologies for Elderly Ubiquitous Healthcare in a Smart Home. Pers. Ubiquitous Comput. 2015, 19, 573–599. [Google Scholar] [CrossRef]
Wang, H.; Choi, H.; Agoulmine, N.; Jamal Deen, M.; Won-Ki Hong, J. Information-based sensor tasking wireless body area networks in U-health systems. In Proceedings of the 2010 International Conference on Network and Service Management, Niagara Falls, ON, Canada, 25–29 October 2010; pp. 517–522. [Google Scholar]
Huynh, N. How the “Big 4” Tech Companies Are Leading Healthcare Innovation. Healthcare Weekly. Available online: https://healthcareweekly.com/how-the-big-4-tech-companies-are-leading-healthcare-innovation/ (accessed on 27 February 2019).
Sodhro, A.; Sangaiah, A.; Sodhro, G.; Lohano, S.; Pirbhulal, S. An Energy-Efficient Algorithm for Wearable Electrocardiogram Signal Processing in Ubiquitous Healthcare Applications. Sensors 2018, 18, 923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Muzammal, M.; Talat, R.; Sodhro, A.H.; Pirbhulal, S. A multi-sensor data fusion enabled ensemble approach for medical data from body sensor networks. Inf. Fusion 2020, 53, 155–164. [Google Scholar] [CrossRef]
Censi, F.; Corazza, I.; Reggiani, E.; Calcagnini, G.; Mattei, E.; Triventi, M.; Boriani, G. P-wave Variability and Atrial Fibrillation. Sci. Rep. 2016, 6, 1–7. [Google Scholar] [CrossRef] [PubMed]
Aytemir, K.; Özer, N.; Atalar, E.; Sade, E.; Aksöyek, S.; Övünç, K.; Oto, A.; Özmen, F.; Kes, S.P. Wave Dispersion on 12-Lead Electrocardiography in Patients with Paroxysmal Atrial Fibrillation. Pacing Clin. Electrophysiol. 2000, 23, 1109–1112. [Google Scholar] [CrossRef]
Andreotti, F.; Carr, O.; Pimentel, M.A.F.; Mahdi, A.; De Vos, M. Comparing feature-based classifiers and convolutional neural networks to detect arrhythmia from short segments of ECG. In Proceedings of the 2017 Computing in Cardiology (CinC), Rennes, France, 24–27 September 2017; pp. 1–4. [Google Scholar] [CrossRef]
Xia, Y.; Wulan, N.; Wang, K.; Zhang, H. Detecting atrial fibrillation by deep convolutional neural networks. Comput. Biol. Med. 2018, 93, 84–92. [Google Scholar] [CrossRef]
Bollepalli, S.C.; Challa, S.S.; Jana, S.; Patidar, S. Atrial Fibrillation Detection Using Convolutional Neural Networks. In Proceedings of the 2017 Computing in Cardiology (CinC), Rennes, France, 24–27 September 2017. [Google Scholar] [CrossRef]
Shashikumar, S.P.; Shah, A.J.; Li, Q.; Clifford, G.D.; Nemati, S. A deep learning approach to monitoring and detecting atrial fibrillation using wearable technology. In Proceedings of the 2017 IEEE EMBS International Conference on Biomedical & Health Informatics (BHI), Orlando, FL, USA, 16–19 February 2017. [Google Scholar] [CrossRef]
Billeci, L.; Costi, M.; Lombardi, D.; Chiarugi, F.; Varanini, M. Automatic Detection of Atrial Fibrillation and Other Arrhythmias in ECG Recordings Acquired by a Smartphone Device. Electronics 2018, 7, 199. [Google Scholar] [CrossRef] [Green Version]
Steinhubl, S.R.; Waalen, J.; Edwards, A.M.; Ariniello, L.M.; Mehta, R.R.; Ebner, G.S.; Carter, C.; Baca-Motes, K.; Felicione, E.; Sarich, T.; et al. Effect of a Home-Based Wearable Continuous ECG Monitoring Patch on Detection of Undiagnosed Atrial Fibrillation. JAMA 2018, 320, 146. [Google Scholar] [CrossRef] [Green Version]
Perez, M.V.; Mahaffey, K.W.; Hedlin, H.; Rumsfeld, J.S.; Garcia, A.; Ferris, T.; Balasubramanian, V.; Russo, A.M.; Rajmane, A.; Cheung, L.; et al. Large-Scale Assessment of a Smartwatch to Identify Atrial Fibrillation. N. Engl. J. Med. 2019, 381, 1909–1917. [Google Scholar] [CrossRef] [PubMed]
Apple. Using Apple Watch for Arrhythmia Detection. 2013. Available online: https://www.apple.com/healthcare/docs/site/Apple_Watch_Arrhythmia_Detection.pdf (accessed on 6 February 2020).
Lian, J.; Wang, L.; Muessig, D. A Simple Method to Detect Atrial Fibrillation Using RR Intervals. Am. J. Cardiol. 2011, 107, 1494–1497. [Google Scholar] [CrossRef] [PubMed]
Tateno, K.; Glass, L. A Method for Detection of Atrial Fibrillation Using RR Intervals. Comput. Cardiol. 2000, 27, 391–394. [Google Scholar] [CrossRef]
Dash, S.; Chon, K.H.; Lu, S.; Raeder, E.A. Automatic Real Time Detection of Atrial Fibrillation. Ann. Biomed. Eng. 2009, 37, 1701–1709. [Google Scholar] [CrossRef] [PubMed]
Apple Watch Series 4. Available online: https://www.apple.com/apple-watch-series-4/ (accessed on 21 June 2019).
MySignals—eHealth and Medical IoT Development Platform. Available online: https://www.the-iot-marketplace.com/solutions/ehealth (accessed on 21 June 2019).
Qardio—State of the Heart Technology. Available online: https://www.getqardio.com/ (accessed on 21 June 2019).
AliveCor. Available online: https://www.alivecor.com/ (accessed on 21 June 2019).
Koltowski, L.; Balsam, P.; Glłowczynska, R.; Rokicki, J.K.; Peller, M.; Blicharz, L.; Maciejewski, K.; Niedziela, M.; Opolski, G.; Grabowski, M. Kardia Mobile applicability in clinical practice: A comparison of Kardia Mobile and standard 12-lead electrocardiogram records in 100 consecutive patients of a tertiary cardiovascular care center. Cardiol. J. 2014. [Google Scholar] [CrossRef] [Green Version]
Evans, G.F.; Shirk, A.; Muturi, P.; Soliman, E.Z. Feasibility of Using Mobile ECG Recording Technology to Detect Atrial Fibrillation in Low-Resource Settings. Glob. Heart 2017, 12, 285–289. [Google Scholar] [CrossRef]
GE Healthcare. Available online: https://www.gehealthcare.com/ (accessed on 21 June 2019).
MyHeartFitnessApp. Available online: https://github.com/JesusUPCT/MyHeartFitnessApp (accessed on 23 January 2020).
Pérez-Valero, J.; Caballero Pintado, M.V.; Melgarejo, F.; García-Sánchez, A.J.; Garcia-Haro, J.; García Córdoba, F.; García Córdoba, J.A.; Pinar, E.; García Alberola, A.; Matilla-García, M.; et al. Symbolic Recurrence Analysis of RR Interval to Detect Atrial Fibrillation. J. Clin. Med. 2019, 8, 1840. [Google Scholar] [CrossRef] [Green Version]
Caballero, M.; Garcia, M.; Ruiz Marín, M. Symbolic recurrence plots to analyze dynamical systems. Chaos Interdisc. J. Nonlinear Sci. 2018, 28, 063112. [Google Scholar] [CrossRef]
Takens, F. Detecting strange attractors in turbulence. In Dynamical Systems and Turbulence, Warwick 1980; Springer: Berlin/Heidelberg, Germany, 1981; pp. 366–381. [Google Scholar]
Moody, G.B.; Mark, R.G. MIT-BIH Arrhythmia Database; Harvard-MIT Division of Health Sciences and Technology, Biomedical Engineering Center: Cambridge, MA, USA, 1992. [Google Scholar] [CrossRef]
Darweesh, M.; Habte, T.; Saleh, H.; Mohammad, B.; Ismail, M. Novel Logarithmic ECG Feature Extraction Algorithm Based on Pan and Tompkins. In Proceedings of the 2016 IEEE 59th International Midwest Symposium on Circuits and Systems (MWSCAS), Abu Dhabi, UAE, 16–19 October 2016. [Google Scholar] [CrossRef]
Silva, I.; Moody, G.B. An Open-source Toolbox for Analysing and Processing PhysioNet Databases in MATLAB and Octave. J. Open Res. Softw. 2014, 2, e27. [Google Scholar] [CrossRef] [PubMed]
Goldberger, A.L.; Amaral, L.A.; Glass, L.; Hausdorff, J.M.; Ivanov, P.C.; Mark, R.G.; Mietus, J.E.; Moody, G.B.; Peng, C.K.; Stanley, H.E. PhysioBank, PhysioToolkit, and PhysioNet: Components of a new research resource for complex physiologic signals. Circulation 2000, 101, e215–e220. [Google Scholar] [CrossRef] [Green Version]
Shrivastava, P.R.K.; Panbude, S.; Narayanan, G.S. Digitization of ECG Paper Records using MATLAB. Int. J. Innov. Explor. Eng. IJITEE 2014, 4. Available online: https://pdfs.semanticscholar.org/61a7/3784741a2f6dc5b2e459f235f695c33aa7e1.pdf (accessed on 11 February 2020).
Papouchado, M.; Walker, P.R.; James, M.A.; Clarke, L.M. Fundamental differences between the standard 12-lead electrocardiograph and the modified (Mason–Likar) exercise lead system. Eur. Heart J. 1987, 8, 725–733. [Google Scholar] [CrossRef] [PubMed]
Clifford, G.D.; Liu, C.; Moody, B.; Li-wei, H.L.; Silva, I.; Li, Q.; Johnson, A.E.; Mark, R.G. AF Classification from a Short Single Lead ECG Recording: The PhysioNet/Computing in Cardiology Challenge 2017. In Proceedings of the 2017 Computing in Cardiology (CinC), Rennes, France, 24–27 September 2017. [Google Scholar]
Behar, J.A.; Rosenberg, A.A.; Yaniv, Y.; Oster, J. Rhythm and quality classification from short ECGs recorded using a mobile device. In Proceedings of the 2017 Computing in Cardiology (CinC), Rennes, France, 24–27 September 2017; pp. 1–4. [Google Scholar] [CrossRef]
Christov, I.; Krasteva, V.; Simova, I.; Neycheva, T.; Schmid, R. Multi-parametric analysis for atrial fibrillation classification in ECG. In Proceedings of the 2017 Computing in Cardiology (CinC), Rennes, France, 24–27 September 2017; pp. 1–4. [Google Scholar] [CrossRef]

Figure 1. SRPs of a 102-length RR interval time series for an embedding dimension m = 3 in a normal sinus (a) and atrial fibrillation (b) patient.

Figure 2. Block diagram showing the main modules that shape the prototype application developed. After the entire procedure, the objective is to determine if a patient is with normal sinus rhythm (NS) or with atrial fibrillation (AF).

Figure 3. Block diagram showing the steps of the digitalization procedure.

Figure 4. Horizontal scanning (green line) vs vertical scanning (red line). Arrows represent the direction of the scanning process.

Figure 5. Transformation process of the signal from PDF or JPEG formats (a), then a binary mask image is obtained (b) and after that, some pixels are removed (c). Finally, a shrink operation is applied to capture only useful information (d), and then the output digitalized signal (e) is shown.

Figure 6. Transformation of the signal before (a) and after (b) applying the filter function

Figure 7. Kardia Mobile device by AliveCorc (a) and MAC800 by General Electric Healthcare (b); both employed for the calibration procedure.

Figure 8. Electrocardiograph (ECG) placement configuration in both devices MAC800 (a) and Kardia Mobile (b)

Figure 9. Different ECG deviations from (a–d) of the Kardia Mobile (blue) with respect to the MAC800 (red) equipment.

Figure 10. RR interval data for 20 samples from 20 different patients measured with Kardia Mobile (blue), MAC800 (red) and corrected Kardia Mobile (green).

Table 1. Sensitivity, specificity, and accuracy of our classification model applied to the raw data set (left) and calibrated data set (right).

Sensitivity	Specificity	Accuracy	Sensitivity	Specificity	Accuracy
0.8765	0.9186	0.9133	0.8765	0.9184	0.9131

Table 2. Comparison of the sensitivity, specificity, and accuracy between some Physionet Challenge 2017 methods and our algorithm.

Reference	Number of Features	Classification Method	Se	Sp	Acc
Billeci et al. [16]	30	Least Square-Support Vector Machine	0.83	0.96	0.95
Behar et al. [41]	42	Feature-Based Machine Learning	0.81	0.91	0.90
Christov et al. [42]	>40	Multi-Parametric Analysis	0.82	0.90	0.89
ReAD-AF [31]	9	Symbolic Recurrence Analysis	0.87	0.91	0.91

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Pérez-Valero, J.; Garcia-Sanchez, A.-J.; Ruiz Marín, M.; Garcia-Haro, J. A Prototype Framework Design for Assisting the Detection of Atrial Fibrillation Using a Generic Low-Cost Biomedical Sensor. Sensors 2020, 20, 896. https://doi.org/10.3390/s20030896

AMA Style

Pérez-Valero J, Garcia-Sanchez A-J, Ruiz Marín M, Garcia-Haro J. A Prototype Framework Design for Assisting the Detection of Atrial Fibrillation Using a Generic Low-Cost Biomedical Sensor. Sensors. 2020; 20(3):896. https://doi.org/10.3390/s20030896

Chicago/Turabian Style

Pérez-Valero, Jesús, Antonio-Javier Garcia-Sanchez, Manuel Ruiz Marín, and Joan Garcia-Haro. 2020. "A Prototype Framework Design for Assisting the Detection of Atrial Fibrillation Using a Generic Low-Cost Biomedical Sensor" Sensors 20, no. 3: 896. https://doi.org/10.3390/s20030896

APA Style

Pérez-Valero, J., Garcia-Sanchez, A.-J., Ruiz Marín, M., & Garcia-Haro, J. (2020). A Prototype Framework Design for Assisting the Detection of Atrial Fibrillation Using a Generic Low-Cost Biomedical Sensor. Sensors, 20(3), 896. https://doi.org/10.3390/s20030896

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

A Prototype Framework Design for Assisting the Detection of Atrial Fibrillation Using a Generic Low-Cost Biomedical Sensor

Abstract

1. Introduction

2. Normal Sinus and Atrial Fibrillation Classification Scheme

3. Prototype Application

3.1. Module 1: ECG Digitization

3.2. Module 2: Signal Processing and Filtering

3.3. Module 3: Calibration Procedure

4. Results

4.1. Validation of the Calibration Procedure

4.2. Validation of the Classification Model

5. Conclusions

Author Contributions

Acknowledgments

Conflicts of Interest

Appendix A

Appendix B

Appendix C

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI