Next Article in Journal
The Use of qPCR Reveals a High Frequency of Phytophthora quercina in Two Spanish Holm Oak Areas
Next Article in Special Issue
Spatial Analysis of Temperate Forest Structure: A Geostatistical Approach to Natural Forest Potential
Previous Article in Journal
Analysis of Environment-Marker Associations in American Chestnut
Previous Article in Special Issue
Effects of Plot Positioning Errors on the Optimality of Harvest Prescriptions When Spatial Forest Planning Relies on ALS Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Precision Individual Tree Diameter and Perimeter Estimation from Close-Range Photogrammetry

1
Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, 165 21 Praha 6–Suchdol, Czech Republic
2
Department of Forest Management and Geodesy, Faculty of Forestry, Technical University in Zvolen, 96053 Zvolen, Slovakia
3
Department of Economics and Management of Forestry, Faculty of Forestry, Technical University in Zvolen, 96053 Zvolen, Slovakia
4
Department of Forest Harvesting, Logistics and Ameliorations, Faculty of Forestry, Technical University in Zvolen, 96053 Zvolen, Slovakia
*
Author to whom correspondence should be addressed.
Forests 2018, 9(11), 696; https://doi.org/10.3390/f9110696
Submission received: 11 October 2018 / Revised: 7 November 2018 / Accepted: 8 November 2018 / Published: 10 November 2018
(This article belongs to the Special Issue Forest Landscape Management: From Data to Decision)

Abstract

:
Close-range photogrammetry (CRP) can be used to provide precise and detailed three-dimensional data of objects. For several years, CRP has been a subject of research in forestry. Several studies have focused on tree reconstruction at the forest stand, plot, and tree levels. In our study, we focused on the reconstruction of trees separately within the forest stand. We investigated the influence of camera lens, tree species, and height of diameter on the accuracy of the tree perimeter and diameter estimation. Furthermore, we investigated the variance of the perimeter and diameter reference measurements. We chose four tree species (Fagus sylvatica L., Quercus petraea (Matt.) Liebl., Picea abies (L.) H. Karst. and Abies alba Mill.). The perimeters and diameters were measured at three height levels (0.8 m, 1.3 m, and 1.8 m) and two types of lenses were used. The data acquisition followed a circle around the tree at a 3 m radius. The highest accuracy of the perimeter estimation was achieved when a fisheye lens was used at a height of 1.3 m for Fagus sylvatica (root mean square error of 0.25 cm). Alternatively, the worst accuracy was achieved when a non-fisheye lens was used at 1.3 m for Quercus petraea (root mean square error of 1.27 cm). The tree species affected the estimation accuracy for both diameters and perimeters.

1. Introduction

Currently, demands for precise information about individual trees and forest stands are increasing in forestry [1]. The most important parameter for the estimation of forest information is the trunk diameter [2]. The diameter is related strongly to other descriptors of tree size such as total height [3] and crown dimensions [4]. Thus, it is a key predictor to estimate the stem and tree volume [5], biomass [6], and other information [7]. In addition to the forestry practice, the high demand for precise and complex three-dimensional (3D) data is evident in multiple research fields; for example, in functional-structural models [8] where it is crucial to obtain accurate temporal 3D data. Furthermore, the demand for a high-quality visualization of the forest environment is high.
Many types of measuring instruments exist that exhibit different properties (e.g., precision, accuracy, operational simplicity, and cost). The methods used can be divided into two categories: contact methods (tape, caliper, electronic tree measuring fork) and non-contact methods (optical calipers, rangefinder dendrometers, optical forks) [9]. The contact methods allow for the highest accuracy. Owing to time constraints and safety issues, they are limited to measurements of the lower bole. Meanwhile, non-contact methods are inaccurate or expensive in terms of time and instrument expense. Hence, alternative methods must be sought [1].
During the last two decades, development methods have rendered 3D point clouds as a useful data source for tree measurements [2]. The 3D data of trees from the ground may be characterized by three primary acquisition methods: (1) magnetic motion tracker, (2) terrestrial laser scanner (TLS), and (3) close-range photogrammetry (CRP) [10]. The magnetic motion tracker allows for the detailed 3D point clouds to be used as a basis for determining the tree diameter. However, this method is expensive and time-consuming [11]. Terrestrial laser scanning produces extremely high point densities and fills the gap between tree-scale manual measurements and large-scale airborne laser scanning measurements by providing a large amount of precise information on various forest structural parameters [12]. The TLS can provide an estimation accuracy of 1–2 cm when the diameter at breast height (DBH) and the diameter on stem curve are estimated [13]. The disadvantage of the TLS is its relatively high cost and lack of personnel training.
CRP is an approach allowing for the automatic reconstruction of 3D models using sets of overlapping two-dimensional digital images [14]. Terrestrial photogrammetry is a perspective method to create 3D point clouds, because hand-held cameras are appropriate for end-users owing to their low-cost, and low-weight sensors that are highly portable and easy to access [15]. Several authors used terrestrial photogrammetry to determine tree parameters [1,2,7,10,15,16,17,18,19]. For example, Liang et al. [15] described a method for terrestrial photogrammetry using a hand-held consumer camera. Forsman et al. [2] developed a method for terrestrial mapping of a tree structure using a multi-camera rig. Surový et al. [10] described the accuracy of point reconstruction on individual stems using a hand-held camera with stop-and-go. Miller et al. [18] designed and described a method for individual tree models. Bauwens et al. [7] used terrestrial photogrammetry as a tool for the non-destructive method of modelling irregularly shaped tropical tree trunks. However, in these studies, the authors did not focus on trunk perimeter estimation from CRP. The perimeter may be more accurate for deriving other parameters such as the diameter. Furthermore, the impact of different types of bark was not investigated. Most of the previous studies focused on trees within a research plot and created images for the whole research area. Capturing the trees individually within the forest stand can result in better estimation accuracies.
The aim of this study was to design an appropriate CRP method in order to create 3D point clouds of individual trees. Furthermore, we estimated the diameter and perimeter of individual trees with different types of bark, at different heights on the tree stem using a digital single-lens reflex (DSLR) camera with a non-fisheye lens and a fisheye lens, separately. The research itself is a first step towards the workflow development that starts from data acquisition to functional-structural forest modelling and visualization.

2. Materials and Methods

2.1. Study Site

The research area is located in the central part of Slovakia, which belongs to the Kremnica Mountains. The forest stands are unique according to the tree species composition, primarily European beech (Fagus sylvatica L.), sessile oak (Quercus petraea (Matt.) Liebl.), Norway spruce (Picea abies (L.) H. Karst.), and European silver fir (Abies alba Mill.). These four species differ significantly by bark type. Therefore, we chose 10 trees from each tree species. However, the light conditions and slope (0%–5%) were the same for all trees. The age varied from 55 to 80 years. Further, the age was equal, and the trees were located in the same forest stand within each tree species (Figure 1).

2.2. Reference Data Acquisition

Tree diameter measurements were performed using a caliper and measuring tape. The diameter was measured at three heights: 0.8 m, 1.3 m, and 1.8 m. At each of those height levels, we marked the trunk on its perimeter with four markings. The diameter was measured with a caliper to be perpendicular to those markings, as shown in Figure 2. The caliper measurement was performed four times for every marking. After measuring with the caliper, we also measured the perimeter by measuring the tape four times starting from each of the markings, and then we calculated the diameter. Altogether, we measured the diameter 20 times at each of the three height levels on the trunk, thus 60 times per tree.

2.3. Image Acquisition

Images of trunks were captured using a DSLR camera, Canon 70D (Canon Inc., Tokyo, Japan), with a CMOS sensor and 20.2-megapixel resolution. Two different lenses were used. The first was the Canon EF 35 mm f/2 IS USM (Canon Inc., Tokyo, Japan) with fixed focal length (non-fisheye). The second lens was the Canon EF 8–15 mm f/4L Fisheye USM (fisheye) (Canon Inc., Tokyo, Japan). The fisheye lens was fixed to 15 mm.
Before the image acquisition, eight A4 papers formatted with 8-bit markers (two on each paper) were placed on the ground, and one was mounted on the reference stick. The reference stick was placed near the trunk (<1 m distance). These markers were used for the scaling and orientation of the point clouds. Altogether, we used nine such papers. For scaling purposes, only one was required. We used multiple papers for verifying the accuracy. The accuracy achieved was under 1 mm in almost all cases. Each tree was captured twice using the non-fisheye lens and subsequently by the fisheye lens. The distance from the tree was maintained at 3 m. The path was a circle around the tree for both lenses. The images were captured in the static position of the operator, and the distance between the positions was approximately 1 m. Two images per position were captured when the non-fisheye lens was used to include the ground and the highest section (1.8 m). When the fisheye lens was used, only one image was required to capture the ground and the highest section. The average number of positions was 35. The average time to collect the images was 4 min for the non-fisheye lens and 2 min for the fisheye lens.

2.4. Data Processing

Images were processed into scaled and oriented point clouds using Agisoft PhotoScan Professional 1.2.6 software (Agisoft LCC, St. Petersburg, Russia). The images were separated into non-fisheye and fisheye. The images were subsequently split into groups. Each grouping corresponded to one individual tree. Subsequently, the images were aligned with automatic camera calibration. When the non-fisheye lens was used, the “frame” camera type was chosen; for the fisheye lens, the “fisheye” camera type was chosen within the calibration options in Agisoft PhotoScan Professional. For all trees, reports were generated. In the reports is information about camera calibration, image overlap, scaling error, etc. Reports are available as Supplementary Materials. An example of the image overlap and camera locations is shown in Figure 3.
Thereafter, we used an option where each image was compared to all other images within the grouping in full resolution with high accuracy, after the alignment markers were detected automatically. Each A4 paper contained two markers that were used for scaling. First, we paired the markers that shared the A4 paper and subsequently set a scale (Figure 4). The distance between the centers of the markers was 14.2 cm. The orientation of the Z axis was set by the markers placed on the reference stick. Subsequently, we generated the dense point cloud and exported the point cloud to the .txt format for further processing.
First, the trunk cross-sections at 0.8 m, 1.3 m, and 1.8 m on the tree trunk were extracted from the point clouds. DendroCloud software was used for this task [20,21]. To cut the cross-sections at different height levels, a digital elevation model (DEM) was generated. We set 0.5 as the grid size where the minimum Z value was assigned from the point cloud. Based on the DEM cross-section at 1.3 m, it was sliced with a 2 cm thickness. The points within the cross-section were grouped spatially, and a circle fitting algorithm was used to calculate the initial diameter and position of the tree. We used the Monte Carlo method with the automatic initial method. More information about the algorithm is reported in [20]. The normalized DEM was calculated based on the position and diameter, and multiple cross-sections at the desired heights were generated. Subsequently, the points in each height were grouped spatially. Each tree contained three grouped cross-sections (0.8 m, 1.3 m, and 1.8 m) for the non-fisheye lens and fisheye lens. The results were exported to the .csv format and subsequently imported to ArcGIS for Desktop software 10.2 (Environmental Systems Research Institute (ESRI), Redlands, CA, USA). Within ArcGIS for Desktop, the convex hull algorithm was applied by the Minimum Bounding Geometry module. Altogether, 240 perimeters were calculated.
The calculated perimeter values were exported by the ArcGIS for Desktop software in a *.txt file, where the name consisted of the following identifiers: tree species with its number, lens type, and height of the tree referred by the calculated perimeter. The file itself contained the value of the calculated perimeter. We created a script in the Python programming language, version 2.7 (Python Software Foundation, Beaverton, OR, USA), to perform the automated acquisition of the above-mentioned parameters from the file name along with the calculated value of the perimeter as a record for each *.txt file. After processing all files, the script saved the array of records (dataset) into the database for further analysis.
To compare the diameters measured using the caliper, we used the position of the markers that were sprayed on the trunks. In Agisoft PhotoScan Professional, these markers were identified as control points. These markers indicated the position where the diameter was measured. Subsequently, the points were exported to ArcGIS for Desktop. Together with points from the cross-section, we measured the diameter in the same position and direction as in the field using the caliper for each tree and cross-section. The entire workflow is shown in Figure 5.

2.5. Data Evaluation

The estimation error was calculated as the difference between the reference measurement and the estimation from the point cloud. The root mean square error (RMSE) (1) of the diameter and perimeter estimation was also calculated. The error and RMSE were calculated for the estimations from the point clouds generated from both lenses that were used:
RMSE % = RMSE i = 1 N R i N × 100 ,
where Ri represents the reference diameter or perimeter.
Furthermore, the factorial analysis of variance (ANOVA) was performed to detect the effects of lens, tree species, and height of measurements (0.8, 1.3, 1.8) on the estimation error. The factorial ANOVA was used separately for the perimeters and diameters. Additionally, the datasets for the perimeters and diameters were separated based on the type of lens, and factorial ANOVA was used again on these datasets.

3. Results

3.1. Reference Data

The perimeters of the same position on the trunk were measured four times each. We calculated the variance of these measurements to determine the error originating from the operator. Figure 6 shows the graphs of the variance diameters and perimeters separated by the tree species and height on the trunk. Fagus sylvatica exhibited the lowest variance with no significant changes between the heights. Quercus petraea and Picea abies exhibited the highest variances. In most cases, we observed that the variance for a height of 1.3 m was the lowest within the same tree species and type of measurement (tape and caliper).

3.2. Perimeter

The RMSE of the perimeter estimation for all tree species when the non-fisheye lens was used varied from 0.27 cm to 1.27 cm. When the fisheye lens was used, the RMSE varied from 0.25 cm to 1.21 cm (Table 1).
The RMSE% varied from 0.21% to 0.99%. Thus, all RMSE% were under 1% (Table 2).

3.3. Diameter

The RMSE of the diameter estimation for all tree species when the non-fisheye lens was used varied from 0.423 cm to 0.71 cm. When the fisheye lens was used, the RMSE varied from 0.386 cm to 0.596 cm (Table 3).
The RMSE% varied from 0.90% to 1.85%. All RMSE% were under 2% (Table 4).

3.4. Impact of Factors

The factorial ANOVA was used to identify if the tree species, height of measurement, and lens affected the errors of estimation significantly for both types of measurements (perimeter and diameter).
Within almost all cases, the tree species was the only factor that affected the estimated error significantly. The p-value ranged from 3.3 × 10−16 to 0.015. Furthermore, in terms of the diameter, the height of the measurement and the lens affected the estimated diameters significantly (p = 0.037 and p = 0.021). The influence of tree species on the errors is illustrated in Figure 7.

4. Discussion

Research on the application of photogrammetry (both terrestrial and airborne) in forestry has increased after the invention of the structure-from-motion (SfM) and multi-view stereo (MVS) techniques, which are closely related to computer vision and learning [14,18]. Forests are measured photogrammetrically at the stand [1], plot [15,17], and tree levels [10,18]. Our study provides the results of the 3D capture and reconstruction of individual trees that can potentially be applied towards quality and quantity estimations [22], the creation of high-precision tree models for growth models and simulations [23], etc. In addition to SfM and MVS, other photographic methods have been reported to be applicable in the evaluation of spatial relations in forests. For example, Dick et al. [24] used 360° panoramic images to measure forest stands using a relation between a known focal length and partially measured field data. Hapca et al. [25] used two converging digital images to reconstruct standing trees up to a height of 12 m.
The accuracy of the photogrammetric 3D reconstruction of trees has already been studied widely, where the DBH is the most frequently estimated parameter. For the DBH of 20 specimens of Cryptomeria japonica, Surový et al. [10] reported an RMSE of 0.59 cm with an SD of 0.72 cm, thus achieving sub-centimeter accuracy as reported in the current study. According to ANOVA, the tree species affected the accuracy significantly; however, these differences were only in the millimeter scale. Liang et al. [15] reported a DBH RMSE of 2.39 cm; however, the data acquisition was plotted based on a 30 m × 30 m rectangular plot. In another study [26], an average DBH error of 1.15 cm was reported for 500 m2 circular plots.
For a fisheye camera, Berveglieri et al. [16] reported an average DBH error of 1.46 cm with a standard deviation of 1.09 cm for a small forest plot consisting of seven trees. They considered these results comparable with manual measurements and values generated from laser scanning. Rodríguez-García et al. [27] studied an application of stereoscopic hemispherical fisheye images and achieved a DBH RMSE of 1.51 cm. Fisheye lenses allow for the coverage of a large area, thereby requiring a lesser number of images. However, their imaging geometry causes large variations in scale, and must be addressed using proper calibration. Some authors [28,29,30] described this calibration process. For example, Berveglieri et al. [16] calibrated a fisheye camera in a terrestrial calibration field composed of coded targets. Our results indicate that, currently, even the “on-the-job” calibration in specialized software can provide results comparable with non-wide-angle cameras. Furthermore, the reprojection errors (Supplementary Materials) ranged from 0.70 pix. to 0.95 pix. for non-fisheye lens and from 0.42 pix. to 0.57 pix. for fisheye lens (Table 5). “Reprojection error is the distance between the point on the image where a reconstructed 3D point can be projected and the original projection of that 3D point detected on the photo and used as a basis for the 3D point reconstruction procedure” [31]. Reprojection errors of fisheye lens are lower for each tree species. Reader should consider that the impact of the automatic camera calibration was not investigated.
In the present study, diameters and perimeters were tested at three heights up to 1.8 m. However, measurements of longer sections of the stem are often required. Bauwens et al. [7] conducted a research where tropical, irregularly shaped trees with massive buttresses were reconstructed up to 12 m in height with an RMSE under 5%. In terms of error variability with increasing height, Surový et al. [10] observed an increased variability in reconstruction errors only on the lowest part of the stem, 10 cm above ground. This likely demonstrates an advantage of the current CRP methods, as a higher error with increasing height was reported for older photogrammetric methods [25].
Laser scanning methods might still be more precise and less laborious, but the high price of the devices is restrictive with respect to a wider application [10]. In addition, if only basic parameters (e.g., tree position and DBH) are required, these can be measured quickly using simple equipment (compass, rangefinder, caliper). The applicability of “classical” methods is also confirmed in our study where references for both perimeters and diameters were measured with a variability of less than one centimeter. The advantage of point cloud applications is related to the estimation of advanced attributes such as the stem curve and stem quality [15].
As suggested by the afore-mentioned studies, the accuracy of the current CRP methods is high, and fully applicable for the most typical tasks related to the trees’ dimensions. Although the computational power of the current computers is extremely high, photogrammetrical processing can be highly time-consuming, as it can potentially continue for days. A photogrammetric assessment of the inevitable data required could be crucial for the optimization of the whole process. Miller et al. [18] measured small trees using 150–180 photos per tree. Surový et al. [10] reported a crucial shift in accuracy for five to eight cameras simultaneously based on the same point. With more than eight cameras, the accuracy did not increase. To reconstruct the stem, 40 photos were captured. This corresponds with our study, where 40 images captured by the fisheye lens were used. The number of necessary images captured using the standard lens was approximately two times higher. The distance between the operator and the tree must also be considered. Bauwens et al. [7] processed the photographs of large tropical trees captured with the focal length of 25 mm; however, the processing of images captured with focal lengths of 35 mm and 50 mm failed. This was likely caused by the small distance between the operator and the tree that led to images with almost the same area and texture as those of the tree. The distance could also be crucial for accuracy. James and Robson [14] suggested a ratio of 1:1000 for the accuracy and distance; however, further testing is required. In our case, the DBH RMSEs ranged from ~4 to ~7 mm at a distance of 3 m. The accuracy could have been influenced by the artificial signs on and around the stems, which could facilitate photo alignment.

5. Conclusions

CRP is suitable for generating 3D point clouds of individual trees with a high accuracy of diameter and perimeter estimations of the tree trunk. We obtained the RMSE of perimeter estimation between 0.27 cm and 1.28 cm. Regarding the diameter, the RMSE varied from 0.39 cm to 0.71 cm. Surprisingly, the camera lens type (non-fisheye vs. fisheye) did not affect the estimation accuracy of the tree trunk perimeter significantly. Meanwhile, the tree species affected the accuracy significantly with regard to the perimeter and diameter estimations. The estimation accuracy of the perimeter and diameter of the trees was highly promising for the monitoring of annual tree trunk increment. Further research will focus on the tree trunk annual increment estimation by CRP, regardless of whether the technique is suitable. The increment in the stem curve should also be investigated.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/9/11/696/s1, reports.

Author Contributions

Conceptualization, M.M., J.V. and J.M.; Data curation, M.M., J.V., A.G., P.V. and M.S.; Formal analysis, M.M.; Funding acquisition, M.M. and J.M.; Methodology, M.M. and J.V.; Supervision, J.M.; Validation, M.M.; Writing—original draft, M.M., J.V. and J.T.; Writing—review and editing, M.M., J.V., J.T., P.V. and J.M.

Funding

This work was supported by the Slovak Research and Development Agency through Grant No. APVV-15-0714 (“Mitigation of climate change risk by optimization of forest harvesting scheduling”), by the Scientific Grant Agency of the Ministry of Education, Science, Research, and Sport of the Slovak Republic under the grant VEGA 1/0881/17, and by Grant No. CZ.02.1.01/0.0/0.0/16_019/0000803 (“Advanced research supporting the forestry and wood-processing sector’s adaptation to global change and the 4th industrial revolution”) financed by OP RDE.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mikita, T.; Janata, P.; Surovỳ, P. Forest stand inventory based on combined aerial and terrestrial close-range photogrammetry. Forests 2016, 7, 165. [Google Scholar] [CrossRef]
  2. Forsman, M.; Börlin, N.; Holmgren, J. Estimation of Tree Stem Attributes Using Terrestrial Photogrammetry with a Camera Rig. Forests 2016, 7, 61. [Google Scholar] [CrossRef]
  3. Zhang, L.; Bi, H.; Cheng, P.; Davis, C.J. Modeling spatial variation in tree diameter—Height relationships. For. Ecol. Manag. 2004, 189, 317–329. [Google Scholar] [CrossRef]
  4. Hemery, G.E.; Savill, P.S.; Pryor, S.N. Applications of the crown diameter—Stem diameter relationship for different species of broadleaved trees. For. Ecol. Manag. 2005, 215, 285–294. [Google Scholar] [CrossRef]
  5. Belyea, H.C. Forest Measurement; John Wiley & Sons, Inc.: New York, NY, USA, 1931. [Google Scholar]
  6. Jenkins, C.; Chojnacky, D.; Heath, L.; Birdsey, R. Comprehensive Database of Diameter-Based Biomass Regressions for North American Tree Species; US Department of Agriculture, Forest Service, Northeastern Research Station: Newtown Square, PA, USA, 2004.
  7. Bauwens, S.; Fayolle, A.; Gourlet-Fleury, S.; Ndjele, L.M.; Mengal, C.; Lejeune, P. Terrestrial photogrammetry: A non-destructive method for modelling irregularly shaped tropical tree trunks. Methods Ecol. Evol. 2017, 8, 460–471. [Google Scholar] [CrossRef]
  8. Sievänen, R.; Godin, C.; de Jong, T.M.; Nikinmaa, E. Functional-structural plant models: A growing paradigm for plant studies. Ann. Bot. 2014, 114, 599–603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Clark, N.; Wynne, R.; Schmoldt, D. A review of past research on dendrometers. For. Sci. 2000, 46, 570–576. [Google Scholar]
  10. Surový, P.; Yoshimoto, A.; Panagiotidis, D. Accuracy of reconstruction of the tree stem surface using terrestrial close-range photogrammetry. Remote Sens. 2016, 8, 123. [Google Scholar] [CrossRef]
  11. Yoshimoto, A.; Surovy, P.; Konoshima, M.; Kurth, W. Constructing tree stem form from digitized surface measurements by a programming approach within discrete mathematics. Trees 2014, 1577–1588. [Google Scholar] [CrossRef]
  12. Liang, X.; Kankare, V.; Hyyppä, J.; Wang, Y.; Kukko, A.; Haggrén, H.; Yu, X.; Kaartinen, H.; Jaakkola, A.; Guan, F.; et al. Terrestrial laser scanning in forest inventories. ISPRS J. Photogramm. Remote Sens. 2016, 115, 63–77. [Google Scholar] [CrossRef] [Green Version]
  13. Liang, X.; Hyyppä, J.; Kaartinen, H.; Lehtomäki, M.; Pyörälä, J.; Pfeifer, N.; Holopainen, M.; Brolly, G.; Francesco, P.; Hackenberg, J.; et al. International benchmarking of terrestrial laser scanning approaches for forest inventories. ISPRS J. Photogramm. Remote Sens. 2018, 144, 137–179. [Google Scholar] [CrossRef]
  14. James, M.R.; Robson, S. Straightforward reconstruction of 3D surfaces and topography with a camera : Accuracy and geoscience application. J. Geophys. Res. Earth Surf. 2012, 117, 1–17. [Google Scholar] [CrossRef]
  15. Liang, X.; Jaakkola, A.; Wang, Y.; Hyyppä, J.; Honkavaara, E.; Liu, J.; Kaartinen, H. The use of a hand-held camera for individual tree 3D mapping in forest sample plots. Remote Sens. 2014, 6, 6587–6603. [Google Scholar] [CrossRef]
  16. Berveglieri, A.; Tommaselli, A.; Liang, X.; Honkavaara, E. Photogrammetric measurement of tree stems from vertical fisheye images. Scand. J. For. Res. 2017, 32, 737–747. [Google Scholar] [CrossRef]
  17. Mokroš, M.; Liang, X.; Surový, P.; Valent, P.; Čerňava, J.; Chudý, F.; Tunák, D.; Saloň, Š.; Merganič, J. Evaluation of Close-Range Photogrammetry Image Collection Methods for Estimating Tree Diameters. ISPRS Int. J. Geo-Inf. 2018, 7, 93. [Google Scholar] [CrossRef]
  18. Miller, J.; Morgenroth, J.; Gomez, C. 3D modelling of individual trees using a handheld camera: Accuracy of height, diameter and volume estimates. Urban For. Urban Green. 2015, 14, 932–940. [Google Scholar] [CrossRef]
  19. Roberts, J.W.; Koeser, A.K.; Abd-Elrahman, A.H.; Hansen, G.; Landry, S.M.; Wilkinson, B.E. Terrestrial Photogrammetric Stem Mensuration for Street Trees. Urban For. Urban Green. 2018, 35, 66–71. [Google Scholar] [CrossRef]
  20. Koreň, M.; Mokroš, M.; Bucha, T. Accuracy of tree diameter estimation from terrestrial laser scanning by circle-fitting methods. Int. J. Appl. Earth Obs. Geoinf. 2017, 63, 122–128. [Google Scholar] [CrossRef]
  21. Koreň, M. DendroCloud: Point Cloud Processing Software for Forestry, Version 1.45. Available online: http://gis.tuzvo.sk/dendrocloud/download/dendrocloud_1_45.pdf (accessed on 21 August 2018).
  22. Cieszewski, C.J.; Strub, M.; Antony, F.; Bettinger, P.; Dahlen, J.; Lowe, R.C. Wood quality assessment of tree trunk from the tree branch sample and auxiliary data based on NIR Spectroscopy and SilviScan. Math. Comput. For. Nat. Sci. 2013, 5, 86–111. [Google Scholar]
  23. Fabrika, M.; Valent, P.; Scheer, Ľ. Thinning trainer based on forest-growth model, virtual reality and computer-aided virtual environment. Environ. Model. Softw. 2018, 100, 11–23. [Google Scholar] [CrossRef]
  24. Dick, A.R.; Kershaw, J.A.; MacLean, D.A. Spatial Tree Mapping Using Photography. North. J. Appl. For. 2010, 27, 68–74. [Google Scholar] [CrossRef]
  25. Hapca, A.I.; Mothe, F.; Leban, J.-M. A digital photographic method for 3D reconstruction of standing tree shape. Ann. For. Sci. 2007, 64, 631–637. [Google Scholar] [CrossRef]
  26. Tomaštík, J.; Saloň, Š.; Tunák, D.; Chudý, F.; Kardoš, M. Tango in forests—An initial experience of the use of the new Google technology in connection with forest inventory tasks. Comput. Electron. Agric. 2017, 141, 109–117. [Google Scholar] [CrossRef]
  27. Rodríguez-García, C.; Montes, F.; Ruiz, F.; Canellas, I.; Pita, P. Stem mapping and estimating standing volume from stereoscopic hemispherical images. Eur. J. For. Res. 2014, 133, 895–904. [Google Scholar] [CrossRef]
  28. Marcato Junior, J.; de Moraes, M.V.A.; Tommaselli, A.M.G. Avaliação experimental de técnicas de calibração de câmaras com objetiva olho de peixe. Bol. Ciencias Geod. 2015, 21, 637–651. [Google Scholar] [CrossRef]
  29. Schneider, D.; Schwalbe, E.; Maas, H.-G. Validation of geometric models for fisheye lenses. ISPRS J. Photogramm. Remote Sens. 2009, 64, 259–266. [Google Scholar] [CrossRef]
  30. Garrido-Jurado, S.; Muñoz-Salinas, R.; Madrid-Cuevas, F.J.; Marín-Jiménez, M.J. Automatic generation and detection of highly reliable fiducial markers under occlusion. Pattern Recognit. 2014, 47, 2280–2292. [Google Scholar] [CrossRef]
  31. Agisoft, L.L.C. Agisoft PhotoScan User Manual: Professional Edition, Version 1.1; Agisoft LLC: St. Petersburg, Russia, 2014. [Google Scholar]
Figure 1. Overview of the trees’ positions within the forest stands, and their positions within the Slovak republic are shown below.
Figure 1. Overview of the trees’ positions within the forest stands, and their positions within the Slovak republic are shown below.
Forests 09 00696 g001
Figure 2. (left) Example of the signs for diameter measurement in the field and (right) diagram of diameter measurement based on the sign on the bark.
Figure 2. (left) Example of the signs for diameter measurement in the field and (right) diagram of diameter measurement based on the sign on the bark.
Forests 09 00696 g002
Figure 3. Camera locations and image overlap from the report generated by Agisoft PhotoScan Professional. The same tree is shown on both images (Abies alba Mill.).
Figure 3. Camera locations and image overlap from the report generated by Agisoft PhotoScan Professional. The same tree is shown on both images (Abies alba Mill.).
Forests 09 00696 g003
Figure 4. Scaling and orientation process using the coding markers.
Figure 4. Scaling and orientation process using the coding markers.
Forests 09 00696 g004
Figure 5. Workflow diagram of the point cloud data processing to the resulted perimeters and diameters of trees.
Figure 5. Workflow diagram of the point cloud data processing to the resulted perimeters and diameters of trees.
Forests 09 00696 g005
Figure 6. Variances of measured reference data separated by tree species. Perimeter is shown on the left and diameter on the right (cm) (Aa: Abies alba Mill.; Fs: Fagus sylvatica L.; Qp: Quercus petraea (Matt.) Liebl.; Pa: Picea abies (L.) H. Karst.).
Figure 6. Variances of measured reference data separated by tree species. Perimeter is shown on the left and diameter on the right (cm) (Aa: Abies alba Mill.; Fs: Fagus sylvatica L.; Qp: Quercus petraea (Matt.) Liebl.; Pa: Picea abies (L.) H. Karst.).
Forests 09 00696 g006
Figure 7. Tree species average error (cm) comparison of (first row) estimated perimeter and (second row) diameter using (first column) non-fisheye and (second column fisheye lens) (Aa: Abies alba Mill.; Fs: Fagus sylvatica L.; Qp: Quercus petraea (Matt.) Liebl.; Pa: Picea abies (L.) H. Karst.).
Figure 7. Tree species average error (cm) comparison of (first row) estimated perimeter and (second row) diameter using (first column) non-fisheye and (second column fisheye lens) (Aa: Abies alba Mill.; Fs: Fagus sylvatica L.; Qp: Quercus petraea (Matt.) Liebl.; Pa: Picea abies (L.) H. Karst.).
Forests 09 00696 g007
Table 1. RMSE (root mean square error) of perimeter estimation (cm) (Aa: Abies alba Mill.; Fs: Fagus sylvatica L.; Qp: Quercus petraea (Matt.) Liebl.; Pa: Picea abies (L.) H. Karst.).
Table 1. RMSE (root mean square error) of perimeter estimation (cm) (Aa: Abies alba Mill.; Fs: Fagus sylvatica L.; Qp: Quercus petraea (Matt.) Liebl.; Pa: Picea abies (L.) H. Karst.).
Non-FisheyeFisheye
80130180All80130180All
Aa0.570.270.480.460.560.310.420.44
Fs0.860.750.760.790.390.250.570.42
Qp0.911.270.771.000.981.211.061.08
Pa1.031.171.101.100.550.750.760.69
Table 2. RMSE% of perimeter estimation (Aa: Abies alba; Fs: Fagus sylvatica; Qp: Quercus petraea; Pa: Picea abies).
Table 2. RMSE% of perimeter estimation (Aa: Abies alba; Fs: Fagus sylvatica; Qp: Quercus petraea; Pa: Picea abies).
Non-FisheyeFisheye
80130180All80130180All
Aa0.420.210.390.360.420.240.350.35
Fs0.990.900.940.950.450.300.710.51
Qp0.590.880.550.690.630.830.760.74
Pa0.590.740.720.680.310.470.500.43
Table 3. RMSE of diameter estimation (cm) (Aa: Abies alba; Fs: Fagus sylvatica; Qp: Quercus petraea; Pa: Picea abies).
Table 3. RMSE of diameter estimation (cm) (Aa: Abies alba; Fs: Fagus sylvatica; Qp: Quercus petraea; Pa: Picea abies).
Non-FisheyeFisheye
80130180All80130180All
Aa0.490.430.430.450.500.490.510.50
Fs0.510.480.430.470.380.390.430.40
Qp0.450.400.410.420.520.530.380.48
Pa0.670.490.460.550.500.450.540.50
Table 4. RMSE% of diameter estimation (Aa: Abies alba; Fs: Fagus sylvatica; Qp: Quercus petraea; Pa: Picea abies).
Table 4. RMSE% of diameter estimation (Aa: Abies alba; Fs: Fagus sylvatica; Qp: Quercus petraea; Pa: Picea abies).
Non-FisheyeFisheye
80130180All80130180All
Aa1.171.101.131.141.191.241.331.25
Fs1.851.831.701.801.391.511.671.52
Qp0.920.900.940.921.081.180.881.06
Pa1.210.980.961.070.920.901.130.98
Table 5. The average of reprojection error of both lenses that were used divided by tree species.
Table 5. The average of reprojection error of both lenses that were used divided by tree species.
Abies AlbaFagus SylvaticaQuercus PetraeaPicea Abies
Non-FisheyeFisheyeNon-FisheyeFisheyeNon-FisheyeFisheyeNon-FisheyeFisheye
Reprojection error (pixel)0.900.570.750.520.700.420.760.44

Share and Cite

MDPI and ACS Style

Mokroš, M.; Výbošťok, J.; Tomaštík, J.; Grznárová, A.; Valent, P.; Slavík, M.; Merganič, J. High Precision Individual Tree Diameter and Perimeter Estimation from Close-Range Photogrammetry. Forests 2018, 9, 696. https://doi.org/10.3390/f9110696

AMA Style

Mokroš M, Výbošťok J, Tomaštík J, Grznárová A, Valent P, Slavík M, Merganič J. High Precision Individual Tree Diameter and Perimeter Estimation from Close-Range Photogrammetry. Forests. 2018; 9(11):696. https://doi.org/10.3390/f9110696

Chicago/Turabian Style

Mokroš, Martin, Jozef Výbošťok, Julián Tomaštík, Alžbeta Grznárová, Peter Valent, Martin Slavík, and Ján Merganič. 2018. "High Precision Individual Tree Diameter and Perimeter Estimation from Close-Range Photogrammetry" Forests 9, no. 11: 696. https://doi.org/10.3390/f9110696

APA Style

Mokroš, M., Výbošťok, J., Tomaštík, J., Grznárová, A., Valent, P., Slavík, M., & Merganič, J. (2018). High Precision Individual Tree Diameter and Perimeter Estimation from Close-Range Photogrammetry. Forests, 9(11), 696. https://doi.org/10.3390/f9110696

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