TAS диаграмма для плутонических пород

Это старая версия документа!

TAS диаграмма для плутонических пород.

Details

Reference	О. А. Богатиков, О. В. Петров, А. Ф. Морозов, Л. П. Шарпенок, В. В. Жданов, А. М. Курчавое, Е. А. Кухаренко, Э. А. Ланда, Б. А. Марковский, В. В. Шатов Петрографический кодекс России. Магматические, метаморфические, метасоматические, импактные образования, 2009 [ref]
Data	bulk
Input	K2O_wtp, Na2O_wtp, SiO2_wtp

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

Код метода

require('ggplot2');

classifydiag.tas_ru_plutonic <- function (wtp_data, ...){

ds <- selectNames(wtp_data, c("Name", "SiO2_wtp", "K2O_wtp", "Na2O_wtp"));

pts <- data.frame(Name=ds$Name,
                  X=ds$SiO2_wtp,
                  Y=ds$K2O_wtp + ds$Na2O_wtp);

marks  <- data.frame(Label=c(""),
                     X=c(0.0),
                     Y=c(0.0));

limits <- c(35,78,0,14);


sec   <- data.frame(X1 = c( 45,  52,  63),
                    Y1 = c(  0,   0,   0),
                    X2 = c( 45,  52,  63),
                    Y2 = c(0.5, 0.5, 0.5));

rads  <- data.frame(X1 = c(41,   41,  37,  35,  37,  41,  55,  69,  73,  45,  69, 73,   41,  45,  49, 61.25,   67, 72.5, 59.2),
                    Y1 = c( 0,    7, 0.8, 0.5, 0.8, 1.5, 5.8, 7.7, 7.9, 1.5, 5.7, 6,     0, 0.2, 0.2,   9.1,   10,   10, 10.5),
                    X2 = c(41, 52.5,  62,  37,  41,  55,  69,  73,  78,  69,  73, 78,   45,  49, 78,     67, 72.5,   78, 66.2),
                    Y2 = c( 7,   14,  14, 0.8, 1.5, 5.8, 7.7, 7.9,   8, 5.7,   6, 5.8, 0.2, 0.2, 4,      10,   10,    9, 12.2));

pls   <- data.frame(X1 = c( 45,  52, 52,  55,  55,  59,  59,  64,   64, 57.7,  69,    69  ,  73,  73, 72.5),
                    Y1 = c(0.2, 0.9,  5, 1.3,   6, 1.8, 6.4, 2.4,    7, 11.7,   3,  7.7 , 3.5, 7.9,   10),
                    X2 = c( 45,  52, 47,  55,  50,  59,  53,  64, 57.7, 52.5,  69,  65.5,  73,72.5, 71.5),
                    Y2 = c(  5,   5,  6,   6, 7.5, 6.4, 9.3,   7, 11.7,   14, 7.7,  12  , 7.9,  10, 11.7));

labels <- c("SiO2", "K2O+Na2O");

ret <- c();

ret$type <- "diagram";
ret$points <- pts;
ret$limits <- limits;
ret$lines <- rbind(rads,sec);
ret$labels  <- labels;

return(ret);
}

classifydiag, wholerock