Solution #3 - MRI of healthy human man

[MATLAB logo] Exploiting MATLAB's Data Structures

Solution #3 - MRI of healthy human man

An MRI scan is an good example of a volume dataset. Here the density of tissues is recorded to create a series of images along slicing planes. In this problem you will learn how to manipulate N-D arrays and display the results of slicing the dataset along various directions.

Goal: Display a multidimensional dataset two ways. Display a single slice as an image and use the slice function to visualize the data in 3-D.

Key concepts: Multidimensional arrays, N-D indexing, slicing.

Key functions: Multidimensional indexing, slice, reshape, whos, size

Solution

The file mri.mat (from the Image Processing Toolbox) contains a 4-D dataset that represents an MRI scan of a healthy male. It is stored as 27 grayscale images. Load the dataset and determine its size with whos or size.
load mri
whos
size(D)
Display the 17th slice as an image.
image(D(:,:,:,17)
colormap(gray)
slice is a volume visualization function that displays slices of a 3-D dataset as planes with the correct orientation. Use slice to display the 17th slice of the dataset as well as slices along the 70th column (x-axis) and 40 and 85th rows (y-axis). Since slice only works on 3-D arrays, you may need to use double and reshape first.
D = double(D);
D = reshape(D,128,128,27);
slice(D,70,[40 85],17)
shading flat
view(-37,45)

[hint] [Show me the solution]