Coil Compression for Improved Phase Image Signal-to-Noise Ratio in Electrical Property Tomography
Kathleen M Ropella1 and Douglas C Noll1

1Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States

Synopsis

The use of multi-channel receivers is essential for acquiring B1+ with sufficient SNR to calculate electrical properties. Combining the individual channel images prior to these calculations typically involves a SENSE-like method or the use of some reference image. In this work we present a modified version of coil compression to provide an automatic and simplified multi-channel array data combination for high SNR phase-based conductivity mapping.

Purpose

Electrical properties tomography (EPT) is a method to calculate the conductivity and permittivity of a material from measured $$$B_1^+$$$ fields using MRI. These calculations result in high noise amplification, warranting the use of high signal-to-noise ratio (SNR) images. Multi-channel receivers can be used to improve SNR of the $$$B_1^+$$$ images, but the individual coil data must be combined while leaving the image phase intact. Combining coil data prior to EPT calculations provides the benefit of eliminating phase artifacts in areas of low signal, which would lead to highly inaccurate electrical properties. Current approaches to combining phase data include SENSE-like methods1 and methods using some sort of reference coil2,3. We present an adaptation of coil compression4,5 to optimally combine the $$$B_1^+$$$ data such that the phase is preserved for EPT. This method requires no reference images or sensitivity maps. We demonstrate the SNR gains and reconstruction efficacy with phase-based conductivity mapping6: $$\sigma = \frac{1}{\omega \mu_0}\nabla^2\phi^+$$ where σ is conductivity, ω is the scanner frequency, µ0 is vacuum permeability, and φ+ is the phase of $$$B_1^+$$$.

Methods

Multi-channel array data was compressed into one virtual coil using the singular value decomposition (SVD), which provides the globally optimal combination of k-space data. The procedure for a fully-sampled slice with dimensions [nx, ny] from a N-channel array is as follows:

1. Perform a 2D IFFT for each coil image and place into a matrix X, which is size [N, nx*ny], such that each row is the vectorized k-space data from a different coil.

2. Calculate the SVD: X = UΣVH. Keep only the first column of U, the first element of Σ, and the first row of VH. Multiply these truncated matrices to form the virtual coil k-space, XV1.

A spatially varying combination can be calculated in a similar manner. This requires an alignment of the singular vectors to ensure the compressed image is smooth, as described by Zhang5. This local compression technique is called geometric decomposition coil compression with alignment, GCC-Align in this work.

The phase from each compressed coil was divided by two to approximate the transmit phase for conductivity mapping, calculated using a model-based method7.

Data was acquired for a saline phantom and a healthy control subject’s brain, each with an 8-channel receiver (GE Healthcare, Waukesha, WI) and a 32-channel receiver (Nova Medical) on a GE Discovery MR750 3.0T scanner. For the phantom, a spin echo (SE) scan was performed with TE/TR = 11/1000 ms, FOV = 24cm x 24cm, 256x256 pixels, 3mm slices. Nominal values of the phantom were measured using the coaxial probe method8. For the brain, a SE scan was performed with TE/TR = 14/1000 ms, FOV = 24cm x 24cm, 256x256 pixels, 3mm slices.

Noise levels in the virtual coils were compared in a uniform spherical phantom using the same acquisition parameters as the saline phantom. The noise variance was calculated as the variance of the residuals after fitting a parabolic surface to phase data in a 15x15 pixel sliding window.

Results

Figure 1 shows that while both compression techniques provide substantial noise reduction in the phase images, the GCC-Align method provides a more uniform noise profile across the object. Figures 2 and 3 show representative coils from the multi-channel arrays and the respective coil combinations for a saline phantom. The arrow in Figure 3 shows an open-ended fringe line in a low signal region. This problematic feature is eliminated in the virtual coil. The phantom conductivity maps are shown in Figure 4, along with the measured values. Human subject conductivity maps are shown in Figure 5.

Discussion

Compressing the SE image data to a single virtual coil decreased the noise variance by at least an order of magnitude. For a small increase in computation, we create a more uniform noise profile with the GCC-Align virtual coil. The conductivity maps shown in Figure 4 show values consistent with those measured, but the GCC-Align images have less variation farther from the coil. While true values for the brain are not available, the higher conductivity of the cerebrospinal fluid is easily identified. The GCC-Align image shows more uniformity in the white matter, as evidenced by fewer black patches, due to reduced phase noise.

Conclusion

We have applied the idea of multi-channel compression to virtual coils for the purpose of high SNR phase maps for EPT. This increases the accuracy of the reconstruction, eliminates the need for reference coils in phase combination, and can be done for the entire FOV or in a locally-varying manner.

Acknowledgements

No acknowledgement found.

References

[1] P. B. Roemer et al. The NMR phased array. Magnetic Resonance in Medicine. 1990; 16:192-225.

[2] D. O. Walsh, A. F. Gmitro, and M. W. Marcellin. Adaptive reconstruction of phased array MR imagery. Magnetic Resonance in Medicine. 2000; 43:682-690.

[3] D. L. Parker et al. Phase Reconstruction from Multiple Coil Data Using a Virtual Reference Coil. Magnetic Resonance in Medicine. 2014; 72:563-569.

[4] F. Huang et al. A software channel compression technique for faster reconstruction with many channels. Magnetic Resonance Imaging. 2008; 26:133-141.

[5] T. Zhang et al. Coil Compression for Accelerated Imaging with Cartesian Sampling. Magnetic Resonance in Medicine. 2013; 69:571-582.

[6] H. Wen. Non-invasive quantitative mapping of conductivity and dielectric distributions using the RF wave propagation effects in high field MRI. in Proceedings of SPIE: Medical Imaging; Physics of Medical Imaging. 2003; 5030:471-477.

[7] K. M. Ropella and D. C. Noll. A Regularized Model-Based Approach to Phase-Based Conductivity Mapping. in Proceedings of the 23rd Annual Meeting of ISMRM. 2015; 3295.

[8] H. Zheng and C. E. Smith. Permittivity Measurements Using a Short Open-Ended Coaxial Line Probe. IEEE Microwave and Guided Wave Letters. 1991; 1:37-339.

Figures

Figure 1: Noise comparison between a single array element, the virtual coil, and the GCC-Align virtual coil from the 32-channel array using a uniform phantom. Images are coil magnitudes and dotted lines denote the profile used in the noise comparison. Right side of phantom is closest to coil.

Figure 2: Representative magnitude and phase images from two coils of the 8-channel array and the virtual coils from the two coil compression techniques. The GCC-Align magnitude is scaled by 1/5. The inner compartment of the phantom is DI water and the outer compartment is saline solution with copper sulfate.

Figure 3: Representative magnitude and phase images from two coils of the 32-channel array and the virtual coils from the two coil compression techniques. The magnitude of the GCC-Align coil has been scaled by a factor of 1/5. The arrow denotes an open-ended fringe line, a problematic feature for reconstruction.

Figure 4: Reconstructed phantom conductivity maps using the phase of the 8-channel and 32-channel receiver virtual coils (top row), and the 8-channel and 32-channel receiver GCC-Align coils (bottom row). Conductivity measured using a coaxial probe is shown at right. The right side was closest to the coil.

Figure 5: Reconstructed human brain conductivity maps using the phase of the 32-channel receiver virtual coil (left) and the 32-channel receiver GCC-Align coil (center). Virtual coil magnitude is shown at right for reference.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1564