MgO-FeO-Al2O3: Crust type, Pearce

Crust type

Использовать Проверить

require('ggtern');

classifydiag.crust_pearce <- function (wtp_data, ...)
{  
ds <- selectNames(wtp_data,c('Name','MgO_wtp', 'FeO_wtp', 'Fe2O3_wtp', 'Al2O3_wtp'))

pts <- data.frame(Name=ds$Name,
                  X = ds$MgO_wtp,
                  Y = ds$FeO_wtp+ds$Fe2O3_wtp,
                  Z = ds$Al2O3_wtp);

labels <- c("MgO", "FeO*", "Al2O3");
 
lines <- data.frame(X1 = c(0.38, 0.14, 0.21, 0.28, 0.21, 0.18, 0.17, 0.15, 0.13, 0.12, 0.15, 0.10),
                    Y1 = c(0.32, 0.51, 0.27, 0.15, 0.27, 0.32, 0.34, 0.35, 0.33, 0.29, 0.35, 0.44),
                    Z1 = c(0.30, 0.35, 0.52, 0.57, 0.52, 0.50, 0.50, 0.52, 0.54, 0.59, 0.52, 0.46),
                    X2 = c(0.23, 0.23, 0.23, 0.21, 0.18, 0.17, 0.15, 0.13, 0.12, 0.10, 0.10, 0.08),
                    Y2 = c(0.33, 0.33, 0.33, 0.27, 0.32, 0.34, 0.35, 0.33, 0.29, 0.21, 0.44, 0.49),
                    Z2 = c(0.46, 0.46, 0.46, 0.52, 0.50, 0.50, 0.52, 0.54, 0.59, 0.69, 0.46, 0.43));
 
marks  <- data.frame(Label=c("Ocean Island", "Ocean Ridge", "Con.", "S.C.I.", "Orogenic"),
                     X=c(0.38, 0.37, 0.10, 0.05, 0.12),
                     Y=c(0.40, 0.25, 0.50, 0.38, 0.18),
                     Z=c(0.22, 0.39, 0.40, 0.58, 0.70));
ret <- c();
ret$type <- "ternogram";
ret$points <- pts;
ret$labels  <- labels;
ret$lines <- lines;
ret$marks <- marks;
 
return(ret);
}