Difference between revisions of "Magnetic Resonance Image Example Based Contrast Synthesis"
(Fig. 1 & 2.) |
|||
| (11 intermediate revisions by 2 users not shown) | |||
| Line 1: | Line 1: | ||
| − | <meta name="title" content="Magnetic Resonance Image Example-Based Contrast Synthesis"/> | + | <!-- <meta name="title" content="Magnetic Resonance Image Example-Based Contrast Synthesis"/> --> |
{{h2| Magnetic Resonance Image Example-Based Contrast Synthesis}} | {{h2| Magnetic Resonance Image Example-Based Contrast Synthesis}} | ||
{{iacl|~roy/|Snehashis Roy}}, {{iacl|~aaron/|Aaron Carass}}, and [[Prince|Jerry L Prince]] | {{iacl|~roy/|Snehashis Roy}}, {{iacl|~aaron/|Aaron Carass}}, and [[Prince|Jerry L Prince]] | ||
{{TOCright}} | {{TOCright}} | ||
| − | |||
| − | |||
{{h3|Introduction}} | {{h3|Introduction}} | ||
| Line 12: | Line 10: | ||
| − | {{h3| | + | {{h3|Method}} |
| + | <div style="background: white; border: 2px solid rgb(150, 150, 150); padding: 2px; text-align: justify;"> | ||
| + | {| | ||
| + | | | ||
| + | {| style="width: 45%" | ||
| + | | | ||
{| | {| | ||
| − | |||
| [[Image:MIMECS_Figure_2.png]] | | [[Image:MIMECS_Figure_2.png]] | ||
| − | |||
| − | |||
|- | |- | ||
| '''Figure 1:''' (a) Atlas T1-weighted SPGR and (b) its corresponding T1-weighted MPRAGE. (c) Subject T1-weighted SPGR scan and (d) its T1-weighted MPRAGE image. Atlas SPGR is deformably registered to the subject SPGR. This deformation is applied to the atlas MPRAGE to obtain (e) a synthetic subject MPRAGE. (f) Synthetic MPRAGE image generated by our algorithm, MIMECS. | | '''Figure 1:''' (a) Atlas T1-weighted SPGR and (b) its corresponding T1-weighted MPRAGE. (c) Subject T1-weighted SPGR scan and (d) its T1-weighted MPRAGE image. Atlas SPGR is deformably registered to the subject SPGR. This deformation is applied to the atlas MPRAGE to obtain (e) a synthetic subject MPRAGE. (f) Synthetic MPRAGE image generated by our algorithm, MIMECS. | ||
| + | |} | ||
| + | |} | ||
| + | |</div> | ||
| + | |<div style="background: white; border: 2px solid rgb(150, 150, 150); padding: 2px; text-align: justify;"> | ||
| | | | ||
| + | {| style="width:45%" | ||
| + | | | ||
| + | {| | ||
| + | | [[Image:MIMECS_FIGURE_3.png]] | ||
| + | |- | ||
| '''Figure 2:''' Left column shows images from four (out of 14) time-points of a normal BLSA subject, where each image was acquired approximately one year apart. Right column shows synthetic SPGR images, synthesized using MIMECS applied to the corresponding MPRAGE images in the left column. | | '''Figure 2:''' Left column shows images from four (out of 14) time-points of a normal BLSA subject, where each image was acquired approximately one year apart. Right column shows synthetic SPGR images, synthesized using MIMECS applied to the corresponding MPRAGE images in the left column. | ||
|- | |- | ||
|} | |} | ||
| − | + | |} | |
| + | |} | ||
| + | </div> | ||
To understand the work on this paper you must understand the idea of sparse priors. Sparse signal reconstruction maintains that because most signals are sparse in some way it is not necessary to observe the full signal in order to accurately reconstruct it. Sparse signal reconstruction has been applied to many image processing algorithms including denoising, image restoration, and super-resolution. Dictionary-based patch matching methods, popular for natural image restoration, are also used. In our case, we build an overcomplete dictionary using patches from human brain images to both normalize a given contrast and synthesize alternate MR tissue contrasts. | To understand the work on this paper you must understand the idea of sparse priors. Sparse signal reconstruction maintains that because most signals are sparse in some way it is not necessary to observe the full signal in order to accurately reconstruct it. Sparse signal reconstruction has been applied to many image processing algorithms including denoising, image restoration, and super-resolution. Dictionary-based patch matching methods, popular for natural image restoration, are also used. In our case, we build an overcomplete dictionary using patches from human brain images to both normalize a given contrast and synthesize alternate MR tissue contrasts. | ||
The MIMECS method is based on analysis of image patches. Patches are typically small and centered on the voxel of interest. The voxels within a patch are always ordered in the same way, using a consistent rasterization. For each voxel, a dictionary is found by casting all patches into one higher dimension and searching for “like” patches. A sparse and nonnegative combination of dictionary patches is found to match the subject patch. Then the patches are combined with patches from another dictionary physically matched to the first dictionary. MIMECS also shares strong similarities with nonlocal means methods, some which have been applied in the medical imaging community for noise reduction, super-resolution, and segmentation. An illustration of how the MIMEC Algorithm works is detailed in Fig. 1. | The MIMECS method is based on analysis of image patches. Patches are typically small and centered on the voxel of interest. The voxels within a patch are always ordered in the same way, using a consistent rasterization. For each voxel, a dictionary is found by casting all patches into one higher dimension and searching for “like” patches. A sparse and nonnegative combination of dictionary patches is found to match the subject patch. Then the patches are combined with patches from another dictionary physically matched to the first dictionary. MIMECS also shares strong similarities with nonlocal means methods, some which have been applied in the medical imaging community for noise reduction, super-resolution, and segmentation. An illustration of how the MIMEC Algorithm works is detailed in Fig. 1. | ||
| − | |||
| − | ''Figure 3' | + | |
| − | + | <div style="background: white; border: 2px solid rgb(150, 150, 150); padding: 2px; text-align: justify;"> | |
| − | ''Illustration of the MIMECS algorithm. Region of the image labeled A, shows the construction of the atlases [[Image:MIMECS_A1.png]] and [[Image:MIMECS_A2.png]]. Region B shows the input subject images converted into patches [[Image:MIMECS_b1(j).png ]]. Region C denotes the estimation of the coefficients [[Image:MIMECS_x(j).png]] which relate [[Image:MIMECS_b1(j).png ]] to the patches of the atlas [[Image:MIMECS_A1.png]]. Finally, region D shows the computation of the synthetic image [[Image:MIMECS_Sn_1.png]] by using the learned coefficients [[Image:MIMECS_x(j).png]] to compute the patch [[Image:MIMECS_b2(j).png]] based on the atlas [[Image:MIMECS_A2.png]]. | + | {| |
| − | + | |- | |
| + | | align=center | [[Image:MIMECS_Figure_1.png]] | ||
| + | |- | ||
| + | | '''Figure 3:''' Illustration of the MIMECS algorithm. Region of the image labeled A, shows the construction of the atlases [[Image:MIMECS_A1.png]] and [[Image:MIMECS_A2.png]]. Region B shows the input subject images converted into patches [[Image:MIMECS_b1(j).png ]]. Region C denotes the estimation of the coefficients [[Image:MIMECS_x(j).png]] which relate [[Image:MIMECS_b1(j).png ]] to the patches of the atlas [[Image:MIMECS_A1.png]]. Finally, region D shows the computation of the synthetic image [[Image:MIMECS_Sn_1.png]] by using the learned coefficients [[Image:MIMECS_x(j).png]] to compute the patch [[Image:MIMECS_b2(j).png]] based on the atlas [[Image:MIMECS_A2.png]]. | ||
| + | |- | ||
| + | |} | ||
| + | </div> | ||
{{h3|Results}} | {{h3|Results}} | ||
| − | The most important results from the set of experiments can be summarized as follows. First, we found that features an nose that are present in the subject image but not in the atlas image (defined in the paper) are suppressed to some degree in the reconstruction. This was true in both recovering the same tissue contrast and in recovering different tissue contrast. This suggests that MIMECS can be used for noise suppression and also shows that MIMECS atlas should contain a patch dictionary that is rich enough to encompass the patches that are expected to be found in the source image. Second, the experiments confirmed that sparsity is important when reconstructing patches with significant detail because otherwise the detail will be smoothed out. Third, MIMECS is robust to the specifics of the subject pulse sequence. In | + | The most important results from the set of experiments can be summarized as follows. First, we found that features an nose that are present in the subject image but not in the atlas image (defined in the paper) are suppressed to some degree in the reconstruction. This was true in both recovering the same tissue contrast and in recovering different tissue contrast. This suggests that MIMECS can be used for noise suppression and also shows that MIMECS atlas should contain a patch dictionary that is rich enough to encompass the patches that are expected to be found in the source image. Second, the experiments confirmed that sparsity is important when reconstructing patches with significant detail because otherwise the detail will be smoothed out. Third, MIMECS is robust to the specifics of the subject pulse sequence. In particular, we found that an atlas comprising both an SPGR and MPRAGE can be used to synthesize either an MPRAGE or SPGR image even if the subject image did not precisely match the SPGR or MPRAGE pulse sequence parameters of the atlas. Fourth, we found that MIMECS normalization is comparable or better than histogram matching when normalizing the same tissue contrast and is far more general in that it can recover missing contrasts. Finally, we found we found that errors in MIMECS synthesis do not take on specific anatomical features, but instead resemble a typical noise pattern of MR images. |
| − | |||
| − | ''Figure 4' | + | |
| − | + | <div style="background: white; border: 2px solid rgb(150, 150, 150); padding: 2px; text-align: justify;"> | |
| − | ''Top row shows a subject T1-w SPGR image, three of the five classification atlases, and fuzzy memberships produced by MIMECS-based tissue classification. Hard segmentations from two leading automatic methods and a manually-corrected method are compared to the hard segmentation of MIMECS. | + | {| |
| − | + | |- | |
| + | | align=center | [[Image:MIMECS_Figure_4.png]] | ||
| + | |- | ||
| + | | '''Figure 4:''' Top row shows a subject T1-w SPGR image, three of the five classification atlases, and fuzzy memberships produced by MIMECS-based tissue classification. Hard segmentations from two leading automatic methods and a manually-corrected method are compared to the hard segmentation of MIMECS. | ||
| + | |- | ||
| + | |} | ||
| + | </div> | ||
| Line 53: | Line 74: | ||
{{h3|Acknowledgements}} | {{h3|Acknowledgements}} | ||
| − | The authors would like to thank | + | The authors would like to thank [[Amod|Amod Jog]] for his help in coding the K-d tree and to D.L. Pham for his insightful comments about this manuscript during its preparation. The authors would also like to thank the BLSA participants and neuroimaging staff for their dedication to these studies. Some data used for this study were downloaded from the Biomedical Informatics Research Network (BIRN) Data Repository, supported by grants to the BIRN Coordinating Center (U24-RR019701) Function BIRN (U24-RR021992, Morphometry BIRN (U24-RR021382), and Mouse BIRN (U24-RR021760) Testbeds, funded by the National Center for Research Resources at the National Institutes of Health, USA |
| − | {{h3| | + | {{h3|Software}} |
| − | + | MIMECS software can be downloaded from the [http://www.nitrc.org/projects/image_synthesis/ Image Synthesis] project page at NITRC. | |
| − | |||
| − | + | {{h3|Publications}} | |
| + | {{iacl-pub|author=S. Roy, A. Carass, and J.L. Prince|title=Magnetic Resonance Image Example Based Contrast Synthesis|jrnl=tmi|number=32(12):2348-2363|when=2013|doi=10.1109/TMI.2013.2282126}} | ||
| + | {{iacl-pub|author=S. Roy, A. Carass, and J.L. Prince|title=Patch Based Intensity Normalization of Brain MR Images|conf=isbi2013}} | ||
| + | {{iacl-pub|author=S. Roy, A. Carass, and J.L. Prince|title=A Compressed Sensing Approach For MR Tissue Contrast Synthesis|conf=ipmi2011|doi=10.1007/978-3-642-22092-0_31|pubmed=21761671|pmcid=3398746}} | ||
| + | {{iacl-pub|author=S. Roy, A. Carass, and J.L. Prince|title=Compressed Sensing Based Intensity Non-Uniformity Correction|conf=isbi2011|doi=10.1109/ISBI.2011.5872364}} | ||
| + | {{iacl-pub|author=S. Roy, A. Carass, P.-L. Bazin, and J.L. Prince|title=Intensity Inhomogeneity Correction of Magnetic Resonance Images using Patches|conf=spie2011|doi=10.1117/12.877466}} | ||
Latest revision as of 00:44, 3 July 2022
Magnetic Resonance Image Example-Based Contrast Synthesis
Snehashis Roy, Aaron Carass, and Jerry L Prince
Introduction
Magnetic resonance (MR) imaging (MRI) is widely used to image the brain. Often times though after the use of an image analysis algorithm the image quality can deteriorate if the image is not acquired according to optimal protocol. There are many practical reasons to disregard protocol such as when making a rapid emergency scan or multi-institutional trials, but it often results with images that are of not the desired tissue contrast or intensity. This paper proposes an image restoration technique called MR image example-based contrast synthesis (MIMECS). The central idea behind MIMECS is to synthesize MR contrasts using patches (small regions of the image). The core idea of MIMECS is that each subject patch can be matched to a combination of a few relevant patches from a dictionary in a nonlocal fashion [43], [45], [46]. It is different from classical histogram matching in the sense that a patch can be thought of as a feature vector of its center voxel, thus including local neighborhood information for that voxel. There are several advantages of MIMECS compared to previous synthesis methods. First, it is a completely automatic preprocessing step that can precede any image processing task. Second, there is no need to use multiple, calibrated pulse sequences in order to estimate underlying tissue parameters. Third, even though MIMECS uses an atlas, there is no need to carry out subject-to-atlas registration, thus avoiding potential errors due to misregistration or missing tissues. Finally, it does not require any segmentation of the subject.
Method
|
|
To understand the work on this paper you must understand the idea of sparse priors. Sparse signal reconstruction maintains that because most signals are sparse in some way it is not necessary to observe the full signal in order to accurately reconstruct it. Sparse signal reconstruction has been applied to many image processing algorithms including denoising, image restoration, and super-resolution. Dictionary-based patch matching methods, popular for natural image restoration, are also used. In our case, we build an overcomplete dictionary using patches from human brain images to both normalize a given contrast and synthesize alternate MR tissue contrasts.
The MIMECS method is based on analysis of image patches. Patches are typically small and centered on the voxel of interest. The voxels within a patch are always ordered in the same way, using a consistent rasterization. For each voxel, a dictionary is found by casting all patches into one higher dimension and searching for “like” patches. A sparse and nonnegative combination of dictionary patches is found to match the subject patch. Then the patches are combined with patches from another dictionary physically matched to the first dictionary. MIMECS also shares strong similarities with nonlocal means methods, some which have been applied in the medical imaging community for noise reduction, super-resolution, and segmentation. An illustration of how the MIMEC Algorithm works is detailed in Fig. 1.
|
Figure 3: Illustration of the MIMECS algorithm. Region of the image labeled A, shows the construction of the atlases  and  . Region B shows the input subject images converted into patches .png){kind=small} . Region C denotes the estimation of the coefficients .png){kind=small} which relate to the patches of the atlas . Finally, region D shows the computation of the synthetic image  by using the learned coefficients to compute the patch .png){kind=small} based on the atlas .
|
Results
The most important results from the set of experiments can be summarized as follows. First, we found that features an nose that are present in the subject image but not in the atlas image (defined in the paper) are suppressed to some degree in the reconstruction. This was true in both recovering the same tissue contrast and in recovering different tissue contrast. This suggests that MIMECS can be used for noise suppression and also shows that MIMECS atlas should contain a patch dictionary that is rich enough to encompass the patches that are expected to be found in the source image. Second, the experiments confirmed that sparsity is important when reconstructing patches with significant detail because otherwise the detail will be smoothed out. Third, MIMECS is robust to the specifics of the subject pulse sequence. In particular, we found that an atlas comprising both an SPGR and MPRAGE can be used to synthesize either an MPRAGE or SPGR image even if the subject image did not precisely match the SPGR or MPRAGE pulse sequence parameters of the atlas. Fourth, we found that MIMECS normalization is comparable or better than histogram matching when normalizing the same tissue contrast and is far more general in that it can recover missing contrasts. Finally, we found we found that errors in MIMECS synthesis do not take on specific anatomical features, but instead resemble a typical noise pattern of MR images.
|
| Figure 4: Top row shows a subject T1-w SPGR image, three of the five classification atlases, and fuzzy memberships produced by MIMECS-based tissue classification. Hard segmentations from two leading automatic methods and a manually-corrected method are compared to the hard segmentation of MIMECS. |
Conclusion
The results demonstrate the potential application of MIMECS, particularly in solving the “missing” or “mis-matched” tissue contrast problem in routine neuroimage processing tasks. The MIMECS framework shares similarities to other patch-based methods, but offers direct synthesis of an alternate tissue contrast. The case study demonstrates the immediate utility of MIMECS for a wide variety of neuroimage processing tasks. MIMECS represents a new class of neuroimage processing methods with a potentially rich future.
Acknowledgements
The authors would like to thank Amod Jog for his help in coding the K-d tree and to D.L. Pham for his insightful comments about this manuscript during its preparation. The authors would also like to thank the BLSA participants and neuroimaging staff for their dedication to these studies. Some data used for this study were downloaded from the Biomedical Informatics Research Network (BIRN) Data Repository, supported by grants to the BIRN Coordinating Center (U24-RR019701) Function BIRN (U24-RR021992, Morphometry BIRN (U24-RR021382), and Mouse BIRN (U24-RR021760) Testbeds, funded by the National Center for Research Resources at the National Institutes of Health, USA
Software
MIMECS software can be downloaded from the Image Synthesis project page at NITRC.
Publications
- S. Roy, A. Carass, and J.L. Prince, "Magnetic Resonance Image Example Based Contrast Synthesis", IEEE Trans. on Medical Imaging, 32(12):2348-2363, 2013. (doi)
- S. Roy, A. Carass, and J.L. Prince, "Patch Based Intensity Normalization of Brain MR Images", Tenth IEEE International Symposium on Biomedical Imaging (ISBI 2013), San Francisco, CA, April 7 - 11, 2013.
- S. Roy, A. Carass, and J.L. Prince, "A Compressed Sensing Approach For MR Tissue Contrast Synthesis", 22nd Conference on Information Processing in Medical Imaging (IPMI 2011), Monastery Irsee, Germany, July 3-8, 2011. (doi) (PubMed) (PMCID 3398746)
- S. Roy, A. Carass, and J.L. Prince, "Compressed Sensing Based Intensity Non-Uniformity Correction", Eighth IEEE International Symposium on Biomedical Imaging (ISBI 2011), Chicago, IL, March 30 - April 2, 2011. (doi)
- S. Roy, A. Carass, P.-L. Bazin, and J.L. Prince, "Intensity Inhomogeneity Correction of Magnetic Resonance Images using Patches", Proceedings of SPIE Medical Imaging (SPIE-MI 2011), Orlando, FL, February 12-17, 2011. (doi)
