 TECHNICAL ADVANCE
 Open Access
 Published:
Sparse reconstruction of compressive sensing MRI using crossdomain stochastically fully connected conditional random fields
BMC Medical Imaging volume 16, Article number: 51 (2016)
Abstract
Background
Magnetic Resonance Imaging (MRI) is a crucial medical imaging technology for the screening and diagnosis of frequently occurring cancers. However, image quality may suffer from long acquisition times for MRIs due to patient motion, which also leads to patient discomfort. Reducing MRI acquisition times can reduce patient discomfort leading to reduced motion artifacts from the acquisition process. Compressive sensing strategies applied to MRI have been demonstrated to be effective in decreasing acquisition times significantly by sparsely sampling the kspace during the acquisition process. However, such a strategy requires advanced reconstruction algorithms to produce high quality and reliable images from compressive sensing MRI.
Methods
This paper proposes a new reconstruction approach based on crossdomain stochastically fully connected conditional random fields (CDSFCRF) for compressive sensing MRI. The CDSFCRF introduces constraints in both kspace and spatial domains within a stochastically fully connected graphical model to produce improved MRI reconstruction.
Results
Experimental results using T2weighted (T2w) imaging and diffusionweighted imaging (DWI) of the prostate show strong performance in preserving fine details and tissue structures in the reconstructed images when compared to other tested methods even at low sampling rates.
Conclusions
The ability to better utilize a limited amount of information to reconstruct T2w and DWI images in a short amount of time while preserving the important details in the images demonstrates the potential of the proposed CDSFCRF framework as a viable reconstruction algorithm for compressive sensing MRI.
Introduction
Magnetic Resonance Imaging (MRI) is a medical imaging technology that is currently used for diagnostic imaging of a wide range of diseases. In particular, since MRI does not use ionizing radiation, it has become a crucial imaging modality for screening frequently occurring cancers such as prostate cancer in men, breast cancer in women, as well as lung and colorectal cancer for both men and women. In 2015, 196,900 new cases of cancer (excluding nonmelanoma skin cancers) were expected, with 51 % of these belonging to the four aforementioned types of cancer in Canada [1]. As such, cancer screening methods with accurate and reliable information such as MRI is highly desired. Of particular interest for cancer screening is multiparametric MRI (MPMRI) since more information can be acquired through different modalities. MPMRI contains different techniques such as diffusion weighted imaging (DWI), correlated diffusion imaging (CDI) [2–4], dynamic contrast enhancement (DCE), T2weighted (T2w) imaging, and T1weighted (T1w) imaging [5]. Although this approach provides a more complete information, acquisition times are significantly longer which causes more patient discomfort and motion artifacts that decrease image quality. As a result, new methods to improve MRI acquisition times are highly desired to facilitate for reliable MPMRI data acquisition.
Compressive sensing has been demonstrated to be an effective strategy for reducing MRI acquisition times by acquiring significantly fewer samples in kspace. A complete signal can then be fully reconstructed through sparse, yet sufficient number of samples [6–8]. In MRI, compressive sampling strategies have been demonstrated to be highly effective at reducing acquisition time while maintaining image quality as different types of tissue structure have been shown to be sparse in certain domains [9]. Furthermore, different techniques have been proposed to improve the imaging process [10] as well as the reconstruction process [11–23] in compressive sensing. Due to the limited amount of data available through compressive sensing, advanced reconstruction algorithms are required to produce high quality and reliable images.
Different methods have been proposed for sparse reconstruction of compressive sensing MRI [11–23]. As a notable example, Block [14] proposed an iterative image reconstruction technique using a modified total variation (TV) constraint [20, 21] for sparse reconstruction of compressive sensing brain MRI. Trzasko [15] introduced a homotopic l _{0} minimization method for the sparse reconstruction of compressive sensing spinal MRI. Wong [12] extended upon this idea and proposed a regional sparsified domain for the sparse reconstruction of breast MRI. A similar technique was also demonstrated by Qu using combined sparsifying transforms and smoothed l_{0} norm minimization [13], where they showed that the use of combined transforms can improve image quality comprised of the reconstructed images from compressive sensing MRI when compared to methods using a single sparsifying transform. However, the downside of the l_{0} norm minimization is the fact that its performance significantly depends on the tuning parameters where these tuning parameters can greatly affect the convergence rate of the algorithm. Other l_{ n } optimization techniques such as the standard l_{2} (least squares) minimization can have high error rates as reported in [24].
An area that is little explored but can reap significant potential benefits is the application of random field modeling for improved sparse reconstruction of compressive sensing MRI. Random field modeling such as Markov random fields (MRF) [25, 26] and conditional random fields (CRF) [27] have long been shown to be powerful tools for incorporating spatial context within a probabilistic graphical modeling framework, which can have significant benefits for reconstructing images from sparse measurements. Despite powerful modeling capabilities and potential benefit to sparse reconstruction, one of the biggest hurdles in leveraging random field models for compressive sensing MRI is the fact that all MRI measurements are made in kspace, whereas the images are reconstructed in spatial domain. As the majority of random field models are typically modeled in a single domain, such models cannot be used directly for the purpose of sparse reconstruction of compressive sensing MRI. This is further complicated by the fact that the MRI measurements in kspace are sparse and incomplete, which make it difficult to leverage existing random field models for this problem. Therefore, a probabilistic graphical modeling framework that can consolidate the fact that partial measurements are made in a domain different than the desired states of the reconstruction images is needed to truly leverage the power of random field modeling for sparse reconstruction of compressed sensing MRI.
This paper proposes a crossdomain Stochastically fully connected conditional random field (CDSFCRF) approach for the reconstruction of compressive sensing MRI at below Nyquist sampling rates [28]. Inspired by the stochastic cliques method presented in [29], the proposed cross domain model is a significant extension upon this method by consolidating two different domains of spatial domain and kspace specifically for MRI compressive sensing. CDSFCRF framework introduces constraints in both kspace and spatial domains within a stochastically fully connected graphical model to produce improved MRI reconstruction. The proposed CDSFCRF framework has the ability to utilize spatial and data driven consistencies in the spatial domain along with data driven consistencies in the kspace domain pertaining to sparse measurements while maintaining edge features and structural details in the reconstructed images. Phantom MRI data as well as prostate MRI data captured using T2w and DWI imaging modalities, which also yields apparent diffusion coefficient (ADC) map images, were used to illustrate the efficacy of the proposed CDSFCRF framework for sparse reconstruction of compressive sensing MRI. To the best of the authors’ knowledge, this is the first time that constraints in both kspace and spatial domains are used in conjunction within a stochastically fully connected graphical model for the sparse reconstruction of compressive sensing MRI, which is the main contribution of this paper.
The paper is formatted as follows. The methodology behind the proposed CDSFCRF framework is described in Section “Methodology”. The experimental setup is described in Section “Experimental setup”. Results and discussions are presented and discussed in Sections “Results” and “Discussion”, respectively. Finally, the conclusion is presented in Section “Conclusions”.
Methodology
In MRI, measurements are made in the kspace [30], with the lower frequency coefficients in the kspace containing coarsegrained contrast information while higher frequency coefficients contain finegrained image detail information. The MRI measurements from the kspace are transformed into the spatial domain to form the reconstructed MRI image. Most compressive sensing strategies [6, 15] sparsely sample the kspace to reduce image acquisition time significantly. Therefore, to fully utilize available information in the reconstruction process, datadriven constraints in the kspace domain and data and spatial driven constraints in the spatial domain would be highly beneficial in improving image reconstruction quality from compressive sensing MRI.
Motivated by this, the proposed crossdomain stochastically fully connected conditional random field (CDSFCRF) introduced here for the purpose of sparse reconstruction of compressive sensing MRI, extends upon the seminal work on stochastically fully connected conditional random fields (SFCRF) first proposed in [29] to facilitate for this crossdomain optimization. SFCRFs are fullyconnected conditional random fields with stochastically defined cliques. Unlike traditional conditional random fields (CRF) where nodal interactions are deterministic and restricted to local neighborhoods, each node in the graph representing a SFCRF is connected to every other node in the graph, with the cliques for each node is stochastically determined based on a distribution probability. Therefore, the number of pairwise cliques might not be the same as the number of neighborhood pairs as in the traditional CRF models. By leveraging longrange nodal interactions in a stochastic manner, SFCRFs facilitate for improved detail preservation while maintaining similar computational complexity as CRFs, which makes SFCRFs particularly enticing for the purpose of improved sparse reconstruction of compressive sensing MRI. However, here the problem is to reconstruct an MRI image in the spatial domain while the available measurements are made in kspace domain. Similar to most CRF models, SFCRFs cannot be leveraged directly for this purpose. Motivated by the significant potential benefits of using SFCRFs in improving reconstruction quality of compressive sensing MRI, we extend the SFCRF model into a crossdomain stochastically fully connected conditional random field (CDSFCRF) model that incorporates crossdomain information and constraints from kspace and spatial domains to reconstruct the desirable MRI image from sparse observations in kspace.
The theory pertaining to sparse reconstruction via a crossdomain stochastically fully connected conditional random field model is detailed in Appendix 1.
Implementation
An implementation of the proposed CDSFCRF framework for the purpose of sparse reconstruction from compressive sensing MRI is illustrated in Fig. 1. Here, an iterative gradient descent optimization approach is employed, and can be described as follows. First, the original compressive sensing MRI data in kspace is transformed to the spatial domain to provide an initial estimate of the reconstructed image. Second, the gradient of the unary and pairwise energy potentials is computed, where the unary data driven consistencies with respect to the original observations are enforced in the kspace, and spatial and data driven consistencies are enforced in the spatial domain. Third, the estimate of the reconstructed image is updated based on the previous estimate and the computed gradient. The second and third steps of this process are repeated until convergence.
Experimental setup
To study the efficacy of the proposed CDSFCRF method for the purpose of sparse reconstruction of compressive sensing MRI, experiments were performed including: i) MRI data acquired of a MRI training phantom, and ii) prostate MPMRI data of 20 patient cases. A detailed description of the phantom data, patient data, and MRI image acquisition procedure to facilitate for the various experiments are described below.
Phantom data
The MRI training phantom used in the experiments, shown in Fig. 2, was a multimodality prostate training phantom from Computerized Imaging Reference Systems Inc (CIRCS MODEL 053). The phantom is composed of a clear acrylic container with dimensions 11.5×7.0×9.5 cm with a front probe opening of 3.2 cm diameter and a rear probe opening of 2.6 cm diameter. The prostate is composed of highscattering Blue Zerdine with dimensions 5×4.5×4.0 cm and is placed in a background gel similar to water with little backscatter attenuation (≤0.07 dB/cm−MHz). Within the prostate, there are 3 randomly placed lesions of sizes between 0.5−1.0 cm placed hypoechoic to the prostate. The urethra and rectal wall are made of low scattering Zerdine with diameter of 0.7 cm with dimensions 6×11×0.5 cm, respectively. This phantom was imaged with an inflatable Medrad eCoil ERC using DWI. The DWI MRI was acquired by a 3T GE Discovery MR750. DWI was collected at b=0 mm^{2}/s at 3NEX^{2}. For the DWI data, the echo time (TE) was 71.70 ms and repetition time (TR) was 10,000.00 ms.
Patient data experiments
To test the efficacy of the proposed CDSFCRF framework within a real clinical scenario, MRI data of 20 patients (17 with cancer and 3 without cancer) were acquired using a Philips Achieva 3.0T machine at Sunnybrook Health Sciences Centre, Toronto, Ontario, Canada. All data was obtained retrospectively under the local institutional research ethics board (Research Ethics Board of Sunnybrook Health Sciences Centre). For each patient, the following MPMRI modalities were obtained (Table 1): T2w and DWI. The patients’ age ranged from 53 to 83. Table 1 summarizes the information about the 20 patients’ datasets used in this study, which includes displayed field of view (DFOV), resolution, echo time (TE), and repetition time (TR).
Compressed sensing configuration
In order to evaluate the efficacy of the proposed CDSFCRF framework at different sample rates, we first acquired MRI measurements at all kspace coefficients. Based on this fullysampled set of kspace measurements, sparse sampling was then conducted using radial sampling patterns with different numbers of radial sampling lines to achieve a desired sampling rate. For example, Fig. 3 shows a radial sampling pattern which corresponds to a sampling rate of 32 % of the kspace. Different sampling rates were tested and evaluated in this study.
Results
In order to evaluate the efficacy of the proposed CDSFCRF framework for sparse reconstruction of compressive MRI sensing, a comparative evaluation analysis was performed alongside a baseline l_{2} minimization (L_{2}) reconstruction method, and a stateoftheart homotopic l_{0} minimization (HL_{0}) [15] reconstruction method. The tested methods were compared quantitatively through peak signaltonoise ratio (PSNR) analysis, and qualitatively via visual assessment. All tested methods were implemented based on the original literature, with optimal parameters used in this study. All tested methods were run until convergence.
Figure 4 shows the PSNR versus sampling rate plots for the tested methods for the phantom MRI data.
Tables 2, 3, and 4 show the PSNR results for the three reconstructed methods for the T2w, DWI, as well as ADC map images for the patient data experiments at different sampling rates.
Figures 5 and 6 shows the visual comparison between the reconstructed images produced using the proposed CDSFCRF framework compared with that produced using the L_{2} and homotopic l_{0} minimization reconstruction methods for three cases for T2w images.
Figures 7 and 8 shows the visual comparison between the reconstructed images produced using the proposed CDSFCRF framework compared with that produced using the L_{2} and HL_{0} methods for three patient cases for DWI (b=100s/mm^{2}) and ADC maps.
Discussion
As it can be observed from Fig. 4, the proposed CDSFCRF framework achieved noticeable PSNR improvements over the other tested methods at all tested sampling rates. The CDSFCRF produced improvements of up to 4 dB over HL_{0} and 7 dB over L_{2} in low sampling conditions. It can also be observed that as sampling rates increase, the performance differences decrease. This is due to the fact that as the sampling rate increases, the amount of available measurements increases, and as such the level of reconstruction quality improvements that can be achieved will naturally decrease given the amount of available information becomes increasingly sufficient for high quality reconstruction. The ability of the CDSFCRF framework to produce high quality reconstruction at very low sampling rates can be demonstrated visually as well.
From additional quantitative analysis of patients MRI data presented in Tables 2, 3, and 4, it can be observed that the proposed CDSFCRF framework achieved the greatest PSNR improvements for the lowest sampling rate (i.e., 17 %) where for T2w, CDSFCRF improved PSNR by 1.78 dB and 1.12 dB over the L_{2} and HL_{0} methods, respectively. For DWI, CDSFCRF improved PSNR by 1.85 dB and 0.28 B over the L_{2} and HL_{0} methods, respectively. Interestingly for ADC maps, the best improvements in PSNR were achieved for the highest sampling rate (47 %) where CDSFCRF improved PSNR by 4.44 dB and 0.21 B over the L_{2} and HL_{0} methods, respectively.
Tables 2, 3 and 4 also show the Pvalues calculated by comparing the proposed CDSFCRF method with L_{2} and HL_{0} methods, respectively. As it can be seen, Pvalues show significant difference between CDSFCRF and the other two methods for T2w and DWI images. For ADC maps, the proposed CDSFCRF was significantly different than L_{2} method as well. The only comparison that did not show significantly different results was CDSFCRF compared to HL_{0} method for ADC maps. This shows that the PSNR improvement for the proposed CDSFCRF framework was meaningful for the majority of cases when compared to other tested methods.
Comparing the results for phantom MRI data (Fig. 4) and patients MRI data shown in Tables 2, 3 and 4 shows that the proposed CDSFCRF framework yields higher performance improvement at 10  20 % sampling range for phantom MRI data compared to patients MRI data. The reason for this difference on PSNR improvement is the fact that the morphological and textural properties of the phantom is significantly less complex than that of real patients’ prostates, and thus the reconstruction problem is a simpler one for the phantom and as a result, greater PSNR gains were achieved using the proposed method.
Qualitative observations from Figs. 5 and 6 show that the L_{2} method resulted in blurry T2w images as well as noticeable radial artifacts at low sampling rates as expected due to the least squares reconstruction being prone to errors. The HL_{0} approach performed better than the L2 minimization and was able to noticeably reduce artifacts and provide a higher quality reconstruction. However, in comparison, the CDSFCRF was able to better restore details and fine tissue structure in the reconstructed image when compared to HL_{0}. This is to be expected as the CDSFCRF takes advantage of more complete data and spatial driven consistencies in a fully connected nature, thus better modeling the underlying tissue detail and structures.
Furthermore, as it can be seen in Figs. 7 and 8, the L_{2} method resulted in blurry DWI and ADC map images again with noticeable radial artifacts. Although the HL_{0} approach performed better than the L_{2} method, it can be observed once again that the proposed CDSFCRF approach was able to preserve more fine tissue structure and detail in the reconstructed image when compared to the HL_{0} method. Nevertheless, an inherent tradeoff exists between preserving fine textural granularity and reducing artifacts due to compressed sensing which can be well utilized in the proposed CDSFCRF framework to achieve a balance between the two competing constraints.
In Figs. 5, 6, 7 and 8, the tumourous regions marked by a radiologist and confirmed by pathology report (biopsy results) are shown by red arrow or white boundary. It can be seen that the proposed CDSFCRF method preserves the separability of the cancerous and healthy tissue in all cases, which is an important measure for usability of the proposed method in practice. As it can be seen, the tumourous regions are blurred in the L_{2} method, which may make it difficult to detect for radiologists.
Both quantitative and qualitative analysis demonstrate the potential of the proposed CDSFCRF framework as a reliable reconstruction approach for compressive sensing in MRI. It demonstrates the ability to produce edge and tissue details at very low sampling rates. The CDSFCRF framework better utilize available information to produce high quality reconstructed images given very limited available information. Preservation of tissue structure and detail enhancement, and noise and artifact mitigation are very important for MRI as the diagnostic quality is directly related to the image quality. This demonstrates that the CDSFCRF framework can be a viable clinical technique as the reduction in acquisition can lead to faster acquisitions and lower patient wait times. With a lower acquisition time and hence lower patient wait time, patients can have access to the necessary treatments in a timely manner, significantly improving the patient outcome and survival rates.
The compressive sensing method used to reconstruct MR images can influence the performance of the computeraided diagnosis (CAD) tools. For example, several radiomicsbased CAD algorithms have been proposed for automatic prostate cancer detection which use T2w and DWI to extract texture and morphological features fed into a classifier [31–36]. These algorithms heavily rely on the quality of regions of interests in similar cases in DWI and therefore, it is expected that a reconstructed MRI with better quality will improve the performance. As future work, we will investigate the effect of the proposed compressive sensing method on the detection accuracies of these radiomicsbased CAD algorithms with respect to the L_{2} and HL_{0} methods. Moreover, recently, computational diffusion MRI (CDMRI) has been introduced which utilizes the wealth of information in DWMRI to computationally construct new sequences of MRI that potentially will help radiologists with more accurate and consistent diagnosis [2, 3]. The proposed CDSFCRF framework will be integrated into CDMRI algorithms [2, 3] to investigate whether CDSFCRF improves the separability of cancerous and healthy tissues in prostate for these computationally generated MR sequences with respect to the L_{2} and HL_{0} methods.
The limitations of the proposed CDSFCRF method that will be addressed in the future direction of this work include the limited sample size. A larger and more diverse dataset will be used to address this limitation. Moreover, in this work, the proposed method was applied only to the prostate. Future work also includes applications of the proposed method to the MRI acquisitions of other organs such as breast or moving organs such as heart. In addition, although the CDSFCRF can significantly decrease MRI acquisition times, because of the fullyconnected nature of this method, the algorithm may require a considerable processing time to complete (although not comparable to original MRI acquisition time). As future work, we will modify the proposed method to improve processing time and the efficiency of the algorithm.
Conclusions
In this paper, a cross domain stochastic fully connected conditional random field (CDSFCRF) framework for sparse reconstruction of compressive sensing MRI was presented. The proposed CDSFCRF framework introduces constraints in both kspace and spatial domains within a stochastically fully connected graphical model to produce improved MRI reconstruction. To test the efficacy of the proposed CDSFCRF framework, quantitative experimentation using peak signaltonoise ratio (PSNR) analysis was performed on phantom MRI data. Quantitative and qualitative experimentations were also performed on prostate MPMRI data of 20 patient cases at different sampling rates. The results show an improvement over other tested sparse reconstruction approaches, especially at low sampling rates. The ability to better utilize available information given very limited information demonstrates the potential of the proposed CDSFCRF framework as a viable reconstruction algorithm for compressive sensing MRI. The proposed CDSFCRF can significantly reduce MRI acquisition times without sacrificing quality and potential reduction in the accuracy of diagnosis. Reducing MRI acquisition time would reduce related cost significantly and lead to less patient discomfort during the MRI acquisition and more importantly, it would reduce the patient wait times considerably. A fast access to MRI would directly translate to better care given to patients who need it the most.
Appendix 1
Sparse reconstruction via crossdomain stochastically fully connected conditional random field
The main goal here is to reconstruct image Y given original sparsely sampled kspace observations X. We model the conditional probability P(YX) of the full state set Y in spatial domain given the set of sparse measurements X in kspace, which can be written as:
where Z(X) is the normalization function and ψ(.) is a combination of unary and pairwise potential functions:
Here y_{ i }∈Y is a single state in the set \(Y=\{y_{i}\}^{n}_{i=1}\), y_{ φ }∈Y encodes a clique structure in the set C, and \(X= \{x_{j}\}^{n}_{j=1}\) is the observations (radially subsampled frequency coefficients) in the frequency domain (kspace). The unary potential ψ_{ u } is enforced in the kspace while the pairwise potential ψ_{ p } is applied in the spatial domain. The unary potential enforces original observations to preserve data fidelity. Since the available observations are captured in kspace in MRI, the model must be formulated in a way to be consistent in both kspace and spatial domain.
The pairwise potential, on the other hand, has to be in the spatial domain to better preserve image detail since neighboring coefficients in the kspace does not contain any meaningful spatial or data consistencies to be utilized by the pairwise potential. Therefore, the optimal way to fully utilize available data within this random field model is to formulate the unary potential in the kspace and the pairwise potential in the spatial domain.
One of main differences between the proposed CDSFCRF framework from conventional CRF models is to incorporate longrange information in the model and preserve boundaries and image structural properties more effectively which is important here due to sparse available observation. To capture longrange information, CDSFCRF assumes fully connected neighboring structure for the underlying graph which each node i has a set of neighbors
where N(i)=n−1 and includes all other nodes in the graph as neighbors of node i. Here the pairwise clique structures are utilized such that:
The active cliques in the inference procedure are determined by the stochastic indicator function \(1^{S}_{\left \{i,j\right \}}=1\). The indicator function decides whether or not nodes can construct a clique, C_{ p }(i) for node i. This stochastic indicator function combines spatial and data driven information to model the probability distribution of informative cliques which informative cliques have higher probability to participate in the inference. The set of active cliques are obtained to extract pairwise potentials in Eq. 2.
As mentioned before, ψ(·) in Eq. 2 is the combination of two potential functions ψ_{ u }(.), the unary potential and ψ_{ p }(.), the pairwise potential. These potential functions are formulated with their corresponding weights λ, respectively as:
where λ controls the importance of each feature function in the energy formulation and it is calculated in the training stages. Although it is possible to provide several arbitrary feature functions to model the conditional probability P(YX), here two feature functions are provided to formulate the image reconstruction for the purpose of sparse reconstruction from compressive sensing MRI. The conditional distribution of Y given X is trained to promote/suppress different features in both the unary and pairwise potentials. Higher \({\lambda _{j}^{u}}\) values promotes a higher reinforcement of original observations while high \({\lambda _{k}^{p}}\) values promotes higher consideration of spatial and data driven neighborhood constraints. In Eq. 6, F refers to the frequency domain potential function. The unary potential is calculated in the kspace while the pairwise remains in the spatial domain. This is the novelty of the CDSFCRF whihc facilitates for better preservation of fine tissue details and contrast in the reconstructed image. The unary potential function F_{ j }(y_{ i },X) can be formulated as:
where \(\mathscr {F}(\cdot,\cdot)\) is the Fourier operator and returns the kspace coefficient corresponding to frequency ω. Based on this formulation, the unary potential is enforced in the kspace and in the inferencing step, the model tries to estimate image Y to be consistent to the original kspace observation \(X = \{x_{\omega }\}_{\omega =\frac {\pi }{2}}^{\frac {\pi }{2}}\).
The pairwise function f_{ k }(y_{ i },y_{ j },X) can be formulated as:
where σ is a control variable for the amount of weighting node pairs in the clique φ={i,j}. Contrary to the unary potential, the pairwise potential is enforced in the spatial domain.
Graph representation
Graph G(V,E) (Fig. 9) is the realization of the CDSFCRF where V is the set of nodes of the graph representing states \(Y=\{y_{i}\}^{n}_{i=1}\), E is the set of edges in the graph. Observations x_{ i }∈X are made in the kspace domain. Our final state estimations Y are in the spatial domain (image). Figure 9 shows the graphical representation of how the spatial and kspace domain are incorporated to model the conditional probability P(YX). x_{ i } comes from sparse measurements in the kspace. In the inference procedure, the kspace observations are transformed into the spatial domain using the Fourier transform to compute the pairwise potentials. Pairwise potentials are calculated in the spatial domain and transformed into the kspace to combine with the unary potential and perform data fidelity. For different types of MRI data, different sparse sampling patterns can be used. Furthermore, pairwise connectivity can be trained for specific types of details and tissue structure.
The proposed CDSFCRF framework utilizes consistencies from the spatial domain through the pairwise potential in conjunction with kspace information through the unary potential. A combination of the two potentials is enforced simultaneously. The unary potential utilizes original observations in the kspace, while the pairwise potential utilizes the spatial domain representation of the observation/state information and calculates pairwise potentials for nodes in the spatial domain. This allows CDSFCRF to take advantage of the lower computational complexity introduced by the stochastically fullyconnected random field model, while leveraging the original kspace observations in improving signal fidelity.
Abbreviations
 MRI:

Magnetic resonance imaging
 CRF:

Conditional random field
 CDSFCRF:

Crossdomain stochastically fully connected conditional random fields
 T2w:

T2weighted
 DWI:

Diffusionweighted imaging
 MPMRI:

Multiparametric MRI
 CDI:

Correlated diffusion imaging
 TV:

Total variation
 ADC:

Apparent diffusion coefficient
 PSNR:

Peak signaltonoise ratio
 stdev:

Standard deviation
 CAD:

Computeraided diagnosis
 CDMRI:

Computational diffusion MRI
References
 1
Canadian Cancer Statistics Special topic : Predictions of the future burden of cancer in Canada. Technical Report 2015 Canadian Cancer Society http://www.cancer.ca/statistics.
 2
Wong A, Glaister J, Cameron A, Haider MA. Correlated diffusion imaging. BMC Med Imaging. 2013; 13(1):26. doi:http://dx.doi.org/10.1186/147123421326.
 3
Wong A, Khalvati F, Haider MA. DualStage Correlated Diffusion Imaging. In: IEEE International Symposium on Biomedical Imaging (ISBI). Brooklyn: 2015. p. 75–8.
 4
Khalvati F, Wong A, Haider MA. Enhanced DualStage Correlated Diffusion Imaging. In: International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Orlando: 2016. p. 5537–40.
 5
Esen T, Turkbey B, Patel A, Futterer J. Multiparametric mri in prostate cancer. BioMed Research International. 2014; 2014:296810. doi:http://dx.doi.org/10.1155/2014/296810.
 6
Donoho DL. Compressed sensing. IEEE Trans Inform Theory. 2006; 52(4):1289–306. doi:http://dx.doi.org/10.1109/TIT.2006.871582.
 7
Candès EJ. Compressive sampling. Int Congr Math. 2006; 3:1433–52.
 8
Baraniuk R. Compressive Sensing. IEEE Signal Process Mag. 2007; 24:1–9.
 9
Lustig M, Donoho D, Pauly JM. Magn Reson Med. 2007; 58(6):1182–95. doi:http://dx.doi.org/10.1002/mrm.21391.
 10
Liang D, Liu B, Wang J, Ying L. Accelerating SENSE using compressed sensing. Magn Reson Med. 2009; 62(6):1574–84. doi:http://dx.doi.org/10.1002/mrm.22161.
 11
Ye JC, Tak S, Han Y, Park HW. Projection reconstruction MR imaging using FOCUSS,. Magn Reson Med. 2007; 57(4):764–5. doi:http://dx.doi.org/10.1002/mrm.21202.
 12
Wong A, Mishra A, Fieguth P, Clausi DA. Sparse reconstruction of breast MRI using homotopic L0 minimization in a regional sparsified domain. IEEE Trans Biomed Eng. 2013; 60(3):743–52. doi:http://dx.doi.org/10.1109/TBME.2010.2089456.
 13
Qu X, Cao X, Guo D, Hu C, Chen Z. Compressed Sensing MRI with Combined Sparsifying transforms and Smooted L0 Norm Minimization. In: IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP). Dallas: 2010. p. 626–9.
 14
Block KT, Uecker M, Frahm J. Magn Reson Med. 2007; 57(6):1086–98. doi:http://dx.doi.org/10.1002/mrm.21236.
 15
Trzasko J, Manduca A. Highly undersampled magnetic resonance image reconstruction via homotopicminimization. Med Imaging IEEE Trans. 2009; 28(1):106–21.
 16
Chartrand R. Fast algorithms for nonconvex compressive sensing: Mri reconstruction from very few data. In: Biomedical Imaging: From Nano to Macro, 2009. ISBI’09. IEEE International Symposium On. Boston: IEEE: 2009. p. 262–5.
 17
Yin W, Osher S, Goldfarb D, Darbon J. Bregman iterative algorithms for ∖ell_1minimization with applications to compressed sensing. SIAM J Imaging Sci. 2008; 1(1):143–68.
 18
Figueiredo MA, Nowak RD, Wright SJ. Gradient projection for sparse reconstruction: Application to compressed sensing and other inverse problems. Selected Topics Signal Process IEEE J. 2007; 1(4):586–97.
 19
Goldstein T, Osher S. The split bregman method for l1regularized problems. SIAM J Imaging Sci. 2009; 2(2):323–43.
 20
Wang Y, Yang J, Yin W, Zhang Y. A new alternating minimization algorithm for total variation image reconstruction. SIAM J Imaging Sci. 2008; 1(3):248–72.
 21
Osher S, Burger M, Goldfarb D, Xu J, Yin W. An iterative regularization method for total variationbased image restoration. Multiscale Model Simul. 2005; 4(2):460–89.
 22
Qu X. Compressed sensing MRI based on nonsubsampled contourlet transform. 2008 IEEE International Symposium on IT in Medicine and Education. 2008:693–696. doi:http://dx.doi.org/10.1109/ITME.2008.4743955.
 23
Yu Y, Hong M, Liu F, Wang H, Crozier S. Compressed sensing MRI using Singular Value Decomposition based sparsity basis. Conf Proc Ann Int Conf IEEE Eng Med Biol Soc. 2011; 2011(1):5734–7. doi:http://dx.doi.org/10.1109/IEMBS.2011.6091419.
 24
Aster R, Borchers B, Thurber C. Parameter estimation and inverse problems. Int Geophys. Series, 90: Elsevier; 2005.
 25
Held K, Kops ER, Krause BJ, Wells III WM, Kikinis R, MullerGartner HW. Markov random field segmentation of brain mr images. Med Imaging, IEEE Trans. 1997; 16(6):878–6.
 26
Blake A, Kohli P, Rother C. Markov Random Fields for Vision and Image Processing. Cambridge: MIT Press; 2011.
 27
Lafferty J, McCallum A, Pereira FCN. Conditional random fields: Probabilistic models for segmenting and labeling sequence data. In: ICML (International Conference on Machine Learning): 2001. p. 282–9.
 28
Nyquist H. Certain topics in telegraph transmission theory. Am Inst Electrical Eng Trans. 1928; 47(2):617–44.
 29
Shafiee MJ, Wong A, Siva P, Fieguth P. IEEE International Conference on Image Processing (ICIP). Paris; 2014, pp. 4289–93.
 30
Ljunggren S. A simple graphical representation of fourierbased imaging methods. J Magn Reson (1969). 1983; 54(2):338–43.
 31
Khalvati F, Wong A, Haider M. Automated Prostate Cancer Detection via Comprehensive MultiParametric Magnetic Resonance Imaging Texture Feature Models. Biomed Central Med Imaging. 2015:15–27. doi:10.1186/s1288001500699.
 32
Khalvati F, Modhafar A, Cameron A, Wong A, Haider M. A multiparametric diffusion magnetic resonance imaging texture feature model for prostate cancer analysis. In: MICCAI 2014 Workshop on Computational Diffusion MRI. Boston: 2014. p. 79–88.
 33
Cameron A, Khalvati F, Haider M, Wong A. MAPS: A Quantitative Radiomics Approach for Prostate Cancer Detection. IEEE Trans Biomed Eng. 2016; 63(6):1145–56.
 34
Cameron A, Modhafar A, Khalvati F, Lui D, Shafiee M, Wong A, Haider MA. Multiparametric MRI Prostate Cancer Analysis via a Hybrid MorphologicalTextural Model. In: International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Chicago: 2014. p. 3357–60.
 35
Chung A, Scharfenberger C, Khalvati F, Wong A, Haider M. Textural Distinctiveness in MultiParametric Prostate MRI for Suspicious Region Detection. In: International Conference on Image Analysis and Recognition (ICIAR). Ontario: 2015. p. 368–76.
 36
Zhang J, Khalvati F, Wong A, Haider M. Superpixelbased Prostate Cancer Detection from Diffusion Magnetic Resonance Imaging. Vision Lett. 2015; 1(1):107.
Funding
This research has been supported by the Canada Research Chairs programs, Natural Sciences and Engineering Research Council of Canada (NSERC), the Ministry of Research and Innovation of Ontario, and Ontario Institute of Cancer Research (OICR).
Availability of data and materials
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request pending the approval of the institute.
Authors’ contributions
EL, FK, MJS, and AW contributed to the design and implementation of the concept. EL, FK, MJS, and AW contributed to the design and implementation of the experiments, and performing statistical analysis. FK and MAH were involved in collecting and reviewing the data. All authors contributed to the writing and reviewing of the paper. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
The institutional research ethics board (Sunnybrook Health Sciences Centre Research Ethics Board) approved this retrospective single institution study and waived the requirement for informed consent.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Li, E., Khalvati, F., Shafiee, M.J. et al. Sparse reconstruction of compressive sensing MRI using crossdomain stochastically fully connected conditional random fields. BMC Med Imaging 16, 51 (2016). https://doi.org/10.1186/s1288001601566
Received:
Accepted:
Published:
Keywords
 Compressive sensing
 Conditional random fields
 Magnetic resonance imaging