import pylab as pl
// Compute the mean EPI: we do the mean along the axis 3, which is time
mean_img = np.mean(fmri_data, axis=3)
// pl.figure() creates a new figure
pl.figure()
// First subplot: coronal view
// subplot: 1 line, 3 columns and use the first subplot
pl.subplot(1, 3, 1)
// Turn off the axes, we don"t need it
pl.axis("off")
// We use pl.imshow to display an image, and use a "gray" colormap
// we also use np.rot90 to rotate the image
pl.imshow(np.rot90(mean_img[:, 32, :]), interpolation="nearest",
cmap=pl.cm.gray)
pl.title("Coronal")