Effects of apodization smoothing and denoising on spectral fitting

A Coupled Statistical and Deterministic Model for Forecasting Climate-Driven Dengue Incidence in Selangor, Malaysia

The mosquito-borne dengue virus remains a major public health concern in Malaysia. Despite various control efforts and measures introduced by the Malaysian Government to combat dengue, the increasing trend of dengue cases persists and shows no sign of decreasing. Currently, early detection and vector control are the main methods employed to curb dengue outbreaks. In this study, a coupled model consisting of the statistical ARIMAX model and the deterministic SI-SIR model was developed and validated using the weekly reported dengue data from year 2014 to 2019 for Selangor, Malaysia. Previous studies have shown that climate variables, especially temperature, humidity, and precipitation, were able to influence dengue incidence and transmission dynamics through their effect on the vector. In this coupled model, climate is linked to dengue disease through mosquito biting rate, allowing real-time forecast of dengue cases using climate variables, namely temperature, rainfall and humidity. For the period chosen for model validation, the coupled model can forecast 1-2 weeks in advance with an average error of less than 6%, three weeks in advance with an average error of 7.06% and four weeks in advance with an average error of 8.01%. Further model simulation analysis suggests that the coupled model generally provides better forecast than the stand-alone ARIMAX model, especially at the onset of the outbreak. Moreover, the coupled model is more robust in the sense that it can be further adapted for investigating the effectiveness of various dengue mitigation measures subject to the changing climate.

Denoising single MR spectra by deep learning: Miracle or mirage?

PURPOSE

The inherently poor SNR of MRS measurements presents a significant hurdle to its clinical application. Denoising by machine or deep learning (DL) was proposed as a remedy. It is investigated whether such denoising leads to lower estimate uncertainties or whether it essentially reduces noise in signal-free areas only.

METHODS

Noise removal based on supervised DL with U-nets was implemented using simulated 1 H MR spectra of human brain in two approaches: (1) via time-frequency domain spectrograms and (2) using 1D spectra as input. Quality of denoising was evaluated in three ways: (1) by an adapted fit quality score, (2) by traditional model fitting, and (3) by quantification via neural networks.

RESULTS

Visually appealing spectra were obtained; hinting that denoising is well-suited for MRS. However, an adapted denoising score showed that noise removal is inhomogeneous and more efficient for signal-free areas. This was confirmed by quantitative analysis of traditional fit results as well as DL quantitation following DL denoising. DL denoising, although apparently successful as judged by mean squared errors, led to substantially biased estimates in both implementations.

CONCLUSION

The implemented DL-based denoising techniques may be useful for display purposes, but do not help quantitative evaluations, confirming expectations based on estimation theory: Cramér Rao lower bounds defined by the original data and the appropriate fitting model cannot be circumvented in an unbiased way for single data sets, unless additional prior knowledge can be incurred in the form of parameter restrictions/relations or applicable substates.

SWADESH: a multimodal multi-disease brain imaging and neuropsychological database and data analytics platform

Multimodal neuroimaging data of various brain disorders provides valuable information to understand brain function in health and disease. Various neuroimaging-based databases have been developed that mainly consist of volumetric magnetic resonance imaging (MRI) data. We present the comprehensive web-based neuroimaging platform "SWADESH" for hosting multi-disease, multimodal neuroimaging, and neuropsychological data along with analytical pipelines. This novel initiative includes neurochemical and magnetic susceptibility data for healthy and diseased conditions, acquired using MR spectroscopy (MRS) and quantitative susceptibility mapping (QSM) respectively. The SWADESH architecture also provides a neuroimaging database which includes MRI, MRS, functional MRI (fMRI), diffusion weighted imaging (DWI), QSM, neuropsychological data and associated data analysis pipelines. Our final objective is to provide a master database of major brain disease states (neurodegenerative, neuropsychiatric, neurodevelopmental, and others) and to identify characteristic features and biomarkers associated with such disorders.

Signal-to-Noise Ratio Enhancement of Single-Voxel In Vivo 31P and 1H Magnetic Resonance Spectroscopy in Mice Brain Data Using Low-Rank Denoising

Magnetic resonance spectroscopy (MRS) is a noninvasive technique for measuring metabolite concentration. It can be used for preclinical small animal brain studies using rodents to provide information about neurodegenerative diseases and metabolic disorders. However, data acquisition from small volumes in a limited scan time is technically challenging due to its inherently low sensitivity. To mitigate this problem, this study investigated the feasibility of a low-rank denoising method in enhancing the quality of single voxel multinuclei (³¹P and ¹H) MRS data at 9.4 T. Performance was evaluated using in vivo MRS data from a normal mouse brain (³¹P and ¹H) and stroke mouse model (¹H) by comparison with signal-to-noise ratios (SNRs), Cramer-Rao lower bounds (CRLBs), and metabolite concentrations of a linear combination of model analysis results. In ³¹P MRS data, low-rank denoising resulted in improved SNRs and reduced metabolite quantification uncertainty compared with the original data. In ¹H MRS data, the method also improved the SNRs, CRLBs, but it performed better for ³¹P MRS data with relatively simpler patterns compared to the ¹H MRS data. Therefore, we suggest that the low-rank denoising method can improve spectra SNR and metabolite quantification uncertainty in single-voxel in vivo ³¹P and ¹H MRS data, and it might be more effective for ³¹P MRS data. The main contribution of this study is that we demonstrated the effectiveness of the low-rank denoising method on small-volume single-voxel MRS data. We anticipate that our results will be useful for the precise quantification of low-concentration metabolites, further reducing data acquisition voxel size, and scan time in preclinical MRS studies.

MP-PCA denoising for diffusion MRS data: promises and pitfalls

Diffusion-weighted (DW) magnetic resonance spectroscopy (MRS) suffers from a lower signal to noise ratio (SNR) compared to conventional MRS owing to the addition of diffusion attenuation. This technique can therefore strongly benefit from noise reduction strategies. In the present work, Marchenko-Pastur principal component analysis (MP-PCA) denoising is tested on Monte Carlo simulations and on in vivo DW-MRS data acquired at 9.4 T in rat brain and at 3 T in human brain. We provide a descriptive study of the effects observed following different MP-PCA denoising strategies (denoising the entire matrix versus using a sliding window), in terms of apparent SNR, rank selection, noise correlation within and across b-values and quantification of metabolite concentrations and fitted diffusion coefficients. MP-PCA denoising yielded an increased apparent SNR, a more accurate B0 drift correction between shots, and similar estimates of metabolite concentrations and diffusivities compared to the raw data. No spectral residuals on individual shots were observed but correlations in the noise level across shells were introduced, an effect which was mitigated using a sliding window, but which should be carefully considered.

PCA denoising and Wiener deconvolution of 31 P 3D CSI data to enhance effective SNR and improve point spread function

Purpose

This study evaluates the performance of 2 processing methods, that is, principal component analysis‐based denoising and Wiener deconvolution, to enhance the quality of phosphorus 3D chemical shift imaging data.

Methods

Principal component analysis‐based denoising increases the SNR while maintaining spectral information. Wiener deconvolution reduces the FWHM of the voxel point spread function, which is increased by Hamming filtering or Hamming‐weighted acquisition. The proposed methods are evaluated using simulated and in vivo 3D phosphorus chemical shift imaging data by 1) visual inspection of the spatial signal distribution; 2) SNR calculation of the PCr peak; and 3) fitting of metabolite basis functions.

Results

With the optimal order of processing steps, we show that the effective SNR of in vivo phosphorus 3D chemical shift imaging data can be increased. In simulations, we show we can preserve phosphorus‐containing metabolite peaks that had an SNR < 1 before denoising. Furthermore, using Wiener deconvolution, we were able to reduce the FWHM of the voxel point spread function with only partially reintroducing Gibb‐ringing artifacts while maintaining the SNR. After data processing, fitting of the phosphorus‐containing metabolite signals improved.

Conclusion

In this study, we have shown that principal component analysis‐based denoising in combination with regularized Wiener deconvolution allows increasing the effective spectral SNR of in vivo phosphorus 3D chemical shift imaging data, with reduction of the FWHM of the voxel point spread function. Processing increased the effective SNR by at least threefold compared to Hamming weighted acquired data and minimized voxel bleeding. With these methods, fitting of metabolite amplitudes became more robust with decreased fitting residuals.