2156

Investigating the Accuracy of PD Measurements at 3T
Zach Lun Pang1, Maria Schmidt2, Evanthia Kousi2, and Julie Hughes3
1Medical Physics, St George’s University Hospitals NHS Foundation Trust, London, United Kingdom, 2Joint Department of Physics, The Institute of Cancer Research and The Royal Marsden NHS Foundation Trust, London, United Kingdom, 3The Royal Marsden NHS Foundation Trust, London, United Kingdom

Synopsis

Keywords: Phantoms, Quantitative Imaging

We investigated the accuracy of clinical proton density (PD) mapping employing the variable flip-angle (VFA) method at 3T, considering the effects of B1+ and B1- fields. Following B1+ correction, the PD measurements were within 5% of the expected value for flip-angle variations up to 18%. Considerable discrepancies were observed for T1<300ms and flip-angle variations >25%. Using image intensity filters, considerable PD variations were found at the edges of the test object. Clinically, the B1+ and B1- corrected PD values of the brain were similar to values reported in the literature, rendering the VFA method a viable tool for PD mapping.

Introduction

Quantitative Magnetic Resonance Imaging (QMRI) methods are increasingly used to study tissue properties non-invasively, mapping relaxation times and proton density (PD) to aid diagnoses and assess response to treatment1. Different methods are available, but long Acquisition Times (Tacq) are impractical2. The Variable Flip Angle (VFA) method offers high resolution imaging with short Tacq3 to map both PD and T14 but requires accurate knowledge of the actual flip angle. Hence, it is vulnerable to errors in the transmitted RF field (B1+), which are significant at 3T due to the conductivity of the subject5.

Objective: To determine whether the VFA method can be used clinically to produce accurate PD maps at 3T when used in conjunction with B1+ mapping tools and post-processing to remove B1+ and receive RF field (B1-) effects. Differences in receiver coil loading were also investigated while focusing on brain tissue mapping.

Methods

Work was undertaken at 3T (Siemens MAGNETOM Vida & Skyra, Erlangen, Germany) at the Royal Marsden Hospital, using the integrated transmit body coil with a 20-channel head receiver coil.

Experiment 1: Initial experiments assessed PD and T1 measurements in a Eurospin (Diagnostic Sonar, Livingston, Scotland) Test Object (TO) with and without B1+ corrections. This required using Siemens’ B1+ mapping tool, which uses the 2D Actual Flip-angle Imaging (AFI) method and implementing it in the calculation for PD and T1. Acquisitions were repeated with and without Image Intensity Filters (IIFs) to assess whether these could correct B1- bias. PD maps were generated using in-house software (MATLAB 2019b), on a voxel-wise basis from the signal intensity and flip angle. B1+ mapping tool performance was assessed by adhering copper to the Eurospin TO, inducing strong B1+ inhomogeneities, and repeating PD measurements.

Experiment 2: For PD mapping to be useful for diagnoses or treatment assessments, individuals’ tissue measurements will need to be evaluated against a range of healthy tissue values. Therefore, it needed to be ascertained whether patient size variations impact PD measurements. Therefore, PD in a uniform TO was measured and repeated with additional TOs in order to load the transmit coil (integrated body) and receive coil (head coil) differently, simulating variations in patient sizes.

Experiment 3: Head coil sensitivity was mapped by measuring PD in an oil phantom (IIFs disabled, B1+ corrected) and assuming PD variations were solely due to B1- bias. This was applied to the previous PD experiment.

Experiment 4: PD was measured in a volunteer with research-ethics approval. VFAs of 3° & 14° were selected as these are optimal for tissues with T1 ~1000ms6, in-between White Matter (WM) and Grey Matter (GM) at 3T. Tacq for two VFA sequences and a B1+ map was ~90s. Coil sensitivity was measured prior to scanning.

Results

Experiment 1: Correcting B1+ errors improved results, but variations in coil sensitivity caused apparent reductions in PD. IIFs caused variations unlike those associated with coil sensitivity. PD was reduced in regions near the copper, indicating a limit of the B1+ mapping tool.

Experiment 2: Frequency adjustments and appropriate B0 shimming resulted in negligible differences in PD measurements with the additional loading. Variations caused by IIFs were more obvious with the uniform TO.

Experiment 3: Correcting PD measurements of B1- bias improved uniformity but was not completely effective. Possibly due to additional post-processing that affects signal intensity.

Experiment 4: B1+ and B1- corrected PD measurements normalised to the PD of CSF were generated, resulting in values like those in the literature (PDWM = 65.9±3.7%, PDGM = 80.5±8.1%). The ventricles, WM and GM were segmented using SPM127 (Functional Imaging Laboratory, Wellcome Centre for Human Neuroimaging, London, UK). Average T1 measurements (T1CSF = 5319±3614ms, T1WM = 957.8±384.1ms, T1GM = 1398±965.8ms) were also found to be like those reported in the literature, providing greater confidence in the results.

Discussion

Results show IIFs do not remove the coil sensitivity weighting completely with TOs showing unexpected variations in PD closest to their edges. Better results were achieved using a B1- map obtained in oil but is only possible in fixed coils. Additionally, B1+ correction is essential for accuracy; the B1+ mapping tool is only accurate for small variations. The remaining areas of concern are tissues with T1 < 300ms and with B1+ inhomogeneities >25%.

Though average PD and T1 values resembled those in the literature8, this study highlights the distribution of brain tissue measurements, meriting further discussion. The distribution can be attributed to noise and partial volume effects, but errors in TOs have also been demonstrated. While the VFA method offers high resolution in 3D, the AFI B1+ mapping produced by Siemens is 2D, and this may be a source of error9. For the future, a 3D AFI should be explored as an alternative.

Conclusion

Measuring PD at 3T is viable with the VFA and B1+ mapping tool. For flip angle variations of up to ~18%, PD measurements were kept within 5% of the correct value. Best results were achieved without IIFs and correcting coil sensitivity using a non-conductive (oil) TO. Unfixed coils are not recommended, and regular QA is needed since changes in coil performance will affect coil sensitivity.

Acknowledgements

This study represents independent research funded by the National Institute for Health and Care Research (NIHR) Biomedical Research Centre and the Clinical Research Facility in Imaging at The Royal Marsden NHS Foundation Trust and The Institute of Cancer Research, London. The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care.

References

1. Baudrexel S, Reitz S, Hof S, et al. Quantitative T1 and proton density mapping with direct calculation of radiofrequency coil transmit and receive profiles from two-point variable flip angle data. NMR in Biomedicine. 2016; 29(3): 349–360.

2. Neeb H, Zilles K, Shah N.J., et al. A new method for fast quantitative mapping of absolute water content in vivo. NeuroImage 2006; 31(3): 1156-1168.

3. Liberman G, Louzoun Y, Bashat D.B., et al. T-1 Mapping Using Variable Flip Angle SPGR Data With Flip Angle Correction. Journal of magnetic resonance imaging 2014; 40.

4. Wang Y, Chen Y, Wu D, et al. STrategically Acquired Gradient Echo (STAGE) imaging, part II: Correcting for RF inhomogeneities in estimating T1 and proton density. Magnetic Resonance Imaging. 2018; 46: 140–150.

5. Vaidya MV, Collins CM, Sodickson DK, et al. Dependence of B1+ and B1- Field Patterns of Surface Coils on the Electrical Properties of the Sample and the MR Operating Frequency. Concepts Magn Reson Part B Magn Reson Eng. 2016; 46(1): 25-40.

6. Wood T.C. Improved formulas for the two optimum VFA flip-angles. Magnetic Resonance in Medicine 2015; 74(1): 1–3.

7. Flandin G, Wellcome Trust Centre for Neuroimaging, University College, London, 12 Queen Square, London, WC1N 3BG

8. Bojorquez J.Z., Bricq S, Acquitter C, et al. What are normal relaxation times of tissues at 3 T? Magnetic Resonance Imaging 2017; 35: 69-80.

9. Malik S.J., Kenny G, Hajnal J.V. Slice profile correction for transmit sensitivity mapping using actual flip angle imaging. Magnetic resonance in medicine 2011; 65(5): 1393–1399.

Figures

Figure 1: Flip angle map generated using the Siemens B1+ mapping sequence. Left is a flip angle map of the test object without the copper tape, right is with the copper tape. Nominal flip angle for both maps was 800. The red arrow indicates the position of the copper tape causing a reduction in the flip angle. T1 ranges from approximately 75 – 2700ms.

Figure 2: Horizontal line profiles showing the effect of the copper in both filtered & unfiltered and B1+ corrected & uncorrected PD maps. Overlapping profiles (with and without copper) are inconsistent near to the copper (indicated by the red arrows). As the distance from the copper tape is increased, profiles become consistent. Additionally, IIFs do not completely remove the coil sensitivity bias, which is noticeable at the edges of the object (green arrows).

Figure 3: PD line profiles with and without the additional coil loading. Left plot shows variations that that cannot be explained by variations in coil sensitivity. Right plot shows variations due to coil sensitivity and closer to what is expected. Red arrow indicates an error due to inappropriate B0 shimming. Green arrow indicates wrap-around artifact. As the intensity of the artifact is reduced, this indicates the IIFs are only partially effective at mitigating coil sensitivity effects.

Figure 4: B1+,- corrected PD map of a uniform TO. Regions of the TO that do not overlap with the oil phantom used to measure coil sensitivity have been cancelled.

Figure 5: In-vivo PD measurements with the ventricles (red), WM (blue) and GM (pink) segmented. PD maps are normalised to the PD of CSF in the ventricles. Voxel size was 1.38 x 1.38mm2.

Proc. Intl. Soc. Mag. Reson. Med. 31 (2023)
2156
DOI: https://doi.org/10.58530/2023/2156