Statistical characterization of residual noise in the low-rank approximation filter framework, general theory and application to hyperpolarized tracer spectroscopy

The use of low-rank approximation filters in the field of NMR is increasing due to their flexibility and effectiveness. Despite their ability to reduce the Mean Square Error between the processed signal and the true signal is well known, the statistical distribution of the residual noise is still undescribed. In this article, we show that low-rank approximation filters are equivalent to linear filters, and we calculate the mean and the covariance matrix of the processed data. We also show how to use this knowledge to build a maximum likelihood estimator, and we test the estimator's performance with a Montecarlo simulation of a 13C pyruvate metabolic tracer. While the article focuses on NMR spectroscopy experiment with hyperpolarized tracer, we also show that the results can be applied to tensorial data (e.g. using HOSVD) or 1D data (e.g. Cadzow filter).Comment: 26 pages, 7 figure

Dual modality optical coherence tomography : Technology development and biomedical applications

Optical coherence tomography (OCT) is a cross-sectional imaging modality that is widely used in clinical ophthalmology and interventional cardiology. It is highly promising for in situ characterization of tumor tissues. OCT has high spatial resolution and high imaging speed to assist clinical decision making in real-time. OCT can be used in both structural imaging and mechanical characterization. Malignant tumor tissue alters morphology. Additionally, structural OCT imaging has limited tissue differentiation capability because of the complex and noisy nature of the OCT signal. Moreover, the contrast of structural OCT signal derived from tissue’s light scattering properties has little chemical specificity. Hence, interrogating additional tissue properties using OCT would improve the outcome of OCT’s clinical applications. In addition to morphological difference, pathological tissue such as cancer breast tissue usually possesses higher stiffness compared to the normal healthy tissue, which indicates a compelling reason for the specific combination of structural OCT imaging with stiffness assessment in the development of dual-modality OCT system for the characterization of the breast cancer diagnosis. This dissertation seeks to integrate the structural OCT imaging and the optical coherence elastography (OCE) for breast cancer tissue characterization. OCE is a functional extension of OCT. OCE measures the mechanical response (deformation, resonant frequency, elastic wave propagation) of biological tissues under external or internal mechanical stimulation and extracts the mechanical properties of tissue related to its pathological and physiological processes. Conventional OCE techniques (i.e., compression, surface acoustic wave, magnetomotive OCE) measure the strain field and the results of OCE measurement are different under different loading conditions. Inconsistency is observed between OCE characterization results from different measurement sessions. Therefore, a robust mechanical characterization is required for force/stress quantification. A quantitative optical coherence elastography (qOCE) that tracks both force and displacement is proposed and developed at NJIT. qOCE instrument is based on a fiber optic probe integrated with a Fabry-Perot force sensor and the miniature probe can be delivered to arbitrary locations within animal or human body. In this dissertation, the principle of qOCE technology is described. Experimental results are acquired to demonstrate the capability of qOCE in characterizing the elasticity of biological tissue. Moreover, a handheld optical instrument is developed to allow in vivo real-time OCE characterization based on an adaptive Doppler analysis algorithm to accurately track the motion of sample under compression. For the development of the dual modality OCT system, the structural OCT images exhibit additive and multiplicative noises that degrade the image quality. To suppress noise in OCT imaging, a noise adaptive wavelet thresholding (NAWT) algorithm is developed to remove the speckle noise in OCT images. NAWT algorithm characterizes the speckle noise in the wavelet domain adaptively and removes the speckle noise while preserving the sample structure. Furthermore, a novel denoising algorithm is also developed that adaptively eliminates the additive noise from the complex OCT using Doppler variation analysis

Integration of magnetic resonance spectroscopic imaging into the radiotherapy treatment planning

L'objectif de cette thèse est de proposer de nouveaux algorithmes pour surmonter les limitations actuelles et de relever les défis ouverts dans le traitement de l'imagerie spectroscopique par résonance magnétique (ISRM). L'ISRM est une modalité non invasive capable de fournir la distribution spatiale des composés biochimiques (métabolites) utilisés comme biomarqueurs de la maladie. Les informations fournies par l'ISRM peuvent être utilisées pour le diagnostic, le traitement et le suivi de plusieurs maladies telles que le cancer ou des troubles neurologiques. Cette modalité se montre utile en routine clinique notamment lorsqu'il est possible d'en extraire des informations précises et fiables. Malgré les nombreuses publications sur le sujet, l'interprétation des données d'ISRM est toujours un problème difficile en raison de différents facteurs tels que le faible rapport signal sur bruit des signaux, le chevauchement des raies spectrales ou la présence de signaux de nuisance. Cette thèse aborde le problème de l'interprétation des données d'ISRM et la caractérisation de la rechute des patients souffrant de tumeurs cérébrales. Ces objectifs sont abordés à travers une approche méthodologique intégrant des connaissances a priori sur les données d'ISRM avec une régularisation spatio-spectrale. Concernant le cadre applicatif, cette thèse contribue à l'intégration de l'ISRM dans le workflow de traitement en radiothérapie dans le cadre du projet européen SUMMER (Software for the Use of Multi-Modality images in External Radiotherapy) financé par la Commission européenne (FP7-PEOPLE-ITN).The aim of this thesis is to propose new algorithms to overcome the current limitations and to address the open challenges in the processing of magnetic resonance spectroscopic imaging (MRSI) data. MRSI is a non-invasive modality able to provide the spatial distribution of relevant biochemical compounds (metabolites) commonly used as biomarkers of disease. Information provided by MRSI can be used as a valuable insight for the diagnosis, treatment and follow-up of several diseases such as cancer or neurological disorders. Obtaining accurate and reliable information from in vivo MRSI signals is a crucial requirement for the clinical utility of this technique. Despite the numerous publications on the topic, the interpretation of MRSI data is still a challenging problem due to different factors such as the low signal-to-noise ratio (SNR) of the signals, the overlap of spectral lines or the presence of nuisance components. This thesis addresses the problem of interpreting MRSI data and characterizing recurrence in tumor brain patients. These objectives are addressed through a methodological approach based on novel processing methods that incorporate prior knowledge on the MRSI data using a spatio-spectral regularization. As an application, the thesis addresses the integration of MRSI into the radiotherapy treatment workflow within the context of the European project SUMMER (Software for the Use of Multi-Modality images in External Radiotherapy) founded by the European Commission (FP7-PEOPLE-ITN framework)

MR-based protein imaging of the human brain by means of dualCEST

Chemical exchange saturation transfer (CEST) is an emerging magnetic resonance imaging (MRI) technique enabling indirect detection of low-concentration cellular compounds in living tissue by their magnetization transfer with water. In particular, protein-attributed CEST signals have been shown to provide valuable diagnostic information for various diseases. While conventional CEST approaches suffer from confounding signals from metabolites and macromolecules, the novel dual-frequency irradiation CEST (dualCEST) technique enables increased protein specificity by selectively detecting the intramolecular spin-diffusion. However, application of this technique has so far been limited to spectroscopic investigations of model solutions at ultrahigh magnetic field strengths. In this thesis, dualCEST was translated to a clinical whole-body MR scanner, enabling protein imaging of the human brain. To this end, several methodological developments were implemented and optimized: (i) improved dual-frequency pulses for signal preparation, (ii) a fast and robust volumetric image readout, (iii) a weighted acquisition scheme, and (iv) an adaptive denoising technique. The resulting improvements are not limited to dualCEST but are relevant for the research field of CEST-MRI in general. Extensive measurements of biochemical model solutions and volunteers demonstrated the protein specificity and reproducibility of dualCEST-MRI. The clinical applicability was verified in pilot studies with tumor and Alzheimer’s patients

Superresolution Reconstruction for Magnetic Resonance Spectroscopic Imaging Exploiting Low-Rank Spatio-Spectral Structure

Magnetic resonance spectroscopic imaging (MRSI) is a rapidly developing medical imaging modality, capable of conferring both spatial and spectral information content, and has become a powerful clinical tool. The ability to non-invasively observe spatial maps of metabolite concentrations, for instance, in the human brain, can offer functional, as well as pathological insights, perhaps even before structural aberrations or behavioral symptoms are evinced. Despite its lofty clinical prospects, MRSI has traditionally remained encumbered by a number of practical limitations. Of primary concern are the vastly reduced concentrations of tissue metabolites when compared to that of water, which forms the basis for conventional MR imaging. Moreover, the protracted exam durations required by MRSI routinely approach the limits for patient compliance. Taken in conjunction, the above considerations effectively circumscribe the data collection process, ultimately translating to coarse image resolutions that are of diminished clinical utility. Such shortcomings are compounded by spectral contamination artifacts due to the system pointspread function, which arise as a natural consequence when reconstructing non-band-limited data by the inverse Fourier transform. These artifacts are especially pronounced near regions characterized by substantial discrepancies in signal intensity, for example, the interface between normal brain and adipose tissue, whereby the metabolite signals are inundated by the dominant lipid resonances. In recent years, concerted efforts have been made to develop alternative, non-Fourier MRSI reconstruction strategies that aim to surmount the aforementioned limitations. In this dissertation, we build upon the burgeoning medley of innovative and promising techniques, proffering a novel superresolution reconstruction framework predicated on the recent interest in low-rank signal modeling, along with state-of-the-art regularization methods. The proposed framework is founded upon a number of key tenets. Firstly, we proclaim that the underlying spatio-spectral distribution of the investigated object admits a bilinear representation, whereby spatial and spectral signal components can be effectively segregated. We further maintain that the dimensionality of the subspace spanned by the components is, in principle, bounded by a modest number of observable metabolites. Secondly, we assume that local susceptibility effects represent the primary sources of signal corruption that tend to disallow such representations. Finally, we assert that the spatial components belong to a class of real-valued, non-negative, and piecewise linear functions, compelled in part through the use of a total variation regularization penalty. After demonstrating superior spatial and spectral localization properties in both numerical and physical phantom data when compared against standard Fourier methods, we proceed to evaluate reconstruction performance in typical in vivo settings, whereby the method is extended in order to promote the recovery of signal variations throughout the MRSI slice thickness. Aside from the various technical obstacles, one of the cardinal prospective challenges for high-resolution MRSI reconstruction is the shortfall of reliable ground truth data prudent for validation, thereby prompting reservations surrounding the resulting experimental outcomes. [...

Advanced Image Acquisition, Processing Techniques and Applications

"Advanced Image Acquisition, Processing Techniques and Applications" is the first book of a series that provides image processing principles and practical software implementation on a broad range of applications. The book integrates material from leading researchers on Applied Digital Image Acquisition and Processing. An important feature of the book is its emphasis on software tools and scientific computing in order to enhance results and arrive at problem solution

A Multiscale Approach for Statistical Characterization of Functional Images

Increasingly, scientific studies yield functional image data, in which the observed data consist of sets of curves recorded on the pixels of the image. Examples include temporal brain response intensities measured by fMRI and NMR frequency spectra measured at each pixel. This article presents a new methodology for improving the characterization of pixels in functional imaging, formulated as a spatial curve clustering problem. Our method operates on curves as a unit. It is nonparametric and involves multiple stages: (i) wavelet thresholding, aggregation, and Neyman truncation to effectively reduce dimensionality; (ii) clustering based on an extended EM algorithm; and (iii) multiscale penalized dyadic partitioning to create a spatial segmentation. We motivate the different stages with theoretical considerations and arguments, and illustrate the overall procedure on simulated and real datasets. Our method appears to offer substantial improvements over monoscale pixel-wise methods. An Appendix which gives some theoretical justifications of the methodology, computer code, documentation and dataset are available in the online supplements

Compressed Sensing Accelerated Magnetic Resonance Spectroscopic Imaging

abstract: Magnetic resonance spectroscopic imaging (MRSI) is a valuable technique for assessing the in vivo spatial profiles of metabolites like N-acetylaspartate (NAA), creatine, choline, and lactate. Changes in metabolite concentrations can help identify tissue heterogeneity, providing prognostic and diagnostic information to the clinician. The increased uptake of glucose by solid tumors as compared to normal tissues and its conversion to lactate can be exploited for tumor diagnostics, anti-cancer therapy, and in the detection of metastasis. Lactate levels in cancer cells are suggestive of altered metabolism, tumor recurrence, and poor outcome. A dedicated technique like MRSI could contribute to an improved assessment of metabolic abnormalities in the clinical setting, and introduce the possibility of employing non-invasive lactate imaging as a powerful prognostic marker. However, the long acquisition time in MRSI is a deterrent to its inclusion in clinical protocols due to associated costs, patient discomfort (especially in pediatric patients under anesthesia), and higher susceptibility to motion artifacts. Acceleration strategies like compressed sensing (CS) permit faithful reconstructions even when the k-space is undersampled well below the Nyquist limit. CS is apt for MRSI as spectroscopic data are inherently sparse in multiple dimensions of space and frequency in an appropriate transform domain, for e.g. the wavelet domain. The objective of this research was three-fold: firstly on the preclinical front, to prospectively speed-up spectrally-edited MRSI using CS for rapid mapping of lactate and capture associated changes in response to therapy. Secondly, to retrospectively evaluate CS-MRSI in pediatric patients scanned for various brain-related concerns. Thirdly, to implement prospective CS-MRSI acquisitions on a clinical magnetic resonance imaging (MRI) scanner for fast spectroscopic imaging studies. Both phantom and in vivo results demonstrated a reduction in the scan time by up to 80%, with the accelerated CS-MRSI reconstructions maintaining high spectral fidelity and statistically insignificant errors as compared to the fully sampled reference dataset. Optimization of CS parameters involved identifying an optimal sampling mask for CS-MRSI at each acceleration factor. It is envisioned that time-efficient MRSI realized with optimized CS acceleration would facilitate the clinical acceptance of routine MRSI exams for a quantitative mapping of important biomarkers.Dissertation/ThesisDoctoral Dissertation Bioengineering 201