Computational and Mathematical Methods in Medicine
Volume 2019 (2019), Article ID 5856970, 10 pages
https://doi.org/10.1155/2019/5856970

Research Article

Trichomonas Detection in Leucorrhea Based on VIBE Method

Xiaohui Du, Lin Liu, Xiangzhou Wang, Jing Zhang, Guangming Ni, Ruqian Hao, Juanxiu Liu, and Yong Liu

University of Electronic Science and Technology of China, MOEMIL Laboratory, School of Optoelectronic Information, No. 4, Section 2, North Jianshe Road, Chengdu 610054, China

Correspondence should be addressed to Xiaohui Du; xiaohuie@126.com

Received 24 October 2018; Revised 20 December 2018; Accepted 25 December 2018; Published 14 January 2019

Academic Editor: Luminita Moraru

Copyright © 2019 Xiaohui Du et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Trichomonas examination is one of the important items in the leucorrhea routine detection. And it cannot be recognized by still images because of the unstable morphology and unfixed focal location caused by motion characteristic. We proposed an improved VIBE algorithm. 6 videos (totally 1414 frames) are collected for testing. In order to compare the effects of the algorithms, we segment each frame artificially as ground truth. Experiments show that percentage of correct classification (PCC) achieves 88%. The proposed improved method can effectively suppress the false detection caused by the formed components such as epithelial cells in the leucorrhea microscopic image and the missed detection caused by the background model update during the movement. At the same time, improvements can effectively suppress smear and ghost areas. The algorithm proposed in this paper can be integrated into the leucorrhea automatic detection system.

1. Introduction

Trichomonas is an important indicator in leucorrhea routine (RT), and the detection of trichomoniasis can prevent a variety of vaginal diseases. At present, the way to detect trichomoniasis in hospitals is under a microscope by manual examination, which is inefficient. So, the microscopic examination algorithm of trichomoniasis is especially necessary for the automation and intelligence of leucorrhea RT. Trichomonas cannot be recognized by still images because of the unstable morphology and unfixed focal location caused by motion characteristic. Trichomonas are morphologically diverse under a microscope. Normally, trichomonas will have flagella and similar texture to white blood cells (WBC) as shown in Figure 1. Its morphological diversity of trichomonas restricts its morphological detection. In most cases, trichomonas may be attached to other cells in the leucorrhea sample, such as epithelial cells, and this diversity poses challenges to ways such as morphologically segmentation.

Figure 1: Trichomonas’ texture: (a) compared to WBC; (b) trichomonas attached to other cells.

Because the trichomonas in the fresh leucorrhea sample is still active, it is reflected in the mirror as a motility, and therefore the detection method of the moving target is no doubt a new idea for the detection of the trichomonas.

The frame difference method is the earliest algorithm used to extract the target area of motion, which has strong robustness to the scene containing the moving target. The operation speed is fast, but the method does not generally detect all pixels of the moving object completely, and the phenomenon of “cavity” is often found inside the detected moving object, so the method is applicable to the case of simple motion detection [1]. Yang and Guo [2] and Cheng and Wang [3] proposed the three-frame difference method to solve the vacancy phenomenon of single-frame difference method, but this method has a great influence on lighting and other factors, which is easy to cause error detection.

Background modeling method is the common detection method used in the moving target detection, which uses the median value of pixel value of continuous K frame image sequence as the background model. The Gaussian mixture model (GMM) is a background modeling method proposed by Stauffer, which can effectively suppress the false detection caused by background changes and accurately detect the moving target area [4]. However, experiments show that the recognition accuracy of GMM model is low in the trichomonas detection. The kernel density estimation method proposed by Elgammal et al. does not need to make any assumptions about the density distribution of the background, but uses the recent image sample information to accurately estimate pixel points by using the standard kernel function to extract moving objects [5]. The experimental results show that the method has good adaptability in complex outdoor scenes, but the performance to trichomonas detection is not good enough.

VIBE is Barnich’s moving target detection method based on pixel points, and it is named the visual background extraction method (VIBE) [6]. This algorithm directly selects a certain number of pixel values randomly according to certain rules for background modeling for each pixel point and then classifies the foreground and background of pixel points using Euclidean distance. The advantage of this algorithm is that it does not need to assume any probability model, and it can detect moving objects in the second frame of video sequence, and the calculation speed is fast. Droogenbroeck and Paquot improved the VIBE algorithm, which improved the PCC (percentages of correct classification) and reduced the PBC (percentage of bad classifications). Although this paper presents a good detection result compared with VIBE, these improvements greatly increase the computation and reduce the real-time performance [7]. Hofmann et al. proposed the pixel-based adaptive segmenter (PBAS) detection method firstly [8]. Based on the advantages of SACON and VIBE, this algorithm combines and improves the detection accuracy of the target further. However, this method is very sensitive to the weak flow of background impurities, resulting in the increase of impurity misdiagnosis.

For the detection of trichomonas, Hao et al. adopted the improved Kalman filtering method to extract trichomonas [9]. Experiments show that this method can adapt to the situation of rapid changes in light. Moreover, when the movement speed of the moving target is slow, or the movement frequency is not high, the accurate extraction of trichomonas cannot be carried out, and the omission ratio is high.

In this paper, we propose an improved VIBE background reconstruction method, with three different background reconstruction optimization methods. Experiments show that PCC achieves 88%, which can filter the misconduct and misdetection caused by impurities effectively. And the phenomena of tailing and ghosts are eliminated. The algorithm proposed in this paper can be integrated into the leucorrhea automatic detection system.

This article is organized as follows. Dataset preparation and the method are described in Section 2. Section 2.1 is the origin VIBE method and Section 2.2 is the improvement to Section 2.1. Section 3 introduces the experiment results and the discussion. Conclusions, limitations, and future work are provided in Section 4.

2. Dataset Preparation and Method

By collecting a fresh sample, smearing it onto a glass slide, and using a 40x objective lens mounted on a Motic BA210 integrated microscope, we can get the video dataset with a frame rate of 16. All sample collections come from the Fourth Affiliated Hospital of Nanchang University. And all methods were carried out in accordance with relevant guidelines and regulations (the Fourth Affiliated Hospital of Nanchang University regulations). The study has been approved by IRB and other agencies (Ethics Committee of drug clinical trial of the Fourth Affiliated Hospital of Nanchang University, SFYLL-PJ-2015-001), and the samples do not contain any identity and diagnostic information of the subjects. While sampling, doctors would verbally inform the patients. None of the samples are obtained from minor patients or without agreement. For the convenience of comparison, we manually labeled 6 video images, a total of 1414 frames for analysis (ground truth).

2.1. VIBE Algorithm

The background point in VIBE is the pixels that do not move, which should be segmented as 0, while the foreground point is the moving pixels, such as trichomonas in leucorrhea images, which should be segmented as 255. In the VIBE model, a background model is created for each pixel in the first frame. In each frame image, calculate the similarity between the pixels to be classified (foreground and background) and the background model. If the pixel value of a point belongs to the background point, it should be close to the sampling value of the background model sample set. The calculation method of similarity degree is the Euclidean distance between the pixel value and the N sample values (, , …, ) in the background model, which is less than the number of given threshold value (R) as shown in Figure 2.

Figure 2: The distribution of pixel value and N sampling values in the background model in two-dimensional color space.

is a color spatial sphere with as the center and radius R as the threshold. The background model with N sample values created for each pixel point is as follows:

Define as the intersection of the background model and , should be classified as background, if satisfies the following condition:

where is the counting function and is the number of samples contained in the intersection . is the threshold for whether it is a background or not. The sensitivity of the model is expressed as

When initializing, VIBE initializes the background model by means of single-frame initialization. The specific method is to select pixel values in the current pixel neighborhood randomly:

The background update of VIBE is random update. According to the update strategy of VIBE, the probability of a sample between time and is . Assuming that time is continuous and there is no memory loss during the selection of background, the probability can be

In the above equation, the sample life cycle in the background model is exponential decay, and the probability of sample retention in time is independent of time . VIBE uses secondary random sampling to realize the possibility of representing an infinite time window with finite samples. The basic idea is to reduce the frequency of background update and extend the life of samples in the background model. According to the principle that neighborhood pixels have similar temporary distributions, VIBE propagates background samples in the neighborhood. This neighborhood propagation mechanism ensures the spatial consistency of the algorithm and can effectively recover the background area covered by the foreground.

2.2. VIBE Improvement

2.2.1. Update of VIBE Background Model

The VIBE model uses three main update strategies for background updates:

(i)The first is the memoryless update, that is, if a pixel x is classified as a background, the algorithm randomly takes a point in the background model and updates the value of with the pixel value of x. This strategy can reduce the probability that the model’s value of the background model will be discarded during the update process.(ii)The second is the time subsampling of the model. The original background point of the model may be background points in consecutive multiple frames; it will reduce the efficiency of the operation if the background model is updated every frame. VIBE uses a random update strategy (adjust φ) to determine whether the background point is updated.(iii)The third is the update of the spatial domain, which can ensure that video does not produce extensive area error detection in jitter. The specific method is that the field set of the background point x is , u is randomly taken from , and the value of the background model of is updated to . This dissemination update mechanism also determines whether the domain points are updated with a certain probability.

The movement of trichomonas in the video is slow, and there is no significant trajectory. When the model is being updated, we adopt the update strategy with memoryless update and time subsampling. For the update of the neighborhood, we update the spatial neighborhood for the background point x on the image of each frame and discard the original dissemination mechanism according to the probability. This is because there are a large number of morphological components such as epithelial cells, bacilli, and cocci in the leucorrhea sample. Since the fresh leucorrhea samples are still active, these morphological components in the video will do a small amount of swimming. Therefore, for each background point x, there will be a high probability of swimming. If there is a certain probability to determine whether the field is updated, many background points on morphological components will be falsely detected as the foreground points.

In this paper, we update the domain point in an absolute update method, that is, if one sample point x is detected as the background point, then background point model of the field is fixedly updated. Randomly take a domain background point y and randomly update its background model to :

The probability that the neighborhood point is selected is , and the probability that the neighborhood point retains its background model is . With the continuation of time, the probability that the neighborhood point background model is not updated during t + dt can be written as

In summary, a point is detected as a background point, and then there are two main situations in which the background model is updated: one is to update with a probability of 1/φ during time subsampling. The other one is that the neighborhood of the point x has a background point y. While performing the update calculation, point x updates with a probability of 1/8, assuming that the number of neighborhood of y is 8:

2.2.2. Update of VIBE Neighborhood Background Model

It is mentioned in Section 2.2.1 that when the point x is detected as a background point, the background model of the neighborhood point will be updated. According to the principle of VIBE, when the point y is the foreground point, there is no need to update the background model of y. Therefore, it is considered that when the neighborhood point y of x is the foreground point at this time, the background model of y does not need to be updated. In judging whether y is a foreground point, the method of Section 2.2.1 will undoubtedly introduce more calculations. In order to simplify the foreground and background judgment of y, this paper proposes a simplified discriminant method.

Since the motion in the leucorrhea sample video changes slowly, it is considered to introduce a frame difference method to prefetch the foreground. The algorithm flow is as follows:

(1)Extract the foreground image by using the frame difference method:

where is the frame, that is, the current frame image, while the is the previous frame image. is the absolute value of the difference of the current frame image minus the previous frame image.(2)Binarize with the Otsu threshold method to obtain (3)Perform morphological hole filling on (4)Perform close operation to the with a circular structural unit with size 11 to connect the discrete binary points of the same component:

where is the circular structural unit.(5)Perform open operation to the with a circular structural unit with size 9 to filter the weak flowing impurities:

The resulting effect picture is shown in Figure 3.

(a)

(b)

Figure 3: The foreground image obtained by the frame difference method: (a) the original image I and (b) the foreground extraction picture .

If the point x is the background point, it is first determined whether the neighborhood point y is the background when updating the neighborhood background model, that is, whether the corresponding pixel point in Figure 3(b) is 0. The update strategy becomes a random selection of the background neighborhood points of x for updating

For the background point x, the probability that the neighborhood y will be updated is

where is the number of neighborhoods that satisfy the condition, that is, in is the background point (pixel value is 0).

2.2.3. Impurity Filtration

We believe that the point should be a background point if a region is detected as a foreground spot for consecutive K frames. Through this improved method, most of the formed components regions can be filtered. This is because the components such as epithelial cells are basically static except for the trichomonas. Although they are affected by the water flow, these impurities undergo a slight movement. If a pixel is detected as a foreground point by successive K frames, then we can force it to be a background point, as shown below:

where indicates the cumulative number of times the point is detected as the foreground point at the frame and is the pixel value of foreground image at the position.

In addition, we believe that the foreground point will not appear out of thin air, and the trichomonas must have existed in the field of vision or moved from outside of the vision to the field of vision. Therefore, if a frame of a pixel point x is detected as a background, the next frame is a foreground point, and there is no foreground point in the neighborhood of the previous frame of x, and then x is considered to appear in the field of view out of thin air as follows:

2.3. Main-Process Stream

According to the improved algorithm based on VIBE, the flow chart of total trichomonas detection is shown in Figure 4.

Figure 4: The flow chart of the improved VIBE method.

During initialization, N = 20, is 2, the radius R is 40 px, and the subsampling probability φ is 3.

3. Experimental Results and Discussion

3.1. Update of VIBE Background Model

The background model update strategy of Section 2.2.1 can effectively filter out regional misdetection caused by weak sway.

As shown in Figure 5, the 27th frame image in the video image and its foreground segmentation are compared. After the update strategy for the neighborhood point is changed to absolute update, it can effectively suppress foreground segmentation misdetection caused by weak disturbances such as epithelial cells, and the updated foreground effect map is closer to the standard foreground image. Experiments have shown that this improved strategy can increase the percentage of correct classification (PCC) by 10%.

(a)

(b)

(c)

(d)

Figure 5: Result of background model update improvement: (a) original image; (b) standard trichomoniasis foreground map; (c) VIBE original algorithm foreground map; (d) improved algorithm (Section 2.2.1) foreground map.

3.2. Update of VIBE Neighborhood Background Model

Only when the neighborhood y of the background point x is the background point, y is possible to be updated. This strategy can effectively avoid some of the foreground point x′ to be updated as a background point in the next frame detection, due to the background model update of the neighborhood y′ as shown in Figure 6.

(a)

(b)

(c)

(d)

Figure 6: Image effect of background model update improvement: (a) original image; (b) standard trichomoniasis foreground map; (c) VIBE original algorithm foreground map; (d) improved algorithm (Section 2.2.2) foreground map of this paper.

In Figure 6, the 45th frame in the video image was used for comparison. It was found that the trichomonas at the box in the original VIBE model produced a missed detection, and after the introduction of the strategy (Section 2.2.2) improvement, the phenomenon was improved. Experiments show that this strategy can increase the detection PCC by about 8%.

3.3. Impurity Filtration Test

Misdetection of the foreground point can be effectively improved by the two strategies in Section 2.2.3. If pixel consecutive K frames are detected as a foreground point, make it as the background point.

(a)

(b)

Figure 7: The impact of changes in K on PCC with (a) and without (b) strategy in Section 2.2.1.

As shown in Figure 7, the effect of filtering the formed component can be effectively achieved, and the change of K has little effect on the final PCC change with the strategy in Section 2.2.1. The average value of the calculated PCC reached more than 80%. In contrast, in the absence of the strategy described in Section 2.2.1, the average value of PCC is smaller, only about 69%, and with the increase of K, PCC decreases, reaching the best when K takes the minimum value of 2.

In addition, the introduction of the strategy described in Section 2.2.3 can effectively eliminate the ghost area, as shown in Figure 8.

(a)

(b)

(c)

Figure 8: The effect of the strategy in Section 2.2.3 on eliminating the ghost region: (a) the original image; (b) K takes infinity; (c) K takes 2.

3.4. Comprehensive Detection Effect

Many metrics can be used to assess the output of a background subtraction algorithm given a series of ground-truth segmentation maps. Barnich and Vanogenbroeck [6] use the percentage of correct classification (PCC) to evaluate

where accounts for the number of correctly detected foreground pixels; accounts for the number of background pixels incorrectly classified as foreground; accounts for the number of correctly classified background pixels; and accounts for the number of foreground pixels incorrectly classified as background.

However, the number of trichomonas in the leucorrhea sample image is generally small, and the area of the whole map is small. Therefore, the introduction of TN in PCC, the black background, will be the majority, so we have changed the PCC as the comparison of PCC is meaningless:

In addition, since the foreground image obtained by segmentation is usually composed of discrete points, in order to detect the target such as trichomoniasis, we connect and fill the points with similar distances; the specific methods are hole filling, closing operation, and opening operation as shown in Figure 9.

(a)

(b)

(c)

Figure 9: Foreground map discrete points filling and connection: (a) original image; (b) foreground image; (c) processing result map.

The final calculation results are as follows.

As can be seen from Table 1, the improved algorithm proposed in this paper has improved by nearly 17% on the original basis. In addition, the effect of the frame difference method and the three-frame difference method is also good because the video shooting conditions are stable, and the test found that when the microscope focal length changes slightly, the effect is very poor, as shown in Figure 10.

Table 1: PCC comparison of different moving target foreground extraction models.

Method
PCC

VIBE [6]
0.7119

Frame difference [1]
0.7557

Three-frame difference [2]
0.7101

GMM [4]
0.3587

KDE [5]
0.0059

CodeBook [10]
0.1336

Kalman [11]
0.3958

Improved Kalman [9]
0.7169

PBAS [8]
0.3087

This paper
0.8803

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 10: Segmentation foreground renderings: (a) original image; (b) standard foreground image (ground truth); (c) algorithm in this paper; (d) VIBE; (e) frame difference; (f) three-frame difference; (g) GMM; (h) KDE; (i) CodeBook; (j) Kalman; (k) improved Kalman; (l) PBAS.

The split effect of frame 157 is shown in Figure 10.

In addition, we presented the intersection of union precision (IoU precision) indicator, which was used to determine the detection effect in target detection:

where represents the intersection of predicted box and ground-truth box. And represents the union. is the function of Area. represents the number of bounding box whose is greater than the threshold. represents the number of bounding boxes that were mistakenly segmented into foreground targets.

The IoU threshold, which is set to 0.5, and the IoU precision are shown in Table 2.

Table 2: Comparison of IoU precision for different moving target foreground extraction models.

Method
IoU precision

VIBE
0.6618

Frame difference
0.8097

Three-frame difference
0.8194

GMM
0.5105

KDE
0.0017

CodeBook
0.2081

Kalman
0.4536

Improved Kalman
0.8554

PBAS
0.3336

This paper
0.9485

As can be seen from Table 2, the proposed algorithm has an IoU precision of 94.8%, which is much higher than other foreground extraction algorithms. This also indicates that the algorithm has a good effect on trichomonas detection in leucorrhea.

The algorithm uses the Intel Core (TM) i7-4810MQ 2.80 GHz PC for detection. The average detection time is 98.08 ms, which meets the requirements of real-time video streaming trichomoniasis detection.

4. Conclusion

In this paper, a detection algorithm of trichomoniasis in leucorrhea microscopy image based on VIBE algorithm is proposed. The experimental results show that the improved algorithm can effectively extract the trichomoniasis region, and the PCC reaches 88%, which is higher than other video motions algorithm.

The proposed improved method can effectively suppress the false detection caused by the formed components such as epithelial cells in the leucorrhea microscopic image and the missed detection caused by the background model update during the movement. At the same time, improvements can effectively suppress smear and ghost areas. The algorithm has good robustness and strong adaptability, and the recognition effect is good under the weak flow of impurities and weak zoom in the video.

Data Availability

The video of trichomonas data used to support the findings of this study is available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

We would like to thank Professor Yutang Ye and all the staff in the MOEMIL Laboratory, who collected, labeled, and counted leukocytes in this study. This research was funded by the National Science Foundation of China (NSFC) (61405028) and the University of Electronic Science and Technology of China (UESTC) (ZYGX2015J046).

References

R. Jain and H. H. Nagel, “On the analysis of accumulative difference pictures from image sequences of real world scenes,” IEEE Transactions on Pattern Analysis & Machine Intelligence, vol. 1, no. 2, pp. 206–214, 1979.

M.-M. Yang and Q.-P. Guo, “A motion detection algorithm based on three frame difference and background difference,” in Proceedings of 8th International Symposium on Distributed Computing and Applications to Business Engineering and Science, pp. 487–490, Wuhan, China, October 2009.

Y. H. Cheng and J. Wang, “A motion image detection method based on the inter-frame difference method,” Applied Mechanics and Materials, vol. 490–491, pp. 1283–1286, 2014.

C. Stauffer and W. E. L. Grimson, “Adaptive background mixture models for real-time tracking,” in Proceedings of 1999 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, vol. 2, p. 257, Fort Collins, CO, USA, June 1999.

A. Elgammal, D. Harwood, and L. Davis, “Non-parametric model for background subtraction,” in Proceedings of International Symposium on Computational Intelligence and Design, pp. 363–367, Hangzhou, China, December 2007.

O. Barnich and M. Vanogenbroeck, “ViBe: a powerful random technique to estimate the background in video sequences,” in Proceedings of IEEE International Conference on Acoustics, Speech & Signal Processing, pp. 945–948, Taipei, Taiwan, April 2009.

M. V. Droogenbroeck and O. Paquot, “Background subtraction: experiments and improvements for ViBe,” in Proceedings of Computer Vision and Pattern Recognition Workshops, pp. 32–37, Providence, RI, USA, June 2012.

M. Hofmann, P. Tiefenbacher, and G. Rigoll, “Background segmentation with feedback: the pixel-based adaptive segmenter,” in Proceedings of Computer Vision and Pattern Recognition Workshops (CVPRW), pp. 38–43, Providence, RI, USA, June 2012.

R. Hao, X. Wang, J Zhang et al., “Automatic detection of trichomonads based on an improved Kalman background reconstruction algorithm,” Journal of the Optical Society of America A, vol. 34, no. 5, pp. 752–759, 2017.

K. Kim, H. Thanarat Chalidabhongse, D. Harwood, and L. Davis, “Real-time foreground–background segmentation using codebook model,” Real-Time Imaging, vol. 11, no. 3, pp. 172–185, 2005.

D. Koller, J. Weber, and J. Malik, Robust Multiple Car Tracking with Occlusion Reasoning, University of California at Berkeley, Berkeley, CA, USA, 1993.