Parallel imaging performance investigation of an 8-channel common-mode differential-mode (CMDM) planar array for 7T MRI
Introduction
Parallel imaging techniques (1-4) using array coils and multichannel receivers (5) in the past decade provide an effective approach to reduce scan time of magnetic resonance imaging (MRI). In parallel imaging, the acquisition time can be dramatically reduced by exploiting the variations in the sensitivity patterns of different elements in an array coil in order to replace some of the spatial encoding normally accomplished by using magnetic field gradients. Numerous data acquisition and image reconstruction strategies have been proposed and adopted for clinical applications in high field MRI such as real time cardiovascular imaging, functional MRI, and contrast-enhanced MR studies (6-9).
In order to perform the parallel imaging techniques, independent multichannel receivers with different sensitivity distribution are demanded to acquire signals simultaneously. The conventional MRI phased arrays with loop structured resonant elements were developed and investigated in previous studies (10-14). Some other types of multichannel receivers were also proposed, such as saddle structure (15) and volume strip arrays (16). It has been demonstrated that high and ultrahigh magnetic fields are fundamentally advantageous in in vivo MR applications due to their inherently high signal-to-noise ratio (SNR), and thus high spatial and spectral resolution (17-20). At high and ultrahigh fields, the interactions between the EM fields and biological samples become distinct. This results in tissue-specific perturbations of EM field patterns, requiring appropriate RF coil designs to improve image quality and to avoid adverse effects in patients. RF coil and coil array designs using microstrip transmission lines were proposed to address the high field RF challenges and electromagnetic coupling issues. Based on the microstrip transmission line theory, the microstrip surface coils were developed at the ultrahigh field of 7T (21,22), which provides the capability of developing a variety of coils such as volume coils (23,24) and phased array coils (25-28). Subsequently, common-mode differential-mode (CMDM) coils implemented by using microstrip transmission line for double-tuned volume coils at 7T were developed (29), showing the capability of intrinsical decoupling between common mode and differential mode resulting from their orthogonal field distribution. Recently, the non-array planar transceiver coil using CMDM method proposed and tested for 7T MR imaging (30), demonstrating high quality factors, high frequency operation capability and high efficiency. The study also suggested that the planar CMDM coil could be a good candidate for design blocks in multichannel RF coil arrays.
In this work, an 8-channel phased array was based on the planar CMDM coil for ultrahigh field MRI. Numerical simulation was performed to investigate electromagnetic field distribution of the CMDM planar phased array. The image of each channel was calculated based on the simulation results. The combined image of the 8-channel CMDM array was obtained by using root-sum-of-squares (rSoS) and covariance weighted root-sum-of-squares (Cov-rSoS) (31). The Cartesian SENSE was employed to investigate imaging acceleration capability. SNR maps, g-factor maps and artifact power at various reduction factors were calculated to evaluate the parallel imaging performance of the proposed array.
Materials and methods
The proposed 8-channel CMDM array consisted of four CMDM blocks as shown in Figure 1. For each CMDM block, the width and length of the block were W2 =43 mm and L2 =94 mm respectively. The width and length of the DM driving loop were W1 =26 mm and L1 =39 mm. The thickness H, width W and length L of the substrate were 12.7, 120 and 223.05 mm respectively. The gap G between adjacent CMDM blocks was 6.35 mm. The width of the copper strips of the driving loops and CMDM blocks were 1.6 and 6.35 mm, respectively. The cubic phantom with 75×75×268 mm3 at 10 mm distance from the top of the substrate was filled with water (εr =78, μr =0.99, σ =1.59 S/m).
Numerical method for full-wave 3D electromagnetic analysis was used to simulate the electromagnetic field distribution of the 8-channel CMDM array. The infinite space was truncated by the perfect matched layers (PML). The distance between the PML and the model was set to 100 mm in all directions. Four layers of PML were exploited to mimic infinite space. Total mesh number of the model was 858,000 with adaptive mesh discretization. The 8-channel CMDM array was tuned to 298 MHz by capacitor Ccm-t1 =91 pF for four CMDM blocks, which corresponds to the proton Larmor frequency at 7T. The value of the capacitor Ccm-t2 was 2.2 pF for CMDM_1 and CMDM_4, and 2.5 pF for CMDM_2 and CMDM_3. The value of tune capacitor Cdm-t for four CMDM blocks was in turn from left to right 1.1, 1.2, 1.2 and 1.0 pF as shown in Figure 1. Eight ports of the array were excited separately by the same RF pulse with equal amplitude and phase. The excitation was set to be current source with input power of 1 W. The electromagnetic field distribution of each channel was acquired for array performance investigation.
Neglecting the effects of the relaxation and susceptibility, the receiving signal intensity S of the gradient echo (GRE) sequence can be written as (32-34):
[1] |
where W is the water content in the voxel, γ and ι are the magnetogyric ratio and the duration, respectively. and are the left-handed rotating field and right-handed rotating field, which correspond to transmit field and receive field of RF coil, respectively. The V is a variable and proportional to coil driving voltage in a given experiment, and thus proportional to flip angle. Assuming that the excitation field was homogenous, the
Two reconstruction methods, rSoS and Cov-rSoS, were employed to generate combined image of the 8-channel CMDM array. Let C and S denote the coil sensitivity vector and aliased image vector respectively. The general expressions of rSoS and Cov-rSoS were described as (31):
[2] |
[3] |
where Ψ denoted the noise covariance matrix, which was calculated by the electric field in phantom as (5):
[4] |
where σ(r) was the conductivity, Ei(r) was the electric field excited by the i-th excitation port, the superscript H indicated conjugate (Hermitian) transpose. If C equaled to S, then
[5] |
[6] |
For the above combination methods, the image SNR was given by:
[7] |
[8] |
In order to evaluate parallel imaging performance of the 8-channel CMDM array, SENSE method was utilized for imaging acceleration. The image intensity of each channel was transformed into k-space as full k-space data. Then, full k-space data were under-sampled by different reduction factors. A non-aliased full-FOV image can be obtained through eight aliased images by implementing the SENSE method [PULSAR toolbox (35)].
When the reduction factor R was 1, the image reconstructed in SENSE was degenerated into Cov-rSoS image, which was an accurate reconstruction image and considered as a reference image in reconstruction error estimation. The SNR for SENSE reconstruction was calculated by:
[9] |
where R and g were reduction factor and g-factor map, respectively. And SNRfull denoted the SNR of full k-space image. SNRfull was replaced by SNRrSoS and SNRCov-rSoS in this study.
The artifact power was used to evaluate the reconstruction error of different reduction factors R. The image of Cov-rSoS reconstruction was considered as reference image. The construction error was evaluated by L2-norm in simulations. The artifact power (AP) was defined as:
[10] |
Results
The simulated distribution excited by each channel separately in transverse plane of the 8-channel CMDM array was shown in Figure 2. The SENSE reconstruction images for 8-channel CMDM array at various reduction factor R were displayed in Figure 3. When reduction factor R =1, the SENSE reconstruction degenerated into Cov-rSoS. For comparisons, the SENSE reconstruction images and corresponding error for 4-channel CM and DM were also given in Figures 4,5, respectively. The SENSE reconstruction images with reduction factor R =2 were almost the same as Cov-rSoS image. The images of 8-channel CMDM array provided better penetration and coverage than the 4-channel CM array and DM array.
The SNR maps for SENSE reconstruction at various reduction factors were demonstrated in Figures 6-8. Figure 6 showed that SNR of 8-channel CMDM array for rSoS and Cov-rSoS reconstruction with R =1, 2, 4. Similarly, the Figures 7,8 displayed that SNR maps of 4-channel CM and DM with R =1, 2, 4. The results demonstrated that 8-channel CMDM array for Cov-rSoS reconstruction with R =1 provided the highest SNR than others. The SNR performance of the 8-channel CMDM array with different reduction factors was better than the 4-channel CM and DM array.
The g-factor for SENSE reconstruction was calculated. The means of g-factor in ROI corresponding to reduction factor R =2, 4 were 1.4355 and 3.3757, respectively. While the maximum of g-factor in ROI corresponding to reduction factor R =2, 4 were 2.408 and 11.6621, respectively. The 1/g-factor maps for SENSE reconstruction at various reduction factors for 8-channel CMDM array were shown in Figure 9. The ROI was selected as shown in Figure 1B.
Finally, the artifact power of SENSE reconstruction for the 8-channel was considered according to Eq. [10]. In Table 1, the artifact power was given for 8-channel CMDM, 4-channel CM, 4-channel DM with various reduction factors. The artifact power of 4-channel CM array was approximately equal to the 4-channel DM array. The artifact power of 8-channel CMDM array was reduced compared with 4-channel CM and DM array.
Full table
Discussion and conclusions
In this study, rSoS and Cov-rSoS are employed to combine the images of 8-channel CMDM array. According to the SNR of combined images in the simulation, Cov-rSoS is superior to rSoS, taking into account the correlated noise among channels.
The reconstructed images acquired using 8-channel CMDM array is compared with that acquired from 4-channel CM and DM array. The SNR of 8-channel CMDM array at various reduction factors is better than the 4-channel CM and DM array. In addition, parallel imaging performance of 8-channel CMDM array based on simulation has been evaluated. The artifact power of 8-channel CMDM array is less than the 4-channel CM and DM array. These indicate the parallel imaging feasibility and better performance of the 8-channel CMDM array.
Practically, the excitation field is not homogeneous at 7T. RF shimming can be employed to eliminate or decrease inhomogeneity of excitation field. Therefore the signal intensity is strongly related to the RF shimming strategies. In this work, in order to focus on the parallel imaging performance evaluation of the proposed array, we only used the in the simulation, which is not indicated the transmit field is homogenous. In our simulation, the transmit field of each channel of the proposed 8-channel CMDM array is shown in Figure 10. The combined transmit field is shown in Figure 11.
In the simulation experiments, there is no noise adding into each channel. The error is zero when reduction factor R is equal to 1. Since the g-factor values demonstrate the noise increase caused by the SENSE reconstruction, it is possible to achieve lower artifact power while the g-factors are high. In order to investigate the artifact power in images with noise, we add Gaussian white noise into image of each channel. The Gaussian white noise of each channel is random and independent from other channels. The signal intensity of the images of each channel is approximate to 5×10–7. We add Gaussian white noise with mean values equal to 0 and standard deviation is 5×10–8 and 1×10–8 respectively. The SNR maps of Cov-rSoS combination and AP at various reduction factors are shown in Figures 12,13 and Tables 2,3. The AP and SNR variation at different noise level could be demonstrated.
Full table
Full table
When the covariance weighted root-sum-of-squares (Cov-rSoS) were employed to combine images of the 8-channel CMDM array in simulation, coupling of each channel was partially compensated by the Cov-rSoS reconstruction method. This is agreed with the result in references (36,37). In practice, the preamplifier decoupling could also be used to further reduce the coupling between the elements. For array fabrication, several decoupling methods such as adjusting the gaps between the adjacent CMDM block, magnetic wall (or metamaterial) decoupling (38) and resonant inductive decoupling (39) might be employed to further reduce the coupling of each channel. Future work includes evaluating the parallel imaging performance of these arrays and decoupling methods with in vivo experiments.
Acknowledgements
This work was supported in part by the National Natural Science Foundation of China under Grant No. 81120108012, 51307171 and 51228702, the Major State Basic Research Development Program of China (973 Program) under Grant No. 2011CB707903, Shenzhen special funds for the development of strategic emerging industries under Contract No. SY294261001 and innovation fund for outstanding youth in Shenzhen Institutes of Advanced Technology of Chinese Academy of Sciences under Grant No. Y3G016.
Disclosure: The authors declare no conflict of interest.
References
- Sodickson DK, Manning WJ. Simultaneous acquisition of spatial harmonics (SMASH): fast imaging with radiofrequency coil arrays. Magn Reson Med 1997;38:591-603. [PubMed]
- Pruessmann KP, Weiger M, Scheidegger MB, et al. SENSE: sensitivity encoding for fast MRI. Magn Reson Med 1999;42:952-62. [PubMed]
- Griswold MA, Jakob PM, Nittka M, et al. Partially parallel imaging with localized sensitivities (PILS). Magn Reson Med 2000;44:602-9. [PubMed]
- Griswold MA, Jakob PM, Heidemann RM, et al. Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med 2002;47:1202-10. [PubMed]
- Roemer PB, Edelstein WA, Hayes CE, et al. The NMR phased array. Magn Reson Med 1990;16:192-225. [PubMed]
- Weiger M, Pruessmann KP, Boesiger P. Cardiac real-time imaging using SENSE. SENSitivity Encoding scheme. Magn Reson Med 2000;43:177-84. [PubMed]
- Huber ME, Kozerke S, Pruessmann KP, et al. Sensitivity-encoded coronary MRA at 3T. Magn Reson Med 2004;52:221-7. [PubMed]
- van den Brink JS, Watanabe Y, Kuhl CK, et al. Implications of SENSE MR in routine clinical practice. Eur J Radiol 2003;46:3-27. [PubMed]
- Park J, McCarthy R, Li D. Feasibility and performance of breath-hold 3D true-FISP coronary MRA using self-calibrating parallel acquisition. Magn Reson Med 2004;52:7-13. [PubMed]
- Fayad ZA, Connick TJ, Axel L. An improved quadrature or phased-array coil for MR cardiac imaging. Magn Reson Med 1995;34:186-93. [PubMed]
- Porter JR, Wright SM, Reykowski A. A 16-element phased-array head coil. Magn Reson Med 1998;40:272-9. [PubMed]
- Zhu Y, Hardy CJ, Sodickson DK, et al. Highly parallel volumetric imaging with a 32-element RF coil array. Magn Reson Med 2004;52:869-77. [PubMed]
- Balu N, Yarnykh VL, Scholnick J, et al. Improvements in carotid plaque imaging using a new eight-element phased array coil at 3T. J Magn Reson Imaging 2009;30:1209-14. [PubMed]
- Li Y, Pang Y, Vigneron D, et al. Investigation of multichannel phased array performance for fetal MR imaging on 1.5T clinical MR system. Quant Imaging Med Surg 2011;1:24-30. [PubMed]
- Zhang X, Webb A. Design of a capacitively decoupled transmit/receive NMR phased array for high field microscopy at 14.1T. J Magn Reson 2004;170:149-55. [PubMed]
- Lee RF, Xue R. A transmit/receive volume strip array and its mode mixing theory in MRI. Magn Reson Imaging 2007;25:1312-32. [PubMed]
- Hoult DI, Richards RE. The signal-to-noise ratio of the nuclear magnetic resonance experiment. J Magn Reson 1976;24:71-85. [PubMed]
- Hoult DI, Phil D. Sensitivity and power deposition in a high-field imaging experiment. J Magn Reson Imaging 2000;12:46-67. [PubMed]
- Vaughan JT, Garwood M, Collins CM, et al. 7T vs. 4T: RF power, homogeneity, and signal-to-noise comparison in head images. Magn Reson Med 2001;46:24-30. [PubMed]
- Adriany G, Auerbach EJ, Snyder CJ, et al. A 32-channel lattice transmission line array for parallel transmit and receive MRI at 7 tesla. Magn Reson Med 2010;63:1478-85. [PubMed]
- Zhang X, Ugurbil K, Chen W. Microstrip RF surface coil design for extremely high-field MRI and spectroscopy. Magn Reson Med 2001;46:443-50. [PubMed]
- Pang Y, Wu B, Wang C, et al. Numerical Analysis of Human Sample Effect on RF Penetration and Liver MR Imaging at Ultrahigh Field. Concepts Magn Reson Part B Magn Reson Eng 2011;39B:206-16. [PubMed]
- Zhang X, Zhu XH, Chen W. Higher-order harmonic transmission-line RF coil design for MR applications. Magn Reson Med 2005;53:1234-9. [PubMed]
- Zhang X, Ugurbil K, Sainati R, et al. An inverted-microstrip resonator for human head proton MR imaging at 7 tesla. IEEE Trans Biomed Eng 2005;52:495-504. [PubMed]
- Adriany G, Van de Moortele PF, Wiesinger F, et al. Transmit and receive transmission line arrays for 7 Tesla parallel imaging. Magn Reson Med 2005;53:434-45. [PubMed]
- Wu B, Wang C, Kelley DA, et al. Shielded microstrip array for 7T human MR imaging. IEEE Trans Med Imaging 2010;29:179-84. [PubMed]
- Wu B, Wang C, Lu J, et al. Multi-channel microstrip transceiver arrays using harmonics for high field MR imaging in humans. IEEE Trans Med Imaging 2012;31:183-91. [PubMed]
- Wu B, Zhang X, Wang C, et al. Flexible transceiver array for ultrahigh field human MR imaging. Magn Reson Med 2012;68:1332-8. [PubMed]
- Pang Y, Zhang X, Xie Z, et al. Common-mode differential-mode (CMDM) method for double-nuclear MR signal excitation and reception at ultrahigh fields. IEEE Trans Med Imaging 2011;30:1965-73. [PubMed]
- Li Y, Yu B, Pang Y, et al. Planar quadrature RF transceiver design using common-mode differential-mode (CMDM) transmission line method for 7T MR imaging. PLoS One 2013;8:e80428. [PubMed]
- Keil B, Wald LL. Massively parallel MRI detector arrays. J Magn Reson 2013;229:75-89. [PubMed]
- Hoult DI, Phil D. Sensitivity and power deposition in a high-field imaging experiment. J Magn Reson Imaging 2000;12:46-67. [PubMed]
- Hoult DI. The principle of reciprocity in signal strength calculations—A mathematical guide. Concepts in Magnetic Resonance 2000;14:173-87.
- Collins CM, Yang QX, Wang JH, et al. Different excitation and reception distributions with a single-loop transmit-receive surface coil near a head-sized spherical phantom at 300 MHz. Magn Reson Med 2002;47:1026-8. [PubMed]
- Ji JX, Son JB, Rane SD. PULSAR: A Matlab toolbox for parallel magnetic resonance imaging using array coils and multiple channel receivers. Concepts in Magnetic Resonance Part B: Magnetic Resonance Engineering 2000;37B:24-36.
- Nordmeyer-Massner JA, De Zanche N, Pruessmann KP. Mechanically adjustable coil array for wrist MRI. Magn Reson Med 2009;61:429-38. [PubMed]
- Ohliger MA, Ledden P, McKenzie CA, et al. Effects of inductive coupling on parallel MR image reconstructions. Magn Reson Med 2004;52:628-39. [PubMed]
- Li Y, Xie Z, Pang Y, et al. ICE decoupling technique for RF coil array designs. Med Phys 2011;38:4086-93. [PubMed]
- Avdievich NI, Pan JW, Hetherington HP. Resonant inductive decoupling (RID) for transceiver arrays to compensate for both reactive and resistive components of the mutual impedance. NMR Biomed 2013;26:1547-54. [PubMed]