Различия
Здесь показаны различия между двумя версиями данной страницы.
classifydiag:tas_ru_plutonic [2016/12/19 14:29] f0ma |
classifydiag:tas_ru_plutonic [2018/11/12 15:07] |
||
---|---|---|---|
Строка 1: | Строка 1: | ||
- | ====== TAS диаграмма для плутонических пород ====== | ||
- | TAS диаграмма для плутонических пород. | ||
- | |||
- | === Details === | ||
- | | Reference | <ref:bogatikov2009> | | ||
- | | Data | bulk | | ||
- | | Input | K2O_wtp, Na2O_wtp, SiO2_wtp | | ||
- | |||
- | <use:classifydiag.tas_ru_plutonic> | ||
- | <test:classifydiag.tas_ru_plutonic> | ||
- | |||
- | === Код метода === | ||
- | |||
- | <Code linenums lang-r> | ||
- | 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.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)); | ||
- | |||
- | 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); | ||
- | } | ||
- | </Code> | ||
- | |||
- | {{tag>classifydiag wholerock}} |