plot.ves: plot.ves

Description Usage Arguments Author(s) See Also Examples

Description

Function to plot the Vertical Electric Sounding data. This function can create two different types of plots: resistivity and interpretation. The resistivity plot includes the apparent resistivity vs electrode spacing. In the interpretation plot, the values of the apparent resistivity vs electrode spacing and the thickness and real resisitivity of the layers are included.

Usage

1
2
3
4
## S3 method for class 'ves'
plot(x, main = NULL, type = c("ves", "transformation",
  "diagnostic"), trans.type = c("direct", "scaling", "zohdy",
  "zohdy.smoothed"), ...)

Arguments

x

A VES object

main

Title of the plot

type

A character string specifying the plot type. Currently only ves (measurements and earth model), transformation (resistivity-depth transformation) and diagnostic are supported.

trans.type

A character string specifying the transformation type. Only direct, scaling, zohdy and zohdy.smoothed are currently supported.

...

Additional parameters to be passed to the plot function

Author(s)

Oscar Garcia-Cabrejo khaors@gmail.com

See Also

Other base functions: print.ves, summary.ves, ves

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
library(grid)
library(gridExtra)
data(ves_data1)
ab2 <- ves_data1$ab2
apprho <- ves_data1$apprho
sev1a <- ves(id= "VES1", ab2 = ab2, apprho = apprho)
# Conventional VES plot
plot(sev1a, type = "ves")
# Apply direct transformation to sev1a
plot(sev1a, type = "transformation", trans.type = "direct")
# Apply scaling transformation to sev1a
plot(sev1a, type = "transformation", trans.type = "scaling")
# Apply Zohdy transformation to sev1a
plot(sev1a, type = "transformation", trans.type = "zohdy")
# Apply Zohdy smoothed transformation to sev1a
plot(sev1a, type = "transformation", trans.type = "scaling")
# Diagnostic plot
rho <- c(40,70,30, 20)
thick <- c(2,10,50,500)
par0 <- c(rho, thick)
res.nls <- calibrate_nls(sev1a, par0, iterations = 30, ireport = 5, trace = FALSE)
sev1a$rhopar <- res.nls$rho
sev1a$thickpar <- res.nls$thickness
sev1a$interpreted <- TRUE
p1 <- plot(sev1a, type = "ves")
print(p1)

khaors/rves documentation built on May 29, 2019, 3:07 a.m.