Converting Siemens MOSAIC
[This article was first published on Rigorous Analytics, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Siemens multi-slice EPI data may be collected as a “mosiac” image; i.e., all slices acquired in a single TR (repitition time) of a functional MRI run are stored in a single DICOM file. The images are stored in an MxN array of images. The function create3D() will try to guess the number of images embedded within the single DICOM file using the AcquisitionMatrix field. If this doesn’t work, one may enter the (M,N) doublet explicitly. Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
The first image is taken from the raw DICOM file, in mosaic format, and displayed with the default margins [(5,4,4,2)+0.1].
The second image is displayed after re-organizing the original DICOM file into a 3D array and then converting it to NIfTI format. Additional information, using the valid NIfTI header fields descrip and aux_file, have been defined for these data.
Remember, all other Siemens DICOM data are stored as a seperate file for each image (as is the case with the other major MR scanner manufacturers).
To leave a comment for the author, please follow the link and comment on their blog: Rigorous Analytics.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.