Planeja Inventários Florestais

planejar() cria um objeto de planejamento de inventário, com informações úteis sobre a suficiência amostral e disposição automática de parcelas em campo.

planejar(
  InvObj,
  processo,
  areapar,
  areatotal,
  pardims = NULL,
  media,
  desvpad,
  cv,
  errotol = 0.1,
  conf = 0.95,
  gerar_parcelas = NULL,
  bounds = NULL,
  seed = 1,
  plottype = c("point", "polygon"),
  complementar = FALSE,
  safe_bounds = TRUE,
  opt_crs = NULL,
  scale_ha = NULL,
  simn = NULL,
  ...
)

Arguments

InvObj: Objeto da classe Inv ou herdeiras.
processo: Caractere contendo o tipo de processo de amostragem. Veja detalhe para opções.
areapar: Numérico contendo a área da parcela em hectares.
areatotal: Numérico contendo a área total do projeto em hectares.
pardims: Vetor de tamanho 2 contendo as dimensões X e Y da parcela, em metros.
media: Numérico contendo a média a ser utilizada para suficiência amostral.
desvpad: Numérico contendo o desvio padrão amostral para cálculo da suficiência amostral.
cv: Numérico contendo o coeficiente de variação em porcentagem para cálculo da suficiência amostral.
errotol: Numérico contendo o erro máximo tolerado em porcentagem para cálculo da suficiência amostral. Padrão é 0.1 (10%).
conf: Numérico contendo o nível de confiança adotado. Padrão é 0.95 (95%).
gerar_parcelas: Lógico (TRUE/FALSE) ou NULL informando se as parcelas devem ser geradas. Veja detalhes.
bounds: Ou caractere com o caminho até arquivo, ou objeto espacial com a delimitação da área para geração de parcelas.
seed: Semente utilizada para sorteio de parcelas. Padrão é 1.
plottype: Caractere informando o tipo de parcela a ser gerada. Opções são 'ponto' ou "polygon".
complementar: Lógico (TRUE/FALSE) informando se apenas parcelas complementares para suficiência devem ser geradas. Usada apenas se InvObj for informado.
safe_bounds: Lógico (TRUE/FALSE) informando se uma margem de segurança deve ser criada para evitar parcelas que extrapolem os limites da área.
opt_crs: Numérico opcional contendo o código epsg do sistema de coordenadas final das parcelas.
scale_ha: Lógico (TRUE/FALSE) informando se media, desvpad etc. estão por hectare. Por padrão se presume que sim.
simn: Numérico opcional que fixa um número de parcelas para cálculo. Eventualmente útil.
...: Demais argumentos corretamente alocados nos métodos de cada processo de amostragem.

Value

Retorna um objeto da família de classes PlanInv*.

Details

planejar() possui dois componentes. Um deles é de processamento, que retorna informações de suficiência amostral. O outro é um componente espacial que trata da geração de parcelas para uso imediato. Esse segundo componente não é obrigatoriamente realizado, pois depende de existir alguma delimitação da floresta/projeto para que se possa lançar as parcelas de acordo com o processo amostral adotado.

A função pode ser utilizada a partir de duas fontes. A primeira é a partir de um objeto Inv ou herdeiro desse, passado para o argumento InvObj=. Se isso ocorrer, presume-se que InvObj= é um inventário piloto, e todas as informações disponíveis em InvObj serão utilizadas por padrão. Isto é, media, desvpad, cv, processo, areapar, areatotal e bounds serão obtidos a partir de InvObj, desde que tais elementos existam no objeto. Se mesmo assim o usuário informar processo, media ou algum já contido em Inv, então esse argumento irá sobrepor o obtido por InvObj.

errotol é o erro relativo tolerado para o inventário. Se for passado um número entre 0 e 1, presume-se que está na escala relativa (e.g. 0.1). Se for um número maior que 1, presume-se que é um valor percentual (e.g. 10%).

Além da área da parcela, também é possível informar as dimensões da parcela em metros com pardims=. Nesse caso, pardims= recebe um vetor de tamanho 2 com os tamanhos do lado de uma parcela retangular. Se for passado o argumento pardims=, então a área da parcela será calculada a partir desse, e areapar será ignorado. Por padrão, pardims= resulta na formação de polígonos das parcelas, independente do que for especificado em plottype=.

plottype= indica o tipo de objeto espacial a ser gerado. As opções são pontos ("point") ou polígonos ("polygon"). Utilizado apenas se o componente espacial da função for executado. Veja que caso pardims= seja especificado, então o informado em plottype= é ignorado e obrigatoriamente serão gerados polígonos em vez do ponto central. Veja ainda que se plottype="polygon" e pardims= for omitido, então será assumida uma parcela quadrada para a área de parcela informada em areapar= ou obtida via InvObj.

gerar_parcelas= tem como padrão NULL, embora receba TRUE/FALSE. Isso ocorre pois caso haja possibilidade de criar as parcelas, a função irá tentará executar o componente espacial. Entretanto, informar TRUE/FALSE irá forçar a construção ou não de parcelas, sujeito a falhas caso não seja possível.

bounds= é a região delimitada onde serão alocadas as parcelas, portanto é elemento crucial para o componente espacial. Se bounds= receber um caractere, será interpretado como o caminho para um arquivo shapefile, como "path/to/file.shp", e será carregado utilizando o pacote {sf}.bounds= também pode receber diretamente um objeto de classe sf ou sfc. Por fim, bounds= pode ser omitido se InvObj já possuir bordadura.

complementar= é utilizada apenas se InvObj for utilizado. Nesse caso, será entendido que InvObj é um inventário piloto, e se complementar=TRUE, então apenas as parcelas que faltam para atingir suficiência amostral serão geradas. Caso um inventário piloto já possua suficiência amostral, então nenhuma parcela será gerada. Utilize complementar=FALSE para ignorar e gerar todas as parcelas.

safe_bounds= cria uma margem de segurança na borda de bounds=, para garantir que as parcelas geradas estejam totalmente contidas na área de interesse. Tal margem é construída independente do tipo de parcela gerada (ponto ou polígono). A margem é calculada como metade do tamanho da diagonal de uma parcela retangular.

opt_crs= recebe o código EPSG de um sistema de coordenadas, caso o usuário queira que as parcelas geradas sejam de algum sistema específico. Todos os objetos espaciais serão transformados para o sistema informado.

scale_ha por padrão é TRUE, e apenas indica se os valores de media, desvpad etc. estão por hectare ou não. Se for omitido e InvObj for informado, então o valor de InvObj será utilizado.

simn recebe um inteiro e fixa o número de parcelas utilizado. Útil para fazer simulações ou testar diferentes números de parcelas e comparar com o erro de amostragem obtido.

A função planejar() retorna um objeto de classe PlanInv*, tal que se for um processo aleatório simples, então será um objeto PlanInvAAS, etc. Cada classe possui métodos personalizados para cada processo de amostragem realizado.

Author

Gabriel Agostini Orso