The main limitation of radiofrequency (RF) ablation numerical simulations reported in the literature is their failure to provide statistical results based on the statistical variability of tissue thermal-electrical parameters. This work developed an efficient probabilistic approach to hepatic RF ablation in order to statistically evaluate the effect of four thermal-electrical properties of liver tissue on the uncertainty of the ablation zone dimensions: thermal conductivity, specific heat, blood perfusion and electrical conductivity. A deterministic thermal-electrical finite element model of a monopolar electrode inserted in the liver was coupled with the unscented transform method in order to obtain coagulation zone confidence intervals, probability and cumulative density functions. The coagulation zone volume, diameter and length were 10.96 cm(3), 2.17 cm and 4.08 cm, respectively (P < 0.01). Furthermore, a probabilistic sensitivity analysis showed that perfusion and thermal conductivity account for >95% of the variability in coagulation zone volume, diameter and length.
One of the current shortcomings of radiofrequency (RF) tumor ablation is its limited performance in regions close to large blood vessels, resulting in high recurrence rates at these locations. Computer models have been used to determine tissue temperatures during tumor ablation procedures. To simulate large vessels, either constant wall temperature or constant convective heat transfer coefficient (h) have been assumed at the vessel surface to simulate convection. However, the actual distribution of the temperature on the vessel wall is non-uniform and time-varying, and this feature makes the convective coefficient variable. This paper presents a realistic time-varying model in which h is a function of the temperature distribution at the vessel wall. The finite-element method (FEM) was employed in order to model RF hepatic ablation. Two geometrical configurations were investigated. The RF electrode was placed at distances of 1 and 5 mm from a large vessel (10 mm diameter). When the ablation procedure takes longer than 1–2 min, the attained coagulation zone obtained with both time-varying h and constant h does not differ significantly. However, for short duration ablation (5–10 s) and when the electrode is 1 mm away from the vessel, the use of constant h can lead to errors as high as 20% in the estimation of the coagulation zone. For tumor ablation procedures typically lasting at least 5 min, this study shows that modeling the heat sink effect of large vessels by applying constant h as a boundary condition will yield precise results while reducing computational complexity. However, for other thermal therapies with shorter treatment using a time-varying h may be necessary.
The surface electromyographic (SEMG) signal is very convenient for prosthesis control because it is non-invasively acquired and intrinsically related to the user's intention. This work presents a feature extraction and pattern classification algorithm for estimation of the intended knee joint angle from SEMG signals acquired using two sets of electrodes placed on the upper leg. The proposed algorithm uses a combination of time-domain and frequency-domain approaches for feature extraction (signal amplitude histogram and auto-regressive coefficients, respectively), a self-organizing map for feature projection and a Levenberg-Marquardt multi-layer perceptron neural network for pattern classification. The new algorithm was quantitatively compared with the method proposed by Wang et al (2006 Med. Biol. Eng. Comput. 44 865-72), which uses wavelet packet feature extraction, principal component analysis and a multi-layer perceptron neural classifier. The proposed method provided lower error-to-signal percentage and peak error amplitudes, higher correlation and fewer error events. The algorithm presented in this work may be useful as part of a myoelectric controller for active leg prostheses designed for transfemoral amputees.
This paper presents a hybrid adaptive algorithm for the compression of surface electromyographic (S-EMG) signals recorded during isometric and/or isotonic contractions. This technique is useful for minimizing data storage and transmission requirements for applications where multiple channels with high bandwidth data are digitized, such as telemedicine applications. The compression algorithm proposed in this work uses a discrete wavelet transform for spectral decomposition and an intelligent dynamic bit allocation scheme implemented by an approach using the Kohonen layer, which improves the bit allocation for sections of the S-EMG with different characteristics. Finally, data and overhead information are packed by entropy coding. The results for the compression of isometric EMG signals showed that this algorithm has a better performance than standard wavelet compression algorithms presented in the literature (presenting a decrease of at least 5% in per cent residual difference (PRD) for the same compression ratio), and a performance that is comparable with the performance of algorithms based on an embedded zero-tree wavelet. For isotonic EMG signals, its performance is better than the performance of the algorithms based on embedded zero-tree wavelets (presenting a decrease in PRD of about 3.6% for the same compression ratios, in the useful compression range).
17 September 2015 in BioMedical Engineering OnLine
In surface electromyography (surface EMG, or S-EMG), conduction velocity (CV) refers to the velocity at which the motor unit action potentials (MUAPs) propagate along the muscle fibers, during contractions. The CV is related to the type and diameter of the muscle fibers, ion concentration, pH, and firing rate of the motor units (MUs). The CV can be used in the evaluation of contractile properties of MUs, and of muscle fatigue. The most popular methods for CV estimation are those based on maximum likelihood estimation (MLE). This work proposes an algorithm for estimating CV from S-EMG signals, using digital image processing techniques. The proposed approach is demonstrated and evaluated, using both simulated and experimentally-acquired multichannel S-EMG signals. We show that the proposed algorithm is as precise and accurate as the MLE method in typical conditions of noise and CV. The proposed method is not susceptible to errors associated with MUAP propagation direction or inadequate initialization parameters, which are common with the MLE algorithm. Image processing -based approaches may be useful in S-EMG analysis to extract different physiological parameters from multichannel S-EMG signals. Other new methods based on image processing could also be developed to help solving other tasks in EMG analysis, such as estimation of the CV for individual MUs, localization and tracking of innervation zones, and study of MU recruitment strategies.