TAS диаграмма для плутонических пород. Сокращения: н/щ - низкощелочной, у/щ - умереннощелочной, щел. - щелочной, кв. - кварцевый.
| 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граниты",
"ультраосновные",
"основные",
"средние",
"кислые",
"н/щ",
"норм.\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);
}