README (original) (raw)

ssebiEF Package

The ssebiEF package provides tools to calculate the Simplified Surface Energy Balance Index (SSEBI) and the Evaporative Fraction (EF) from raster data. It utilizes albedo and surface-air temperature difference (TS-TA) to generate the SSEBI or EF raster.

Installation

You can install the development version of ssebiEFdirectly from your local folder:

Requirements

• The input rasters (TS and TA) must have the same unit of measurement (Kelvin preferably) to ensure accurate calculations of TS-TA and Simplified Surface Energy Balance Index (SSEBI) or the Evaporative Fraction (EF).
• The input rasters (albedo and TS-TA) must:
  • Have the same Coordinate Reference System (CRS).
  • Cover the same geographic extent.
• If these conditions are not met, the function attempt to align the rasters automatically (with warnings).

Dependencies

The ssebiEF package automatically manages the following dependencies: - terra for raster manipulation. -dplyr for data manipulation. - stats for statistical calculations.

You do not need to load these packages manually; they are handled internally by ssebiEF

Examples Usage of thessebiEF Package

To use the ssebiEF package, start by loading the necessary library in R.

• Define the paths to the input files (albedo and TS-TA rasters) and the output path for the calculated SSEBI raster.
  # Load the necessary package library(ssebEF)
  # Input raster paths albedo_path <- “path/to/albedo.tif” TS_TA_path <- “path/to/TS_TA.tif”
  # Output raster path output_path <- “path/to/SSEBI.tif”
  # Calculate SSEBI calculate_SSEBI( albedo_path = albedo_path, TS_TA_path = TS_TA_path, output_path = output_path, n_intervals = 20, percentile = 0.01 )
  # Check the output SSEBI <- terra::rast(output_path) plot(SSEBI)
• Define the paths to the input files (albedo and TS-TA rasters) and the output path for the calculated EF raster.
  # Load the necessary package library(ssebiEF)
  # Input raster paths albedo_path <- “path/to/albedo.tif” TS_TA_path <- “path/to/TS_TA.tif”
  # Output raster path output_path <- “path/to/EF.tif”
  # Calculate EF calculate_EF( albedo_path = albedo_path, TS_TA_path = TS_TA_path, output_path = output_path, n_intervals = 20, percentile = 0.01 )
  # Check the output EF <- terra::rast(output_path) plot(EF)