Conduction velocity of the uterine contraction in serial magnetomyogram (MMG) data: Event based simulation and validation

We propose a novel approach to calculate the conduction velocity (CV) of the uterine contraction bursts in magnetomyogram (MMG) signals measured using a multichannel SQUID array. For this purpose, we partition the sensor coordinates into four different quadrants and identify the contractile bursts using a previously proposed Hilbert-wavelet transform approach. If contractile burst is identified in more than one quadrant, we calculate the center of gravity (CoG) in each quadrant for each time point as the sum of the product of the sensor coordinates with the Hilbert amplitude of the MMG signals normalized by the sum of the Hilbert amplitude of the signals over all sensors. Following this we compute the delay between the CoGs of all (six) possible quadrant pairs combinations. As a first step, we validate this approach by simulating a stochastic model based on independent second-order autoregressive processes (AR2) and we divide them into 30 second disjoint windows and insert burst activity at specific time instances in preselected sensors. Also we introduce a lag of 5 ± 1 seconds between different quadrants. Using our approach we calculate the CoG of the signals in a quadrant. To this end, we compute the delay between CoGs obtained from different quadrants and show that our approach is able to reliably capture the delay incorporated in the model. We apply the proposed approach to 19 serial MMG data obtained from two subjects and show an increase in the CV as the subjects approached labor.