Next Article in Journal
Exposing Manipulated Photos and Videos in Digital Forensics Analysis
Next Article in Special Issue
Customized Efficient Neural Network for COVID-19 Infected Region Identification in CT Images
Previous Article in Journal
Comparing Stacking Ensemble Techniques to Improve Musculoskeletal Fracture Image Classification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Super Resolution of Magnetic Resonance Images

1
Indian Institute of Technology Mandi, Mandi, Himachal Pradesh 175005, India
2
Post Graduate Institute of Medical Education & Research, Chandigarh 160012, India
*
Author to whom correspondence should be addressed.
Submission received: 8 April 2021 / Revised: 4 June 2021 / Accepted: 7 June 2021 / Published: 21 June 2021
(This article belongs to the Special Issue Intelligent Strategies for Medical Image Analysis)

Abstract

:
In this work, novel denoising and super resolution (SR) approaches for magnetic resonance (MR) images are addressed, and are integrated in a unified framework, which do not require example low resolution (LR)/high resolution (HR)/cross-modality/noise-free images and prior information of noise–noise variance. The proposed method categorizes the patches as either smooth or textured and then denoises them by deploying different denoising strategies for efficient denoising. The denoising algorithm is integrated into the SR approach, which uses a gradient profile-based constraint in a sparse representation-based framework to improve the resolution of MR images with reduced smearing of image details. This constraint regularizes the estimation of HR images such that the estimated HR image has gradient profiles similar to the gradient profiles of the original HR image. For this, the gradient profile sharpness (GPS) values of an unknown HR image are estimated using an approximated piece-wise linear relation among GPS values of LR and upsampled LR images. The experiments are performed on three different publicly available datasets. The proposed SR approach outperforms the existing unsupervised SR approach addressed for real MR images that exploits low rank and total variation (LRTV) regularization, by an average peak signal to noise ratio (PSNR) of 0.73 dB and 0.38 dB for upsampling factors 2 and 3, respectively. For the super resolution of noisy real MR images (degraded with 2% noise), the proposed approach outperforms the LRTV approach by an average PSNR of 0.54 dB and 0.46 dB for upsampling factors 2 and 3, respectively. The qualitative analysis is shown for real MR images from healthy subjects and subjects with Alzheimer’s disease and structural deformity, i.e., cavernoma.

1. Introduction

High quality magnetic resonance (MR) images are generally desired for precise medical diagnosis and analysis, and are typically characterized by high spatial resolution and high signal to noise ratio (SNR). However, the acquisition of such high quality MR images is limited due to various constraints such as acquisition time and requirement of expensive hardware. For example, an increase in the spatial resolution is limited as it leads to a decrease in SNR in MR images for a fixed acquisition time. It has been experimentally demonstrated that the post-acquisition MR image processing can increase the spatial resolution of MR images along with improved SNR [1,2]. The focus of this work is to improve the resolution (using a super resolution (SR) approach) and SNR (using a denoising approach) of MR images.
Different paradigms of machine learning approaches: The limited availability of example/training MR images makes it crucial to classify the existing denoising and SR approaches based on the requirement of paired/labeled example data as (i) supervised and (ii) unsupervised. The supervised approaches require paired example images for learning the mapping between input and ground truth images. The unsupervised approaches, on the other hand, are not supervised by the paired/labeled data. Several unsupervised approaches have been addressed in the literature which require unpaired/unlabeled data to learn the mapping between input and ground truth images. For example, the cycle-generative adversarial networks (GAN)-based unsupervised SR approach exploits transductive learning and requires example low resolution (LR) and high resolution (HR) images without any correspondences among them [3]. Many existing unsupervised SR approaches require either example LR or example HR or example cross-modality MR images [4,5,6,7,8,9,10,11]. Another class of unsupervised approaches assumes a mathematical model between input and desired output images, and poses the problem as ill-posed inverse problems [12,13,14,15,16,17]. The main advantage of these approaches is that these approaches do not require any example images. However, there might be other requirements in these approaches such as several denoising approaches from this class which obviate the requirement of example images and require the prior knowledge of noise: noise variance or distribution of noise [14,15,16,17]. The denoising approaches which exploit prior knowledge of noise have been observed to outperform the approaches which do not require the prior knowledge of noise, but such requirements of noise are generally not met because the variance of noise is often not known a priori. Several hybrid versions for supervised and unsupervised frameworks also exist such as a semi-supervised framework which requires a few paired and many unpaired example MR images, for example cycle GAN-based methods [18].
In the last decades, several deep learning-based approaches have been addressed to super resolve the MR images [7,19,20]. These methods use three dimensional convolutional neural networks in [7] to learn the mapping among given LR and HR images to estimate HR images. Such methods are further investigated to estimate an HR image along with learning the residual/difference in estimated HR and ground truth images to improve the estimation of HR images [19]. A similar residual learning-based strategy is used for super resolving the musculoskeletal-knee MR images in the slice-select direction [20]. Most of the addressed deep learning-based methods are supervised and the recent research is inclined towards reducing the requirement of example images [18].
The performances of supervised, semi-supervised, and unsupervised approaches that exploit example images are better than the latter class of unsupervised approaches without the requirement of any example images, but the applicability of these approaches is highly restricted due to the limited availability of example images in real scenarios [3,18]. Further, the performance of supervised approaches is highly dependent on the training database [21]. These issues can be addressed by unsupervised approaches, which obviate all requirements of example data, but there are very few in the literature [12,13]. Therefore, in this direction, we focus on SR and denoising approaches in an unsupervised framework, which obviates the additional requirements such as prior-knowledge or estimation of noise variance and the distribution of noise, along with the requirement of example images. For simplicity, the latter class of unsupervised approaches, which do not require any example images, are here onwards referred to as unsupervised approaches.
Challenges in unsupervised image denoising: Image denoising specifically for MR images is well explored using Bayesian approaches [22], partial differential equation [23], adaptive smoothing and filtering [24], collaborative filtering [25,26,27], etc. The low-rank approximation (LRA) of a given image patch estimated by either a nuclear norm or singular value decomposition (SVD) or principal component analysis (PCA) has been extensively studied for denoising of MR images [28,29,30,31]. The LRA-based image denoising methods, which use PCA to provide the first few (l) eigenvectors with significant eigenvalues that are used to represent denoised image patches, can be categorized as unsupervised approaches. The performance of these methods depends upon the criterion to choose the appropriate l to discard l + 1 onward eigenvectors, which are mainly associated with noise present in the given image patches. Using the same criterion (e.g., the same l) to select the eigenvectors for denoising different kinds of image patches (such as smooth patches, and patches with varying strengths of edges) may not be suitable to improve the trade-off between suppression of noise and smearing of edges in images. The issue of the same criterion to choose l for different patch structures has been addressed by choosing l adaptive to the patch structure as in [14,15,16,17]. These approaches address the computation of weights that are adaptive to the patch structure. These weights are used to weigh the projection coefficients of image patches on eigenvectors to reconstruct the denoised image patches. Hence, it leads to efficient denoising of an image. For smooth/flat image patches, the approaches in [14,17] provide exactly zero weight, and when multiplied with projection coefficients (producing l = 0 ), they lead to better denoising of smooth patches. However, these approaches require prior knowledge of variance of noise for the computation of adaptive weights. Other approaches [15,16] do not require prior knowledge of noise but in these approaches, the value computed for l is usually greater than or equal to one even for smooth patches and thus may result in noisy smooth patches. This concern is addressed by the proposed denoising approach in the unsupervised framework by categorizing the smooth and textured patches without the requirements of prior knowledge of noise variance and any example image. The categorized patches are denoised using different denoising strategies. This leads to effective denoising of smooth patches similar to [14,17] and better than [15,16] but without the requirement of noise variance in the proposed denoising approach.
Challenges in unsupervised SR: The initial attempt made in unsupervised SR of MR images is based on the non-local means (NLM3D) [13]. In this approach, each small 3D cube volume is represented using weighted averaging of similar 3D cubes present non-locally in interpolated MR image volume. Here, the weighted averaging of image patches generally tends to blur the image details. Though weighted averaging aids in its robustness to noise while super-resolving the noisy images, it can be effective for only a small strength of noise. Feng Shi et al. exploit the low-rank structure of brain MR image patches along with the total variation regularization in MR images low rank and total variation (LRTV) [12]. This approach reduces the blur in reconstructed image details but is observed to have a tendency to produce staircase artifacts in reconstructed MR images (a staircase artifact is generally associated with TV prior [32]). On the other hand, the degradation of gradient profiles in MR images differs from natural images. Thus, such a difference shall be considered while deploying the gradient-based prior in the SR of MR images. The proposed SR approach addresses this difference and utilizes it to regularize the estimation of HR images such that the gradient profiles of an HR image (depicted by tissue boundaries) are restored with better accuracy and without the significant generation of artifacts like the staircase artifact.
Scope of this work: This paper addresses the SR of noisy brain MR images by proposing novel unsupervised denoising and a super resolution approach, which are embedded in a unified framework. The proposed denoising approach is based on the hypothesis that the trade-off between suppression of noise and smearing of image details can be improved if different structures of patches are denoised with different strategies. The categorization of different patch structures, such as smooth patches and textured/edge patches, from noisy image patches is a challenging task, especially without the prior knowledge of noise variance. The proposed work addresses such patch categorization by exploiting the relation between different structures of MR image patches and their respective low rank approximation (LRA)-based reconstruction error. The LRA-based reconstruction error of a patch here denotes the mean square error between the patch and its low-rank approximated version. Further, we propose to use a statistical measure of LRA-based reconstruction error values of an ensemble of image patches to indicate the strength of noise, and is used to categorize the patch structures adaptive to the noise strength. The eigenvectors used to compute LRA of the image patches are estimated from the ensemble of given noisy image patches. Thus, the estimated low-rank structure of image patches possesses inaccuracy. To address this, we progressively reduce such inaccuracy by re-estimating the LRA of patches from the denoised image patches in an iterative manner. This aids in the accurate reconstruction of image details. On the other hand, the gradient profiles of image details convey important semantic information and are generally degraded while improving the resolution. The gradient information thus shall be preserved to aid the reconstruction of clear tissue boundaries in MR images. To do so, we propose to regularize the estimation of the HR image in the SR approach with approximated reference gradient profile sharpness (GPS) values. The unavailable reference GPS values, i.e., GPS values of an HR image, are approximated by estimating a piecewise linear relation (for natural images the relation is linear) between GPS values of available LR and upsampled LR images. The super resolution approach addressed in this work exploits the sparse representation framework and uses the self-similarity concept to utilize similarity among patches of the test image itself for enriching the dictionary construction.
The novelty of this work is as follows: (i) a novel denoising approach is proposed that addresses the categorization of noisy image patches as smooth or textured/edges without the estimated/prior knowledge of noise such as noise variance, and (ii) a novel gradient profile sharpness-based constraint is proposed specific to the structure of MR images in the super resolution approach.
The preliminary work by us for in-plane (2D) super resolution of MR images with GPS constraint has been reported in [33]. The work in this paper has been further extended by defining a novel denoising approach that addresses the inhomogeneity among image patches and formulates a noise adaptive approach for categorization of different patch structures and for improving the resolution of noisy three dimensional (3D) MR volume. The extensive experimental analysis is reported in this manuscript including the analysis of the parameters in the proposed approach, and the performance on both healthy and unhealthy subjects.
Contributions of our work: The important contributory aspects of the proposed work are as follows: (i) The relation between structure of a patch and its LRA version is exploited to categorize the patch structure without the requirements of prior information about noise, (ii) the statistical mean of LRA-based error values of patches is used to indicate the strength of noise and is used to categorize the patches adaptive to the strength of noise, (iii) a progressive improvement in the estimated low-rank structure of image patches is addressed by re-estimating the eigenvectors from denoised image patches for better reconstruction of image details, (iv) the piecewise linear relation among GPS values of LR and an upsampled LR image is formulated to regularize the estimation of a denoised HR image that aids in estimation of clearer tissue boundaries, (v) extensive experimental analysis is performed for healthy and unhealthy subjects to validate the accurate estimation of MR image details.
The organization of the paper is as follows: Section 2 is an overview of related SR and denoising approaches. The proposed denoising approach and SR approach are discussed in Section 3. The experimental results are demonstrated and discussed in Section 4. The work is summarized in Section 5.

2. Related Works

The SR of noisy MR images has been addressed in the literature in two ways: (i) using a disjoint combination of SR and denoising approaches (first denoise and then super resolve the images or vice-versa) [12,13,25,34,35], and (ii) using the embedded framework for SR and denoising approaches [21,36]. It has been demonstrated that the image details are restored with better accuracy in the case of the embedded framework than the disjoint combination [21,36]. The existing works to improve the quality of MR images are explained in three categories: approaches that address (i) both SR and denoising in a single framework, (ii) denoising only, and (iii) SR only, in the following Section 2.1, Section 2.2 and Section 2.3, respectively. Here, we only discuss the closely related approaches in detail.

2.1. Embedded Approaches for SR and Denoising

The embedded frameworks integrate the denoising approach in the SR approach using the framework of either sparse representation of image patches [15,21,36,37] or neural networks, and thus omit the drawbacks of the disjoint combination of SR and denoising approaches. Numerous ways to shrink the coefficients in the wavelet domain [24] and PCA [15] have been reported in the literature to denoise while super resolving the images, but require prior knowledge of noise variance. In the case of medical images, especially for MR images, D. Trinh et al. proposed a single framework for SR and denoising [21]. Here, each HR image patch is estimated from the linear combination of a few selected example HR image patches, for the given LR image patch. However, this method requires example HR images and is thus supervised in nature with limited applicability in the real scenario. The approach in [36] addressed a sparse-representation framework for SR of noisy MR images without the requirement of example images. However, it follows the denoising approach as in [14] and thus requires prior knowledge of the strength of noise, i.e., noise variance.

2.2. Image Denoising Approaches

In magnitude MR images (such as T1 and T2 weighted MR images), the noise follows a Rician distribution [38]. However, for high SNR regions (pixels/patches with brain tissues of high signal intensities) the distribution of noise can be well approximated by a white Gaussian distribution [38,39]. Hence, the denoising of MR images for both Gaussian distributed noise (for high SNR regions) [28,29,40,41] and Rician distributed noise along with mean removal [25,42] is addressed in the literature. Further, many image denoising techniques in the MR images’ literature process the 3D image’s volume [27,28,35,41]. That means while denoising one MR image, the details from its adjacent slices are also considered. Apart from 3D level processing, there exist MR image denoising methods which process each image independently (2D) and do not consider the information from other slices in the same volume [34]. In LRA-based unsupervised denoising approaches, the low rank structure estimated by eigenvectors possess inaccuracy. This is because the eigenvectors, employed to reconstruct the denoised image, are estimated from the given noisy image patches (especially in unsupervised framework) and thus may not represent the reliable low-rank structure of image patches. This is addressed by an LRA-based denoising approach which is closely related to our method, and is addressed in computer vision literature, i.e., local patch group-based PCA [14]. In this method, similar patches are grouped and eigenvectors are estimated from the covariance matrix of patches in each cluster. The eigenvectors are re-estimated from the denoised images to improve the accuracy of estimated low rank structure of image patches. The projection coefficients for a given patch on the eigenvectors are adaptively weighted according to the structure of the patch for efficient denoising. However, this approach requires the prior knowledge of noise variance for the computation of these adaptive weights, in contrast to the proposed denoising approach.

2.3. Super Resolution Approaches

The SR approaches to improve the resolution of MR images have been well explored using different methods [4,5,6,7,9,10,12,13,43,44,45,46,47,48,49,50], and in different frameworks like sparse representation framework, deep neural networks, etc. [2,6,7,8,11,49]. Various regularizers are explored in these frameworks such as sparse derivative priors [50], gradient [49], and structural priors [47] for better restoration of HR image details. The low rank and total variation priors used in [12] aid in reconstruction of sharp image details but have been observed to induce the staircase effect [32]. The SR approach is closely related to the proposed SR approach [51] in that it addresses GPS-based approximation of the gradient field for natural images and uses the approximated gradient field to constrain the estimation of the HR image. The relation between GPS values of LR and upsampled LR natural images is approximated by linear relation and is used to estimate GPS values of the HR image. However, this approach addresses the GPS values obtained using fixed shape edge models like the triangle and mixed Gaussian edge models. These models may not fit the edges/gradient profiles of MR images. Thus, the proposed work utilizes the similar formulation of GPS from [51], but is defined using a more generalized edge model [52] for MR images instead of fixed shape edge models. Further, it has been observed that the GPS values for LR and upsampled LR MR images are not linearly related as in natural images but are piecewise linearly related [33]. We have focused on the gradient-based prior specific to MR images which addresses the above-mentioned difference between MR and natural images.

3. Method

For the evaluation of the proposed method, both simulated as well as real T1 weighted brain MR images are used. The simulated images are obtained from MR images dataset-I (https://brainweb.bic.mni.mcgill.ca/ accessed on 1 December 2019) and eight real MR image volumes are randomly selected from the publicly available dataset-II (https://www.humanconnectome.org/study/hcp-young-adult/data-releases accessed on on 1 December 2019). The MR images are obtained from dataset-I using a normal brain database with slice thickness 1 mm and 0% intensity non-uniformity. The selected MR images from dataset-II are acquired from a 3T MRI scanner using Magnetization Prepared Rapid Acquisition Gradient Echo (MPRAGE) sequence, time to repeat (TR) = 2400 ms, time to echo (TE) = 2.14, time inversion (TI) = 1000 ms, and flip angle (FA) = 8°. The spatial resolutions of each volume in dataset-I and dataset-II are 181 × 217 × 181 and 173 × 173 × 207 , respectively, which are used to simulate noisy LR images, and then super resolved and/or denoised using different algorithms. The performance analysis is also evaluated for real MR images with uncommon image details such as structural deformity (cavernoma is found in one subject from dataset-II), and another subject with Alzheimer’s disease from dataset-III which is acquired for the Alzheimer’s Disease Neuroimaging Initiative [53]. The selected MR images in dataset-III are acquired from a 3T MRI scanner with an MPRAGE sequence with a spatial resolution of 196 × 256 × 256 .
In the proposed method, an embedded framework for SR and denoising of MR images in an unsupervised framework is addressed. The proposed denoising method is built upon an LRA-based unsupervised framework for image denoising [15,29,30]. The novelty of the proposed denoising work lies in the method addressed to categorize the patch as smooth or textured, without the prior/estimated knowledge of noise variance. Further, the proposed idea to differently denoise the categorized patches leads to efficient denoising, and the re-estimation of LRA leads to better representation of image patches. This aids in improving the trade-off between noise suppression and smearing of edges. The proposed SR approach is developed upon the existing sparse representation-based framework for MR image patches [6,12]. The GPS values of brain MR images convey important semantic information such as tissue boundaries and should be preserved well and close to the gradient profiles in HR brain MR images. In the case of the proposed SR approach, the novelty lies in formulating a piecewise relation between GPS values of LR and upsampled LR images to estimate GPS values of unknown HR images. The estimated GPS values are used to regularize the estimation of HR images.
The initial HR image estimate is approximated by interpolating the noisy LR image which is denoised using the proposed denoising approach. Each iteration of the proposed SR approach is integrated with the proposed denoising approach, consisting of three consecutive steps in each iteration: (i) estimate the HR image by estimating its high-frequency component in a sparse representation framework and estimate its low-frequency component using a non-local means approach, (ii) apply the GPS-based regularization to enhance the image edge details, i.e., tissue boundaries, and (iii) the regularized image is denoised using the proposed denoising method to suppress the leftover noise. The SR approach provides better restoration of fine image details and the denoised image patches possess improved SNR, and are thus averaged for the next iteration of the SR algorithm. These three steps are followed by back projection and are repeated until the changes in the estimated HR image are insignificant. The flow of these steps in the proposed SR approach for noisy images and detailed steps for both SR and denoising are illustrated in Figure 1. The details of the proposed denoising and SR approaches are given in Section 3.1 and Section 3.2, respectively.

3.1. Proposed Denoising Method

The MR brain image consists of different kinds of patches such as patches with edges/texture depicting tissue boundaries, referred to as signal enriched (SE) patches (high signal related information) and patches with the relatively smooth region due to white matter and gray matter, referred to as noise enriched (NE) patches in the case of presence of noise.
It has been conjectured that noise belongs to the high-frequency range of the image spectrum, so the high-frequency (HF) component of an image patch is more contaminated with noise as compared to its low-frequency (LF) component. Thus, the HF component ( x h ) of a noisy image patch is processed for denoising and the respective LF component is added later in the estimated denoised HF component of image patch. The decomposition of LF ( x l ) and HF ( x h ) components of an image patch ( x ) is performed as x = x h + x l , where x l is the vector with each element as μ = i = 1 m x [ i ] m . Here, x = [ x [ i ] ] , i = 1 to m , and m represents the dimensionality of a patch vector.
In the proposed denoising method, the HF components of patches are extracted and grouped into clusters using the Euclidean distance-based k-means algorithm (for simplicity, we will omit the notion of an HF component of the patch from here onward). The grouping of patches helps in improving the homogeneity among patches in a cluster. However, the clustering may not be accurate due to the presence of noise, thus leading to certain inhomogeneity left in the clusters. We categorize the image patches in each cluster further into smooth/NE and textured/edge/SE patches using the LRA-based reconstruction error of a patch. The LRA of a patch x h is obtained as x ^ h = i = 1 l c i u i , l < m . The LRA-based reconstruction error for the jth patch is represented by e j and computed as
e j = | | x h j x ^ h j | | 2 2 i = l + 1 m | | c i | | 2
where x h j represents the HF component of the jth patch, x ^ h j represents the LRA of the jth image patch, and the respective projection coefficient of x h j on ith eigenvector is represented by c i . The ith unit norm eigenvector is represented by u i , after arranging eigenvectors according to descending order of eigenvalues. The inequality in Equation (1) holds true in general by keeping the eigenvectors of unit norm [54]. The LRA-based reconstruction error of a patch is upper bounded by its projection coefficients on the last eigenvectors. The different behavior of projection coefficients for different patch structures is used to categorize the patches. The details for categorization of the patch structure for each cluster and the denoising strategies for different patch structures and re-estimation of eigenvectors are given in Section 3.1.1 and Section 3.1.2.

3.1.1. Categorization of Smooth and Textured/Edge Patches

It has been observed that the smooth (NE) patches have smaller magnitude variations (due to noise) than the textured/edge (SE) patches (due to edges of varying strengths). Thus, the projection coefficients ( c i ) in smooth patches are constrained to small values in contrast to textured/edge patches. Following this and from Equation (1), small projection coefficients in smooth patches lead to less LRA-based reconstruction error ( e j ) for smooth patches than textured/edge patches, and thus can be used to categorize the patch as smooth or textured/edge.
This has been experimentally demonstrated in Figure 2 by computing the eigenvalues for textured/edge and smooth patches, separately. The patches extracted from a randomly selected noisy MR image (degraded with 2% noise) are categorized as smooth and textured/edge. A threshold is defined to divide the reconstruction error values into two ranges, each corresponding to either smooth or textured/edge patches and is discussed in the following subsection. We have considered patches with less reconstruction error ( e j ) as smooth patches and other patches as textured/edge patches, i.e., patches with varying strength of edges. The eigenvalues computed for smooth and textured/edge patches are shown using green and cyan color, respectively, in Figure 2. It can be observed that the textured/edge patches tend to obey the power law, i.e., the eigenvalues decay drastically from the first eigenvector to last eigenvector. However, the last eigenvalues of textured/edge patches are still higher than the eigenvalues of smooth patches, as shown in Figure 2, leading to higher reconstruction error for textured patches than smooth patches (followed from Equation (1)).
Estimation of threshold:
It shall be noted that the choice of threshold for reconstruction error of a patch ( e j ) to categorize the patch as smooth or textured is crucial, and is challenging in the presence of noise. Since LRA-based reconstruction error of a patch changes with the strength of noise, the threshold should also be adaptive to the strength of noise.
Adaptation of threshold to the strength of noise in image patches: The LRA-based reconstruction error for each patch (represented by e j for the jth patch) is concatenated to form reconstruction error vector for ensemble of patches as e = [ e 1 e 2 e 3 e N ] . The LRA-based reconstruction error for an ensemble of patches, i.e., e , is used to characterize the strength of noise present in image patches. For fixed l, an increase in noise variance leads to an increase in each eigenvalue and thus an increase in the reconstruction error of each patch. To demonstrate this, a randomly selected MR image is synthetically corrupted with Rician noise with different noise variance values and the histogram of reconstruction error values of all image patches ( e ) for different noise variance values is shown in Figure 3a. It can be observed that the histogram of reconstruction error values shifts right with increase in noise variance. The same has been observed for image corrupted with the Gaussian noise. Hence, a statistical measure of reconstruction error values obtained for image patches like mean of e computed as κ = j = 1 N e j N , can be used to correlate with strength of noise. Here, N represents the number of patches. The experimental demonstration of relation between mean of reconstruction error values ( κ ) and noise variance is shown in Figure 3b, and it can be observed that κ follows the linear relation with noise variance.
Hence, the adaptation of the threshold to strength of noise can be taken care of by choosing the threshold which is linearly related with κ , mean of reconstruction error (from Figure 3b). Following this, the categorization of an image patch ( x h j ) is done as
patch kind = textured / edges , i f τ 1 < e j < τ 2 smooth , o t h e r w i s e .
Here, τ 1 = ζ 1 j = 1 n e j n , and τ 2 = ζ 2 j = 1 n e j n such that ζ 2 > ζ 1 . The values for ζ 1 and ζ 2 are empirically chosen, and are described in the experimental section. The e j = | | x h j x ^ h j | | 2 2 represents reconstruction error of the jth patch, x h j represents the HF component of the jth patch and n represents the number of patches in the selected cluster. Here, τ 1 and τ 1 represent the noise adaptive thresholds used to categorize the image patch structure. The τ 2 represents the threshold for patches having higher projection coefficients only along the last few eigen directions. Such patches are usually rare but should also be denoised in the same manner as smooth patches. This is because the reconstruction of a patch using only the last few eigen directions conveys very little semantic meaning about image details and is usually considered to possess noise like characteristics.
The reason for employing Equation (2) in the categorization of image patches is experimentally illustrated in Figure 4. It shows the histogram of reconstruction error values obtained using LRA of smooth and textured patches of a randomly selected MR image in Figure 4a. The reconstruction error values (the x-axis of Figure 4a) are split into low and high values by using a threshold ( τ 1 ), defined using the mean of the reconstruction error of the ensemble of patches and considering ζ 1 = 0.7 , i.e., τ 1 = 0.7 j = 1 n e j n . The patches corresponding to low ( < τ 1 ) reconstruction error values are shown in the red box, and patches corresponding to high ( τ 1 ) reconstruction error values are shown in the green box. These patches are placed on their corresponding locations in 2D to form an image as shown in Figure 4b,c, respectively. It can be observed that the patches corresponding to low reconstruction error values are mostly smooth patches corrupted with noise, and high reconstruction error values correspond to patches dominated with texture/edges.
Different denoising strategies for different kinds of patches: The LRA for both smooth and textured/edge patches is computed using the considerate number of eigenvectors such that the edges are present in textured/edge patches even though the noise in these patches is not suppressed completely. The denoised version of textured/edge patches is represented by their estimated LRA, i.e., x ^ h . The pixels in each categorized smooth patch, are replaced by the average of pixel values in that patch. The averaging acts as a low pass filter and thus denoises the smooth patches. As these are smooth patches, averaging does not lead to any significant loss of image details, thus leading to efficient denoising.

3.1.2. Re-Estimation of LRA in an Iterative Manner

Eigenvectors play an essential role in estimation of textured/edge patches. Since the estimation of eigenvectors is done from noisy image patches, the approximation of a low-rank structure possessed by textured/edge patches may be noisy. To improve the reliability of the low-rank structure estimated by eigenvectors and for better reconstruction of image details, the proposed denoising procedure is repeated in an iterative manner. The eigenvectors are estimated in each iteration from the denoised patches obtained in the previous iteration. The patches, categorized as smooth patches in the previous iteration and denoised by an averaging operator, lie near the origin, i.e., are with insignificant values after mean removal (for extracting HF component). Thus, these patches do not play any role in the computation of eigenvectors in the next iteration. However, the averaging of categorized smooth patches in the previous iteration affects the overall distribution of patches, which leads to a change in estimation of eigenvectors in the next iteration such that the effect of noise on the computation of eigenvectors is reduced. Thus, repeating this process iteratively helps in increasing the reliability of estimated eigenvectors which capture the low-rank structure of noiseless patches, thus aiding in better reconstruction of image details. The summarized flow of steps followed in the proposed denoising approach is shown in Figure 5.

3.2. Super Resolution of Noisy MR Images

The above-mentioned denoising procedure of image patches is performed in each iteration of the SR algorithm (refer to Figure 1), which is based on the framework of sparse representation of image patches. In the SR approach, the LF components of HR image patches are reconstructed using a non-local means approach. On the contrary, the HF components of HR image patches are estimated using the framework of sparse representation framework. The unavailability of example HR MR images for forming the dictionary for sparse representation is addressed using the concept of self-similarity. The image is up and down scaled to exploit the self-similarity among patches of the image by forming an image pyramid. The HF components of image patches and patches from the up-down scaled version of the same image are clustered using the Euclidean distance-based k-means algorithm. The patches of an image and the up-down scaled version of the image aid in enrichment of dictionary construction. The k-means approach is a typical unsupervised method for clustering, and clusters the input image patches without any requirement of example data [54]. For each cluster, the eigenvectors are estimated from the covariance matrix of an ensemble of patches in that cluster and are concatenated column wise in the matrix, named the PCA-based dictionary (for more details about dictionary building, please refer to [33]). Please note that the estimated HR image is denoised using the proposed denoising approach in each iteration of the SR approach but the steps carried out in both are independent of each other.
The mathematical model for degradation of HR MR images for super resolution is adapted using a point spread function ( H ) (PSF, generally considered Gaussian kernel) and a downsampling operator ( B ) such as y = BHx , similar to the adaption of SR for MR images in [2,12]. Here, y R p × 1 represents the noisy LR patch vector obtained by blurring and downsampling the HR patch vector x R q × 1 , q > p using H R q × q and B R p × q , respectively. Thus, the objective function for the sparse representation-based SR algorithm with the constraint for preserving gradient profiles sharpness, proposed in this work, is as follows:
m i n α , η | | y BHA α | | 2 2 + | | α | | 1 + λ | | E x ^ { η } E x { η H } | | 2 2 .
The interpolated version of denoised y using the proposed denoising approach is considered as the initial estimate of x . The sparse representation of an image patch is computed as x ^ = A α where A is the PCA-based compact dictionary estimated for a cluster to which the initial estimate of x belongs and α denotes the coefficient vector. Since the sparsity is implicitly induced by thresholding the insignificant projection coefficients ( α ), the weight for term | | α | | 1 is not mentioned. Here, λ represents the weight for the regularizer. The η and η h represent the images with GPS values for the estimated HR image (its patches are represented by x ^ ) and the original HR image (its patches are represented by x ), respectively. Here, the GPS values corresponding to patch x ^ are extracted from η from the same location as patch x ^ using operator E x ^ { η } . Similarly, the E x { . } operator extracts the corresponding GPS values for patch x from η ^ H . Since the reference GPS values η H are unavailable, the constraint becomes | | E x ^ { η } E x { η ^ H } | | 2 2 , where η ^ H represents the approximated η H . The computation of GPS values and the estimation of η ^ H are explained below.
The GPS value indicates the ratio of height (h) and width (w) of the edge centered on an edge pixel in image, i.e., h / w as defined in [51]. Higher values of h indicate higher contrast and high values of w indicate more blur in the corresponding edge. Thus, higher GPS values correspond to sharper image details. The GPS value for each pixel in an image is computed to obtain GPS valued image η . It has been observed in the case of natural images [51] that the GPS values of an HR image ( η H ) possess a linear relation with GPS values in an upsampled LR image (UR) ( η U ) as η H = β η U . However, we observed that the relations among GPS values of HR and UR MR images are not exactly linear especially when the super-resolution factor is increased (for more details, refer to [33]). It can be due to the fact that there exist many small details in MR images which are degraded severely, and due to image details with varying sharpness in MR images which are degraded differently while upscaling. Thus, there exists a non-linearity between GPS values of UR and HR MR images. This non-linear behavior is approximated with a piecewise linear relation by dividing the range of η U values into four regions. For each ith region of η U values, GPS values of HR can be related as
η H i = β i η U i , here i = 1 to 4 .
The η H i and η U i represent the GPS values for the ith region of η U values. Here, β i , i = 1 to 4 , is unknown and is estimated to approximate η H i . It has been observed that most of the GPS values for a UR image lie in the small valued region, whereas there exist higher GPS values for both LR and HR images. Such behavior of the UR image can be due to the basic interpolation used for upsampling LR images that lead to blurred details and hence lead to small GPS values. Further, the percentage of edge pixels in LR and HR images is also generally considered as similar. Following which, the histogram of GPS values of LR images ( h i s t ( η L ) ) is assumed to be similar to h i s t ( η H ) . Hence, from Equation (4), we can estimate β i by minimizing the Chi-squared distance between h i s t ( η L i ) and h i s t ( β i η U i ) , for each region.
We thus estimate the GPS values for an HR image as
η ^ H i = β i η U i , i = 1 to 4 .
using LR and UR images only. This is used to constrain the solution space for the estimate of the HR image ( I ) by using η H η ^ H in Equation (3). Such a piecewise linear relation-based approximation for GPS values of an HR image aids in the improvement of the linear correlation of η H and η U for each region of η U , thus leading to a better approximation of the underlying relation. Further, it also helps in considering the low GPS values (obtained from relatively smoother regions) and high GPS values (from regions with edges) differently, which aids better preservation of tissue boundaries. It has to be noted that the four regions for a piecewise linear relation are empirically chosen as choosing three or five regions provides lesser overall linear correlation among η H and η U values. The step wise flow of the proposed approach is shown in Algorithm 1.
Algorithm 1 Estimate denoised HR image.
Input: Noisy LR image patches y
Output: Denoised HR image patches x ^
Variables: M: Number of iterations
k = 1 to M
Step 1: Denoise y as discussed in Section 3.1. Interpolate the denoised LR image patches to provide initial denoised HR patches x , represented by x ^ .
Step 2: Estimate GPS values η ^ H of x ^ as explained in Equation (5).
Step 3: Apply a low pass filter to x ^ to get the LF component and subtract it from x ^ to get the HF component of x ^ .
Step 4: Represent the LF component of each patch of x ^ using a non-local mean approach.
Step 5: Apply k-means clustering to HF components of x ^ , and up-down scaled versions of image patches x ^ . Construct a PCA-based dictionary for each cluster.
Step 6: Represent the HF component of x ^ using one of the dictionaries, i.e., the dictionary constructed from a cluster to which the HF component of x ^ belongs.
Step 7: Update the estimate of the HR patch x ^ by adding estimates of the LF (from Step 4) and HF (from Step 6) components.
Step 8: Regularize and update the estimated HR image patch x ^ obtained in Step 7 using GPS constraint.
Step 9: Denoise the estimate of HR patch x ^ obtained in Step 8 as discussed in Section 3.1.
If k > 0 and I ^ k I ^ k 1 < ϵ
      break;
else
Update x ^ as average of outputs from Step 8 and Step 9, i.e., I ^ k , in Step 3 and repeat the steps from Step 3 to Step 9.
end If
end while loop

4. Results and Discussion

The performance of the proposed algorithm is evaluated for SR of noisy MR images as well as individually for the denoising approach and super resolution along three directions, i.e., in-plane and slice-select direction. The reconstruction quality of the resultant images is reported using peak signal to noise ratio (PSNR), structural similarity index (SSIM) [55], and feature similarity index metric (FSIM) [56] values. Higher PSNR and SSIM values indicate the overall quality improvement of the reconstructed image. Higher FSIM indicates the higher subjective evaluation based on image gradient magnitude and phase congruency, and is generally used to indicate better features related to edge information [56]. The weight for GPS regularizer ( λ ) is empirically chosen as 0.001 for the best reconstruction of image patches and ζ 1 and ζ 2 are empirically chosen as 0.7 and 1.8, respectively. The patch size used in the experiments is 5 × 5 . The number of iterations for the denoising algorithm is fixed to 3 because it has been observed that further increase in iteration does not improve the performance. The maximum number of iterations in the proposed SR approach is fixed at 320 and is otherwise stopped early if the change in the estimated HR image is insignificant.

4.1. Denoising of MR Images

The experiments are performed for images corrupted with Rician and Gaussian distributed noise. The noise is synthetically added in simulated images as well as real MR images to simulate the real scenario. Here, the variance of noise is decided as a percentage of brightest tissue in an MR image, followed from the literature [34,41,57]. For example, 2% noise here means the noise variance is 0.02 × 255 = 5.1 for the image with 255 as its maximum intensity value. The proposed denoising approach is compared with conventional as well as state-of-the-art denoising methods addressed for MR images. These methods require prior knowledge of noise variance. In addition to noise variance, few approaches also require prior knowledge of noise distribution. The existing methods that process 2D images include non local means (NLM), universal NLM (UNLM) [34] and variance stabilization transform (VST) [25] with block-matching and 3D filtering (BM3D) [27]. Other existing methods denoise 3D image volume by considering the adjacent slices information. These methods include optimized blockwise NLM denoising filter optimized blockwise non local means (ONLM) [57] and its adapted version ONLM (AONLM) [40], multi-resolution-based ONLM multi-resolution based ONLM (MRONLM) [41], oracle-based discrete cosine transform filter oracle based discrete cosine transform (ODCT) explained in [57], PCA-based denoising over ODCT method PCA-based denoising over ODCT (PRINLM) [35], and VST in conjunction with block matching with 4D filtering (BM4D) and AONLM [25].
One real MR image is randomly selected, degraded with 2% Gaussian noise, and is denoised using different algorithms. The denoised images are shown in Figure 6. It can be observed that the proposed approach has reduced smearing of image details (in the yellow box), provides efficient denoising in the smooth region, and estimates edges with comparable quality to existing methods in the cyan box in Figure 6.
Table 1 gives the quantitative comparison of the proposed approach with existing approaches. Here, the highest PSNR and SSIM values are highlighted in red color. PSNR and SSIM values obtained for proposed work are made bold. The blue colored values indicate the performance of existing algorithms that are comparable to the proposed work. The estimation of noise variance is crucial and is observed to be inaccurate as the level of noise increases, leading to decreased performance of the existing algorithms [58]. Hence, for a fair comparison, the existing works with the requirement of prior knowledge of noise variance are given the true variance value in the experiments.
It can be observed that the proposed work, without the requirement of knowledge of variance and distribution of noise, performs comparable to the existing classical methods such as NLM and UNLM [34], and the recent methods like ODCT and PRINLM [35] that, on contrary, have such requirements. The proposed method provides PSNR values within the 0.5 dB range of these existing approaches. This can be explained by the fact that in NLM and UNLM approaches, the images are denoised using weighted averaging and thus several image details get blurred, leading to its comparable performance to the proposed approach.
The existing supervised 3D volume-based approaches addressing denoising of MR images perform better than the proposed approach. These approaches assume the prior distribution of noise and such assumption is the reason that VST-based methods, i.e., VST+BM4D and VST+AONLM, perform better for images with Rician distributed noise than with Gaussian distributed noise. These approaches also require prior knowledge of noise variance. Hence, such prior knowledge aids in better performance for these approaches as compared to the proposed approach, but is not generally available in the real scenario.

4.1.1. Significance of the Proposed Patch Categorization and Re-estimation of Eigenvectors

The significance of the proposed patch categorization is demonstrated by incorporating it in an existing LRA-based denoising method [15]. For this, the image patches are categorized as smooth and textured/edge patches using the proposed method. The smooth patches are denoised using the proposed method and the textured/edge patches are denoised using the existing method [15]. The eigenvectors in the existing method are re-estimated in the progressive manner as mentioned in the proposed denoising approach. A randomly selected denoised image (3% Rician noise) obtained by this procedure is shown in Figure 7d. The noisy input and noise-free images are shown in Figure 7a,e, respectively. The images denoised using the method described in [15] is shown in Figure 7b. It can be observed that smooth patches are not efficiently denoised in the existing LRA-based approach [15]. The denoising of such smooth patches is improvised by adapting the existing approach to the proposed denoising method indicating the advantages of patch categorization. On the other hand, the tissue boundaries represented by the proposed denoising approach are not denoised well but improvise by using patch adaptive l for textured/edge patches from the existing approach [15]. This can be due to the empirically chosen fixed l in the proposed denoising approach for representation textured/edge patches.

4.1.2. Parameter Analysis for Optimal Selection of τ 1 in the Proposed Denoising Approach

In the proposed denoising approach, τ 1 is a crucial parameter as it is challenging to separate the smooth and textured/edge patches in the presence of noise. In addition to the noise adaptive threshold τ 1 , the optimal threshold has to be chosen such that no texture/edge patch is categorized as smooth patch because it may lead to loss of image details. To analyze the change in denoising with the change in τ 1 = ζ 1 j = 1 n e j n , a randomly selected noisy MR image (2% noise) is denoised with varying values of ζ 1 . The obtained PSNR values for denoised images are plotted against different values for ζ 1 in Figure 8. It can be observed that there is a slight increase in PSNR values for smaller values of ζ 1 . The increase in PSNR value at ζ 1 = 0.7 compared to the image denoised with ζ 1 = 0.25 can be explained by further reduction in noise at ζ 1 = 0.7 , as shown in Figure 8. The PSNR value is reduced significantly as the value for ζ 1 is further increased. This can be explained by the corresponding image at ζ 1 = 1 , i.e., τ 1 = 1 j = 1 n e j n , where even the textured/edge patches are considered as smooth and are denoised accordingly. This lead to huge image detail loss and thus the PSNR value drops. Hence, ζ 1 = 0.7 is used for the experiments in this work. Similarly, the value of ζ 2 is used to govern the minimum loss of image details while suppressing the noise and it has been observed that the value of 1.8 gives the best results.

4.2. Super Resolution of MR Image Volumes

The LR image volumes are simulated by blurring the HR image volumes using a Gaussian kernel of the standard deviation of one voxel size and followed by downsampling. One MR image volume is randomly selected from dataset-II and super-resolved with super-resolution factor (SRF) 2 using the proposed SR approach and is compared with the existing state-of-the-art unsupervised methods NLM3D [13] and LRTV [12]. The reconstructed images using these approaches are shown in Figure 9. The highlighted region in the coronal slice is zoomed and shown in the red box. The pink and green arrows show the location for tissue boundary between cerebrospinal fluid (CSF) and gray matter, and white matter and gray matter, respectively. It can be observed that the proposed approach reconstructs the image details with relatively sharper tissue boundaries as well as better preservation of image details (near green arrow) as compared to the existing methods. Table 2 shows the mean PSNR, SSIM and FSIM values computed for different subjects for SRFs 2 and 3 reconstructed by different methods. The graphical representation of Table 2 is shown in Figure 10. It can be observed that the proposed approach provides comparably higher PSNR, SSIM, and FSIM values than existing methods, and thus illustrates the advantages of the proposed work in pixel-intensity, and structural and gradient-feature-based similarities, respectively.
A popular supervised approach in neural network super-resolution convolutional neural network (SR-CNN) [9] (after training using paired example LR-HR images) is also tested for SR of the MR image volumes used for evaluation of performance of proposed approach. It provides the maximum SSIM value 0.9012 and the average SSIM value reaches 0.8666 for (SRF = 2) which is comparable to the proposed approach.
To emphasize the significance of reconstruction of uncommon image details in the proposed SR method, we chose a real MR image volume with a cavernoma and improved its resolution by factor 2, as shown in Figure 11. It can be observed that the LRTV [12] approach in Figure 11d reduces the blur as compared to interpolation and NLM3D [13] in Figure 11b,c, but tends to produce staircase effect. It can be observed that the tissue boundaries and image details inside the cavernoma, shown in the red box, are more clear and distinct for the proposed method. In addition, the skull outlines can be seen to be well defined (in the sagittal plane) in the case of the proposed algorithm (see Figure 11e) as compared to existing methods.

Parameter Analysis for Optimal Selection of λ in the Proposed Super Resolution Approach

In the proposed approach, the weight to GPS-based regularizer is denoted by λ . It is a crucial parameter that decides the balance between the data fidelity term and the regularizer for gradient profiles. The optimal value of λ is chosen empirically based on the maximization of PSNR for estimated HR images. For this, an example MR image is upscaled by factor 2 with different values of λ , and PSNR values are plotted in Figure 12. It can be observed that with the high value of λ , i.e., beyond 0.005, the PSNR values dropped. Similarly for lesser values of λ , the PSNR value is less. The optimal of λ is thus chosen as 0.001 for which the PSNR value obtained is maximum, and holds true for another ensemble of MR images.

4.3. Resolution Improvement of Noisy MR Image Volumes

The performance of the proposed single framework for the super resolution of noisy MR images has been evaluated here for 2% Gaussian noise and different upsampling factors. For existing SR methods, the images are first denoised using AONLM [40]. The qualitative results using different algorithms for upsampling a randomly selected noisy image by factor 2 are shown in Figure 13. To illustrate the significance of the proposed approach for efficient denoising as well as the preservation of image details, different kinds of regions are zoomed and shown separately. The green box in Figure 13 comprises the region with strong edges, the red box shows the region with minor edges and texture, and the cyan box shows the smooth region. It can be observed that the reconstructed image using spline interpolation and the non local means in three dimensions (NLM3D) approach in Figure 13b,c blur the image details in every kind of region and hence provide poorer contrast but with better denoising in each region. The LRTV method, however, tends to provide better contrast among tissues and preservation of image details, yet it is not able to provide efficiently denoised smooth regions (see red and cyan boxes in Figure 13d). The proposed work, on the contrary, can be observed to reconstruct strong edges while preserving the image details with limited noise (see green box Figure 13e), it can relatively better suppress the noise in patches with smooth and texture/minor edges as compared to LRTV [12] (can be seen in red and cyan boxes), and it can provide improved contrast as compared to NLM3D [13]. It demonstrates that the proposed method tends to improve the trade-off between the preservation of image details and suppression of noise. The obtained PSNR, SSIM, and FSIM values for different subjects are summarized in Table 3 and show that the proposed work performs relatively better than existing unsupervised methods when combined with supervised denoising method.
To see the effect of different noise levels, a volume is randomly selected followed by degradation with downsampling factor 2 and varying noise levels (1% to 4%). The obtained PSNR values using different algorithms are mentioned in Table 4. It can be observed that the increase in PSNR values of images reconstructed by the proposed work are decreased as the noise level is increased and this is expected because restoring the image details for a high level of noise is a challenging problem. In addition, the PSNR metric represents an average measure of similarity and thus can be higher even for the blurred image as compared to the image with left out noise. This is observed to be true in case of 4% noise, where NLM3D [13] provides the blurred image with comparable PSNR. The FSIM and SSIM values are obtained for the estimated HR image using the proposed approach for varying variance of noise and upscaling factors, and are plotted in Figure 14. The metrics are obtained for one randomly chosen real MR image volume from dataset-II. It can be observed that the FSIM values decrease at a lesser rate when strength of noise is increased as compared to the case where the upscale factor is increased.

Analysis for Alzheimer Subjects

The clinical applications of the proposed approach can be valid only if no artifacts are introduced while denoising or super resolving. To experimentally verify this, a noisy LR image of an Alzheimer’s subject is restored by different algorithms and compared in the first row of Figure 15. The zoomed version of the region in the blue rectangle (highlighted in the first row) is shown in the second row. It can be observed that interpolated images and NLM3D [13] provide blurred image details. The LRTV [12] approach highlights the CSF flow artifact significantly (which is not as present in the noise-free HR image), may compromise the detection of pathology in the ventricular system, thus compromising the diagnostic capability. The image reconstructed by the proposed approach looks relatively noisy but clearer than existing methods and does not highlight any such artifact; hence, it can aid in accurate diagnosis. The same values of ζ 1 = 0.7 , ζ 2 = 1.8 , and λ = 0.001 were used to process the MR images of an Alzheimer’s subject obtained from dataset-III as for processing dataset-II and dataset-I. It indicates the generalization of the proposed approach that uses the same values of parameters for different datasets.

5. Summary

This paper addresses the embedded framework for the super resolution of the noisy MR images. The smooth and textured/edge patches are categorized, without the prior/estimated knowledge of noise, and are denoised using different denoising strategies for efficient denoising. Further, the degradation of gradient profiles is addressed specifically for MR images to regularize the estimation of the HR image with clear tissue boundaries. The experimental section demonstrates the significance of proposed different denoising strategies for categorized patches to achieve efficient denoising. The importance of GPS constraint is demonstrated in the results section by clear tissue boundaries, i.e., reduced blur and reduced staircase effect, and higher FSIM values in estimated HR images. The experimental section demonstrates that the proposed unsupervised denoising approach performs comparably to the conventional [34] and a few recent supervised denoising approaches [25,35,40], which require prior knowledge of noise. The proposed integrated framework performs better than the existing unsupervised approaches [12,13] for super resolution of real MR images (with and without synthetically added noise). The proposed unsupervised SR method performs better than existing unsupervised methods for each kind of region (smooth/textured/strong edges) in noisy MR images, and thus improves the trade-off between denoising and smearing of image details. In addition, the proposed work performs comparable to the supervised deep learning-based SR approach [9] indicating the potential of unsupervised approaches.

Author Contributions

Conceptualization, P.K. and A.K.S.; methodology, P.K. and A.K.S.; validation, P.K., A.K.S. and C.K.A.; formal analysis, P.K.; investigation, P.K.; writing—original draft preparation, P.K.; writing—review and editing, P.K., A.K.S. and C.K.A.; visualization, P.K.; supervision, A.K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Minstry of Electronics and Information Technology (MietY) India under grant VISPHD-MEITY-750.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Plenge, E.; Poot, D.H.J.; Bernsen, M.; Kotek, G.; Houston, G.; Wielopolski, P.; van der Weerd, L.; Niessen, W.J.; Meijering, E. Super-resolution methods in MRI: Can they improve the trade-off between resolution, signal-to-noise ratio, and acquisition time? Magn. Reson. Med. 2012, 68, 1983–1993. [Google Scholar] [CrossRef]
  2. Van Reeth, E.; Tham, I.W.K.; Tan, C.H.; Poh, C.L. Super-resolution in magnetic resonance imaging: A review. Concepts Magn. Reson. Part A 2012, 40A, 306–325. [Google Scholar] [CrossRef]
  3. Kim, G.; Park, J.; Lee, K.; Lee, J.; Min, J.; Lee, B.; Han, D.K.; Ko, H. Unsupervised Real-World Super Resolution with Cycle Generative Adversarial Network and Domain Discriminator. In Proceedings of the 2020 IEEE/CVF Conference on Computer Vision and Pattern Recognition Workshops (CVPRW), Seattle, WA, USA, 14–19 June 2020; pp. 1862–1871. [Google Scholar] [CrossRef]
  4. Rudin, L.I.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D Nonlinear Phenom. 1992, 60, 259–268. [Google Scholar] [CrossRef]
  5. Shi, J.; Li, Z.; Ying, S.; Wang, C.; Liu, Q.; Zhang, Q.; Yan, P. MR Image Super-Resolution via Wide Residual Networks With Fixed Skip Connection. IEEE J. Biomed. Health Inform. 2019, 23, 1129–1140. [Google Scholar] [CrossRef]
  6. Rueda, A.; Malpica, N.; Romero, E. Single-image super-resolution of brain MR images using overcomplete dictionaries. Med. Image Anal. 2013, 17, 113–132. [Google Scholar] [CrossRef]
  7. Pham, C.H.; Ducournau, A.; Fablet, R.; Rousseau, F. Brain MRI super-resolution using deep 3D convolutional networks. In Proceedings of the IEEE International Symposium on Biomedical Imaging (ISBI), Melbourne, Australia, 18–21 April 2017; pp. 197–200. [Google Scholar]
  8. Chen, Y.; Xie, Y.; Zhou, Z.; Shi, F.; Christodoulou, A.G.; Li, D. Brain MRI super resolution using 3D deep densely connected neural networks. In Proceedings of the IEEE International Symposium on Biomedical Imaging (ISBI), Washington, DC, USA, 4–7 April 2018; pp. 739–742. [Google Scholar]
  9. Dong, C.; Loy, C.C.; He, K.; Tang, X. Image Super-Resolution Using Deep Convolutional Networks. IEEE Trans. Pattern Anal. Mach. Intell. 2016, 38, 295–307. [Google Scholar] [CrossRef] [Green Version]
  10. Yang, J.; Wright, J.; Huang, T.S.; Ma, Y. Image Super-Resolution Via Sparse Representation. IEEE Trans. Image Process. 2010, 19, 2861–2873. [Google Scholar] [CrossRef] [PubMed]
  11. Zhu, Q.; Ren, Y.; Qiu, Z.; Wang, W. Robust MR image super-resolution reconstruction with cross-modal edge-preserving regularization. Int. J. Imaging Syst. Technol. 2019, 29, 491–500. [Google Scholar] [CrossRef]
  12. Shi, F.; Cheng, J.; Wang, L.; Yap, P.T.; Shen, D. LRTV: MR Image Super-Resolution With Low-Rank and Total Variation Regularizations. IEEE Trans. Med. Imaging 2015, 34, 2459–2466. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Manjón, J.V.; Coupé, P.; Buades, A.; Fonov, V.; Louis Collins, D.; Robles, M. Non-local MRI upsampling. Med. Image Anal. 2010, 14, 784–792. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, L.; Dong, W.; Zhang, D.; Shi, G. Two-stage image denoising by principal component analysis with local pixel grouping. Pattern Recognit. 2010, 43, 1531–1549. [Google Scholar] [CrossRef] [Green Version]
  15. Mandal, S.; Bhavsar, A.; Sao, A.K. Noise adaptive super-resolution from single image via non-local mean and sparse representation. Signal Process. 2017, 132, 134–149. [Google Scholar] [CrossRef]
  16. Manjón, J.V.; Coupé, P.; Buades, A. MRI noise estimation and denoising using non-local PCA. Med. Image Anal. 2015, 22, 35–47. [Google Scholar] [CrossRef]
  17. He, Y.; Gan, T.; Chen, W.; Wang, H. Adaptive Denoising by Singular Value Decomposition. IEEE Signal Process. Lett. 2011, 18, 215–218. [Google Scholar] [CrossRef]
  18. Qu, L.; Wang, S.; Yap, P.T.; Shen, D. Wavelet-based Semi-supervised Adversarial Learning for Synthesizing Realistic 7T from 3T MRI. In Medical Image Computing and Computer Assisted Intervention—MICCAI 2019; Shen, D., Liu, T., Peters, T.M., Staib, L.H., Essert, C., Zhou, S., Yap, P.T., Khan, A., Eds.; Springer International Publishing: Cham, Switzerland, 2019; pp. 786–794. [Google Scholar]
  19. Xue, X.; Wang, Y.; Li, J.; Jiao, Z.; Ren, Z.; Gao, X. Progressive Sub-Band Residual-Learning Network for MR Image Super Resolution. IEEE J. Biomed. Health Inform. 2020, 24, 377–386. [Google Scholar] [CrossRef]
  20. Chaudhari, A.S.; Fang, Z.; Kogan, F.; Wood, J.; Stevens, K.J.; Gibbons, E.K.; Lee, J.H.; Gold, G.E.; Hargreaves, B.A. Super-resolution musculoskeletal MRI using deep learning. Magn. Reson. Med. 2018, 80, 2139–2154. [Google Scholar] [CrossRef]
  21. Trinh, D.; Luong, M.; Dibos, F.; Rocchisani, J.; Pham, C.; Nguyen, T.Q. Novel Example-Based Method for Super-Resolution and Denoising of Medical Images. IEEE Trans. Image Process. 2014, 23, 1882–1895. [Google Scholar] [CrossRef]
  22. Awate, S.P.; Whitaker, R.T. Feature-Preserving MRI Denoising: A Nonparametric Empirical Bayes Approach. IEEE Trans. Med. Imaging 2007, 26, 1242–1255. [Google Scholar] [CrossRef] [Green Version]
  23. Mumford, D.; Shah, J. Optimal approximations by piecewise smooth functions and associated variational problems. Commun. Pure Appl. Math. 1989, 42, 577–685. [Google Scholar] [CrossRef] [Green Version]
  24. Donoho, D.L.; Johnstone, I.M. Ideal spatial adaptation by wavelet shrinkage. Biometrika 1994, 81, 425–455. [Google Scholar] [CrossRef]
  25. Foi, A. Noise estimation and removal in MR imaging: The variance-stabilization approach. In Proceedings of the 2011 IEEE International Symposium on Biomedical Imaging: FromNano to Macro, Chicago, IL, USA, 30 March–2 April 2011; pp. 1809–1814. [Google Scholar]
  26. Ahmad, R.; Bouman, C.A.; Buzzard, G.T.; Chan, S.; Liu, S.; Reehorst, E.T.; Schniter, P. Plug-and-Play Methods for Magnetic Resonance Imaging: Using Denoisers for Image Recovery. IEEE Signal Process. Mag. 2020, 37, 105–116. [Google Scholar] [CrossRef] [Green Version]
  27. Dabov, K.; Foi, A.; Katkovnik, V.; Egiazarian, K. Image Denoising by Sparse 3-D Transform-Domain Collaborative Filtering. IEEE Trans. Image Process. 2007, 16, 2080–2095. [Google Scholar] [CrossRef]
  28. Parekh, A.; Selesnick, I.W. Enhanced Low-Rank Matrix Approximation. IEEE Signal Process. Lett. 2016, 23, 493–497. [Google Scholar] [CrossRef]
  29. Does, M.D.; Olesen, J.L.; Harkins, K.D.; Serradas-Duarte, T.; Gochberg, D.F.; Jespersen, S.N.; Shemesh, N. Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry. Magn. Reson. Med. 2019, 81, 3503–3514. [Google Scholar] [CrossRef] [PubMed]
  30. Fu, Y.; Dong, W. 3D magnetic resonance image denoising using low-rank tensor approximation. Neurocomputing 2016, 195, 30–39. [Google Scholar] [CrossRef]
  31. Nguyen, H.M.; Peng, X.; Do, M.N.; Liang, Z. Denoising MR Spectroscopic Imaging Data With Low-Rank Approximations. IEEE Trans. Biomed. Eng. 2013, 60, 78–89. [Google Scholar] [CrossRef] [Green Version]
  32. Wei, S.; Xu, S. Staircasing reduction model applied to total variation based image reconstruction. In Proceedings of the 2009 17th European Signal Processing Conference, Glasgow, UK, 24–28 August 2009; pp. 2579–2583. [Google Scholar]
  33. Kaur, P.; Sao, A. Gradient Profile Based Super Resolution of MR Images. In Proceedings of the Medical Image Computing and Computer Assisted Intervention, Granada, Spain, 16–20 September 2018; pp. 109–117. [Google Scholar]
  34. Manjón, J.V.; Carbonell-Caballero, J.; Lull, J.J.; García-Martí, G.; Martí-Bonmatí, L.; Robles, M. MRI denoising using Non-Local Means. Med. Image Anal. 2008, 12, 514–523. [Google Scholar] [CrossRef]
  35. Manjón, J.V.; Coupé, P.; Buades, A.; Collins, D.L.; Robles, M. New methods for MRI denoising based on sparseness and self-similarity. Med. Image Anal. 2012, 16, 18–27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lee, O.Y.; Lee, J.W.; Kim, J.O. Combining self-learning based super-resolution with denoising for noisy images. J. Vis. Commun. Image Represent. 2017, 48, 66–76. [Google Scholar] [CrossRef]
  37. Hu, J.; Wu, X.; Zhou, J. Noise robust single image super-resolution using a multiscale image pyramid. Signal Process. 2018, 148, 157–171. [Google Scholar] [CrossRef]
  38. Gudbjartsson, H.; Patz, S. The rician distribution of noisy MRI data. Magn. Reson. Med. 1995, 34, 910–914. [Google Scholar] [CrossRef]
  39. Nowak, R.D. Wavelet-based Rician noise removal for magnetic resonance imaging. IEEE Trans. Image Process. 1999, 8, 1408–1419. [Google Scholar] [CrossRef] [Green Version]
  40. Manjón, J.V.; Coupé, P.; Martí-Bonmatí, L.; Collins, D.L.; Robles, M. Adaptive non-local means denoising of MR images with spatially varying noise levels. J. Magn. Reson. Imaging 2010, 31, 192–203. [Google Scholar] [CrossRef] [Green Version]
  41. Coupe, P.; Manjon, J.V.; Robles, M.; Collins, D.L. Adaptive multiresolution non-local means filter for three-dimensional magnetic resonance image denoising. IET Image Process. 2012, 6, 558–568. [Google Scholar] [CrossRef]
  42. Liu, R.W.; Shi, L.; Huang, W.; Xu, J.; Yu, S.C.H.; Wang, D. Generalized total variation-based MRI Rician denoising model with spatially adaptive regularization parameters. Magn. Reson. Imaging 2014, 32, 702–720. [Google Scholar] [CrossRef]
  43. Peled, S.; Yeshurun, Y. Superresolution in MRI: Application to human white matter fiber tract visualization by diffusion tensor imaging. Magn. Reson. Med. 2001, 45, 29–35. [Google Scholar] [CrossRef]
  44. Shilling, R.Z.; Robbie, T.Q.; Bailloeul, T.; Mewes, K.; Mersereau, R.M.; Brummer, M.E. A Super-Resolution Framework for 3-D High-Resolution and High-Contrast Imaging Using 2-D Multislice MRI. IEEE Trans. Med. Imaging 2009, 28, 633–644. [Google Scholar] [CrossRef] [PubMed]
  45. Tamez-Pena, J.G.; Totterman, S.; Parker, K.J. MRI isotropic resolution reconstruction from two orthogonal scans. In Medical Imaging 2001: Image Processing, International Society for Optics and Photonics; SPIE: Bellingham, WA, USA, 2001; Volume 4322, pp. 87–97. [Google Scholar]
  46. Zhao, X.; Zhang, Y.; Zhang, T.; Zou, X. Channel Splitting Network for Single MR Image Super-Resolution. IEEE Trans. Image Process. 2019, 28, 5649–5662. [Google Scholar] [CrossRef] [Green Version]
  47. Cherukuri, V.; Guo, T.; Schiff, S.J.; Monga, V. Deep MR Brain Image Super-Resolution Using Spatio-Structural Priors. IEEE Trans. Image Process. 2020, 29, 1368–1383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Rousseau, F.; Initiative, T.A.D.N. A non-local approach for image super-resolution using intermodality priors. Med. Image Anal. 2010, 14, 594–605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Zheng, H.; Zeng, K.; Guo, D.; Ying, J.; Yang, Y.; Peng, X.; Huang, F.; Chen, Z.; Qu, X. Multi-Contrast Brain MRI Image Super-Resolution With Gradient-Guided Edge Enhancement. IEEE Access 2018, 6, 57856–57867. [Google Scholar] [CrossRef]
  50. Zhang, D.; He, J.; Zhao, Y.; Du, M. MR image super-resolution reconstruction using sparse representation, nonlocal similarity and sparse derivative prior. Comput. Biol. Med. 2015, 58, 130–145. [Google Scholar] [CrossRef] [PubMed]
  51. Yan, Q.; Xu, Y.; Yang, X.; Nguyen, T.Q. Single Image Superresolution Based on Gradient Profile Sharpness. IEEE Trans. Image Process. 2015, 24, 3187–3202. [Google Scholar] [PubMed]
  52. Guan, J.; Zhang, W.; Gu, J.; Ren, H. No-reference blur assessment based on edge modeling. J. Vis. Commun. Image Represent. 2015, 29, 1–7. [Google Scholar] [CrossRef]
  53. Jack, C.R., Jr.; Bernstein, M.A.; Fox, N.C.; Thompson, P.; Alexander, G.; Harvey, D.; Borowski, B.; Britson, P.J.; Whitwell, J.L.; Ward, C.; et al. The Alzheimer’s disease neuroimaging initiative (ADNI): MRI methods. J. Magn. Reson. Imaging 2008, 27, 685–691. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Bishop, C.M. Pattern Recognition and Machine Learning; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  55. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: from error visibility to structural similarity. IEEE Trans. Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [Green Version]
  56. Zhang, L.; Zhang, L.; Mou, X.; Zhang, D. FSIM: A Feature Similarity Index for Image Quality Assessment. IEEE Trans. Image Process. 2011, 20, 2378–2386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Coupe, P.; Yger, P.; Prima, S.; Hellier, P.; Kervrann, C.; Barillot, C. An Optimized Blockwise Nonlocal Means Denoising Filter for 3-D Magnetic Resonance Images. IEEE Trans. Med. Imaging 2008, 27, 425–441. [Google Scholar] [CrossRef] [Green Version]
  58. Coupé, P.; Manjón, J.V.; Gedamu, E.; Arnold, D.; Robles, M.; Collins, D.L. Robust Rician noise estimation for MR images. Med. Image Anal. 2010, 14, 483–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Detailed illustration of the proposed approach.
Figure 1. Detailed illustration of the proposed approach.
Jimaging 07 00101 g001
Figure 2. Demonstration of eigenvalues obtained from an ensemble of smooth patches (green) and textured/edge patches (cyan) present in a randomly selected magnetic resonance (MR) image degraded with 2% Rician distributed noise. The eigenvalues for patches with only noise are shown in red color. The x-axis denotes the number of eigenvector and the y-axis denotes the corresponding eigenvalue. The plot is zoomed for the 15th to the 25th eigenvectors and is shown to better visualize the small differences among eigenvalues. Here, 2% noise denotes that noise variance is 0.02 × 255 = 5.1 for an image with 255 as its maximum pixel intensity.
Figure 2. Demonstration of eigenvalues obtained from an ensemble of smooth patches (green) and textured/edge patches (cyan) present in a randomly selected magnetic resonance (MR) image degraded with 2% Rician distributed noise. The eigenvalues for patches with only noise are shown in red color. The x-axis denotes the number of eigenvector and the y-axis denotes the corresponding eigenvalue. The plot is zoomed for the 15th to the 25th eigenvectors and is shown to better visualize the small differences among eigenvalues. Here, 2% noise denotes that noise variance is 0.02 × 255 = 5.1 for an image with 255 as its maximum pixel intensity.
Jimaging 07 00101 g002
Figure 3. (a) Illustration of the relation between distribution of LRA-based reconstruction error values ( e ) of an ensemble of image patches and strength of noise, i.e., noise variance. (b) Illustration of the relation between mean of reconstruction error values with strength of noise, i.e., noise variance.
Figure 3. (a) Illustration of the relation between distribution of LRA-based reconstruction error values ( e ) of an ensemble of image patches and strength of noise, i.e., noise variance. (b) Illustration of the relation between mean of reconstruction error values with strength of noise, i.e., noise variance.
Jimaging 07 00101 g003
Figure 4. Demonstration of the relation between reconstruction error vector e and patch structure, i.e., smooth and textured/edge, using the noise adaptive threshold τ 1 , for (mean subtracted patches of) a randomly selected image corrupted with 2% Rician distributed noise: (a) Histogram of reconstruction error values e and the corresponding image, (b) histogram of reconstruction error values less than τ 1 and the corresponding patches in the image, and (c) histogram of reconstruction error values greater or equal to τ 1 and the corresponding patches in the image.
Figure 4. Demonstration of the relation between reconstruction error vector e and patch structure, i.e., smooth and textured/edge, using the noise adaptive threshold τ 1 , for (mean subtracted patches of) a randomly selected image corrupted with 2% Rician distributed noise: (a) Histogram of reconstruction error values e and the corresponding image, (b) histogram of reconstruction error values less than τ 1 and the corresponding patches in the image, and (c) histogram of reconstruction error values greater or equal to τ 1 and the corresponding patches in the image.
Jimaging 07 00101 g004
Figure 5. Illustration of major steps in the proposed denoising approach.
Figure 5. Illustration of major steps in the proposed denoising approach.
Jimaging 07 00101 g005
Figure 6. Demonstration of image denoising achieved by different denoising algorithms on real MR images with synthetically added 2% Gaussian noise. From left to right—noisy image, optimized blockwise non local means (ONLM) [57], multi-resolution based ONLM (MRONLM) [41], oracle based discrete cosine transform (ODCT), PCA-based denoising over ODCT (PRINLM) [35], variance stabilized transform with blockwise matching and 4D filtering (VST-BM4D) [25], adaptive ONLM (AONLM) [40], the proposed approach, and the original noiseless image.
Figure 6. Demonstration of image denoising achieved by different denoising algorithms on real MR images with synthetically added 2% Gaussian noise. From left to right—noisy image, optimized blockwise non local means (ONLM) [57], multi-resolution based ONLM (MRONLM) [41], oracle based discrete cosine transform (ODCT), PCA-based denoising over ODCT (PRINLM) [35], variance stabilized transform with blockwise matching and 4D filtering (VST-BM4D) [25], adaptive ONLM (AONLM) [40], the proposed approach, and the original noiseless image.
Jimaging 07 00101 g006
Figure 7. Illustration of the significance of the proposed categorization of smooth and textured/edge patches and progressive estimation of eigenvectors by adapting it with the existing low-rank approximation (LRA)-based denoised approach [15]. (a) Noisy image (3% Rician Noise), (b) denoised using the existing LRA-based denoising approach [15], (c) the proposed denoising approach, (d) the existing approach [15] adapted with the proposed patch categorization and re-estimation of eigenvectors, (e) noise-free image.
Figure 7. Illustration of the significance of the proposed categorization of smooth and textured/edge patches and progressive estimation of eigenvectors by adapting it with the existing low-rank approximation (LRA)-based denoised approach [15]. (a) Noisy image (3% Rician Noise), (b) denoised using the existing LRA-based denoising approach [15], (c) the proposed denoising approach, (d) the existing approach [15] adapted with the proposed patch categorization and re-estimation of eigenvectors, (e) noise-free image.
Jimaging 07 00101 g007
Figure 8. Analysis of parameter ζ 1 for optimal selection of τ 1 = ζ 1 i = 1 N e i , e i denotes the low rank-based reconstruction error of the ith patch and N denotes the number of patches, for categorization of the patch as smooth or textured/edge in the proposed denoising approach. The peak signal to noise ratio (PSNR) values are computed for an MR image denoised using different values of ζ 1 , and plotted in this Figure. The denoised images using ζ 1 as 0.25, 0.7 and 1.0 are also shown to indicate the reason for rise and drop in PSNR values.
Figure 8. Analysis of parameter ζ 1 for optimal selection of τ 1 = ζ 1 i = 1 N e i , e i denotes the low rank-based reconstruction error of the ith patch and N denotes the number of patches, for categorization of the patch as smooth or textured/edge in the proposed denoising approach. The peak signal to noise ratio (PSNR) values are computed for an MR image denoised using different values of ζ 1 , and plotted in this Figure. The denoised images using ζ 1 as 0.25, 0.7 and 1.0 are also shown to indicate the reason for rise and drop in PSNR values.
Jimaging 07 00101 g008
Figure 9. Illustration of reconstruction results for real MR image volume for super-resolution factor 2 using different super-resolution algorithms: (a) nearest neighbor interpolation, (b) spline interpolation, (c) non-local means [13], (d) low rank and total variation (LRTV) [12], (e) the proposed method, and (f) the original denoised high resolution (HR) image. Zoomed version of the red box shown in the axial slice is shown to demonstrate the difference (specifically in tissue boundaries indicated with arrows).
Figure 9. Illustration of reconstruction results for real MR image volume for super-resolution factor 2 using different super-resolution algorithms: (a) nearest neighbor interpolation, (b) spline interpolation, (c) non-local means [13], (d) low rank and total variation (LRTV) [12], (e) the proposed method, and (f) the original denoised high resolution (HR) image. Zoomed version of the red box shown in the axial slice is shown to demonstrate the difference (specifically in tissue boundaries indicated with arrows).
Jimaging 07 00101 g009
Figure 10. Demonstration of PSNR, structural similarity index (SSIM) and feature similarity index metric (FSIM) values obtained for super resolution of eight MR image volumes by upscale factors 2 and 3.
Figure 10. Demonstration of PSNR, structural similarity index (SSIM) and feature similarity index metric (FSIM) values obtained for super resolution of eight MR image volumes by upscale factors 2 and 3.
Jimaging 07 00101 g010
Figure 11. Illustration of super resolution results for structural deformity cavernoma in real MR images, by different algorithms: (a) nearest neighbor, (b) spline interpolation, (c) non local means in three dimensions (NLM3D) [13], (d) low rank total variation based method (LRTV) [12], (e) the proposed approach, and (f) the original HR image. Each slice in axial, sagittal, and coronal planes is shown. The zoomed version of the cavernoma region from the coronal slice is highlighted in red rectangle. Please zoom for better visualization.
Figure 11. Illustration of super resolution results for structural deformity cavernoma in real MR images, by different algorithms: (a) nearest neighbor, (b) spline interpolation, (c) non local means in three dimensions (NLM3D) [13], (d) low rank total variation based method (LRTV) [12], (e) the proposed approach, and (f) the original HR image. Each slice in axial, sagittal, and coronal planes is shown. The zoomed version of the cavernoma region from the coronal slice is highlighted in red rectangle. Please zoom for better visualization.
Jimaging 07 00101 g011
Figure 12. Optimal selection of parameter λ in the super resolution approach. A random MR image is selected and upscaled by factor 2 using the proposed super-resolution (SR) approach. The HR image is obtained with different values of λ in Equation (3), and the obtained PSNR values are plotted.
Figure 12. Optimal selection of parameter λ in the super resolution approach. A random MR image is selected and upscaled by factor 2 using the proposed super-resolution (SR) approach. The HR image is obtained with different values of λ in Equation (3), and the obtained PSNR values are plotted.
Jimaging 07 00101 g012
Figure 13. Demonstration of reconstruction quality of different kinds of regions/patches, region with edges (green boxes), region with texture (red boxes), and smooth region (cyan boxes), after super resolving real MR images degraded with downsampling factor 2 and 2% noise, using different algorithms. (a) NN interpolation of noisy LR image, (b) spline interpolation of denoised LR image, (c) NLM3D [13] applied on denoised LR image, (d) LRTV [12] applied on denoised LR images, (e) the proposed work applied on a noisy LR image, and (f) the original noise-free HR image.
Figure 13. Demonstration of reconstruction quality of different kinds of regions/patches, region with edges (green boxes), region with texture (red boxes), and smooth region (cyan boxes), after super resolving real MR images degraded with downsampling factor 2 and 2% noise, using different algorithms. (a) NN interpolation of noisy LR image, (b) spline interpolation of denoised LR image, (c) NLM3D [13] applied on denoised LR image, (d) LRTV [12] applied on denoised LR images, (e) the proposed work applied on a noisy LR image, and (f) the original noise-free HR image.
Jimaging 07 00101 g013
Figure 14. Analysis of FSIM and SSIM values for different noise variances and upscale factors.
Figure 14. Analysis of FSIM and SSIM values for different noise variances and upscale factors.
Jimaging 07 00101 g014
Figure 15. Illustration of super resolving the real MR images of an Alzheimer’s subject degraded with downsampling factor 2 and 2% noise. From Left: Interpolated noisy LR image, spline interpolation of denoised image, NLM3D-based super-resolved denoised image, LRTV-based super resolution of the denoised image, the proposed approach super resolving the noisy LR image, and the original noiseless HR image.
Figure 15. Illustration of super resolving the real MR images of an Alzheimer’s subject degraded with downsampling factor 2 and 2% noise. From Left: Interpolated noisy LR image, spline interpolation of denoised image, NLM3D-based super-resolved denoised image, LRTV-based super resolution of the denoised image, the proposed approach super resolving the noisy LR image, and the original noiseless HR image.
Jimaging 07 00101 g015
Table 1. Quantitative comparison (peak signal to noise ratio (PSNR) and structural similarity index (SSIM) values) of the proposed denoising approach with existing approaches. Each row indicates a different method and each column indicates a different noise level (in %, i.e., 2% means noise variance is 255 × 0.02 = 5.1 for image with 255 maximum intensity value) with mentioned distribution (Rician/Gaussian) and mentioned datasets (real/simulated dataset).
Table 1. Quantitative comparison (peak signal to noise ratio (PSNR) and structural similarity index (SSIM) values) of the proposed denoising approach with existing approaches. Each row indicates a different method and each column indicates a different noise level (in %, i.e., 2% means noise variance is 255 × 0.02 = 5.1 for image with 255 maximum intensity value) with mentioned distribution (Rician/Gaussian) and mentioned datasets (real/simulated dataset).
I/P Rician NoiseGaussian Noise
MethodAssumed
Noise
Need
Var.
Simulated DatasetReal DatasetSimulated DatasetReal Dataset
3%5%7%9%11%3%5%7%9%11%3%5%7%9%11%3%5%7%9%11%
2D Image based MethodsNoisy
Image
--30.4126.0023.1220.9819.2820.4519.8719.1418.3217.4730.3925.9523.0320.8519.1020.3819.7018.8417.9116.96
0.84010.76230.69460.63690.58740.50030.49140.48280.47470.46690.83970.75980.68860.62710.57460.49840.48620.47300.46000.4476
NLMRiceYes35.8832.6230.2328.3726.7434.5431.4529.0627.1825.6435.8732.6630.3628.7227.4834.6631.8829.9328.5527.55
0.92390.88660.84690.80820.76940.89270.83010.77050.72040.68040.92370.88570.84450.80730.77630.89420.83720.78750.74890.7192
UNLMRiceYes36.3533.2130.9029.1327.7734.4431.6729.7828.4227.3836.2632.9430.3028.1026.1734.4931.6629.4627.4625.50
0.93070.89900.86390.82960.79840.89550.84990.80180.75690.71910.93040.89660.85220.80200.74940.89670.85370.80830.76380.7197
VST
BM3D
RiceYes36.2833.5431.7530.3629.2435.1932.8731.3330.0728.9036.1232.9130.2827.9325.8635.0832.1629.3126.4924.06
0.93040.90700.88580.86350.84230.90180.87010.84090.81120.78100.92970.89990.86070.81360.76340.90240.87010.83340.78610.7317
ProposedNillNo34.2132.5830.1928.2125.8034.5732.0429.0827.0825.6633.8632.1430.1428.2625.9334.1032.0329.7426.8024.08
0.92350.90100.85640.80450.76230.89030.85600.79890.72470.68620.91920.89530.85850.81820.75590.88280.85930.81070.73500.5677
3D Volume based MethodsORNLMRiceYes36.3433.4531.5329.7427.6335.1332.3130.5329.2028.2836.0033.0331.3730.2329.3535.5433.1531.5530.2929.32
0.93120.90670.88190.84920.80160.91300.87190.84000.81320.79020.93080.90510.88430.86640.85060.91640.88210.85360.82950.8080
AONLMRiceYes36.2833.6431.9530.6129.4834.7432.7131.1229.9528.7736.1033.4331.7430.4729.4434.6732.6831.2730.1629.24
0.92950.90560.88420.86280.84210.90610.87940.84980.82420.79640.92930.90410.88220.86150.84220.90440.87800.85340.83040.8090
MRNLMRiceYes36.0533.2631.3829.6427.5935.4932.8431.0029.5828.6035.4732.5730.9629.8629.0035.5233.0531.3629.9828.95
0.93290.90910.88520.85430.80870.91600.87950.84840.82120.79910.93140.90500.88370.86550.84940.91620.88220.85310.82740.8057
ODCTRiceYes35.8832.6930.3528.2626.1331.2028.9627.4727.9427.7835.8732.4629.9428.0526.5932.1230.8330.2828.5027.18
0.92950.90170.86940.82890.77540.85750.81030.77790.78610.77650.93030.90010.86540.83270.80190.87420.85470.84110.79320.7496
PRINLMRiceYes36.2632.8730.2928.1125.7930.9327.3426.8327.6427.1836.0632.2829.5127.5526.0331.7030.8730.1328.0026.61
0.93730.91060.87570.82810.76160.85330.80910.77290.79560.77140.93630.90350.86340.82290.77930.86770.85940.83720.77510.7129
VST
BM4D
RiceYes36.6533.9032.1630.8029.7135.4833.1431.5430.2529.2936.5733.4230.7528.3526.1535.4332.5429.6326.7224.21
0.93530.91420.89450.87460.85550.90570.87650.84960.82270.79850.93580.90890.87150.82710.77770.90670.87720.84090.79360.7400
VST
NLM3D
RiceYes36.0433.2131.3429.8828.6235.0032.5730.9129.5728.4235.9332.6529.9627.5225.2934.8631.8928.9325.9123.49
0.92780.90070.87540.85080.82520.90610.87300.84240.81300.78380.92720.89290.85100.80240.74850.90660.87200.82960.77310.7146
Here, red color denotes the highest PSNR/SSIM values obtained, bold numbers denote the PSNR and SSIM values obtained using the proposed work. Blue color denotes the PSNR/SSIM values comparable (within 0.5 dB PSNR range) to the proposed work obtained using different approaches.
Table 2. PSNR, SSIM, and FSIM values obtained using different algorithms for super resolution of real MR images in all three directions by factors 2 and 3.
Table 2. PSNR, SSIM, and FSIM values obtained using different algorithms for super resolution of real MR images in all three directions by factors 2 and 3.
SubjectMetricNNSplineNLM3DLRTVProposed
SRF = 21PSNR29.098129.940631.210333.416034.63
SSIM0.72750.75750.80790.87590.8859
FSIM0.88310.88560.95190.95970.9484
2PSNR29.2129.9631.1033.1233.56
SSIM0.75460.77500.81090.86370.8638
FSIM0.89990.90120.92970.95380.9697
3PSNR27.7328.4329.5831.6132.21
SSIM0.76820.79000.82410.87270.8765
FSIM0.90070.90560.91980.94670.9668
4PSNR27.9228.4729.5631.3132.20
SSIM0.82690.84370.87240.90890.9161
FSIM0.90510.90830.92170.94820.9639
5PSNR29.7130.1531.4233.4033.79
SSIM0.73890.78210.84130.87860.8819
FSIM0.89840.90840.91020.91940.9248
6PSNR28.6729.1530.5532.0032.82
SSIM0.70090.71040.82770.86950.8762
FSIM0.88720.89720.90990.91850.9274
7PSNR27.4428.7732.4434.2435.40
SSIM0.71030.74970.84360.90830.9293
FSIM0.89800.90690.91030.92910.9413
8PSNR28.2328.7029.9531.7332.01
SSIM0.75110.76530.79070.84070.8683
FSIM0.87200.90140.91040.92480.9236
SRF = 31PSNR26.2327.1228.3028.4128.89
SSIM0.69410.73080.77870.79860.8066
FSIM0.87490.88400.90300.90960.9182
2PSNR27.4428.3129.6029.7630.26
SSIM0.71630.75340.79630.81080.8180
FSIM0.89010.90910.92850.93960.9390
3PSNR26.1027.7928.7928.8529.11
SSIM0.65570.69430.74410.76470.7699
FSIM0.84980.87400.89190.90380.9148
4PSNR27.0328.2029.3029.2429.57
SSIM0.62710.66350.71920.74420.7419
FSIM0.85710.87200.89040.90810.9117
5PSNR27.9728.9030.3030.3230.67
SSIM0.58380.63720.71210.73420.7448
FSIM0.83010.85390.89490.90380.9175
6PSNR26.3427.8928.9729.1329.45
SSIM0.61030.76120.77310.78120.7865
FSIM0.85340.86410.88370.90880.9130
7PSNR26.7727.6228.9028.9229.32
SSIM0.66680.70590.76120.78070.7875
FSIM0.86930.88410.89220.90760.9127
8PSNR27.1428.0429.1829.2729.64
SSIM0.61150.65510.71210.73410.7421
FSIM0.86100.88010.90670.91350.9195
Bold numbers indicate the best performance of respective metrics.
Table 3. Quantitative results for different algorithms for super resolution of real MR images with upscaling factors 2 and 3, degraded with synthetically added 2% noise.
Table 3. Quantitative results for different algorithms for super resolution of real MR images with upscaling factors 2 and 3, degraded with synthetically added 2% noise.
SubjectsMetricNNSplineNML3DLRTVProposed
SRF = 21PSNR29.0129.5030.1731.2931.95
SSIM0.68930.70990.75620.79720.7983
FSIM0.88200.88410.91020.94040.9546
2PSNR29.4730.2530.9732.2032.81
SSIM0.72340.78020.79710.80680.8114
FSIM0.89870.89810.91780.94560.9537
3PSNR27.9828.7029.4430.6131.26
SSIM0.70120.72890.73130.80290.8032
FSIM0.89790.89860.91860.94430.9530
4PSNR31.6432.9633.5934.6434.52
SSIM0.79010.80010.81620.85120.8524
FSIM0.90200.90280.91980.94180.9507
5PSNR28.3129.0029.6830.6631.23
SSIM0.71730.72140.74450.79910.8041
FSIM0.89040.90080.90380.90480.9121
6PSNR29.1329.5430.2631.4331.97
SSIM0.75290.77120.78300.80900.8095
FSIM0.86940.90270.90710.91030.9174
7PSNR29.2129.8030.4931.5432.32
SSIM0.73740.74420.78340.81860.8291
FSIM0.89430.90410.90940.91940.9255
8PSNR28.4729.1229.8330.8631.47
SSIM0.74040.76190.79100.81380.8108
FSIM0.86320.89260.90880.91220.9198
SRF = 3 1PSNR26.4526.9727.8727.6228.08
SSIM0.71920.72310.74500.74010.7532
FSIM0.86110.88140.90060.90590.9134
2PSNR27.1927.8629.1728.9929.26
SSIM0.72670.73420.76620.75800.7667
FSIM0.88540.90720.92230.92320.9280
3PSNR26.9727.2428.4828.2828.49
SSIM0.71420.72080.72210.72410.7354
FSIM0.84430.86540.89070.90010.9062
4PSNR26.9527.3528.7228.2928.93
SSIM0.67070.67440.68190.68040.6880
FSIM0.84860.86620.88880.90000.9026
5PSNR28.1628.8829.5228.9329.44
SSIM0.64230.64870.65600.64460.6689
FSIM0.82180.85010.88450.88750.8974
6PSNR26.9127.1228.3327.9028.43
SSIM0.71280.71460.72570.72250.7351
FSIM0.84240.86310.87340.90750.9047
7PSNR27.1927.5728.6228.2928.83
SSIM0.63670.64560.67010.66470.6807
FSIM0.85340.87320.89820.90410.9105
8PSNR27.0527.6728.4929.0929.63
SSIM0.68480.68960.70150.69710.7035
FSIM0.85170.87600.90280.90640.9142
Bold numbers indicate the best performance of respective metrics.
Table 4. PSNR values for reconstructed images of randomly selected subject degraded with downsampling factor 2 and different levels of noise.
Table 4. PSNR values for reconstructed images of randomly selected subject degraded with downsampling factor 2 and different levels of noise.
Noise LevelInterpolationNLM3D [13]LRTV [12]Proposed
1%30.1230.531.6732.19
2%30.2630.5831.4331.78
3%29.7530.0030.3430.45
4%29.5429.6229.2829.34
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kaur, P.; Sao, A.K.; Ahuja, C.K. Super Resolution of Magnetic Resonance Images. J. Imaging 2021, 7, 101. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7060101

AMA Style

Kaur P, Sao AK, Ahuja CK. Super Resolution of Magnetic Resonance Images. Journal of Imaging. 2021; 7(6):101. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7060101

Chicago/Turabian Style

Kaur, Prabhjot, Anil Kumar Sao, and Chirag Kamal Ahuja. 2021. "Super Resolution of Magnetic Resonance Images" Journal of Imaging 7, no. 6: 101. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7060101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop