Automated breast-region segmentation in the axial breast MR images

https://doi.org/10.1016/j.compbiomed.2015.04.001Get rights and content

Highlights

  • The method for the breast-region segmentation in the axial MR images is proposed.

  • It uses the shortest-path search incorporating information from the adjacent slice.

  • The cost function incorporates the edges obtained using a tunable Gabor filter.

  • The method is applicable for cases with no visible contrast at breast-region border.

  • The method requires no training procedure.

Abstract

Purpose

The purpose of this study was to develop a robust breast-region segmentation method independent from the visible contrast between the breast region and surrounding chest wall and skin.

Materials and methods

A fully-automated method for segmentation of the breast region in the axial MR images is presented relying on the edge map (EM) obtained by applying a tunable Gabor filter which sets its parameters according to the local MR image characteristics to detect non-visible transitions between different tissues having a similar MRI signal intensity. The method applies the shortest-path search technique by incorporating a novel cost function using the EM information within the border-search area obtained based on the border information from the adjacent slice. It is validated on 52 MRI scans covering the full American College of Radiology Breast Imaging-Reporting and Data System (BI-RADS) breast-density range.

Results

The obtained results indicate that the method is robust and applicable for the challenging cases where a part of the fibroglandular tissue is connected to the chest wall and/or skin with no visible contrast, i.e. no fat presence, between them compared to the literature methods proposed for the axial MR images. The overall agreement between automatically- and manually-obtained breast-region segmentations is 96.1% in terms of the Dice Similarity Coefficient, and for the breast-chest wall and breast-skin border delineations it is 1.9 mm and 1.2 mm, respectively, in terms of the Mean-Deviation Distance.

Conclusion

The accuracy, robustness and applicability for the challenging cases of the proposed method show its potential to be incorporated into computer-aided analysis systems to support physicians in their decision making.

Introduction

Magnetic Resonance Imaging (MRI) is an invaluable tool in the clinical work-up of patients suspected of having breast cancer [1], [2]. Moreover, the breast MRI is gaining popularity as a screening modality for patients with dense breasts or high-risk patients with BRCA 1 & 2 gene mutations (screened at a younger age when their breasts are dense) [3], [4]. As the diagnostic efficiency is highly dependent on the level of experience of the radiologist [5], [6], recent researches have focused on developing computer-aided analysis methods aimed at helping the radiologists in their diagnostic tasks, particularly those radiologists experienced in mammography reading but less experienced in interpreting the breast MRI [7].

The automated breast-region segmentation is required for performing a fully-automated computer-aided analysis of the breast MR images overcoming drawbacks of the manual and user-assisted segmentation which are impractical for processing large amounts of the MRI data being time consuming and biased by both the intra-observer and inter-observer variability. The breast-region segmentation refers to separation of the breast as an organ from other body parts in the MR images. The importance of the automated breast-region segmentation has recently been recognized in a number of applications for the breast MRI, such as the breast-density measurement [8] and tumour detection for the subsequent tumour classification [9]. It can also facilitate the pharmacokinetic-model calibration with respect to the reference tissues for improving the diagnostic performance of the dynamic contrast-enhancement breast tumours [10] where, for the purpose of calibration, the pectoral muscle (the most anterior chest-wall muscle) can be used given its properties [11]. This would require, among others, a precise breast-chest wall border determination which is a part of the breast-region segmentation.

To automatically segment the breast-region in the MR images, the segmentation algorithm needs to determine the breast-chest wall and breast-skin borders as well as the lateral-posterior breast ends. In the T1-weighted MR images – which are almost always included in the typical clinical breast MRI protocols for tumour analysis and used for the breast-density measurement – the fibroglandular tissue, skin and chest wall have a similar signal intensity. Therefore, for the challenging cases where a part of the fibroglandular tissue is connected to the chest wall and/or skin resulting in no visible contrast between them, automated breast-region segmentation is a highly demanding task. Moreover, the MRI signal intensity is usually affected by inhomogeneity, partial volume effect, aliasing, ghosting artefacts, etc., furtherly aggravating the breast-region segmentation.

There have been a few methods proposed in the literature to automatically segment the breast region from the MR images. The simplest ones are based on the intensity threshold usually followed by morphological operations [12], [13]. The more advanced ones, such as fuzzy C-mean (FCM) clustering [14], sign of gradients [9], [15], region growing [16], [17], and Markov random field [18], use various connectivity rules between pixels. These methods are highly dependent on the visible contrast between the breast region and chest wall/skin, i.e. on the presence of the fat along the breast-region border, and therefore tend to fail in the challenging cases.

The performance of the model-based [19] and atlas-based [20], [21] methods for the breast-region segmentation depends on the size and variety of the training database which is critical in achieving a reasonable accuracy. Moreover, the atlas-based segmentation techniques rely on an accurate intensity-based registration of the image to be segmented with the population atlas where the intensity-based registration assumes a similar grey-level distribution in the image compared to the atlas which it is impossible to account for the wide range of the clinical acquisition protocols (e.g., T1-weighted, T1-fat suppressed, T2-weighted), while it is impractical to create high-quality, expert-annotated atlases for every different acquisition protocol [19]. The absence of one or both breasts, as a result of the mastectomy surgery, might also challenge the model-based methods for the breast-region segmentation. Moreover, the previously proposed atlas-based methods have not been validated for the challenging cases while the proposed model-based method has failed in the cases when the contrast between the breast and chest wall is lower than the contrast between the chest wall and chest.

The edge-based methods [17], [22], [23] for the breast-chest wall border determination incorporate classical edge-detectors, such as Canny and Hessian, that are highly sensitive to the noise, low contrast and weak-textural edges affecting the accuracy of the breast-region segmentation. Moreover, the edge-based methods developed for the axial MR images [17], [23] rely on the presence of the fat along the anterior side of the chest wall (according to their authors) and, therefore, they are not applicable for segmentation of the challenging cases.

In this paper we present a novel fully-automated method for segmenting the breast region in the axial MR images using a tunable Gabor filter which sets its parameters according to the local MR image characteristics to detect non-visible transitions between different tissues having a similar MRI signal intensity. The breast-chest wall border is determined using a multi-slice approach where the border information from the adjacent slice is used to specify the border-search area in the current slice due to the high similarity in the chest-wall morphology (shape and inner structure) between two adjacent breast MRI slices. The method applies the shortest-path search technique within the border-search area by incorporating a novel cost function using the edges obtained with a tunable Gabor filter and subsequent extension of the resulting path along the fat up to the lateral-posterior breast ends. The lateral-posterior breast ends are determined using the body-contour landmarks, while the breast-skin border is determined by incorporating the edge information of the tunable Gabor filter in a slice-wise manner. Validated on 52 MRI scans covering the full range of American College of Radiology Breast Imaging-Reporting and Data System (BI-RADS) breast-density categories, our method shows to be independent from the visible contrast between the breast-region and surrounding chest wall and skin, i.e., independent from the presence of the fat along the breast-region border.

Our paper is organized as follows. After presenting the aim of our research and overview of the works related to the breast-region segmentation in Section 1, Section 2 describes the clinical breast MRI scans used in our research. Section 3 provides a detailed explanation of the proposed method for the breast-region segmentation. Section 4 introduces the validation metrics used to measure the segmentation accuracy. Section 5 reports our experimental results. Section 6 discusses the main findings of our research, outlines its limitations and determines the areas of our future work. Finally, Section 7 draws conclusions of our research.

Section snippets

Materials

The used database of 52 pre-contrast MRI breast scans (i.e. volumes) was obtained from the Institute of Oncology, Ljubljana, Slovenia. These were all clinical cases of different patients where the screening MRI revealed a lesion suspected of being malignant. The age of the patients included in our database ranged from 28 to 67 years with an average of 46 years. The patients were scanned in a prone position.

The axial T1-weighted images were acquired using a 3D fast low-angle shot-pulse sequence

Method

The proposed method for breast-region segmentation revolves around an edge-enhancement algorithm (see Section 3.1) originally designed for the fingerprint edge-enhancement in automatic fingerprint matching [25]. This algorithm performs a tunable Gabor filtering by setting the filter parameters according to the local MR image characteristics which is particularly suitable for detecting non-visible transitions between different tissues having a similar MRI signal intensity.

Validation metrics

To measure the breast-region segmentation accuracy of the proposed method, the agreement between the automatically- and manually-generated segmentation is assessed slice-wise within ROI using two validation metrics: the Dice Similarity Coefficient (DSC) [20], [22], [28], [29], [30] and the Mean-Deviation Distance (MDD) [17], [22]. It should be noted that the manual segmentation (described in Section 2) is performed within ROI to have the consistency in the lateral-posterior breast ends between

Results

The proposed method was implemented using MATLAB (v. R2009b; Mathworks, Natick, MA) requiring ~4.1 min running on a PC desktop (Intel Core CPU 3.67 GHz, 8 GB RAM) for processing one breast MRI scan (having a resolution of 448×448×144) compared to ~53 min required for the manual segmentation. This is a considerable increase in the time-efficiency requiring no user involvement. The speed of our method can be further increased by its implementation in C/C++. The edge-enhancement algorithm is performed

Discussion

In this paper we present a novel fully-automated method for the breast-region segmentation in the axial MR images. The method first determines the lateral-posterior breast ends and then the breast-chest wall and breast-skin borders where each step relies on the edge map obtained by the edge-enhancement algorithm based on a tunable Gabor filter. This filter preserves the edges in the local pixel neighbourhood having the same orientation as the orientation of the tuned filter and attenuating the

Conclusion

In this paper we present a fully-automated method for the breast-region segmentation in the axial MR images. The main advantage of our method compared to other methods proposed in the literature for the axial MR images is in its applicability for the challenging cases where a part of the fibroglandular tissue is connected to the chest wall and/or skin with no visible contrast, i.e. no fat presence, between them. The important advantage is also in not requiring any training set for learning the

Conflict of interest statements

None declared.

Acknowledgement

The authors are grateful to Dr. Bram Platel from Department of Radiology at the Radboud University Nijmegen Medical Centre, M.Sc. Kristijana Hertl MD and radiographer Gasper Podobnik from the Radiology department at the Institute of Oncology, Ljubljana, Slovenia, for their generous assistance in realization of this paper and for valuable discussions. The project was supported by the Slovenian Research Agency under the PR-03116-1 grant.

References (32)

  • T.E. Yankeelov et al.

    Quantitative pharmacokinetic analysis of DCE-MRI data without an arterial input function: a reference region model

    Magn. Reson. Imaging

    (2005)
  • G. Ertaş et al.

    Breast MR segmentation and lesion detection with cellular neural networks and 3D template matching

    Comput. Biol. Med.

    (2008)
  • Breast Imaging Reporting and Data System (BI-RADS), Breast Imaging Atlas, American College of Radiology,...
  • J.G.M. Klijn

    Early diagnosis of hereditary breast cancer by magnetic resonance imaging: what is realistic?

    J. Clin. Oncol.

    (2010)
  • S.G. Orel et al.

    MR imaging of the breast for the detection, diagnosis, and staging of breast cancer

    Radiology

    (2001)
  • A.J. Rijnsburger et al.

    BRCA1-associated breast cancers present differently from brca2-associated and familial cases: long-term follow-up of the Dutch MRISC screening study

    J. Clin. Oncol.

    (2010)
  • C.W. Piccoli

    Contrast-enhanced breast MRI: factors affecting sensitivity and specificity

    Eur. Radiol.

    (1997)
  • M.J. Stoutjesdijk et al.

    Variability in the description of morphologic and contrast enhancement characteristics of breast lesions on magnetic resonance imaging

    Invest. Radiol.

    (2005)
  • H.P. Chan et al.

    Improvement of radiologists’ characterization of mammographic masses by using computer-aided diagnosis: an ROC study

    Radiology

    (1999)
  • A. Gubern-Mérida et al.

    Volumetric breast density estimation from full-field digital mammograms: a validation study

    PLoS ONE

    (2014)
  • J. Yao et al.

    Breast tumor analysis in dynamic contrast enhanced MRI using texture features and wavelet transform

    IEEE J. Select. Top. Signal Process.

    (2009)
  • D.A. Kovar et al.

    A new method for imaging perfusion and contrast extraction fraction: input functions derived from reference tissues

    J. Magn. Reson. Imag.

    (1998)
  • T. Twellmann et al.

    An adaptive tissue characterization network for model-free visualization of dynamic contrast-enhanced magnetic resonance image data

    IEEE Trans. Med. Imaging

    (2005)
  • K. Nie et al.

    Development of a quantitative method for analysis of breast density based on three-dimensional breast MRI

    Med. Phys.

    (2008)
  • V. Giannini, A. Vignati, L. Morra, D. Persano, D. Brizzi, L. Carbonaro, A. Bert, V. Sardanelli, D.A. Regge, A fully...
  • W. Qi, D. Hui, W. Guang-zhi, A multistage processing procedure for 4D breast MRI segmentation, in: Proceeding of IEEE...
  • Cited by (28)

    • Segmentation of breast MR images using a generalised 2D mathematical model with inflation and deflation forces of active contours

      2019, Artificial Intelligence in Medicine
      Citation Excerpt :

      However, before further analysis can be done, it is crucial for the computer algorithm itself to capture regions of interest accurately. In breast CAD, breast segmentation is an important first pre-processing step to speed up the subsequent processes without losing any important anatomical information [8–16]. For example, breast segmentation (estimating the breast boundary and chest wall boundaries) removes unnecessary regions such as the heart, lung and liver.

    • Automatic and fast segmentation of breast region-of-interest (ROI) and density in MRIs

      2018, Heliyon
      Citation Excerpt :

      However, the accuracy depends upon the size and variability of the database and requires an atlas that is representative of the population, which is computationally expensive. An edge based approach was proposed by [28] that is independent from the visible contrast between the breast ROI and chest wall. This method calculates cost function using edge information obtained from tunable Gabor filter.

    • Clinical applicability and relevance of fibroglandular tissue segmentation on routine T1 weighted breast MRI

      2017, Clinical Imaging
      Citation Excerpt :

      Comparison between quantified FGT and radiologist qualitative FGT assessment suggests a need for incorporation of FGT quantification into clinical practice. The FGT segmentation method described in this study includes only one non-automatic step, namely the use of semi-automated chest wall contouring to first isolate the breasts from the chest, as in several other studies [10–13]. Fully automated chest wall detection algorithms are being evaluated to automate this step, but with noteworthy limitations.

    View all citing articles on Scopus
    View full text