RLum.OSL_decomposition {OSLdecomposition} | R Documentation |
Separate CW-OSL components in RLum.Analysis data sets
Description
Calculates the CW-OSL signal component intensities for each CW-OSL measurement under the requirement that the decay rates are already given. The signal decomposition process uses an analytical approach described in detail in Mittelstrass (2019) and Mittelstrass et al. (in preparation). This function processes RLum.Analysis data sets created within the Luminescence-package (Kreutzer et al. 2012).
Usage
RLum.OSL_decomposition(
object,
record_type = "OSL",
K = 3,
decay_rates = NULL,
report = FALSE,
report_dir = NULL,
image_format = "pdf",
open_report = TRUE,
rmd_path = NULL,
verbose = TRUE
)
Arguments
object |
RLum.Analysis or list of RLum.Analysis (required):
Data set of one or multiple CW-OSL measured aliquots. The data set must either
contain a list element |
record_type |
character (with default):
Type of records, selected by the RLum.Analysis attribute |
K |
numeric (with default):
Number of components. Selects the according result table in the |
decay_rates |
numeric vector or data.frame (optional):
User-defined component decay rates. If this parameter is defined, the parameter |
report |
logical (with default):
Creates a |
report_dir |
character (optional):
Path of output directory if |
image_format |
character (with default):
Image format of the automatically saved graphs if |
open_report |
logical (with default):
If set to |
rmd_path |
character (with default): For advanced users: File path to the rmarkdown source code file of the report. This allows to execute a manipulated version of the report. |
verbose |
logical (with default): Enables console text output. |
Details
The workflow of this function is as follows:
-
optimise_OSLintervals: Approximates the optimal integration intervals. Uses the global average curve as time axis template. If none global average curve is given, one is automatically created using sum_OSLcurves.
-
decompose_OSLcurve: Calculates component intensities for all
record_type
measurements. Uses the"det"
algorithm if a background correction was performed with RLum.OSL_correction or the"det+nls"
algorithm if no background correction was performed. For error estimation, the"empiric"
approach is used. Creates a
html
report to summarize the results (optional).
Data sets must be formatted as RLum.Analysis objects and should have been processed with RLum.OSL_correction and RLum.OSL_global_fitting beforehand. Output objects are also RLum.Analysis objects and are meant for equivalent dose determination with Luminescence::analyse_SAR.CWOSL.
If report = TRUE
, a html
report of the results is rendered by the rmarkdown-package
and saved in the working directory, which is usually the directory of the data file.
This report can be displayed, shared and published online without any requirements regarding
the operation system or installed software. However, an internet connection is needed to display
the MathJax encoded equations and special characters.
The Rmarkdown source code of the report can be found with the following command:
system.file("rmd", "report_Step2.Rmd", package = "OSLdecomposition")
Value
The input object
, a list of RLum.Analysis objects is returned but with
a new list element object[["DECOMPOSITION"]]
, containing:
-
$decompositon.input
data.frame: Set of input components. Relevant is just the column$lambda
-
$results
data.frame: Overview table of decomposition -
$parameters
list: Input and algorithm parameters
The RLum.Data.Curve attribute @info
of each CW-OSL record contains the
new entry $COMPONENTS
with the curve-individual signal component parameters.
It can be read for example by:
object[[i]]@records[[j]]@info[["COMPONENTS"]]
Last updates
2022-05-02, DM: Added new parameter open_report
to give control over automatic browser opening
Author(s)
Dirk Mittelstrass, dirk.mittelstrass@luminescence.de
Please cite the package the following way:
Mittelstraß, D., Schmidt, C., Beyer, J., Heitmann, J. and Straessner, A.: R package OSLdecomposition: Automated identification and separation of quartz CW-OSL signal components, in preparation.
References
Kreutzer, S., Schmidt, C., Fuchs, M.C., Dietze, M., Fischer, M., Fuchs, M., 2012. Introducing an R package for luminescence dating analysis. Ancient TL, 30 (1), 1-8.
Mittelstraß, D., 2019. Decomposition of weak optically stimulated luminescence signals and its application in retrospective dosimetry at quartz (Master thesis). TU Dresden, Dresden.
See Also
RLum.OSL_global_fitting, decompose_OSLcurve, optimise_OSLintervals, Luminescence::analyse_SAR.CWOSL
Examples
#'FB_10Gy' is a dose recovery test with the Fontainebleau quartz
# measured in a lexsyg research with green LED stimulation
data_path <- system.file("examples", "FB_10Gy_SAR.bin", package = "OSLdecomposition")
data_set <- Luminescence::read_BIN2R(data_path, fastForward = TRUE)
# Separate components
data_set_decomposed <- RLum.OSL_decomposition(
data_set, decay_rates = c(0.8, 0.05))