Ciências Florestais

# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Processamento de dados LiDAR ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
# Autor: Luiz Carlos Estraviz Rodriguez
#        Departamento de Ciências Florestais
#        ESALQ/USP - 05/Jun/2022
#
# Linguagem de programação: R
#                           pacote lidR (c) Jean Romain Roussel
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Define diretório e nome da nuvem de pontos LiDAR ~~~~~~~~~~~~~~~~~~~~
lidarDir <- "C:/LiDAR/"
lidarNuv <- "interpine.laz"

# Define nome completo do aquivo com a nuvem de pontos LiDAR ~~~~~~~~~~
nuvemArq <- paste0(lidarDir, lidarNuv)

# Instala (se ainda não estiver instalado) e carrega o pacote lidR ~~~~
if(!require(lidR)){
  install.packages('lidR')
}
library(lidR)

# Leitura parcial dos dados LiDAR ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
las <- readLAS(nuvemArq,   # lê apenas XYZ, intensidade e classificação
               select = "xyzic")  

# Confere a 'sanidade" da nuvem original de dados LiDAR ~~~~~~~~~~~~~~~
las_check(las)

# Características originais da nuvem de dados LiDAR ~~~~~~~~~~~~~~~~~~~
summary(las)

# Colunas de dados extraídas da nuvem original ~~~~~~~~~~~~~~~~~~~~~~~~
print(las@data)                                         # head(las@data)

# Estatísticas básicas dos dados extraídos da nuvem original ~~~~~~~~~~
summary(las@data)

# Plot da nuvem colorida por altura ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
plot(las,
     color  = "Z",
     bg     = "black")

# Plot da nuvem colorida por classificação ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# os pontos azuis representam pontos de solo ~~~~~~~~~~~~~~~~~~~~~~~~~~
plot(las,
     size   = 3,
     color  = "Classification", 
     bg     = "gray", axis = TRUE)
Para mais detalhes sobre o uso do pacote lidR no processamento de dados LiDAR, consulte as seguintes referências: