Document Type : Research Paper

Authors

1 Postdoctoral Researcher in Forestry, Faculty of Natural Resources, Sari University of Agricultural Sciences, Mazandaran, Sari, Iran

2 Professor, Department of Forestry, Faculty of Natural Resources, Sari University of Agricultural Sciences and Natural Resources, Sari, Iran

Abstract

  Extended Abstract
 Introduction
 Estimation of forest habitat characteristics is a necessary issue in order to collect information for sustainable forest management (Ahmadi et al., 2020). Data collection methods require a lot of time and money. Therefore, it is always tried to use complementary methods, with lower costs and acceptable accuracy, using the achievements obtained in various scientific fields (Sivanpillai et al., 2006). Sentinel 2 is a new generation optical satellite for Earth monitoring developed by the European Space Agency with new spectral capabilities, wide coverage and good spatial and temporal resolution for data continuity and enhanced Landsat and Spot missions (Wang et al., 2017). When the size of the population is not very large, the application of each of the simple random, classification and systematic methods leads to a more or less similar result. But when the size of the community increases, these methods are associated with problems such as: preparing a sampling framework, high cost of surveying sample units with high dispersion and preparing a sampling plan from units far from each other (Zubair, 2007). The cluster method is one of the recommended methods for large areas in which instead of one sample plot, several sample plots are harvested in one part of the study area (Yim et al., 2015). Among the researches done on the mentioned subjects are the research of Kleinn (1994), Ismaili et al. (1396), Behera et al. (2021), Sibanda et al. (2021), Praticò et al. (2021), Nazariani et al. (1400) and Dabija et al. (2021). Although studies on estimating quantitative forest characteristics using distance measurement data and nonparametric algorithms in Zagros forests may have been done extensively, the effect of main and artificial bands to estimate canopy characteristics and density (number Per hectare) using Sentinel 2 images in the forests of Watershed Orfi Olad Ghobad Koohdasht with the aim of selecting the optimal cluster design to save time and money to achieve forest inventory has not been reported, so in this study, we tried to investigate this issue.

Materials and methods

 In order to conduct the present study, a part of the Zagros forests located 35 km north of Koohdasht city, named Watershed Olad Ghobad was selected. Sampling points were determined in a regular-random manner using a grid with dimensions of 600 × 500 meters. Then, at each sampling point, 16 different cluster sampling designs with four circular and square subplots were designed and implemented. The radius of the circular subplots was 15 meters, the diameter of the square sample was 37 meters and the distance between the subplots was 60 meters. Then, the information on the characteristics of the number per hectare and canopy of trees including the number, of two large and small canopy diameters per sample was measured. In this study, Sentinel 2 sensor images related to August 6, 2021, equivalent to summer 1400, were used at the L1C correction level. This level of correction is geometrically error-free due to the reference ground and because their reflection is at the upper level of the atmosphere. In the present study, four bands (2-blue band, 3-green band, 4-red band, and 8-near-infrared band) of this sensor with a resolution of 10 meters were used. In general, Sentinel 2 image preprocessing operations involve radiometric and geometric correction. The image processing also includes various operations such as grading, texture analysis, band integration, and fabrication of plant features (Naghavi, 2014). In addition to the main bands, artificial bands were created by applying appropriate processing, which was used in the modeling process. Spectral values ​​equivalent to ground plots were extracted from the main and artificial bands and used as an independent variable in the models. In order to evaluate and fit the regression models, 25% of the data were randomly selected (Lu et al, 2004) and excluded from the evaluation data set. The validity of statistical models was evaluated using the coefficient of determination of the mean squared error squared, bias, mean squared error, and squared percentage. In total, ArcGIS software was used to implement the sample parts on the image, ENVI software was used for image processing and STATISTICA software was used for modeling.
Results
In this method, during data validation, the results showed the characteristic of number per hectare of cluster 16 and the characteristic of canopy cover of cluster 15 with a coefficient of explanation (0.66) and (0.59), respectively, it has the highest accuracy. The results obtained from the application of the nearest neighbor algorithm with four criteria of Euclidean distance, Euclidean square, Manhattan, and Chapichev showed that for the number of characteristics per hectare, the Euclidean distance criterion with cluster 16 and for the canopy characteristic of the Euclidean distance criterion with cluster three, respectively (R2 = 0.59 and RMSE=5.70%) and (R2 = 0.62 and RMSE= 12.30%). The accuracy and efficiency of the support vector machine algorithm are influenced by the type of kernel used. The results of different kernels by considering different cluster sampling designs in the backup vector machine method showed for the characteristic number of linear kernel trees and 13 cluster sampling designs with an explanation coefficient of 0.72 and for the canopy characteristic. The linear kernel and the cluster sampling design of seven with a coefficient of determination of 0.65 have the best results. Evaluation of the artificial neural network model showed that the MLP algorithm is more suitable than the RBF algorithm in estimating the studied characteristics with its high accuracy and average squared percentage. Based on this, among the 16 designs used with the MLP algorithm, they showed the most suitable results for the number of characteristics per hectare of cluster six with a coefficient of reflection of 0.86 and for the canopy characteristic of cluster 10 with a coefficient of reflection of 0.76, respectively. Based on the values ​​of the coefficient of explanation and the lowest squared percentage of the mean squares of error, the most appropriate model was selected from the four types of algorithms studied in modeling and the results showed both characteristics of the artificial neural network model respectively (with MLP algorithms MLP 80-20-1 and MLP 80-11-1) presented optimal results with explanation coefficients of 0.86 and 0.76.
Discussion and conclusion
The modeling results with four studied algorithms for the canopy characteristic showed that the artificial neural network model algorithm with a cluster sampling design of 10 with an explanation coefficient of 0.76 was the most suitable method. The results are consistent with the study (Yim et al., 2015;) and show the superiority of using cluster sampling, nonparametric modeling of the artificial neural networks and Sentinel 2 images in the structure of the forest ecosystem. Yim et al. (2015) acknowledged that in natural environments, the correlation between sub-plots and habitat conditions in terms of their shape and size should be more sensitive to forest structure. According to the study of Sivanpillai et al. (2006) in poorer masses, due to the presence of more gaps in the canopy, absorption and distribution occur. In contrast, Dabija et al. (2021) compared support vector machine and stochastic forest algorithms for canopy mapping using Sentinel-2 and Landsat 8 satellite imagery to evaluate regional and spatial classification and development in three different regions. Catalonia, Poland, and Romania paid. The results showed that Sentinel-2 satellite images were better than Landsat 8 data inaccuracy (8-10%) in land cover classification and radial-based support vector algorithm than in random forest with accuracy (6-7%). Function. Nazariani et al. (1400) also had the stochastic forest algorithm as the most suitable model for estimating the canopy characteristic, which is not consistent with the results of the present study. The reason for the difference can be found in the type of algorithm obtained and the accuracy achieved.

Keywords

Main Subjects

  1. Ahmadi, K., Kalantar, B., Saeidi, V., Harandi, E. K., Janizadeh, S., & Ueda, N. (2020). Comparison of machine learning methods for mapping the stand characteristics of temperate forests using multi-spectral sentinel-2 data. Remote Sensing, 12(18), 3019.
  2. Behera, M. D., Barnwal, S., Paramanik, S., Das, P., Bhattyacharya, B. K., Jagadish, B. & Behera, S. K. (2021). Species-Level Classification and Mapping of a Mangrove Forest Using Random Forest—Utilisation of AVIRIS-NG and Sentinel Data. Remote Sensing, 13(11), 2027.
  3. Bell, E., Howell, B.M., Johnson, G.V., Solie, J.B., Raun, W.R., & Stone, M.L. (2004). A comparison of measurements obtained using optical sensing with turf growth, chlorophyll content, and tissue nitrogen. Horticultural Science, 39(5): 1130-1132.
  4. Breiman, L. (2001). Random forests. Machine learning, 45(1), 5-32.
  5. Castleman, K. R. (1996). Digital image processing. Prentice Hall Press.
  6. Chen, J. M. (1996). Evaluation of vegetation indices and a modified simple ratio for boreal applications. Canadian Journal of Remote Sensing, 22(3), 229-242.
  7. Dabija, A., Kluczek, M., Zagajewski, B., Raczko, E., Kycko, M., Al-Sulttani, A. H., & Corbera, J. (2021). Comparison of support vector machines and random forests for corine land cover mapping. Remote Sensing, 13(4), 777.
  8. Egbers, R. (2016). Sentinel-2 data processing and identifying glacial features in Sentinel-2 imagery. TU Delft, University of Technology, The Netherlands.
  9. Esmaeeli, Z., Pilehvar, B., Kaboodi, A., & Mirazadi, Z. (2017). The Appropriate Sampling Method for Estimating Density and Crown Canopy of Declined Oak Stands in Dinarkooh Protected Forest, Abdanan, Ilam. Ecology of Iranian Forest, 5 (10) :53-60.
  10. Fatolahi, M., Fallah, A., Hojati, S. M., & Kalbi, S. (2013). Determining the capability of SPOT-HRG measuring data in estimating the number of trees per hectare. Journal of Wood and Forest Science and Technology, 20(4): 117-133.
  11. Hamidi, S. K., Weiskittel, A., Bayat, M., & Fallah, A. (2021). Development of individual tree growth and yield model across multiple contrasting species using nonparametric and parametric methods in the Hyrcanian forests of northern Iran. European Journal of Forest Research, 140(2), 421-434.
  12. Holmström, H., & Fransson, J. E. (2003). Combining remotely sensed optical and radar data in k NN-estimation of forest variables. Forest Science, 49(3), 409-418.
  13. Khorami, R. (2004). Investigating the capability of Landsat7ETM+ satellite data in estimating the standing volume of beech stands (case study: Sengdeh forests). Master's thesis. Forestry Department. Department of natural resources. University of Tehran.
  14. Kleinn, C. (1994). Comparison of the performance of line sampling to other forms of cluster sampling. Forest ecology and management, 68(2-3), 365-373.
  15. Kumar, M., Singh, M. P., Singh, H., Dhakate, P. M., & Ravindranath, N. H. (2020). Forest working plan for the sustainable management of forest and biodiversity in India. Journal of Sustainable Forestry, 39(1), 1-22.
  16. Lu, D., Mausel, P., Brondızio, E., & Moran, E. (2004). Relationships between forest stand parameters and Landsat TM spectral responses in the Brazilian Amazon Basin. Forest ecology and management, 198(1-3), 149-167.
  17. Masihpoor, M., Darvishsefat, A.A., Rahmani, R., & Fatehi, P. (2021). Phenological parameters trend of the southern Zagros forests based on MODIS-NDVI time series during 2000-2017. Iranian Journal of Forest, 12 (4): 577-590.
  18. Mohammadi, J., Shataei, Sh., & Namirania, M. (2014). Comparison of quantitative and qualitative characteristics of forestsstructure and composition in natural and managed forest stands
    (Case study: Shast Kalate forests of Gorga). Journal of Wood & Forest Science and Technolog, 21 (1): 65-83.
  19. Naghavi, H. (2014). Application of Quickbird satellite images in estimating the canopy level of Zagros forests (case study: Gol Khorramabad Castle area). Dissertation of the doctoral course in Forestry, Faculty of Natural Resources, Sari University of Agricultural Sciences and Natural Resources.
  20. , Fallah, A., Ramezani, H., Naghavi, H., & Jalilvand, H. (2021). The Effect of Different Cluster Sampling Schemes in Estimating the Quantitative Characteristics of Zagros Forests Using Sentinel 2 Sensor Images. Remote Sensing and GIS Scientific-Research Journal of Iran, Published online on November 15, 2022. 10.52547/gisj.2021.222017.1024.
  21. Norian, N., Shatai Joibari, Sh., Mohammadi, J., & Yazdani, S. (2014). Estimation of forest quantitative characteristics using ASTER sensor data and CART algorithm (Study area: Shast Kolate forest of Gorgan). Iranian Journal of Forest and Poplar Research, 22(3): 434-446.
  22. Norodini, A. & Porshkuri, F. (2011). Forest canopy classification on aerial photographs using textural analysis (case study: Taf Lorestan forest). Iran Remote Sensing and GIS, 3 (4): 33-46.
  23. Praticò, S., Solano, F., Di Fazio, S., & Modica, G. (2021). Machine learning classification of mediterranean forest habitats in google earth engine based on seasonal sentinel-2 time-series and input image composition optimisation. Remote Sensing, 13(4), 586.
  24. Roujean, J. L., & Breon, F. M. (1995). Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote sensing of Environment, 51(3), 375-384.
  25. Rouse J.W., Haas, R.H., Schell, J.A., & Deering, D.W. (1973). Monitoring vegetation system in the Great Plains with ERTS: Proceedings of the Third Earth Resources Technology Satellite-1 Symposium, Washington DC, 309-317.
  26. Shakeri, Z., Mohammadi-Samani, K., Maarofi, H., Khoonsiavashan, S., & Sharifi, K. (2020). Species diversity, life form, and chorotypes of plant species in sacred groves and
    surrounding silvopastoral woodlands of Northern Zagros, Iran. Iranian Journal of Forest and Poplar Research, 29(2): 101-113.
  27. Shataee, S., Kalbi, S., Fallah, A., & Pelz, D. (2012). Forest attributes imputation using machine-learning methods and ASTER data: comparison of k-NN, SVR and random forest regression algorithms. International journal of remote sensing, 33(19), 6254-6280.
  28. Sibanda, M., Gumede, N., & Mutanga, O. (2021). Estimating leaf area index of the Yellowwood tree (Podocarpus spp.) in an indigenous Southern African Forest, using Sentinel 2 Multispectral Instrument data and the Random Forest regression ensemble. Geocarto International, 1-22.
  29. Sivanpillai, R., Smith, C. T., Srinivasan, R., Messina, M. G., & Wu, X. B. (2006). Estimation of managed loblolly pine stands age and density with Landsat ETM+ data. Forest Ecology and Management, 223(1-3), 247-254.
  30. Smith, F. (1996). Biological diversity, ecosystem stability and economic development. Ecological Economics, 16(3), 191-203.
  31. Townsend, P. A. (2002). Estimating forest structure in wetlands using multitemporal SAR. Remote sensing of environment, 79(2-3), 288-304.
  32. Tucker, C. J. (1979). Red and photographic infrared linear combinations for monitoring vegetation. Remote sensing of Environment, 8(2), 127-150.
  33. Wang, Q., Blackburn, G. A., Onojeghuo, A. O., Dash, J., Zhou, L., Zhang, Y., & Atkinson, P. M. (2017). Fusion of Landsat 8 OLI and Sentinel-2 MSI data. IEEE Transactions on Geoscience and Remote Sensing, 55(7), 3885-3899.
  34. Yim, J. S., Shin, M. Y., Son, Y., & Kleinn, C. (2015). Cluster plot optimization for a large area forest resource inventory in Korea. Forest science and technology, 11(3), 139-146.
  35. Zobeiri, M. (2007). Forest biometrics. Tehran: Tehran University Press.