Systems and methods for maternal uterine activity detection

Abstract

A method includes receiving acoustic inputs; generating signal channels from the acoustic inputs; pre-processing data in the signal channels; extracting S 1 -S 2 peaks from the pre-processed data; removing artifacts and outliers from the S 1 -S 2 peaks; generating S 1 -S 2 signal channels based on the S 1 -S 2 peaks in the pre-processed signal channels; selecting two or more of the S 1 -S 2 signal channels; and combining the selected two or more S 1 -S 2 signal channels to produce an acoustic uterine monitoring signal.

What is claimed is: 
     
         1 . A computer-implemented method, comprising:
 receiving, by at least one computer processor, a plurality of raw acoustic inputs,
 wherein each of the raw acoustic inputs being received from a corresponding one of a plurality of acoustic sensors, and 
 wherein each of the plurality of acoustic sensors is positioned so as to measure a respective one of the raw acoustic inputs of a pregnant human subject; 
   generating, by the at least one computer processor, a plurality of signal channels from the plurality of raw acoustic inputs,
 wherein the plurality of signal channels comprises at least three signal channels; 
   pre-processing, by the at least one computer processor, respective signal channel data of each of the signal channels to produce a plurality of pre-processed signal channels,
 wherein each of the pre-processed signal channels comprises respective pre-processed signal channel data; 
   extracting, by the at least one computer processor, a respective plurality of S 1 -S 2  peaks from the pre-processed signal channel data of each of the pre-processed signal channels to produce a plurality of S 1 -S 2  peak data sets,   wherein each of the S 1 -S 2  peak data sets comprises a respective plurality of S 1 -S 2  peaks;   removing, by the at least one computer processor, from the plurality of S 1 -S 2  peak data sets, at least one of: (a) at least one signal artifact or (b) at least one outlier data point,
 wherein the at least one signal artifact is one of a movement-related artifact or a baseline artifact; 
   replacing, by the at least one computer processor, the at least one signal artifact, the at least one outlier data point, or both, with at least one statistical value determined based on a corresponding one of the S 1 -S 2  peak data sets from which the at least one signal artifact, the at least one outlier data point, or both was removed, to produce a plurality of interpolated S 1 -S 2  peak data sets;   generating, by the at least one computer processor, a respective S 1 -S 2  signal data set for a respective S 1 -S 2  signal channel at a predetermined sampling rate based on each respective interpolated S 1 -S 2  peak data set to produce a plurality of S 1 -S 2  signal channels;   selecting, by the at least one computer processor, at least one first selected S 1 -S 2  signal channel and at least one second selected S 1 -S 2  signal channel from the plurality of S 1 -S 2  signal channels based on at least one correlation between (a) a respective S 1 -S 2  signal data set of at least one first particular S 1 -S 2  signal channel and (b) a respective S 1 -S 2  signal data set of at least one second particular S 1 -S 2  signal channel; and   generating, by the at least one computer processor, acoustic uterine monitoring data representative of an acoustic uterine monitoring signal based on at least the respective S 1 -S 2  signal data set of the first selected S 1 -S 2  signal channel and the respective S 1 -S 2  signal data set of the second selected S 1 -S 2  signal channel.   
     
     
         2 . The computer-implemented method of  claim 1 , further comprising:
 sharpening, by the at least one computer processor, the acoustic uterine monitoring data to produce a sharpened acoustic uterine monitoring signal.   
     
     
         3 . The computer-implemented method of  claim 2 , further comprising:
 determining, by the at least one computer processor, whether the acoustic uterine monitoring data is calculated based on based on a selected one of the acoustic uterine monitoring signal channels that is a corrupted acoustic uterine monitoring signal channel,   wherein the step of sharpening, by the at least one computer processor, the acoustic uterine monitoring data is omitted if the acoustic uterine monitoring data is calculated based on the selected one of the acoustic uterine monitoring signal channels that is the corrupted acoustic uterine signal monitoring channel.   
     
     
         4 . The computer-implemented method of  claim 2 , further comprising:
 post-processing the sharpened acoustic monitoring signal data to produce a post-processed acoustic uterine monitoring signal.   
     
     
         5 . The computer-implemented method of  claim 2 , wherein the step of sharpening, by the at least one computer processor, the acoustic uterine monitoring data comprises:
 identifying a set of peaks in the acoustic uterine monitoring signal data;   determining a prominence of each of the peaks;   removing, from the set of peaks, peaks having a prominence that is less than at least one threshold prominence value;   calculating a mask based on remaining peaks of the set of peaks;   smoothing the mask based on a moving average window to produce a smoothed mask; and   adding the smoothed mask to the acoustic uterine monitoring signal data to produce the sharpened acoustic uterine monitoring signal data.   
     
     
         6 . The computer-implemented method of  claim 5 , wherein the at least one threshold prominence value includes at least one threshold prominence value selected from a group consisting of an absolute prominence value and a relative prominence value calculated based on a maximal prominence of the peaks in the set of peaks. 
     
     
         7 . The computer-implemented method of  claim 5 , wherein the mask includes zero values outside areas of the remaining peaks and nonzero values inside areas of the remaining peaks,
 wherein the nonzero values are calculated based on a Gaussian function.   
     
     
         8 . The computer-implemented method of  claim 1 , wherein the at least one filtering step of the pre-processing step includes applying at least one filter selected from a group consisting of a DC removal filter, a powerline filter, and a high pass filter. 
     
     
         9 . The computer-implemented method of  claim 1 , wherein the step of extracting, by the at least one computer processor, the respective plurality of S 1 -S 2  peaks comprises:
 receiving a set of maternal ECG peaks for the pregnant human subject; and 
 identifying S 1 -S 2  peaks in each of the pre-processed signal channels within a predetermined time window before and after each of the maternal ECG peaks in the set of maternal ECG peaks as a maximum absolute value in each of the pre-processed signal channels within the predetermined time window. 
 
     
     
         10 . The computer-implemented method of  claim 1 , wherein the step of removing at least one of a signal artifact or an outlier data point comprises removing at least one movement artifact by a process comprising:
 identifying at least one corrupted peak in one of the plurality of S 1 -S 2  peak data sets based on the at least one corrupted peak having an inter-peaks root mean square value that is greater than a threshold; and   replacing the corrupted peak with a median value, wherein the median value is either a local median or a global median.   
     
     
         11 . The computer-implemented method of  claim 1 , wherein the step of removing at least one of a signal artifact or an outlier data point comprises removing at least one baseline artifact by a process comprising:
 identifying a change point in S 1 -S 2  peaks in one of the plurality of S 1 -S 2  peak data sets;   subdividing the one of the plurality of S 1 -S 2  peak data sets into a first portion located prior to the change point and a second portion located subsequent to the change point;   determining a first root-mean-square value for the first portion;   determining a second root-mean-square value for the second portion;   determining an equalization factor based on the first root-mean-square value and the second root-mean-square value; and   modifying the first portion by multiplying S 1 -S 2  peaks in the first portion by the equalization factor.   
     
     
         12 . The computer-implemented method of  claim 1 , wherein the step of removing at least one of a signal artifact or an outlier point comprises removing at least one outlier in accordance with a Grubbs test for outliers. 
     
     
         13 . The computer-implemented method of  claim 1 , wherein the step of generating a respective S 1 -S 2  data set based on each respective S 1 -S 2  peak data set comprises interpolating between the S 1 -S 2  peaks of each respective S 1 -S 2  peak data set, and wherein the interpolating between the S 1 -S 2  peaks comprises interpolating using an interpolation algorithm that is selected from a group consisting of a cubic spline interpolation algorithm and a shape-preserving piecewise cubic interpolation algorithm. 
     
     
         14 . The computer-implemented method of  claim 1 , wherein the step of selecting at least one first one of the S 1 -S 2  signal channels and at least one second one of the S 1 -S 2  signal channels from the plurality of S 1 -S 2  signal channels comprises:
 identifying contact issues in each of the plurality of S 1 -S 2  signal channels; 
 determining a percentage of prior intervals in which each of the plurality of S 1 -S 2  signal channels experienced contact issues; 
 selecting candidate S 1 -S 2  signal channels from the plurality of S 1 -S 2  signal channels based on the percentage of prior intervals in which each of the plurality of S 1 -S 2  signal channels experienced contact issues; 
 grouping the selected candidate S 1 -S 2  signal channels into a plurality of couples, wherein each of the couples includes two of the selected candidate S 1 -S 2  channels that are independent from one another; 
 calculating a correlation value of each of the couples; 
 identifying at least one of the couples having a correlation value that exceeds a threshold correlation value; 
 selecting, as the selected at least one first one of the S 1 -S 2  signal channels, the first candidate S 1 -S 2  signal channel of each of the identified at least one of the couples; and 
 selecting, as the selected at least one second one of the S 1 -S 2  signal channels, the second candidate S 1 -S 2  signal channel of each of the identified at least one of the couples. 
 
     
     
         15 . The computer-implemented method of  claim 1 , wherein the step of calculating the acoustic uterine monitoring signal comprises calculating the acoustic uterine monitoring signal that is, at each point in time, a predetermined percentile of all of the selected at least one first one of the S 1 -S 2  signal channels and the selected at least one second one of the S 1 -S 2  signal channels at that point in time. 
     
     
         16 . The computer-implemented method of  claim 15 , wherein the predetermined percentile is an 80th percentile. 
     
     
         17 . The computer-implemented method of  claim 1 , wherein the statistical value is one of a local median, a global median, or a mean.