Calculate Human Footprint Tolerance Index (HFTI)

Description

This function calculates the HFTI for each species occurrence based on the average human footprint within a 1 km buffer around each occurrence.

Usage

calculate_HFTI(occurrences, footprint, divisor = 50, buffer_m = 1000)

Arguments

Value

A list with:

data

The sf object of occurrences with a new column HFTI_value.

HFTI

A data frame with the median HFTI per species.

Examples

# Note: simple example using the internal "toy" data
footprint <- load_human_footprint()
# We need dummy occurrence data for the example to run
occ_data <- data.frame(species = "SpA", longitude = -38.4, latitude = -12.6)
processed <- process_occurrences(occ_data)
hfti_result <- calculate_HFTI(processed, footprint, divisor = 50, buffer_m = 1000)

Calculate Road Tolerance Index (RTI)

Description

Esta função calcula o RTI para cada ocorrência de espécie com base na distância mínima até a estrada mais próxima. A fórmula usada é:

RTI = \frac{ref\_dist}{dist + ref\_dist}

Usage

calculate_RTI(occurrences, roads, ref_dist = 3.5)

Arguments

Details

Se dist = 0, então RTI = 1. Conforme dist aumenta, RTI aproxima-se de 0.

Value

Uma lista com:

data

O objeto sf de ocorrências, porém com uma nova coluna RTI_value com o índice calculado.

RTI

Um data frame com a mediana de RTI por espécie.

Examples

# Carregar rodovias de exemplo (Toy data)
roads <- load_roads()

# Criar dados de ocorrencia ficticios na area do Toy Data
occ_data <- data.frame(
  species = "SpA", 
  longitude = -38.4, 
  latitude = -12.6
)

# Processar ocorrencias
processed <- process_occurrences(occ_data)

# Calcular RTI
rti_res <- calculate_RTI(processed, roads, ref_dist = 3.5)

Compute both RTI and HFTI

Description

This main function processes occurrence data (removing duplicates, filtering close points), then calculates the Highway Tolerance Index (RTI) and the Human Footprint Tolerance Index (HFTI).

Usage

compute_indices(
  data,
  remove_duplicates = TRUE,
  buffer_km = 1,
  ref_dist = 3.5,
  divisor = 50,
  data_type = "toy"
)

Arguments

Value

A list containing:

RTI

A data frame with median RTI values per species.

HFTI

A data frame with median HFTI values per species.

indices

A data frame combining both RTI and HFTI by species.

processed_data

The sf object of occurrences with RTI and HFTI values.

Examples

# Example usage with internal toy data
# Using coordinates that fall within the provided toy dataset (Bahia, Brazil)
sample_data <- data.frame(
   species = c("SpA", "SpA", "SpB", "SpB"),
   longitude = c(-38.40, -38.41, -38.42, -38.43),
   latitude = c(-12.60, -12.61, -12.62, -12.63)
)

result <- compute_indices(
   data = sample_data,
   remove_duplicates = TRUE,
   buffer_km = 1,
   ref_dist = 3.5,
   divisor = 50
)

print(result$RTI)
print(result$HFTI)
print(result$indices)

Generate and save maps & density plots for RTI and HFTI of all species

Description

Esta funcao recebe o resultado de compute_indices(), um shapefile de area de estudo e o diretorio de saida. Para cada especie presente em result$processed_data, ela gera:

1. Mapa de RTI (points coloridos por RTI_value)

2. Mapa de HFTI (points coloridos por HFTI_value)

3. Grafico de densidade de RTI (com linha vertical na mediana)

4. Grafico de densidade de HFTI (com linha vertical na mediana) Alem disso, salva um arquivo CSV com RTI e HFTI medios por especie.

Usage

generate_all_reports(
  result,
  area_shapefile,
  out_dir,
  map_crs = 4326,
  point_size = 1.5,
  density_fill = "steelblue",
  density_alpha = 0.35
)

Arguments

Value

Invisivel. A funcao retorna NULL apos salvar todos os arquivos em out_dir.

Load Human Footprint GeoTIFF

Description

This function loads the human footprint data. By default, it loads a small "toy" raster included in the package. To use the full Brazil raster, set type = "full".

Usage

load_human_footprint(type = "toy")

Arguments

Value

A RasterLayer object containing human footprint data.

Load Brazilian Roads Data

Description

This function loads the roads data. By default, it loads a small "toy" dataset included in the package for testing and examples. To use the full dataset (all Brazilian roads), set type = "full". The full dataset will be downloaded from Zenodo upon the first request and cached locally.

Usage

load_roads(type = "toy")

Arguments

Value

An sf object containing roads data.

Process occurrence data

Description

This function processes species occurrence data by optionally removing duplicates and filtering out points that are closer than a specified distance (buffer).

Usage

process_occurrences(data, remove_duplicates = TRUE, buffer_km = 1)

Arguments

Value

An sf object with the processed occurrences.

Examples

# Example data
occ_data <- data.frame(
  species = c("SpA", "SpA", "SpB", "SpB"),
  longitude = c(-50, -50.0005, -51, -51.001),
  latitude = c(-10, -10.0005, -11, -11.001)
)

processed <- process_occurrences(
  data = occ_data,
  remove_duplicates = TRUE,
  buffer_km = 1
)
processed