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/12/03 03:32] – 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>
lidar/projetos/rcode_invconve.1701574339.txt.gz · Última modificação: 2024/03/23 20:17 (edição externa)