Acessibilidade / Reportar erro

Determinação do volume de madeira em povoamento de eucalipto por escâner a laser aerotransportado

Determining timber volume of eucalyptus stands by airborne laser scanning

Resumos

O objetivo deste trabalho foi avaliar a possibilidade de se estimar o diâmetro à altura do peito (DAP) com os dados de altura e de número de árvores derivados do escâner a laser aerotransportado (LiDAR, "light detection and ranging"), e determinar o volume de madeira de talhão de Eucalyptus sp. a partir dessas variáveis. O número total de árvores detectadas foi obtido com uso da filtragem de máxima local. A altura de plantas estimada pelo LiDAR apresentou tendência não significativa à subestimativa. A estimativa do DAP foi coerente com os valores encontrados no inventário florestal; porém, também mostrou tendência à subestimativa, em razão do comportamento observado quanto à altura. A variável número de fustes apresentou valores próximos aos observados nas parcelas do inventário. O LiDAR subestimou o volume total de madeira do talhão em 11,4%, em comparação ao volume posto na fábrica. A tendência de subestimação da altura das árvores (em média, cerca de 5%) impactou a estimativa do volume individual de árvores e, consequentemente, a do volume do talhão. No entanto, é possível gerar equações de regressão que estimam o DAP com boa precisão, a partir de dados de altura de plantas obtidos pelo LiDAR. O modelo parabólico é o que possibilita as melhores estimativas da produção volumétrica dos talhões de eucalipto.

Eucalyptus; inventário florestal; LiDAR; modelos biométricos; sensoriamento remoto.


The objective of this work was to evaluate the possibility of estimating the diameter at breast height (DBH) with tree height and number data derived from airborne laser scanning (LiDAR, light detection and ranging) dataset, and to determine the timber volume of an Eucalyptus sp. stand from these variables. The total number of detected trees was obtained using a local maxima filtering. Plant height estimated by LiDAR showed a nonsignificant tendency to underestimation. The estimate for DBH was coherent with the results found in the forest inventory; however, it also showed a tendency towards underestimation due to the observed behavior for height. The variable number of stems showed values close to the ones observed in the inventory plots. LiDAR underestimated the total timber volume in the stand in 11.4%, compared to the total volume delivered to the industry. The underestimation tendency of tree height (5% mean value) impacted the individual tree volume estimate and, consequently, the stand volume estimate. However, it is possible to obtain regression equations that estimate DBH with good precision, from the LiDAR plant height derived data. The parabolic model is the one that provides the best estimates for timber volumetric yield of eucalyptus stands.

Eucalyptus; forest inventory; LiDAR; biometric models; remote sensing.


Introdução

A crescente necessidade de aprimorar as estimativas de inventários florestais tem motivado estudos com dados de sensoriamento remoto, como auxílio à quantificação volumétrica de povoamentos florestais.

O uso do "light detection and ranging" (LiDAR), que é um escâner a laser aerotransportado, para a determinação da elevação de terrenos iniciou-se na segunda metade da década de 1970 (Lillesand & Kiefer, 2000 LILLESAND, T.M.; KIEFER, R.W. Remote sensing and image interpretation. Madison: John Wiley, 2000. 724p.) e, desde os anos 1980, seu emprego para aplicações florestais tem avançado significativamente. A maior disponibilidade de imagens de alta resolução espacial torna viável a obtenção de informações de sensoriamento remoto relativas a árvores individuais, e não mais apenas para talhões ou povoamentos (Disperati et al., 2007DISPERATI, A.A.; SANTOS, J.R. dos; OLIVEIRA FILHO, P.C. de; NEEFF, T. Aplicação da técnica "filtragem de locais máximas" em fotografia aérea digital para a contagem de copas em reflorestamento de Pinus elliottii. Scientia Forestalis, n.76, p.45-55, 2007.).

Atualmente, o LiDAR tem uma gama de aplicações em diversas áreas, como arqueologia (Harmon et al., 2006HARMON, J.M.; LEONE, M.P.; PRINCE, S.D.; SNYDER, M. LiDAR for archaeological landscape analysis: a case study of two eighteenth-century Maryland plantation sites. American Antiquity, v.71, p.649-670, 2006. DOI: 10.2307/40035883.
https://doi.org/10.2307/40035883....
), planejamento costeiro, avaliação de risco de inundações, telecomunicações, redes de transmissão de energia, florestas, agricultura, petróleo, transportes, planejamento urbano e mineração, entre outras (Giongo et al., 2010GIONGO, M.; KOEHLER, H.S.; MACHADO, S. do A.; KIRCHNER, F.F.; MACHETTI, M. LiDAR: princípios e aplicações florestais. Pesquisa Florestal Brasileira, v.30, p.231-244, 2010. DOI: 10.4336/2010.pfb.30.63.231.
https://doi.org/10.4336/2010.pfb.30.63.2...
). No Brasil, no entanto, a linha de pesquisa com o equipamento ainda é incipiente, com destaque para os trabalhos de Castro & Centeno (2005), Castro (2006)CASTRO, F.C. e; CENTENO, T.M. Segmentação de imagens geradas por perfilamento a laser para delimitação de árvores individuais em uma área de reflorestamento de eucaliptos. In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 12., 2005, Goiânia. Anais. Goiânia: INPE, 2005. p.737-744., Zandoná et al. (2008) ZANDONÁ, D.F.; LINGNAU, C.; NAKAJIMA, N.Y. Varredura a laser aerotransportado para estimativa de variáveis dendrométricas., Scientia Forestalis v.36, p.295-306, 2008., Macedo (2009)MACEDO, R. de C. Estimativa volumétrica de povoamento clonal de Eucalyptus sp. através de laser scanner aerotransportado. 2009. 143p. Dissertação (Mestrado) - Instituto Nacional de Pesquisas Espaciais, São José dos Campos., Rodriguez et al. (2010)RODRIGUEZ, L.C.E.; POLIZEL, J.L.; FERRAZ, S.F.B.; ZONETE, M.F.; FERREIRA, M.Z.Inventário florestal com tecnologia laser aerotransportada de plantios de Eucalyptus spp. no Brasil. Ambiência, v.6, p.67-80, 2010., Zonete (2009)ZONETE, M.F. Análise do uso da tecnologia laser aerotransportado para inventários florestais em plantios clonais de Eucalyptus sp. no sul da Bahia2009. 95p. Dissertação (Mestrado) - Escola Superior de Agricultura Luiz de Queiroz, Piracicaba. e Zonete et al. (2010)ZONETE, M.F.; RODRIGUEZ, R.C.E.; PACKALÉN, P. Estimação de parâmetros biométricos de plantios clonais de eucalipto no sul da Bahia: uma aplicação da tecnologia laser aerotransportada., Scientia Forestalis v.38, p.225-235, 2010.. No mundo, há vários estudos para estimar parâmetros de parcelas com o uso do LiDAR, para a diminuição da intensidade amostral. Esses estudos têm sido desenvolvidos com contagem de árvores individuais em campo, com a quantificação do estoque florestal e da biomassa, como nos trabalhos de Naesset & Bjerknes (2001)NAESSET, E.; BJERKNES, K.O. Estimating tree heights and number of stems in young forest stands using airborne laser scanner data., Remote Sensing of Environment v.78, p.328-340, 2001. DOI: 10.1016/S0034-4257(01)00228-0.
https://doi.org/10.1016/S0034-4257(01)00...
, Popescu et al. (2002, 2003)POPESCU, S.C.;; WYNNE, R.H. NELSON, R.F.Measuring individual tree crown diameter with LiDAR and assessing its influence on estimating forest volume and biomass. Canadian Journal of Remote Sensing, v.29, p.564-577, 2003. DOI: 10.5589/m03-027
https://doi.org/10.5589/m03-027...
, Maltamo et al. (2004)MALTAMO, M.; EERIKÄINEN, K.; PITKÄNEN, J.; HYYPPÄ, J.; VEHMAS, M. Estimation of timber volume and stem density based on scanning laser altimetry and expected tree size distribution functions. Remote Sensing of Environment, v.90, p.319-330, 2004. DOI: 10.1016/j.rse.2004.01.006.
https://doi.org/10.1016/j.rse.2004.01.00...
, Naesset (2004)NAESSET, E. Practical large-scale forest stand inventory using a small-footprint airborne scanning laser., Scandinavian Journal of Forest Research v.19, p.164-179, 2004. DOI: 10.1080/02827580310019257.
https://doi.org/10.1080/0282758031001925...
, Yu et al. (2004)YU, X.W.; HYYPPÄ, J. KAARTINEN, H.; MALTAMO, M.Automatic detection of harvested trees and determination of forest growth using airborne laser scanning., Remote Sensing of Environment v.90, p.451-462, 2004. DOI: 10.1016/j.rse.2004.02.001.
https://doi.org/10.1016/j.rse.2004.02.00...
, Tiede et al. (2005)TIEDE, D.; HOCHLEITNER, G.; BLASCHKE, T. A full GIS-based workflow for tree identification and tree crown delineation using laser scanning. In: STILLA, U.; ROTTENSTEINER, F.; HINZ, S. (Ed.). 3D city models, road databases, and traffic monitoring: concepts, algorithms, and evaluation. Vienna: IAPRS, 2005. p.9-14., Crow et al. (2007)CROW, P.; BENHAM, S.; DEVEREUX, B.J.; AMABLE, G.S. Woodland vegetation and its implications for archaeological survey using LiDAR. Forestry, v.80, p.241-252, 2007. DOI: 10.1093/forestry/cpm018.
https://doi.org/10.1093/forestry/cpm018....
, Jupp et al. (2007) JUPP, D.; CULVENOR, D.; LOVELL, G.; NEWHAM, G.; STRAHLER, A.; WOODCOCK, C. Estimating forest LAI profiles and structural parameters using a ground-based laser called Echidna. Tree Physiology, v.29, p.171-181, 2007. DOI: 10.1093/treephys/tpn022.
https://doi.org/10.1093/treephys/tpn022...
, Naesset (2007)NAESSET, E. Airbone laser scanning as method in operational foresty inventory: status of accuracy assessments accomplished in Scandinavia. Scandinavian Journal of Forest Research, v.22, p.433-442, 2007. DOI: 10.1080/02827580701672147.
https://doi.org/10.1080/0282758070167214...
, Heurich & Thoma (2008)HEURICH, M.; THOMA, F. Estimation of forestry stand parameters using laser scanning data in temperate, structurally rich natural European beech (Fagus sylvatica) and Norway spruce (Picea abies) forests., Forestry v.81, p.645-661, 2008. DOI: 10.1093/forestry/cpn038.
https://doi.org/10.1093/forestry/cpn038....
, Lauri et al. (2008)LAURI, K.; JUSSI, P.; JUKKA, M.; AKI, S.; MATTI, M.; PETTERI, P.; JYRKI, K. The use of airborne laser scanning to estimate sawlog volumes., Forestry v.81, p.499-509, 2008. DOI: 10.1093/forestry/cpn018.
https://doi.org/10.1093/forestry/cpn018....
Sasaki et al. (2008)SASAKI, T.; IMANISHI, J. IOKI, K.; MORIMOTO, Y. KITADA, K.Estimation of leaf area index and canopy openness in broad-leaved forest using airborne laser scanner in comparison with high-resolution near-infrared digital photography. Landscape and Ecological Engineeringv.2, p.47-55, 2008. DOI: 10.1007/s11355-008-0041-8.
https://doi.org/10.1007/s11355-008-0041-...
, e Ioki et al. (2010)IOKI, K.; JUNICHI, I.; SASAKI, T.; MORIMOTO, Y.; KITADA, K. Estimating stand volume in broad-leaved forest using discrete-return LiDAR: plot-based approach. Landscape and Ecological Engineering v.6, p.29-36, 2010. DOI: 10.1007/s11355-009-0077-4.
https://doi.org/10.1007/s11355-009-0077-...
.

O desenvolvimento de novas técnicas, bem como a validação das já existentes para condições brasileiras, é importante para consolidar o uso desta tecnologia no setor florestal. Entre as variáveis que podem ser obtidas com dados do LiDAR estão altura, densidade de plantio, diâmetro de copa, distribuição espacial de plantas, índices de heterogeneidade dos talhões, ocorrências de falhas em reboleiras, entre outras. O volume de madeira do povoamento é a variável mais valiosa para a gestão florestal, tanto no âmbito estratégico como no operacional. A altura das árvores é uma variável dendrométrica de difícil aquisição no inventário florestal, mas apresenta forte correlação com o volume de madeira de um talhão e, portanto, é de grande importância para a determinação do incremento volumétrico em talhões homogêneos (Skovsgaard & Vanclay, 2008SKOVSGAARD, J.P.; VANCLAY, J.K. Forest site productivity: a review of the evolution of dendrometric concepts for even-aged stands., Forestry v.81, p.13-31, 2008. DOI: 10.1093/forestry/cpm041.
https://doi.org/10.1093/forestry/cpm041....
).

Zandoná et al. (2008) ZANDONÁ, D.F.; LINGNAU, C.; NAKAJIMA, N.Y. Varredura a laser aerotransportado para estimativa de variáveis dendrométricas., Scientia Forestalis v.36, p.295-306, 2008. utilizaram informações de altura derivadas de LiDAR para estimar o diâmetro de copa por meio de regressão linear. Posteriormente, uma segunda regressão, que utilizava o diâmetro de copa e a altura como parâmetros, foi utilizada pelos autores para estimar a variável dendrométrica diâmetro à altura do peito (DAP). Os autores então utilizaram os percentis 10, 30, 50, 70 e 90, para classificar as estimativas das variáveis biométricas, que foram utilizadas em modelos para estimar o volume de árvores individuais. Modelos de regressão permitem estimar o DAP a partir da altura das árvores. Esses modelos são tradicionalmente usados de forma inversa; assim, para a determinação da altura a partir do volume, Ribeiro et al. (2010) RIBEIRO, A.; FERRAZ FILHO, A.C.; MELLO, J.M.;; FERREIRA, M.Z. LISBOA, P.M.M.; SCOLFORO, J.R.S. Estratégias e metodologias de ajuste de modelos hipsométricos em plantios de Eucalyptus sp., Cerne v.16, p.22-31, 2010. utilizaram o modelo de Curtis e obtiveram bons resultados na determinação da altura a partir do DAP.

O objetivo deste trabalho foi avaliar a possibilidade de se estimar o diâmetro à altura do peito (DAP), com os dados de altura e de número de árvores derivados do escâner a laser aerotransportado (LiDAR, "light detection and ranging"), e determinar o volume de talhão de Eucalyptus sp. a partir dessas variáveis.

Material e Métodos

A área experimental está situada no extremo sul da Bahia, entre as coordenadas 16º05'10"S, 39º25'30"W e 16º05'40"S, 39º25'00"W (Figura 1). O clima da região, conforme classificação de Köppen, é do tipo tropical equatorial, com chuvas distribuídas durante o ano inteiro e média anual de 1.200 mm. A temperatura média é de 24ºC, com pequena amplitude, por se tratar de região costeira. A área compreende um talhão de Eucalyptus sp. de 23,5 ha, plantado no espaçamento de 5,00 x 2,40 m, com sete anos de idade à época do voo.

Figura 1.
Localização geográfica do talhão escaneado com o LiDAR e parcelas circulares usadas no inventário florestal (em vermelho).

O talhão foi inventariado na mesma semana do voo com o LiDAR, tendo-se considerado uma suficiência amostral de quatro parcelas. O número de parcelas foi definido com vistas à obtenção de um erro no inventário inferior a 10%, para a variável volume, e uma superfície de análise de pelo menos 2.000 m2. Assim, quatro parcelas circulares de 530,93 m2, com o centro georreferenciado com DGPS, foram distribuídas aleatoriamente pelo talhão, conforme os círculos em vermelho da Figura 1. Nessas parcelas, foi feita a medição do DAP, em todas as árvores, e da altura das seis primeiras árvores. Para a estimativa das alturas totais das demais árvores, utilizou-se a relação hipsométrica proposta por Campos et al. (1984CAMPOS, J.C.C.; RIBEIRO, J.O.; PAULA NETO, F. Inventário florestal nacional, reflorestamento: Minas Gerais. Brasília: Instituto Brasileiro de Desenvolvimento Florestal, 1984. 126p.), ajustada para o talhão em estudo, pela equação HTest = exp(β0 + β1/DAP + β2Ln Hdom) + ε, em que: HTest é a altura estimada; DAP é o diâmetro à altura do peito (1,30 m do solo); Hdom é a altura dominante; βi representa os parâmetros estimados; e ε é o erro.

Um teste t pareado foi realizado para averiguar se havia diferenças significativas entre os valores médios de altura em cada parcela, obtidos pelo LiDAR e pelo inventário florestal.

A estimativa de volume foi realizada pelo modelo de Schumacher & Hall por Vi = β0 DAPβ1 HTβ2 + ε, em que: Vi é o volume total estimado na iésima árvore; DAP é o diâmetro à altura do peito; HT é a altura total, βi representa os parâmetros estimados; e ε é o erro. A equação foi ajustada com informações de 65 árvores, cubadas rigorosamente para o clone existente no talhão, com a mesma idade e em diferentes classes diamétricas.

Os dados do LiDAR foram coletados com o sistema ALTM 3100, acoplado a um avião bimotor Piper Seneca II. O voo foi realizado à velocidade de 150 km h-1, à altura de 1.000 m, diâmetro do feixe de 20 cm, divergência do feixe de 0,3 mrad, ângulo de varredura de 15º, frequência de varredura de 58,7 Hz e intensidade de 1,5 ponto m-2. O erro planimétrico foi de 0,5 m (1,0s), o altimétrico de 0,15 m (1,0 s), e a área de cobertura do laser foi de 0,32 m. Os dados dos sinais do LiDAR referentes às áreas dos talhões foram organizados e calibrados com base nas informações de latitude, longitude e altitude (x, y, z), pulsos (1º, 2º e 3º retornos) e intensidade do sinal. Esses fatores podem influenciar a densidade da nuvem de pontos resultante, que pode variar de 0,5 a 5 pontos m-2. A variação é decorrente do tipo de produto encomendado, da altura e da velocidade de voo, e de características do sensor laser. A precisão altimétrica do LiDAR é de aproximadamente 15 cm, segundo Crow et al. (2007)CROW, P.; BENHAM, S.; DEVEREUX, B.J.; AMABLE, G.S. Woodland vegetation and its implications for archaeological survey using LiDAR. Forestry, v.80, p.241-252, 2007. DOI: 10.1093/forestry/cpm018.
https://doi.org/10.1093/forestry/cpm018....
.

Após a estruturação das informações, analisou-se a sobreposição lateral de 30% entre as linhas de voo e a rota de voo da aeronave, para separar os pontos sobrepostos. As informações de sobreposição não foram consideradas no presente estudo. O ruído advindo da sobreposição pode ter sido causado pelo provável deslocamento do dossel pelo vento, nas diferentes passagens da aeronave, em que a interporlação pelo método IDW ("inverse distance weighted interpolation") são mais realçadas (Oliveira et al., 2012)OLIVEIRA, L.T. de; CARVALHO, L.M.T. de; FERREIRA, M.Z.; OLIVEIRA, T.C. de A.; ACERBI JUNIOR, F.W. Application of LIDAR to forest inventory for tree count in stands of Eucalyptus sp. Cerne, v.18, p.75-184, 2012. DOI: 10.1590/S0104-77602012000200001.
https://doi.org/10.1590/S0104-7760201200...
. Os dados foram processados conforme Oliveira et al. (2012)OLIVEIRA, L.T. de; CARVALHO, L.M.T. de; FERREIRA, M.Z.; OLIVEIRA, T.C. de A.; ACERBI JUNIOR, F.W. Application of LIDAR to forest inventory for tree count in stands of Eucalyptus sp. Cerne, v.18, p.75-184, 2012. DOI: 10.1590/S0104-77602012000200001.
https://doi.org/10.1590/S0104-7760201200...
, tendo-se analisado duas classes de elevação, segundo o modelo digital de superfície (MDS): na primeira classe, todas as informações referentes ao primeiro pulso acima de 10 m de altura; e, na segunda, apenas as informações referentes ao primeiro pulso acima de 1 m.

A altura das árvores foi obtida por meio da diferença entre os dados fornecidos pelo modelo digital de terreno (MDT) e os do MDS. Para isso, utilizaram-se regras de lógica binária, tendo-se avaliado as informações de pulso, altitude e angulação vertical e horizontal entre os pontos vizinhos, a partir do desvio-padrão das altitudes de cada pulso, do número de árvores e de suas respectivas alturas, obtidas pelo método proposto por Oliveira et al. (2012)OLIVEIRA, L.T. de; CARVALHO, L.M.T. de; FERREIRA, M.Z.; OLIVEIRA, T.C. de A.; ACERBI JUNIOR, F.W. Application of LIDAR to forest inventory for tree count in stands of Eucalyptus sp. Cerne, v.18, p.75-184, 2012. DOI: 10.1590/S0104-77602012000200001.
https://doi.org/10.1590/S0104-7760201200...
. De acordo com esse método, os dados brutos do LiDAR são organizados em função das informações x, y, z e da separação dos retornos (1º., 2º. e 3º.), e classificados como solo, vegetação acima de 1 m de altura e vegetação acima de 10 m, considerando-se os dados de sobreposição do voo, que são retirados do processamento. A altura é calculada para todos os pontos, em relação ao MDT, e utilizada para geração do MDS considerando-se apenas o primeiro retorno; este modelo é rasterizado pelo método IDW. Este raster é submetido a um filtro de local máxima, com uso de janela móvel de 5x5, de onde são extraídos o número e a altura das árvores detectadas.

Para a determinação do volume, ajustou-se o modelo de regressão para estimação do DAP a partir da altura das árvores obtida pelo LiDAR. Para isso, foram testados quatro modelos tradicionais para a relação hipsométrica, tomados em sua forma inversa: modelo parabólico, DAP = β0 + β1HT + β2HT2 + εi; modelo da linha reta, DAP = β0 + β1HT + εi; modelo de Stofel, Log(DAP) = β0 + β1 Log(HT) + Log(εi); e modelo de Curtis, DAP = β0 + β1 1/Log(HT) + εi.

Para a seleção do melhor modelo, avaliou-se o gráfico de resíduos, o erro-padrão da estimativa (% e cm) e, por último, o coeficiente de determinação R2. Todos os modelos foram ajustados conforme os dados do inventário florestal, tendo-se utilizado a melhor equação para a estimativa do DAP, nas duas classes de dados do LiDAR (classe 1, primeiro pulso acima de 10 m; e classe 2, primeiro pulso acima de 1 m), interpolados pelo IDW por meio de uma janela móvel 5x5 (Oliveira et al., 2012OLIVEIRA, L.T. de; CARVALHO, L.M.T. de; FERREIRA, M.Z.; OLIVEIRA, T.C. de A.; ACERBI JUNIOR, F.W. Application of LIDAR to forest inventory for tree count in stands of Eucalyptus sp. Cerne, v.18, p.75-184, 2012. DOI: 10.1590/S0104-77602012000200001.
https://doi.org/10.1590/S0104-7760201200...
).

Após a estimativa da altura, aplicou-se a equação ajustada para o modelo volumétrico de Schumacher & Hall.

Os resultados das duas classes obtidas com LiDAR foram comparadas com as informações obtidas pelo inventário, para todo o talhão, com o registro da contagem do número de árvores por hectare, cortadas pelo harvester, e com o volume entregue na fábrica. Neste último caso, a madeira foi pesada e sua densidade foi estimada para a obtenção do volume; a medição ocorreu 35 dias após a colheita. A Figura 2 apresenta um resumo da sequência metodológica empregada neste estudo.

Figura 2.
Fluxograma para a execução da metodologia proposta.

Resultados e Discussão

As médias de altura obtidas pelo inventário florestal e pelo LiDAR, para a classe 1 (média = 1,65, DP = 1,738, N = 4) e para a classe 2 (média = 1,70, DP = 1,676, N = 4), não diferiram significativamente [t(3) = 1,90 e p bicaudal = 0,154, para a classe 1; e t(3) = 2,03 e p bicaudal = 0,135, para a classe 2)]. As Tabelas 1 e 2 apresentam uma síntese das informações obtidas no inventário florestal, realizado no talhão à época do voo. Os modelos de relação hipsométrica apresentaram boas estatísticas de ajuste, com alto coeficiente de determinação (média de 92%) e baixos valores de erro (média do erro-padrão residual de 1,1 cm). Estes valores estão de acordo com estudos com relação hipsométrica em eucalipto, tais como os de Soares et al. (2004) SOARES, T.S.;SCOLFORO, J.R.S. FERREIRA, S.O.; MELLO, J.M.Uso de diferentes alternativas para viabilizar a relação hipsométrica no povoamento florestal. Revista Árvore, v.28, p.845-854, 2004. DOI: 10.1590/S0100-67622004000600009.
https://doi.org/10.1590/S0100-6762200400...
e Ribeiro et al. (2010) RIBEIRO, A.; FERRAZ FILHO, A.C.; MELLO, J.M.;; FERREIRA, M.Z. LISBOA, P.M.M.; SCOLFORO, J.R.S. Estratégias e metodologias de ajuste de modelos hipsométricos em plantios de Eucalyptus sp., Cerne v.16, p.22-31, 2010..

Tabela 1.
Síntese das informações obtidas pelo inventário de um talhão de Eucalyptus sp.

Tabela 2.
Síntese das medidas de precisão, obtidas pelo inventário, para a variável volume de madeira.

Os dados de R2 e do erro-padrão residual apontaram para o modelo parabólico como o de melhor ajuste (Tabela 3). O gráfico de dispersão (Figura 3) também indica que houve melhor ajuste do modelo parabólico. Assim, o modelo foi selecionado para estimar o DAP das árvores a partir de seus dados de altura obtidos pelo LiDAR. Esse modelo, por ser essencialmente matemático, tem aplicação restrita ao intervalo de dados ao qual foi ajustado.

Tabela 3.
Parâmetros para o ajuste dos modelos analisados.

Figura 3.
Dispersão dos resíduos e inclinação da reta para os modelos: A e B, parabólico; C e D, linha reta; E e F, Stofel; e G e H, Curtis.

Os valores de DAP estimados apresentaram boa aderência aos valores obtidos no inventário florestal, nas parcelas 1 e 2 (Tabela 4), com diferença média de -4%. Já os valores de DAP para as parcelas 3 e 4 foram subestimados em 17,8%. Esta subestimativa do DAP ocorreu em consequência da subestimativa da altura média pelo LiDAR, nas parcelas 3 e 4, que foi, em média, 9,8% menor do que a altura obtida no inventário. Assim, se a variável altura for corretamente estimada pelo LiDAR, espera-se um bom acerto na variável DAP.

Tabela 4.
Valores médios de altura, obtidos com o LiDAR, e DAP, estimado com o modelo parabólico, e valores médios de altura e DAP obtidos pelo inventário florestal.

O modelo para estimar o volume de madeira, a partir do DAP (V = 0,000041DAP1,76383 × HT1,165339), apresentou os seguintes parâmetros: R2 ajustado, 98,40; erro-padrão residual (m3), 0,020851; e CV, 5,3%. Os dados do LiDAR na classe 2 foram os que apresentaram menor diferença quanto à estimativa da produção volumétrica, em comparação à produção medida na fábrica (Tabela 5). Todos os dados volumétricos para o talhão, excetuando-se os da classe 1, estiveram dentro do intervalo de confiança do inventário (Tabela 2); porém, a classe 1 excedeu o limite inferior em 256,28 m3 e gerou um erro de subestimativa de 14% em relação ao volume posto na fábrica.

Tabela 5.
Número de árvores e volume, encontrados com contagem visual, LiDAR, inventário, inventário pré-corte e "harvester".

O volume obtido com "harvester" foi um pouco mais alto do que o posto em fábrica (Tabela 5). Isso pode ocorrer, pois este equipamento, quando em operação, tende a contar árvores mais de uma vez, quando elas se quebram durante a operação. A tendência da contagem individual, de subestimar o número de árvores encontradas nas parcelas do inventário florestal, não era esperada e deve ter ocorrido em razão da supressão de árvores muito próximas na contagem individual.

A subestimativa para o volume estimado com os dados do LiDAR, em relação aos dados do inventário florestal, pode ter sido provocada pela subestimativa das alturas detectadas nos pontos de máxima. Assim, em razão da altura subestimada, o volume também foi subestimado em relação ao posto na fábrica. No entanto, o volume total da classe 2 ficou dentro do intervalo de confiança do inventário. A altura é uma variável dendrométrica de grande importância para o inventário florestal (Skovsgaard & Vanclay, 2008SKOVSGAARD, J.P.; VANCLAY, J.K. Forest site productivity: a review of the evolution of dendrometric concepts for even-aged stands., Forestry v.81, p.13-31, 2008. DOI: 10.1093/forestry/cpm041.
https://doi.org/10.1093/forestry/cpm041....
) que, se não tivesse sido subestimada pela contagem automática, teria proporcionado estimativas de volume mais próximas do real. A existência de sobreposição ou de árvores muito próximas foi responsável pela subestimativa do número de árvores captado pelo LiDAR, o que também foi relatado por Oliveira et al. (2012)OLIVEIRA, L.T. de; CARVALHO, L.M.T. de; FERREIRA, M.Z.; OLIVEIRA, T.C. de A.; ACERBI JUNIOR, F.W. Application of LIDAR to forest inventory for tree count in stands of Eucalyptus sp. Cerne, v.18, p.75-184, 2012. DOI: 10.1590/S0104-77602012000200001.
https://doi.org/10.1590/S0104-7760201200...
, e pelo fato de os pontos de máxima não terem captado os picos das árvores nas classes 1 e 2. Assim, é possível que o LiDAR acerte o número total das árvores no talhão, em razão de erros compensativos no decorrer da análise. Macedo (2009)MACEDO, R. de C. Estimativa volumétrica de povoamento clonal de Eucalyptus sp. através de laser scanner aerotransportado. 2009. 143p. Dissertação (Mestrado) - Instituto Nacional de Pesquisas Espaciais, São José dos Campos. verificou que, apesar de os algoritmos para estimativa volumétrica que utilizam dados LiDAR subestimarem a quantidade de árvores e a altura, e superestimarem o DAP, a estimativa volumétrica se mantém dentro de um intervalo aceitável, tanto em termos de intervalo de confiança quanto de erro amostral.

Em geral, o deslocamento da copa das árvores, em áreas limítrofes das parcelas, causado pelo vento nas diferentes passagens da aeronave, faz com que elas sejam suprimidas ou acrescentadas ao cômputo total de árvores em cada parcela, o que seria responsável por erros, ainda que pequenos. Esses erros são responsáveis pela diferença entre o valor das informações observadas no inventário, na contagem visual e na detecção automática. Apesar disso, os dados indicam boa precisão quanto ao volume total do talhão para a classe 2. Portanto, o LiDAR é uma ferramenta promissora para o setor florestal, embora ainda demande estudos para sua correta implementação.

Conclusões

  1. É possível gerar equações de regressão para estimar o DAP com boa precisão a partir de dados de altura obtidos pelo LiDAR, e o modelo parabólico é o que possibilita as melhores estimativas.

  2. A subestimativa da altura das árvores detectadas pelo LiDAR é a principal responsável pela diferença entre o volume de madeira estimado com o uso equipamento e o volume posto em fábrica.

Agradecimentos

À Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes) e à Fundação de Amparo à Pesquisa do Estado de Minas Gerais (Fapemig), pelas bolsas de estudo.

  • CAMPOS, J.C.C.; RIBEIRO, J.O.; PAULA NETO, F. Inventário florestal nacional, reflorestamento: Minas Gerais. Brasília: Instituto Brasileiro de Desenvolvimento Florestal, 1984. 126p.
  • CASTRO, F. da C. e. Segmentação de imagens geradas por laser scanning aerotransportado para delimitação de árvores individuais em áreas de reflorestamento de eucaliptos. 2006. 76p. Dissertação (Mestrado) - Universidade Tecnológica Federal do Paraná, Curitiba.
  • CASTRO, F.C. e; CENTENO, T.M. Segmentação de imagens geradas por perfilamento a laser para delimitação de árvores individuais em uma área de reflorestamento de eucaliptos. In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 12., 2005, Goiânia. Anais. Goiânia: INPE, 2005. p.737-744.
  • CROW, P.; BENHAM, S.; DEVEREUX, B.J.; AMABLE, G.S. Woodland vegetation and its implications for archaeological survey using LiDAR. Forestry, v.80, p.241-252, 2007. DOI: 10.1093/forestry/cpm018.
    » https://doi.org/10.1093/forestry/cpm018.
  • DISPERATI, A.A.; SANTOS, J.R. dos; OLIVEIRA FILHO, P.C. de; NEEFF, T. Aplicação da técnica "filtragem de locais máximas" em fotografia aérea digital para a contagem de copas em reflorestamento de Pinus elliottii. Scientia Forestalis, n.76, p.45-55, 2007.
  • GIONGO, M.; KOEHLER, H.S.; MACHADO, S. do A.; KIRCHNER, F.F.; MACHETTI, M. LiDAR: princípios e aplicações florestais. Pesquisa Florestal Brasileira, v.30, p.231-244, 2010. DOI: 10.4336/2010.pfb.30.63.231.
    » https://doi.org/10.4336/2010.pfb.30.63.231.
  • HARMON, J.M.; LEONE, M.P.; PRINCE, S.D.; SNYDER, M. LiDAR for archaeological landscape analysis: a case study of two eighteenth-century Maryland plantation sites. American Antiquity, v.71, p.649-670, 2006. DOI: 10.2307/40035883.
    » https://doi.org/10.2307/40035883.
  • HEURICH, M.; THOMA, F. Estimation of forestry stand parameters using laser scanning data in temperate, structurally rich natural European beech (Fagus sylvatica) and Norway spruce (Picea abies) forests., Forestry v.81, p.645-661, 2008. DOI: 10.1093/forestry/cpn038.
    » https://doi.org/10.1093/forestry/cpn038.
  • IOKI, K.; JUNICHI, I.; SASAKI, T.; MORIMOTO, Y.; KITADA, K. Estimating stand volume in broad-leaved forest using discrete-return LiDAR: plot-based approach. Landscape and Ecological Engineering v.6, p.29-36, 2010. DOI: 10.1007/s11355-009-0077-4.
    » https://doi.org/10.1007/s11355-009-0077-4.
  • JUPP, D.; CULVENOR, D.; LOVELL, G.; NEWHAM, G.; STRAHLER, A.; WOODCOCK, C. Estimating forest LAI profiles and structural parameters using a ground-based laser called Echidna. Tree Physiology, v.29, p.171-181, 2007. DOI: 10.1093/treephys/tpn022.
    » https://doi.org/10.1093/treephys/tpn022
  • LAURI, K.; JUSSI, P.; JUKKA, M.; AKI, S.; MATTI, M.; PETTERI, P.; JYRKI, K. The use of airborne laser scanning to estimate sawlog volumes., Forestry v.81, p.499-509, 2008. DOI: 10.1093/forestry/cpn018.
    » https://doi.org/10.1093/forestry/cpn018.
  • LILLESAND, T.M.; KIEFER, R.W. Remote sensing and image interpretation. Madison: John Wiley, 2000. 724p.
  • MACEDO, R. de C. Estimativa volumétrica de povoamento clonal de Eucalyptus sp. através de laser scanner aerotransportado. 2009. 143p. Dissertação (Mestrado) - Instituto Nacional de Pesquisas Espaciais, São José dos Campos.
  • MALTAMO, M.; EERIKÄINEN, K.; PITKÄNEN, J.; HYYPPÄ, J.; VEHMAS, M. Estimation of timber volume and stem density based on scanning laser altimetry and expected tree size distribution functions. Remote Sensing of Environment, v.90, p.319-330, 2004. DOI: 10.1016/j.rse.2004.01.006.
    » https://doi.org/10.1016/j.rse.2004.01.006.
  • NAESSET, E. Airbone laser scanning as method in operational foresty inventory: status of accuracy assessments accomplished in Scandinavia. Scandinavian Journal of Forest Research, v.22, p.433-442, 2007. DOI: 10.1080/02827580701672147.
    » https://doi.org/10.1080/02827580701672147.
  • NAESSET, E. Practical large-scale forest stand inventory using a small-footprint airborne scanning laser., Scandinavian Journal of Forest Research v.19, p.164-179, 2004. DOI: 10.1080/02827580310019257.
    » https://doi.org/10.1080/02827580310019257.
  • NAESSET, E.; BJERKNES, K.O. Estimating tree heights and number of stems in young forest stands using airborne laser scanner data., Remote Sensing of Environment v.78, p.328-340, 2001. DOI: 10.1016/S0034-4257(01)00228-0.
    » https://doi.org/10.1016/S0034-4257(01)00228-0.
  • OLIVEIRA, L.T. de; CARVALHO, L.M.T. de; FERREIRA, M.Z.; OLIVEIRA, T.C. de A.; ACERBI JUNIOR, F.W. Application of LIDAR to forest inventory for tree count in stands of Eucalyptus sp. Cerne, v.18, p.75-184, 2012. DOI: 10.1590/S0104-77602012000200001.
    » https://doi.org/10.1590/S0104-77602012000200001.
  • POPESCU, S.C.; WYNNE, R.H.; NELSON, R.F. Estimating plot-level tree heights with lidar: local filtering with a canopy-height based variable window size. Computers and Electronics in Agriculture, v.37, p.71-95, 2002. DOI: 10.1016/S0168-1699(02)00121-7.
    » https://doi.org/10.1016/S0168-1699(02)00121-7.
  • POPESCU, S.C.;; WYNNE, R.H. NELSON, R.F.Measuring individual tree crown diameter with LiDAR and assessing its influence on estimating forest volume and biomass. Canadian Journal of Remote Sensing, v.29, p.564-577, 2003. DOI: 10.5589/m03-027
    » https://doi.org/10.5589/m03-027
  • RIBEIRO, A.; FERRAZ FILHO, A.C.; MELLO, J.M.;; FERREIRA, M.Z. LISBOA, P.M.M.; SCOLFORO, J.R.S. Estratégias e metodologias de ajuste de modelos hipsométricos em plantios de Eucalyptus sp., Cerne v.16, p.22-31, 2010.
  • RODRIGUEZ, L.C.E.; POLIZEL, J.L.; FERRAZ, S.F.B.; ZONETE, M.F.; FERREIRA, M.Z.Inventário florestal com tecnologia laser aerotransportada de plantios de Eucalyptus spp. no Brasil. Ambiência, v.6, p.67-80, 2010.
  • SASAKI, T.; IMANISHI, J. IOKI, K.; MORIMOTO, Y. KITADA, K.Estimation of leaf area index and canopy openness in broad-leaved forest using airborne laser scanner in comparison with high-resolution near-infrared digital photography. Landscape and Ecological Engineeringv.2, p.47-55, 2008. DOI: 10.1007/s11355-008-0041-8.
    » https://doi.org/10.1007/s11355-008-0041-8
  • SKOVSGAARD, J.P.; VANCLAY, J.K. Forest site productivity: a review of the evolution of dendrometric concepts for even-aged stands., Forestry v.81, p.13-31, 2008. DOI: 10.1093/forestry/cpm041.
    » https://doi.org/10.1093/forestry/cpm041.
  • SOARES, T.S.;SCOLFORO, J.R.S. FERREIRA, S.O.; MELLO, J.M.Uso de diferentes alternativas para viabilizar a relação hipsométrica no povoamento florestal. Revista Árvore, v.28, p.845-854, 2004. DOI: 10.1590/S0100-67622004000600009.
    » https://doi.org/10.1590/S0100-67622004000600009
  • TIEDE, D.; HOCHLEITNER, G.; BLASCHKE, T. A full GIS-based workflow for tree identification and tree crown delineation using laser scanning. In: STILLA, U.; ROTTENSTEINER, F.; HINZ, S. (Ed.). 3D city models, road databases, and traffic monitoring: concepts, algorithms, and evaluation. Vienna: IAPRS, 2005. p.9-14.
  • YU, X.W.; HYYPPÄ, J. KAARTINEN, H.; MALTAMO, M.Automatic detection of harvested trees and determination of forest growth using airborne laser scanning., Remote Sensing of Environment v.90, p.451-462, 2004. DOI: 10.1016/j.rse.2004.02.001.
    » https://doi.org/10.1016/j.rse.2004.02.001
  • ZANDONÁ, D.F.; LINGNAU, C.; NAKAJIMA, N.Y. Varredura a laser aerotransportado para estimativa de variáveis dendrométricas., Scientia Forestalis v.36, p.295-306, 2008.
  • ZONETE, M.F. Análise do uso da tecnologia laser aerotransportado para inventários florestais em plantios clonais de Eucalyptus sp. no sul da Bahia2009. 95p. Dissertação (Mestrado) - Escola Superior de Agricultura Luiz de Queiroz, Piracicaba.
  • ZONETE, M.F.; RODRIGUEZ, R.C.E.; PACKALÉN, P. Estimação de parâmetros biométricos de plantios clonais de eucalipto no sul da Bahia: uma aplicação da tecnologia laser aerotransportada., Scientia Forestalis v.38, p.225-235, 2010.

Datas de Publicação

  • Publicação nesta coleção
    Set 2014

Histórico

  • Recebido
    07 Abr 2014
  • Aceito
    19 Ago 2014
Embrapa Secretaria de Pesquisa e Desenvolvimento; Pesquisa Agropecuária Brasileira Caixa Postal 040315, 70770-901 Brasília DF Brazil, Tel. +55 61 3448-1813, Fax +55 61 3340-5483 - Brasília - DF - Brazil
E-mail: pab@embrapa.br