1786

Evaluation of Multicoil SENSE Reconstruction of Undersampled Data using a Human Observer Model of Signal Detection
Alexandra G O'Neill1, Tavianne M Kemp1, Sajan G Lingala2, and Angel R Pineda1
1Mathematics Department, Manhattan College, Riverdale, NY, United States, 2Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, IA, United States

Synopsis

We evaluated undersampling in MRI using a multicoil SENSE reconstruction with no regularization based the detection of signals by humans. We used a sparse difference-of-Gaussians (S-DOG) model to predict human performance in the detection of a small and large signal in anatomical backgrounds. The prediction was then validated using human observer two-alternative forced choice (2-AFC) tasks. Our model predicted a decrease in performance for both the small and large signal from 4X to 5X acceleration. Our observer study validated that prediction. This approach may lead to a way of assessing image quality that predicts human performance with fewer observer studies.

Purpose

Undersampling of magnetic resonance imaging (MRI) can be used to reduce the time required for a patient to be in the scanner. Ideally the undersampling would not result in a loss of diagnostic information. A surrogate for diagnostic information is the performance of human observers in simple detection tasks. We use task-based assessment of image quality1 by conducting human observer studies and implementing a sparse difference-of-Gaussians (S-DOG)2 observer model to evaluate the detectability of small and large signals in accelerated MRI using multicoil SENSE. The objective is to be able to predict human performance with fewer time-consuming human observer studies.

Methods

The multi-coil brain data used in this study were from fluid-attenuated inversion recovery (FLAIR) data sets from the fastMRI open source dataset3. All the data sets used had nky = 320, nkx = 640 and nc = 20. We carried out one-dimensional undersampling where we kept 16 k-space lines around the zero frequency (5% of the data) and everywhere else collected every kth line. For simplicity, we refer to this acquisition as a kX undersampling. The SENSE reconstructions were done in the BART4 environment with coil maps estimated via the sum of squares approach using the central 16 k-space lines5. In order to generate the images with the signals, the multi-coil k-space data of the small signal (disk with radius = 0.25 pixels blurred by a Gaussian kernel, σ=1 pixel) and the large signal (disk with radius = 3.5 pixels blurred by a Gaussian kernel, σ=1 pixel) was added to the background image data before undersampling and reconstruction. The amplitude of the signals was chosen to be on the threshold of being visible to better evaluate the effect of undersampling.

Three human observers were trained in performing two-alternative forced-choice (2-AFC) studies (Figures 1 and 2), in which the observer determines which sub-image (128 x 128) contains the signal. From these studies, the percentage of images the user identifies as correctly (PC) containing the signal is obtained and used for quantifying image quality. The 2-AFC trials were performed by all observers on a Barco MDRC 2321 monitor in a dark room. Observers were instructed to stay approximately 50 cm away from the monitor, a marker at 50 cm was provided for reference. For each condition every observer did 200 2-AFC trials from which their percent correct was obtained.

The S-DOG model uses frequency channels (Figure 3) and an internal noise model to model detection by humans2. The S-DOG model has also been used to predict the detection performance in human observer studies due to total variation regularization6. The level of internal noise was chosen to match the average human performance at 4X acceleration for both the small and large signal. For the small signal, a small amount of internal noise was needed (0.1% of background variability) but a larger amount of internal noise was needed for the large signal (13.6% of the background variability). All other accelerations were estimated with those resulting models.

Results & Discussion

The results for the small signal (Figure 4) show a prediction from the S-DOG model observer that the first large drop in performance in the detection task would occur between the 4X and 5X accelerations. This prediction was validated by a subsequent observer study. The predicted drop by the S-DOG model was 0.04 and the drop in average human performance was 0.07. Similarly, for the large signal (Figure 5), the model observer predicted the first large drop in performance between the 4X and 5X acceleration which was verified by an observer study. The predicted drop by the S-DOG model was 0.03 and the drop in average human performance was 0.04. While in both cases the S-DOG underestimated the drop in human performance, it provided an indication of the change occurring at those acceleration factors. As expected the acceleration affected the detection of the small signal more than the detection of the large signal. Future work will include the use of more human observers and test conditions to better understand the limits of both this model observer and the 2-AFC task used.

Conclusion

The S-DOG channel was able to predict a drop in human detection performance from the 4X to 5X acceleration which was validated by human observers.

Acknowledgements

We acknowledge support from NIH R15-EB029172, the Manhattan College Faculty Development Grant and the Kakos Center for Scientific Computing. The authors also thank Dr. Krishna S. Nayak, and Dr. Craig K. Abbey for their helpful insights.

References

1. HH Barrett and KJ Myers 2004 Foundations of Image Science (Hoboken, NJ: Wiley)

2. CK Abbey, HH Barrett, "Human-and model-observer performance in ramp-spectrum noise: effects of regularization and object variability", JOSA A, 18: 1237-1242 (1994)

3. J Zbontar, F Knoll, A Sriram, et al, "fastMRI: An Open Dataset and Benchmarks for Accelerated MRI", arXiv: 1811.08839v1 (2018)

4. M Uecker, F Ong, JI Tamir, et al, "Berkeley Advanced Reconstruction Toolbox", Proc. Intl. Soc. Mag. Reson. Med, 23, 2486 (2015)

5. K Pruessmann, M Weiger, M Scheidegger, P Boesiger, “SENSE:Sensitivity Encoding for Fast MRI”, Magn Reson Med, 42:952-962 (1999)

6. AG O’Neill, EL Valdez, SG Lingala, AR Pineda, “Modeling human observer detection in undersampled magnetic resonance imaging (MRI)”, Proc. of SPIE Medical Imaging 2021, 115990. 115990H.

Figures

Figure 1. Sample 2AFC trial with a small signal detection task where observers choose which image contains the signal (always in the center of the image). The signal is in the right image in this trial.

Figure 2. Sample 2-AFC trial with a large signal detection task and observers choose which image contains the signal (always in the center of the image). The signal is in the left image for this trial. Even though both images had the same level of undersampling, they may have come from different volumetric acquisitions. The image to the left has more ghosting artifacts.

Figure 3. The cross section of S-DOG channels in the frequency domain. The channels are inspired by the frequency sensitivity of the human visual system which requires higher contrast for low and high frequencies2.

Figure 4. The percent correct (PC) for the 2-AFC task with the small signal. The SDOG model observer was calibrated to match the average human performance at acceleration 4X and then used to predict all other accelerations. A following observer study at 5X is consistent with the predictions of the model observer.

Figure 5. The percent correct (PC) for the 2-AFC task with the large signal. This signal was chosen to be subtler than the small signal. The large signal is also less affected by acceleration than the small signal. The SDOG model observer was calibrated to match the average human performance at acceleration 4X and then used to predict all other accelerations. A following observer study at 5X is consistent with the predictions of the model observer.

Proc. Intl. Soc. Mag. Reson. Med. 30 (2022)
1786
DOI: https://doi.org/10.58530/2022/1786