Aerial photo-interpretation using Z/I DMC images for estimation of forest variables

Abstract

Interpretation and tree height measurements in aerial photographs using photogrammetric workstations are frequently performed in standwise forest inventory. Images acquired by digital aerial cameras are now replacing the traditional film-based aerial photographs. In this study, digital images from the airborne Z/I DMC system for standwise estimation of stem volume, tree height and tree species composition were investigated at a 1200 ha forest area located in southern Sweden (58°30′N, 13°40′E). The 56 selected stands were dominated by Norway spruce [Picea abies (L.) Karst.] and Scots pine (Pinus sylvestris L.) with stem volume in the range of 30–630 m³ ha^?1 (average 300 m³ ha^?1) and tree height in the range of 6–28 m (average 20 m). The large-format pansharpened colour infrared images were captured at a flight altitude of 4800 m above ground level corresponding to a pixel size of 0.48 m. The photo-interpretation was conducted by four professional interpreters, independently. In particular, two different base-to-height ratios (i.e. the ratio between the ground distance between image centres at the time of exposure and the flight altitude above ground level) of 0.26 and 0.39 were evaluated, but no significant difference in the estimation accuracy for stem volume and tree height was found. The accuracy for stem volume estimation in terms of relative root mean square error, corrected for systematic errors, was on average 24% (in the range of 17–39%). The corresponding accuracy for tree height estimation was on average 1.4 m (in the range of 0.9–1.6 m). The tree species composition accuracy assessment using a fuzzy set evaluation procedure showed that 95% of the stands were correctly classified. The estimation accuracies are in agreement with previous results using conventional film-based aerial panchromatic photographs.     

Tree Species Classification using Semi-automatic Delineation of Trees on Aerial Images

The purpose of this study was to develop a method for classifying tree species from remote sensing images by combining a semi-automatic pattern recognition technique and spectral properties of trees. Five stands in southern Finland were studied. Individual trees in the digital colour infrared (CIR) aerial photographs were segmented by a method based on the recognition of tree crown patterns at subpixel accuracy. The images were filtered with the Gaussian N-by-N smoothing operator and local maxima above a threshold level were segmented. The segments were classified into three tree species classes. The kappa coefficients for stands varied from 0.43 to 0.86 when the training data and test data were from the same aerial photograph. When training data from other photographs were used as reference data, the kappa coefficients ranged from 0.40 to 0.75. The method described provides an interesting approach for detecting tree species semi-automatically in digital aerial data.     

Deciduous forest mapping using change detection of multi-temporal canopy height models from aerial images acquired at leaf-on and leaf-off conditions

Discrimination of deciduous trees using spectral information from aerial images has only been partly successfully due to the complexity of the reflectance at different view angles, times of acquisition, phenology of the trees and inter-tree radiance. Therefore, the objective was to evaluate the accuracy of estimating the proportion of deciduous stem volume (P) utilizing change detection between canopy height models (CHMs) generated by digital photogrammetry from leaf-on and leaf-off aerial images instead of using spectral information. The study was conducted at a hemi-boreal study area in Sweden. Using aerial images from three seasons, CHMs with a resolution of approximately 0.5?m were generated using semi-global matching. For training plots, metrics describing the change between leaf-on and leaf-off conditions were calculated and used to model the continuous variable P, using the Random Forest approach. Validated at sub-stands, the estimation accuracy of P in terms of root mean square error and bias was found to be 18% and ?6%, respectively. The overall classification accuracy, using four equally wide classes, was 83% with a kappa value of 0.68. The validation plots in classes of high proportion of coniferous or deciduous stem volume were well classified, whereas the mixed forest classes showed lower classification accuracies.     

Stem Form and Taper Changes After Thinning and Nitrogen Fertilization in Picea abies and Pinus sylvestris Stands

Stem form and taper changes after thinning and thinning combined with N fertilization were studied in 23 Norway spruce [Picea abies (L.) Karst.] and 46 Scots pine (Pinus sylvestris L.) stands in Sweden observed for 5-28 yrs. Average taper was calculated for the stem section 1.3-6.0 m above ground level and as the ratio between the diameter at breast height and total tree height for periods of 5-10 yrs. At the start of the experiment, before the first thinning, the dominant height was 12-15 m. Once 65% of the basal area had been removed in a single thinning from below, the remaining trees showed a strong increase in taper compared with trees in unthinned stands. Stems in stands treated with heavily recurrent thinnings from below also increased in taper, although the increases were not always statistically significant. Trees in thinned, N-fertilized Scots pine stands in middle and northern Sweden developed a more pronounced taper compared with stems in equally thinned, unfertilized stands.     

Patterns of Stem Damage by Moose (Alces alces) in Young Pinus sylvestris Stands in Sweden

The incidence of recent moose browsing of the main stem on young Scots pine (Pinus sylvestris L.) trees, and correlations with tree or stand characteristics, were examined using data from a large-scale survey in Sweden. On average, 10% of all Scots pine main stems showed one or more types of recent damage, with browsing of the apical leader accounting for about 75% of all damage recorded. Stripping of bark off the main stem occurred most frequently in the southern area, where site productivity is highest. Apical leader browsing was negatively correlated with height of the main stems, length of the apical leader and incidence of pre-commercial thinning, whereas bark stripping and stem breaking were positively correlated with the same variables. All types of recent stem damage correlated positively with the extent of previous stem damage and also with the extent of recent browsing of lateral shoots.     

Relationships of stem surface area to other stem dimensions for Japanese cedar (<Emphasis Type="Italic">Cryptomeria japonica</Emphasis> D. Don) and Japanese cypress (<Emphasis Type="Italic">Chamaecyparis obtusa</Emphasis> Endl.) trees

I investigated the relationships between stem surface area outside bark and other stem dimensions for Japanese cedar (Cryptomeria japonica D. Don) and Japanese cypress (Chamaecyparis obtusa Endl.) trees. The stem dimensions used here were the basal area and the product of diameter at breast height (dbh) and total tree height. The regression equation of the stem surface area s against basal area g was s = 184.216g for the cedar trees and s = 156.878g for the cypress trees. The slope of the equation was significantly different between the two species. For the same dbh, the cedar trees tended to have a higher total tree height than the cypress trees. The cedar trees also had a larger surface area of relative stem form than the cypress trees. These results indicated that the difference in the slope was produced by the differences in both the stem slenderness and tapering between the two species. On the other hand, the regression equation between the stem surface area and the product of dbh and total tree height dh was s = 1.937dh for the cedar trees and s = 1.921dh for the cypress trees, whereas no significant difference in the slope was found. The obtained slopes for the cedar and cypress trees seemed to be in accord with that for other coniferous species reported in earlier studies, suggesting that the variation in the slope among coniferous species would be small. The estimation from the basal area would provide a simpler means for estimating the stem surface area and would be useful in obtaining an approximation of the surface area. By contrast, the estimation from the product of dbh and total tree height would provide a more accurate and precise estimate as well as a wider applicable range, i.e., a parameter for physiological growth models. In conclusion, it could not be judged which regression equation examined in the present study was superior to the other, and thus it was important to select an appropriate equation depending both on the purpose and on the time and labor available.     

Mixed-effect non-linear modelling for diameter estimation along the stem of Tectona grandis in mid-western Brazil

This study evaluated the efficiency of taper functions and the application of mixed-effect modelling for diameter estimation along the stems of Tectona grandis. We sampled 266 trees of Tectona grandis, measuring the diameter at relative heights for volume determination, grouping the data according to three form-factor classes. Six taper functions were fitted, selecting the function with better fit performance. Six taper functions were fitted, selecting the function with better fit performance. The selected function was fitted in its basic formulation, and with the mixed non-linear modelling technique in different scenarios, and for the stem stratified in three portions of the total height. The precision and selection of the adjusted models were evaluated regarding the coefficient of determination, standard error of estimate, the Akaike information criterion, bias, quadratic error and absolute bias. According to the statistical criteria used, the model of Kozak was selected for the adjustments. For diameter estimation, the scenario with two coefficients as random effects provided an accuracy increase of 11.91%, and the mixed non-linear modelling better estimated the stem diameter for the stratified stems. In conclusion, the model of Kozak can be used to describe the stem shape of Tectona grandis, and the mixed-effect non-linear model approach was the best technique to estimate diameter along the stem of Tectona grandis.     

Forest variable estimation using photogrammetric matching of digital aerial images in combination with a high-resolution DEM

Abstract

The rapid development in aerial digital cameras in combination with the increased availability of high-resolution Digital Elevation Models (DEMs) provides a renaissance for photogrammetry in forest management planning. Tree height, stem volume, and basal area were estimated for forest stands using canopy height, density, and texture metrics derived from photogrammetric matching of digital aerial images and a high-resolution DEM. The study was conducted at a coniferous hemi-boreal site in southern Sweden. Three different data-sets of digital aerial images were used to test the effects of flight altitude and stereo overlap on an area-based estimation of forest variables. Metrics were calculated for 344 field plots (10 m radius) from point cloud data and used in regression analysis. Stand level accuracy was evaluated using leave-one-out cross validation of 24 stands. For these stands the tree height ranged from 4.8 to 26.9 m (17.8 m mean), stem volume 13.3 to 455 m³ ha^?1 (250 m³ ha^?1 mean), and basal area from 4.1 to 42.9 m² ha^?1 (27.1 m² ha^?1 mean) with mean stand size of 2.8 ha. The results showed small differences in estimation accuracy of forest variables between the data-sets. The data-set of digital aerial images corresponding to the standard acquisition of the Swedish National Land Survey (Lantmäteriet), showed Root Mean Square Errors (in percent of the surveyed stand mean) of 8.8% for tree height, 13.1% for stem volume and 14.9% for basal area. The results imply that photogrammetric matching of digital aerial images has significant potential for operational use in forestry.     

Automatic Stem Mapping in Three Dimensions by Template Matching from Aerial Photographs

Recent developments in remote sensing have shown promising results for the detection and species classification of individual trees from high-resolution imagery of forests in leaf-on condition. Existing image-analysis algorithms produce a two-dimensional stem map in the observed imagery. This paper develops a new algorithm that generalizes existing template-matching algorithms to produce a stem map in the three dimensions of a forest in leaf-off condition. A self-contained algorithm to produce artificial template images is introduced. The image-analysis algorithm was applied to oak ( Quercus robur L.) grown under standard silvicultural treatment in Denmark. One analysis gave detection rates in the range 85-98% and root mean squared errors on the stem base coordinates below 8.1, 7.6 and 20.5 cm in the x -, y - and z -directions, respectively. It is concluded that it is indeed feasible to produce a quite accurate three-dimensional stem map.     

Stem form variations in the natural stands of major commercial softwoods in eastern Canada

Stem form is often used as a sawlog assortment criterion and has an important effect on lumber recovery and mechanical properties. Based on 7018 stems collected from the natural stands in 3 regions in eastern Canada, this study quantified the variations in stem form (taper, sweep and eccentricity) for five major commercial softwood species, viz., jack pine (Pinus banksiana Lamb.), black spruce (Picea mariana), white spruce (Picea glauca), red spruce (Picea rubens) and balsam fir (Abies balsamea). The majority of the stems had a diameter at breast height (DBH) ranging from 10 to 32 cm and total tree height from 10 to 20 m. Stem taper and butt taper increased steadily from 0.38 to 1.75 cm/m and from 0.74 to 4.23 cm/m, respectively, with DBH increasing from 8 to 48 cm. Stem sweep and basal sweep ranged from 0.12 to 7.28 cm/m and from 0.12 to 9.85 cm/m, respectively, with mean values of 0.87 and 0.95 cm/m. There were no significant differences in stem sweep and basal sweep between DBH classes and both sweeps distributed over a wide range within DBH classes, species and regions. About 0.77% of the stems and 1.2% of the butt logs had seriously curved shapes (sweep >3.0 cm/m). Eccentricity (ratio of smaller diameter to larger diameter) at the breast height (0.96) was very close to 1. Balsam fir and jack pine had the best and worst stem forms, respectively, while stem forms in black spruce were in-between. Butt logs were most seriously curved in red spruce and most tapered in white spruce. Region C (mainly eastern Quebec and the Atlantic Provinces) produced stems with smaller butt taper, basal sweep and stem sweep and larger eccentricity compared with other regions in eastern Canada. Overall stems in eastern Canada have a good stem form in terms of taper, sweep and eccentricity. The evaluation of stem form provides valuable timber quality information for the wood industry to optimize wood processing as well as for forest managers to define appropriate silvicultural regimes to improve stem form and economic value.     

Determinants of stem form: Application to Tectona grandis (Linn. F) stands

ABSTRACT

Stem form information can enhance rapid estimation of stem volume. Studies on important determinants of stem form in pure stands are scanty. This study aimed at identifying determinants of stem form using Tectona grandis stands in Omo Forest Reserve, Nigeria. Twenty 0.04 ha square plots were randomly selected across four stand ages. Crown diameter, total and merchantable heights (MH), diameters at breast height (DBH), at 5.27 m tree height and at midpoint from breast height to tip of tree were measured in each plot. Measured trees were classified into four canopy layers (CL)—dominant, codominant, intermediate, and suppressed. Merchantable volume, basal area (BA), form factor (F_f), normal form quotient (Q_n), Girard form class (Q_G), absolute form quotient (Q_a), ratio of diameter at merchantable limit to diameter at the base (D_t/D_b), and tree slenderness coefficient were derived from the data. Analysis of variance, correlation, and regression analyses were used to analyze the data. Significant differences were observed in stem form and tree size variables under different canopy layers. The DBH, BA, MHT, and CL were identified as important variables associated with stem form. Inclusion of Q_a improved most the predictive ability of the selected stem volume equation.     

Possibilities for Harvester-based Forest Inventory in Thinnings

This study investigated the possible estimation of forest characteristics using the information collected by the harvester in first thinnings. For the analysis a complete forest inventory was carried out in a stand, which was subsequently thinned. The global mean values of tree diameter, tree height, basal area and stem density were estimated, and further, a spatial analysis was carried out to investigate whether the data collected by the harvester could be used to generate a continuous spatial model of the forest. The results indicated that the global mean diameter and height may be estimated, whereas area-related properties, such as basal area and stem density, are more difficult to estimate. The spatial distribution of the diameter and height remained similar after the thinning, whereas the basal area and stem density had become more homogeneous after the thinning. From the trees removed in the thinning a continuous spatial model of tree diameter was developed. It reproduced the spatial structure of the original trees to some extent ( R 2 = 0.27, RMSE = 14.3 mm).     

Models for Predicting Dry Matter Content of Norway spruce

Models were developed for predicting the dry matter content of single Norway spruce [Picea abies (L.) Karst.] stems from non-destructive field measurement data such as tree height, diameter and age. The material under study comprised 1612 stem discs originating from 235 trees grown at nine locations in Denmark and southern Sweden. In total, 153 trees (1054 discs) originating from a thinning experiment in southern Denmark were used for establishing the models, which were later validated on independent material. Results of the validation were that the 95% confidence interval for single stem density was ±9 to ±16% when based on average ring width at breast height, but within ±7% when based on wood samples taken at breast height. The density models were combined with standard volume models in order to predict dry matter content. Validation on independent material showed that in general the dry matter content of single stems can be estimated with less than 10% error (RMSE).     

A stem taper equation compatible to volume equation for <Emphasis Type="Italic">Paraserianthes falcataria </Emphasis>in Pare,East Java,Indonesia: its implications for the plantation management

A stem taper equation compatible to the volume equation was derived from the data of 172 trees sampled in the three plantations of Paraserianthes falcataria in Pare, East Java. The stem taper equation in this study was based on a conventional logarithmic volume equation, with the addition of one free parameter that minimized the standard errors of estimate to the observed diameter of stems at 1.0-m intervals. The derived equation was able to predict the stem taper of P. falcataria reasonably well, while it was not flexible enough to describe delicate changes of the taper in the upper portion of stem. The equation was able to predict bole length and the volume up to the top diameter 20 cm, the diameter size merchantable for sawn timber in Indonesia, with reasonable accuracy. An appropriate rotation age for sawn-timber production was found to be 8–9 years old, at which the mean annual increment of sawn timber would be maximized.     

Xylem cavitation in roots and stems of Douglas-fir and white fir

Roots of hardwoods have been shown to be more vulnerable to xylem cavitation than stems. This study examined whether this pattern is also observed in a conifer species. Vulnerability to cavitation was determined from the pressure required to inject air into the vascular system of hydrated roots and stems, and reduce hydraulic conductance of the xylem. According to the air-seeding hypothesis for the cavitation mechanism, these air pressures predict the negative xylem pressure causing cavitation in dehydrating stems. This was evaluated for stems of Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) and white fir (Abies concolor (Gord. & Glend.) Lindl.). The air-injection method was applied to roots and stems of different sizes and positions in Douglas-fir trees. Roots, especially smaller roots with a xylem diameter < 5 mm, were more vulnerable to cavitation than stems. Mean cavitation pressure for smaller roots was -2.09 +/- 0.42 versus -3.80 +/- 0.19 MPa for larger roots (> 8 mm diameter). Within the shoot system, smaller stems (< 5 mm diameter) were most vulnerable to cavitation, having a mean cavitation pressure of -4.23 +/- 0.565 versus -5.27 +/- 0.513 MPa for large stems (> 8 mm diameter). There was no correlation between tracheid diameter and mean cavitation pressure within root or stem systems, despite larger tracheid diameters in roots (23.3 +/- 3.9 micro m) than in stems (9.2 +/- 1.6 micro m). Smaller safety margins from cavitation in roots may be beneficial in limiting water use during mild drought, and in protecting the stem from low xylem pressures during extreme drought.     

Early growth and quality of radiata pine in a silvopastoral system in New Zealand

The effect of different understoreys on the early growth and quality of radiata pine (Pinus radiata) in a silvopastoral system was evaluated in Canterbury, New Zealand. Fitted models successfully represented the effect of understorey on the growth (basal area and height) and a quality of trees (diameter over stubs). Basal area, diameter, height and volume growth were affected by understorey competition. Trees growing with no-understorey had 34%and 29% higher volume than trees growing with lucerne (Medicagosativa) and cocksfoot (Dactylis glomerata)understoreys, respectively. However, the height of trees at age 10 grown without understorey competition was only 9% and 15% more than trees grown with grass and lucerne understorey, respectively. Understorey competition effects on growth became more evident during the summer and early autumn months when water deficits occurred. The quality of trees was also affected by understorey competition. While trees grown in bare ground had higher growth than trees grown with understorey competition, they had the poorest stem quality with higher diameters over stubs, branch size, sweep and proportion of stem defects. This will have an impact on the agroforestry system profitability by affecting the net harvest revenue. This revised version was published online in June 2006 with corrections to the Cover Date.     

基于大比例尺航片的针叶树种冠幅的提取

基于凉水国家级自然保护区2009年拍摄的1：2000航空像片和同期森林资源二类调查的固定样地数据,采用子像元分类方法分别提取出红松、落叶松和云冷杉的专题影像图。在此基础上,将栅格专题影像图转换为矢量图形,采用目视解译的方法提取上层针叶林的树冠信息。通过将针叶树冠形似为圆形提取出各树种的冠幅,用固定样地实测数据进行对比分析和精度评价,并建立航片上提取冠幅与实测冠幅之间的一元线性回归模型。结果表明：红松、落叶松和云冷杉冠幅的提取精度分别达到83．50％、84．35％和82．26％,其预测精度分别达到83．60％、81．46％和83．57％。     

Basal sweep and compression wood in young Scots pine trees

This study demonstrates the correlation between stem form and compression wood content in 36 sampled trees from a 6-yr-old container grown Scots pine ( Pinus sylvestris L.) plantation in central Sweden. Root analyses were carried out to examine causes for the high incidences of basal sweep. On the sampled trees basal sweep, stem sweep and maximum bow-height were measured. Discs were cut at different heights and the compression wood content was analysed using digital image analysis. Pith eccentricity and out-of-roundness were calculated for each disc as well as per stem. Compression wood and pith eccentricity was most pronounced near the stem base, where sweep tends to be high. In general, pith eccentricity increased with degree of basal sweep. As a conclusion, expression of basal sweep can be used to predict compression wood content in young Scots pine trees. This study also shows the importance of anchoring of a tree to prevent basal sweep formation due to instability.     

Functions for Biomass Estimation of Young Pinus sylvestris,Picea abies and Betula spp. from Stands in Northern Sweden with High Stand Densities

New silvicultural regimes with high within-stand competition require new functions for estimation of standing stock and growth of biomass components, since the allometry of trees is changed by light competition. This paper presents functions for estimation of the aboveground biomass dry weights for stem wood, stem bark, branches and leaves of young (diameter at breast height <10 cm) Scots pine (Pinus sylvestris L.), Norway spruce [Picea abies (L.) Karst.] and birch (Betula pendula Roth. and Betula pubescens Ehrh.) trees growing in dense mixed stands. The functions were derived from a sample consisting of 84 Scots pine, 43 Norway spruce and 66 birch trees from six stands in northern Sweden with high stand densities (>10000 st ha^-1). The logarithmically transformed power function displayed a good ability to stabilize the variance of dry weights and showed a good fit to the material (0.37< R ² <0.99). A comparison with the most commonly used biomass functions in Sweden today showed that they overestimated the weight of stem wood and branches, while the weight of foliage was underestimated. The nature of these discrepancies suggested that the precision of biomass estimations might also be improved for young trees at wider spacing.     

A forest biomass yield table based on an empirical model

We report an empirical model for estimating unutilized wood biomass, and its application to Cryptomeria japonica D. Don and Larix kaempferi in Tohno City, Iwate Prefecture, northeast Japan. Outputs from the model are the quantity of unutilized wood biomass and merchantable volume produced by timber harvest. The unutilized wood biomass is divided into stumps, tops, branches, foliages, small trees, and unutilized stems due to their defects. Inputs to the model are mean diameter at breast height (DBH), mean tree height, trees per unit area, and timber utilization standards. DBH distribution, DBH–height curve, stem form, bark thickness, and relationship of stem biomass to foliage and branch biomass could be described by the proposed model, indicating its validity. The proposed model enables us to develop the forest biomass yield tables modified from the existing stem volume yield tables. The developed forest biomass yield tables indicated that the unutilized wood biomass due to defects accounted for the largest part of the whole unutilized wood biomass, and that the ratio of unutilized parts in stem volume to total stem volume could vary with stand age and site productivity class. Based on a comparison of the developed forest biomass yield tables with those reported previously, we concluded that the proposed model-based forest biomass yield table would be useful for estimating the quantity of unutilized wood biomass.