# Location and Recognition Fruit Trees Based on Binocular Stereo Vision

Xueguan Zhao, Yuanyuan Gao, Songlin Wang, Xiu Wang, Pengfei Fan, Qingcun Feng

► To cite this version:

Xueguan Zhao, Yuanyuan Gao, Songlin Wang, Xiu Wang, Pengfei Fan, et al. Location and Recognition Fruit Trees Based on Binocular Stereo Vision. 11th International Conference on Computer and Computing Technologies in Agriculture (CCTA), Aug 2017, Jilin, China. pp.460-472, 10.1007/978-3-030-06179-1_46. hal-02111527

HAL Id: hal-02111527

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

Submitted on 26 Apr 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
Location and Recognition Fruit Trees Based on Binocular Stereo

Vision

Xueguan Zhao ¹,², Yuanyuan Gao ¹,², Songlin Wang ¹,², Xiu Wang ¹,² (✉),
pengfei Fan ¹,², Qingcun Feng ¹,²

¹Beijing Research Center of Intelligent Equipment for Agriculture, Beijing, China

²National Research Center of Intelligent Equipment for Agriculture, Beijing, China
{zhaoxg, gaoyy, wangsl, wangx, fanpf, fengqc}@nercita.org.cn

**Abstract.** In order to improve pesticide utilization rate and reduce the environmental pollution caused by pesticide ground loss, this paper proposes to use binocular vision to recognize the contour and distance information of fruit trees. To improve the recognition accuracy and speed, focusing on the optimization of SIFT stereo matching algorithm, A method for matching the feature points of left and right images based on cosine distance and the vector modulus is proposed. On this basis, two stereo matching algorithms are compared. The accuracy of the Improved SIFT stereo matching algorithm is improved by 1.53%. With this method, the recognition time is almost unchanged, and the stability of depth measurement is analyzed. When the target distance sensor is 180cm-220cm, the standard deviation is 1.3592cm, can meet the requirements of the work.

**Keywords:** Image recognition, binocular vision, stereo match, visual location

# 1 Introduction

The misuse of pesticide application is very common in the actual process of agricultural production. Excessive application of will not only cause the waste of pesticide, but also bring pollution to the environment. Automatic target spray technology of fruit trees is an important way to achieve high efficiency and low pollution. The realization of target spray based on machine vision technology has become a research hotspot in the field of precision spray technology at home and abroad due to its high flexibility and low equipment development cost.

Target technology mainly uses sensors for target detection, such as infrared sensors, ultrasonic sensors, laser radar, image sensors (CCD) and so on. Li Li et al. designed a target spraying control system, in which the infrared sensor was used to judge whether the target was available or not, and a green sensor was used to judge whether the target was green, so as to achieve the goal of applying pesticide only to green crops, thereby further reducing the waste of pesticides. However, because of the small detection area of infrared sensor, some switch jitter may be caused by the space between the branches and leaves. As the ultrasonic sensor could detect surfaces, it could solve the similar problem effectively. Gil et al. designed a multi nozzle sprayer consists of 3 ultrasonic sensors and 3 solenoid valves in early 2007. It could realize variable spraying according to the variation of grape leaves and save 58% pesticide compared with traditional spraying method.

Changyuan Zhai et al. built a target contour detection platform based on the ultrasonic sensor, and some target detection experiments were carried out on the regular crown and cherry trees. Good results

Beijing research center for intelligent equipment technology in agricultural; Beijing Nongke Mansion, No. 11, Shuguang huayuan Middle Road, Haidian District, Beijing, China; wangxiu@nercita.org.cn; zhaoxg@nercita.org.cn;
were obtained, which proved the feasibility of the ultrasonic target detection. Solanelles et al. applied ultrasonic sensors and proportional control valves to the air blower and controlled the spray flow according to the measured tree width, thereby greatly saving the pesticide. Although the cost of CCD image sensor is higher and the processing speed is slower, it can also detect crop diseases and insect pests with a certain image processing technology. Honghui Rao et al. used CCD image sensors to collect target information, and sprayed the target by controlling motor movement after images processing from CCD sensors. In 2010, Tianxiang Hu et al. studied the application of binocular vision technology in intelligent target spray, and further improved the accuracy and efficiency of spray.

Zacharie, Doerr used two-dimensional laser radar and GPS as sensors, and successfully developed a tractor autonomous navigation system for orchard operations. The accuracy of the system was higher

We can see that the application and research of machine vision in agriculture and forestry have entered the stage of development, and various theoretical and practical results emerge in an endless stream. But the accuracy of ultrasonic sensors and infrared sensors is poor, and laser ranging is expensive, at the same time target detection of machine vision is mainly used in weed identification in the field, and the research of target detection in orchard based on machine vision is still very few. At the same time, as the nerve center of target spraying system, the vision detection decision system needs to be further improved at the speed and accuracy of the detection algorithm.

## 2 Real-time precision target system

The hardware of Precision target spray decision system mainly include: Daheng CCD camera, PC (host computer), indoor fruit tree test stand. Image analysis development based on Microsoft Visual Studio 2012, It mainly to achieve image acquisition, segmentation, measurement, ranging, intelligent decision-making and sending results and other functions. The image acquisition function is completed by the CCD camera. The image of the fruit tree collected by the CCD camera is stored in the form of an image sequence to the computer for subsequent image processing.

The image acquisition function is completed by CCD camera, and the fruit images collected by the CCD camera are stored in the form of image sequences for the computer to be used for subsequent image processing. The applicable image segmentation, recognition and measurement algorithm is sought and develop software used for image segmentation and processing, realize description of fruit tree by gravity, perimeter, shape complexity, depth and other parameters. And can recover the real object of fruit tree from the image information to lay the foundation for the precise application of pesticides, functional flow chart shown in Figure 1.
<figure><img src="image_2.png" /><figcaption>Fig1.Functional flow chart</figcaption></figure>

## 3 System calibration

MATLAB camera calibration toolbox has higher calibration accuracy, and the calibration method is simple, therefore this paper uses MATLAB to about camera calibration, the mercury series Daheng camera MER-500-7UC, 2592*1944 resolution, spatial position and binocular vision sensor calibration plate as shown in figure 2:

<figure><img src="image_3.png" /><figcaption>Fig. 2 spatial position of calibration board and binocular vision sensor</figcaption></figure>

MATLAB's camera calibration toolbox has a high calibration accuracy, and calibration method is simple, so use MATLAB to calibrate the left and right cameras in this paper, and the calibration results are as follows, the internal parameters of the left and right camera are:

$$A = \begin{bmatrix} 3503.12 & 0 & 1207.31 \\ 0 & 3504.61 & 967.84 \\ 0 & 0 & 1 \end{bmatrix} \quad B = \begin{bmatrix} 3511.25 & 0 & 1198.93 \\ 0 & 3511.71 & 960.13 \\ 0 & 0 & 1 \end{bmatrix}$$
The equivalent focal length of the camera in the X direction is fx, the inner parameter matrix of the left camera obtained from the test shows that the fx of the left camera is 3503.12, This article adopts Computar brand 8mm focal length lens, the physical dimensions of a camera pixel in the X direction have been given by the camera: dx = 2.2um. According to fx = f / dx, can get f1 = 7.706mm, the same can get f r = 7.724mm, A computer series of 8mm focal length lenses are used. Which can be obtained left and right camera calibration error were 0.37% and 0.34% respectively. Then the rotation matrix and the translation matrix are calculated according to the common feature points, and R and T are

$$
R = \begin{bmatrix}
0.0220 & 0.6674 & 0.7446 \\
0.9964 & -0.0759 & 0.0386 \\
0.0822 & 0.7408 & -0.0666
\end{bmatrix}
$$

$$
T = \begin{bmatrix} -166.1646 & -84.2903 & 955.0699 \end{bmatrix}
$$

The obtained R and T results represent the matrix required for the left camera and the right camera to achieve coplanar, write R and T as XML files for cvstereoRectify function calls in Opencv for stereo correction, the results of the stereo correction are shown in fig3:

<figure><img src="image_4.png" /><figcaption>Fig3. image before and after the rectification</figcaption></figure>

4 Location

According to the principle of stereo vision, the simplest binocular stereo vision model is a stereo
camera made up of two parallel lenses to capture the same scene image, as shown in figure 4. Since there
is a distance between the two lenses (Figure 4, B), the target perceived through these two lenses will
produce bias in the captured image (Figure 4, dl ,dr). According to the triangulation principle, these
deviations are proportional to the distance between the camera and the target (Z in Fig. 4), so that these
deviations can be used to calculate the depth information of the target.
<figure><img src="image_5.png" /><figcaption>Fig4. Parallel optical axis geometric model of fruit trees Stereoscopic vision</figcaption></figure>

b represents the horizontal distance between the two cameras, the baseline of the stereo vision system, the F is focal length, and the Z is depth. The parallax d can be calculated according to the formula (2-1) once the deviation (dl and dr) of the horizontal direction of the two images is determined.

$$d = d' l - d' r \quad (1)$$

See from the coordinates in Figure 4, the dr is negative, so in fact the parallax d is the sum of dl and dr. Where the relevant parameters from the triangular similarity relationship can be obtained from the baseline mathematical expression:

$$b = (dl + dr) \cdot \frac{Z}{f} = d \cdot \frac{Z}{f} \quad (2)$$

Considering the relationship between length and pixels, the unit unity of parallax and depth calculation results is:

$$Z_{[mm]} = \frac{b_{[mm]} \cdot f_{[mm]}}{k_{[mm/pixel]} \cdot d_{[pixels]}} \quad (3)$$

In equation (3), Z, b, f units are mm, d is represented by pixels, and k is the size of each pixel.

## 5 Binocular matching

Identification of fruit trees in the environment in real time and accurately, it is the key to the target spraying vision system, and the essence of object recognition is image segmentation. In this study, indoor simulation experiments were conducted to collect images of fruit trees. Because there is a big difference between the surface color and the background color of fruit trees, there are different distribution characteristics in color space, so that image segmentation can be used to extract the target fruits and vegetables from the background. In this paper, 2G-R-B is used as a partition factor of the super green method. The picture of the fruit tree captured by the binocular vision system is shown in figure 5.

<figure><img src="image_6.png" /><figcaption>Fig5. A Pair of original stereo images</figcaption></figure>

### 5.1 SIFT feature point matching algorithm
The appearance of fruit trees is complicated and prone to mismatching. Due to a hypothetical of region matching is the space plane is the plane, which is parallel to the camera plane, in the actual orchard scene there are a lot of non-positive planes, so people began to consider the use of some of the more significant feature points (points of interest) to match, this method can also be called feature matching method.

SIFT feature matching algorithm is an image local feature description algorithm based on scale space, invariance to image scaling, rotation, and even affine transformations, therefore, it is widely used in matching technology. Because the number of feature points detected by SIFT feature point detection algorithm is larger, and each feature point descriptor is a vector of 128 dimensions. The feature utilization rate of the algorithm is not very high, it takes longer time, and there are some problems such as matching error or repetition.

## 5.2 Improved stereo matching strategy

As the most common distance measurement standard, “euclidean distance “is also used as a matching criterion for SIFT based matching, and good results are obtained. However, there are still problems such as matching errors or repetitions. From the analysis of the eigenvector itself: if the direction of the two vectors is the same, the smaller the angle is, on the basis of that, the modulus of the vector is taken into account. If the modulus is equal or close, the two vectors are considered equal. To determine whether the two vectors are in the same direction, it is necessary to use the cosine theorem to compute the angle of the vector. The cosine similarity of vector a and B is calculated as follows

$$ \cos \theta = \frac{a \cdot b}{|a||b|} \qquad (4) $$

From the above SIFT feature descriptor, we can see that the dimension of vectors a and b is a=(a1,a2, ...an)T、b=(b1,b2, ...bn)T, then the cosine of vectors a and b is:

$$ \cos \theta = \frac{\sum_{i=1}^{n} a_i y_i}{\sqrt{\sum_{i=1}^{n} a_i^2 \sum_{i=1}^{n} b_i^2}} \qquad (5) $$

Compared to Euclidean distance, cosine distance pays more attention to the difference of the two vectors in the direction. If we compare the length of the two vectors on the basis of the direction, we can accurately extract the feature matching points, the length of the vector, that is, the norm of the vector:

$$ \| x \| = \sqrt{x_1^2 + x_2^2 + L x_n^2} \qquad (6) $$

The specific implementation steps of the proposed stereo matching strategy are as follows: Firstly, the key feature points of the left / right image are extracted, and the image with the key feature points is used as the reference image, and the image with less key feature points is used as the image to be matched. The vector matrix of the K1 critical feature points of the baseline image is X, the vector matrix of the K2 critical feature points of the matched image is Y, X and Y are shown below:

$$ X = \begin{bmatrix} x_{1,1} & x_{1,2} & x_{1,3} & x_{1,4} \\ x_{2,1} & x_{2,2} & x_{2,3} & x_{2,4} \\ \mathbf{M} & \mathbf{M} & \mathbf{M} & \mathbf{M} \\ x_{k1,1} & x_{k1,2} & x_{k1,3} & x_{k1,4} \end{bmatrix} \quad Y = \begin{bmatrix} y_{1,1} & y_{1,2} & y_{1,3} & y_{1,4} \\ y_{2,1} & y_{2,2} & y_{2,3} & y_{2,4} \\ \mathbf{M} & \mathbf{M} & \mathbf{M} & \mathbf{M} \\ y_{k2,1} & y_{k2,2} & y_{k2,3} & y_{k2,4} \end{bmatrix} \qquad (7) $$
Take the first row vector of A $[x_{1,1} \ x_{1,2} \ x_{1,3} \ x_{1,1}]$, that is, the first critical feature point vector of the reference image, calculates cosine distance from all row vectors in B and takes inverse cosine, obtain an angle sequence of two vectors {$\theta_{1,1}, \theta_{1,2}, \theta_{1,3}, \mathbf{L} \theta_{1,k}$}. To increase robustness,

Take the minimum value of the sequence of five values corresponding to the vector for norm comparison, take X first row feature vector and Y line k feature vector as an example:

$$ \delta = abs(\| a \| - \| b \|) = \sqrt{x_{1,1}^2 + x_{1,2}^2 + \mathbf{L} x_{1,n}^2} - \sqrt{y_{k,1}^2 + y_{k,2}^2 + \mathbf{L} y_{k,n}^2} \quad (8) $$

Take the vector corresponding to the minimum of the δas the feature point. And so on, until all of the row vectors in A have been completed in turn to calculate the cosine distance of all row vectors in B, and the comparison of the norm. That is, all the key K1 feature points of the reference image are matched with all the K2 key feature points of the image to be matched.

### 5.3 key feature point matching contrast experiment

According to the research scheme described in this paper, the application of improved SIFT algorithm in stereo matching is studied, and simulation experiments are carried out to verify and analyze the feasibility of the improved stereo matching strategy. Computer configuration in the test are: CPU core i5-2410M, Memory is 2 GB, the operating system is Win7, and simulation platform is Matlab 2013. According to this method, the key feature points of left and right images are generated based on the improved SIFT algorithm, and the feature vectors extracted from the key feature points of left and right pictures are shown in Figure 6 (a) and (B),

<figure><img src="image_7.png" /><figcaption>Figure 6 Feature extraction result of the key feature points of the left / right image</figcaption></figure>

In this paper, the indoor environment image is selected, and the improved SIFT algorithm is tested under the illumination changes, scale changes and rotation changes. The matching results are shown in figure 7.

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

Control group

light is reduced by 30%.
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>Fig7. Matching results of improved SIFT algorithm</figcaption></figure>

In order to compare the difference in matching time between the improved SIFT algorithm and the SIFT algorithm, about 20 of the image acquisition in contrast, rotation, zoom, change of illumination conditions, using SIFT algorithm and improved SIFT matching algorithm for feature matching results statistics, statistical results from 20 on average of the image, as shown in the table1. The improved rule matching, stereo matching method increases the matching number, improve the matching efficiency, matching repetition or error reduce 1.53%, it is more advantageous to 3D reconstruction and localization of robot vision system.

<table><caption>Table 1 Statistical results of feature point matching</caption><thead><tr><th rowspan="2">Matching method</th><th colspan="2">Number of key<br>feature points</th><th rowspan="2">Number of<br>feature points<br>matching</th><th rowspan="2">Number of error<br>matching points</th><th rowspan="2">Correct matching<br>rate /%</th><th rowspan="2">Matching time<br>/ms</th></tr><tr><th>Left</th><th>Right</th></tr></thead><tbody><tr><td>SIFT matching<br>algorithm</td><td>490</td><td>498</td><td>206</td><td>7</td><td>96.60</td><td>697.21</td></tr><tr><td>Improved SIFT<br>matching algorithm</td><td>490</td><td>498</td><td>214</td><td>4</td><td>98.13</td><td>702.63</td></tr></tbody></table>

## 5.4 Binocular stereo vision experiment

The distance between the fruit tree and the camera is between 1.5 and 2.5m, the length of the baseline of the two camera is 50mm, and the algorithm is used to match. The matching effect is as Fig 8:
<figure><img src="image_11.png" /><figcaption>Fig8. Matching result of fruit tree binocular image</figcaption></figure>

## 5.5 Test results and discussions

In order to further observe the influence of matching region selection on the region matching effect, the SIFT matching algorithm and the improved SIFT matching algorithm are tested respectively. The binocular vision sensor is used to reconstruct the fruit tree, and then the width, height and depth of the fruit tree are measured in space, and the comparison is made. In order to test the accuracy of the binocular vision sensor in space reconstruction at different depths from the target, the distance between binocular camera and fruit tree from 1.5 m to 2.5 m is measured every 0.05 m. The measurement results are shown in Fig 9.

According to the three-dimensional reconstruction method, combined with the camera's internal and external parameters and parallax images, the 3D point cloud is measured at different angles, to get the height, width of the fruit trees and the depth in the vertical direction of the image plane. Where AB is the connection between vertex and diameter center of stem, The CD is the corners of the vertical line at 1/8 of the AB line, and the outer contour, 1 is the center of the AB line, 2 is the 1/8 at one end of the AB line, 3,4, respectively, at both ends 1/6 of the CD, the results shown in Table 1.

<figure><img src="image_12.png" /><figcaption>Fig9. Test point distribution</figcaption></figure>

<table><thead><tr><th rowspan="2">The actual distance between the lens and the target 1 points(cm)</th><th colspan="6">measured value</th></tr><tr><th>height (cm)</th><th>width</th><th>Target point A</th><th>Target point B</th><th>Target point D</th><th>Target point E</th></tr></thead><tbody><tr><td>150</td><td>120.0</td><td>45.5</td><td>136.4</td><td>147.6</td><td>168.6</td><td>173.3</td></tr><tr><td>160</td><td>119.9</td><td>46.7</td><td>147.3</td><td>148.7</td><td>176.1</td><td>179.7</td></tr></tbody></table>
<table><tr><td>170</td><td>120.4</td><td>45.2</td><td>159.4</td><td>164.5</td><td>190.7</td><td>194.5</td></tr><tr><td>180</td><td>121.8</td><td>47.1</td><td>177.8</td><td>182.0</td><td>208.2</td><td>207.5</td></tr><tr><td>190</td><td>122.1</td><td>48.4</td><td>186.6</td><td>192.4</td><td>219.6</td><td>219.2</td></tr><tr><td>200</td><td>121.4</td><td>48.1</td><td>210.4</td><td>212.8</td><td>227.1</td><td>241.6</td></tr><tr><td>210</td><td>120.9</td><td>47.8</td><td>221.2</td><td>224.3</td><td>244.0</td><td>252.6</td></tr><tr><td>220</td><td>122.0</td><td>48.8</td><td>228.9</td><td>230.8</td><td>260.3</td><td>261.9</td></tr><tr><td>230</td><td>122.1</td><td>49.1</td><td>236.6</td><td>246.8</td><td>274.4</td><td>271.9</td></tr><tr><td>240</td><td>124.5</td><td>49.6</td><td>258.9</td><td>261.0</td><td>287.5</td><td>287.1</td></tr><tr><td>250</td><td>123.2</td><td>50.3</td><td>273.8</td><td>274.9</td><td>302.6</td><td>302.4</td></tr><tr><td>True value</td><td>121.5</td><td>48.0</td><td>-</td><td>-5.4</td><td>+20.3</td><td>+22.7</td></tr></table>

The true value of the target 2-4 points is the relative depth of the target 1, and the "+" represents the direction of great depth.

In order to analyze the stability of the depth measurement, the difference between the target A and the mean value of target D,E as the length, the difference between the target D and the target E is half of the width, from the three-dimensional reconstruction of the measurement results can be found, compared to its height and width, length and manual measurement results, the standard deviation was 1.3805, 1.6224, 3.6081. And length value corresponds to the depth value in space, and its beating is greater compared with the width and height, the, and the width and height than the beating larger. When the sensor is away from the target at 180cm-220cm, the standard deviation is 1.174, relatively stable, so in the actual positioning of the apple space, the Deviation range of distance from visual sensor to the target 1 point range should be measured within the 40cm.

<figure><img src="image_13.png" /><figcaption>Fig10. Real-time measurement results</figcaption></figure>

Through a lot of test analysis, the error is mainly due to the accuracy of the depth information acquisition and several factors related to: the distance between the trees and the camera; the uniformity of light on the regional matching is also greater, due to the uneven light resulting in the match is not accurate, easy to cause the depth of information error, coupled with the impact of random noise, image segmentation results will be biased, affecting the results of distance measurement.

## Acknowledgement

This research was financially supported by the National Key R&D Program of China (2016YFD0200604).

## References

1. Li Li, Li Heng, He Xiongkui, et al. Development and experiment of automatic detection device for infrared target[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012(12): 159-163.

2. Zhai Changyuan, Zhao Chunjiang, Wang Xiu, et al. Probing method of tree spray target profile[J]. Transactions of
the CSAE, 2010(12): 173-177.

3. Gile, Escoaa, Roselljr, et al. Variable rate application of plant protection products in vineyard using ultrasonic sensors [J]. Crop Protection, 2007, 26(8): 1287 - 1297.

4. Stajnkod,Berkp, LESNIK M, et al. Programmable ultrasonic sensing system for targeted spraying in orchards[J]. Sensors, 2012, 12(11): 15500 - 15519

5. Solanelles F, Escola A, Planas S, et al. An electronic control system for pesticide application proportional to the canopy width of tree crops[J].Biosystems Engineering, 2006, 95(4): 473-481.

6. Rao honghui. Study on crop spraying control based on machine vision[J].Journal of Nanjing Agricultural University, 2007 (1): 120-123

7. Hu Tianxiang, Zheng Jiaqiang, Zhou Hongping. Distance measurement of tree images based on binocular vision. Proceedings of the agricultural mechanical society, 2010, 41 (11): 158-162

8. Cai Jianrong, Sun Haibo, Li Yongping. Three dimensional information acquisition and reconstruction based on binocular stereo vision. Chinese Journal of agricultural machinery, 2012, 43 (3): 152-156

9. Zacharie Doerr, Claude Lague. 2009 Evaluating the ability to detect foreign objects in crops using range scanners mounted on agricultural vehicles[C]. ASABE Annual International Meeting, Nevada, June 21-24.

10. ZOU L, LI Y. A method of stereo vision matching based on Open CV [C] //Proceedings of 2010 International Conference on Audio Language and Image Processing (ICALIP). Piscataway, NJ: IEEE, 2010: 185 – 190.

11. SPAMPINATO G, LIDHOLM J, AHLBERG C, et al. An embedded tereo vision module for industrial vehicles automation [C] //2013 EEE International Conference on Industrial Technology (ICIT). Pis- ataway, NJ: IEEE, 2013: 52-5

12. GANGANATH N, LEUNG H. Mobile robot localization using dometry and Kinect sensor[C] //2010 IEEE International onference on Emerging Signal Processing Applications (ES-PA). NJ, USA, 2010: 91-94.

13. Bradley D,Heidrich W. Binocular camera calibration using ectification error [C] //Computer and Robot Vision (CRV), 2010 Canadian Conference on. IEEE, 2010: 183-190.

14. Zhi Shuai. Research on the Error Theory of Binocular Vision easurement system on Non-cooperative Targe[D]. Changchun: Changchun Institute of Optics,fine Mechanics and Physics, Chinese Academy of Sciences, 2014.

15. Ravari A R N,Taghirad H D,Tamjidi A H.Vision-based Fuzzy Navigation of Mobile Robots in Grassland Environments[C]//Proc. Interational Conference on Advanced Intelligent Mechatronics.[S. I.]: IEEE/ASME Press,2009

16. SONG Tao, TANG Baoping, ZHAO Minghang, et al. Anaccurate 3-D fire location method based on sub-pixel edge detection and non-parametric stereo matching[J]. Measure-ment, 2014, 50: 160-171.