CANOPY HEIGHT MODELS TO MONITOR AREAS UNDER DEGRADATION IN MIXED OMBROPHILE

This research was aimed at providing a study of digital elevation models from Airborne Laser Scanner (ALS) data. The goal was to improve the forest monitoring, or more specifically, arboreal canopy monitoring, which contains the most elevated vegetation of the forest. It is also described as the surface above the canopy top. The region of rain forest between the states of São Paulo and Paraná in Brazil was analyzed. The proposed method employs the ALS system that provides information related to different vegetative cover bedding, to analyze it and reach to estimate for the most elevated stage. The methodology is based on differentiating the identified surfaces in the original set of data. It uses these surfaces with a raster, which was obtained from an interpolation of different sets of points from the surfaces of interest. The raster structure, unlike vector data, allows for processes linked to the digital processing of images. It enables implementation of digital filters over the data. The method helped in reducing the influence of errors during classification and achieving the theoretical surface. The results revealed a coherent relationship with anthropogenic influence such as the introduction of exotic vegetation. We concluded that arboreal vegetation could be an index of conservation for environmental registration and licensing organizations.


INTRODUCTION
Arboreal canopy height has been described in the literature as an important parameter for characterizing the regeneration or degradation stage of a region.As such, it is highly relevant in studies and environmental impact reports required before engaging in activities that could be environmentally destructive.
In Brazil, studies such as Schaaf et al (2006) and Oliveira et al (2003) discuss the importance and special characteristics of tropical forest monitoring, particularly for Mixed and Dense Ombrophile Forests.
Dense Ombrophile Forests span the forest formations distributed on the intermediary portion of the slopes of Serra do Mar, specifically at altitudes greater than 400 m.Regions with descending slopes and mild temperatures, influenced by the flow of cold plateau air, exhibit flora that is completely different from those at lower levels.Nonetheless, montane and submontane ecosystems are physiognomically and structurally very similar.The Montane Dense Ombrophile Forest has an average height of 20 m.
Various species of the Lauraceae family are predominant in the upper stratum of this forest type, particularly Ocotea catharinensis, Ocotea odorifera, and Cryptocarya aschersoniana.Other species are also common, such as Copaifera trapezifolia, Aspidosperma olivaceum, Pouteria torta, and Lamanonia speciosa.Meanwhile, the lower strata are predominantly characterized by the Myrtaceae and Rubiaceae families.This formation is also characterized by the abundance of pteridophytes (ferns) present in the understory, particularly Dicksonia sellowiana and several species of the Cyatheacea family (SEMA, 2003).
This study aims to characterize the arboreal canopy height of a residual forest located in the surrounding area of the Governador Parigot de Souza Hydroelectric Power Plant in Paraná, Brazil, using data derived from Airborne Laser Scanner (ALS) and treated with image processing techniques.

Characteristics of the study area
The evaluated data correspond to the area surrounding the Governador Parigot de Souza Hydroelectric Power Plant dam located on the Capivari River.The site in question is in the municipality of Campina Grande do Sul, Paraná, Brazil alongside Highway BR-116, 50 km north of the southern end of the portion between Curitiba and São Paulo.
ALS data was collected by LACTEC -Institute of Technology for Development, based in Paraná.This center was created by COPEL, the state energy management company of Paraná, making the electrical sector one of its main research foci.
LACTEC scans were carried out using an ATLM 2050, manufactured by the Canadian company Optech.The system has a frequency of up to 50,000 Hz, records the first and last return pulses, has a variable scanning angle of +/-20°, and can operate at a maximum nominal altitude of 2,000 m.Table 1 shows the configurations used for data collection.
The results obtained showed an average error of 10 cm and a standard deviation of 50 cm in altitude.The larger error values were found in areas where the mesh of points generated by the topographic survey was not dense enough.During data collection, there was difficulty in determining what points corresponded to terrain in some areas, due to the large amount of litter and tripping vegetation.The ALS data present a series of constraints for their manipulation, seeing how the large volume of points presents two operating possibilities.
The first approach would consist of interpolating the data for processing in a regular grid or grayscale image.In this case, tools used in digital image processing  Canopy height models to monitor areas under... would be applied, since data would be aggregated in a raster grid.The downside of this approach is in the necessary interpolation and resampling, as this may heavily alter the original values.
An alternative approach is the use of coordinates as vectors, which would require changing the environment to a Computer Aid Design.In such a case, the work algorithms, necessary functions for the assortment, selection, division, and operation of the surfaces, increase in complexity.However, unlike the first method, data are not interpolated, i.e., they are less influenced, depending upon the available hardware capacity.There is an almost linear relationship between RAM and the number of processed laser points.

Method
Figure 1 shows the profile of a Dense Ombrophile Forest.Stratum A contains the crowns of dominant or emergent trees, which usually stand out above the continuous canopy of the forest.Just below, stratum B contains codominant trees, whose crowns touch and, as with the dominant trees, receive direct sunlight.
Eventually, stratum B might open up due to the loss of one or more trees, thereby allowing the penetration of light directly to stratum C.This stratum is composed of trees that do not typically receive sunlight otherwise, meaning that they are usually shorter.However, when a forest clearing opens, these trees quickly respond to the stimulus and tend to occupy the upper strata.These strata comprise the arboreal canopy.Under this canopy lies vegetation that constitutes the understory -namely, stratum D, which is composed of bushes, and stratum E, an herbaceous layer composed of prostrate weeds between 30 and 70 cm in height (POGGIANI, 1989).Samples with different average heights were taken from both steep and flat regions within the available area.
The most common technique for obtaining the digital terrain model (DSM) of the crowns of the trees is to calculate the triangulated irregular network (TIN) of the highest reflection (with the highest point of a preset neighborhood removed) and then interpolate the points using Delaunay triangulation.Canopy models are obtained with the subtraction of the digital terrain model (DTM) from the corresponding DSM.
In order to detect points in the forest, a sufficient density of laser pulse coverage in the region is necessary.This allows for the detection of return pulses that penetrated the vegetation through the spaces left between leaves and twigs.
After the selection, individual archives were separated for each sample and, as in the original archive, ASCII features were maintained.
Samples were separated into two classes of interest: a)Terrain Points: all points classified as terrestrial and usually derived from the last return mode.b) Vegetation Points: all points classified as vegetation (corresponding to the canopy, plants, bushes, buildings) and usually derived from the first return mode.
In this step, data was used in its vector form.As such, vegetation and terrain points were classified based on a visual assessment of sectional points, without the use of algorithms and automated segmentation techniques.
These two classes were established in order to produce two elevation models, the first corresponding to the terrain and the second corresponding to the canopy height.This separation of points resulted in the models shown in Figure 1a.In the following step, altimetric data for each of the models were correlated.It was not possible to do so between the surfaces without first obtaining the values corresponding to the altitude in the whole area of interest.Therefore, a regular grid model was used to compare.ALIXANDRINI JUNIOR, M.J. et al.
Subsequently, the subtraction of the terrain altitude from the altitude of the trees resulted in a third model corresponding to tree height, as shown in Figure 1b.
From this model, it was possible to obtain parameters such as average tree height, the number of trees of specific species, and crown volume.

RESULTS
A total of 30 equally sized rectangles, 50 x 50 m in area, were sampled in regions where the relief showed smooth and mountainous terrain, as shown in Figure 2.
The profiles were generated through the previously discussed classification of ALS points in regions of the profile line.The profiles were obtained through the vectorization of the union of those points.A large variation in altitude, directly associated with the proximity of the formation of Serra do Capivari, which is located east of the dam, can be observed.
Each sample contained between 13,000 and 30,000 ALS points derived from the archive containing the MDE, i.e., grouped points from the first and second laser pulses.From here, the areas were individually processed for the operations proposed in the methodology.
Table 2 shows the main characteristics of each area, specifically declivity, altitude, the total number of points, and total penetration of the vegetation.
The percent declivity values correspond to the ratio between the altimetric variation of the terrain and planimetric variation.Therefore, a declivity of 100% would correspond to a 45° angle between the average land surface and a horizontal surface.

DISCUSSION
The estimates correlated with different responses that could be directly associated with the regeneration of vegetation in the sampled areas.This was apparent when the responses were located in the area surrounding the reservoir.
The higher values obtained are compatible with the expected values for a Montane Mixed Ombrophile Forest.However, it is possible that the lower values  Canopy height models to monitor areas under... that do not correspond indicate the insertion of an exotic species into the region.In this way, vegetation height estimation could act as strong indication of environmental change.
Associating this indicating factor with a more careful description of the testing environment can be of great value for both environmental records and environmental licensing, as it is a much more detailed survey compared to the one currently used for both cases.
This research showed one potential application of ALS, specifically in acquiring altimetric data.Such an application could highlight relationships with surrounding properties, as well as help in evaluating and monitoring areas in the hydrographic basin of the reservoir.

CONCLUSIONS
Exotic species used in traditional Brazilian silviculture, such as Pinus elliot, Pinus taeda, Eucalyptus dunni, and Eucalyptus bethami, can at least have identified.This is because the pinus and eucalyptus species, which are mostly cultivated in the south and southeast, have an adult height of up to 40 and 50 m, respectively.This is in comparison to the average height of the Montane Mixed Ombrophile Forest, which is between 20 and 25 m.This opens up the possibility of monitoring the migration movement of these species through the observation of exotic adults, mainly in preservation parks of the original biome.A specific evaluation of the influence of litter and tripping vegetation on the final accuracy of ALS data collection is necessary.Automated classifications that give a better response to tropical vegetation in soil recognition for MDT derivation are needed.The major problem in this case is that most of the algorithms developed internationally are based on problems generated in temperate climates; for example, vegetation and the corresponding canopy density are typically assumed to vary by season.Testing this approach by obtaining corroborating data in homogenous regions is recommended.This will also prove useful in verifying differences in estimates for coverage on different system configurations.
The resulting reliance on winter and summer data does not accurately reflect the problems caused by various types of flora in Brazil.Applying TIN techniques to this methodology, instead of using regular grids, is also recommended.This will decrease iteration with the original data, opening up the possibility of using it for comparison, without samples, in different areas with similar vegetation.

ACKNOWLEDGMENTS
Thanks to the LACTEC team, represented by the engineer and cartographer Marlo Martins, for support and access to ALS data.

Figure 1 -
Figure 1 -a) Digital Terrain Model and Digital Surface Model and b) Vegetation Height Model.Figura 1 -a) Modelo digital de terreno e modelo digital de superfície b) Modelo da altura da Vegetação.