Estimation of Crop Height Distribution for Mature Rice Based on a Moving Surface and 3D Point Cloud Elevation
Abstract
:1. Introduction
2. Materials and Methods
2.1. LiDAR Measuring System and Point Cloud Data Acquisition
2.1.1. LiDAR Measuring System Construction
2.1.2. Data Acquisition of Field-Mature Rice Point Cloud
2.2. Point Cloud Classification Based on Moving Surface Fitting Elevation
2.2.1. D Point Cloud Denoising
- (1)
- Determine neighborhood radius r and the standard deviation-multiple threshold m.
- (2)
- Traverse all point cloud data xi within the neighborhood radius r, and calculate the average distance di between all point cloud data; the calculation of d is as Equation (1):
- (3)
- Calculate the mean u and standard deviation σ of each point; the calculation equation is shown in Equations (2) and (3)
- (4)
- Define the point not belonging to the noise point if μ-m∙σ ≤ di ≤ μ-m∙σ, delete it otherwise;
- (5)
- Traverse all the point cloud data, save the undeleted point cloud data.
2.2.2. Ground and Crop Point Cloud Classification
- All point cloud data have 3D coordinates (Xi, Yi, Zi). In order to speed up the process of finding adjacent points and reduce the time spent in subsequent classification, all point cloud data are sorted based on the K–D tree algorithm.
- Among all point cloud data, the three points with the smallest height value are taken as the initial ground points. The x and y values of these three points are input into Equation (5), and the surface equation of the initial ground is obtained by fitting; the initial ground height is Zo.
- In order to reduce the amount of data processing, a threshold T1 is set, which is determined according to the fluctuation of the ground to remove the points that do not belong to the ground points, determined according to the following rules:
- For the remaining points that need to be further classified (that is, the target area), find the adjacent points of the initial ground point according to the K–D tree and substitute the x and y values of the adjacent points into the initial ground equation obtained in step 2 to obtain the fitting height Zi’. If the difference between the fitted value and the actual value is within the set threshold range, it is considered that the point belongs to the ground point cloud data. Otherwise, it does not belong to the ground point data.
- The newly determined ground points are fitted with the initial ground points. When the number of ground points is greater than six, quadratic surface fitting is performed according to Equation (4). Ground points and discriminate new data are added continuously if the adjacent new point cloud data are ground points. Once the number of fitted ground points reaches the set number, the number of ground point cloud data would be retained. As the new ground point is being determined, the old one would be removed, and the surface equation would be updated as well.
- Step 5 would be repeated until all point cloud data are determined; then the ground and non-ground point cloud data would be saved separately.
2.3. Plant Height Distribution Estimation
2.3.1. Mesh Generation Point Cloud Elevation
- (1)
- Calculate the boundaries of point cloud data (L1, L2) = (Xmax − Xmin, Ymax − Ymin); then, the coordinate XOY is established according to the data boundaries.
- (2)
- Set the grid side length l according to the requirements of plant height estimation resolution, and the number of grids a and b are determined based on Equations (8) and (9).
- (3)
- The ground and crop data point clouds are divided into each grid area ([Xmin + (m − 1)l, Xmin + ml], [Ymin + (n − 1)l, Ymin + nl]) of the grid coordinate system; among them, m ∈ [1,2,3……a], n ∈ [1,2,3……b].
2.3.2. Estimation of Rice Plant Height Distribution
3. Results
4. Discussion
5. Conclusions
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Liang, Z.W.; Li, Y.M.; Xu, L.Z. Grain sieve loss fuzzy control system in rice combine harvesters. Appl. Sci. 2018, 9, 114. [Google Scholar] [CrossRef] [Green Version]
- Maertens, K.; Baerdemaeker, J.D. Design of a virtual combine harvester. Math. Comput. Simulat. 2004, 65, 49–57. [Google Scholar] [CrossRef]
- Chai, X.Y.; Xu, L.Z.; Li, Y.; Qiu, J.; Li, Y.M.; Lv, L.Y.; Zhu, Y.H. Development and experimental analysis of a fuzzy grey control system on rapeseed cleaning loss. Elecrtronics 2020, 9, 1764. [Google Scholar] [CrossRef]
- Zhang, Y.W.; Du, C.; Yin, Y.X.; Wang, X.; Wang, S.M. Experimental study of feed rate related factors of combine harvester based on grey correlation. IFAC-PapersOnLine 2018, 51, 402–407. [Google Scholar] [CrossRef]
- Zhang, Z.Q.; Sun, Y.F.; Liu, R.J.; Zhang, M.; Li, H.; Li, M.Z. Design and test of feed rate monitoring system for combine harvester. Nongye Jixie Xuebao/Trans. Chin. Soc. Agric. Mach. 2019, 50, 85–92. [Google Scholar]
- Long, D.S.; McCallum, J.D. On-combine, multi-sensor data collection for post-harvest assessment of environmental stress in wheat. Precis. Agric. 2015, 16, 492–504. [Google Scholar] [CrossRef]
- Chang, A.J.; Jung, J.H.; Maeda, M.M.; Landivar, J. Crop height monitoring with digital imagery from Unmanned Aerial System (UAS). Comput. Electron. Agr. 2017, 141, 232–237. [Google Scholar] [CrossRef]
- Maertens, K.; Reyns, P.; DeClippel, J.; DeBaerdemaeker, J. First experiments on ultrasonic crop density measurement. J. Sound Vib. 2003, 266, 655–665. [Google Scholar] [CrossRef]
- Munch, P. Predictive feed-rate control for combine harvesters. LandTechnik AgEng 2019, 2361, 545–552. [Google Scholar] [CrossRef]
- Yuan, W.; Li, J.; Bhatta, M.; Shi, Y.; Baenziger, P.S.; Ge, Y. Wheat height estimation using LiDAR in comparison to ultrasonic sensor and UAS. Sensors 2018, 18, 3731. [Google Scholar] [CrossRef] [Green Version]
- Saeys, W.; Lenaerts, B.; Craessaerts, G.; Baerdemaeker, J. Estimation of the crop density of small grains using LiDAR sensors. Biosyst. Eng. 2009, 102, 22–30. [Google Scholar] [CrossRef]
- Blanquart, J.E.; Sirignano, E.; Lenaerts, B.; Saeys, W. Online crop height and density estimation in grain fields using LiDAR. Biosyst. Eng. 2020, 198, 1–14. [Google Scholar] [CrossRef]
- Zhang, L.; Grift, T. A Lidar-based crop height measurement system for miscanthus giganteus. Comput. Electron. Agr. 2012, 85, 70–76. [Google Scholar] [CrossRef]
- Hämmerle, M.; Höfle, B. Effects of reduced terrestrial Lidar point density on high-resolution grain crop surface models in precision agriculture. Mol. Divers. Preserv. Int. 2014, 14, 24212–24230. [Google Scholar] [CrossRef] [Green Version]
- Zhu, H. Development of UAV-Based Lidar Crop Height Mapping System. Master’s Thesis, University of Illinois at Urbana-Champaign, Urbana, IL, USA, 2017. [Google Scholar]
- Phan, A.; Takahashi, K.; Rikimaru, A. Method for estimating rice plant height without ground surface detection using Laser scanner measurement. J. Appl. Remote Sens. 2016, 10, 046018. [Google Scholar] [CrossRef] [Green Version]
- Anice, C.; Olivier, N.; Mathilde, C.; Sebastien, C.; Frederic, L.; James, T. An algorithm to automate the filtering and classifying of 2D LiDAR data for site-specific estimations of canopy height and width in vineyards. Biosyst. Eng. 2020, 200, 450–465. [Google Scholar] [CrossRef]
- Tilly, N.; Hoffmeister, D.; Cao, Q.; Huang, S.Y.; Lenz-Wiedemann, V.; Miao, Y.X.; Bareth, G. Multitemporal crop surface models: Accurate plant height measurement and biomass estimation with terrestrial laser scanning in paddy rice. J. Appl. Remote Sens. 2014, 8, 83671. [Google Scholar] [CrossRef]
- Jimenez-Berni, J.A.; Deery, D.M.; Rozas-Larraondo, P.; Condon, A.G.; Rebetzke, G.J.; James, R.A.; Bovill, W.D.; Furbank, R.T.; Sirault, X.R.R. High Throughput Determination of Plant Height, Ground Cover, and Above-Ground Biomass in Wheat with LiDAR. Front. Plant Sci. 2018, 9, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Sritarapipat, T.; Rakwatin, P.; Kasetkasem, T. Automatic Rice Crop Height Measurement Using a Field Server and Digital Image Processing. Sensors 2014, 14, 900–926. [Google Scholar] [CrossRef] [Green Version]
- Rogers, D.R., III. Unmanned Aerial System for Monitoring Crop Status. Ph.D. Thesis, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 2014. [Google Scholar]
- Liu, K. Research on the Application of LiDAR Technology in Crop Geometric Feature Detection; JiangSu University: Zhenjiang, China, 2019. [Google Scholar]
- Dhami, H.; Yu, K.; Xu, T.; Zhu, Q.; Dhakal, K.; Friel, J.; Li, S.; Tokekar, P. Crop height and plot estimation from Unmanned Aerial Vehicles using 3D LiDAR. arXiv 2019, arXiv:1910.14031. [Google Scholar]
- Hall, D.S.; Kerstens, P.J.; Rekow, M.N. Integrated Illumination and Detection for LIDAR Based 3-D Imaging. U.S. Patent No. 030780, 10 July 2018. [Google Scholar]
- Xu, J.Z.; Lai, X.D.; Zhang, S.W. Filtering of LIDAR Points Clouds Based on Multi-Resolution Directional Prediction; Geomatics and Information Science of Wuhan University: Wuhan, China, 2007. [Google Scholar]
- Dimitrievski, M.; Veelaert, P.; Philips, W. Semantically aware multilateral filter for depth up sampling in automotive LiDAR point clouds. In Proceedings of the 2017 IEEE Intelligent Vehicles Symposium (IV), Los Angeles, CA, USA, 11–14 June 2017. [Google Scholar] [CrossRef] [Green Version]
- Haris, B.; Jasmin, V.; Walter, B.; Geert, D.C.; Bruno, S. Fast Statistical Outlier Removal Based Method for Large 3D Point Clouds of Outdoor Environments. IFAC-PapersOnLine 2018, 51, 348–353. [Google Scholar] [CrossRef]
- Rusu, R.B. Semantic 3D Object Maps for Every-Day Manipulation in Human Living Environments. Ph.D. Thesis, University of Munich, Munich, Germany, 2009. [Google Scholar]
- Wang, Q. The volatility forecast based on the leverage effect CARR model. Appl. Stat. Manag. 2016, 1, 51–58. [Google Scholar]
- James, W.; James, E.; Glenn, M.; Haydn, K. Photogrammetry for the estimation of wheat biomass and harvest index. Field Crop. Res. 2018, 216, 165–174. [Google Scholar] [CrossRef]
- Chen, B.; Cao, Q.X. Grain height and boundary detection based on laser and IMU. Mechatronics 2020, 1, 6. [Google Scholar]
Performance Target | Parameter | Performance Target | Parameter |
---|---|---|---|
Number of laser lines | 16 | Horizontal angle resolution (°) | 0.1–0.4 |
Range (mm) | <100 | Frequency (Hz) | 5–20 |
Accuracy (cm) | ±3 | Working voltage (V) | 9–32 |
Vertical angle resolution (°) | 2 | Weight (g) | 830 |
Horizontal measurement angle range (°) | 360 | Dimensions (mm) | 103 × 72 |
Wavelength (nm) | 903 | Power (W) | 8 |
Laser beam size (mm) | 9.5 × 12.7 | Divergence angle (°) | 0.18 (3.0 mrad) |
Serial Number | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
1 | 775 | 700 | 780 | 790 | 795 | 793 | 765 | 790 | 780 | 795 |
2 | 750 | 750 | 760 | 800 | 820 | 800 | 780 | 785 | 790 | 790 |
3 | 790 | 810 | 810 | 810 | 800 | 805 | 800 | 810 | 805 | 810 |
4 | 810 | 815 | 810 | 830 | 845 | 845 | 845 | 830 | 830 | 820 |
5 | 820 | 860 | 860 | 875 | 890 | 890 | 890 | 830 | 820 | 800 |
6 | 830 | 835 | 840 | 820 | 820 | 885 | 895 | 840 | 830 | 815 |
7 | 840 | 845 | 845 | 845 | 850 | 850 | 850 | 850 | 850 | 815 |
8 | 830 | 850 | 835 | 825 | 800 | 810 | 810 | 860 | 760 | 800 |
9 | 800 | 800 | 820 | 810 | 800 | 830 | 825 | 815 | 815 | 790 |
10 | 785 | 770 | 775 | 770 | 800 | 760 | 780 | 765 | 780 | 790 |
Group Number | RMSE | MAPE/% | R2 |
---|---|---|---|
Group 1 | 12.18 | 13.36 | 0.974 |
Group 2 | 9.26 | 11.05 | 0.886 |
Group 3 | 9.42 | 9.73 | 0.929 |
Group 4 | 7.37 | 8.51 | 0.853 |
Group 5 | 5.27 | 5.96 | 0.858 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Sun, Y.; Luo, Y.; Zhang, Q.; Xu, L.; Wang, L.; Zhang, P. Estimation of Crop Height Distribution for Mature Rice Based on a Moving Surface and 3D Point Cloud Elevation. Agronomy 2022, 12, 836. https://doi.org/10.3390/agronomy12040836
Sun Y, Luo Y, Zhang Q, Xu L, Wang L, Zhang P. Estimation of Crop Height Distribution for Mature Rice Based on a Moving Surface and 3D Point Cloud Elevation. Agronomy. 2022; 12(4):836. https://doi.org/10.3390/agronomy12040836
Chicago/Turabian StyleSun, Yixin, Yusen Luo, Qian Zhang, Lizhang Xu, Liying Wang, and Pengpeng Zhang. 2022. "Estimation of Crop Height Distribution for Mature Rice Based on a Moving Surface and 3D Point Cloud Elevation" Agronomy 12, no. 4: 836. https://doi.org/10.3390/agronomy12040836
APA StyleSun, Y., Luo, Y., Zhang, Q., Xu, L., Wang, L., & Zhang, P. (2022). Estimation of Crop Height Distribution for Mature Rice Based on a Moving Surface and 3D Point Cloud Elevation. Agronomy, 12(4), 836. https://doi.org/10.3390/agronomy12040836