Профили грунта с цветными объемными фракциями с «aqp» в R - PullRequest
Новая
       41

Профили грунта с цветными объемными фракциями с «aqp» в R

0 голосов
/ 11 октября 2019

Я пытаюсь построить профиль почвы в R, используя пакет aqp: алгоритмы количественной педологии. Профиль должен представлять цвет матрицы, плюс цвет пятен и процент. Для этой цели я использую функцию addVolumeFraction, которая в некоторой степени хорошо работает: она наносит точки на профиль, соответствующие правильному проценту пятнистости для каждого горизонта, но не назначает соответствующие цвета. Вот пример: 

#Variables for the soil profile
id <- rep(1, 4)
hor <- c("H1", "H2", "H3", "H4")
tops <- c(0,15,35,60)
bottoms <- c(15, 35, 60, 95)
mx_Hex <- c("#695F59FF", "#A59181FF", "#9E9388FF", "#A59181FF")
mot_Hex <- c("#EEB422","#EEB422", "#CD4F39", "#CD4F39")
mot_perc <- c(5, 10, 40, 8)

#Soil profile df
soildf <- data.frame(id,hor,tops,bottoms, mx_Hex, mot_Hex, mot_perc)
soildf$mx_Hex <- as.character(mx_Hex) #the class "SoilProfileCollection" needs colors as characters
soildf$mot_Hex <- as.character(mot_Hex)

# Transform df to "SoilProfileCollection"
depths(soildf) <- id ~ tops + bottoms

#Plot
plot(soildf, name = "hor", color = "mx_Hex", divide.hz = TRUE)
addVolumeFraction(soildf, "mot_perc",pch = 19, cex.min = 0.4, cex.max = 0.5, col = soildf$mot_Hex)

График профиля почвы

Как видно на графике, цвета пятен смешаны вдоль профиля. Я бы хотел, чтобы вместо этого были пятна определенного цвета для соответствующего горизонта. Кто-нибудь может помочь мне решить эту проблему?

Спасибо !!

...