MIXANDMIX: numerical techniques for the computation of empirical spectral distributions of population mixtures

The MIXANDMIX (mixtures by Anderson mixing) tool for the computation of the empirical spectral distribution of random matrices generated by mixtures of populations is described. Within the population mixture model the mapping between the population distributions and the limiting spectral distribution can be obtained by solving a set of systems of non-linear equations, for which an efficient implementation is provided. The contributions include a method for accelerated fixed point convergence, a homotopy continuation strategy to prevent convergence to non-admissible solutions, a blind non-uniform grid construction for effective distribution support detection and approximation, and a parallel computing architecture. Comparisons are performed with available packages for the single population case and with results obtained by simulation for the more general model implemented here. Results show competitive performance and improved flexibility.Comment: 17 pages, 6 figure

Non-Rigid Groupwise Registration for Motion Estimation and Compensation in Compressed Sensing Reconstruc- tion of Breath-Hold Cardiac Cine MRI

Purpose: Compressed sensing methods with motion estimation and compensation techniques have been proposed for the reconstruction of accelerated dynamic MRI. However, artifacts that naturally arise in compressed sensing reconstruction procedures hinder the estimation of motion from reconstructed images, especially at high acceleration factors. This work introduces a robust groupwise non-rigid motion estimation technique applied to the compressed sensing reconstruction of dynamic cardiac cine MRI sequences. Theory and Methods: A spatio-temporal regularized, groupwise, non-rigid registration method based on a B-splines deformation model and a least squares metric is used to estimate and to compensate the movement of the heart in breath-hold cine acquisitions and to obtain a quasi-static sequence with highly sparse representation in temporally transformed domains. Results: Short axis in vivo datasets are used for validation, both original multi-coil as well as DICOM data. Fully sampled data were retrospectively undersampled with various acceleration factors and reconstructions were compared with the two well-known methods k-t FOCUSS and MASTeR. The proposed method achieves higher signal to error ratio and structure similarity index for medium to high acceleration factors. Conclusions: Reconstruction methods based on groupwise registration show higher quality recon- structions for cardiac cine images than the pairwise counterparts tested

Autocalibrated cardiac tissue phase mapping with multiband imaging and k-t acceleration

PURPOSE: To develop an autocalibrated multiband (MB) CAIPIRINHA acquisition scheme with in-plane k-t acceleration enabling multislice three-directional tissue phase mapping in one breath-hold. METHODS: A k-t undersampling scheme was integrated into a time-resolved electrocardiographic-triggered autocalibrated MB gradient-echo sequence. The sequence was used to acquire data on 4 healthy volunteers with MB factors of two (MB2) and three (MB3), which were reconstructed using a joint reconstruction algorithm that tackles both k-t and MB acceleration. Forward simulations of the imaging process were used to tune the reconstruction model hyperparameters. Direct comparisons between MB and single-band tissue phase-mapping measurements were performed. RESULTS: Simulations showed that the velocities could be accurately reproduced with MB2 k-t (average ± twice the SD of the RMS error of 0.08 ± 0.22 cm/s and velocity peak reduction of 1.03% ± 6.47% compared with fully sampled velocities), whereas acceptable results were obtained with MB3 k-t (RMS error of 0.13 ± 0.58 cm/s and peak reduction of 2.21% ± 13.45%). When applied to tissue phase-mapping data, the proposed technique allowed three-directional velocity encoding to be simultaneously acquired at two/three slices in a single breath-hold of 18 heartbeats. No statistically significant differences were detected between MB2/MB3 k-t and single-band k-t motion traces averaged over the myocardium. Regional differences were found, however, when using the American Heart Association model for segmentation. CONCLUSION: An autocalibrated MB k-t acquisition / reconstruction framework is presented that allows three-directional velocity encoding of the myocardial velocities at multiple slices in one breath-hold

Multi-Stencil Streamline Fast Marching: a general 3D Framework to determine Myocardial Thickness and Transmurality in Late Enhancement Images

We propose a fully three-dimensional methodology for the computation of myocardial non-viable tissue transmurality in contrast enhanced magnetic resonance images. The outcome is a continuous map defined within the myocardium where not only current state-of-the-art measures of transmurality can be calculated, but also information on the location of non-viable tissue is preserved. The computation is done by means of a partial differential equation framework we have called Multi- Stencil Streamline Fast Marching (MSSFM). Using it, the myocardial and scarred tissue thickness is simultaneously computed. Experimental results show that the proposed 3D method allows for the computation of transmurality in myocardial regions where current 2D methods are not able to as conceived, and it also provides more robust and accurate results in situations where the assumptions on which current 2D methods are based —i.e., there is a visible endocardial contour and its corresponding epicardial points lie on the same slice—, are not met