st_integration_index: Calculate integration index

View source: R/st_integration_index.R

st_integration_indexR Documentation

Calculate integration index

Description

This function calculates an index defining the intergration of new land use area into conisting land use area of the same land use class conisdering two different time steps.

Usage

st_integration_index(
  tool = "sf",
  geom.old,
  geom.new,
  geom.boundary = NULL,
  tol = 0.1,
  precision = 0,
  env.rsaga = NULL,
  snap.rgrass = 1e-05,
  dist.new = 0,
  override = TRUE,
  return.geom = FALSE,
  quiet = TRUE
)

Arguments

tool

tool to perform erase-vector-operation. Default is "sf". However, for complex polygons "grass" is highly recommended. Fot the use of "grass" a valid GRASS GIS-session mus be initiated. "saga" is also supported, then env.rsaga must be properly set

geom.old

object of class sf representing a land use class of a previous time step

geom.new

object of class sf representing a land use class of a following time step

geom.boundary

polygon of class sf representing subregions, e.g. administrative boundaries

tol

tolerance value for overlapping area meter square

precision

precision for process using sf::st_set_precision. See st_precision. Default: 0

env.rsaga

environment of SAGA GIS. If st_erase fails then SAGA GIS erase is used. Default: NULL, but in function call if not set: rsaga.env

snap.rgrass

similar to sf::st_set_precision. Snapping threshold for boundaries. Default: 1e-04

dist.new

buffer distance in mapping units around new geometry (result of erase) using sf::st_buffer. Default: 0

override

If TRUE projection of geometry-operation is replaced by projection of geom.old. This flag prevents errors when the original datum is not found by the tool. Dafault: FALSE

return.geom

If set to TRUE, intermediate geometries are returned as well. Default: FALSE

quiet

If set to FALSE actual state is printed to console. Default: TRUE.

Value

Vector with integration index (completly integrated: 2/3 < R < 1; good integrated: 1/3 < R < 2/3; low integrated: 0 < R < 1/3; not integrated: 0)

Note

If the unique border line has incomprehensible gabs(can happen for "complicated" geometries), it can help to use a small buffer distance (< 0.05, e.g. 1e-04) by dist.new to get a complete edge line Code is based on the following references:

  • Meinel, G., & Winkler, M. (2002). Spatial analysis of settlement and open land trends in urban areas on basis of RS data œ studies of five European cities over a 50-year period.

  • Siedentop, S., Heiland, S., Lehmann, I., & Schauerte-Lüke, N. (2007). Regionale Schlüsselindikatoren nachhaltiger Flächennutzung für die Fortschrittsberichte der Nationalen Nachhaltigkeitsstrategie–Flächenziele (Nachhaltigkeitsbarometer Fläche). Abschlussbericht, BBR-Forschungen, H, 130.

Depending on the selected tool, the result of the vector-operation can differ significantly!


raff-k/VLSM documentation built on Oct. 13, 2023, 11:13 a.m.