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

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

require('ggplot2');

classifydiag.wholerock_type_rus <- function (wtp_data, annotate=T, labels=F, ...)
{  
  curve_1 <- data.frame(si=c(45,78), nak=c(0,4));
  curve_2 <- data.frame(si=c(39,53,64,78), nak=c(1,5,7.5,8));
  curve_3 <- data.frame(si=c(35.5,44,48,56.5,61,65,69,72,78), nak=c(1.5,4,5.5,11.5,12.5,12,10,9.5,9));
  curve_4st <- data.frame(si=c(31,32,33.2), nak=c(5,5.8,6.7));
  curve_4 <- data.frame(si=c(33.2,33.7,34.7,36,38,42,46,53,57,62,66,72,78), nak=c(6.7,9.5,12,14,16.5,19,20,19,17,13.5,12,11,10));
  curve_5 <- data.frame(si=c(54,59,64,66,67,69.5,71,75,78), nak=c(7.75,9,10,10.1,10,9,8.5,8.2,8.1));
  
  font_size=3;
  line_height=0.7;
  
  p <- ggplot(wtp_data, aes(x=SiO2_wtp,y=Na2O_wtp+K2O_wtp)) + geom_point() +
  geom_vline(xintercept = 44) + geom_vline(xintercept = 53) + geom_vline(xintercept = 64) +
  geom_line(data=data.frame(spline(curve_1, n = 20, method = "natural")), aes(x=x,y=y))+
  geom_line(data=data.frame(spline(curve_2, n = 20, method = "natural")), aes(x=x,y=y))+
  geom_line(data=data.frame(spline(curve_3, n = 100, method = "natural")), aes(x=x,y=y))+
  geom_line(data=data.frame(spline(curve_4st, n = 100, method = "natural")), aes(x=x,y=y))+
  geom_line(data=data.frame(spline(curve_4, n = 100, method = "natural")), aes(x=x,y=y))+
  geom_line(data=data.frame(spline(curve_5, n = 100, method = "natural")), aes(x=x,y=y))+
  theme_bw()+
  xlim(c(30,80)) +
  ylim(c(0,21)) +
  xlab('SiO2, wt%')+
  ylab('Na2O+K2O, wt%');
  
  if (labels)
  {
      p <- p + geom_text(aes(label=Name),hjust=0, vjust=0, size=font_size, lineheight=line_height);
  }
  
  if (annotate)
  {
      ants <- 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лейкограниты")
      );

      p <- p+ geom_text(data=ants,aes(x=x,y=y,label=label),size=font_size, lineheight=line_height);
  }
  
  return (p);
}