The present application generally relates to biomedical signal processing. Particularly, the application provides a method and system for removing corruption in photoplethysmogram (PPG) signals captured using a mobile communication device. More particularly, the application provides a method and system for removing corruption in photoplethysmogram (PPG) signals for monitoring cardiac health of patients.
An intriguing problem in healthcare analytics, particularly in photoplethysmogram (PPG) signals based proactive cardiac health monitoring, is the frequent presence of corruption in PPG signal from multiple noise sources like motion artifacts, ambient noise. Data preprocessing through denoising results in improved interpretation of the cardiovascular systems and removal of false alarms. In fact, the occurrence of high false alarms (precisely false positives due to the presence of noise) leads to ‘alarm fatigue’. In contrary, high false negatives would be fatal to the patients.
In all practicality, it is infeasible to capture ECG or ABP signals without using extra sensors or through invasive procedures. It is observed that majority of computational techniques do not find widespread clinical use because of their lack of capability in reducing the artifacts and other noises. Prior art illustrates 80 different artifact detection techniques and it is perceived that most of the techniques are highly explicit, hard coded to suit for certain device settings, usage and thus limiting their applicability in almost all practical purposes. Current trend is towards multi-signal corruption analysis and detection. The underlying assumption is that all of the cardiological signals like viz. PPG, ECG, ABP cannot be corrupted simultaneously while source of corruption is independent. Stand-alone method of PPG signal denoising and quality assessment is presented in prior art, which employs computationally heavy machine learning algorithms and multivariate ‘voting’ threshold mechanism. Such schemes are not suitable for smartphones to deliver real-time performance. Another problem is the single-stage error detection is itself erroneous and prone to have high false negatives.
Thereby, removing corruption in photoplethysmogram (PPG) signals for monitoring cardiac health of patient's negatives is still considered to be one of the biggest challenges of the technical domain.
Embodiments of the present disclosure present technological improvements as solutions to one or more of the above-mentioned technical problems recognized by the inventors in conventional systems. For example, in one embodiment, a method for removing corruption in photoplethysmogram (PPG) signals for monitoring cardiac health of patients is provided. The method comprises extracting photoplethysmogram signals from the patient, detecting and eliminating corruption caused by larger and transient disturbances in the extracted photoplethysmogram (PPG) signals, segmenting photoplethysmogram (PPG) signals post detection and elimination of corruption caused by larger and transient disturbances, identifying of inconsistent segments from the segmented photoplethysmogram (PPG) signals, detecting anomalies from the identified inconsistent segments of the photoplethysmogram (PPG) signals, analysing the detected anomalies of the photoplethysmogram (PPG) signals and identifying photoplethysmogram (PPG) signal segments corrupted by smaller and prolonged disturbances.
In another embodiment, a system for removing corruption in photoplethysmogram (PPG) signals for monitoring cardiac health of patients is provided. The system (200) comprises of an image capturing device (202) coupled with a mobile communication device (204), an extrema elimination module (206), a segmentation module (208), an inconsistency identification module (210), an anomaly detection module (212), an anomaly analytics module (214) and a decision module (216)..
In yet another embodiment, a non-transitory computer readable medium storing instructions is disclosed wherein when executed by a possessor on a system, the instructions cause the processor to perform method for monitoring cardiac health of patients comprising extracting photoplethysmogram signals from the patient, detecting and eliminating corruption caused by larger and transient disturbances in the extracted photoplethysmogram (PPG) signals, segmenting photoplethysmogram (PPG) signals post detection and elimination of corruption caused by larger and transient disturbances, identifying of inconsistent segments from the segmented photoplethysmogram (PPG) signals, detecting anomalies from the identified inconsistent segments of the photoplethysmogram (PPG) signals, analysing the detected anomalies of the photoplethysmogram (PPG) signals and identifying photoplethysmogram (PPG) signal segments corrupted by smaller and prolonged disturbances.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention, as claimed.
The accompanying drawings, which are incorporated in and constitute a part of this disclosure, illustrate exemplary embodiments and, together with the description, serve to explain the disclosed principles.
Exemplary embodiments are described with reference to the accompanying drawings. In the figures, the left-most digit(s) of a reference number identifies the figure in which the reference number first appears. Wherever convenient, the same reference numbers are used throughout the drawings to refer to the same or like parts. While examples and features of disclosed principles are described herein, modifications, adaptations, and other implementations are possible without departing from the spirit and scope of the disclosed embodiments. It is intended that the following detailed description be considered as exemplary only, with the true scope and spirit being indicated by the following claims.
The present disclosure provides a method and system for removing corruption in photoplethysmogram (PPG) signals for monitoring cardiac health of patients.
In an embodiment of the present invention, a system (200) for removing corruption in photoplethysmogram (PPG) signals for monitoring cardiac health of patients comprises of an image capturing device (202) coupled with a mobile communication device (204), an extrema elimination module (206), a segmentation module (208), an inconsistency identification module (210), an anomaly detection module (212), an anomaly analytics module (214) and a decision module (216).
Referring to
The process starts at step 102, photoplethysmogram signals are extracted from the user using the image capturing device (202) coupled with the mobile communication device (204). At the step 104, corruption caused by larger and transient disturbances in the extracted photoplethysmogram signals are detected and eliminated. At the step 106, photoplethysmogram (PPG) signals post detection and elimination of corruption caused by larger and transient disturbances are segmented. At step 108, inconsistent segments from the segmented photoplethysmogram signals are identified. At step 110, anomalies from the identified inconsistent segments of the photoplethysmogram signals are detected. At step 112, the detected anomalies of the photoplethysmogram signals are analyzed and at step 114, photoplethysmogram (PPG) signal segments corrupted by smaller and prolonged disturbances are identified.
Referring to
In another embodiment of the present invention, the image capturing device (202) coupled with the mobile communication device (204) is adapted for extracting photoplethysmogram signals from the patient. The photoplethysmogram signals are extracted from patient's peripheral body parts selected from a group comprising but not limited to finger, ear, and toe. In a specific embodiment, the photoplethysmogram signals are extracted from user's forehead. The mobile communication device (204) captures photoplethysmogram signal in reflective mode. The mobile communication device (204) is selected from a group comprising of smart phone, mobile phone, laptop, tablet, and personal digital assistant.
The image capturing device (202) coupled with the mobile communication device (204) is a camera and have a light emitting source for extracting photoplethysmogram signals from the patient's peripheral body parts selected from a group comprising but not limited to finger, ear, toe; forehead, thereby, obtaining a video sequence of the light, reflected from patient's peripheral body parts.
Corruption in PPG signal induces erroneous feature extraction and eventually leads to wrong diagnosis. The corruption in PPG is classified as: 1. Extremas and 2. Intricates, where extremas are the corruptions due to larger, transient disturbance and intricates are due to smaller and mostly prolonged disturbances as shown in
In another embodiment of the present invention, the extremas are first processed for—to detect the larger variation of corruption, which would help in unbiased analysis.
In another embodiment of the present invention, let P={pn}, n=1,2, . . . , N be the captured PPG signal. Firstly, PPG is filtered to remove the residual invalide sections because human Heart Rate (HR) is a bounded variable and PPG cannot have frequency components beyond that bound. PPG signal P={pn} is passed through a low pass filter (LPF) with transfer function:
In another embodiment of the present invention, a modified Thompson Tau technique is considered to statistically decide the inconsistencies on the filtered PPG signal for finding the extremas—
where
is the critical student t's value with a=0.05 with single degree of freedom.
are the extremas, Ex ∈ P.
In an exemplary embodiment of the present invention, the extremas detection of PPG signal of patient #v748s is depicted (data taken from publicly available MIT-Physionet challenge 2015) in
In another embodiment of the present invention, the below table depicts the notational meanings of the symbols used—
With successful extremas detection, decorrupting the entire PPG becomes less error prone. In the next steps, intrcicate detection takes place.
PPG signal is a time series records sampled at a predefined frequency. It consists of series of continuous segments and each segment start is indicated by cardiac cycle onset time (
In another embodiment of the present invention, most probable onset duration is computed using a method to estimate individual's cardiac characteristics. In the said method, two parameters
In another embodiment of the present invention, PPG being a deterministic signal (at least the valid PPG signal would always maintain certain pattern) and each segment length carries that signature or pattern. The most obvious PPG cardiac segment sample is considered as the template :={t1, . . . , tM} of length M corresponds to most common human heart rate of 72 beats per minute.
To counter the nonlinearity in the PPG segments, each segment is restricted to most probable segment length ιp, i.e. Ωk={ω1, ω2, . . . , ωι
In order to compute the dissimilarity of given PPG segment, distance is measured with the computed ‘Morphologically Valid PPG Template’ (MVPT). Let P80 :={p1, p2, . . . , pL} be the PPG signal of length L ∈ for segment λ and :={t1, t2, . . . , tM} be the template MVPT of length M ∈. In order to accommodate the non-stationary nature and nonlinearity in each of the segments due to dynamicity of cardiac output, heart beat rate for single person and person to person, normalized dynamic window-adaptive dynamic time warping (NDWADTW) is introduced here. Instead of computing on the entire length of each PPG segment at cardiac cycle {ok−1−ok} (which is variable in length due to change of heart rate), the temporal-length mismatch bias is eliminated by fixing the analysis segment length of a person to its most probable segment length, i.e. L=ιp.
With two stages of decorruption performed through extremas and segment length variability detection, finally normalized dynamic window-adaptive dynamic time warping (NDWADTW) distance is computed to address the amplitude-abnormality and is applied for PPG signal cleaning. Each segment is normalized as
that minimizes the effect of nonlinear changes in the amplitude of each segment due to changes in cardiac output. Thus the similarity measurement is individual-adaptive and is not biased with the non-stationary nature and nonlinearity of the segment lengths. Normalized dynamic window-adaptive dynamic time warping (NDWADTW) distance is computed between template and each of the segments
k∀K of the extracted PPG signal of lengths M, ιp respectively as:
Where
In
To classify the PPG segments as corrupt or non-corrupt, robust extreme value analysis is performed, which requires masking effect and swamping effect are to be minimized. Masking effect causes one or more anomalous or corrupt points to get undetected, while normal or non-corrupt points are depicted as outlier by swamping effect. As it is intuitive as well as evident from
For minimizing masking effect, masking breakdown point ξM is to be close to 0.5. For a large class of identifier Hampel identifier, a nonlinear local outlier detector based on Median Absolute Deviation (MAD) scale estimation provides outlier detection with ξM→0.5 as follows:
MAD(εN)=median(|ε1−median(εN)|, . . . , |εN−median(εN)|)
An observation is identified as corrupt when: . Typically Thr=3 is considered. The Hampel identifier is applied to decide whether a DTW distance is corrupt. When Hampel identifier declares a NDWADTW distance outlier, the corresponding PPG segment is declared as corrupt.
In an exemplary embodiment of the present invention, Total 750 patients from four hospitals in the USA and Europe are chosen at random and publicly available as MIT-BIH arrhythmia Physionet CinC Challenge 2015 datasets. Corrupt-segments are manually annotated through crowd-sourcing. Random 25 PPG signals out of 750 are taken and number of annotators annotates the corrupt regions and voting result is used as ground truth. Exemplary case of de-corruption of PPG signal for patient # v748s is shown in
In another embodiment of the present invention, the dynamic multi-level cluster-based anomaly analytics implemented herein is described below:
In another embodiment of the present invention, multi-level cluster-based anomaly analytics implemented herein is dynamic in nature because of the system's ability to adapt to the changing nature of the PPG signals.
In another embodiment of the present invention, as prediction problem is evaluated by confusion matrix parameters: False Positives (FP), False Negatives (FN), True Positives (TP) and True Negatives (TN) are shown in the table below. It is observed that have achieved very low false negatives (<=4%), indicating that the very less number of non-detections. With such mechanism of PPG signal de-corruption, better heart condition status and lesser false alarms are expected.
In an exemplary embodiment of the present invention, 20 subjects with 10 lab-generated and 10 from MIMIC 2 datasets from physionet have been considered—the performance in terms of precision=78%, recall=80% and specificity=96%.
The illustrated steps are set out to explain the exemplary embodiments shown, and it should be anticipated that ongoing technological development will change the manner in which particular functions are performed. These examples are presented herein for purposes of illustration, and not limitation. Further, the boundaries of the functional building blocks have been arbitrarily defined herein for the convenience of the description. Alternative boundaries can be defined so long as the specified functions and relationships thereof are appropriately performed. Alternatives (including equivalents, extensions, variations, deviations, etc., of those described herein) will be apparent to persons skilled in the relevant art(s) based on the teachings contained herein. Such alternatives fall within the scope and spirit of the disclosed embodiments. Also, the words “comprising,” “having,” “containing,” and “including,” and other similar forms are intended to be equivalent in meaning and be open ended in that an item or items following any one of these words is not meant to be an exhaustive listing of such item or items, or meant to be limited to only the listed item or items. It must also be noted that as used herein and in the appended claims, the singular forms “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise.
Furthermore, one or more computer-readable storage media may be utilized in implementing embodiments consistent with the present disclosure. A computer-readable storage medium refers to any type of physical memory on which information or data readable by a processor may be stored. Thus, a computer-readable storage medium may store instructions for execution by one or more processors, including instructions for causing the processor(s) to perform steps or stages consistent with the embodiments described herein. The term “computer-readable medium” should be understood to include tangible items and exclude carrier waves and transient signals, i.e., be non-transitory. Examples include random access memory (RAM), read-only memory (ROM), volatile memory, nonvolatile memory, hard drives, CD ROMs, DVDs, flash drives, disks, and any other known physical storage media.
It is intended that the disclosure and examples be considered as exemplary only, with a true scope and spirit of disclosed embodiments being indicated by the following claims.
Number | Date | Country | Kind |
---|---|---|---|
201621008875 | Mar 2016 | IN | national |
This U.S. patent application claims priority under 35 U.S.C. §119 to: India Application No. 201621008875, filed on 14 Mar. 2016. The entire contents of the aforementioned application are incorporated herein by reference.