Learn R Programming

Luminescence (version 0.2.4)

CW2pLM: 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

values
data.frame (required): measured curve data of type stimulation time (t) (values[,1]) and measured counts (cts) (values[,2])

Value

  • Returns a data.frame of the same form as the input of the transformed curve data.

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=\sqrt{(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)

References

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

See Also

CW2pLM The output of the function can be further used for LM-OSL fitting:

CW2pLMi, CW2pHMi, CW2pPMi, fit_LMCurve

Examples

Run this code
##read curve from BINfileData transform curve an plot values
data(ExampleData.BINfileData)

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

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

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

##plot
plot(values.transformed)

Run the code above in your browser using DataLab