Towards a Parameter-Free ESPIRiT: Soft-Weighting for Robust Coil Sensitivity Estimation
Siddharth Srinivasan Iyer1, Frank Ong1, and Michael Lustig1

1Electrical Engineering and Computer Sciences, University of California, Berkeley, Berkeley, CA, United States

Synopsis

ESPIRiT is a robust, auto-calibrating approach to parallel MR image reconstruction that estimates the subspace of sensitivity maps using an eigenvalue-based method. While it is robust to a range of parameter choices, having parameters that result in a tight subspace yields the best performance. We propose an automatic, parameter free method that appropriately weights the subspace using a shrinkage operator derived from Stein's Unbiased Risk Estimate. We demonstrate the efficacy of our technique by showing superior map estimation without user intervention in simulation and in-vivo data compared to the current default method of subspace estimation.

Introduction

ESPIRiT is a robust, auto-calibrating approach to parallel MR image reconstruction that estimates the subspace of coil sensitivity maps using an eigenvalue-based method. While it is robust to a range of parameter choices, having parameters that result in a tight subspace yields the best performance. We propose an automatic, parameter free method that appropriately weights the subspace using a shrinkage operator derived from Stein's Unbiased Risk Estimate. We demonstrate the efficacy of our technique by showing superior map estimation without user intervention in simulation and in-vivo data compared to the current default method of subspace estimation.

Theory

In ESPIRiT, a calibration matrix $$$(A)$$$ is constructed from auto-calibration signal (ACS) data. Its singular value decomposition (SVD) is used to characterize the signal subspace $$$(V_{||})$$$ and the noise subspace $$$(V_\perp)$$$.1

$$\text{SVD}(A) = USV^*, \text{ where } V = \left[V_{||},V_{\perp}\right]$$

A self-consistency operator $$$\mathcal{G}_q$$$ is derived from the signal subspace. Self consistency implies $$$\mathcal{G}_q \, x = x$$$ where $$$x$$$ is the image. The point-wise eigenvalue decomposition of $$$\mathcal{G}_q$$$ is done in the image domain and eigenvectors with eigenvalues close to one with a tolerance threshold $$$(\gamma)$$$ are said to span sensitivity maps.1

A hard threshold ($$$\lambda$$$) on the singular values is used to estimate the signal subspace $$$(V_{||})$$$. If $$$s_1 \geq s_2 \geq \dots \geq s_{\min \dim A}$$$ are the singular values of $$$A$$$,

$$V_{||} = VW, \text{ where } W = \text{diag}\left[\frac{\mathbb{1}(s_i > \lambda s_1)}{s_i}\right]$$

While the hard threshold works well in most cases, it can be sensitive to the choice of $$$\lambda$$$. It might include vectors from the noise subspace or null vectors from the signal subspace. Finding the right $$$\lambda$$$ is therefore crucial, but this problem is not convex. To overcome this, we replace the hard threshold with a soft threshold that weighs down the noise subspace. Soft thresholding is more robust to choice of threshold and it is possible to determine the optimal soft threshold $$$(\hat \lambda)$$$ using Stein's Unbiased Risk Estimate.2

Once the optimal soft threshold $$$(\hat \lambda)$$$ is determined, we define a new weighted subspace estimate.

$$\hat V_{||} =V\hat W, \text{ where } \hat W = \text{diag}\left[\frac{(s_i - \hat \lambda)_+}{s_i}\right]$$

This weighting operator $$$\hat W$$$ appropriately shrinks the norms of the singular vectors. This weighting of $$$V$$$ to estimate $$$V_{||}$$$ results in $$$\hat{\mathcal{G}}_q$$$ being related to $$$\mathcal{G}_q$$$ as follows,

$$\hat{\mathcal{G}}_q = W^* \mathcal{G}_q W$$

We assume the eigenvectors of $$$\hat{\mathcal{G}}_q$$$ span "weighted" sensitivity maps and hence the eigenvalues do not necessarily meet the "=1" eigenvalue condition.1 To accommodate this, the second threshold $$$(\gamma)$$$ is set to $$$0.98$$$ times the maximum eigenvalue of the eigenvalue decomposition of $$$\hat{\mathcal{G}}_q$$$.

There are different ways to estimate the noise standard deviation. Since we hypothesize that $$$A$$$ is low-rank, we expect that the last few singular values should be contributions from noise only. In order to estimate noise standard deviation ($$$\sigma$$$), we generate a zero-mean, unit-variance Gaussian noise calibration matrix and fit its last one-fourth singular values to the last one-fourth singular values of $$$A$$$.

Methods

We look at the ESPIRiT maps generated from the two methods and qualitatively compare them. We use two fully sampled knee-datasets acquired from a Discovery MR 750 GE Scanner.3 A kernel of dimension $$$[6 \times 6]$$$ and a calibration region of dimension $$$[24 \times 24]$$$ are used. We use ESPIRiT's default threshold parameters $$$(\lambda = 0.001, \gamma = 0.8)$$$ to generate one set of ESPIRiT maps. The weighting method described above is used to generate the other set. We project the data on the above two sets of maps and compare the results and their difference images. (The difference image is the null projection. It should have only noise and no signal.)

Results

Figures 1, 2, 3 and 4 exemplify our results. Even with data being acquired from the same scanner, we see that the default parameters can result in loose sensitivity maps as in the Figure 2 and tight maps as in Figure 4. Contrast this with the soft-weighting method which results in tight sensitivity maps in both cases. Looking at Figure 1 specifically, we see the tight sensitivity maps result in more noise in the difference image (1e) compared to difference image (1c), which is desirable. This is especially noticeable in regions surrounding the knee, where we expect no signal and only noise.

Discussion and Conclusion

The effects of weighting the signal subspace on the second threshold need to be further studied. Nevertheless, the soft-weighting of singular vectors is a promising step towards having an ESPIRiT implementation that is both robust to noise and completely auto-calibrating.

Acknowledgements

No acknowledgement found.

References

1. Uecker, Martin, Peng Lai, Mark J. Murphy, Patrick Virtue, Michael Elad, John M. Pauly, Shreyas S. Vasanawala, and Michael Lustig. "ESPIRiT-an Eigenvalue Approach to Autocalibrating Parallel MRI: Where SENSE Meets GRAPPA." Magnetic Resonance in Medicine 71, no. 3 (2013): 990-1001.

2. Candès, Emmanuel J., Carlos A. Sing-Long, and Joshua D. Trzasko. "Unbiased Risk Estimates for Singular Value Thresholding and Spectral Estimators." IEEE Transactions on Signal Processing 61, no. 19 (2013): 4643-657.

3. Epperson, Kevin, Anne Marie Sawyer, Michael Lustig, Marcus Alley, Martin Uecker, Patrick Virtue, Peng Lai and Shreyas Vasanawala. "Creation of Fully Sampled MR Data Repository for Compressed Sensing of the Knee" Paper presented at SMRT Conference, Salt Lake City, UT, 2013. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.402.206&rep=rep1&type=pdf.

Figures

Figure 1: (1a) is the root-sum-of-squares reconstruction of the original fully sampled dataset. (1b) is the resulting projection of the image onto the ESPIRiT maps generated using default parameters and (1c) is the difference image of the (1a) and (1b) times 5. (1d) is the resulting projection onto maps generated using the soft-weighting method and (1e) is the difference image between (1a) and (1d) times 5. Note that (1e) has more noise than (1c), which is desirable.

Figure 2: These are the sensitivity maps generated by the two methods from the dataset illustrated in Figure 1. (2a) shows ESPIRiT maps generated from using the default threshold parameters and (2b) shows maps generated from using the soft-weighting method. Notice how much more tighter the maps are when using the soft-weighting method.


Figure 3: Similar to Figure 1, (3a) depicts the root-sum-of-squares reconstruction of the original fully sampled dataset. (3b) is the resulting projection from using the default parameters and (3c) is the difference image times 5. (3d) is the resulting projection from using the soft-weighting method and (3e) is the difference image times 5. Note that now, both projections are similar.

Figure 4: Similar to Figure 2, these are the sensitivity maps estimated by the two methods from the dataset in Figure 3. (4a) shows ESPIRiT maps generated from using the default threshold parameters and (4b) shows maps generated from using the soft-weighting method. Notice how both sets of maps look similar and are tight.




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