Description: We propose a novel, physics-based method for detecting multi-scale tubular features in ultrasound images. The detector is based on a Hessian-matrix eigenvalue method, but unlike previous work, our detector is guided by an optimal model of vessel-like structures with respect to the ultrasound-image formation process. Our method provides a voxel-wise probability map, along with estimates of the radii and orientations of the detected tubes. These results can then be used for further processing, including segmentation and enhanced volume visualization. Most Hessian-based algorithms, including the well-known Frangi filter, were developed for CTA or MRA; they implicitly assume symmetry about the vessel centerline. This is not consistent with ultrasound data. We overcome this limitation by introducing a novel filter that allows multi-scale estimation both with respect to the vessel's centerline and with respect to the vessel's border. We use manually-segmented ultrasound imagery from 35 patients to show that our method is superior to standard Hessian-based methods. We evaluate the performance of the proposed methods based on the sensitivity and specificity like measures, and finally demonstrate further applicability of our method to vascular ultrasound images of the carotid artery, as well as ultrasound data for abdominal aortic aneurysms.

Description: A single-task functional magnetic resonance imaging (fMRI) experiment may only partially highlight alterations to functional brain networks affected by a particular disorder. Multivariate analysis across multiple fMRI tasks may increase the sensitivity of fMRI-based diagnosis. Prior research using multi-task analysis in fMRI, such as those that use joint independent component analysis (jICA), has mainly assumed that brain activity patterns evoked by different tasks are independent. This may not be valid in practice. Here, we use sparsity, which is a natural characteristic of fMRI data in the spatial domain, and propose a joint sparse representation analysis (jSRA) method to identify common information across different functional subtraction (contrast) images in data from a multi-task fMRI experiment. Sparse representation methods do not require independence, or that the brain activity patterns be nonoverlapping. We use functional subtraction images within the joint sparse representation analysis to generate joint activation sources and their corresponding sparse modulation profiles. We evaluate the use of sparse representation analysis to capture individual differences with simulated fMRI data and with experimental fMRI data. The experimental fMRI data was acquired from 16 young (age: 19–26) and 16 older (age: 57–73) adults obtained from multiple speech comprehension tasks within subjects, where an independent measure (namely, age in years) can be used to differentiate between groups. Simulation results show that this method yields greater sensitivity, precision, and higher Jaccard indexes (which measures similarity and diversity of the true and estimated brain activation sources) than does the jICA method. Moreover, superiority of the jSRA method in capturing individual differences was successfully demonstrated using experimental fMRI data.

Description: Balanced steady-state free precession (SSFP) imaging suffers from irrecoverable signal losses, known as banding artifacts, in regions of large B0 field inhomogeneity. A common solution is to acquire multiple phase-cycled images each with a different frequency sensitivity, such that the location of banding artifacts are shifted in space. These images are then combined to alleviate signal loss across the entire field-of-view. Although high levels of artifact suppression are viable using a large number of images, this is a time costly process that limits clinical utility. Here, we propose to accelerate individual acquisitions such that the overall scan time is equal to that of a single SSFP acquisition. Aliasing artifacts and noise are minimized by using a variable-density random sampling pattern in k-space, and by generating disjoint sampling patterns for separate acquisitions. A sparsity-enforcing method is then used for image reconstruction. Demonstrations on realistic brain phantom images, and in vivo brain and knee images are provided. In all cases, the proposed technique enables robust SSFP imaging in the presence of field inhomogeneities without prolonging scan times.

Description: Image guided interventions have seen growing interest in recent years. The use of X-rays for the procedure impels limiting the dose over time. Image sequences obtained thereby exhibit high levels of noise and very low contrasts. Hence, the development of efficient methods to enable optimal visualization of these sequences is crucial. We propose an original denoising method based on the curvelet transform. First, we apply a recursive temporal filter to the curvelet coefficients. As some residual noise remains, a spatial filtering is performed in the second step, which uses a magnitude-based classification and a contextual comparison of curvelet coefficients. This procedure allows to denoise the sequence while preserving low-contrasted structures, but does not improve their contrast. Finally, a third step is carried out to enhance the features of interest. For this, we propose a line enhancement technique in the curvelet domain. Indeed, thin structures are sparsely represented in that domain, allowing a fast and efficient detection. Quantitative and qualitative evaluations performed on synthetic and real low-dose sequences demonstrate that the proposed method enables a 50% dose reduction.

Description: Fluorescence lifetime imaging (FLIM) has demonstrated potential for robust assessment of atherosclerotic plaques biochemical composition and for complementing conventional intravascular ultrasound (IVUS), which provides information on plaque morphology. The success of such a bi-modal imaging modality depends on accurate segmentation of the IVUS images and proper angular registration between these two modalities. This paper reports a novel IVUS segmentation methodology addressing this issue. The image preprocessing consisted of denoising, using the Wiener filter, followed by image smoothing, implemented through the application of the alternating sequential filter on the edge separability metric images. Extraction of the lumen/intima and media/adventitia boundaries was achieved by tracing the gray-scale peaks over the A-lines of the IVUS preprocessed images. Cubic spline interpolation, in both cross-sectional and longitudinal directions, ensured boundary smoothness and continuity. The detection of the guide-wire artifact in both modalities is used for angular registration. Intraluminal studies were conducted in 13 ex vivo segments of human coronaries. The IVUS segmentation accuracy was assessed against independent manual tracings, providing 91.82% sensitivity and 97.55% specificity. The proposed methodology makes the bi-modal FLIM and IVUS approach feasible for comprehensive intravascular diagnosis by providing co-registered biochemical and morphological information of atherosclerotic plaques.

Description: Statistical X-ray computed tomography (CT) reconstruction can improve image quality from reduced dose scans, but requires very long computation time. Ordered subsets (OS) methods have been widely used for research in X-ray CT statistical image reconstruction (and are used in clinical PET and SPECT reconstruction). In particular, OS methods based on separable quadratic surrogates (OS-SQS) are massively parallelizable and are well suited to modern computing architectures, but the number of iterations required for convergence should be reduced for better practical use. This paper introduces OS-SQS-momentum algorithms that combine Nesterov's momentum techniques with OS-SQS methods, greatly improving convergence speed in early iterations. If the number of subsets is too large, the OS-SQS-momentum methods can be unstable, so we propose diminishing step sizes that stabilize the method while preserving the very fast convergence behavior. Experiments with simulated and real 3D CT scan data illustrate the performance of the proposed algorithms.

Description: A novel adaptation of the unscented Kalman filter (UKF) was recently introduced in literature for simultaneous multitensor estimation and fiber tractography from diffusion MRI. This technique has the advantage over other tractography methods in terms of computational efficiency, due to the fact that the UKF simultaneously estimates the diffusion tensors and propagates the most consistent direction to track along. This UKF and its variants reported later in literature however are not intrinsic to the space of diffusion tensors. Lack of this key property can possibly lead to inaccuracies in the multitensor estimation as well as in the tractography. In this paper, we propose a novel intrinsic unscented Kalman filter (IUKF) in the space of diffusion tensors which are symmetric positive definite matrices, that can be used for simultaneous recursive estimation of multitensors and propagation of directional information for use in fiber tractography from diffusion weighted MR data. In addition to being more accurate, IUKF retains all the advantages of UKF mentioned above. We demonstrate the accuracy and effectiveness of the proposed method via experiments publicly available phantom data from the fiber cup-challenge (MICCAI 2009) and diffusion weighted MR scans acquired from human brains and rat spinal cords.

Description: The detection of newly appearing and changing pigmented skin lesions (PSLs) is essential for timely diagnosis of cutaneous melanoma. Total body skin examination (TBSE) procedures, currently practiced for this purpose, can be extremely time-consuming for patients with numerous lesions. In addition, these procedures are prone to subjectivity when selecting PSLs for baseline image comparison, increasing the risk of missing a developing cancer. To address this issue, we propose a new photogrammetry-based total body scanning system allowing for skin surface image acquisition using cross-polarized light. Equipped with 21 high-resolution cameras and a turntable, this scanner automatically acquires a set of overlapping images, covering 85%–90% of the patient's skin surface. These images are used for the automated mapping of PSLs and their change estimation between explorations. The maps produced relate images of individual lesions with their locations on the patient's body, solving the body-to-image and image-to-image correspondence problem in TBSEs. Currently, the scanner is limited to patients with sparse body hair and, for a complete skin examination, the scalp, palms, soles and inner arms should be photographed manually. The initial tests of the scanner showed that it can be successfully applied for automated mapping and temporal monitoring of multiple lesions: PSLs relevant for follow-up were repeatedly mapped in several explorations. Moreover, during the baseline image comparison, all lesions with artificially induced changes were correctly identified as “evolved.”