# Non-invasive Edge Detection of Leaves Based on Order Morphology

Yanlei Xu, Qi Zhang, Chenxiao Li, Xindong Wang, Xiaotian Meng

► To cite this version:

Yanlei Xu, Qi Zhang, Chenxiao Li, Xindong Wang, Xiaotian Meng. Non-invasive Edge Detection of Leaves Based on Order Morphology. 11th International Conference on Computer and Computing Technologies in Agriculture (CCTA), Aug 2017, Jilin, China. pp.428-439, 10.1007/978-3-030-06137-1_39. hal-02124217

HAL Id: hal-02124217

https://inria.hal.science/hal-02124217

Submitted on 9 May 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Non-Invasive Edge Detection of Leaves based on Order Morphology

Yanlei Xu¹(✉), Qi Zhang¹, Chenxiao Li¹, Xindong Wang¹, and Xiaotian Meng¹

1. College of Information and Technology, Jilin Agricultural University, Changchun, China
(yanleixu,z163zhangqi,lichenxiao2000,dong__2017,mengxiaotian11}@163.com

**Abstract.** Non-invasive edge detection of leaves is the key step of leaf feature extraction. Traditional algorithms for edge detection of leaves are usually invasive, since the detection is carried out after the leaves are picked. We apply order morphology in this study to non-invasive edge detection of leaves. First we analyze the algorithm for order morphology edge detection of leaves. In particular, the impact of structural elements and percentile on the detection is discussed. Based on the theory of order morphology transform of leaf images, the operator for detecting the leaf edge is constructed. Finally simulation experiment is carried out on leaf images under the conditions of natural illumination and artificial noise, respectively. Results show that the proposed algorithm is accurate and fast in leaf edge extraction and not sensitive to noise.

**Keywords:** order morphology; plant leaf; edge detection.

## 1 Introduction

In agricultural field, non-contact detection and measurement of plants as a way of monitoring plant growth status is an emerging topic. Leaves are the organs where photosynthesis takes place, and the morphological features of leaves represent an important aspect of plant growth status[1]. Shape is one major morphological feature of leaves[2-7], and edge extraction of leaves is one step towards leaf shape extraction.

Traditional edge detection algorithms are invasive and involve the use of picked leaves. In non-invasive detection and measurement, leaf images may present the leaves as overlapped with complex shape features. It is difficult to determine the curvature at different positions along the edge. Moreover, illumination also has a great impact on the clarity of the extracted edge. It is now urgent to establish an algorithm for non-invasive edge detection of leaves in open field.

Edge detection algorithms generally depend on the discontinuity of pixel features in different regions. Edges exist between adjacent regions that have different pixel features, and the grayscale values of the pixels along the edge are discontinuous. Considering this fact, edge detection can be carried out using differential operators. Good segmentation can be achieved with Roberts operator, Sobel operator, Prewitt operator and Canny operator where the edges between different regions are clear and the grayscale values of the pixels vary more obviously. But these operators are sensitive to noise[8-9]. Leaf images collected under natural light usually contain considerable noise and the edge information extracted is obscure. Conditional differential operators are not applicable in that case.

More studies are concerned with edge detection of picked leaves than with non-invasive edge detection of leaves. Z. Husin et al. [10] proposed threshold-based segmentation for edge detection of leaves, but the algorithm was applied to the entire leaf rather than to the edge. Therefore, this method is not conducive to

* Yanlei Xu, College of Information and Technology, Jilin Agricultural University, Xincheng Avenue2888, Changchun,130118, CHINA., E-mail:yanleixu@163.com
shape extraction, and only specific types of crops can be detected. More importantly, the leaves were detected only after being picked. Li et al. [11]proposed the algorithm for edge detection of leaves in cotton images based on Mean-shift and lifted wavelet transform. This algorithm reduced the non-edge noise and more effectively extracted the edge of overlapping leaves. But this algorithm was more fit for images collected under laboratory conditions than under natural light. Lin et al. [12] described the method of edge detection of leaves based on fuzzy logics, which reduced the pseudo-edge of leaves. But this method was confined to leaf images collected in greenhouse environment. Xu et al. [13] put forward the method of edge detection of leaf veins based on fuzzy order morphology. Hu et al. [14]established the method of edge detection of seedlings of broad-leaved trees based on intuitionistic fuzzy set. The above algorithms are based on fuzzy theory and they are capable of extracting clear edges, though the computing load is very large and the complexity is high.

Order morphology theory is the combination of mathematical morphology and order statistics. Order morphology transform of digital images is the ordering of finite data of grayscale values of pixels within the range of structural elements[15]. Order morphology provides good filtering effect and noise inhibiting effect, especially in the open field. We first analyze the order morphology theory and the order morphology transform of leaf images. Structural elements and percentile are identified as the influence factors of the detection result. Finally, the leaf edge is extracted by order morphology transform. Simulation experiments are carried out under natural light condition, and the extracted leaf edges are clear and accurate and insensitive to noise.

# 2 Theory of order morphology transform

## 2.1 Basic concepts

The concept and property of order morphology transform are first defined[16-17].

Definition 1: Suppose the set of discrete points is A, B is structural element, $0 \le \mu(B) = k < +\infty$ (metric $\mu(\cdot)$ is the count of B points). Thus the order morphology transform of A with respect to B is: $A \oplus B$ ($p=0,1/(k-1)$, $2/(k-1), \dots, 1$), defined as

$$x = A \oplus B = \{x : \mu(A \cap B) = (k \ 1)p+1\}, p = 0, \frac{1}{k-1}, \frac{2}{k-1}, \dots, 1 \quad (1)$$

where p is the percentile for order morphology transform. Discrete order morphology transform can be based either on morphology sum and morphology difference.

$$A \oplus B = \{x : \mu(A \cap B) = 1\} = A \odot B \quad (2)$$

$$A \oplus B = \{x : \mu(A \cap B) = k\} = A \oplus B \quad (3)$$

The geometrical meaning of discrete is: $A \oplus B$ is the x set of $(k-1)p+1$ points which B include A

The properties of order morphological transformation $A \oplus B$ are: monotonicity, expansibility and inverse expansibility.

Property 1 $A \oplus B \subset A \oplus B \subset A \oplus B$
$A \textcircled{1} B \subset A \textcircled{2} B$, $0 \le p_1 \le p_2 \le 1$

Property 2 $A \subsetneq C \Rightarrow A \oplus B \subsetneq C \oplus B$

Property 3 $(A \oplus B)^c = A^c \oplus B, p+q=1$

Property 3 show that the morphology transformation $A \oplus B$ and $A \oplus B$, $p+q=1$, are the dual transformation.

Definition 2: Suppose $0 \le p, q \le 1$, $A(p,q)B = A \oplus B \oplus$ is called as multiplex order morphology transformation, obviously:

$$A(0,1)B = A \circ B \quad (4)$$

$$A(1,0)B = A \bullet B \quad (5)$$

$$A(0,1)B \subset A \subset A(1,0)B \quad (6)$$
$$A^c(p,q)B = [A(1-p,1-q)B]^c \qquad (7)$$

Definition 2: n multiplex order morphology is defined as:

$$A^{\oplus n}B = (\dots(A^{\oplus p}B) \oplus B)\dots B \dots B \qquad (8)$$

## 2.2 Order morphology transformation of plan leaf

The order morphology is applied to the plant leaf, which is the order transformation of plan leaf.

(1) Z(n) is the n-D digital space

(2) f or g is the gray image of plant leaf

(3) B is the structure element

(4) $nB=B \oplus B \oplus B \oplus B \dots \oplus B$

(5) $f \le g$ is $f(x) \le g(x)$, $f-g$ is $f(x)-g(x)$, $x \in Z(n)$

Definition 1 Supposed plant leaf image f: $0 \le f(x) \le m(m$ is the gray value), $B=\{x_1,x_2, \dots,x_k\}$ is the structure element, $0 \le \mu(B)=k<+\infty$ $\mu(\cdot)$ is the count of B. The order of k values in structure element B is: $f(x_1^*) \le f(x_2^*) \le \dots \le f(x_k^*)$, define the d rank order variable in structure element B:

$$\text{order}(d) = \frac{d!}{d!(k-d)!} B \quad (d \le k, 2, \dots, k) \qquad (9)$$

Definition 2 The d rank order variable of f(x) in B ($d=(k-1)p+1$ is rank) is defined as the order morphological transformation of f about B, marked $f^{\oplus}B$, that is

$$g_{i,j}^{(d)} = f^{\oplus}B = \text{order}(d: f/B) = f(x_d^*)(p=0,1/(k-1),2/(k-1),\dots,1) \qquad (10)$$

Definition 3 $f(p,q)B = (f^{\oplus}B)\oplus B$, ($p,q=0,1/(k-1),\dots,1$) is named duplicate symmetrical mix order morphological transformation of f about structure element B.

# 3 Analysis of edge features in order morphology transform

As shown above, order morphology transform can act as a filter using different percentile p on the grayscale images. Thus the grayscale images can be treated by various mathematical operations [18-28].

## 3.1 Filtering and edge detection by order morphology transform

The leaf images contain much noise, which can be filtered by order morphology transform using different percentile. When $0 < p < 1/2$, the positive pulse noise in the image can be inhibited; when $1/2 < p < 1$, the negative pulse noise is inhibited[29]. Two transforms can be performed separately, that is, using p value of 0-1/2 and p of 1/2-1. Thus both positive and negative pulse noise is eliminated. Suppose $0 < p < 1/2$ and $1/2 < q < 1$. Let $f(p,q)$ B be open operation and $f(q,p)$ B closed operation, both of which can act as filter. The open operation can filter out the high-luminance noise and the closed operation the smaller noise. Therefore, all types of noise can be removed by setting different p and q values.

The noise in the image can be removed and the edge be detected using the above operators. Fig. 1 is the schematic diagram of edge detection based on order morphology transform. When structural element B is located in the flat region, the grayscale values of pixels within B do not vary much. So the difference is not great before and after the transform. When B is located in the region of abrupt change of grayscale values, the grayscale values of the pixels within B vary considerably. The output image will be very different from the original image. When $0<p<1/2$, the output image $f^{\oplus}B$ is subtracted from the original image f. Given the anti-extensibility of order morphology transform, the difference $f-f^{\oplus}B$ ($0<p<1/2$) is the edge information of the original image. When $1/2<p<1$, the original image f is subtracted from the output image $f^{\oplus}B$, and the difference $f^{\oplus}B-f(1/2<p<1)$ is the edge information of the original image.
<figure><img src="image_2.png" /><figcaption>Fig.1 Schematic diagram of edge detection</figcaption></figure>

The working principle of edge detection of leaves based on order morphology transform can be described briefly as follows: When the structural element B is located within the range of the leaf, the grayscale values of the pixels within B vary slightly. Therefore, the input and output grayscale values differ little. When B is located over the edge of the leaves, the grayscale values within B vary greatly, so the input and output grayscale values also show large difference. By amplifying the changes of grayscale values of the edges through the transform, the leaf edges can be accurately recognized. Thus the edge detection operator based on order morphology has the following form:

Order morphology has unparallel advantages in edge detection by choosing proper structural elements and percentile. We build the operator based on order morphology for edge detection of leaf images.

$$G(f) = f - f^{\oplus} B \qquad (p < 1/2)$$

$$G(f) = f^{\oplus} B - f \qquad (p > 1/2) \qquad (11)$$

## 3.2 Building the edge detection operator

If only one structural element is used, then the output image will contain only one type of geometrical information, which means a loss of leaf edge details. To detect the edge more accurately, several structural elements may used. These structural elements may differ in size and shape. We present the edge detection operator based on several structural elements.

$$G(f) = f^{\oplus} n_1 B_1 \oplus n_2 B_2 \oplus n_3 B_3 - f^{\oplus} n_4 B_4 \oplus n_5 B_5 \qquad (12)$$

$$G(f) = f^{\oplus} n_1 B_1 \oplus n_2 B_2 - f^{\oplus} n_3 B_3 \oplus n_4 B_4 \oplus n_5 B_5 \qquad (13)$$

$$G(f) = f^{\oplus} n_1 B_1 \oplus n_2 B_2 \oplus n_3 B_3 - f^{\oplus} n_4 B_4 \oplus n_5 B_5 \oplus n_6 B_6 \qquad (14)$$

In the operator, nᵢ represents different structural element; Bᵢ indicates the shape of the structural element. The outer edge of the image is extracted in formula (12), and the edge spanning over the Euclidean border is extracted in formula (14). Different shape of the structural elements corresponds to order morphology operation in different directions of the leaf. The final leaf edge is obtained by combining together the edges extracted in all directions. For a given structural element, different size corresponds to different scale of extraction. Small-size structural elements have lower noise-removing capacity, but they are more sensitive to edge details. On the contrary, large-size structural elements have stronger noise-removing capacity, but the more details may be lost. So the best way is to use several structural elements with varying size and shape. An illustration is carried out using formula (14). This formula extracts the Euclidean edge of the leaf. Suppose the structural elements n₁B₁ and n₂B₂, n₄B₄ and n₅B₅ have identical shape and size, respectively. If n₃B₃ and n₆B₆ have the same shape but different size, then the larger the n₃ and n₆, the wider and the brighter the extracted edge; otherwise, the edge will be narrower and darker. If n₃B₃ and n₆B₆ have the same size but
different shape, it is necessary to choose appropriate shape of the structural element depending on specific
leaf so as to extract the edge more accurately.

3.3 Impact of structural element and percentile on leaf edge detection

In mathematical morphology, the choice of structural elements has a large impact on the extracted edge[30]. This is also true in order morphology transform, but no relevant studies have been reported yet. We discuss the impact of structural elements and percentile on leaf edge detection.

The thickness and direction of the extracted edge can be controlled by controlling the shape and size of the structural elements. The structural element is expressed as the morphology sum of several structural element Bj, which is nᵢBj (n is the size and B is the shape). Fig. 2 shows the edges of different thickness when the size of the structural element is 3×3, 4×4 and 5×5, respectively, with the constant shape of a rectangle.

<figure><img src="image_3.png" /><figcaption>Fig.2 Leaf edge of different thickness</figcaption></figure>

It can be seen from Fig. 2 that the larger the size of the structural element, the thicker the extracted edge is. But thicker edges are usually obscure with the generation of some pseudo-edges, which increases the difficulty of further processing. To overcome this defect, the size of the structural element can be varied depending on the needs.

Similarly, by changing the shape of the structural elements, the edges extracted will show different features. Given the same size of the structural elements, different shape of the structural elements will result in different features of edges. Here we choose the constant size of 3×3 and vary the shape of the structural elements. The results are shown in Fig. 3.

<figure><img src="image_4.png" /><figcaption>Fig.3 Edge of different shape</figcaption></figure>

As shown in Fig. 3, edges with different features are extracted by setting different shape of the structural
elements. For example, a cross-like shape results in satisfactory edges in transverse and vertical direction,
but poor edges in other directions. In Fig. 3(b), both the cross-like shape and the diagonal shape result in
clear edges. Because the leaf edges are approximately elliptical or circular or oblique lines but rarely
rectangular, the use of rectangular structural element results in less clear edges (See Fig. 3(a)). It is important
to choose the appropriate shape of structural elements for specific leaves, as discussed below.
In order morphology transform, the values of percentile (p and q) are important to edge detection. We discuss the impact of percentile on edge detection by choosing different percentile values. Let $0 < p < 1/2$ and $1/2 < q < 1$ and the structural element is rectangular with size of $6 \times 6$. The leaf vein edges are extracted at different p and q values, as shown in Fig. 4. When $p=1/3$ and $q=35/36$, the extracted leaf edges are shown in Fig. 4(a); when $p=1/36$ and $q=35/36$, the extracted leaf edges are shown in Fig. 4(b); when $p=1/3$ and $q=35/36$, the leaf edges are shown in Fig. 4(c).

<figure><img src="image_5.png" /><figcaption>Fig. 4 Edge of different percent</figcaption></figure>

It is easy to see that if p and q differ greatly, the extracted edges are more clear and thicker, though some details are lost (comparing (a) and (b), the leaf veins are more clear). But pseudo-edges appear in (a). After several experiments, it is found that a large difference existing between p and q is more conducive to edge detection.

## 4 Simulation experiment on order morphology edge detection of leaves

We carry out simulation experiment on order morphology detection of leaves under natural light condition. Fig. 5(a) is the original leaf image, and the leaf edges are extracted using traditional Sobel operator and Canny operator. The results are shown in Fig. 5(b) and (c). We apply order morphology edge detection operator in formula (11) to extract leaf edges, as shown in Fig. 5(d). Using the proposed operator based on several structural elements of different size and shape, the extracted edges are shown in Fig. 5(e) and (f).

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>(d) classical order morphology</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(e) operator1 of the paper</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(f) operator2 of the paper</figcaption></figure>

<caption>**Fig.5** Edge detection of plant leaf</caption>

The leaf edges extracted by Sobel operator are less clear and discontinuous with considerable loss of details. Canny operator fails to extract leaf edges and the resulting image is obscure and does not show the general contour of the leaves. Though the traditional order morphology operators can detect the leaf edges, the edges have low contrast and do not present a clear contour with loss of some of the edge information. Formula (14) is used as the edge detection operator based on several structural elements of varying size and shape. The edges in Fig. (e) and Fig. (f) are extracted using different structural elements. In Formula (14), $n_1$ is 5×5 and B1 is rectangular; $n_2$ is 3×3 and B2 is rectangular; $n_3$ is 7×7 and B3 is cross-like shape. $n_4B_4$ and $n_5B_5$, $n_1B_1$ and $n_2B_2$ have identical shape and size, respectively. When $n_6$ is 5×5 and $B_6$ is rectangular, the detection result is shown in Fig. (e). Then $n_1$ is 5×5 and $B_1$ is rectangular; $n_2$ is 3×3 and $B_2$ is rectangular; $n_3$ is 5×5 and $B_3$ is diagonal; $n_4B_4$ and $n_5B_5$, $n_1B_1$ and $n_2B_2$ have identical shape and size, respectively. When $n_6$ is 5×5 and $B_6$ is cross-like shape, the detection result is shown in Fig. (f). All edges are clear and intact using different structural elements. The edges in Fig. (3) have the highest clarity and contrast, because the structural elements used are larger and have more variable sizes than those in Fig. f.

One major advantage of order morphology edge detection is that the noise is filtered simultaneously. So we perform another experiment where 10% salt and pepper noise is added into the leaf images. The edge detection of leaves is carried out using Sobel operator, Canny operator and the proposed operator in Formula (11), respectively. The results are shown in Fig. 6.

<figure><img src="image_10.png" /><figcaption>(a) leaf image (10% S&P noise)</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(b) Sobel operator</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(c) Canny operator</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>(d) classical order morphology</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>(e) operator1 of the paper</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>(f) operator2 of the paper</figcaption></figure>

<caption>**Fig.6** Edge detection of noised plant leaf</caption>
It can be seen from the figure that the traditional operators are sensitive to noises and they cannot filter the
noise; the edges extracted are obscure and contain a large amount of noise. The traditional order morphology
transform operator only extracts the general contour of the leaves contaminated by the noise. Our operator
obtains very clear leaf edges. The structural elements in Fig. (e) and Fig. (f) are the same as those in Fig. 5. It
can be seen that the operator 1 is little affected by noise and the extracted leaf edges are clear and intact, with
noise basically filtered. Operator 2 is also effective in removing a major part of noise. As shown in Fig. (f), the
extracted leaf edges are basically not affected by noise, and compared with Fig. 5, the edges are clear and intact.
This proves the ability of the proposed operator in filtering noise and its high applicability to natural light
condition. Since the proposed operator allows for the adjustment of size and shape of structural elements and
the percentile value, it can adapt to a variety of leaf features perfectly.

5 Conclusion

We propose a method of non-invasive order morphology edge detection of leaves. Starting with thorough
analysis of the order morphology theory and order morphology transform of leaf images, we proceed to the
noise filtering and edge detection ability of order morphology. Structural elements and percentile have large
impact on the detection results. We then build the order morphology edge detection operator based on several
structural elements with varying size and shape. The algorithm is verified by simulation experiment.

The experiment consists of two parts. One is the edge extraction from leaf images collected under natural
light condition. The other is edge extraction in the presence of artificial noise. Results show that the leaf edges
extracted by the proposed operator are more clear and accurate and insensitive to noise. This operator can
adapt to natural light condition and can be used for shape recognition of plant leaves and plant growth
monitoring in the future.

References

1. Allen J Coombes, Zsolt Debreczy. The Book of Leaves [M]. London, The University of Chicago Press, Ltd, 10-15(2010)

2. Gheng L Y, Zhang J T, Wang Q Y. Mean-shift-based color segmentation of images containing green vegetation[J]. Computer and Electronics in Agriculture, 65 (1):93-98(2009)

3. Lopez-Molina C, De Baets B, Bustince H. Generating fuzzy edge images from gradient magnitudes[J]. Computer Vision and Image Understanding, 115 (11):1571-1580(2011)

4. M.Wen and C.Zhong. Application of sobel algorithm in edge detectionof images[J].hina High-tech Enterprise, 6(5):57-62(2008)

5. Wang Dong, Zhou Shisheng. Color image recognition method based on the prewitt operator [C] //International Conference on Computer Science and Software Engineering, 170-173(2008)

6. Cheng Jinyong, Xue Ruojuan.Segmentation of medical images with canny operator and GVF snake model[C] // Proceedings of the 7th World Congress on Intelligent Control and Automation, 899-903(2008)

7. Tizhoosh H R. Fast fuzzy edge detection [C] // Proceedings of Annual Meeting of the North American on Fuzzy Information Processing Society, 239-242(2002)

8. Ma Xianmin, Jiang Yong. A fast edge detection roberts algorithm of coal gangue binary image[J]. Chinese Journal of Scientific Instrument, 26(4):595-597(2005)

9. M.Wen and C.Zhong. Application of sobel algorithm in edge detectionof images[J].hina High-tech Enterprise, 6(5):57-62(2008)

10. Z. Husin, A.Y.M. Shakaff, A.H.A. Aziz. Embedded portable device for herb leaves recognition using image processing techniques and neural network algorithm[J]. Computers and Electronics in Agriculture, 89(1): 18-29(2012)

11. Li Han, Wang Ku, Bian Haoyi. Cotton leaf image edge detection using Mean-shift algorithm and lifting wavelet transform[J]. Transactions of the CSAE, 26(Supp.1): 182-186(2010)

12. Lin Kaiyan, Si Huiping, Zhou Qiang, et al. Plant leaf edge detection based on fuzzy logic[J]. Transactions of the Chinese Society for Agricultural Machinery, 46 (6) : 227-231. (in Chinese with English abstract)( 2013)

13. Xu Yanlei, Jia Honglei, Bao Jialin. Plant leaf vein edge detection based on fuzzy order morphology [J]. Transactions of the Chinese Society of Agricultural Engineering, 31(13): 193-198(2015)

14. Hu Chunhua, Li Pingping. Edge detection of hardwood seedlings leaves based on intuitionistic fuzzy set [J]. Journal of Nanjin Forestry University(Natural Sciences Edition), 38(6):193-198(2014)

15. Arce G R, Foster R E. Detail preserving ranked-order based filters for image processing[J]. IEEE Trans. on ASSP,
ASSP-37(1):83-98(1989)

16. YAN Xue-qiang. Applications of Order Morphology to Image Edges Detection[J]. Signal and Processing. 13(4):357-362(1997)

17. Dong Hu, Xianzhong Tian. A Multi-directions Algorithm for Edge Detection Based on Fuzzy Mathematical Morphology[C]. Proceedings of the 16th International Conference on Artificial Reality and Telexistence Workshops (ICAT'06), 2133-2136(2006)

18. Jean Serra. Morphological filtering: An overview[J]. Signal Processing, 38(1): 3-11(1994)

19. Jean Cousty, Laurent Najman, Fabio Dias, Jean Serra. Morphological filtering on graphs[J]. Computer Vision and Image Understanding, 117(4): 370-385(2013)

20. Bart Wilburn. Theory of ranked-order filters with applications to feature extraction and interpretive transforms[J]. Advances in Imaging and Electron Physics, 112:233-332(2000)

21. Ye Bin, Peng Jiaxiong, Lu Hanqing. Study on Application of Order Morphology Filtering for Detection of Small Target and Point Target[J]. Acta Automatic Sinic, 28(6): 990-994(2002)

22. ZHAO Chun-hui, HUI Jun-ying, WANG Wei, etc. A Class of Adaptive Ranked-Order Morphological Filters[J]. Journal of Image and Graphics, 5(8):674-677(2000)

23. Yanlei Xu, Jiyin Zhao, Yubin Jiao. Gray-scale Image Edge Detection Based on Order Morphology Transformation[C]. The 7th World Congress on Intelligent Control and Automation (WCICA'08), 5970-5975(2008)

24. Chi Jian-nan, Fang Shuai, etc. Gray-scale Image Edge Detection Based on Multi-structuring Elements Order Morphology Transformation[J]. Journal of Image and Graphics, 11(1):41-46(2006)

25. Chi Jian-nan, Fang Shuai, etc. Infrared Image Edge Detection Based on Multiplex Order Morphology Transformation[J]. Journal of Northeastern University ( Natural Science), 26(2): 103-106(2005)

26. Yan Haixia, Liu Yanjun. Edge detection method based on adaptive order morphology filter[J]. Application Research of Computers, 28(5): 1979-1980(2011)

27. Yan Haixia, Zhao Xiaohui. Edge detection method based on Tsallis entropy difference of remotesensing image.[J]. Application Research of Computers, 26(9): 3598-3600(2009)

28. Ye Bin, Peng Jiaxiong. Mov ing Small Target Detect ion Based on Order Morphology Filtering in Infrared Image Sequences[J]. Journal of Data Acquisition & Processing, 16(3): 315-319(2001)

29. Jiyin Zhao, Yanlei Xu, Yubin Jiao. A Kind of Fast Arithmetic of Gray-scale Image Edge Detection Based on the Order Morphology[J]. Chinese Journal of Electronics, 36(11): 2195-2199(2008)

30. GUO Jingfeng, LIN Xudong. The Analysis and Study of Structure Element of Mathematical Morphology[J]. Computer Science, 29(7):113-115(2002)