Automated Breast MRI Segmentation Method for Background Parenchymal Enhancement
Vignesh A Arasu1, Roy Harnish1, Cody McHargue1, Wen Li1, Lisa J Wilmes1, David Newitt1, Ella Jones1, Laura J Esserman2, Bonnie N Joe1, and Nola M Hylton1

1Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States, 2Surgery, University of California, San Francisco, San Francisco, CA, United States

Synopsis

Automated measurements of whole breast segmentation are becoming an essential process to the development of quantitative and reproducible imaging biomarkers. We have developed a method for automated whole breast tissue segmentation and assess its performance using a test dataset, and found approximately 75% of cases had satisfactory segmentation requiring none to minimal manual modification. The current method can likely provide accurate assessment of mean background parenchymal enhancement, but further refinement of breast-chest wall boundary identification is required for other measurements (e.g. breast density).

Purpose

Automated measurements of breast segmentation are becoming an essential process to the development of quantitative and reproducible imaging biomarkers. Background parenchymal enhancement has been suggested as a possible biomarker to assess breast cancer risk and for early prediction of neoadjuvant response. Manual processing is time consuming, prone to inter- and intra-rater reliability limitations, and therefore not practical for use as a regular clinical metric or for large-scale data imaging research. Background parenchymal enhancement represents the mean value of the initial enhancement of normal fibroglandular tissue, and because it represents an averaged value, requires only a large portion of tissue be segmented for accurate measurement. We have tested a method for automated breast tissue segmentation and assessed its performance using a pilot dataset.

Methods

In the IRB/HIPAA compliant study, breast tissue segmentation (Figure 1) was obtained using a coronal reformat of a pre-contrast MRI axial series resampled to a 1.0 mm3 isotropic voxel size. Image processing algorithms were applied using MATLAB and Statistics Toolbox Release 2012a (The MathWorks, Inc., Natick, Massachusetts, United States). Parallel imaging artifacts were eliminated and an initial tissue/air mask was obtained. The mask was traversed slice by slice from anterior to posterior and the number of 2D connected components in each slice was tracked. When the number of connected components dropped from 2 to 1 signaling that the sternal region was reached, an additional 1 cm band in the A/P direction was included and mask slices beyond that point were eliminated. The mask was then resampled into the original image matrix and b-splines were applied for manual editing using in-house software programmed in IDL (ITT Visual Information Solutions, Boulder, CO, USA).

The automated breast segmentation method was then applied to a test dataset of 49 patients with bilateral axial breast MRIs. These examinations were originally obtained in patients who received breast MRI to monitor neoadjuvant response. Two reviewers evaluated the segmented examinations qualitatively against the reference standard of manually drawn breast contours. Method accuracy was assessed if manual modification was required based on volume of breast tissue covered and extent of skin/chest wall inclusion to be sufficient for background parenchymal enhancement measurement.

Results

Of the 49 examinations tested, 61% of patients had fibroglandular tissue successfully segmented that required no further modification in order to measure background parenchymal enhancement (e.g. Figure 2a). A further 13% had satisfactory breast segmentation that only required minimal manual modification. The final 26% required significant manual adjustment either due to segmenting <50% of breast tissue or including skin/chest wall (e.g. Figure 2b). Of those that required significant manual adjustment, most were due to inadequate breast volume covered and required adjustment of the b-splines to extend posteriorly to the chest wall. A single case with breast reconstruction and implants required the most significant adjustment. Overall, segmentation cases captured more of the extreme craniocaudal sections of the breast than was obtained by manual segmentation, but performed worse along the posterior segmentation boundary.

Conclusion

Initial results suggests that this automated segmentation method was successful in approximately 75% of bilateral breast MRI examinations in our pilot dataset requiring minimal to no manual correction. The method in its current form would likely characterize background parenchymal enhancement accurately, and improved upon manual segmentation by including the extreme craniocaudal sections of the breast. Future versions will focus on better characterization of the breast-chest wall boundary and optimization with non-symmetric breasts such as post-reconstruction patients.

Acknowledgements

This work was supported by an NIH / NIBIB T32 EB001631 training grant.

References

King, V., Brooks, J. D., Bernstein, J. L., Reiner, A. S., Pike, M. C., & Morris, E. A. (2011). Background Parenchymal Enhancement at Breast MR Imaging and Breast Cancer Risk. Radiology, 260(1), 50–60. doi:10.1148/radiol.11102156

Preibsch, H., Wanner, L., Bahrs, S. D., Wietek, B. M., Siegmann-Luz, K. C., Oberlecher, E., et al. (2015). Background parenchymal enhancement in breast MRI before and after neoadjuvant chemotherapy: correlation with tumour response. European Radiology. doi:10.1007/s00330-015-4011-x

van der Velden, B. H. M., Dmitriev, I., Loo, C. E., Pijnappel, R. M., & Gilhuijs, K. G. A. (2015). Association between Parenchymal Enhancement of the Contralateral Breast in Dynamic Contrast-enhanced MR Imaging and Outcome of Patients with Unilateral Invasive Breast Cancer. Radiology, 276(3), 675–685. doi:10.1148/radiol.15142192

Figures

Figure 1. Coronal reformats were obtained and processed first to a) eliminate parallel imaging artifacts using distance map from a rough breast-air boundary. Subsequently, the b) 2D image mask is created and c) connected components added together. Finally, d) the 3D mask is created as the number of connected components vs. slice index going from anterior to posterior through the coronal slices.

Figure 2. Breast segmentation results were qualitatively assessed in a test dataset. An example is shown of a) successful segmentation of fibroglandular breast tissue would accurately measure background parenchymal enhancement and not require manual modification. b) Unsuccessful segmentation of < 50% of the fibroglandular tissue.



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