In this study we classify short echo-time brain magnetic resonance spectroscopic imaging (MRSI) data by applying a model-based canonical correlation analyses algorithm and by using, as prior knowledge, multimodal sources of information coming from high-resolution magic angle spinning (HR-MAS), MRSI and magnetic resonance imaging. The potential and limitations of fusing in vivo and ex vivo nuclear magnetic resonance sources to detect brain tumors is investigated. We present various modalities for multimodal data fusion, study the effect and the impact of using multimodal information for classifying MRSI brain glial tumors data and analyze which parameters influence the classification results by means of extensive simulation and in vivo studies. Special attention is drawn to the possibility of considering HR-MAS data as a complementary dataset when dealing with a lack of MRSI data needed to build a classifier. Results show that HR-MAS information can have added value in the process of classifying MRSI data.

A new technique is presented to create nosologic images of the brain based on magnetic resonance imaging (MRI) and magnetic resonance spectroscopic imaging (MRSI). A nosologic image summarizes the presence of different tissues and lesions in a single image by color coding each voxel or pixel according to the histopathological class it is assigned to. The proposed technique applies advanced methods from image processing as well as pattern recognition to segment and classify brain tumours. First, a registered brain atlas and a subject-specific abnormal tissue prior, obtained from MRSI data, are used for the segmentation. Next, the detected abnormal tissue is classified based on supervised pattern recognition methods. Class probabilities are also calculated for the segmented abnormal region. Compared to previous approaches, the new framework is more flexible and able to better exploit spatial information leading to improved nosologic images. The combined scheme offers a new way to produce high-resolution nosologic images, representing tumour heterogeneity and class probabilities, which may help clinicians in decision making.

A new fast and accurate tissue typing technique has recently been successfully applied to prostate MR spectroscopic imaging (MRSI) data. This technique is based on canonical correlation analysis (CCA), a statistical method able to simultaneously exploit the spectral and spatial information characterizing the MRSI data. Here, the performance of CCA is further investigated by using brain data obtained by two-dimensional turbo spectroscopic imaging (2DTSI) from patients affected by glioblastoma. The purpose of this study is to investigate the applicability of CCA when typing tissues of heterogeneous tumors. The performance of CCA is also compared with that of ordinary correlation analysis on simulated as well as in vivo data. The results show that CCA outperforms ordinary correlation analysis in terms of robustness and accuracy.

A novel method is described that combines high-resolution scanning microdiffraction techniques, Rietveld quantitative phase analysis and a statistical method known as canonical correlation analysis (CCA). The method has been applied to a sample taken from a bone-tissue-engineered bioceramic porous scaffold implanted in a mouse for six months. The CCA technique allows the detection of those pixels throughout the investigated sample that best correlate with signal models. Besides the standard usage of this approach, which requires theoretical profiles as signal models, a novel application is presented here, which consists of picking the model spectra out of the experimental data set. Patterns representative of a reasonable range of phase compositions were selected among the huge number of two-dimensional patterns ( folded in onedimensional profiles) to extract quantitative phase fractions. At this stage, the CCA approach was also used to overcome the low Poisson statistic of signal models, so making Rietveld quantitative analysis more reliable. These patterns have been used as profile models for CCA. The final classification map, obtained by assigning the considered pixel to the model spectrum with the highest canonical coefficient, provides the spatial variation of phase concentration.

Precise knowledge of X-ray diffraction profile shape is crucial in the investigation of the properties of matter in crystals powder. Line-broadening analysis is the fourth pre-processing step in most of the full powder pattern fitting softwares. The final result of line-broadening analysis strongly depends on three further steps: noise filtering, removal of background signal, and peak fitting. In this work a new model independent procedure for two of the aforementioned steps (background suppression and peak fitting) is presented. The former is dealt with by using morphological mathematics, while the latter relies on the Hankel–Lanczos singular value decomposition technique. Real X-ray powder diffraction (XRPD) intensity profiles of Ceria samples are used to test the performance of the proposed procedure. Results show the robustness of this approach and its capability of efficiently improving the disentangling of instrumental broadening. These features make the proposed approach an interesting and user-friendly tool for the pre-processing of XRPD data.

A filter based on the Hankel-Lanczos singular value decomposition (HLSVD) technique is presented and applied for the first time to X-ray diffraction (XRD) data. Synthetic and real powder XRD intensity profiles of nanocrystals are used to study the filter performances with different noise levels. Results show the robustness of the HLSVD filter and its capability to extract easily and effciently the useful crystallographic information. These characteristics make the filter an interesting and user-friendly tool for processing of XRD data.

A reliable and automatic method is applied to crystallographic data for tissue typing. The technique is based on canonical correlation analysis, a statistical method which makes use of the spectral-spatial information characterizing X-ray diffraction data measured from bone samples with implanted tissues. The performance has been compared with a standard crystallographic technique in terms of accuracy and automation. The proposed approach is able to provide reliable tissue classification with a direct tissue visualization without requiring any user interaction.

Magnetic resonance spectroscopic imaging (MRSI) provides information about the spatial metabolic heterogeneity of an organ in the human body. In this way, MRSI can be used to detect tissue regions with abnormal metabolism, e.g. tumor tissue. The main drawback of MRSI in clinical practice is that the analysis of the data requires a lot of expertise from the radiologists. In this article, we present an automatic method that assigns each voxel of a spectroscopic image of the brain to a histopathological class. The method is based on Canonical Correlation Analysis (CCA), which has recently been shown to be a robust technique for tissue typing. In CCA, the spectral as well as the spatial information about the voxel is used to assign it to a class. This has advantages over other methods that only use spectral information since histopathological classes are normally covering neighbouring voxels. In this paper, a new CCA-based method is introduced in which MRSI and MR imaging information is integrated. The performance of tissue typing is compared for CCA applied to the whole MR spectra and to sets of features obtained from the spectra. Tests on simulated and in vivo MRSI data show that the new method is very accurate in terms of classification and segmentation. The results also show the advantage of combining spectroscopic and imaging data.

In this article an accurate and efficient technique for tissue typing is presented. The proposed technique is based on Canonical Correlation Analysis, a statistical method able to simultaneously exploit the spectral and spatial information characterizing the Magnetic Resonance Spectroscopic Imaging (MRSI) data. Recently, Canonical Correlation Analysis has been successfully applied to other types of biomedical data, such as functional MRI data. Here, Canonical Correlation Analysis is adapted for MRSI data processing in order to retrieve in an accurate and efficient way the possible tissue types that characterize the organ under investigation. The potential and limitations of the new technique have been investigated by using simulated as well as in vivo prostate MRSI data, and extensive studies demonstrate a high accuracy, robustness, and efficiency. Moreover, the performance of Canonical Correlation Analysis has been compared to that of ordinary correlation analysis. The test results show that Canonical Correlation Analysis performs best in terms of accuracy and robustness

Accurate quantitation of Magnetic Resonance Spectroscopy (MRS) signals is an essential step before converting the estimated signal parameters, such as frequencies, damping factors, and amplitudes, into biochemical quantities (concentration, pH). Several subspace-based parameter estimators have been developed for this task, which are efficient and accurate time-domain algorithms. However, they suffer from a serious drawback: they allow only a limited inclusion of prior knowledge which is important for accuracy and resolution. In this paper, a new method is presented: KNOB-SVD and its improved variant KNOB-TLS. KNOB-SVD is a recently proposed method, based on the Singular Value Decomposition (SVD), which allows the use of more prior knowledge about the signal parameters than previously published subspace-based methods. We compare its performance in terms of robustness and accuracy with the performance of three commonly used methods for signal parameter estimation: HTLS, a subspace-based method which does not allow any inclusion of prior knowledge, except for the model order; HTLSPK(Dfdeq), a subspace-based method obtained by incorporating in HTLS the prior information that the frequency differences between doublet components are known and the damping factors are equal; and AMARES, an interactive maximum likelihood method that allows the inclusion of a variety of prior knowledge. Extensive simulation and in vivo studies, using 31P as well as proton MRS signals, show that the new method outperforms HTLS and HTLSPK(Dfdeq) in robustness, accuracy, and resolution, and that it provides parameter estimates comparable to the AMARES ones.

Some inverse eigenvalue problems for matrices with Toeplitz-related structure are considered in this paper. In particular, the solutions of the inverse eigenvalue problems for Toeplitz-plus-Hankel matrices and for Toeplitz matrices having all double eigenvalues are characterized, respectively, in close form. Being centrosymmetric itself, the Toeplitz-plus-Hankel solution can be used as an initial value in a continuation method to solve the more difficult inverse eigenvalue problem for symmetric Toeplitz matrices. Numerical testing results show a clear advantage of such an application.

Magnetic resonance spectroscopy (MRS) has been shown to be a potentially important medical diagnostic tool. The success of MRS depends on the quantitative data analysis, i.e., the interpretation of the signal in terms of relevant physical parameters, such as frequencies, decay constants, and amplitudes. A variety of time-domain algorithms to extract parameters have been developed. On the one hand, there are so-called blackbox methods. Minimal user interaction and limited incorporation of prior knowledge are inherent to this type of method. On the other hand, interactive methods exist that are iterative, require user involvement, and allow inclusion of prior knowledge. We focus on blackbox methods. The computationally most intensive part of these blackbox methods is the computation of the singular value decomposition (SVD) of a Hankel matrix. Our goal is to reduce the needed computational time without affecting the accuracy of the parameters of interest. To this end, algorithms based on the Lanczos method are suitable because the main computation at each step, a matrix-vector product, can be efficiently performed by means of the fast Fourier transform exploiting the structure of the involved matrix. We compare the performance in terms of accuracy and efficiency of four algorithms: the classical SVD algorithm based on the QR decomposition, the Lanczos algorithm, the Lanczos algorithm with partial reorthogonalization, and the implicitly restarted Lanczos algorithm. Extensive simulation studies show that the latter two algorithms perform best. © 2002 Elsevier Science (USA).