The magnetic susceptibility of tissue can be determined in gradient echo MRI by deconvolving the local magnetic field with the magnetic field generated by a unit dipole. This Quantitative Susceptibility Mapping (QSM) problem is unfortunately ill-posed. By transforming the problem to the Fourier domain, the susceptibility appears to be undersampled only at points where the dipole kernel is zero, suggesting that a modest amount of additional information may be sufficient for uniquely resolving susceptibility. A Morphology Enabled Dipole Inversion (MEDI) approach is developed that exploits the structural consistency between the susceptibility map and the magnitude image reconstructed from the same gradient echo MRI. Specifically, voxels that are part of edges in the susceptibility map but not in the edges of the magnitude image are considered to be sparse. In this approach an L1 norm minimization is used to express this sparsity property. Numerical simulations and phantom experiments are performed to demonstrate the superiority of this L1 minimization approach over the previous L2 minimization method. Preliminary brain imaging results in healthy subjects and in patients with intracerebral hemorrhages illustrate that QSM is feasible in practice.

For optimal image quality in susceptibility-weighted imaging and accurate quantification of susceptibility, it is necessary to isolate the local field generated by local magnetic sources (such as iron) from the background field that arises from imperfect shimming and variations in magnetic susceptibility of surrounding tissues (including air). Previous background removal techniques have limited effectiveness depending on the accuracy of model assumptions or information input. In this article, we report an observation that the magnetic field for a dipole outside a given region of interest (ROI) is approximately orthogonal to the magnetic field of a dipole inside the ROI. Accordingly, we propose a nonparametric background field removal technique based on projection onto dipole fields (PDF). In this PDF technique, the background field inside an ROI is decomposed into a field originating from dipoles outside the ROI using the projection theorem in Hilbert space. This novel PDF background removal technique was validated on a numerical simulation and a phantom experiment and was applied in human brain imaging, demonstrating substantial improvement in background field removal compared with the commonly used high-pass filtering method.

Magnetic susceptibility varies among brain structures and provides insights into the chemical and molecular composition of brain tissues. However, the determination of an arbitrary susceptibility distribution from the measured MR signal phase is a challenging, ill-conditioned inverse problem. Although a previous method named calculation of susceptibility through multiple orientation sampling (COSMOS) has solved this inverse problem both theoretically and experimentally using multiple angle acquisitions, it is often impractical to carry out on human subjects. Recently, the feasibility of calculating the brain susceptibility distribution from a singleangle acquisition was demonstrated using morphology enabled dipole inversion (MEDI). In this study, we further improved the original MEDI method by sparsifying the edges in the quantitative susceptibility map that do not have a corresponding edge in the magnitude image. Quantitative susceptibility maps generated by the improved MEDI were compared qualitatively and quantitatively with those generated by calculation of susceptibility through multiple orientation sampling. The results show a high degree of agreement between MEDI and calculation of susceptibility through multiple orientation sampling, and the practicality of MEDI allows many potential clinical applications. Magn Reson Med 66:777-783,

a b s t r a c tWe propose a new framework to extract the activity-related component in the BOLD functional magnetic resonance imaging (fMRI) signal. As opposed to traditional fMRI signal analysis techniques, we do not impose any prior knowledge of the event timing. Instead, our basic assumption is that the activation pattern is a sequence of short and sparsely distributed stimuli, as is the case in slow event-related fMRI.We introduce new wavelet bases, termed ''activelets'', which sparsify the activityrelated BOLD signal. These wavelets mimic the behavior of the differential operator underlying the hemodynamic system. To recover the sparse representation, we deploy a sparse-solution search algorithm.The feasibility of the method is evaluated using both synthetic and experimental fMRI data. The importance of the activelet basis and the non-linear sparse recovery algorithm is demonstrated by comparison against classical B-spline wavelets and linear regularization, respectively.

Abstract-In this paper, an approach is introduced based on differential operators to construct wavelet-like basis functions. Given a differential operator L with rational transfer function, elementary building blocks are obtained that are shifted replicates of the Green's function of L. It is shown that these can be used to specify a sequence of embedded spline spaces that admit a hierarchical exponential B-spline representation. The corresponding B-splines are entirely specified by their poles and zeros; they are compactly supported, have an explicit analytical form, and generate multiresolution Riesz bases. Moreover, they satisfy generalized refinement equations with a scale-dependent filter and lead to a representation that is dense in 2 . This allows us to specify a corresponding family of semi-orthogonal exponential spline wavelets, which provides a major extension of earlier polynomial spline constructions. These wavelets are completely characterized, and it is proven that they satisfy the following remarkable properties: 1) they are orthogonal across scales and generate Riesz bases at each resolution level; 2) they yield unconditional bases of 2 -either compactly supported (B-spline-type) or with exponential decay (orthogonal or dual-type); 3) they have vanishing exponential moments, where is the order of the differential operator; 4) they behave like multiresolution versions of the operator L from which they are derived; and 5) their order of approximation is ( ), where and give the number of poles and zeros, respectively. Last but not least, the new wavelet-like decompositions are as computationally efficient as the classical ones. They are computed using an adapted version of Mallat's filter bank algorithm, where the filters depend on the decomposition level.

Dysregulated host inflammatory response causes many diseases, including cardiovascular and neurodegenerative diseases, cancer, and sepsis. Sensitive detection of the site of inflammation will, therefore, produce a wide-ranging impact on disease diagnosis and treatment. We hypothesized that nanoprobes designed to mimic the molecular interactions occurring between inflamed leukocytes and endothelium may possess selectivity toward diverse host inflammatory responses. To incorporate inflammation-sensitive molecular interactions, superparamagnetic iron oxide nanoparticles were conjugated with integrin lymphocyte function-associated antigen (LFA)-1 I domain, engineered to mimic activated leukocytes in physiology. Whole body optical and magnetic resonance imaging in vivo revealed that leukocyte-mimetic nanoparticles localized preferentially to the vasculature within and in the invasive front of the tumor, as well as to the site of acute inflammation. This study explored in vivo detection of tumor-associated vasculature with systemically injected inflammation-specific nanoparticles, presenting a possibility of tumor detection by inflamed tumor microenvironment.

Abstract-Tomographic reconstruction from positron emission tomography (PET) data is an ill-posed problem that requires regularization. An attractive approach is to impose an 1 -regularization constraint, which favors sparse solutions in the wavelet domain. This can be achieved quite efficiently thanks to the iterative algorithm developed by Daubechies et al., 2004. In this paper, we apply this technique and extend it for the reconstruction of dynamic (spatio-temporal) PET data. Moreover, instead of using classical wavelets in the temporal dimension, we introduce exponential-spline wavelets (E-spline wavelets) that are specially tailored to model time activity curves (TACs) in PET. We show that the exponential-spline wavelets naturally arise from the compartmental description of the dynamics of the tracer distribution. We address the issue of the selection of the "optimal" E-spline parameters (poles and zeros) and we investigate their effect on reconstruction quality. We demonstrate the usefulness of spatio-temporal regularization and the superior performance of E-spline wavelets over conventional Battle-Lemarié wavelets in a series of experiments: the 1-D fitting of TACs, and the tomographic reconstruction of both simulated and clinical data. We find that the E-spline wavelets outperform the conventional wavelets in terms of the reconstructed signal-to-noise ratio (SNR) and the sparsity of the wavelet coefficients. Based on our simulations, we conclude that replacing the conventional wavelets with E-spline wavelets leads to equal reconstruction quality for a 40% reduction of detected coincidences, meaning an improved image quality for the same number of counts or equivalently a reduced exposure to the patient for the same image quality.

No abstract