Next Article in Journal
An Investigation on Pervasive Technologies for IoT-based Thermal Monitoring
Previous Article in Journal
Single-Equipment with Multiple-Application for an Automated Robot-Car Control System
Previous Article in Special Issue
Spectral Data Collection by Dual Field-of-View System under Changing Atmospheric Conditions—A Case Study of Estimating Early Season Soybean Populations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LiDAR-Based 3D Scans of Soil Surfaces and Furrows in Two Soil Types

by
Frederik F. Foldager
1,2,*,
Johanna Maria Pedersen
1,
Esben Haubro Skov
2,
Alevtina Evgrafova
2 and
Ole Green
2
1
Department of Engineering, Aarhus University, Inge Lehmanns Gade 10, 8000 Aarhus C, Denmark
2
Agro Intelligence ApS, Agro Food Park 13, 8200 Aarhus N, Denmark
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(3), 661; https://doi.org/10.3390/s19030661
Submission received: 30 November 2018 / Revised: 23 January 2019 / Accepted: 24 January 2019 / Published: 6 February 2019
(This article belongs to the Special Issue Sensors in Agriculture 2018)

Abstract

:
Soil surface measurements play an important role in the performance assessment of tillage operations and are relevant in both academic and industrial settings. Manual soil surface measurements are time-consuming and laborious, which often limits the amount of data collected. An experiment was conducted to compare two approaches for measuring and analysing the cross-sectional area and geometry of a furrow after a trailing shoe sweep. The compared approaches in this study were a manual pinboard and a Light Detection and Ranging (LiDAR) sensor. The experiments were conducted in coarse sand and loamy sand soil bins exposed to three levels of irrigation. Using the LiDAR, a system for generating 3D scans of the soil surface was obtained and a mean furrow geometry was introduced to study the geometrical variations along the furrows. A comparison of the cross-sectional area measurements by the pinboard and the LiDAR showed up to 41% difference between the two methods. The relation between irrigation and the resulting furrow area of a trailing shoe sweep was investigated using the LiDAR measurements. The furrow cross-sectional area increased by 11% and 34% under 20 mm and 40 mm irrigation compared to non-irrigated in the coarse sand experiment. In the loamy sand, the cross-sectional area increased by 17% and 15% by irrigation of 20 mm and 40 mm compared to non-irrigated measured using the LiDAR.

1. Introduction

A wide range of sensors are used to describe various properties in arable ecosystems. Weed detection can be performed using high precision cameras [1], soil organic carbon stocks and many soil properties can be obtained using spectroscopic techniques [2], as well as soil surface roughness can be evaluated using scanners [3] and digital cameras [4]. Soil surface roughness is an important measure in agricultural research and engineering as it affects properties such as soil water interaction [5] and the risk of soil erosion [6,7]. The soil roughness and cross-sectional shapes of furrows have been determined to develop and assess the performance of tillage tools by measuring the resulting soil disturbance of tillage operations [8,9]. Several methods exist for evaluating the cross-sectional area and geometry of a furrow such as a pinboard [10] and a chain method [11]. These methods are manually operated and thus limited to assess only a finite number of cross-sectional profiles along a furrow. Digital and automated measurements of soil surface roughness illustrated a better description of the spatial variations in both natural and agro-ecosystems [12]. The automated measurement methods also allow for on-the-go measurements of soil roughness as shown in [3] for site-specific cultivation using a Light Detection and Ranging (LiDAR) sensor. The LiDAR has the advantage of being a non-contact, high-resolution technique for spatial data acquisition [13]. Another method for surface detection is the camera-based Structure from Motion (SfM) method [14], where maps are created from a series of pictures. The method has shown its applicability in different scales from mapping landscapes [15] as well as in measuring soil roughness using a commercial grade camera [16].
The aim of this study was to use the 2D LiDAR sensor to assess the furrow geometry and area of a trailing shoe sweep. The trailing shoe is applied to establish a soil furrow to contain the liquid slurry during a slurry application. The use of the trailing shoes is expected to lower emissions in comparison to the slurry application by trailing hoses. The former reduces the surface contact between the slurry and the surrounding atmosphere, which is recognised as a significant factor [17,18,19]. Therefore, accurate measurements of how furrow geometry varies with soil texture and soil moisture content can support predictions of when application by trailing shoe is an efficient emission reduction technique. The objectives were (i) to compare two methods for evaluating soil surface and furrow cross-sectional profiles after a trailing shoe operation using a digital LiDAR sensor and an analogue pinboard and (ii) to measure the changes in furrow geometry and the spatial variations along furrows in coarse sand and loamy sand soils with three different amounts of irrigation.

1.1. Background

1.1.1. Manual Soil Surface Measurements

The relief meter, which was an early version of the pinboard was presented for measuring the soil roughness between cultivation and seedbed preparation [10]. The pinboards vary in size, numbers of pins, spacing between those, materials used, and method of collecting data, but the principle of the measurements is the same. The method is commonly used as it is cost effective and easy to operate. Since the early development, the pinboard method has been used to measure soil surface variations for different purposes [20,21]. Further development of the pinboard for measuring soil surface roughness and soil variations at the microscale was presented [22]. The pinboard has also been used for evaluating soil loss from quarter-drains by quantifying the changes of the cross-sectional area before and after rain events [23].

1.1.2. LiDAR and Camera-Based Soil Surface Sensing

In agricultural engineering, LiDAR-based systems have been applied for various applications such as obstacle detection for autonomous field-robots [24], 3D imaging of crop development [25], and soil surface analysis [6]. The automated systems for soil surface detection have been presented since the 1980s [26]. Moreover, sensor-based systems have been applied for studying microtopography and soil roughness on different scales, from point samples [27,28] to larger areas [29] as well as for airborne applications [30]. The laser-based method performs within the accuracy of the pin displacement unit [31], which made the laser a feasible tool for researching erosion [32]. In addition, microtopography effects of rainfall has been studied using both a photogrammetric technique and a high-resolution laser scanner [33] as well as 1D distance sensors [34,35] and 2D scanners [36] have been applied for soil roughness estimations [37,38]. Besides the use of scanners, camera-based techniques have also been applied to determine soil surface roughness, from hand-held cameras [39] to modern commercial grade cameras [40]. In addition, stereo photogrammetry was also used for mapping soil surfaces [41] and estimating microtopography and soil roughness [42]. Many stationary LiDAR and camera-based systems have been presented for automated soil surface measurements. However, LiDAR was mentioned as a feasible choice for stationary use as well as for on-the-go sensing of soil roughness in the field [3].

2. Materials and Methods

2.1. Experimental Settings

In order to compare two methods such as the pinboard and LiDAR, an experiment consisting of 18 trailing shoe furrows which were obtained in two semi-field soil bins at three irrigation amounts was conducted (Figure 1). For each furrow, a 3D scan and a pinboard measurement of the geometry and area were obtained.
The experiment was conducted in the semi-field facility at Aarhus University Foulum Research Center, Tjele. The semi-field facility was established in 1993 [43], consisting of three soil bins under a moveable roof that allows for controlling the precipitation. Three weeks prior to and during the experiment, the roof was kept over the studied soil bins and the soil was kept bare during the experiments. Two out of three semi-field soil bins were used in this research, coarse sand and loamy sand [44,45,46,47]. The coarse sand (Ortic Haplohumod) had 4% clay (<0.002 mm), 5% silt (0.002–0.063 mm), and 2.0% soil organic matter, while loamy sand (Typic Hapludult) was characterized by 9% clay, 24% silt, and 2.5% soil organic matter [48]. Approximately 24 h prior to the experiment, six plots (1 × 1 m) were irrigated with three different amounts of tap water: 0, 20 and 40 mm. The soil water contents were determined gravimetrically [49] using a soil core of 100 cm 3 (Table 1). The dry bulk densities were 1.28 ± 0.09 g cm 3 for loamy sand and 1.39 ± 0.07 g cm 3 for coarse sand. Due to an unusually dry and hot period prior to the experiments, no significantly different soil water contents were obtained except for the coarse sand soils with 40 mm of irrigation compared to the non-irrigated coarse sand soils (Table 1).
In order to study soil displacements of the furrows, three Bomech trailing shoes (Bomech B.V., Albergen, The Netherlands) were attached to a metal frame on wheels with a distance of 25 cm between each trailing shoe (Figure 2a,b). A vertical load of 117.7 N (12.0 kg) was exerted on each trailing shoe. Three sweeps were conducted for each soil bin at a constant speed of approximately 2.0 km h 1 . In total, 18 furrows were obtained and measured using the analogue pinboard and the digital LiDAR approaches.

2.2. LiDAR-Based Soil Surface Measurement Unit

A soil surface measurement unit was assembled to obtain high-resolution and continuous scans of the soil surface. The 2D LiDAR sensor (SICK LMS511-20100 PRO) was chosen based on the indicated applicability for roughness measurements in the field [50]. The LiDAR operates from a minimum distance to the target at 0.70 m. According to the operating instructions, the statistical error is ±7 mm at 1–10 m distance. The LiDAR was mounted on a linear rail system (Figure 3a) that allowed to obtain the 3D scans of the soil surface by moving the scanner at an elevation of 0.75 m above the surface and parallel to the direction of the furrow while monitoring the position of the scanner. The scans were conducted with an angular resolution of 0.167 and an angular frequency of 25 Hz, while the linear speed of the scanner ( y ˙ ) (Figure 3b) was constant at 0.003 m s 1 , driven by a stepper motor. The distance between two line scans along the furrow (y-direction) was 0.33 mm. Each line scan was 37 cm wide and consisted of 166 data points with a maximum distance between two points of 2.3 mm. The beam diameter increases from 13 mm at the front screen of the sensor to approximately 16.5 mm in diameter at the soil surface of this elevation using the high resolution setting of the LiDAR.
The raw sensor data was collected using the software, SOPAS Engineering Tool ver. 2018.2 (SICK, Germany). No automatic filtering was applied during data sampling. However, the data points were transformed from polar coordinates into cartesian coordinates to study the vertical distances from the scanner to the soil surface in the post-processing and to assemble the 3D scans. A low-pass Gaussian filter ( σ = 1.5 ) was applied in order to remove sensor noise before generating 3D scans using the package: scipy.ndimage.filters.gaussian_filter of SciPy ver. 2.0 in Python 3.6.5.
The x and z coordinates were obtained using the LiDAR sensor, whereas the y coordinates were monitored through the stepper drive (Figure 3b).

2.3. Furrow Cross-Sectional Geometry and Area Measurements

2.3.1. LiDAR Measurements

3D scans of the furrows were obtained using LiDAR measurements and generated based on the filtered data. In order to compare the results of the 3D scanned furrows, a reference frame was introduced to each measured furrow profile. A coordinate system was located in the xz-plane of the furrow. The horizontal xy-plane was located at the elevation of the undisturbed soil. The yz-plane was intersecting with the deepest point of the furrow (Figure 3b). A mean furrow geometry was introduced to express 3D scans as 2D profiles and, hence, to visualise the spatial variations along furrows. The mean furrow was defined for each scan as the mean of the vertical distance measurements (z) for each value of (x) in the relative coordinate systems along the furrow (y). The mean furrow geometries were generated based on 3D scans of which the area was calculated using the package sklearn.metrics.auc ver. 0.19.1 in Python 3.6.5 by integrating the elevations measurements (z) over the distance x l to x r using the trapezoidal rule (Figure 3b). The 2D mean furrow and standard deviations were obtained based on one scan per centimetre of the furrows.

2.3.2. Pinboard Measurements

The analogue measurements of the furrow geometries and areas were performed using pinboard (Figure 4). The pinboard consists of pins (width of 3 mm) that moves freely in the vertical direction. The device was operated by levelling the pins on the soil surface. The contour of the surface was then drawn on paper based on the position of the pins. Subsequent to this operation, the areas were obtained by moving a Digitizing Area-line Meter (Super PLANIX b, Tamaya Technics Inc., Tokyo, Japan) along the drawn contour of the furrow profile, hence referred to as the analogue area measurements. This procedure was done three times for each analogue cross-sectional measurement, and the mean value of these were used for further analysis. In this study, the Computer Aided Design (CAD) software, SOLIDWORKS 2018 SP1 (Dassault Systémes, Paris, France) was applied to digitalise the drawn pinboard-based contours.

3. Results and Discussion

3.1. 2D and 3D Soil Surface Profiles

The 3D scans allowed to evaluate furrow geometries as a result of different irrigation amounts as well as to estimate the variation along furrows within a larger area. The LiDAR gave a more representative description of the furrow compared to analogue measurements as only a limited number of manual measurements can be conducted along the furrow due to the time-consuming process [6]. The cross-sectional geometry, such as width, depth and shape, varied in the two soil types, where the loamy sand plots with no irrigation were characterized with a more well-defined edge and a narrower furrow-width compared to the non-irrigated coarse sand plots (Figure 5 and Figure 6). The furrow geometry in the coarse sand was affected by the increased soil water content. However, in the loamy sand it was not possible to conclude on furrow geometry changes as a factor of soil water content due to the fact that the latter was not significantly different (Table 1). Although, the experiments were performed in semi-field conditions under controlled irrigation, the LiDAR will provide similar results under field conditions. The soil surface in the semi-field is comparable to soil surface conditions in fields, however, a limitation of the LiDAR is the need for bare soil to obtain an uninterrupted line of sight between the sensor and the soil surface.
The furrow profiles of the two soil types under the three irrigation amounts were illustrated in Figure 5 and Figure 6. The subfigures, (a–c) present the 2D profiles of the pinboard measurements, three repetitions each, digitalised using CAD software. The subfigures (d–f) present a subset of the furrow 3D scans and was generated based on the stepper position and LiDAR measurements. In the coarse sand plots, a visual effect of the irrigation was observed by leaving a deeper furrow compared to the non-irrigated soil (Figure 5). This indicates that the resulting furrow geometry was dependent on the water content. Hence, water content and soil type should be considered when the trailing shoe is used for the liquid slurry application.
The mean furrows and the standard deviations of the elevation measurements in the direction of the furrows were determined as the spatial variation along the furrows was observed. Using the mean furrow, it was possible to present 3D scans as 2D cross-sectional plots, which allowed to compare the digital mean furrow plots with the analogue pinboard measurements of the furrow (Figure 7) as the common reference frame was applied (Figure 3b). The location of the common reference frame was a feasible choice for the soil bin study. However, in a field setting, another choice of reference frame could be introduced, e.g., the deepest point of a furrow.
The ability to visualize and take into account the variation along the furrow is needed for research and development purposes related to soil and tillage tools [51]. The magnitude of the standard deviation provides an insight into the variation of the furrow geometry. By evaluating the 3D scans, the three-dimensional soil displacement become available in the assessment of tillage tools [52]. The pinboard has previously been applied for comparing results of simulations-based and measured soil-tool interactions [48]. By considering the variations along the furrow, it is possible to include a physics-based error tolerance when validating numerical models of corresponding soil-tool interactions [53].
The non-contact and automated approach of the LiDAR ensures consistent test results independent from soil conditions, which is an advantage for research purposes. The LiDAR has previously shown applicability in a field setting performing on-the-go measurements of soil surface [50]. This indicates that this technology can be further developed for applications in a field setting as a method for increasing the site-specific treatments based on LiDAR measurements. With respect to the slurry application using trailing shoes, on-board measurements of the furrow geometry and furrow area could act as an input for controlling the pressure on the trailing shoe to obtain the desired furrow area, which is shown to be dependent on the soil type and water content.

3.2. Furrow Area Measurements

The cross-sectional area of the furrow is important due to the application of the trailing shoe, namely to establish a furrow in which the liquid slurry is contained. Furthermore, high variations in the reported effect of trailing shoe on emission reduction were observed [54,55,56,57]. The inconsistencies in reported emission reduction by trailing shoe are often attributed to the soil properties, particularly, soil moisture content [17,56], as this affects the furrow cross-sectional area. Hereby, accurate and consistent methods for assessing the cross-sectional area under different conditions are important as well as being able to correlate the volume of the furrow with the volume of the slurry. The cross-sectional area was a subject for the comparison between the digital and analogue soil surface measurement techniques. Other parameters as width and depth could in a similar way have been considered as the measure for comparison. The cross-sectional areas were measured in every soil plot (Figure 1) in order to determine the change in area due to various irrigation amounts in the two soils. The digital cross-sectional areas were determined using the mean furrow geometry of the 3D scans (Figure 7). A comparison between the analogue profiles and the closest line scan as well as an example of the areas applied for the pinboard measurements are shown in Figure 8.
The analogue cross-sectional area measurements and the LiDAR-based mean furrow areas were shown in Figure 9. The values are based on the mean of the three repetitions for each plot, i.e., combination of soil type and irrigation amount. As an example, for the coarse sand at 40 mm irrigation, the analogue area was calculated as the mean of the area above the red curves of Figure 7, whereas the digital area corresponds to the average area above the three solid blue curves.
There was no significant difference in the areas measured by the digital and analogue methods, except for the non-irrigated coarse sand and coarse sand irrigated at 40 mm (Tukey’s HSD, p < 0.05). However, by comparing the areas of the digital and analogue methods in the coarse sand, it was found that the analogue pinboard predicted lower than the estimations of the digital approach by 7%, 15% and 8% for the 0, 20, and 40 mm irrigation amounts. In the loamy sand plots, all the analogue measurements provided larger cross-sectional areas compared to the digital measurements of 41%, 3% and 16%, respectively. This inconsistency indicated that the results of the analogue measurements was affected by the process of first measuring the surface contour using the pinboard and then calculating the area. By comparing the magnitude of the standard deviations of the analogue area measurements and the standard deviation based on the scanned area, it was found that the standard deviation was lower when LiDAR scans were used, compared to the analogue approach. Hence, the scanner provided more consistent results than the analogue method due to the non-contact and automated process. It was also expected as that the number of data points collected for 3D scans and the mean furrows was higher than the data points collected for the pinboard measurements. By comparing the digital cross-sectional area estimations for the dry coarse sand soil with the irrigated coarse sand of 20 and 40 mm, an increase in furrow cross-sectional area was obtained by 11% and 34%. For the loamy sand, the cross-sectional area increased 17% and 15% by irrigation of 20 and 40 mm. Nevertheless, more experiments under different soil conditions and irrigation amounts are needed to fully evaluate these effects.

3.3. Impact Analysis of Pinboard Measurement

The disturbance of the surface in loose soils during the application of the pinboard [58] was studied by comparing the digital 2D mean furrow geometries and the analogue pinboard measurements. Differences in the furrow geometry between the analogue and digital measurements were observed in the coarse sand plots as a result of the analogue measurement. The impacts were measured by scanning the surface prior to and after the pinboard measurements. In Figure 10, a 3D scan of the soil after the analogue measurement is shown. The positions where the pinboard measurements were taken are marked with arrows (Figure 10). 2D line scans of the furrows in coarse sand with 40 mm irrigation were compared to the corresponding pinboard measurement.
It was observed that the analogue measurements did not accurately capture the geometry of the furrow. Particularly, the features of the furrow for values z > 0, using the notation of Figure 3b, were not captured. Furthermore, it was observed that the analogue measurements provide a wider representation of the furrow cross-sectional geometry compared to the scan conducted before the pinboard measurements (Figure 8).
By the means of the LiDAR, an assessment of the pinboard performance was accessible and an issue of the pinboard was identified, namely, the pinboard approach was challenged in capturing small cross-sectional geometries in loose sand due to disturbance of the soil. This issue has previously been reported as an uncertainty when using the pin-based methods in sandy soil [31].

4. Conclusions

Two methods were compared to evaluate the soil surface and furrow cross-sectional area after a trailing shoe sweep in two soils. A 2D mean profile of the furrow was constructed to compare the digital LiDAR-based results with the results of the analogue pinboard. From the experiments it was concluded:
(i)
The geometric variations in the direction of the furrow were observed by increasing the resolution from 2D pinboard measurements to 3D scans.
(ii)
The results indicated the importance of applying a non-contact method for accurate soil surface measurements in loose soils. The analogue furrow areas were below the LiDAR-based areas in coarse sand by up to 15% and for loamy sand, the analogue areas were above the LiDAR-based areas by up to 41%.
(iii)
An increase in the cross-sectional areas was measured using the LiDAR-based method, in coarse sand the area increased by 11% and 34% for the 20 and 40 mm irrigation compared to the dry coarse sand. In the loamy sand, the cross-sectional area of the furrow increased 17% and 15% by irrigation of 20 and 40 mm. However, more experiments are needed to fully evaluate these effects.

Author Contributions

Conceptualization, F.F.F. and J.M.P.; Methodology, F.F.F. and J.M.P.; Software, E.H.S. and F.F.F.; Experiment design and execution, E.H.S., J.M.P. and F.F.F.; Investigation, F.F.F., J.M.P. and A.E.; Resources, O.G. and J.M.P.; Data Curation, F.F.F., E.H.S. and J.M.P.; Writing—Original Draft Preparation, F.F.F., J.M.P. and A.E.; Writing—Review & Editing, F.F.F., J.M.P., A.E., E.H.S., and O.G.; Visualization, E.H.S.; Supervision, F.F.F. and O.G.; Project Administration, F.F.F.; Funding Acquisition, F.F.F. and O.G.

Funding

This work is partly funded by the Danish Innovation Fund under grant number 7038-00231B and by the Ministry of Environment and Food of Denmark as a green development and demonstration program (GUDP) with the project title New application method for slurry in growing crops.

Acknowledgments

We would like to acknowledge all the technical staff for the support in and contributions to the performance of the tests and Henri Bonnart for assisting in assembling the test setup.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tellaeche, A.; BurgosArtizzu, X.P.; Pajares, G.; Ribeiro, A.; Fernández-Quintanilla, C. A new vision-based approach to differential spraying in precision agriculture. Comput. Electron. Agric. 2008, 60, 144–155. [Google Scholar] [CrossRef]
  2. Viscarra Rossel, R.A.; Lobsey, C.R.; Sharman, C.; Flick, P.; McLachlan, G. Novel proximal sensing for monitoring soil organic C stocks and condition. Environ. Sci. Technol. 2017, 51, 5630–5641. [Google Scholar] [CrossRef] [PubMed]
  3. Jensen, T.; Karstoft, H.; Green, O.; Munkholm, L.J. Assessing the effect of the seedbed cultivator leveling tines on soil surface properties using laser range scanners. Soil Tillage Res. 2017, 167, 54–60. [Google Scholar] [CrossRef]
  4. Moreno, R.G.; Álvarez, M.C.; Alonso, A.T.; Barrington, S.; Requejo, A.S. Tillage and soil type effects on soil surface roughness at semiarid climatic conditions. Soil Tillage Res. 2008, 98, 35–44. [Google Scholar] [CrossRef]
  5. Bauer, T.; Strauss, P.; Grims, M.; Kamptner, E.; Mansberger, R.; Spiegel, H. Long-term agricultural management effects on surface roughness and consolidation of soils. Soil Tillage Res. 2015, 151, 28–38. [Google Scholar] [CrossRef]
  6. Thomsen, L.M.; Baartman, J.E.M.; Barneveld, R.J.; Starkloff, T.; Stolte, J. Soil surface roughness: comparing old and new measuring methods and application in a soil erosion model. Soil 2015, 1, 399–410. [Google Scholar] [CrossRef]
  7. Jester, W.; Klik, A. Soil surface roughness measurement - Methods, applicability, and surface representation. Catena 2005, 64, 174–192. [Google Scholar] [CrossRef]
  8. Milkevych, V.; Munkholm, L.J.; Chen, Y.; Nyord, T. Modelling approach for soil displacement in tillage using discrete element method. Soil Tillage Res. 2018, 183, 60–71. [Google Scholar] [CrossRef]
  9. Zhang, X.; Chen, Y. Soil disturbance and cutting forces of four different sweeps for mechanical weeding. Soil Tillage Res. 2017, 168, 167–175. [Google Scholar] [CrossRef]
  10. Kuipers, H. A relief meter for soil cultivation studies. Neth. J. Agric. Sci. 1957, 5, 255–262. [Google Scholar]
  11. Saleh, A. Soil roughness measurement: chain method. J. Soil Water Conserv. 1993, 48, 527–529. [Google Scholar]
  12. Wolock, D.M.; Price, C.V. Effects of digital elevation model map scale and data resolution on a topography-based watershed model. Water Resour. Res. 1994, 30, 3041–3052. [Google Scholar] [CrossRef]
  13. Oz, I.; Arav, R.; Filin, S.; Assouline, S.; Furman, A. High-Resolution Measurement of Topographic Changes in Agricultural Soils. Vadose Zone J. 2017, 16. [Google Scholar] [CrossRef]
  14. Westoby, M.J.; Brasington, J.; Glasser, N.F.; Hambrey, M.J.; Reynolds, J.M. ‘Structure-from-Motion’ photogrammetry: A low-cost, effective tool for geoscience applications. Geomorphology 2012, 179, 300–314. [Google Scholar] [CrossRef]
  15. Eltner, A.; Maas, H.G.; Faust, D. Soil micro-topography change detection at hillslopes in fragile Mediterranean landscapes. Geoderma 2018, 313, 217–232. [Google Scholar] [CrossRef]
  16. Snapir, B.; Hobbs, S.; Waine, T.W. Roughness measurements over an agricultural soil surface with Structure from Motion. ISPRS J. Photogramm. Remote Sens. 2014, 96, 210–223. [Google Scholar] [CrossRef]
  17. Smith, K.A.; Jackson, D.R.; Misselbrook, T.H.; Pain, B.F.; Johnson, R.A. Reduction of ammonia emission by slurry application techniques. J. Agric. Eng. Res. 2000, 77, 277–287. [Google Scholar] [CrossRef]
  18. Wulf, S.; Maeting, M.; Clemens, J. Application technique and slurry co-fermentation effects on ammonia, nitrous oxide, and methane emissions after spreading. J. Environ. Qual. 2002, 31, 1795–1801. [Google Scholar] [CrossRef]
  19. Webb, J.; Sørensen, P.; Velthof, G.; Amon, B.; Pinto, M.; Rodhe, L.; Salomon, E.; Hutchings, N.; Burczyk, P.; Reid, J. An assessment of the variation of manure nitrogen efficiency throughout Europe and an appraisal of means to increase manure-N efficiency. In Advances in Agronomy; Elsevier: Amsterdam, The Netherlands, 2013; Volume 119, pp. 371–442. [Google Scholar]
  20. Ram, D.N.; Zwerman, P. A convenient soil compaction integrator. Agron. J. 1960, 52. [Google Scholar] [CrossRef]
  21. Burwell, R.E.; Allmaras, R.R.; Amemiya, M. A field measurement of total porosity and surface microrelief of soils. Soil Sci. Soc. Am. Proc. 1963, 27, 697–700. [Google Scholar] [CrossRef]
  22. Moreno, R.G.; Álvarez, M.C.D.; Requejo, A.S.; Tarquis, A.M. Multifractal Analysis of Soil Surface Roughness. Vadose Zone J. 2008, 7, 512. [Google Scholar] [CrossRef]
  23. Kornecki, T.S.; Fouss, J.L.; Prior, S.A. A portable device to measure soil erosion/deposition in quarter-drains. Soil Use Manag. 2008, 24, 401–408. [Google Scholar] [CrossRef]
  24. Foldager, F.F.; Larsen, P.G.; Green, O. Development of a Driverless Lawn Mower Using Co-Simulation; Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics); Springer: Berlin/Heidelberg, Germany, 2018; Volume 10729 LNCS, pp. 330–344. [Google Scholar]
  25. Garrido, M.; Paraforos, D.S.; Reiser, D.; Arellano, M.V.; Griepentrog, H.W.; Valero, C. 3D maize plant reconstruction based on georeferenced overlapping lidar point clouds. Remote Sens. 2015, 7, 17077–17096. [Google Scholar] [CrossRef]
  26. Harral, B.B.; Cove, C.A. Development of an optical displacement transducer for the measurement of soil surface profiles. J. Agric. Eng. Res. 1982, 27, 421–429. [Google Scholar] [CrossRef]
  27. Frede, H.G.; Gäth, S. Soil surface roughness as the result of aggregate size distribution 1. report: Measuring and evaluation method. Z. Pflanzenernähr. Bodenkd. 1995, 158, 31–35. [Google Scholar] [CrossRef]
  28. Darboux, F.; Huang, C.H. An Instantaneous-Profile Laser Scanner to Measure Soil Surface Microtopography. Soil Sci. Soc. Am. J. 2003, 67, 92. [Google Scholar] [CrossRef]
  29. Verhoest, N.E.; Lievens, H.; Wagner, W.; Álvarez-Mozos, J.; Moran, M.S.; Mattia, F. On the soil roughness parameterization problem in soil moisture retrieval of bare surfaces from synthetic aperture radar. Sensors 2008, 8, 4213–4248. [Google Scholar] [CrossRef]
  30. Turner, R.; Panciera, R.; Tanase, M.A.; Lowell, K.; Hacker, J.M.; Walker, J.P. Estimation of soil surface roughness of agricultural soils using airborne LiDAR. Remote Sens. Environ. 2014, 140, 107–117. [Google Scholar] [CrossRef]
  31. Rice, C.; Wilson, B.N.; Appleman, M. Soil topography measurements using image processing techniques. Comput. Electron. Agric. 1988, 3, 97–107. [Google Scholar] [CrossRef]
  32. Eltz, F.L.F.; Norton, L.D. Surface roughness changes as affected by rainfall erosivity, tillage, and canopy cover. Soil Sci. Soc. Am. J. 1997, 61, 1746–1755. [Google Scholar] [CrossRef]
  33. Nouwakpo, S.K.; Huang, C.H. A Simplified Close-Range Photogrammetric Technique for Soil Erosion Assessment. Soil Sci. Soc. Am. J. 2012, 76, 70. [Google Scholar] [CrossRef]
  34. Raper, R.L.; Grift, T.E.; Tekeste, M.Z. A Portable Tillage Profiler for Measuring Subsoiling Disruption. Trans. ASAE 2004, 47, 23–27. [Google Scholar] [CrossRef]
  35. Arvidsson, J.; Bölenius, E. Effects of soil water content during primary tillage—laser measurements of soil surface changes. Soil Tillage Res. 2006, 90, 222–229. [Google Scholar] [CrossRef]
  36. Jensen, T.; Green, O.; Munkholm, L.J.; Karstoft, H. Fourier and granulometry methods on 3D images of soil surfaces for evaluating soil aggregate size distribution. Appl. Eng. Agric. 2016, 32, 609–615. [Google Scholar] [CrossRef]
  37. Huang, C.H.; Bradford, J.M. Portable Laser Scanner for Measuring Soil Surface Roughness. Soil Sci. Soc. Am. J. 1990, 54, 1402. [Google Scholar] [CrossRef]
  38. Huang, C.; Bradford, J. Applications of a laser scanner to quantify soil microtopography. Soil Sci. Soc. Am. J. 1992, 56, 14–21. [Google Scholar] [CrossRef]
  39. Warner, W.S. Mapping a three-dimensional soil surface with hand-held 35 mm photography. Soil Tillage Res. 1995, 34, 187–197. [Google Scholar] [CrossRef]
  40. Rieke-Zapp, D.H.; Nearing, M.A. Digital close range photogrammetry for measurement of soil erosion. Photogramm. Rec. 2005, 20, 69–87. [Google Scholar] [CrossRef]
  41. Taconet, O.; Ciarletti, V. Estimating soil roughness indices on a ridge-and-furrow surface using stereo photogrammetry. Soil Tillage Res. 2007, 93, 64–76. [Google Scholar] [CrossRef]
  42. Gilliot, J.M.; Vaudour, E.; Michelin, J. Soil surface roughness measurement: A new fully automatic photogrammetric approach applied to agricultural bare fields. Comput. Electron. Agric. 2017, 134, 63–78. [Google Scholar] [CrossRef]
  43. Nyord, T.; Kristensen, E.F.; Munkholm, L.J.; Jørgensen, M.H. Design of a slurry injector for use in a growing cereal crop. Soil Tillage Res. 2010, 107, 26–35. [Google Scholar] [CrossRef]
  44. Manevski, K.; Børgesen, C.D.; Andersen, M.N.; Kristensen, I.S. Reduced nitrogen leaching by intercropping maize with red fescue on sandy soils in North Europe: A combined field and modeling study. Plant Soil 2015, 388, 67–85. [Google Scholar] [CrossRef]
  45. Olesen, J.E.; Askegaard, M.; Rasmussen, I.A. Design of an organic farming crop-rotation experiment. Acta Agric. Scand. Sect. B Soil Plant Sci. 2000, 50, 13–21. [Google Scholar] [CrossRef]
  46. Pandey, A.; Li, F.; Askegaard, M.; Rasmussen, I.A.; Olesen, J.E. Nitrogen balances in organic and conventional arable crop rotations and their relations to nitrogen yield and nitrate leaching losses. Agric. Ecosyst. Environ. 2018, 265, 350–362. [Google Scholar] [CrossRef]
  47. De Jonge, L.W.; Sommer, S.G.; Jacobsen, O.H.; Djurhuus, J. Infiltration of slurry liquid and ammonia volatilization from pig and cattle slurry applied to harrowed and stubble soils. Soil Sci. 2004, 169, 729–736. [Google Scholar] [CrossRef]
  48. Chen, Y.; Munkholm, L.J.; Nyord, T. A discrete element model for soil–sweep interaction in three different soils. Soil Tillage Res. 2013, 126, 34–41. [Google Scholar] [CrossRef]
  49. Burt, R. Soil Survey Field and Laboratory Methods Manual; Soil Survey Investigations Report No. 51, Version 1.0; United States Department of Agriculture: Washington, DC, USA, 2009; p. 407.
  50. Jensen, T.; Munkholm, L.; Green, O.; Karstoft, H. A mobile surface scanner for soil studies. In Proceedings of the Second International Conference on Robotics, Associated High-Technologies and Equipment for Agriculture and Forestry—RHEA 2014, Madrid, Spain, 21–23 May 2014; pp. 187–194. [Google Scholar]
  51. Znova, L.; Melander, B.; Lisowski, A.; Klonowski, J.; Chlebowski, J.; Edwards, G.T.; Kirkegaard Nielsen, S.; Green, O. A new hoe share design for weed control: measurements of soil movement and draught forces during operation. Acta Agric. Scand. Sect. B Soil Plant Sci. 2018, 68, 139–148. [Google Scholar] [CrossRef]
  52. Solhjou, A.; Fielke, J.M.; Desbiolles, J.M. Soil translocation by narrow openers with various rake angles. Biosyst. Eng. 2012, 112, 65–73. [Google Scholar] [CrossRef]
  53. Shmulevich, I. State of the art modeling of soil-tillage interaction using discrete element method. Soil Tillage Res. 2010, 111, 41–53. [Google Scholar] [CrossRef]
  54. Hafner, S.D.; Pacholski, A.; Bittman, S.; Burchill, W.; Bussink, W.; Chantigny, M.; Carozzi, M.; Génermont, S.; Häni, C.; Hansen, M.N.; et al. The ALFAM2 database on ammonia emission from field-applied manure: Description and illustrative analysis. Agric. Forest Meteorol. 2018, 258, 66–79. [Google Scholar] [CrossRef]
  55. Webb, J.; Pain, B.; Bittman, S.; Morgan, J. The impacts of manure application methods on emissions of ammonia, nitrous oxide and on crop response—A review. Agric. Ecosyst. Environ. 2010, 137, 39–46. [Google Scholar] [CrossRef]
  56. Rochette, P.; Guilmette, D.; Chantigny, M.H.; Angers, D.A.; MacDonald, J.D.; Bertrand, N.; Parent, L.É.; Côté, D.; Gasser, M.O. Ammonia volatilization following application of pig slurry increases with slurry interception by grass foliage. Can. J. Soil Sci. 2008, 88, 585–593. [Google Scholar] [CrossRef]
  57. Häni, C.; Sintermann, J.; Kupper, T.; Jocher, M.; Neftel, A. Ammonia emission after slurry application to grassland in Switzerland. Atmos. Environ. 2016, 125, 92–99. [Google Scholar] [CrossRef]
  58. Khorashahi, J.; Byler, R.K.; Dillaha, T.A. An opto-electronic soil profile meter. Comput. Electron. Agric. 1987, 2, 145–155. [Google Scholar] [CrossRef]
Figure 1. Experimental design sketch. 18 furrows were obtained and measured using the pinboard and LiDAR within six plots of two soil types and three irrigation amounts.
Figure 1. Experimental design sketch. 18 furrows were obtained and measured using the pinboard and LiDAR within six plots of two soil types and three irrigation amounts.
Sensors 19 00661 g001
Figure 2. Frame with three of the Bomech trailing shoes (a) and three furrows in coarse sand at 40 mm irrigation (b).
Figure 2. Frame with three of the Bomech trailing shoes (a) and three furrows in coarse sand at 40 mm irrigation (b).
Sensors 19 00661 g002
Figure 3. The LiDAR-based soil surface measurement unit in the semi-field facility (a); Sketch of the furrow geometry that includes the reference frame of the cross-sectional profile (b), solid line indicates the soil surface and dashed lines indicate the vertical and horizontal axis of the origin.
Figure 3. The LiDAR-based soil surface measurement unit in the semi-field facility (a); Sketch of the furrow geometry that includes the reference frame of the cross-sectional profile (b), solid line indicates the soil surface and dashed lines indicate the vertical and horizontal axis of the origin.
Sensors 19 00661 g003
Figure 4. The analogue measurement using pinboard of the furrow in loamy sand at 40 mm irrigation.
Figure 4. The analogue measurement using pinboard of the furrow in loamy sand at 40 mm irrigation.
Sensors 19 00661 g004
Figure 5. The cross-sectional measurements in coarse sand after irrigation of 0 mm(a,d), 20 mm (b,e), and 40 mm (c,f). Pinboard measurements (ac), where three repetitions are shown as red, blue and green lines and LiDAR measurements (df), where the colourbar represent elevation (z) in mm.
Figure 5. The cross-sectional measurements in coarse sand after irrigation of 0 mm(a,d), 20 mm (b,e), and 40 mm (c,f). Pinboard measurements (ac), where three repetitions are shown as red, blue and green lines and LiDAR measurements (df), where the colourbar represent elevation (z) in mm.
Sensors 19 00661 g005
Figure 6. The cross-sectional measurements in loamy sand after irrigation of 0 mm (a,d), 20 mm (b,e), and 40 mm (c,f). Pinboard measurements (ac), where three repetitions are shown as red, blue and green lines and LiDAR measurements (df), the colourbar represent elevation (z) in mm.
Figure 6. The cross-sectional measurements in loamy sand after irrigation of 0 mm (a,d), 20 mm (b,e), and 40 mm (c,f). Pinboard measurements (ac), where three repetitions are shown as red, blue and green lines and LiDAR measurements (df), the colourbar represent elevation (z) in mm.
Sensors 19 00661 g006
Figure 7. Mean profiles (solid blue lines), standard deviation on the elevation measurement (light blue area), analogue pinboard measurement (red line). Results are shown for all repetitions in coarse sand at 40 mm irrigation.
Figure 7. Mean profiles (solid blue lines), standard deviation on the elevation measurement (light blue area), analogue pinboard measurement (red line). Results are shown for all repetitions in coarse sand at 40 mm irrigation.
Sensors 19 00661 g007
Figure 8. Cross-sectional geometries in coarse sand at irrigation amount of 40 mm. Line scans prior to the pinboard measurement (green line), line scan subsequent to the pinboard measurement (blue), and analogue pinboard measurement (red).
Figure 8. Cross-sectional geometries in coarse sand at irrigation amount of 40 mm. Line scans prior to the pinboard measurement (green line), line scan subsequent to the pinboard measurement (blue), and analogue pinboard measurement (red).
Sensors 19 00661 g008
Figure 9. Cross-sectional area measurements based on pinboard (Analogue) and LiDAR scans (Digital) in the two soils.
Figure 9. Cross-sectional area measurements based on pinboard (Analogue) and LiDAR scans (Digital) in the two soils.
Sensors 19 00661 g009
Figure 10. 3D scan subsequent to a pinboard measurement. The arrows indicate the location of two pinboard measurements. Experiments conducted in coarse sand, irrigation 40 mm.
Figure 10. 3D scan subsequent to a pinboard measurement. The arrows indicate the location of two pinboard measurements. Experiments conducted in coarse sand, irrigation 40 mm.
Sensors 19 00661 g010
Table 1. Soil water content [g/g] (n = 3) with the standard deviation after three irrigation amounts such as 0, 20, and 40 mm were measured in the topsoils (0–10 cm). The letters indicate significance based on the Tukey’s HDS test, p ≤ 0.05.
Table 1. Soil water content [g/g] (n = 3) with the standard deviation after three irrigation amounts such as 0, 20, and 40 mm were measured in the topsoils (0–10 cm). The letters indicate significance based on the Tukey’s HDS test, p ≤ 0.05.
Irrigation Level, mmLoamy Sand Coarse Sand
00.13 (0.02)a0.06 (0.00)a
200.16 (0.00)a0.08 (0.01)ab
400.15 (0.03)a0.10 (0.02)b

Share and Cite

MDPI and ACS Style

Foldager, F.F.; Pedersen, J.M.; Haubro Skov, E.; Evgrafova, A.; Green, O. LiDAR-Based 3D Scans of Soil Surfaces and Furrows in Two Soil Types. Sensors 2019, 19, 661. https://doi.org/10.3390/s19030661

AMA Style

Foldager FF, Pedersen JM, Haubro Skov E, Evgrafova A, Green O. LiDAR-Based 3D Scans of Soil Surfaces and Furrows in Two Soil Types. Sensors. 2019; 19(3):661. https://doi.org/10.3390/s19030661

Chicago/Turabian Style

Foldager, Frederik F., Johanna Maria Pedersen, Esben Haubro Skov, Alevtina Evgrafova, and Ole Green. 2019. "LiDAR-Based 3D Scans of Soil Surfaces and Furrows in Two Soil Types" Sensors 19, no. 3: 661. https://doi.org/10.3390/s19030661

APA Style

Foldager, F. F., Pedersen, J. M., Haubro Skov, E., Evgrafova, A., & Green, O. (2019). LiDAR-Based 3D Scans of Soil Surfaces and Furrows in Two Soil Types. Sensors, 19(3), 661. https://doi.org/10.3390/s19030661

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop