The present invention is directed to systems and methods for detecting an arrhythmic cardiac event from a cardiac signal obtained from a subject being monitored for cardiac function.
Early detection of cardiac arrhythmias can be critical for patient recovery. Increasingly sophisticated systems and methods for monitoring for various cardiac events are needed to improve diagnosis and treatment. The present invention is directed toward detecting an arrhythmic cardiac event from a cardiac signal obtained from a subject being monitored for cardiac function.
The following U.S. patents, U.S. patent applications, and Publications are incorporated herein in their entirety by reference.
“Identifying A Type Of Cardiac Event From A Cardiac Signal Segment”, U.S. patent application Ser. No. 14/492,948, by Xu et al. (Attorney Docket No. 20140386 US01).
“Estimating Cardiac Pulse Recovery From Multi-Channel Source Data Via Constrained Source Separation”, U.S. patent application Ser. No. 13/247,683, by Mestha et al.
“Continuous Cardiac Pulse Rate Estimation From Multi-Channel Source Video Data”, U.S. patent application Ser. No. 13/528,307, by Kyal et al.
“Continuous Cardiac Pulse Rate Estimation From Multi-Channel Source Video Data With Mid-Point Stitching”, U.S. patent application Ser. No. 13/871,728, by Kyal et al.
“Continuous Cardiac Signal Generation From A Video Of A Subject Being Monitored For Cardiac Function”, U.S. patent application Ser. No. 13/871,766, by Kyal et al.
“Determining Cardiac Arrhythmia From A Video Of A Subject Being Monitored For Cardiac Function”, U.S. patent application Ser. No. 13/532,128, by Mestha et al.
What is disclosed is a system and method for detecting an arrhythmic cardiac event from a cardiac signal obtained from a subject being monitored for cardiac function. One embodiment of the present method involves the following. First, a plurality of different cardiac signals are received. The cardiac signals can be received from a plurality of cardiac function monitoring devices or may be a videoplethysmographic signal obtained from processing overlapping batches of video image frames of the subject where a plethysmographic signal can be registered by at least one channel of the video device used to capture those video images. The received cardiac signals are then transformed into frequency domain signals using, for example, a Fast Fourier Transform (FFT) or a Discrete Cosine Transform (DCT) method. The frequency domain signals are then changed such that dominant frequencies in each of the signals are substantially aligned. This alignment forms a matrix of feature vectors, whose dimension is then reduced by performing Principal Component Analysis. The reduced matrix is an eigen feature matrix. The eigen feature matrix is used to train a classifier. The classifier can be, for instance, a Support Vector Machine (SVM) classifier, or a neural network, as are commonly understood. Once the classifier has been trained, one or more cardiac signals from a subject are received. The subject's cardiac signal is transformed to a frequency domain signal which is changed in a manner such that a dominant frequency in that cardiac signal are substantially aligned with a dominant frequency of the signals used to train the classifier. If the subject's frequency domain signal has a frequency bandwidth which is different than that of the signals used to train the classifier then the subject's frequency domain signal is changed to have a same frequency bandwidth and a same signal length as those training signals. The subject's frequency domain signal is then provided as a new feature vector to the classifier. The classifier then uses the new feature vector to classify the subject as having one of: an arrhythmic event, and a non-arrhythmic event. Many features and advantages of the present method will become readily apparent from the following detailed description and accompanying drawings.
The foregoing and other features and advantages of the subject matter disclosed herein will be made apparent from the following detailed description taken in conjunction with the accompanying drawings, in which:
What is disclosed is a system and method for detecting an arrhythmic cardiac event from a cardiac signal obtained from a subject being monitored for cardiac function.
A “subject” refers to a living being. Although the term “person” or “patient” may be used throughout this disclosure, it should be appreciated that the subject may be something other than a human such as, for example, a primate. Therefore, the use of such terms is not to be viewed as limiting the scope of the appended claims strictly to humans.
A “cardiac signal” is a signal which relates to the function of the subject's heart. Cardiac signals can be, for instance, a videoplethysmographic (VPG) signal extracted from a time-series signal obtained from processing overlapping batches of image frames captured by an imaging device such as, a video camera, a RGB camera, a multi-spectral or hyperspectral imaging system, and a hybrid device comprising any combination thereof. Such imaging devices typically have a plurality of outputs where the captured images are obtained on a per-channel basis.
A “feature vector” contains features of interest obtained by analyzing the cardiac signal. A feature of interest may be one or more aspects of, for instance, a frequency domain version of the cardiac signal. In addition to this, features of interest may be one or more higher order statistical quantities obtained by analyzing a set of peak-to-peak intervals of a cardiac signal with respect to any of: a mean, standard deviation, skewness, and kurtosis. Features of interest may be one or more heart rate variability metrics obtained by analyzing a cardiac signal with respect to any of: a Standard Deviation of RR Intervals (SDRR), Root Mean Square of Successive RR Difference (RMSSD), Proportion of NN or RR interval exceeding 50 milliseconds (pNN50), Shannon Entropy (ShE), Standard Deviation 1 (SD1), Standard Deviation 2 (SD2), Pulse Harmonic Strength (PHS), and Normalized Pulse Harmonic Strength (NPHS). Patient information and medical histories may further be associated with various cardiac signals and features of interest.
A “matrix of feature vectors” is a dimensionality reduced matrix comprising an eigen feature matrix. In one embodiment, the feature vector is obtained by having performed a Principal Component Analysis (PCA) method. The matrix of features vectors is used to train a classifier.
A “classifier” is an artificial intelligence system which functions to map feature spaces to labels which define each space. Once trained, the classifier receives a new set of features for an unclassified event and assigns a label to those features. Classifiers can take any of a variety of forms including a Support Vector Machine (SVM), a neural network, a Bayesian network, a Logistic regression, Naïve Bayes, Randomized Forests, Decision Trees and Boosted Decision Trees, K-nearest neighbor, and a Restricted Boltzmann Machine (RBM), as are understood in the machine learning arts. For an in-depth discussion, the reader is directed to any of a wide variety of texts on classifiers, including: “Foundations of Machine Learning”, by Mehryar Mohri, Afshin Rostamizadeh, Ameet Talwalkar, MIT Press (2012), ISBN-13: 978-0262018258, and “Design and Analysis of Learning Classifier Systems: A Probabilistic Approach”, by Jan Drugowitsch, Springer (2008), ISBN-13: 978-3540798651, both of which are incorporated herein in their entirety by reference. The classifier uses the subject's feature vector to classify the subject as having an arrhythmic or non-arrhythmic event.
An “arrhythmic event” refers to a cardiac arrhythmia, (also known as cardiac dysrhythmia), which is an irregular heartbeat. One common cardiac arrhythmia is atrial fibrillation. Other arrhythmias include ventricular tachycardia, sinus tachycardia, sinus bradycardia, and ventricular premature contraction, as are generally understood in the medical arts. The reader is referred to any of a variety of medical texts for a thorough discussion on the different cardiac arrhythmias which comprise an arrhythmic event.
A “non-arrhythmic event” refers to the heart's normal sinus rhythm (SR).
A “frequency domain signal” is a signal obtained by performing any of: a non-parametric spectral density estimation on a cardiac signal, a parametric spectral density estimation on a cardiac signal, a Fourier Transform, a Wavelet Transform, and/or a Discrete Cosine Transform. Methods for transforming a signal into a frequency domain signal are well established in the signal processing arts.
“Receiving cardiac signals” is intended to be widely construed and includes retrieving, capturing, acquiring, or otherwise obtaining cardiac signals for processing in accordance with the methods disclosed herein. The cardiac signals can be retrieved from a memory or storage of the device used to acquire those signals. Cardiac signals can be retrieved from a media such as a CDROM or DVD, or can be received from a remote device over a network. Cardiac signals may be downloaded from a web-based system or application which makes such signals available for processing.
“Changing frequency domain signals” means to shift, scale, or otherwise modify the frequency domain signals such that a dominant frequency in each of those signals are substantially aligned. In one embodiment, the frequency domain signals are changed to have the same sampling frequency and the same signal length as the signals used to train the classifier. The aligned signals form feature vectors which, in turn, are used to form the matrix of feature vectors.
It should be appreciated that the steps of “transforming”, “changing”, “training”, “providing”, “processing”, and the like, as used herein, include the application of various signal processing and mathematical operations as applied to data and signals, according to any specific context or for any specific purpose. The terms in this Detailed Description and claims include any activity, in hardware or software, having the substantial effect of the mathematical or signal-processing action. It should be appreciated that such steps may be facilitated or otherwise effectuated by a microprocessor executing machine readable program instructions retrieved from a memory or storage device.
Reference is now being made to the flow diagram of
At step 302, receive a plurality of different cardiac signals.
At step 304, transform each of the cardiac signals to frequency domain signals.
At step 306, change the frequency domain signals such that a dominant frequency in each of the signals are substantially aligned to form a matrix of feature vectors.
At step 308, train a classifier with the matrix of feature vectors.
At step 310, receive a cardiac signal of a subject. One embodiment of cardiac signals being obtained for a subject is shown in
At step 312, transform the subject's cardiac signal to a frequency domain signal.
At step 314, change the subject's frequency domain signal such that a dominant frequency in the signal is substantially aligned with the dominant frequency of the signals used to train the classifier.
Reference is now being made to the flow diagram of
At step 316, provide the subject's changed frequency domain signal as a new feature vector to the trained classifier, the classifier using the new feature vector to classify the subject as having one of: an arrhythmic event, and a non-arrhythmic event.
At step 318, communicate the subject's classification to a display device. An alert signal may be initiated and/or a signal sent to a medical professional.
At step 320, a determination is made whether to receive and process another cardiac signal from the subject. If so, then processing repeats with respect to node B wherein, at step 310, a next cardiac signal of the subject is received. Processing repeats in a similar manner for this next received cardiac signal. Processing can repeat for each next received cardiac signal on a real-time and on a continuous basis. If, at step 320, no more of the subject's cardiac signals are to be processed then further processing stops.
It should be appreciated that the flow diagrams depicted herein are illustrative. One or more of the operations illustrated in the flow diagrams may be performed in a differing order. Other operations may be added, modified, enhanced, or consolidated. Variations thereof are intended to fall within the scope of the appended claims.
Reference is now being made to
Signal Receiver 501 acts as a buffer and receives n≧2 cardiac signals (collectively at 502) and stores the received cardiac signals to storage device 503. Signal Receiver 501 may be used for queuing additional signals, as needed including information about the received signals signal length, detected peaks, time/date information, and the like, and may further be configured to also store data, mathematical formulas and other representations to facilitate processing of received cardiac signals in accordance with the teachings hereof. Signal Transform Module 504 retrieves the stored cardiac signals and transforms each of the received cardiac signals to frequency domain signals. The frequency domain signals are then stored to storage device 503. Signal Alignment Module 505 retrieves the frequency domain signals and proceeds to change the frequency domain signals such that a dominant frequency in each of the training signals is substantially aligned. Alternatively, the dimensionality of the aligned signals is reduced. One method for reducing dimensionality is by Principle Component Analysis (PCA). The aligned signals are stored to storage device 503. Module 505 also provides as output a matrix of feature vectors 506 to Training Module 507. The Training Module uses the matrix of feature vectors to train Classifier 508. Although the Classifier is shown as being external to the Signal Processing System 500, in other embodiments, the Classifier is integrated with the other modules and processors of the system 500. Central Processor (CPU) 509 is provided generally to facilitate the functionality of any of the modules of the signal processing system 500. Microprocessor 509, operating alone or in conjunction with other processors, retrieves machine readable program instructions from Memory 510 and is configured to assist or otherwise facilitate the performance of any of the functionality of any of the components of system 500. Processor 509 further effectuates communication between various modules of system 500 and with workstation 521.
Video imaging device 511 acquires streaming video of the subject 100 of
Workstation 520 has a computer case 521 which houses various components such as a motherboard with a processor and memory, a network card, a video card, a hard drive capable of reading/writing to machine readable media 522 such as a floppy disk, optical disk, CD-ROM, DVD, magnetic tape, and the like, and other software and hardware needed to perform the functionality of a computer workstation. The workstation further includes a display device 523, such as a CRT, LCD, or touchscreen device, for displaying information, video, distances, clusters, features of interest, computed values, medical information, results, and the like, which are produced or are otherwise generated by any of the block modules of system 500. A user can view any of that information and make a selection from menu options displayed thereon. Keyboard 524 and mouse 525 effectuate a user input or selection as needed. The workstation implements a database 526 wherein records (collectively at 527) are stored, manipulated, and retrieved in response to a query. Such records, in various embodiments, take the form of patient medical history stored in association with information identifying the patient along with medical information. Although the database is shown as an external device, the database may be internal to the workstation mounted, for example, on a hard disk therein.
The workstation has an operating system and other specialized software configured to display alphanumeric values, menus, scroll bars, dials, slideable bars, pull-down options, selectable buttons, and the like, for entering, selecting, modifying, and accepting information needed for training a classifier and using that classifier in accordance with the methods disclosed herein. A user or technician may use the workstation to analyze the cardiac signals, identify features of interest, associate various features of interest with different cardiac events, train the classifier, set various parameters, select cardiac signals for processing, and/or use the workstation to facilitate the functionality of any of the modules and processing units of Signal Processing System 500. User input and selections may be stored/retrieved in either of storage devices 503 and 526. Default settings and initial parameters can also be retrieved from any of these storage devices. A user may adjust various parameters being utilized or dynamically adjust settings in real-time. The alert signal generated may be received and viewed by the workstation 520 and/or communicated to one or more remote devices over network 528.
Although shown as a desktop computer, it should be appreciated that the workstation 520 can be a laptop, mainframe, or a special purpose computer such as an ASIC, circuit, or the like. The embodiment of the workstation is illustrative and may include other functionality known in the arts. Any of the components of the workstation may be placed in communication with any of the modules of system 500 or any devices placed in communication therewith. Moreover, any of the modules of system 500 can be placed in communication with storage devices 503 and 526 and/or computer readable media 522 and may store/retrieve therefrom data, variables, parameters, signals, records, mathematical functions, and machine readable/executable program instructions, as needed to perform their intended functions. Moreover, any of the modules of system 500 may be placed in communication with one or more remote devices over network 528.
It should be appreciated that some or all of the functionality performed by any of the modules or processing units of the system of
The methods disclosed herein were applied on 11 patients using data obtained through a Xerox-University of Rochester Medical Center (URMC) Research Agreement. Patients with atrial fibrillation (AF) had undergone electrocardioversion to recover back to normal sinus rhythm (SR). The setup included ECG recording via a Holter monitor and video recording with a standard RGB camera at the same time. We were also provided with the manual annotations i.e. whether the beat is normal, AF beat or Ventricular Premature Contraction (VPC), etc. Initially a facial region of interest (ROI) was selected manually from the first video frame. Motion tracking of the ROI was used to automatically select multiple video segments of a duration of 15 seconds with minimum or no motion. These video segments were then pre-processed to extract a VPG signal for the subject. This was followed by calculating the power spectral density (PSD) for all signals. We detected a total of 408 segments [143 (AF)+265 (SR)] from the 11 patients followed by PCA and then training an SVM classifier. ECG RR-intervals and manual annotations for the corresponding video segments were found at the same time.
The teachings hereof can be implemented in hardware or software using any known or later developed systems, structures, devices, and/or software by those skilled in the applicable art without undue experimentation from the functional description provided herein with a general knowledge of the relevant arts. One or more aspects of the methods described herein are intended to be incorporated in an article of manufacture. The article of manufacture may be shipped, sold, leased, or otherwise provided separately either alone or as part of a product suite or a service.
The above-disclosed and other features and functions, or alternatives thereof, may be desirably combined into other different systems or applications. Presently unforeseen or unanticipated alternatives, modifications, variations, or improvements may become apparent and/or subsequently made by those skilled in this art which are also intended to be encompassed by the following claims. The teachings of any publications referenced herein are hereby incorporated by reference in their entirety.