Education, tips and tricks to help you conduct better fMRI experiments.
Sure, you can try to fix it during data processing, but you're usually better off fixing the acquisition!

Friday, December 30, 2016

Use of split slice GRAPPA (aka Leak Block) for SMS-EPI reconstruction

Accurate separation of the simultaneously acquired slices is one of the bigger limitations of the SMS-EPI method, compared to the processing used for conventional multislice EPI. The default SMS reconstruction, as used in my two introductory posts on the SMS sequences from CMRR (MB-EPI) and MGH (Blipped CAIPI), is a slice dimension adaptation of the GeneRalized Autocalibrating Partial PArallel (GRAPPA) method that was originally applied in-plane to acceleration of the phase encoding direction. It's not essential to understand the GRAPPA method applied in-plane for the purposes of understanding this post or for SMS reconstruction more generally. But if you're curious I wrote a brief introduction to in-plane GRAPPA in 2011. That post was specifically concerned with motion sensitivity of (in-plane) GRAPPA. I'll be looking in more detail at the motion sensitivity of SMS in a future post. In this post I want to compare the standard SMS reconstruction - what is generally termed Slice GRAPPA - with an alternative known as Split Slice GRAPPA. The latter option is termed "Leak Block" in the CMRR pulse sequence, MB-EPI.

What's the concern?

CMRR's parameter nomenclature offers a strong clue to the problem. In conventional EPI reconstruction we use a 2D Fourier transform (FT) which produces some amount of ringing. We also use slices that have some degree of cross-talk to neighboring slices, arising out of the limitations of frequency selectivity. So, while we think of voxels as perfect little rectangles or cubes, in reality they are blurry beasts that spread their signal into adjoining voxels because of a non-rectangular point-spread function (PSF). The dimensions we assign a voxel are entirely nominal.

With SMS we have a broader spatial problem than just non-cubic PSF. Separation of the simultaneous slices can leave signal in an incorrect position that is quite some distance from where it is supposed to be. It's a longer length scale error than the simple PSF of a voxel. Let's suppose we acquire four 2 mm slices simultaneously, 84 total slices. In one SMS acquisition we will have four slices separated by one quarter of the total slice dimension extent of 168 mm, or about 42 mm (assuming no additional inter-slice gap). Do a quick thought experiment. Imagine that in the first slice there is a very strong activation and nothing in the other three. If there is a large residual spatial error arising from poor SMS separation then we might start seeing this activation projected 4.2, 8.4 or even 12.6 cm from where it should be! And how would we know that the distant activation sites were erroneous?

This slice leakage, as it's usually called in the literature, may be strongest for simultaneously acquired neighbors but may extend throughout the slice dimension, between simultaneously acquired slices that might be quite far apart in anatomical space. And, as the thought experiment illustrates, one might assume that distant leakage would be harder to spot than the conventional cross-talk between successively acquired slices in conventional multislice EPI, or errors arising from the PSF more generally. The PSF can usually be interpreted as a local phenomenon, with errors decreasing monotonically from a voxel. Not so with SMS slice separation, meaning there is more risk of interpreting a false positive remote from the true activation site.

At this point we can recognize that reducing leakage is a noble, perhaps essential, goal. As usual with MRI, however, there's a catch. Reducing leakage using the Split Slice GRAPPA reconstruction may come at the cost of increasing in-plane artifacts. The overall (total) artifact level might be higher, too. I'll go into these issues in some detail below. The goal of this post is to perform a rudimentary assessment of the artifacts and determine the circumstances when Split Slice GRAPPA might be preferred over the conventional Slice GRAPPA reconstruction. For the CMRR sequence this amounts to whether or not to enable the Leak Block option.

Thursday, October 13, 2016

Motion traces for the respiratory oscillations in EPI and SMS-EPI

This is a follow-up post to Respiratory oscillations in EPI and SMS-EPI. Thanks to Jo Etzel at WashU, you may view here the apparent head motion reported by the realignment algorithm in SPM12 for the experiments described in the previous post. Each time series is 200 volumes long, TR=1000 ms per volume. The realignment algorithm uses the first volume in each series as the template. The motion is plotted in the laboratory frame, where Z is the magnet bore axis (head-to-foot for a supine subject), X is left-right and Y is anterior-posterior for a supine subject.

In the last post I said that there were five total episodes of a deep breath followed by sigh-like exhale, but actually the subject produced a breath-exhale on average every 30 seconds throughout the runs. (This was a self-paced exercise.) Thus, what you see below (and in the prior post) has a rather large degree of behavioral variability. Still, the main points I made previously are confirmed in the motion traces. I'll begin with the axial scan comparison. Here are the motion parameters for the MB=6 axial acquisition with standard foam head restraint (left) versus the custom printed restraint (right):

MB=6, axial slices. Left: foam restraint. Right: custom 3D printed headcase restraint

The effect of the custom restraint is quite clear. The deep breath-then-sigh episodes are especially apparent when using only foam restraint. Note the rather similar appearance of the high frequency oscillations, particularly apparent in the blue (Y axis) traces between the two restraint systems, suggesting that the origin of these fluctuations is B0 modulation from chest motion rather than direct mechanical motion of the head. We cannot yet be sure of this explanation, however, and I am keeping an open mind just in case there are small movements that the custom head restraint doesn't fix.

Friday, October 7, 2016

Respiratory oscillations in EPI and SMS-EPI

tl;dr   When using SMS there is a tendency to acquire smaller voxels as well as use shorter TR. There are three mechanisms contributing to the visibility of respiratory motion with SMS-EPI compared to conventional EPI. Firstly, smaller voxels exhibit higher apparent motion sensitivity than larger voxels. What was intra-voxel motion becomes inter-voxel motion, and you see/detect it. Secondly, higher in-plane resolution means greater distortion via the extended EPI readout echo train, and therefore greater sensitivity to changes in B0. Finally, shorter TR tends to enhance the fine structure in motion parameters, often revealing oscillations that were smoothed at longer TR. Hence, it's not the SMS method itself but the voxel dimensions, in-plane EPI parameters and TR that are driving the apparent sensitivity to respiration. Similar respiration sensitivity is obtained with conventional single-shot EPI as for SMS-EPI when spatial and temporal parameters are matched.


The effects of chest motion on the main magnetic field, B0, are well-known. Even so, I was somewhat surprised when I began receiving reports of likely respiratory oscillations in simultaneous multi-slice (SMS) EPI data acquired across a number of projects, centers and scanner manufacturers. (See Note 1.) Was it simply a case of a new method getting extra attention, revealing an issue that had been present but largely overlooked in regular EPI scans? Or was the SMS scheme exhibiting a new, or exacerbated, problem?

Upper section of Fig. 4 from Power,, showing the relationship between apparent head motion (red trace) reported from a realignment algorithm and chest motion (blue trace) recorded by a respiratory belt. See the paper for an explanation of the bottom B&W panel.

Wednesday, June 29, 2016

Starting points for SMS-EPI at 3 T: Part II

In an earlier post I presented three starting protocols for the CMRR version of SMS-EPI, referred to as the MB-EPI sequence here. I'll use italics to indicate a specific pulse sequence whereas SMS-EPI, no italics, refers to the family of simultaneous multi-slice methods. In this post I'll develop a similar set of three starting protocols for the Massachusetts General Hospital (MGH) version of SMS-EPI, called Blipped-CAIPI. I'm going to build upon the explanations of the last post so please cross reference for parameter explanations and background.

As for the previous post there are several things to bear in mind. This series is Siemens-centric, specifically Trio-centric. While many of the concepts and parameter options may apply to other platforms there will be minor differences in parameter naming conventions and, perhaps, major differences in implementation that you will need to consider before you proceed. For Siemens users, I am running aging software, syngoMR version B17A. The age of the software and the old reconstruction board on the scanner means that you can expect to see much faster reconstruction on a newer system. I hope, but cannot guarantee, that the actual image quality and artifact level won't differ massively from a Trio running VB17A to a new Prisma running VE11C. I'll keep you updated as I learn more.


As before, for this post I am going to be using a 32-channel receive-only head coil. The SMS-EPI sequences can be made to work with a 12-channel coil but only in a reduced fashion because the 12-channel coil has minimal receive field heterogeneity along the magnet z axis - the struts run parallel with the magnet axis except at the coil's rear, where they converge - and generally we want to do axial slices (along z) for fMRI. I don't yet know whether SMS-EPI would work well on the 20-channel head/neck coil on a Prisma, it's something I hope to investigate in the near future. But a 64-channel head/neck coil on a Prisma will definitely work for SMS-EPI. Better or worse than a 32-channel coil on a Prisma? I have no idea yet.

The Blipped-CAIPI sequence version 2.2 was obtained through a C2P (Core Competence Partnership) with MGH. Installation was a breeze: a single executable to port to the scanner and one click, done. The development team offers an informative but brief 7-page manual which will be useful to anyone who has read the SMS-EPI literature and has a basic understanding of how SMS works. It's not a starting point for everyday neuroscience, however. The manual mentions a .edx (protocol) file as a starting point for 2, 2.5 and 3 mm resolution scans, but in the file I downloaded for VB17A the contents didn't include it. Perhaps contact MGH if you are on another software version and you'd like a .edx file rather than building your own protocol, e.g. by recreating what you see here.

Tuesday, June 28, 2016

Exploiting Tanzania

So a massive helium reserve may have been found in Tanzania's Rift Valley. Wonderful. All the western headlines this morning have put a typically western spin on it. Hurrah! We are saved! We get to go plunder a foreign place again for what we need to save our own lives! Before we get too carried away with ourselves, let's take a few seconds to think about a few things. Like, say, how many MRI scanners are there in Tanzania right now? How many Tanzanian lives will be saved? Anyone care to estimate? This scanner in Dar es Salaam makes headlines when it breaks down.

What about the wildlife in Tanzania? Will lives be saved there, too? Note the concentration of national parks and game reserves in and around the Rukwa region of Tanzania. Now, I'm not intimately familiar with how helium gas is extracted, concentrated or liquefied but I'm going to guess that some of it has to be done where the gas is found. Even if the gas doesn't just float conveniently into collection chambers instead of needing some sort of gas forcing process (We love fracking, right?) and miles and miles of pipelines, it's a fair assumption that there will be massive energy needs to liquefy it. Then the cryogenic liquid helium must be transported. So we'll need roads, maybe an airport for the suits to get in and out quickly, and perhaps a railway to move the product to a sea port. Or we could just push the gas down a long pipe to the coast where it could be liquefied, then transported abroad. This is all going to be great news for African nature, I'm sure of it!

Tuesday, February 9, 2016

Starting points for SMS-EPI at 3 T

Several people have approached me for advice on using simultaneous multi-slice (SMS) EPI for fMRI experiments. This is the sequence also known as multiband (MB) EPI. I'll come back to nomenclature in a moment. First, though, a brief introduction to what may become a lengthy series of posts. I'm going to focus on BOLD-based fMRI exclusively for the time being - sorry diffusion and ASL folks - and because I presently only have a Siemens Trio at my disposal, everything I write will have strong bias in that direction. That said, I do anticipate writing later posts dealing with SMS-EPI (for fMRI) on a Siemens Prisma at least, and I can already envisage a need for posts dealing with receive field normalization, in-plane parallel imaging, distortion correction options, reconstruction options and multi-echo SMS, to name just a few advanced topics. But first things first - to get going!

Options for SMS-EPI on a Siemens 3 T scanner

I am aware of three SMS-EPI pulse sequences for a Siemens Trio. One comes from the University of Minnesota's Center for Magnetic Resonance Research (hereafter CMRR), one comes from the Martinos Center for Biomedical Imaging at Massachusetts General Hospital (hereafter MGH), and one comes from Siemens as a work-in-progress (WIP) aftermarket sequence. For this post I'm going to be using the sequence provided by CMRR. Since CMRR refer to their sequence as multiband (MB) EPI I shall stick to this nomenclature here, and reserve the term SMS-EPI to apply to the broader family of pulse sequences. I may do posts on the MGH and WIP sequences in the future, but the CMRR sequence has been used the most broadly to date (e.g. the Human Connectome Project, which I'll discuss at length below) and so it offers the most immediate, road-tested place to start.

Friday, February 5, 2016

Corrective lenses for tight head coils

Most labs have plastic goggles for correcting a subject's vision in the MRI. Here's our box of tricks:

These goggles work pretty well for a standard sized head coil, such as the 12-channel TIM coil on my Siemens Trio. But for a tighter fitting coil, such as the 32-channel head coil, there is simply no way to cram a subject wearing goggles into the space available. For a start the goggles' frame prevents the subject's nose from penetrating the appropriate gap in the front of the coil.

A simple solution is to relocate the corrective lenses on the outer surface of the head coil. All that's required is a different way to hold the lenses in place. My ace engineer, Rick, made this pair of holders for our 32-channel coil: