New functionalities

Feel free to contact me if you have any troubles, questions, comments, wishes, etc.

emanuel.huber@alumni.ethz.ch

Check the companion website for tutorials:

http://emanuelhuber.github.io/RGPR/

Time-zero, first wave break, etc.

data(frenkeLine00)
tfb <- firstBreak(frenkeLine00)
t0 <- firstBreakToTime0(tfb, frenkeLine00, c0 = 0.299)
time0(frenkeLine00) <- t0
frenkeLine00_2 <- time0Cor(frenkeLine00, method = "pchip")

Recent functions for the class GPRsurvey

1. create an object of the class GPRsurvey

SU <- GPRsurvey(LINES)

2. the function linkCoordFid() relates the fiducials markers to the recorded trace positions (x,y,z)

fids <- linkCoordFid(y = SU, xyz = TO[, c("E", "N", "Z")], pcode = TO$PCODE)
which(sapply(fids, is.null))     # test it worked well

3. interpolate trace positions

SU <- interpPos( SU , fids , r = r)
plot(SU)

4. reverse the GPR lines s.t. they are all oriented in the same direction

SUr <- reverse(SU)
plot(SUr)

you can also only reverse some specific lines

SUr <- reverse(SU, id = c("YLINE01" , "YLINE30"))
SUr <- reverse(SU, id = seq(1, 30, by = 2))

5. rotate survey

SUrot <- rotate(SUr, alpha = pi/75)
plot(SUrot)

6. sequentially apply processing step on an object of the class GPRsurvey

   Create a list of lists defining the processing steps. The name of each sublist (e.g., agc) correspond to the name of the processing function, each sublist corresponds to the argument = values as defined by the processing function.

prc <- list("gain"    = list(type = "agc", w = 10),
            "dewow"   = list(type = "Gaussian", w = 15),
            "fFilter" = list(f = c(10, 40, 150, 200), type = "bandpass",
                             plotSpec = FALSE))

Apply all the processing step with the function papply()

SUr1 <- RGPR::papply(SUr, prc)