LS_band_weighting: Distance-weighting of Landsat variable

In STBrinkmann/DRIGLUCoSE: DRI-GLUCoSE

Distance-weighting of Landsat variable

Description

To account for the diminishing effect of spatial variables (e.g. NDVI exposure) as distance increases, we fitted a logit function to weight each incremental isochrone.

Usage

LS_band_weighting(
  isochrones,
  tag = "tag",
  time = "time",
  landsat_list,
  band = "NDVI",
  b = 8,
  m = 0.5,
  stats = list("sd", "median", list("percentile", 0.05), list("percentile", 0.95),
    "skew"),
  cores = 1
)

Arguments

isochrones	object of class sf containing road features, derived from the isodistances function.
tag	character or NA; string containing the column name, that indicates the unique tag column.
time	character; string containing the column name, that indicates the time column.
landsat_list	list of landsat products, derived from the LS_L1C function. See Details.
band	character; spectral band or index from the landsat object, to which distance-weighting should be applied.
b	numeric; slope of distance decay function
m	numeric; x for g(x) = 0.5
stats	The function to be applied. See Details
cores	the number of cores to use.

Details

All unique tags from the isochrones shapefile will be overlapped with the landsat_list, and only the first overlapping landsat image will be used. This allows you to analyze multiple areas at once.

Supported functions for the stat object:
"sum", "mean", "min", "max", "sd", "rms", "skew", "median" and "percentile". When using "percentile", provide a list as follows:
list("percentile", upper_lim), where upper_lim is the upper limit of the percentile (e.g. 0.95).

For a detailed explanation of this method, see documentation on GitHub.

Value

A tibble containing the spatially weighted statistics of the intersects of the isochrone and Landsat objects.

STBrinkmann/DRIGLUCoSE documentation built on Oct. 14, 2022, 11:26 a.m.