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

TAS диаграмма для плутонических пород. Сокращения: н/щ - низкощелочной, у/щ - умереннощелочной, щел. - щелочной, кв. - кварцевый.

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

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("фойдовые\nпороды",
                             "щелочные\nпороды разной\nосновности",
                             "монцо-\nгаббро",
                             "габбро",
                             "н/щ\nгаббро",
                             "монцо-\nгаббро\nдиориты",
                             "габбро\nдиориты",
                             "н/щ\nгаббро\nдиориты",
                             "монцо-\nдиориты",
                             "диориты",
                             "н/щ\nдиориты",
                             "монцониты",
                             "кв. диориты",
                             "н/щ кв.\nдиориты",
                             "фоидовые\nсиениты",
                             "щел.\nсиениты",
                             "сиениты",
                             "грано-\nсиениты",
                             "грано-\nдиориты",
                             "н/щ грано-\nдиориты",
                             "щел.\nграниты",
                             "у/щ\nграниты",
                             "граниты",
                             "н/щ\nграниты",
                             "щел.\nлейко-\nграниты",
                             "у/щ\nлейко-\nграниты",
                             "лейко-\nграниты",
                             "н/щ\nлейко-\nграниты",
                             "ультраосновные",
                             "основные",
                             "средние",
                             "кислые",
                             "н/щ",
                             "норм.\nщел.",
                             "у/щ",
                             "щел."
                             ),
                     X=c( 42,47.5,47.5, 48.5,  49, 50.8, 53.5, 53.5,  54,  57,  57 ,  58, 61.5, 61.5, 57.5, 61.5,   63,  65, 66.5, 66.5, 69, 70.5,  71,  71,   75,  75 , 75.5, 75.5, 40, 48.5, 57.5, 70, 79, 79, 79, 79),
                     Y=c(9.5, 8.5, 4.7,  2.8, 1.2,   6.3,    4,    2, 7.4, 4.8, 2.5 , 8.5,  5.5,    3, 13.5, 12.2, 10.3, 8.5,  6.5,  4.1, 11,    9, 6.8, 4.5, 10.7, 8.7,  7,    5, -0.5, -0.5, -0.5, -0.5, 5, 7, 8.5, 9.5));

limits <- c(35,79,-1,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,pls);
ret$labels  <- labels;
ret$marks <- marks;

return(ret);
}