 Research
 Open Access
 Published:
Lesion segmentation in breast ultrasound images using the optimized marked watershed method
BioMedical Engineering OnLine volume 20, Article number: 57 (2021)
Abstract
Background
Breast cancer is one of the most serious diseases threatening women’s health. Early screening based on ultrasound can help to detect and treat tumours in the early stage. However, due to the lack of radiologists with professional skills, ultrasoundbased breast cancer screening has not been widely used in rural areas. Computeraided diagnosis (CAD) technology can effectively alleviate this problem. Since breast ultrasound (BUS) images have low resolution and speckle noise, lesion segmentation, which is an important step in CAD systems, is challenging.
Results
Two datasets were used for evaluation. Dataset A comprises 500 BUS images from local hospitals, while dataset B comprises 205 opensource BUS images. The experimental results show that the proposed method outperformed its related classic segmentation methods and the stateoftheart deep learning model RDAUNET. Its accuracy (Acc), Dice similarity coefficient (DSC) and Jaccard index (JI) reached 96.25%, 78.4% and 65.34% on dataset A, and its Acc, DSC and sensitivity reached 97.96%, 86.25% and 88.79% on dataset B, respectively.
Conclusions
We proposed an adaptive morphological snake based on marked watershed (AMSMW) algorithm for BUS image segmentation. It was proven to be robust, efficient and effective. In addition, it was found to be more sensitive to malignant lesions than benign lesions.
Methods
The proposed method consists of two steps. In the first step, contrast limited adaptive histogram equalization (CLAHE) and a side window filter (SWF) are used to preprocess BUS images. Lesion contours can be effectively highlighted, and the influence of noise can be eliminated to a great extent. In the second step, we propose adaptive morphological snake (AMS). It can adjust the working parameters adaptively according to the size of the lesion. Its segmentation results are combined with those of the morphological method. Then, we determine the marked area and obtain candidate contours with a marked watershed (MW). Finally, the best lesion contour is chosen by the maximum average radial derivative (ARD).
Background
According to the 2020 global cancer data report, breast cancer ranks first among the three most common cancers in women, indicating that it has become a serious threat to the health of women worldwide [1]. Studies show that early detection and diagnosis of breast cancer can effectively increase the cure rate [2]. At present, digital mammography (DM) and breast ultrasound are the two main tools used in breast screening in China. However, ultrasound has no ionizing radiation and can show the anatomy and pathology of dense breast tissue, which DM cannot achieve. Therefore, ultrasound is more suitable for detecting breast lesions in Asian women with high density than DM and is becoming a popular screening tool for breast cancer [3]. Geisel et al. also demonstrated the effectiveness, practicability and feasibility of breast ultrasound as a screening tool for the early detection of occult breast cancer [4]. However, in the process of breast ultrasound imaging, the speckle noise generated by coherent waves greatly reduces the image quality, which requires a high degree of professionalism for radiologists to address. Due to the lack of radiologists in remote areas, ultrasoundbased breast cancer screening cannot truly be popularized.
With the development of artificial intelligence technology, computeraided diagnosis (CAD) systems based on medical images have made great achievements in cancer detection. In particular, the development of an ultrasoundbased breast cancer CAD system is impressive. It can realize intelligent screening and diagnosis. When the system receives realtime images, it can perform lesion detection, segmentation and diagnosis automatically. Its application will greatly alleviate the lack of radiologists. However, due to the inherent problems of breast ultrasound (BUS) images such as speckle noise and low contrast, the accuracy of lesion segmentation has not been effectively improved, which greatly affects the reliability of the diagnosis results. Thus, finding a stable and effective BUS image segmentation method is of great significance to promote the application and popularization of ultrasoundbased breast cancer CAD systems.
Therefore, we conducted this research. Focusing on solving the inherent problems of BUS images quality, we are committed to designing a stable and efficient image segmentation method. In recent years, many excellent image segmentation algorithms have emerged. Level set, first introduced in 1994 [5] and improved in 1995 [6], 2005 [7], 2012 [8] and so on, has proven to be very effective in image segmentation. However, it needs a great deal of time to solve partial deferential equations (PDEs), which is not very practical. To solve this problem, morphological snake (MS) was proposed. It uses morphological operations on a binary level set to approach the differential operators of a standard PDE [9]. It needs only numerical calculations, so MS is simple and fast. In terms of the field of BUS image segmentation, many scholars have used parameter deformable models and geometric deformable model technology [8, 10, 11]. However, to achieve ideal segmentation results, an appropriate initial tumour boundary or a precise edgebased stop function should be set in advance. Other researchers have used and improved graphbased segmentation methods, such as [12] and [13]. Boukerroui attempted to overcome the biggest drawbacks of the MRF model, i.e., a low optimization speed and local optimization [14]. In 2013, Zhao proposed the generalized fuzzy cluster method (FCM) with spatial information, which performs well in segmentation and has a rapid convergence speed [15]. FCM and the improved FCM algorithm were applied to lesion detection in BUS images [16, 17]. Since 2013, there have been an increasing number of segmentation methods with supervised and semisupervised learning, especially with the increasing popularity of deep learning, which has made great progress in solving the problem of BUS image segmentation. Supervised learning includes support vector machines, artificial neural networks (ANNs), and convolutional neural networks (CNNs), which have been applied to BUS image segmentation and have made great progress [18,19,20,21]. Zhuang proposed the RDAUNET model [22], which performs best on BUS image segmentation compared to other models. To date, deep learning models have proven to be the best way to perform image segmentation. However, they face some major problems, which are also the main bottleneck for further development. For example, the prediction result is not sufficiently robust. Robustness is the basic performance metric determining whether a model can be widely used [23,24,25,26]. Additionally, the model is not explainable, and training data are not sufficient. To solve these problems, a new approach has integrated visual saliency into a deep learning model for BUS image segmentation [27]. Attention blocks were introduced into a UNet architecture, which learns feature representations that prioritize spatial regions with high saliency levels and achieved a Dice similarity coefficient (DSC) of 90.5% on a data set of 510 images. However, this method relies greatly on the quality of the saliency maps, which is also not sufficiently robust.
Therefore, how can we have an efficient and robust BUS image segmentation method? We again turned to some excellent classic segmentation methods. It has been reported that tomography watersheds have a certain effect on solving complex segmentation problems and are more stable than other existing methods, but they are sensitive to noise and might cause oversegmentation. In view of this, many scholars have improved watersheds. Huang and Chen [10] combined a watershed with the active contour model to obtain a relatively accurate tumour boundary. In 2009, Gomez used a marked watershed (MW) algorithm incorporating morphological techniques and an average radial derivative function [28]. The method was improved by using an anisotropic diffusion filter guided by texture descriptors derived from a set of Gabor filters and creating segmenting functions generated by Newton filters to facilitate more precise segmentation [29]. However, the anisotropic diffusion filter requires many iterations to obtain a good preprocessing result, which takes a long time. In addition, the acquisition of the marker function is slightly complex, which reduces the efficiency of the algorithm. In view of these two problems, we made improvements in a previous study [30]. We combined contrastlimited adaptive histogram equalization (CLAHE) and curvature filtering to preprocess the images and used a morphological method to obtain the marker function, which is simple and efficient. However, this method not only improves the segmentation accuracy and DSC but also brings a higher false positive rate (FPR), which means that many false positive tissues are also segmented. Therefore, to further improve the performance, this paper makes technical contributions that are summarized as follows:
1. We use CLAHE and a side window filter (SWF) to enhance the lesion contour and eliminate the influence of noise. Compared with some other preprocessing methods, it is the most beneficial to BUS image segmentation. We propose an embedded segmentation method, adaptive morphological snake (AMS). It is more robust and stable than MS when processing complex datasets with different sizes of lesions collected from different types of ultrasound equipment.
2. We propose an optimized marked watershed segmentation method, adaptive morphological snake based on marked watershed (AMSMW). Its marker region is corrected by AMS. Taking full consideration of the advantages of classical segmentation algorithms, such as the level set method [5], morphological snake (MS) [9] and MW [31], we find that AMSMW has higher segmentation precision and is 3–4 times faster than other existing methods.
Results
Evaluation metrics
We used both area and contour error metrics, which include the accuracy (Acc), true positive ratio (TPR), false positive ratio (FPR), Jaccard index (JI), Dice similarity coefficient (DSC), area error ratio (AER), Hausdorff error (HE), and mean absolute error (MAE), to evaluate dataset A. The calculation formulas of these indicators are listed below. In addition, we used the Dice coefficient (DC), areaundercurve (AUC), precision (PC), sensitivity (Sen), specificity (Sp), F1score (F1) and meanintersectionoverunion (MIOU) to evaluate the proposed method on dataset B. The calculation formulas of these indicators can be found in Zhuang’s paper [22]:
where \(A_{*}\) is the number of pixels in *, A is the number of pixels contained in the image, and the subscripts G and S represent the ground truth and segmentation result, respectively. C represents the contour of the ROI. z and k are the points in the contour.
Xian has noted how important these metrics are [12]. Large JI values and small AER, HE and MEA values indicate good performance. Supposing that JI is small, when AER, HE, and MAE are large, if TPR and FPR are both large, then the lesion was overestimated. If the TPR and FPR are both small, the lesion was underestimated.
Experiment details
First, we obtain 100 BUS images and 50 BUS images from dataset A and dataset B respectively as the research object of the research experiment on preprocessing method. And the remaining 500 images in dataset A and the remaining 205 images in dataset B are used as the test set for the comparison experiment. Then, we use Python to implement the algorithm and calculate the evaluation metrics. The parameter setting of the AMSMW method is shown in Table 1. Next, we introduce the process of the three experiments in detail.
Finding the most suitable BUS image preprocessing method
We explored the effect of preprocessing methods on the segmentation results by using four preprocessing schemes: SWF, CLAHE&SWF, CLAHE&CF&SWF, and CLAHE&CF to preprocess the 100 images from dataset A and the 50 images from dataset B. In addition, we set up a control group without preprocessing operations to determine the effectiveness of the preprocessing method on the segmentation result.
Comparing AMSMW with other classical image segmentation methods and deep learning methods
First of all, we used the best preprocessing method obtained in the previous step to preprocess the test set. Then, we compared the proposed method with some related segmentation methods and RDAUNET on the test set of dataset A. Finally, we compared the proposed method with the typical deep learning segmentation method on the test set of dataset B.
In terms of traditional segmentation methods, we implemented some related and classical methods that include level set [5], MS [9], MW [31], and FSMW [30]. For the MS method, we set its initialization position and radius to be the centre of the RROI and 70% of the smallest length and width of the RROI. Moreover, it should be noted that we used the same preprocessing method when performing the comparison experiments except for FSMW.
In terms of deep learning methods, many excellent image segmentation models have been borrowed, improved and used. In [22], several typical deep learning segmentation models were compared on dataset B. The results showed that the RADUNET model performs best. To make an objective comparison between AMSMW and RADUNET, we performed the following two experiments. In the first experiment, we used fivefold crossvalidation method. First of all, we redivided the test set of dataset A, according to the ratio of training set:validation set:test set is 6:2:2 and then used the angle transformation method to expand each training set by four times, and then used them to train and test RDAUNET. Finally, we obtained five segmentation results. In order to achieve a scientific and fair comparison, we tested AMSMW on the five test sets to obtain five results too. The second experiment is that we used AMSMW to conduct a segmentation experiment on the test set of dataset B, and compared the quantitative results with those published in the paper[22].
Studying the sensitivity of AMSMW to benign and malignant lesions
Benign and malignant tumours are very different in size, morphology, margins, and internal state, which may greatly affect the algorithm’s performance. If a relationship can be found, it will be of great significance to design a more adaptive BUS image segmentation algorithm. Therefore, we conducted an exploratory experiment on the algorithm’s sensitivity in segmenting benign and malignant lesions. The specific operation was to first group dataset A into 250 benign and 250 malignant groups. Then, a quantitative segmentation experiment was performed.
Experimental results
Combined with CLAHE, SWF can enhance the edge of the lesion and contribute to better BUS segmentation results
Quantitative results and some examples of qualitative results are shown in Table 2 and Fig. 1, respectively. As shown in Table 2, the “\(\surd \)” means “used” and “\(\times \)” means “not used”. Separate (A) and separate (B) respectively represents the 100 images from dataset A and the 50 images from dataset B used for the preprocessing method experiment. The “Overall” column lists the average values of Dice. Observed form Table 2, we can draw the following three conclusions. First, the segmentation results using preprocessing methods are much better than those without preprocessing methods, and different preprocessing methods improve the segmentation performance to different degrees, which shows that it is very necessary to choose a suitable preprocessing method. Second, the Dice of the CLAHE&SWF method is similar to that of CLAHE&CF&SWF, but the average value of the CLAHE&SWF method is slightly higher, indicating that the method is more stable on different source BUS image datasets. Third, it also shows that CLAHE can improve the contrast and SWF can smooth the noise and well preserve the lesion boundary. It can also be observed directly from Fig. 1. Compared with the images in the last two columns, the contrast of the images in the first three columns on the left which were preprocessed with CLAHE is significantly more balanced and stronger. Images in the third and fourth columns were treated with SWF and their lesion boundaries were clearly highlighted. However, although the images in the first column were also preprocessed by SWF, their lesion boundaries became blurred after using CF.
AMSMW performs best on both quantitative results and qualitative results
First, some relevant and excellent traditional segmentation methods were tested on dataset A, and the quantitative and qualitative results are shown later. It can be observed from Table 3 that even without the preprocessing method, MW still performs the best on TPR, indicating that MW can segment the entire lesion area more sensitively and comprehensively, whereas the error rate is also the highest, causing the FPR to be too high. This means that a large part that does not belong to the lesion area is also segmented, which can be seen intuitively from the qualitative result in Fig. 2. Taking the images in the third and fourth columns as examples, many normal tissues were segmented by MW. However, level set has a lower FPR than MW, indicating that it can effectively improve the ability to identify lesion contours. However, the values of the other indicators of level set were relatively low, indicating that it cannot find the entire lesion area. In addition, it can be seen from the list of standard deviations in Table 3 that the dispersions are generally small, which shows that the data distribution is appropriate and that the experimental results are credible.
Compared with level set, MS has the greatest advantage in that it uses morphological methods to replace the process of solving numerical differential equations, greatly improving the efficiency. Experiments on ordinary notebooks show that it takes 6 s for MS and 17 s for level set to segment an image. At the same time, it can be inferred from the quantitative results that MS is better than level set in most metrics, indicating that it can segment tumours more precisely than level set. However, as shown in Fig. 2, taking the third and fourth columns as examples, there are some lesions with many calcification points inside, which cannot be segmented completely by the MS method. Compared with MS, although AMS is slightly lower than MS on Acc, FPR, AER, HE and MAE, it has obvious advantages on TPR, DSC and JI, showing that AMS is more stable and can obtain more complete lesion. As shown in Fig. 2, images in the fourth row are segmentation results of AMS. It can be observed clearly that although parts of the surrounding normal tissue were mistaken for part of the tumour by AMS, all parts of the tumour were completely included, which is of great significance for AMSMW to obtain a complete marked area later.
By comparing MS+MW with AMSMW, we find that AMS has great effects on improving performance. The level set+MW algorithm is a method in which level set is used instead of AMS as the embedded segmentation method. By comparing the quantitative results of level set+MW with other methods, it can be found that AMSMW has obvious advantages on indicators other than TPR. In addition, as observed from Fig. 2, we find that the qualitative result of AMSMW is much closer to that of GT. Taking the images in the sixth row and fourth column as an example, AMSMW can not only perfectly resist the interference of calcification points and speckle noise in the lesion area, but also resist the interference of back echo of the lesion, so as to accurately identify tumour boundaries. And reduce FPR as much as possible. Moreover, the AMSMW method runs fastest. In summary, we believe that AMSMW has the highest efficiency and effectiveness.
Second, we performed segmentation on the 205 images from dataset B, using AMSMW. Furthermore, We performed fivefold crossvalidation experiments on the 500 images from dataset A, using RDAUNET and AMSMW respectively. The quantitative and qualitative results are shown in Table 4 and Fig. 3 and in Table 5 and Fig. 4, respectively. As observed from Table 4, AMSMW is slightly lower than RDAUNET on SP, PC and MIOU but obviously higher on the other five metrics. This indicates that AMSMW has good adaptability and can segment lesion areas more precisely than RDAUNET, which can also be seen more intuitively from Fig. 3. We take images in the first, eighth and tenth column for example. They either have calcification points inside the tumor, severe speckle noise, or strong echoes behind the tumor. Facing with so much interference, AMSMW can still accurately identify tumour boundaries without causing excessive segmentation. RDAUNET can also find the lesion area, whereas it segments more normal tissue area, increasing the FPR. It can be observed from Fig. 4 clearly, there are relatively more serious problems of oversegmentation and false positive in the segmentation results of RDAUNET. In addition, observed from Table 5, RDAUNET does not show a strong generalization ability because it outperforms AMSMW only on Sp and Pc. And it is far worse on the other seven metrics. Overall, it is obvious that if there are not enough training data and excellent hardware resources, even the best deep learning model is not available in regard to new or more complex data. Therefore, there is still a long way for deep learning to go to improve its generalization performance. In other words, although traditional algorithms cannot be fully automated, the semiautomated capability is sufficient to alleviate the burden on doctors, and excellent traditional segmentation algorithms still have good performance when processing complex data.
AMSMW is more sensitive to malignant tumours than benign tumours
As shown in Table 6, the performance of AMSMW in segmenting benign and malignant tumours is comparable. From the four indicators, TPR, DSC, JI and AER, we can conclude that AMSMW is more sensitive to discriminating malignant tumours. However, due to the strong echo behind the lesion and the possible strong inner calcification points, AMSMW has a high FPR when segmenting malignant lesions, causing the algorithm to perform poorly on the other four indicators. Therefore, AMSMW is more stable when segmenting benign lesions.
Discussion
The proposed method consists of two parts: an image preprocessing scheme and an image segmentation algorithm. These two parts are complementary and inseparable. The goal of the preprocessing scheme is to highlight the boundary of the lesion to obtain more accurate segmentation results. Some preprocessing methods can suppress noise and improve contrast but cannot highlight the boundary. However, some preprocessing methods are the opposite. Therefore, it is necessary to find a suitable image preprocessing scheme. By exploring the effect of preprocessing methods on the segmentation results, we find that for BUS images, CLAHE & SWF is a better image processing method. Theoretically, SWF can preserve boundaries well, while CLAHE can suppress noise and improve contrast. The combination of the two methods is especially suitable for breast cancer ultrasound images with low contrast and speckle noise. In image segmentation, considering the complexity and particularity of BUS images, our goal is to find a robust and efficient lesion segmentation method. MW is very robust in solving complex image segmentation problems. However, its accuracy largely depends on the accuracy of the “marked area”. Theoretically, the “marked area” is the known lesion area. The more accurate the “marked area” is, the higher the accuracy of the algorithm will be. Thus, we mainly optimized MW by improving the method of obtaining the “marked area”. Our idea is to find an excellent algorithm as the acquisition method for the “marked area”. AMS is an improved MS method that can adaptively change the working parameters without tedious calculation of the PED equation. It is proven to be very suitable for acquiring a “marked area” for MW. The experimental results show that the proposed algorithm achieves better segmentation results. Moreover, by comparing it with some other classic traditional segmentation methods, we find that the proposed method is the most efficient and effective algorithm.
In addition, we compared the proposed method with the stateoftheart deep learning models RDAUNET and UNET. It still performed well on most of the metrics on both dataset A and dataset B. Because it does not need a training set, it does not depend too greatly on data sets with different data distributions. Therefore, theoretically, its generalization performance should be better than that of deep learning algorithms. Therefore, in the present era of deep learning, which has attracted much attention and is widely sought after, most deep learning models do not have an ideal generalization ability, which leads to a bottleneck in its continued development. However, the traditional segmentation method is stable and efficient, which provides a way to solve the bottleneck problem of deep learning to a certain extent. Therefore, in future work, we should not neglect classical segmentation methods. Most likely, it would be a good solution to integrate the efficient and stable traditional segmentation method with a deep learning model to complete image segmentation work, and the results would certainly be greatly improved.
By evaluating the sensitivity of the algorithm in segmenting benign and malignant tumours, we find that the proposed method has high sensitivity in the delineation of malignant tumour boundaries and is relatively stable for benign tumours. Therefore, in future work, we could take the strong echo and characteristics of malignant tumours into account and set up an adaptive ideal segmentation method.
Moreover, the current study has great potential for further studies. It will result in better and faster precise diagnosis and treatment of oncological diseases. There are many more medical specialities and diseases [32,33,34] in which there are known and applicable diagnostic imaging methods, but there are still few predictive modelling bases. Noninvasive diagnostic methods can be used not only in oncology but also in other medical specialities. Therefore, this study can also be applied to computeraided diagnosis of these diseases.
Conclusions
In this paper, an efficient semiautomatic BUS image segmentation method was proposed and evaluated quantitatively. It was proven to be the most robust and effective BUS image segmentation method compared with classic traditional segmentation methods and a stateoftheart deep learning model. In addition, by studying the sensitivity of AMSMW in segmenting benign and malignant lesions, we found that it is more sensitive to malignant lesions and more stable to benign lesions, which is of great significance for algorithm research in precision medicine in the future. Moreover, since the RROI in the proposed method is drawn manually, we are considering adding a deep learning network to automatically identify RROIs and completely liberate radiologists from this task in our future work.
Methods
The flowchart of the proposed method is shown in Fig. 5. It consists of five main parts: data acquisition, rectangular region of interest (RROI) acquisition, image preprocessing, marked area acquisition and final contour acquisition. The first three parts are image preparation and preprocessing. The last three parts are the process of image segmentation.
Data acquisition and difference analysis of the two datasets
Dataset A was collected by us. It has 600 BUS images, which include 300 benign solid cysts and 300 malignant solid cysts. They are captured from different devices, such as GE LOGIQ E9 and PHILIPS EPIQ 5, in a local hospital. The patient information in all images was hidden. An experienced radiologist sketched the lesion boundary for each image as the ground truth (GT). Dataset B is open source. It contains a total of 255 images. Among them, 213 images are from [22], and 42 images are from [35]. To study the generalization ability of the algorithm on different datasets, we analysed the significant differences between dataset A and dataset B. We used a grey level cooccurrence matrix to extract the following statistics for each image: the difference entropy, sum entropy, correlation, angular second moment, sum average, contrast, difference variance, entropy, homogeneity, sum variance, variance and information measure of correlation. Then, we used the MannWhitney U test to perform statistical analysis on datasets A and B to obtain the pvalue of each statistic, as shown in Table 7. We can see that the pvalues of all statistics are less than 0.05. Therefore, we find that datasets A and B have significant differences.
RROI acquisition
The RROI was obtained manually by the following steps: first, a point was selected as a starting point by leftclicking the mouse and holding down the left mouse button to move diagonally until the end position was found. Here, we define the RROI’s four vertices as \(w_{1}\),\(w_{2}\), \(h_{1}\), and \(h_{2}\), where w and h represent points along the tumour’s width and height, respectively, and the subindices 1 and 2 represent the lower and upper limits of the tumour’s width and height, respectively. The geometric centre of the RROI, which will be used later, can be defined as
Image preprocessing
Contrast enhancement
BUS images are characterized by low contrast and considerable noise, which can be improved by applying CLAHE, an optimization method based on adaptive histogram equalization (AHE) that limits the increase in contrast. It effectively overcomes the problem of overamplifying noise in the AHE algorithm.
Edge highlighting
Local windows, whose centres align with the pixels being processed, usually cause blurred edges. To avoid this, [36] proposed SWF, which can significantly preserve edges. Thus, we introduced SWF to highlight the edges of lesions in BUS images. We give a brief introduction to SWF, and more information can be found in [36].
As shown in Fig. 6, eight side windows are defined only in a discrete case, where (x, y) are the coordinates of target pixel i and r and \(\theta \) are the radius and angle of the window, respectively. \(\rho \in \{\hbox {o},\hbox {r}\}\) , \(\theta =\hbox {k}\times \pi /2\), and k\(\in \)[0,3]. Thus, we can obtain four side windows, \(W_{Di}\), \(W_{Ri}\), \(W_{Ui}\) and \(W_{Li}\), by setting \(\rho \)=r, which aligns i with their sides. While \(\rho \)=0, we have \(W_{SWi}\), \(W_{SEi}\), \(W_{NEi}\) and \(W_{NWi}\), which align i with their corners. For each pixel, the process of filtering can be regarded as the process of finding the \(I_{am}\) value, which satisfies
where,
\(W_{ij}\) is the weight of pixel j, which neighbours pixel i, based on the filtering kernel F; \(q_{j}\) is the intensity of image q at location i; and S=L, R, U, D, NW, NE, SW, SE is the set of side window indices. The result of filtering by SWF is defined as
where
\(I_{i}^{'\theta ,\rho ,\gamma }\) is the result for the eight side windows when \(\rho \in \{\hbox {o},\hbox {r}\}\), \(\theta =\hbox {k}\times \pi /2\), and k\(\in \)[0,3].
Constrained Gaussian kernel set
Similar to the method proposed in [30], we multiply Gaussian functions with the filtered image \(I_{CF}\) to obtain the region of interest (ROI). However, the difference is that we use a union of five constrained Gaussian distributions that have the same variances to make the lesion area more prominent:
The only difference between the five constrained Gaussian functions is the centre position. One is centred at the geometric centre of the RROI, and the other four are translated by half of the diagonal lengths in the four diagonal directions of the RROI. Hence, taking the Gaussian function centred at the geometric centre of the RROI as an example, its function can be expressed as
where \(\hat{p}\)(m,n) represents the pixel’s location and \(s_{\alpha }\) is the diagonal covariance matrix. This can be expressed as
We superimpose these five Gaussian functions to obtain their union \(G_{T}\) and then multiply it with \(I_{CF}\), which was negative previously:
Therefore, a specific highlighted ROI, whose surrounding tissue is greatly darkened, is obtained, as shown in Fig. 7. The experiments show that the illuminated ROIs obtained by these five Gaussian function sets are more complete than before, which is significant for determining accurate tumour boundaries and performing efficient segmentation.
Then, the input image of the marker function, JmnM, and the input image of the segmentation function, JmnN, are obtained separately by performing the opening operation on J(m,n) using a 9 and a 15pixel radius disk, respectively.
Marked area acquisition
The MW algorithm depends greatly on the marked area. The proposed method mainly improves the method of obtaining the marked area. We obtain the marked area by taking the intersection of the marker function and segmentation function.
Marker function
Similar to [30], we obtain the marker function through a series of morphological operations. First, we binarize JmnM with 1255 as the threshold. The 255 binarized images are denoted as \(f_{p}^{th}\)(m,n), (th=1,2,...,255), corresponding to 255 marker functions. Referring to Eq. 20, the marker function can be obtained by performing morphology operations:
where \(f_{p}^{th}\) represents the marker function and \(f_{ex}^{th}\) and \(f_{in}^{th}\) represent the external and internal markers, respectively. \(\delta \) and \(\varepsilon \) are morphological dilation and erosion, respectively, and \(B_{1}\) and \(B_{2}\) are two structural elements with a 15pixelradius disk and a 15pixelwide square, respectively.
Segmentation function
In [30], we discussed and proved that the segmentation function plays a large role in whether we can obtain accurate markers and makes a great contribution to obtaining good segmentation results. Therefore, to obtain more precise segmentation results, we evaluate the existing segmentation methods and propose an optimized method to obtain the segmentation function.
(1) MS: Let \(u_{e}\): \(R^{+} \times R^{2} \rightarrow \)R be an implicit representation of C such that C(t)=(x,y); u(t,(x,y))=0. MS uses a combination of binary morphological operators whose infinitesimal behaviour is equivalent to the flow expressed by the active contour PDE. Therefore, the curve is given as the zero level set of a binary piecewise constant function u: \(R^{2}\) \(\rightarrow \) \(\{\)0,1\(\}\). We take u(x)=1 for every point x inside the curve and u(x)=0 for every point x outside the curve. The morphological operators act on u and implicitly evolve the curve:
where v\(\in \) R is the balloon force parameter and g(I) selects which regions of I attract the curve. In the MS model, we use two common morphological operators: erosion and dilation. The dilation of a function is defined as
The erosion is defined as
The balloon force PDE can be expressed as
Given that the snake evolution at iteration n is \(u^{n}\): \(R^{2}\) \(\rightarrow \) \(\{\)0,1\(\}\), it can be solved using the following morphological approach:
where \(D_{d}\) and \(E_{d}\) are the discrete versions of dilation and erosion. Therefore, the morphological implementation of Eq. 27 can be expressed as
where \(SI_{d}\) and \(IS_{d}\) are smoothing operators. In a binary image u, \(SI_{d}\) works only on white pixels, and \(IS_{d}\) works only on black pixels. Taking \(SI_{d}\) as an example, for every white pixel \(x_{1}\) in a binary image, the \(SI_{d}\) operator looks for small (3pixellong) straight lines of white pixels that contain \(x_{1}\). This search is done in the four possible orientations corresponding to the four segments in P, where P is a collection of four discretized segments centred at the origin:
If no straight line exists, the pixel is made black (see Fig. 8). Sharp edges (Fig. 8b and d) are detected and removed as pixels that are not part of a straight line. White pixels on smooth edges (Fig. 8a and c) remain unchanged.
If no straight line exists, the pixel is made black (see Fig. 8). Sharp edges (Fig. 8b and d) are detected and removed as pixels that are not part of a straight line. White pixels on smooth edges (Fig. 8a and c) remain unchanged.
(2) AMS: Considering the different sizes of tumours, MS is not sensitive to especially large or small tumours. Therefore, we propose AMS, which is an optimized model based on the MS model, by applying adjustments to choose appropriate parameters.
In the AMS model, different shapes and types of tumour are considered. We use the geometric centre of the manually acquired RROI as the initial point and adjust the radius and iterations of the circle level set in real time according to the aspect ratio of the tumour. In the MS model, these parameters are fixed. Table 8 lists the relevant adjustable parameters of MS and AMS.
Finally, we can find the minimum boundary \(f_{sm}th(m, n)\), referring to Eq. 32. Then, we obtain the marked area by performing a closing operation with a 25pixelradius disk after binarization:
Final contour acquisition
First, we obtain 255 candidate contours by setting \(f_{label}^{th}\)(m, n) as the input of MW, referring to (33).
Second, we take the contour corresponding to the largest average radial derivative (ARD) value as the final contour. After the ARD is calculated for some sample images, 96 is determined as the average threshold value corresponding to the maximum ARD (for more details, please refer to [30]). To improve the efficiency of the algorithm and ensure that the selected boundary line is close to the ideal boundary line, we directly take the candidate boundary corresponding to the threshold of 96 as the final contour, thus avoiding the calculation of the ARD for 255 candidate boundaries in an image.
Availability of data and materials
The datasets analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 BUS::

Breast ultrasound
 CLAHE::

Contrast limited adaptive histogram equalization
 SWF::

Side window filter
 RROI::

Rectangular region of interest
 AMS::

Adaptive morphological snake
 MW::

Marked watershed
References
Siegel R.L., Miller K.D, Jemal A. Cancer statistics, 2020. Cancer J Clin. 2020;70(1):7–30. https://doi.org/10.3322/caac.21590.
Jemal A, Bray F, Center M.M, Ferlay J, Ward E, Forman D. Global cancer statistics. Cancer J Clin. 2011;61(2):69–90.
Drukker K, Giger ML, Horsch K, Kupinski MA, Vyborny CJ, Mendelson EB. Computerized lesion detection on breast ultrasound. Med Phys. 2002;29(7):1438–46. https://doi.org/10.1118/1.1485995.
Geisel J, Raghu M, Hooley R. The role of ultrasound in breast cancer screening: the case for and against ultrasound. Semin Ultrasound CT MR. 2018;39(1):25–34. https://doi.org/10.1053/j.sult.2017.09.006.
Sussman M. A level set approach for computing solutions to incompressible twophase flow. J Comput Phys. 1994;114(1):146–59. https://doi.org/10.1006/jcph.1994.1155.
Adalsteinsson D, Sethian JA. A fast level set method for propagating interfaces. J Comput Phys. 1995;118(2):269–77.
Shi Y, Karl WC. Realtime tracking using level sets. In: 2005 IEEE computer society conference on computer vision and pattern recognition (CVPR’05). USA: IEEE; 2005. vol. 2, p. 34–41. https://doi.org/10.1109/CVPR.2005.294.
Gao L, Liu X, Chen W. Phase and GVFbased level set segmentation of ultrasonic breast tumors. J Appl Math. 2012. https://doi.org/10.1155/2012/810805.
Álvarez L, Baumela L, Henríquez P, MárquezNeila P. Morphological snakes. In: Proceedings of the IEEE computer society conference on computer vision and pattern recognition. 2010. p. 2197–202. https://doi.org/10.1109/CVPR.2010.5539900.
Huang YL, Jiang YR, Chen DR, Moon WK. Watershed segmentation for breast tumor in 2D sonography. Int J Comput Assist Radiol Surg. 2006;1(SUPPL. 7):63–5. https://doi.org/10.1007/s1154800600446.
Gómez W, Infantosi AFC, Leija L, Pereira WCA. Active contours without edges applied to breast lesions on ultrasound. IFMBE Proc. 2010;29(2):292–5. https://doi.org/10.1007/978364213039773.
Xian M, Zhang Y, Cheng HD. Fully automatic segmentation of breast ultrasound images based on breast characteristics in space and frequency domains. Pattern Recognit. 2015;48(2):485–97. https://doi.org/10.1016/j.patcog.2014.07.026.
Chiang HH, Cheng JZ, Hung PK, Liu CY, Chung CH, Chen CM. Cellbased graph cut for segmentation of 2D/3D sonographic breast images. In: 2010 7th IEEE international symposium on biomedical imaging: from Nano to Macro, Netherlands: IEEE; 2010. p. 177–80. https://doi.org/10.1109/ISBI.2010.5490384.
Boukerroui D, Basset O, Guérin N, Baskurt A. Multiresolution texture based adaptive clustering algorithm fo breast lesion segmentation. Eur J Ultrasound. 1998;8(2):135–44. https://doi.org/10.1016/S09298266(98)000627.
Zhao F, Jiao L, Liu H. Kernel generalized fuzzy cmeans clustering with spatial information for image segmentation. Digit Signal Process. 2013;23(1):184–99. https://doi.org/10.1016/j.dsp.2012.09.016.
Lo C, Shen YW, Huang CS, Chang RF. Computeraided multiview tumor detection for automated whole breast ultrasound. Ultrasonic Imaging. 2014;36(1):3–17. https://doi.org/10.1177/0161734613507240.
Moon WK, Lo CM, Chen RT, Shen YW, Chang JM, Huang CS, Chen JH, Hsu WW, Chang RF. Tumor detection in automated breast ultrasound images using quantitative tissue clustering. Med Phy. 2014. https://doi.org/10.1118/1.4869264.
Liu B, Cheng HD, Huang J, Tian J, Tang X, Liu J. Fully automatic and segmentationrobust classification of breast tumors based on local texture analysis of ultrasound images. Pattern Recognit. 2010;43(1):280–98. https://doi.org/10.1016/j.patcog.2009.06.002.
Huang SF, Chen YC, Woo KM. Neural network analysis applied to tumor segmentation on 3D breast ultrasound images. In: 2008 5th IEEE international symposium on biomedical imaging: from nano to macro. France: IEEE; 2008. p.1303–06. https://doi.org/10.1109/ISBI.2008.4541243.
Xu Y, Wang Y, Yuan J, Cheng Q, Wang X, Carson P.L. Medical breast ultrasound image segmentation by machine learning. Ultrasonics. 2019;91(March 2018):1–9. https://doi.org/10.1016/j.ultras.2018.07.006.
Ronneberger O, Fischer P, Brox T. Unet: Convolutional networks for biomedical image segmentation. In: International conference on medical image computing and computerassisted intervention. Springer; 2015. p. 234–41.
Zhuang Z, Li N, Joseph Raj AN, Mahesh VG, Qiu S. An rdaunet model for lesion segmentation in breast ultrasound images. PloS ONE. 2019;14(8):0221535.
Abedinia O, Zareinejad M, Doranehgard MH, Fathi G, Ghadimi N. Optimal offering and bidding strategies of renewable energy based large consumer using a novel hybrid robuststochastic approach. J Clean Prod. 2019;215:878–89.
Saeedi M, Moradi M, Hosseini M, Emamifar A, Ghadimi N. Robust optimization based optimal chiller loading under cooling demand uncertainty. Appl Thermal Eng. 2019;148:1081–91.
Gao W, Darvishan A, Toghani M, Mohammadi M, Abedinia O, Ghadimi N. Different states of multiblock based forecast engine for price and load prediction. Int J Electr Power Energy Syst. 2019;104:423–35.
Khodaei H, Hajiali M, Darvishan A, Sepehr M, Ghadimi N. Fuzzybased heat and power hub models for costemission operation of an industrial consumer using compromise programming. Appl Thermal Eng. 2018;137:395–405.
Vakanski A, Xian M, Freer PE. Attentionenriched deep learning model for breast tumor segmentation in ultrasound images. Ultrasound Med Biol. 2020;46(10):2819–33.
Gómez W, Leija L, Pereira WCA, Infantosi AFC. Morphological operators on the segmentation of breast ultrasound images. In: 2009 Pan American health care exchanges. Mexico: IEEE; 2009. p. 67–71. https://doi.org/10.1109/PAHCE.2009.5158367.
Gómez W, Leija L, Alvarenga AV, Infantosi AFC, Pereira WCA. Computerized lesion segmentation of breast ultrasound based on markercontrolled watershed transformation. Med Phys. 2010;37:82–95. https://doi.org/10.1118/1.3265959.
Shen X, Liu J, Li H, Sun H, Ma H. A novel lesion segmentation method based on breast ultrasound images. In: ACM international conference proceeding series. 2019; p. 32–8. https://doi.org/10.1145/3366174.3366176.
Lewis SH, Dong A. Detection of breast tumor candidates using markercontrolled watershed segmentation and morphological analysis. In: 2012 IEEE southwest symposium on image analysis and interpretation. IEEE; 2012. p. 1–4.
LyssekBoroń A, Wylęgała A, Polanowska K, Krysik K, Dobrowolski D. Longitudinal changes in retinal nerve fiber layer thickness evaluated using avanti rtvuexr optical coherence tomography after 23g vitrectomy for epiretinal membrane in patients with openangle glaucoma. J Healthc Eng. 2017;2017:4673714.
Chatterjee A, He D, Fan , Antic T, Yulei J. Diagnosis of prostate cancer by use of mriderived quantitative risk maps: a feasibility study. Am J Roentgenol. 2019;213(2):W66W75.
Krysik K, Dobrowolski D, Polanowska K, LyssekBoroń A, Edward AW. Measurements of corneal thickness in eyes with pseudoexfoliation syndrome: comparative study of different image processing protocols. J Healthc Eng. 2017;2017:1–6.
Yap MH, Pons G, Marti J, Ganau S, Sentis M, Zwiggelaar R, Davison AK, Marti R. Automated breast ultrasound lesions detection using convolutional neural networks. IEEE J Biomed Health Inform. 2017;22(4):1218–26.
Yin H, Gong Y, Qiu G. Side window filtering. In: Proceedings of the IEEE conference on computer vision and pattern recognition. 2019. p. 8758–66.
Acknowledgements
Not applicable
Funding
This research is supported by the Guizhou Province Science and Technology Project under Grant Qiankehezhicheng [2019] 2794.
Author information
Authors and Affiliations
Contributions
XS carried out the algorithm design and implementation and drafted the manuscript. HM participated in the design of the study and in coordination and contributed suggestions to complete the manuscript. RL implemented RDAUNET. HL contributed to discussions. JH participated in collecting data and delineating tumours. XW participated in the discussion, organizing the data and writing the formulas with LaTeX code. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of the First Hospital of China Medical University and was in accordance with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. All subjects gave written informed consent in accordance with the Declaration of Helsinki.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Shen, X., Ma, H., Liu, R. et al. Lesion segmentation in breast ultrasound images using the optimized marked watershed method. BioMed Eng OnLine 20, 57 (2021). https://doi.org/10.1186/s12938021008917
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12938021008917
Keywords
 Breast cancer
 Morphological snake
 Side window
 Ultrasound
 Watershed