gfGapLength: Julendat function to identify gap lengths in measurement...

View source: R/gfGapLength.R

gfGapLengthR Documentation

Julendat function to identify gap lengths in measurement series

Description

This is a function taken from Julendat to identify gap lengths in (eco-)climatological measurement series of a given parameter.

Usage

gfGapLength(
  data.dep,
  pos.na,
  gap.limit,
  end.datetime = Sys.Date(),
  units = "days",
  ...
)

Arguments

data.dep

Object of class ki.data. See as.ki.data, gsod2ki for further information.

pos.na

Numeric. Indices of missing data points.

gap.limit

Numeric. Maximum length of a measurement gap. All gaps exceeding this threshold will not be considered.

end.datetime

Object of class Date, default is Sys.Date(). Not required for GSOD data processing.

units

Character. Measurement interval, typically "days" for GSOD data.

...

Additional arguments. Currently not in use.

Value

A list containing start, end, and length of each data gap.

Author(s)

Florian Detsch

Examples

## Not run: 
moshi <- subset(gsodstations, `STATION NAME` == "MOSHI")

gsod_moshi <- dlGsodStations(usaf = moshi$USAF,
                             start_year = 1990, end_year = 1995,
                             dsn = tempdir(),
                             unzip = TRUE)

# Conversion to KiLi SP1 `ki.data` object
ki_moshi <- gsod2ki(data = gsod_moshi,
                    prm_col = c("TEMP", "MIN", "MAX"),
                    df2ki = TRUE)

# Identify length per data gap
gfGapLength(data.dep = ki_moshi, 
            pos.na = which(is.na(methods::slot(ki_moshi, "Parameter")$TEMP)), 
            gap.limit = 365,
            units = "days", 
            end.datetime = Sys.Date())

## End(Not run)


environmentalinformatics-marburg/GSODTools documentation built on Jan. 5, 2024, 12:19 a.m.