MRI denoising using nonlocal neutrosophic set approach of Wiener filtering

In this paper, a new filtering method is presented to remove the Rician noise from magnetic resonance images (MRI) acquired using single coil MRI acquisition system. This filter is based on nonlocal neutrosophic set (NLNS) approach of Wiener filtering. A neutrosophic set (NS), a part of neutrosophy theory, studies the origin, nature, and scope of neutralities, as well as their interactions with different ideational spectra. Now, we apply the neutrosophic set into image domain and define some concepts and operators for image denoising. First, the nonlocal mean is applied to the noisy MRI. The resultant image is transformed into NS domain, described using three membership sets: true (T), indeterminacy (I) and false (F). The entropy of the neutrosophic set is defined and employed to measure the indeterminacy. The ω-Wiener filtering operation is used on T and F to decrease the set indeterminacy and to remove the noise. The experiments have been conducted on simulated MR images from Brainweb database and clinical MR images. The results show that the NLNS Wiener filter produces better denoising results in terms of qualitative and quantitative measures compared with other denoising methods, such as classical Wiener filter, the anisotropic diffusion filter, the total variation minimization and the nonlocal means filter. The visual and the diagnostic quality of the denoised image are well preserved.     

A comparison of denoising methods for X-ray fluoroscopic images

Fluoroscopic images exhibit severe signal-dependent quantum noise, due to the reduced X-ray dose involved in image formation, that is generally modelled as Poisson-distributed. However, image gray-level transformations, commonly applied by fluoroscopic device to enhance contrast, modify the noise statistics and the relationship between image noise variance and expected pixel intensity. Image denoising is essential to improve quality of fluoroscopic images and their clinical information content. Simple average filters are commonly employed in real-time processing, but they tend to blur edges and details. An extensive comparison of advanced denoising algorithms specifically designed for both signal-dependent noise (AAS, BM3Dc, HHM, TLS) and independent additive noise (AV, BM3D, K-SVD) was presented. Simulated test images degraded by various levels of Poisson quantum noise and real clinical fluoroscopic images were considered. Typical gray-level transformations (e.g. white compression) were also applied in order to evaluate their effect on the denoising algorithms. Performances of the algorithms were evaluated in terms of peak-signal-to-noise ratio (PSNR), signal-to-noise ratio (SNR), mean square error (MSE), structural similarity index (SSIM) and computational time. On average, the filters designed for signal-dependent noise provided better image restorations than those assuming additive white Gaussian noise (AWGN). Collaborative denoising strategy was found to be the most effective in denoising of both simulated and real data, also in the presence of image gray-level transformations. White compression, by inherently reducing the greater noise variance of brighter pixels, appeared to support denoising algorithms in performing more effectively.     

Characterizing Venous Vasculatures of Hepatocellular Carcinoma Using a Multi-Breath-Hold Two-Dimensional Susceptibility Weighted Imaging

The aim of our study is to characterize the venous vasculatures of hepatocellular carcinoma (HCC) using a multi-breath-hold two-dimensional (2D) susceptibility weighted imaging (SWI) in comparison with conventional Magnetic Resonance Imaging (MRI) sequences. Twenty-nine patients with pathologically confirmed HCC underwent MR examination at a 3.0 T scanner. The number of venous vascularity in or around the lesion was counted and the image quality was subjectively evaluated by two experienced radiologists independently based on four image sets: 1) SWI, 2) T1-weighted sequence, 3) T2-weighted sequence, and 4) T1-weighted dynamic contrast-enhanced (DCE) sequence. Of the 29 patients, a total of 33 liver lesions were detected by both SWI and conventional MR sequences. In the evaluation of the conspicuity of venous vascularity, a mean of 10.7 tumor venous vessels per mass was detected by the SWI and 3.9 tumor vasculatures were detected by T1-weighted DCE (P<0.0001), while none was detected by T1-, T2-weighted sequences. The Pearson correlation coefficients between the lesion sizes and the number of tumor vasculatures detected by T1-weighted DCE was 0.708 (P<0.001), and 0.883 by SWI (P<0.001). Our data suggest that SWI appears to be a more sensitive tool compared to T1-weighted DCE sequence to characterize venous vasculature in liver lesions.     

N2NSR-OCT: Simultaneous denoising and super-resolution in optical coherence tomography images using semisupervised deep learning

Optical coherence tomography (OCT) imaging shows a significant potential in clinical routines due to its noninvasive property. However, the quality of OCT images is generally limited by inherent speckle noise of OCT imaging and low sampling rate. To obtain high signal-to-noise ratio (SNR) and high-resolution (HR) OCT images within a short scanning time, we presented a learning-based method to recover high-quality OCT images from noisy and low-resolution OCT images. We proposed a semisupervised learning approach named N2NSR-OCT, to generate denoised and super-resolved OCT images simultaneously using up- and down-sampling networks (U-Net (Semi) and DBPN (Semi)). Additionally, two different super-resolution and denoising models with different upscale factors (2× and 4× ) were trained to recover the high-quality OCT image of the corresponding down-sampling rates. The new semisupervised learning approach is able to achieve results comparable with those of supervised learning using up- and down-sampling networks, and can produce better performance than other related state-of-the-art methods in the aspects of maintaining subtle fine retinal structures.     

Evaluation of radiomic texture feature error due to MRI acquisition and reconstruction: A simulation study utilizing ground truth

The purpose of this study was to examine the dependence of image texture features on MR acquisition parameters and reconstruction using a digital MR imaging phantom. MR signal was simulated in a parallel imaging radiofrequency coil setting as well as a single element volume coil setting, with varying levels of acquisition noise, three acceleration factors, and four image reconstruction algorithms. Twenty-six texture features were measured on the simulated images, ground truth images, and clinical brain images. Subtle algorithm-dependent errors were observed on reconstructed phantom images, even in the absence of added noise. Sources of image error include Gibbs ringing at image edge gradients (tissue interfaces) and well-known artifacts due to high acceleration; two of the iterative reconstruction algorithms studied were able to mitigate these image errors. The difference of the texture features from ground truth, and their variance over reconstruction algorithm and parallel imaging acceleration factor, were compared to the clinical “effect size”, i.e., the feature difference between high- and low-grade tumors on T1- and T2-weighted brain MR images of twenty glioma patients. The measured feature error (difference from ground truth) was small for some features, but substantial for others. The feature variance due to reconstruction algorithm and acceleration factor were generally smaller than the clinical effect size. Certain texture features may be preserved by MR imaging, but adequate precautions need to be taken regarding their validity and reliability. We present a general simulation framework for assessing the robustness and accuracy of radiomic textural features under various MR acquisition/reconstruction scenarios.     

Elimination of biological and physical artifacts in abdomen and brain computed tomography procedures using filtering techniques

IntroductionMedical images are usually affected by biological and physical artifacts or noise, which reduces image quality and hence poses difficulties in visual analysis, interpretation and thus requires higher doses and increased radiographs repetition rate.ObjectivesThis study aims at assessing image quality during CT abdomen and brain examinations using filtering techniques as well as estimating the radiogenic risk associated with CT abdomen and brain examinations.Materials and MethodsThe data were collected from the Radiology Department at Royal Care International (RCI) Hospital, Khartoum, Sudan. The study included 100 abdominal CT images and 100 brain CT images selected from adult patients. Filters applied are namely: Mean filter, Gaussian filter, Median filter and Minimum filter. In this study, image quality after denoising is measured based on the Mean Squared Error (MSE), Peak Signal-to-Noise Ratio (PSNR), and the Structural Similarity Index Metric (SSIM).ResultsThe results show that the images quality parameters become higher after applications of filters. Median filter showed improved image quality as interpreted by the measured parameters: PSNR and SSIM, and it is thus considered as a better filter for removing the noise from all other applied filters.DiscussionThe noise removed by the different filters applied to the CT images resulted in enhancing high quality images thereby effectively revealing the important details of the images without increasing the patients’ risks from higher doses.ConclusionsFiltering and image reconstruction techniques not only reduce the dose and thus the radiation risks, but also enhances high quality imaging which allows better diagnosis.     

The discriminative bilateral filter: an enhanced denoising filter for electron microscopy data

Advances in three-dimensional (3D) electron microscopy (EM) and image processing are providing considerable improvements in the resolution of subcellular volumes, macromolecular assemblies and individual proteins. However, the recovery of high-frequency information from biological samples is hindered by specimen sensitivity to beam damage. Low dose electron cryo-microscopy conditions afford reduced beam damage but typically yield images with reduced contrast and low signal-to-noise ratios (SNRs). Here, we describe the properties of a new discriminative bilateral (DBL) filter that is based upon the bilateral filter implementation of Jiang et al. (Jiang, W., Baker, M.L., Wu, Q., Bajaj, C., Chiu, W., 2003. Applications of a bilateral denoising filter in biological electron microscopy. J. Struc. Biol. 128, 82-97.). In contrast to the latter, the DBL filter can distinguish between object edges and high-frequency noise pixels through the use of an additional photometric exclusion function. As a result, high frequency noise pixels are smoothed, yet object edge detail is preserved. In the present study, we show that the DBL filter effectively reduces noise in low SNR single particle data as well as cellular tomograms of stained plastic sections. The properties of the DBL filter are discussed in terms of its usefulness for single particle analysis and for pre-processing cellular tomograms ahead of image segmentation.     

Imaging of Reperfused Intramyocardial Hemorrhage with Cardiovascular Magnetic Resonance Susceptibility Weighted Imaging (SWI)

PurposeTo report initial experience with TE-averaged susceptibility weighted imaging (SWI) in normal subjects and acute myocardial infarction (AMI) patients for the detection of intramyocardial hemorrhage (IMH).ResultsThere were six patients with microvascular obstruction (MVO) and four patients with IMH detected by TE-averaged SWI imaging. All patients with IMH on SWI scans had MVO on late gadolinium-enhanced (LGE) imaging. There was a three-fold increase in IMH contrast with SWI compared to magnitude images. IMH contrast decreased and signal-to-noise increased with increased TE averages.ConclusionsTE-averaged SWI imaging is a promising method for myocardial tissue characterization in the setting of AMI for the detection of IMH. Along with gray-scale colormap inversion, it combines not only magnitude and phase information, but also images across TEs to provide a single image sensitive to IMH with characteristics similar to LGE imaging.     

Applications of a bilateral denoising filter in biological electron microscopy

Due to the sensitivity of biological sample to the radiation damage, the low dose imaging conditions used for electron microscopy result in extremely noisy images. The processes of digitization, image alignment, and 3D reconstruction also introduce additional sources of noise in the final 3D structure. In this paper, we investigate the effectiveness of a bilateral denoising filter in various biological electron microscopy applications. In contrast to the conventional low pass filters, which inevitably smooth out both noise and structural features simultaneously, we found that bilateral filter holds a distinct advantage in being capable of effectively suppressing noise without blurring the high resolution details. In as much, we have applied this technique to individual micrographs, entire 3D reconstructions, segmented proteins, and tomographic reconstructions.     

Quantitative Susceptibility Mapping-Based Microscopy of Magnetic Resonance Venography (QSM-mMRV) for In Vivo Morphologically and Functionally Assessing Cerebromicrovasculature in Rat Stroke Model

Abnormal cerebral oxygenation and vessel structure is a crucial feature of stroke. An imaging method with structural and functional information is necessary for diagnosis of stroke. This study applies QSM-mMRV (quantitative susceptibility mapping-based microscopic magnetic resonance venography) for noninvasively detecting small cerebral venous vessels in rat stroke model. First, susceptibility mapping is optimized and calculated from magnetic resonance (MR) phase images of a rat brain. Subsequently, QSM-mMRV is used to simultaneously provide information on microvascular architecture and venous oxygen saturation (SvO₂), both of which can be used to evaluate the physiological and functional characteristics of microvascular changes for longitudinally monitoring and therapeutically evaluating a disease model. Morphologically, the quantification of vessel sizes using QSM-mMRV was 30% smaller than that of susceptibility-weighted imaging (SWI), which eliminated the overestimation of conventional SWI. Functionally, QSM-mMRV estimated an average SvO₂ ranging from 73% to 85% for healthy rats. Finally, we also applied QSM to monitor the revascularization of post-stroke vessels from 3 to 10 days after reperfusion. QSM estimations of SvO₂ were comparable to those calculated using the pulse oximeter standard metric. We conclude that QSM-mMRV is useful for longitudinally monitoring blood oxygen and might become clinically useful for assessing cerebrovascular diseases.     

A Feasibility Study of an Intravascular Imaging Antenna to Image Atherosclerotic Plaques in Swine Using 3.0 T MRI

Purpose

To investigate the feasibility of an intravascular imaging antenna to image abdominal aorta atherosclerotic plaque in swine using 3.0T magnetic resonance imaging (MRI).

Methods

Atherosclerotic model was established in 6 swine. After 8 months, swine underwent an MR examination, which was performed using an intravascular imaging guide-wire, and images of the common iliac artery and the abdominal aorta were acquired. Intravascular ultrasound (IVUS) was performed in the right femoral artery; images at the same position as for the MR examination were obtained. The luminal border and external elastic membrane of the targeted arteries were individually drawn in the MR and IVUS images. After co-registering these images, the vessel, lumen, and vessel wall areas and the plaque burden in the same lesions imaged using different modalities were calculated and compared. The diagnostic accuracy of intravascular MR examination in delineating the vessel wall and detecting plaques were analyzed and compared using IVUS.

Results

Compared with IVUS, good agreement was found between MRI and IVUS for delineating vessel, lumen, and vessel wall areas and plaque burden (r value: 0.98, 0.95, 0.96 and 0.91, respectively; P<0.001).

Conclusion

Compared with IVUS, using an intravascular imaging guide-wire to image deep seated arteries allowed determination of the vessel, lumen and vessel wall areas and plaque size and burden. This may provide an alternative method for detecting atherosclerotic plaques in the future.     

一种滤除医学影像噪声的混合滤波算法

目的：医学影像在获取、存储、传输过程中会不同程度地受到噪声污染,这极大影像了其在临床诊疗中的应用。为了有效地滤除医学影像噪声,提出了一种混合滤波算法。方法：该算法首先将含有高斯和椒盐噪声的图像进行形态学开运算,然后对开运算后的图像进行二维小波分解,得到高频和低频小波分解系数。保留低频系数不变,将高频系数经过维纳滤波器进行滤波,最后进行小波系数重构。结果：采用该混合滤波算法、小波阚值去噪、中值滤波、维纳滤波分别对含有混合噪声的医学影像分别进行滤除噪声处理,该滤波算法去噪后影像的PSNR值明显高于其他三种方法。结论：该混合滤波算法是一种较为有效的医学影像噪声滤除方法。     

Development of a dynamic imaging method for gravitropism in pea sprouts using clinical magnetic resonance imaging system

Although magnetic resonance imaging (MRI) is a useful technique, only a few studies have investigated the dynamic behavior of small subjects using MRI owing to constraints such as experimental space and signal amount. In this study, to acquire high-resolution continuous three-dimensional gravitropism data of pea (Pisum sativum) sprouts, we developed a small-bore MRI signal receiver coil that can be used in a clinical MRI and adjusted the imaging sequence. It was expected that such an arrangement would improve signal sensitivity and improve the signal-to-noise ratio (SNR) of the acquired image. All MRI experiments were performed using a 3.0-T clinical MRI scanner. An SNR comparison using an agarose gel phantom to confirm the improved performance of the small-bore receiver coil and an imaging experiment of pea sprouts exhibiting gravitropism were performed. The SNRs of the images acquired with a standard 32-channel head coil and the new small-bore receiver coil were 5.23±0.90 and 57.75±12.53, respectively. The SNR of the images recorded using the new coil was approximately 11-fold higher than that of the standard coil. In addition, when the accuracy of MR imaging that captures the movement of pea sprout was verified, the difference in position information from the optical image was found to be small and could be used for measurements. These results of this study enable the application of a clinical MRI system for dynamic plant MRI. We believe that this study is a significant first step in the development of plant MRI technique.     

Performance of deep learning restoration methods for the extraction of particle dynamics in noisy microscopy image sequences

Particle tracking in living systems requires low light exposure and short exposure times to avoid phototoxicity and photobleaching and to fully capture particle motion with high-speed imaging. Low-excitation light comes at the expense of tracking accuracy. Image restoration methods based on deep learning dramatically improve the signal-to-noise ratio in low-exposure data sets, qualitatively improving the images. However, it is not clear whether images generated by these methods yield accurate quantitative measurements such as diffusion parameters in (single) particle tracking experiments. Here, we evaluate the performance of two popular deep learning denoising software packages for particle tracking, using synthetic data sets and movies of diffusing chromatin as biological examples. With synthetic data, both supervised and unsupervised deep learning restored particle motions with high accuracy in two-dimensional data sets, whereas artifacts were introduced by the denoisers in three-dimensional data sets. Experimentally, we found that, while both supervised and unsupervised approaches improved tracking results compared with the original noisy images, supervised learning generally outperformed the unsupervised approach. We find that nicer-looking image sequences are not synonymous with more precise tracking results and highlight that deep learning algorithms can produce deceiving artifacts with extremely noisy images. Finally, we address the challenge of selecting parameters to train convolutional neural networks by implementing a frugal Bayesian optimizer that rapidly explores multidimensional parameter spaces, identifying networks yielding optimal particle tracking accuracy. Our study provides quantitative outcome measures of image restoration using deep learning. We anticipate broad application of this approach to critically evaluate artificial intelligence solutions for quantitative microscopy.     

An iterative deconvolution algorithm for image recovery in clinical CT: A phantom study

PurposeTo study the feasibility of using an iterative reconstruction algorithm to improve previously reconstructed CT images which are judged to be non-diagnostic on clinical review. A novel rapidly converging, iterative algorithm (RSEMD) to reduce noise as compared with standard filtered back-projection algorithm has been developed.Materials and methodsThe RSEMD method was tested on in-silico, Catphan^®500, and anthropomorphic 4D XCAT phantoms. The method was applied to noisy CT images previously reconstructed with FBP to determine improvements in SNR and CNR. To test the potential improvement in clinically relevant CT images, 4D XCAT phantom images were used to simulate a small, low contrast lesion placed in the liver.ResultsIn all of the phantom studies the images proved to have higher resolution and lower noise as compared with images reconstructed by conventional FBP. In general, the values of SNR and CNR reached a plateau at around 20 iterations with an improvement factor of about 1.5 for in noisy CT images. Improvements in lesion conspicuity after the application of RSEMD have also been demonstrated. The results obtained with the RSEMD method are in agreement with other iterative algorithms employed either in image space or with hybrid reconstruction algorithms.ConclusionsIn this proof of concept work, a rapidly converging, iterative deconvolution algorithm with a novel resolution subsets-based approach that operates on DICOM CT images has been demonstrated. The RSEMD method can be applied to sub-optimal routine-dose clinical CT images to improve image quality to potentially diagnostically acceptable levels.     

A deep‐learning‐based approach for noise reduction in high‐speed optical coherence Doppler tomography

Optical coherence Doppler tomography (ODT) increasingly attracts attention because of its unprecedented advantages with respect to high contrast, capillary‐level resolution and flow speed quantification. However, the trade‐off between the signal‐to‐noise ratio of ODT images and A‐scan sampling density significantly slows down the imaging speed, constraining its clinical applications. To accelerate ODT imaging, a deep‐learning‐based approach is proposed to suppress the overwhelming phase noise from low‐sampling density. To handle the issue of limited paired training datasets, a generative adversarial network is performed to implicitly learn the distribution underlying Doppler phase noise and to generate the synthetic data. Then a 3D based convolutional neural network is trained and applied for the image denoising. We demonstrate this approach outperforms traditional denoise methods in noise reduction and image details preservation, enabling high speed ODT imaging with low A‐scan sampling density.     

An Objective Method to Optimize the MR Sequence Set for Plaque Classification in Carotid Vessel Wall Images Using Automated Image Segmentation

A typical MR imaging protocol to study the status of atherosclerosis in the carotid artery consists of the application of multiple MR sequences. Since scanner time is limited, a balance has to be reached between the duration of the applied MR protocol and the quantity and quality of the resulting images which are needed to assess the disease. In this study an objective method to optimize the MR sequence set for classification of soft plaque in vessel wall images of the carotid artery using automated image segmentation was developed. The automated method employs statistical pattern recognition techniques and was developed based on an extensive set of MR contrast weightings and corresponding manual segmentations of the vessel wall and soft plaque components, which were validated by histological sections. Evaluation of the results from nine contrast weightings showed the tradeoff between scan duration and automated image segmentation performance. For our dataset the best segmentation performance was achieved by selecting five contrast weightings. Similar performance was achieved with a set of three contrast weightings, which resulted in a reduction of scan time by more than 60%. The presented approach can help others to optimize MR imaging protocols by investigating the tradeoff between scan duration and automated image segmentation performance possibly leading to shorter scanning times and better image interpretation. This approach can potentially also be applied to other research fields focusing on different diseases and anatomical regions.     

A C-GAN Denoising Algorithm in Projection Domain for Micro-CT

Micro-CT provides a high-resolution 3D imaging of micro-architecture in a non-invasive way, which becomes a significant tool in biomedical research and preclinical applications. Due to the limited power of micro-focus X-ray tube, photon starving occurs and noise is inevitable for the projection images, resulting in the degradation of spatial resolution, contrast and image details. In this paper, we propose a C-GAN (Conditional Generative Adversarial Nets) denoising algorithm in projection domain for Micro-CT imaging. The noise statistic property is utilized directly and a novel variance loss is developed to suppress the blurry effects during denoising procedure. Conditional Generative Adversarial Networks (C-GAN) is employed as a framework to implement the denoising task. To guarantee the pixelwised accuracy, fully convolutional network is served as the generator structure. During the alternative training of the generator and the discriminator, the network is able to learn noise distribution automatically. Moreover, residual learning and skip connection architecture are applied for faster network training and further feature fusion. To evaluate the denoising performance, mouse lung, milkvetch root and bamboo stick are imaged by micro-CT in the experiments. Compared with BM3D, CNN-MSE and CNN-VGG, the proposed method can suppress noise effectively and recover image details without introducing any artifacts or blurry effect. The result proves that our method is feasible, efficient and practical.     

早期胚胎全场OCT图像的边界去模糊研究

全场光学相干层析成像技术(全场OCT)是研究早期胚胎形态发育的最理想成像设备,然而所采集图像难免受噪声干扰.这些噪声可模糊早期胚胎内不同组织结构的边界,从而给基于图像边界的结构划分带来干扰.为解决这一问题,本文运用中值滤波、维纳滤波、各向异性扩散算法处理全场OCT获得的早期胚胎图像,并运用信噪比、均方误差、峰值信噪比和边缘保留等指标评价图像处理效果.结果表明:经各向异性扩散算法处理的早期胚胎图像,可完整地保留原始图像信息,且边界最清晰,视觉效果最好.     

Fast denoising and lossless spectrum extraction in stimulated Raman scattering microscopy

Stimulated Raman scattering (SRS) microscopy is a nonlinear optical imaging method for visualizing chemical content based on molecular vibrational bonds. However, the imaging speed and sensitivity are currently limited by the noise of the light beam probing the Raman process. In this paper, we present a fast non-average denoising and high-precision Raman shift extraction method, based on a self-reinforcing signal-to-noise ratio (SNR) enhancement algorithm, for SRS spectroscopy and microscopy. We compare the results of this method with the filtering methods and the reported experimental methods to demonstrate its high efficiency and high precision in spectral denoising, Raman peak extraction and image quality improvement. We demonstrate a maximum SNR enhancement of 10.3 dB in fixed tissue imaging and 11.9 dB in vivo imaging. This method reduces the cost and complexity of the SRS system and allows for high-quality SRS imaging without use of special laser, complicated system design and Raman tags.