Transform a CW-OSL curve into a pLM-OSL curve

Description

Transforms a conventionally measured continuous-wave (CW) curve into a pseudo linearly modulated (pLM) curve using the equations given in Bulur (2000).

Usage

CW2pLM(values)

Arguments

Details

According to Bulur (2000) the curve data are transformed by introducing two new parameters P (stimulation period) and u (transformed time):

P=2*max(t)

u=√{(2*t*P)}

The new count values are then calculated by

ctsNEW = cts(u/P)

and the returned data.frame is produced by: data.frame(u,ctsNEW)

Value

The function returns the same data type as the input data type with the transformed curve values.

Function version

0.4.1 (2015-11-29 17:27:48)

Note

The transformation is recommended for curves recorded with a channel resolution of at least 0.05 s/channel.

Author(s)

Sebastian Kreutzer, IRAMAT-CRP2A, Universite Bordeaux Montaigne (France)
R Luminescence Package Team

References

Bulur, E., 2000. A simple transformation for converting CW-OSL curves to LM-OSL curves. Radiation Measurements, 32, 141-145.

Further Reading

Bulur, E., 1996. An Alternative Technique For Optically Stimulated Luminescence (OSL) Experiment. Radiation Measurements, 26, 701-709.

See Also

CW2pHMi, CW2pLMi, CW2pPMi, fit_LMCurve, lm, RLum.Data.Curve

The output of the function can be further used for LM-OSL fitting: CW2pLMi, CW2pHMi, CW2pPMi, fit_LMCurve, RLum.Data.Curve, plot_RLum

Examples

##read curve from CWOSL.SAR.Data transform curve and plot values
data(ExampleData.BINfileData, envir = environment())

##read id for the 1st OSL curve
id.OSL <- CWOSL.SAR.Data@METADATA[CWOSL.SAR.Data@METADATA[,"LTYPE"] == "OSL","ID"]

##produce x and y (time and count data for the data set)
x<-seq(CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"]/CWOSL.SAR.Data@METADATA[id.OSL[1],"NPOINTS"],
       CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"],
       by = CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"]/CWOSL.SAR.Data@METADATA[id.OSL[1],"NPOINTS"])
y <- unlist(CWOSL.SAR.Data@DATA[id.OSL[1]])
values <- data.frame(x,y)

##transform values
values.transformed <- CW2pLM(values)

##plot
plot(values.transformed)

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(Luminescence)
Welcome to the R package Luminescence version 0.6.0 [Built: 2016-05-30 16:47:30 UTC]
Push Pin: 'Made of used sample carriers'
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Luminescence/CW2pLM.Rd_%03d_medium.png", width=480, height=480)
> ### Name: CW2pLM
> ### Title: Transform a CW-OSL curve into a pLM-OSL curve
> ### Aliases: CW2pLM
> ### Keywords: manip
> 
> ### ** Examples
> 
> 
> 
> ##read curve from CWOSL.SAR.Data transform curve and plot values
> data(ExampleData.BINfileData, envir = environment())
> 
> ##read id for the 1st OSL curve
> id.OSL <- CWOSL.SAR.Data@METADATA[CWOSL.SAR.Data@METADATA[,"LTYPE"] == "OSL","ID"]
> 
> ##produce x and y (time and count data for the data set)
> x<-seq(CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"]/CWOSL.SAR.Data@METADATA[id.OSL[1],"NPOINTS"],
+        CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"],
+        by = CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"]/CWOSL.SAR.Data@METADATA[id.OSL[1],"NPOINTS"])
> y <- unlist(CWOSL.SAR.Data@DATA[id.OSL[1]])
> values <- data.frame(x,y)
> 
> ##transform values
> values.transformed <- CW2pLM(values)
> 
> ##plot
> plot(values.transformed)
> 
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>