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("фойдовые\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),
                     Y=c(9.5, 8.5, 4.7,  2.8, 1.2,   6.1,    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.3, 8.7,  6.5,    5));
 
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,pls);
ret$labels  <- labels;
ret$marks <- marks;
 
return(ret);
}