Ciências Florestais

Ferramentas do usuário

Ferramentas do site

Visitou:
lidar:projetos:rcode_invconve

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

Ambos lados da revisão anteriorRevisão anterior
Próxima revisão		Revisão anterior
lidar:projetos:rcode_invconve [2023/11/26 12:48] – lcerlidar:projetos:rcode_invconve [2024/03/23 20:17] (atual) – edição externa 127.0.0.1
Linha 8: Linha 8:
 # Autor: Luiz Carlos Estraviz Rodriguez # Autor: Luiz Carlos Estraviz Rodriguez
 #        Departamento de Ciências Florestais #        Departamento de Ciências Florestais
-#        ESALQ/USP - 25/Nov/2023+#        ESALQ/USP - 16/Jan/2024
 # #
 # Inventário florestal da Fazenda Modelo # Inventário florestal da Fazenda Modelo
Linha 19: Linha 19:
 # Linguagem de programação: # Linguagem de programação:
 #       R (v 4.3) #       R (v 4.3)
-#       pacote lidR* (v 4.0.3 March 11th, 2023) Jean Romain Roussel 
 # #
 # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Linha 86: Linha 85:
 parcelas <- parcelas <-
   parcelasSemGeo %>%   parcelasSemGeo %>%
-  select(SUBTALHAO, CHAVE2, DATAREALIZ, IDINV, AREAPARCEL, MHDOM, VTCC, AB) %>% +  select(SUBTALHAO, CHAVE2, DATAREALIZ, IDINV, AREAPARCEL, MHDOM, VTCC) %>% 
   arrange(SUBTALHAO) %>% as.data.frame   arrange(SUBTALHAO) %>% as.data.frame
  
Linha 134: Linha 133:
 # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 dirNome <- paste0('C:/LiDAR/', prjNome) dirNome <- paste0('C:/LiDAR/', prjNome)
-  dir.create(dirNome, showWarnings = F)+dir.create(dirNome, showWarnings = F)
 dirNome <- paste0(dirNome, '/GRAPH/') dirNome <- paste0(dirNome, '/GRAPH/')
-  dir.create(dirNome, showWarnings = F)+dir.create(dirNome, showWarnings = F)
  
 # Cria mapa dos talhões com localização das parcelas (EPSG: 31983) # Cria mapa dos talhões com localização das parcelas (EPSG: 31983)
Linha 184: Linha 183:
  
 # Cria lista com variáveis de interesse # Cria lista com variáveis de interesse
-variaveisDeInteresse <- c("MHDOM", "VTCC", "AB")+variaveisDeInteresse <- c("MHDOM", "VTCC")
 # Cria dataframe com valores observados das variáveis de interesse # Cria dataframe com valores observados das variáveis de interesse
 df <- parcelas %>% select(all_of(variaveisDeInteresse)) df <- parcelas %>% select(all_of(variaveisDeInteresse))
Linha 319: Linha 318:
   B    = errDesired * mean(y)   B    = errDesired * mean(y)
   n    = sum( Nh^2 * vars / Wh ) /    n    = sum( Nh^2 * vars / Wh ) / 
-         ((B^2 * (sum(Nh)^2))/4 + sum(Nh * vars))+    ((B^2 * (sum(Nh)^2))/4 + sum(Nh * vars))
   return(n)   return(n)
 } }
Linha 409: Linha 408:
 # save.image(file='Modelo_inventario.RData')       # faz backup dos dados # save.image(file='Modelo_inventario.RData')       # faz backup dos dados
 # load('Modelo_inventario.RData')                       # recupera backup # load('Modelo_inventario.RData')                       # recupera backup
- 
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
-# Este script foi atualizado em Novembro/2023 depois de instalar, 
-# na seguinte sequência, as mais recentes versões do pacote lidR*: 
-# 
-#   devtools::install_github("Jean-Romain/rlas", dependencies=TRUE) 
-#   devtools::install_github("Jean-Romain/lidR") 
-# 
-#   * lidR latest development version 
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
 #  ---- #  ----
 </code> </code>
Linha 427: Linha 416:
   * [[https://r-lidar.github.io/lidRbook/ | The lidR package (eBook)]]   * [[https://r-lidar.github.io/lidRbook/ | The lidR package (eBook)]]
   * [[https://hal.inrae.fr/hal-02972284/file/roussel_2020.pdf | Roussel et al.(2020) lidR: An R package for analysis of ALS data]]   * [[https://hal.inrae.fr/hal-02972284/file/roussel_2020.pdf | Roussel et al.(2020) lidR: An R package for analysis of ALS data]]
 +  * [[https://github.com/r-lidar/lidR/wiki/Area-based-approach-from-A-to-Z | Inventário passo a passo com lidR]] 
   * [[https://www.rdocumentation.org/packages/lidR/versions/4.0.1 | RDocumentation (lidR)]]   * [[https://www.rdocumentation.org/packages/lidR/versions/4.0.1 | RDocumentation (lidR)]]
   * [[https://cran.rstudio.com/web/packages/lidR/vignettes/ | Vignettes]]   * [[https://cran.rstudio.com/web/packages/lidR/vignettes/ | Vignettes]]
 +
 +e sobre os pacotes //terra// e //sf// para processamento de dados espaciais:
 +  * [[https://rspatial.org/terra/pkg/index.html | The terra package]]
 +  * [[https://r-spatial.github.io/sf/articles/sf1.html | Simple Features for R ]]
lidar/projetos/rcode_invconve.1701002927.txt.gz · Última modificação: 2024/03/23 20:17 (edição externa)