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: The chaotic vascular network surrounding malignant tumors leads to pulsatile blood flow patterns that differ from those in benign regions of the breast. This study aimed to determine if high-speed electrical impedance tomography (EIT) is able to detect conductivity changes associated with cyclic blood-volume changes and to gauge the potential of using these signatures to differentiate malignant from benign regions within the breast. EIT imaging of pulsating latex membranes submerged in saline baths provided initial validation of its use for tracking temporally varying conductivities. Nineteen women (10 with cancer, nine without) were imaged with EIT over the course of several heartbeats in synchrony with pulse-oximetry acquisition. Eight parameters ( $r_{s}$ , $phi(r_{t,max})$ , $r_{t,max}$ , $P_{rm low:full}$ , $P_{rm high:full}$ , $P_{rm low:high}$ ) relating the conductivity images and pulse-oximeter signatures were extracted and used as a means of comparing malignant and benign regions of the breast. Significant differences $({ p}〈0.01)$ between malignant and benign regions of interest were noted in seven of the eight parameters. The maximum correlation between conductivity and pulse-oximeter signals, $r_{t,max}$ , was observed to be the optimal discriminating parameter with a receiver operating characteristic area under the curve of 0.8 and a specificity of 81% at a sensitivity of 77%. Assessing the dynamic conductivity of breast ma- provide additional clinical utility to that of standard imaging modalities, but further investigation is necessary to better understand the biophysical mechanisms leading to the observed conductivity changes.

Description: Classical surgery is being overtaken by minimally invasive and transcatheter procedures. As there is no direct view or access to the affected anatomy, advanced imaging techniques such as 3D C-arm computed tomography (CT) and C-arm fluoroscopy are routinely used in clinical practice for intraoperative guidance. However, due to constraints regarding acquisition time and device configuration, intraoperative modalities have limited soft tissue image quality and reliable assessment of the cardiac anatomy typically requires contrast agent, which is harmful to the patient and requires complex acquisition protocols. We propose a probabilistic sparse matching approach to fuse high-quality preoperative CT images and nongated, noncontrast intraoperative C-arm CT images by utilizing robust machine learning and numerical optimization techniques. Thus, high-quality patient-specific models can be extracted from the preoperative CT and mapped to the intraoperative imaging environment to guide minimally invasive procedures. Extensive quantitative experiments on 95 clinical datasets demonstrate that our model-based fusion approach has an average execution time of 1.56 s, while the accuracy of 5.48 mm between the anchor anatomy in both images lies within expert user confidence intervals. In direct comparison with image-to-image registration based on an open-source state-of-the-art medical imaging library and a recently proposed quasi-global, knowledge-driven multi-modal fusion approach for thoracic-abdominal images, our model-based method exhibits superior performance in terms of registration accuracy and robustness with respect to both target anatomy and anchor anatomy alignment errors.

Description: Image reconstruction from low-count positron emission tomography (PET) projection data is challenging because the inverse problem is ill-posed. Prior information can be used to improve image quality. Inspired by the kernel methods in machine learning, this paper proposes a kernel based method that models PET image intensity in each pixel as a function of a set of features obtained from prior information. The kernel-based image model is incorporated into the forward model of PET projection data and the coefficients can be readily estimated by the maximum likelihood (ML) or penalized likelihood image reconstruction. A kernelized expectation-maximization algorithm is presented to obtain the ML estimate. Computer simulations show that the proposed approach can achieve better bias versus variance trade-off and higher contrast recovery for dynamic PET image reconstruction than the conventional maximum likelihood method with and without post-reconstruction denoising. Compared with other regularization-based methods, the kernel method is easier to implement and provides better image quality for low-count data. Application of the proposed kernel method to a 4-D dynamic PET patient dataset showed promising results.

Description: Advances in magnetic resonance imaging (MRI) allow to gain critical insight into the structure of neural networks and their functional dynamics. To relate structural connectivity [as quantified by diffusion-weighted imaging (DWI) tractography] and functional connectivity [as obtained from functional MRI (fMRI)], increasing emphasis has been put on computational models of brain activity. In the present study, we use structural equation modeling (SEM) with structural connectivity to predict functional connectivity. The resulting model takes the simple form of a spatial simultaneous autoregressive model (sSAR), whose parameters can be estimated in a Bayesian framework. On synthetic data, results showed very good accuracy and reliability of the inference process. On real data, we found that the sSAR performed significantly better than two other reference models as well as than structural connectivity alone, but that the Bayesian procedure did not bring significant improvement in fit compared to two simpler approaches. Nonetheless, we also found that the values of the region-specific parameters inferred using Bayesian inference differed significantly across resting-state networks. These results demonstrate 1) that a simple abstract model is able to perform better that more complex models based on more realistic assumptions, 2) that the parameters of the sSAR can be estimated and can potentially be used as biomarkers, but also 3) that the sSAR, while being the best-performing model, is at best still a very crude model of the relationship between structure and function in MRI.

Description: We propose a novel deformation corrected compressed sensing (DC-CS) framework to recover contrast enhanced dynamic magnetic resonance images from undersampled measurements. We introduce a formulation that is capable of handling a wide class of sparsity/compactness priors on the deformation corrected dynamic signal. In this work, we consider example compactness priors such as sparsity in temporal Fourier domain, sparsity in temporal finite difference domain, and nuclear norm penalty to exploit low rank structure. Using variable splitting, we decouple the complex optimization problem to simpler and well understood sub problems; the resulting algorithm alternates between simple steps of shrinkage-based denoising, deformable registration, and a quadratic optimization step. Additionally, we employ efficient continuation strategies to reduce the risk of convergence to local minima. The decoupling enabled by the proposed scheme enables us to apply this scheme to contrast enhanced MRI applications. Through experiments on numerical phantom and in vivo myocardial perfusion MRI datasets, we observe superior image quality of the proposed DC-CS scheme in comparison to the classical k-t FOCUSS with motion estimation/correction scheme, and demonstrate reduced motion artifacts over classical compressed sensing schemes that utilize the compact priors on the original deformation uncorrected signal.

Description: Identification of error in nonrigid registration is a critical problem in the medical image processing community. We recently proposed an algorithm that we call “Assessing Quality Using Image Registration Circuits” (AQUIRC) to identify nonrigid registration errors and have tested its performance using simulated cases. In this paper, we extend our previous work to assess AQUIRC's ability to detect local nonrigid registration errors and validate it quantitatively at specific clinical landmarks, namely the anterior commissure and the posterior commissure. To test our approach on a representative range of error we utilize five different registration methods and use 100 target images and nine atlas images. Our results show that AQUIRC's measure of registration quality correlates with the true target registration error (TRE) at these selected landmarks with an ${rm R}^{2}=0.542$ . To compare our method to a more conventional approach, we compute local normalized correlation coefficient (LNCC) and show that AQUIRC performs similarly. However, a multi-linear regression performed with both AQUIRC's measure and LNCC shows a higher correlation with TRE than correlations obtained with either measure alone, thus showing the complementarity of these quality measures. We conclude the paper by showing that the AQUIRC algorithm can be used to reduce registration errors for all five algorithms.

Description: This paper presents a robust atlas-based segmentation (ABS) algorithm for segmentation of the breast boundary in 3-D MR images. The proposed algorithm combines the well-known methodologies of ABS namely probabilistic atlas and atlas selection approaches into a single framework where two configurations are realized. The algorithm uses phase congruency maps to create an atlas which is robust to intensity variations. This allows an atlas derived from images acquired with one MR imaging sequence to be used to segment images acquired with a different MR imaging sequence and eliminates the need for intensity-based registration. Images acquired using a Dixon sequence were used to create an atlas which was used to segment both Dixon images (intra-sequence) and T1-weighted images (inter-sequence). In both cases, highly accurate results were achieved with the median Dice similarity coefficient values of $94% pm 4%$ and $87 pm 6.5%$ , respectively.