Диаграмма классификации горных пород

Диаграмма классификации горных пород.

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

require('ggplot2');

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

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

  pts <- data.frame (Name = ds$Name,
                     X = ds$SiO2_wtp,
                     Y = ds$Na2O_wtp + ds$K2O_wtp);
 
  labels <- c("SiO2", "Na2O+K2O");
                  
  curve_1 <- data.frame(X1=c(45),
                        Y1=c(0),
                        X2=c(78),
                        Y2=c(4));
                        
  curve_2 <- data.frame(X1=c(39,53,64),
                        Y1=c(1,5,7.5),
                        X2=c(53,64,78),
                        Y2=c(5,7.5,8)
                        );
  
  
  curve_3 <- data.frame(X1=c(35.5,44,48,56.5,61,65,69,72),
                        Y1=c(1.5,4,5.5,11.5,12.5,12,10,9.5),
                        X2=c(44,48,56.5,61,65,69,72,78),
                        Y2=c(4,5.5,11.5,12.5,12,10,9.5,9));
                        
  curve_4st <- data.frame(X1=c(31,32),
                          Y1=c(5,5.8),
                          X2=c(32,33.2),
                          Y2=c(5.8,6.7));
                          
  curve_4 <- data.frame(X1=c(33.2,33.7,34.7,36,38,42,46,53,57,62,66,72),
                        Y1=c(6.7,9.5,12,14,16.5,19,20,19,17,13.5,12,11),
                        X2=c(33.7,34.7,36,38,42,46,53,57,62,66,72,78),
                        Y2=c(9.5,12,14,16.5,19,20,19,17,13.5,12,11,10));
  
  curve_5 <- data.frame(X1=c(54,59,64,66,67,69.5,71,75),
                        Y1=c(7.75,9,10,10.1,10,9,8.5,8.2),
                        X2=c(59,64,66,67,69.5,71,75,78),
                        Y2=c(9,10,10.1,10,9,8.5,8.2,8.1));
  
  ret <- c();
  
  brbline <- data.frame(X1=c(44,53,64),
                        Y1=c(0,0,0),
                        X2=c(44,53,64),
                        Y2=c(21,21,21));

  limits <- c(30,80,0,21);
  
  marks <- data.frame(
          X=c(37,48.5,58,71,35,39,42,41.5,40,36.5,39.5,46.5,47.5,48,47,51,50,55.5,57,61,61,66,56,60.5,56,61,59,56,60.75,70,66,73.5,73.5,68.5,71.5,57,66,70,76),
          Y=c(21,21,21,21,1,1.7,1.4,1,0.5,6.5,14.5,17.5,11.5,9.5,4.25,6.5,3,15.25,14.5,11.5,10.0,10.0,8.2,9.25,7.5,8.5,7.25,4.5,5.5,10.75,9,10.25,9.5,8.5,9,8.5,5.5,6,6),
          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граниты","Риолиты,\nлейкограниты")
      );

  ret$type <- "diagram";
  ret$points <- pts;
  ret$limits <- limits;
  ret$marks <- marks;
  ret$lines <- rbind(brbline, curve_1, curve_2, curve_3, curve_4st, curve_4, curve_5);
  ret$labels  <- labels;
  ret$markscale <- 0.5;
  
  return(ret);
}