We introduce an optimized six-minute high-resolution whole-brain protocol of 3D, mainly isotropic acquisitions for T2w SPACE, FLAIR (3D), T1w MPRAGE and T2*w SWI; all acquired at 9-fold acceleration using Wave-CAIPI encoding. Design strategies were developed to optimize wave-encoding parameter selection, coil sensitivity computation, and data-driven automated gradient trajectory error estimation. We validated our exam on five volunteers on a 3T MRI scanner and achieved good SNR with negligible g-factor noise and good image contrast despite the 9-fold acceleration. Our fast high-resolution neurological exam should benefit several clinical and research applications and help mitigate patient compliance issues and motion artifacts.
The Wave-CAIPI acquisition scheme was implemented into T2w SPACE, FLAIR (3D), T1w MPRAGE, and T2*w SWI prototype sequences. Auto-calibrated reconstruction [4] was implemented for all sequences to correct for wave-gradient trajectory errors during the image reconstruction.
To ensure robust, high-quality imaging, the following design strategies were investigated:
In compliance with IRB requirements, five healthy volunteers were scanned on a 3T MRI Scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany) using a 32-channel head coil.
Fig.1 shows results from a six-minute Wave protocol optimized based on findings from Fig.3-5 (see details below). Detailed acquisition parameters are provided in Fig.2. Despite high acceleration (R=3x3), all scans achieve good SNR and contrast, which is reflected in the low g-factor maps (see Fig.1). Here, gavg is close to 1 for all acquisitions, with minor hotspots in the bottom part of the brain. Hotspots are highest for the SPACE sequences (T2w and FLAIR) where larger bandwidth (592 Hz/px) was used in comparison to MPRAGE (200 Hz/px) and SWI (100 Hz/px). This limits the wave-spreading and g-factor performance.
Fig.3a illustrates that the extent of Wave spreading (i.e. controlled-aliasing along x, which correlates to g-factor performance) is mainly determined by the Wave’s gradient amplitude and not the number of wave cycles. However, using a high #cycles will reduce the Wave’s k-space radius for a given amplitude, which results in much less mixing of signal from different partition/phase-encoding steps experiencing different signal modulation due to T1 recovery or T2 decay (e.g. in MPRAGE, SPACE). This in turn leads to a sharp point-spread-function with very small side-lobes (Fig.3b) and reduced artifacts as shown in Fig.4b,c.
In Fig.3c the wave frequency of the shown PSF does not exactly agree with the ADC duration due to differences in raster times. This removes the sparse nature of the PSF in frequency space as additional side-lobes appear besides the center frequency and DC peak, which are not estimated in the auto-calibrated reconstruction. We mitigated associated artefacts by enforcing the Wave cycle duration to be an exact integer multiple of the ADC duration, as experimentally validated in Fig.4a,b.
We found that ESPIRiT can robustly obtain accurate coil sensitivities using just 24x24 k-space lines of our 3D-GRE. To mitigate potential artifacts from motion during the exam the fast pre-scan (2s) was performed prior to each Wave acquisition.
Typically, missing coil sensitivity information from the shoulder causes FOV wraps in the corresponding phase/partition-encoding plane. However, for Wave, such wraps can cause artefacts which spread along the readout (head-food) direction (due to incorrect deblurring) into the area of interest, e.g. the brain (see Fig.5a,b). The use of a second set of sensitivities from ESPIRiT was able to mitigate this issue which has been incorporated into all reconstructions.
This work was supported in part by NIH research grants: R01EB020613, R01EB019437, R24MH106096, P41EB015896, and the shared instrumentation grants: S10RR023401, S10RR019307, S10RR019254, S10RR023043