Convert 2D siemens mosaic MRI image to multi slices

I have a siemens mosaic MRI image, and it's 512*512 unit-16 matrix after I use dicomread to read it. It contains 60 slices and the size for each slice is 64*64. I want to convert the montage image to slices. Is there some toolbox to do that? Thanks!

Answer 1

[Walter Roberson]

slices = mat2cell(MosaicArray, 64 * ones(1,size(MosaicArray,1)/64), 64 * ones(1,size(MosaicArray,2)/64) );

Now you would need to get rid of 4 of the items, but at the moment I do not know which 4 are the empty ones and I do not know whether the unraveling should be along the rows or down the columns. Probably it should be

slices = transpose(slices);    %rows become priority
slices(61:end) = [];   %remove the 4 empty slices

[Hui WANG]

It's right! Thank you very much!

[Makis]

I am having the same problem. I have some mosaic 2D dicom data with dimensions 1024x1024x1 and I want to write them as multi-slice mosiac dicom data with dimensions 128x128x64.

After using:

slices = mat2cell(MosaicArray, 128 * ones(1,size(MosaicArray,1)/128), 128 * ones(1,size(MosaicArray,2)/128) );

how can I save these slices into the final desired mosaic dicom format?