Jizhou Li

Dynamic measurement of pennation angle of gastrocnemius muscles during contractions based on ultrasound imaging

BackgroundMuscle fascicle pennation angle (PA) is an important parameter related to musculoskeletal functions, and ultrasound imaging has been widely used for measuring PA, but manually and frame by frame in most cases. We have earlier reported an automatic method to estimate aponeurosis orientation based on Gabor transform and Revoting Hough Transform (RVHT).MethodsIn this paper, we proposed a method to estimate the overall orientation of muscle fascicles in a region of interest, in order to complete computing the orientation of the other side of the pennation angle, but the side found by RVHT. The measurements for orientations of both fascicles and aponeurosis were conducted in each frame of ultrasound images, and then the dynamic change of pennation angle during muscle contraction was obtained automatically. The method for fascicle orientation estimation was evaluated using synthetic images with different noise levels and later on 500 ultrasound images of human gastrocnemius muscles during isometric plantarflexion.ResultsThe muscle fascicle orientations were also estimated manually by two operators. From the results it’s found that the proposed automatic method demonstrated a comparable performance to the manual method.ConclusionsWith the proposed methods, ultrasound measurement for muscle pennation angles can be more widely used for functional assessment of muscles.

The Sensitive and Efficient Detection of Quadriceps Muscle Thickness Changes in Cross-Sectional Plane Using Ultrasonography: A Feasibility Investigation

As a direct determinant parameter to quantify muscle activity, the muscle thickness (MT) has been investigated in many aspects and for various purposes. Ultrasonography (US) is a promising modality to detect muscle morphological changes during contractions since it is portable, noninvasive, and real time. However, there are few reports on sensitive and efficient estimation of changes of MT in a cross-sectional plane. In this feasibility investigation, we proposed a coarse-to-fine method based on a compressive-tracking algorithm for estimation of MT changes during an example task of isometric knee extension using ultrasound images. The sensitivity and efficiency are evaluated with 1920 US images from quadriceps muscle (QM) in eight subjects. The detection results were compared with those obtained from both traditional manual measurement and the well known normalized cross-correlation method, and the effect of the size of tracking window on detection performance was evaluated as well. It is demonstrated that the proposed method agrees well with the manual measurement. Meanwhile, it is not only sensitive to relatively small changes of MT but also computationally efficient.

Estimation and visualization of longitudinal muscle motion using ultrasonography: A feasibility study

Ultrasonography is a convenient and widely used technique to look into the longitudinal muscle motion as it is radiation-free and real-time. The motion of localized parts of the muscle, disclosed by ultrasonography, spatially reflects contraction activities of the corresponding muscles. However, little attention was paid to the estimation of longitudinal muscle motion, especially towards estimation of dense deformation field at different depths under the skin. Yet fewer studies on the visualization of such muscle motion or further clinical applications were reported in the literature. A primal-dual algorithm was used to estimate the motion of gastrocnemius muscle (GM) in longitudinal direction in this study. To provide insights into the rules of longitudinal muscle motion, we proposed a novel framework including motion estimation, visualization and quantitative analysis to interpret synchronous activities of collaborating muscles with spatial details. The proposed methods were evaluated on ultrasound image sequences, captured at a rate of 25 frames per second from eight healthy subjects. In order to estimate and visualize the GM motion in longitudinal direction, each subject was asked to perform isometric plantar flexion twice. Preliminary results show that the proposed visualization methods provide both spatial and temporal details and they are helpful to study muscle contractions. One of the proposed quantitative measures was also tested on a patient with unilateral limb dysfunction caused by cerebral infarction. The measure revealed distinct patterns between the normal and the dysfunctional lower limb. The proposed framework and its associated quantitative measures could potentially be used to complement electromyography (EMG) and torque signals in functional assessment of skeletal muscles.

Deconvolution of poissonian images with the PURE-LET approach

We propose a non-iterative image deconvolution algorithm for data corrupted by Poisson noise. Many applications involve such a problem, ranging from astronomical to biological imaging. We parametrize the deconvolution process as a linear combination of elementary functions, termed as linear expansion of thresholds (LET). This parametrization is then optimized by minimizing a robust estimate of the mean squared error, the “Poisson unbiased risk estimate (PURE)”. Each elementary function consists of a Wiener filtering followed by a pointwise thresholding of undecimated Haar wavelet coefficients. In contrast to existing approaches, the proposed algorithm merely amounts to solving a linear system of equations which has a fast and exact solution. Simulation experiments over various noise levels indicate that the proposed method outperforms current state-of-the-art techniques, in terms of both restoration quality and computational time.

Fast and accurate three-dimensional point spread function computation for fluorescence microscopy

The point spread function (PSF) plays a fundamental role in fluorescence microscopy. A realistic and accurately calculated PSF model can significantly improve the performance in 3D deconvolution microscopy and also the localization accuracy in single-molecule microscopy. In this work, we propose a fast and accurate approximation of the Gibson-Lanni model, which has been shown to represent the PSF suitably under a variety of imaging conditions. We express the Kirchhoffs integral in this model as a linear combination of rescaled Bessel functions, thus providing an integral-free way for the calculation. The explicit approximation error in terms of parameters is given numerically. Experiments demonstrate that the proposed approach results in a significantly smaller computational time compared with current state-of-the-art techniques to achieve the same accuracy. This approach can also be extended to other microscopy PSF models.

An attempt to bridge muscle architecture dynamics and its instantaneous rate of force development using ultrasonography.

Muscle force output is an essential index in rehabilitation assessment or physical exams, and could provide considerable insights for various applications such as load monitoring and muscle assessment in sports science or rehabilitation therapy. Besides direct measurement of force output using a dynamometer, electromyography has earlier been used in several studies to quantify muscle force as an indirect means. However, its spatial resolution is easily compromised as a summation of the action potentials from neighboring motor units of electrode site. To explore an alternative method to indirectly estimate the muscle force output, and with better muscle specificity, we started with an investigation on the relationship between architecture dynamics and force output of triceps surae. The muscular architecture dynamics is captured in ultrasonography sequences and estimated using a previously reported motion estimation method. Then an indicator named as the dorsoventrally averaged motion profile (DAMP) is employed. The performance of force output is represented by an instantaneous version of the rate of force development (RFD), namely I-RFD. From experimental results on ten normal subjects, there were significant correlations between the I-RFD and DAMP for triceps surae, both normalized between 0 and 1, with the sum of squares error at 0.0516±0.0224, R-square at 0.7929±0.0931 and root mean squared error at 0.0159±0.0033. The statistical significance results were less than 0.01. The present study suggested that muscle architecture dynamics extracted from ultrasonography during contraction is well correlated to the I-RFD and it can be a promising option for indirect estimation of muscle force output.

Real-time detection of muscle thickness changes during isometric contraction from ultrasonography: A feasibility study

Ultrasonography is a convenient and widely-used technique to look into the muscle contraction as it is non-invasive and real-time. Muscle thickness has been investigated in many aspects and for various purposes. However, there are few studies on automatic estimation of changes of muscle thickness in cross-sectional plane during contraction. In this study we proposed a coarse-to-fine method based on compressive tracking algorithm for real-time estimation of muscle thickness during isometric contraction on ultrasound images. The real-time performance and precision are investigated with data from quadriceps muscle in five subjects. The detection results were compared with those obtained from both cross-correlation algorithm and manual measurement. It is shown that the proposed method agrees well with the manual measurement and outperforms the cross-correlation method in the sense of both accuracy and computation cost.

PURE-LET Image Deconvolution

We propose a non-iterative image deconvolution algorithm for data corrupted by Poisson or mixed Poisson-Gaussian noise. Many applications involve such a problem, ranging from astronomical to biological imaging. We parameterize the deconvolution process as a linear combination of elementary functions, termed as linear expansion of thresholds. This parameterization is then optimized by minimizing a robust estimate of the true mean squared error, the Poisson unbiased risk estimate. Each elementary function consists of a Wiener filtering followed by a pointwise thresholding of undecimated Haar wavelet coefficients. In contrast to existing approaches, the proposed algorithm merely amounts to solving a linear system of equations, which has a fast and exact solution. Simulation experiments over different types of convolution kernels and various noise levels indicate that the proposed method outperforms the state-of-the-art techniques, in terms of both restoration quality and computational complexity. Finally, we present some results on real confocal fluorescence microscopy images and demonstrate the potential applicability of the proposed method for improving the quality of these images.

Pure-let deconvolution of 3D fluorescence microscopy images

Three-dimensional (3D) deconvolution microscopy is very effective in improving the quality of fluorescence microscopy images. In this work, we present an efficient approach for the deconvolution of 3D fluorescence microscopy images based on the recently developed PURE-LET algorithm. By combining multiple Wiener filtering and wavelet denoising, we parametrize the deconvolution process as a linear combination of elementary functions. Then the Poisson unbiased risk estimate (PURE) is used to obtain the optimal coefficients. The proposed approach is non-iterative and outperforms existing techniques (usually, variants of Richardson-Lucy algorithm) both in terms of computational efficiency and quality. We illustrate its effectiveness on both synthetic and real data.

A multi-frame optical flow spot tracker

Accurate and robust spot tracking is a necessary tool for quantitative motion analysis in fluorescence microscopy images. Few trackers however consider the underlying dynamics present in biological systems. For example, the collective motion of cells often exhibits both fast dynamics, i.e. Brownian motion, and slow dynamics, i.e. time-invariant stationary motion. In this paper, we propose a novel, multi-frame, tracker that exploits this stationary motion. More precisely, we first estimate the stationary motion and then use it to guide the spot tracker. We obtain the stationary motion by adapting a recent optical flow algorithm that relates one image to another locally using an all-pass filter. We perform this operation over all the image frames simultaneously and estimate a single, stationary optical flow. We compare the proposed tracker with two existing techniques and show that our approach is more robust to high noise and varying structure. In addition, we also show initial experiments on real microscopy images.