Antonio Fratini

Relevance of motion artifact in electromyography recordings during vibration treatment

Electromyography readings (EMGs) from quadriceps of fifteen subjects were recorded during whole body vibration treatment at different frequencies (10-50 Hz). Additional electrodes were placed on the patella to monitor the occurrence of motion artifact, triaxial accelerometers were placed onto quadriceps to monitor motion. Signal spectra revealed sharp peaks corresponding to vibration frequency and its harmonics, in accordance with the accelerometer data. EMG total power was compared to that associated with vibration harmonics narrow bands, before and during vibration. On average, vibration associated power resulted in only 3% (+/-0.9%) of the total power prior to vibration and 29% (+/-13.4%) during vibration. Often, studies employ surface EMG to quantitatively evaluate vibration evoked muscular activity and to set stimulation frequency. However, previous research has not accounted for motion artifacts. The data presented in this study emphasize the need for the removal of motion artifacts, as they consistently affect RMS estimation, which is often used as a concise muscle activity index during vibrations. Such artifacts, rather unpredictable in amplitude, might be the cause of large inter-study differences and must be eliminated before analysis. Motion artifact filtering will contribute to thorough and precise interpretation of neuromuscular response to vibration treatment.

Muscle motion and EMG activity in vibration treatment

The aim of this study is to highlight the relationship between muscle motion, generated by whole body vibration, and the correspondent electromyographic (EMG) activity and to suggest a new method to customize the stimulation frequency. Simultaneous recordings of EMG and tri-axial accelerations of quadriceps rectus femoris from fifteen subjects undergoing vibration treatments were collected. Vibrations were delivered via a sinusoidal oscillating platform at different frequencies (10-45 Hz). Muscle motion was estimated by processing the accelerometer data. Large EMG motion artifacts were removed using sharp notch filters centred at the vibration frequency and its superior harmonics. EMG-RMS values were computed and analyzed before and after artifact suppression to assess muscular activity. Muscles acceleration amplitude increased with frequency. Muscle displacements revealed a mechanical resonant-like behaviour of the muscle. Resonance frequencies and dumping factors depended on subject. Moreover, RMS of artifact-free EMG was found well correlated (R(2)=0.82) to the actual muscle displacement, while the maximum of the EMG response was found related to the mechanical resonance frequency of muscle. Results showed that maximum muscular activity was found in correspondence to the mechanical resonance of the muscle itself. Assuming the hypothesis that muscle activation is proportional to muscle displacement, treatment optimization (i.e. to choose the best stimulation frequency) could be obtained by simply monitoring local acceleration (resonance), leading to a more effective muscle stimulation. Motion artifact produced an overestimation of muscle activity, therefore its removal was essential.

Bluetooth portable device for continuous ECG and patient motion monitoring during daily life

Continuous patient monitoring during daily life can provide valuable information to different medical specialties. Indeed, long recording of related cardiac signals such as ECG, respiration and also other information such as body motion can improve diagnosis and monitor the evolution of many widespread diseases. Key-issues for portable or even wearable biomedical devices are: power consumption, longterm sensors, comfortable wearing, easy and wireless connectivity. Within this scenario, is valuable to realize prototypes making use of novel electronic technologies and common available communication technologies to assess practical use of long-term personal monitoring and foster new ways to provide healthcare services. We realized a small, battery powered, portable monitor capable to record ECG and body three-axes acceleration and continuously wireless transmit to any Bluetooth device including PDA or cellular phone. The ECG front end offers ultra-high input impedance allowing use of dry, long-lasting electrodes such conductive rubbers or novel textile electrodes that can be embedded in clothes. A small size MEMS 3-axes accelerometer was also integrated. Patient monitor incorporate a microprocessor that controls 12-bit ADC of signals at programmable sampling frequencies (e.g. 100 Hz) and drives a Bluetooth module capable to reliable transmit real-time signals within 10 m range. All circuitry can be powered by a standard mobile phone like Ni-MH 3.6V battery that can sustain more than seven day continuous functioning utilizing the Bluetooth Sniff mode to reduce TX power. At the moment we are developing dedicated software to process data and to extract concise parameters valuable for medical studies.

2D-3D registration of CT vertebra volume to fluoroscopy projection: a calibrationmodel assessment

This study extends a previous research concerning intervertebral motion registration by means of 2D dynamic fluoroscopy to obtain a more comprehensive 3D description of vertebral kinematics. The problem of estimating the 3D rigid pose of a CT volume of a vertebra from its 2D X-ray fluoroscopy projection is addressed. 2D-3D registration is obtained maximising a measure of similarity between Digitally Reconstructed Radiographs (obtained from the CT volume) and real fluoroscopic projection. X-ray energy correction was performed. To assess the method a calibration model was realised a sheep dry vertebra was rigidly fixed to a frame of reference including metallic markers. Accurate measurement of 3D orientation was obtained via single-camera calibration of the markers and held as true 3D vertebra position; then, vertebra 3D pose was estimated and results compared. Error analysis revealed accuracy of the order of 0.1 degree for the rotation angles of about 1 mm for displacements parallel to the fluoroscopic plane, and of order of 10 mm for the orthogonal displacement.

An algorithm for FHR estimation from foetal phonocardiographic signals

The long-term foetal surveillance is often to be recommended. Hence, the fully non-invasive acoustic recording, through maternal abdomen, represents a valuable alternative to the ultrasonic cardiotocography. Unfortunately, the recorded heart sound signal is heavily loaded by noise, thus the determination of the foetal heart rate raises serious signal processing issues. In this paper, we present a new algorithm for foetal heart rate estimation from foetal phonocardiographic recordings. A filtering is employed as a first step of the algorithm to reduce the background noise. A block for first heart sounds enhancing is then used to further reduce other components of foetal heart sound signals. A complex logic block, guided by a number of rules concerning foetal heart beat regularity, is proposed as a successive block, for the detection of most probable first heart sounds from several candidates. A final block is used for exact first heart sound timing and in turn foetal heart rate estimation. Filtering and enhancing blocks are actually implemented by means of different techniques, so that different processing paths are proposed. Furthermore, a reliability index is introduced to quantify the consistency of the estimated foetal heart rate and, based on statistic parameters; [,] a software quality index is designed to indicate the most reliable analysis procedure (that is, combining the best processing path and the most accurate time mark of the first heart sound, provides the lowest estimation errors). The algorithm performances have been tested on phonocardiographic signals recorded in a local gynaecology private practice from a sample group of about 50 pregnant women. Phonocardiographic signals have been recorded simultaneously to ultrasonic cardiotocographic signals in order to compare the two foetal heart rate series (the one estimated by our algorithm and the other provided by cardiotocographic device). Our results show that the proposed algorithm, in particular some analysis procedures, provides reliable foetal heart rate signals, very close to the reference cardiotocographic recordings.

A comparison of denoising methods for X-ray fluoroscopic images

Fluoroscopic images exhibit severe signal-dependent quantum noise, due to the reduced X-ray dose involved in image formation, that is generally modelled as Poisson-distributed. However, image gray-level transformations, commonly applied by fluoroscopic device to enhance contrast, modify the noise statistics and the relationship between image noise variance and expected pixel intensity. Image denoising is essential to improve quality of fluoroscopic images and their clinical information content. Simple average filters are commonly employed in real-time processing, but they tend to blur edges and details. An extensive comparison of advanced denoising algorithms specifically designed for both signal-dependent noise (AAS, BM3Dc, HHM, TLS) and independent additive noise (AV, BM3D, K-SVD) was presented. Simulated test images degraded by various levels of Poisson quantum noise and real clinical fluoroscopic images were considered. Typical gray-level transformations (e.g. white compression) were also applied in order to evaluate their effect on the denoising algorithms. Performances of the algorithms were evaluated in terms of peak-signal-to-noise ratio (PSNR), signal-to-noise ratio (SNR), mean square error (MSE), structural similarity index (SSIM) and computational time. On average, the filters designed for signal-dependent noise provided better image restorations than those assuming additive white Gaussian noise (AWGN). Collaborative denoising strategy was found to be the most effective in denoising of both simulated and real data, also in the presence of image gray-level transformations. White compression, by inherently reducing the greater noise variance of brighter pixels, appeared to support denoising algorithms in performing more effectively.

Analysis and modelling of muscles motion during whole body vibration

The aim of the study is to characterize the local muscles motion in individuals undergoing whole body mechanical stimulation. In this study we aim also to evaluate how subject positioning modifies vibration dumping, altering local mechanical stimulus. Vibrations were delivered to subjects by the use of a vibrating platform, while stimulation frequency was increased linearly from 15 to 60 Hz. Two different subject postures were here analysed. Platform and muscles motion were monitored using tiny MEMS accelerometers; a contra lateral analysis was also presented. Muscle motion analysis revealed typical displacement trajectories: motion components were found not to be purely sinusoidal neither in phase to each other. Results also revealed a mechanical resonant-like behaviour at some muscles, similar to a second-order system response. Resonance frequencies and dumping factors depended on subject and his positioning. Proper mechanical stimulation can maximize muscle spindle solicitation, which may produce a more effective muscle activation.

PSD modifications of FHRV due to interpolation and CTG storage rate

Cardiotocographic data provide physicians information about foetal development and permit to assess conditions such as foetal distress. An incorrect evaluation of the foetal status can be of course very dangerous. To improve interpretation of cardiotocographic recordings, great interest has been dedicated to foetal heart rate variability spectral analysis. It is worth reminding, however, that foetal heart rate is intrinsically an uneven series, so in order to produce an evenly sampled series a zero-order, linear or cubic spline interpolation can be employed. This is not suitable for frequency analyses because interpolation introduces alterations in the foetal heart rate power spectrum. In particular, interpolation process can produce alterations of the power spectral density that, for example, affects the estimation of the sympatho-vagal balance (computed as low-frequency/high-frequency ratio), which represents an important clinical parameter. In order to estimate the frequency spectrum alterations of the foetal heart rate variability signal due to interpolation and cardiotocographic storage rates, in this work, we simulated uneven foetal heart rate series with set characteristics, their evenly spaced versions (with different orders of interpolation and storage rates) and computed the sympatho-vagal balance values by power spectral density. For power spectral density estimation, we chose the Lomb method, as suggested by other authors to study the uneven heart rate series in adults. Summarising, the obtained results show that the evaluation of SVB values on the evenly spaced FHR series provides its overestimation due to the interpolation process and to the storage rate. However, cubic spline interpolation produces more robust and accurate results.

Waveform type evaluation in congenital nystagmus

Congenital nystagmus is an ocular-motor disorder characterised by involuntary, conjugated and bilateral to and fro ocular oscillations. In this study a method to recognise automatically jerk waveform inside a congenital nystagmus recording and to compute foveation time and foveation position variability is presented. The recordings were performed with subjects looking at visual targets, presented in nine eye gaze positions; data were segmented into blocks corresponding to each gaze position. The nystagmus cycles were identified searching for local minima and maxima (SpEp sequence) in intervals centred on each slope change of the eye position signal (position criterion). The SpEp sequence was then refined using an adaptive threshold applied to the eye velocity signal; the outcome is a robust detection of each slow phase start point, fundamental to accurately compute some nystagmus parameters. A total of 1206 slow phases was used to compute the specificity in waveform recognition applying only the position criterion or adding the adaptive threshold; results showed an increase in negative predictive value of 25.1% using both features. The duration of each foveation window was measured on raw data or using an interpolating function of the congenital nystagmus slow phases; foveation time estimation less sensitive to noise was obtained in the second case.

Subject identification via ECG fiducial-based systems: influence of the type of QT interval correction.

Electrocardiography (ECG) has been recently proposed as biometric trait for identification purposes. Intra-individual variations of ECG might affect identification performance. These variations are mainly due to Heart Rate Variability (HRV). In particular, HRV causes changes in the QT intervals along the ECG waveforms. This work is aimed at analysing the influence of seven QT interval correction methods (based on population models) on the performance of ECG-fiducial-based identification systems. In addition, we have also considered the influence of training set size, classifier, classifier ensemble as well as the number of consecutive heartbeats in a majority voting scheme. The ECG signals used in this study were collected from thirty-nine subjects within the Physionet open access database. Public domain software was used for fiducial points detection. Results suggested that QT correction is indeed required to improve the performance. However, there is no clear choice among the seven explored approaches for QT correction (identification rate between 0.97 and 0.99). MultiLayer Perceptron and Support Vector Machine seemed to have better generalization capabilities, in terms of classification performance, with respect to Decision Tree-based classifiers. No such strong influence of the training-set size and the number of consecutive heartbeats has been observed on the majority voting scheme.