Impact of number of channels on RF shimming at 3T

Object At high-field strengths (≥3T) inhomogeneity of the radio frequency (RF) field and RF power deposition become increasingly problematic. Parallel Transmission (PTx)—the use of segmented transmission arrays with independently driven elements—affords the ability to combat both of these issues. There are a variety of existing designs for PTx coils, ranging from systems with two channels to systems with eight or more. In this work, we have investigated the impact of the number of independent channels on the achievable results for both homogeneity improvement and power reduction in vivo. Materials and methods A 3T Philips Achieva MRI system fitted with an 8-channel PTx body coil was driven so as to emulate configurations with 1, 2 4 and 8 independent channels. RF shimming was used in two different anatomies in order to assess improvements in RF homogeneity. Results Significant homogeneity improvements were observed when increasing from 1 to 2, 2 to 4, and 4 to 8 channel configurations. Reductions in RF power requirements and local SAR were predicted for increasing numbers of channels. Conclusion Increasing the number of RF transmit channels adds extra degrees of freedom which can be used to benefit homogeneity improvement or power reduction for body imaging at 3T.


Introduction
At high field strengths (3T and greater), the excitation radio-frequency (RF) wavelength becomes comparable to the size of the human body, causing strong coil-tissue interactions which are subject-specific, depending on geometry and tissue electrical properties. The result of these interactions is an inhomogeneous RF transmit (B 1 ? ) field, which leads to spatially varying image contrast. Furthermore, conventional single-channel RF systems generally have a one-time adjustment to match the RF coil to a nominal load presented by an ''average patient'', whereas in fact the loading of the RF coil is subjectdependent, which can result in variable performance. The use of multiple transmit channels offers a means for mitigating these effects by allowing for greater control over the RF electromagnetic fields. RF shimming [1], in which a subject-specific adjustment is made to the relative amplitudes and phases of each transmit channel, can be used to improve B 1 ? homogeneity and/or reduce SAR. While such techniques have been technically demonstrated for transmit arrays with eight or more channels [2,3], systematic clinical studies of RF shimming have primarily used commercially available two channel systems [4]. RF shimming at 3T has been used to stabilise and improve B 1 ? homogeneity while achieving shorter image acquisition times compared with single-channel quadrature excitation [4], and to also improve diagnostic quality of resulting images [5]. These improvements come directly from the ability to manipulate the spatial B 1 ? distribution as a linear sum of contributing fields from individual coil elements: In the above expression, r i (r) is the transmit sensitivity of the ith coil, N c is the total number of coils, S(r) is the overall transmit sensitivity, and b i are complex weights applied to each channel. Transmit sensitivity is defined as the ratio of the actual B 1 ? field to the nominal value; it is a dimensionless value which describes spatial variation of the field amplitude and phase. Equation 1 shows that the ability to obtain uniform overall sensitivity depends on the sensitivities of the different coils in the array. It might be expected, therefore, that more coils, or perhaps different coil designs, would yield even larger improvements than a two-channel configuration. In this investigation, a 3T MRI scanner with eight physical transmit channels has been used to emulate different 2-, 4-and 8-channel configurations in vivo. Possible benefits of RF shimming were assessed in a multiple-volunteer study focusing on B 1 ? homogeneity, image quality and RF power reduction.

Methods
Experiments were performed using an 3T Achieva MRI system (Philips Healthcare, Best, The Netherlands) equipped with an 8-channel parallel transmit (PTx) body coil. The coil [3] consists of an array of eight elements, arranged around the bore of the magnet to give maximum spatial variation in transmit sensitivity in the axial plane.

RF shim framework
The overall transmit sensitivity of the array coil is a linear sum depending on the complex weighting applied to each channel (Eq. 1). This can be rewritten in matrix form as S = rb, where S is an N 9 1 vector (N is the number of voxels), r is an N 9 N c sensitivity matrix and b is an N c 9 1 vector containing the complex coil weightings for each channel. For the system used, N c = 8 and b = [1 1 1 1 1 1 1 1] T is defined as quadrature mode (single-channel) transmission. This definition of ''quadrature'' mode depends on the system being calibrated; this will be discussed later. Emulation of systems with a reduced number of channels can be performed by introducing an N c 9 N r channel reduction matrix, R. This reduction matrix allows for emulation of systems with N r \ N c channels via linear combinations of the scanner's physical channels. The transmit sensitivity is then written S = rRb, and b represents the weightings applied now not to individual physical coils, but to the N r emulated channels. Full control of the physical coils implies R is the N c 9 N c identity matrix. The phase/amplitude relationship between all channels that defines quadrature mode is a property of the system hardware that must be determined by calibration. For our system, this is done using a standard phantom. Choice of the calibration conditions does not affect full 8-channel operation. However, it adds a bias to reduced channel combinations by locking in some aspects of the calibration. In order to reduce the effect that this has, a new system calibration matrix was computed by performing full 8-channel shimming on 10 previously acquired pelvic B 1 ? maps (10 different subjects) individually, then taking the complex average of all results for each channel. The resulting weights c i ,which provide an optimised average starting shim, can be incorporated into a diagonal matrix, C ii = c i which is included in the overall expression for the combined system sensitivity: S = rCRb. The calibration matrix was explicitly defined prior to any experimental work and not altered subsequently.
RF Shim weights were calculated using Magnitude Least Squares (MLS) optimisation [6,7] using the local variable exchange method S target is the target combined transmit sensitivity; in this work S target = 1.0 for all space. The regularisation parameter k is designed to control total RF power. Note that the solution, b, has N r entries; the actual weights applied to the N c physical coils are determined from the product Rb.
Hence the second term in Eq. 2, which enumerates the total RF power, considers this product. The overall scaling of R is not important in the term k||Rb|| 2 since scalar factors in R can be subsumed into the parameter k. Several candidate reduction matrices for emulating 2-and 4-channel systems were tested on the pilot data. Candidates included phased combinations, Fourier modes [2] and pair-wise coil combinations. These were evaluated using Eq. 2 with a range of values of k to find the combinations that produced the optimal combination of B 1 ? homogeneity and lowest total RF power when applied to the pilot data. The optimised 4-channel reduction matrix used in this study is shown below; in this case adjacent coils were paired together, corresponding to creation of an enveloping array with fewer, larger coil elements placed adjacently.
Two 2-channel combinations were used. The first used the following reduction matrix which was found to yield the best results in the pilot data. The second was chosen to emulate orthogonal, nominally linear modes of a 2-channel birdcage system and was defined as: , calculated all shim settings for a subject in less than 10 s, allowing for straightforward implementation in a scanning session. The curves were used to select two solutions for each case: the first maximised the homogeneity improvement obtainable whilst using the same RF power level as quadrature mode excitation, the second minimised the RF power level while maintaining approximately the same homogeneity as quadrature excitation. Figure 1 shows how these solutions were obtained using L-curves calculated for one of the study subjects. Imaging experiments were performed for the first solutions (maximised homogeneity) where measurable differences in B 1 ? were expected; in order to keep the examination duration within 1 h for each subject, the minimized power solutions were not used for imaging.

Imaging experiments
RF shimming experiments were carried out on healthy volunteers, with the focus being thigh and pelvis imaging. A six-element phased array coil was used for signal reception. Research ethics committee approval was obtained for the study, and all participants gave written informed consent prior to enrolment. In total, seven healthy volunteers (two male, five female) underwent pelvic and/or thigh scanning (five volunteers underwent both).
Axial B 1 ? field maps were acquired using the actual flip angle imaging sequence [8] including a modified spoiling regime [9], slice profile correction [10] and linear combination mapping to improve SNR [11][12][13]. Imaging parameters were FOV = 440 mm 9 220 mm, 64 9 32 matrix, slice thickness 10 mm, nominal flip angle = 808 and pulse repetition times of 30 ms and 150 ms. A T1-weighted turbo spin echo (T1w-TSE) sequence (440 mm 9 220 mm FOV, 400 9 318 matrix, TR = 605 ms, TE = 10 ms, TSE turbo factor = 3) and a T2-weighted turbo spin echo (T2w-TSE) sequence (440 mm 9 220 mm FOV, 400 9 318 matrix, TR = 2,000 ms, TE = 100 ms, TSE turbo factor = 10) were chosen for body imaging so the effect of varying the number of channels had on image quality could be assessed. All imaging was performed in the transverse plane using a single slice acquisition. The TSE sequences outlined above are often used to acquire multiple interleaved slices; use of a single slice allows the SAR to be kept as low as possible while maintaining the same contrast parameters as used in clinical examinations. A full examination, consisting of B 1 ? mapping, RF shim calculation and acquisition of images for all channel combinations, could be completed within 45 min. B 1 ? maps were also obtained for each optimized shim setting to allow comparison between the predicted and achieved shimmed B 1 ? distributions, independent of the receiver coil.

RF safety and SAR analysis
Accurate calculation of SAR, particularly the local SAR, requires subject specific E-field models based on the individual's anatomy. Obtaining this information and completing the required calculations remains a challenge to achieve during an examination and is a current research topic in its own right [14,15]. In this study whole body SAR was dealt with simplistically using the ||Rb|| 2 term in the RF shimming cost function (Eq. 2). This method allowed for a time-efficient, indirect minimisation of global SAR through the limiting of total RF power. A singlesubject, single-position E-field model was also used to calculate the local SAR prior to every scan performed using the method from Ref. [16], allowing for quantification of the SAR reduction potential of multiple transmit channels. The E-field model is for a large male subject, positioned with torso at isocentre, therefore the results acquired are unlikely to be accurate for different body shapes and different body positions relative to the RF coils. Due to the uncertainty in the SAR taken from this model, a maximum estimated SAR level of 10 % of the relevant regulatory limit [17] was allowed for any imaging sequence, to allow a tenfold margin for error. . RMSE quantifies the deviation in absolute terms from the desired field pattern, and is equivalent to the square root of the first term in Eq. 2. The coefficient of variation (standard deviation divided by mean) is purely a measure of spatial homogeneity; unlike RMSE, C v is not altered by scaling the B 1 ? fields. We hypothesised that increasing the number of channels would improve homogeneity, and therefore decrease the C V . Statistical analysis was carried out on the C V data using SPSS (IBM). Kolmogorov-Smirnov tests were used to test the normality of the data (there was no evidence for non-normality), then paired sample t-tests were used to examine differences in C V for incremental increases in the number of channels controlled (2 vs. 1 (quad), 4 vs. 2, 8 vs. 4). The requirement for statistical significance was p \ 0.05 for one-tailed comparisons (i.e. the alternative hypothesis is that C V for the configuration with more channels is less). A blind grading system was used to assess changes in imaging performance. An experienced radiologist was given all images in a randomised order and asked to grade them according to a five-point scoring system, previously described by Willinek et al. [4], but modified as all volunteers were healthy and no pathology was assessed. Perceived image quality was rated on a scale of 1-5, with 1 being an image of non-diagnostic quality and 5 demonstrating optimal contrast with no B 1 ? inhomogeneity artefacts.

B 1
? homogeneity and Image quality improvements A consistent pattern of substantial inhomogeneity was found in the B 1 ? field generated by quadrature excitation, in both the pelvis and thighs of all volunteers. The inhomogeneity was characterised by diagonally offset areas of low B 1 ? intensity in the anterior and posterior regions. This resulted in areas of very low signal intensity in the resulting T 1 w-TSE images; the T 2 w-TSE images were more robust to the inhomogeneity. The B 1 ? fields and the corresponding T 1 w and T 2 w images acquired in quadrature mode for all volunteers are shown in Fig. 2 for the pelvis and Fig. 3 for the thighs. Note that the variability seen in the average B 1 ? value between different subjects in Fig. 2 is the result of inaccuracy in the scanner's standard global poweroptimisation procedure [18]. Figure 4 contains RF homogeneity, power and SAR estimates for the pelvis and thigh experiments. All shim solutions had an approximately constant RF forward drive, relative to quadrature mode, as can be seen in parts e, f. The maximum local SAR as indicated by the E-field model based SAR calculator is displayed on parts g, h; note that the local SAR measurements are quoted for the AFI B 1 ? mapping sequence. This figure shows that by increasing the number of transmit channels, the B 1 ? field can be made progressively more homogeneous for both the pelvis and thighs. The results of paired sample t-tests comparing C V for configurations with successively increasing numbers of channels showed that improvements were significant (Table 1) for both anatomies studied and for all increases in channel number. Examples of the effect of these improvements in single volunteers can be seen in Figs. 5 and 6, which contain the B 1 ? maps, T1w-TSE and T2w-TSE images obtained for each configuration from a single volunteer.
The mean maximum local SAR measurements from the online SAR model decreased with increasing numbers of channels for the pelvis experiment, although 4-and 8-channel configurations performed similarly. In the thighs there was also a gradual reduction with the exception of the 2-channel (non-birdcage) configuration which achieved a much lower calculated average local SAR than the other configurations. When transmitting in single-channel (nominal quadrature) mode, the maximum local SAR for the AFI sequence is always 3 W/kg according to the experimental online electric field model as no phase/ amplitude modifications have been made and the system is calibrated.
The blind grading results from the anatomical images for all coil configurations averaged over all subjects are displayed in Table 2. In the case of the T1w images, increasing channel numbers also increases perceived image performance for both anatomies tested. For pelvic imaging there were no perceived improvements in the T2w images when increasing the number of transmit channels; however, a slight improvement was noted when imaging the thighs.

RF power reduction with increasing channel numbers
For all experiments, optimized RF shims were also computed for minimizing RF power (quantified by ||Rb|| 2 ) while maintaining a fixed homogeneity level (constant RMSE). No imaging was performed for these solutions; B 1 ? homogeneity was instead estimated from the predicted field maps for the solutions found. The results are summarised by Fig. 7. In both anatomies, the total RF power reduces as the number of channels increases, with the 8-channel system having the lowest power requirement. The predicted local SAR is also reduced for channel counts higher than one, however these do not follow the same pattern as total RF power. In addition the C V decreased with increasing channel number in both cases.

Discussion
In this study, conducted at 3T, pronounced areas of B 1 ? inhomogeneity were found in the body with the areas of low B 1 ? intensity following a distinctive diagonal characteristic of 'elliptic' shaped loads, which the human body loosely resembles [19]. Our results show that increasing the number of RF channels, and hence the number of effective degrees of freedom, leads to successive improvement in the achievable B 1 ? homogeneity. This is visible from the diagram in Fig. 1, which was generated using one of the datasets included in the study. The solutions found for configurations with more degrees of freedom are nested within those found with fewer degrees of freedom, indicating that a more favourable trade-off between power and inhomogeneity can be achieved by increasing the channel Fig. 2 Nominal quadrature mode axial pelvis results for multiple volunteers. Axial AFI transmit sensitivity maps (S(r)) (top), T 1 w-TSE images (middle) and T 2 w-TSE images (bottom). Homogeneity is variable in all volunteers. Shading artefacts are observed particularly in the T 1 w images. Note that the T2w image is missing for one volunteer because technical problems led to the examination being terminated before all data had been acquired number. A systematic study using multiple volunteers confirmed that increasing the number of channels from one to two, two to four, and four to eight all resulted in statistically significant improvements in homogeneity. Improved RF homogeneity also led to a quantifiable improvement in image quality as measured by a blinded Taking 0.05 as the threshold for significance, increasing from one to two, two to four channels and four to eight channels yielded significant improvements in homogeneity in all cases. Note that only the optimal 2-channel configuration (not the 'birdcage' configuration) is shown here to reduce the number of comparisons evaluation, especially in T1 weighted images. The T2 weighted TSE sequence used was found to be more robust to B 1 ? homogeneity. The amplitude of a single spin echo produced by refocusing flip angle a is proportional to sin 2 (a/2), whereas the pseudo steady-state echo amplitude produced by a TSE sequence with multiple refocusing is (birdcage), 2-channel (optimal), 4-channel and 8-channel. The T 2 w images appear more robust to B 1 ? inhomogeneity than the T1w sequence Fig. 6 Example data from thigh scans of one volunteer. Top AFI axial thigh S(r) maps, middle T 1 w-TSE and bottom T 2 w-TSE axial pelvic images; columns from left to right: quadrature, 2-channel (birdcage), 2-channel (optimal), 4-channel and 8-channel. The T 2 w images appear more robust to B 1 ? inhomogeneity than the T1w sequence proportional to sin(a/2) [20]. For both the T1 and T2 weighted sequences, a was 130°, however the former used centric phase encoding, while the latter used linear. As a result it would be expected that the T1 weighted images would be more influenced by the first echo (proportional to sin 2 (a/2)) while the T2 images would be closer to the pseudo steady-state (sin(a/2)), which is in accordance with the finding that the T2 weighted sequence is more robust to variations in B 1 ? . As well as B 1 ? field homogeneity, parallel transmission also allows greater control over both global and local SAR. RF shims were calculated for maximal homogeneity at a fixed RF drive power level, where drive power is the sum of squares of the RF drive levels supplied to each channel, quantified by ||Rb|| 2 . A single station E-field model for a single subject was available for the estimation of peak local SAR for these solutions (Fig. 4g, h for pelvis and thigh respectively). These results indicate that there was some reduction in local SAR for the shimmed solutions, but importantly none of the shimmed solutions produced a higher local SAR than quadrature mode, according to the model. It should be noted that this optimisation did not aim to minimise local SAR, and in fact held relative power constant. A second set of experiments which held homogeneity constant and aimed to reduce total power did indeed show consecutive reductions in total power as channel number increased. The local SAR, as predicted by the model, also decreased when compared with quadrature Fig. 7 Comparisons of homogeneity, power requirements and local SAR with different numbers of channels when performing RF shimming at a fixed level of homogeneity. Unlike Fig. 4, homogeneity values (RMSE/C V ) are predicted rather than measured. Pelvis results are shown in the left column, thigh results on the right. a, b RMSE is held the same for all configurations; variability comes from inter-subject differences. c, d C v falls with increasing channel count, indicating that although the RMSE is held constant the actual field patterns do vary. e, f RF power requirements fall as the number of channels increases. g, h Local SAR prediction of AFI B 1 ? mapping sequence acquired with power reduction shims mode, but again decreases were erratic. This is partly to be expected because the SAR model did not accurately describe the imaging scenarios, either in the positioning or the size and build of the subject. Personalised models for each subject might allow a subject-specific estimation of the local SAR within the current framework. These could be used to much greater effect by allowing the replacement of the regularization term in Eq. 2 with a term based directly on a true SAR estimate. A recent study of this type by Homann et al. [15] investigated shimming with an 8-channel coil of a similar design to the one used in this study. In that study, it was found that including a subjectspecific SAR model when determining shim settings produced substantial SAR reductions, which is in line with our findings using a single standard model and in relation to total input power. The method employed in the current study, in which an L-curve is calculated to explore balancing B 1 ? homogeneity with total drive power, or ideally with peak local SAR, allows flexible choice of operating point. In this study, constant power or constant homogeneity solutions were examined, but in practice any point along the curves could be selected, depending on the situation. As indicated by the imaging examples presented in Figs. 2, 3, 5 and 6, different sequences have variable sensitivity to B 1 ? inhomogeneity. The T2w sequence employed here is more robust, and so a lower power solution might be selected to increase efficiency, whereas other parts of an examination may be much more vulnerable and so warrant a different homogeneity-power tradeoff. This allows flexible balancing of SAR load across a whole examination. At the other extreme, Fig. 1 shows that it is possible to achieve highly homogeneous B 1 ? fields over a large (axial) field of view with a multi-transmit system at 3T if the total power is not a concern. However, useful solutions emerge as soon as any power considerations are introduced, in the case of the method used here, by setting k to any non-zero value.
The lack of suitability of the SAR model for the subjects imaged meant that the true SAR values were uncertain. For this reason, a conservative approach was taken, and imaging was limited to run at 10 % of the relevant SAR limit. This limitation is the subject of a considerable research effort in the wider community: the use of generic models [21], simplified models generated in real time [22] or even actual in situ measurement of the SAR distribution [23] are all current possibilities. More accurate models would allow subject-specific SAR reductions to be calculated; if predictions could be treated with confidence, then reductions in quoted SAR would be translatable into an increase in permissible RF duty cycle, which can be used to reduce TR or increase the maximum number of slices within a fixed TR. If the number of slices were increased, extra B 1 ? maps would also need to be acquired, potentially slowing the acquisition process. New methods such as that in Ref. [24] provide the possibility of rapid B 1 ? measurement over large volumes, alleviating this issue. Similarly, while computation times were acceptable for this study, they would not scale well to multiple slices. The code used was, however, not optimized for speed since it was not an issue here, though it could be accelerated if necessary.
The human body, with its approximately elliptic cylindrical geometry, produces field variation mainly in the axial plane [19]. As a result, prevailing coil designs consist of elements distributed azimuthally around the magnet, giving maximal axial control. Consequently, in this study all imaging and field measurements were obtained in axial planes, where the effects of field inhomogeneity are generally greatest and the existing coil design offers its maximum flexibility in field control. We have found qualitatively similar performance differences between the different channel combinations for coronal orientation scanning. Given the asymmetry in the coil design, with its lack of any head-foot segmentation, the manner in which changing the orientation of the scanning plane modulates the magnitude of the effect of the number of degrees of freedom exploited is likely to be strongly influenced by the properties of the coil. A systematic investigation of these effects in the context of an array configuration that has multiple elements arrayed in all three directions around the body will be of further interest and the current work suggests that such additional degrees of freedom would result in a powerful capacity for homogeneity control. Similarly, while this study was conducted purely at 3T, B 1 ? inhomogeneity issues are greater at higher field strengths. The results of this work cannot directly be translated to these other scenarios, but they do indicate that transmit arrays, perhaps with more than 8 channels, could mitigate related issues at field strengths greater than 3T.

Conclusion
The impact the number of RF transmit channels has on B 1 ? homogeneity, RF power requirements/SAR reduction and image quality when performing subject-specific RF shimming was investigated. Increasing the number of transmit channels was found to yield progressive improvements in B 1 ? homogeneity achieved under approximately constant total drive power, or reduction in total drive required to achieve the same level of B 1 ? homogeneity. The regularised shimming approach adopted in this study allows B 1 ? homogeneity and total RF power to be traded against one another and the balance adjusted throughout an examination to optimise imaging performance and control RF power deposition. RF shimming with large coil arrays allows for significant increases in homogeneity, beyond those found with 2-channels. In this study, the maximum number of channels was 8. Potentially, further performance increases could be achieved in some anatomies with a greater number of transmit channels, with the use of suitable array coils.