对于我的专长,我想绘制一些土壤样本。一切正常,但是颜色。 我用Munsell颜色获得它们,将它们更改为RGB代码,但现在看来R只是从中制造随机颜色。 怎么了?我检查了RGB颜色,它们是正确的。它应该全部是褐色。但是R只是给我胡说八道的颜色作为输出。
代码:
install.packages ("aqp")
install.packages("RColorBrewer")
install.packages("latticeExtra")
install.packages("plyr")
install.packages("reshape")
require(aqp)
require(RColorBrewer)
require(latticeExtra)
require(plyr)
require(reshape)
setwd("C:/R_Working_Directory")
dat <- read.csv2("Boden_R.csv", header = TRUE, dec=".", stringsAsFactors = FALSE)
dat <- within (dat, label <- paste(name, kurz, sep=" - "))
print(dat)
str(dat)
print(dat)
class(dat)
summary(dat)
depths(dat) <- id ~ top + bottom
str(dat)
class(dat)
summary(dat)
plot(dat, name = "name", cex.name=1.1)
数据:
structure(list(depth_units = "cm"), class = "data.frame", row.names = c(NA,
-1L)), horizons = structure(list(id = c("GOWP1", "GOWP1",
"GOWP1", "GOWP1", "GOWP1", "GOWP2", "GOWP2", "GOWP2", "GOWP2",
"GOWP2", "GOWP2", "GOWP3", "GOWP3", "GOWP3", "GOWP3", "GOWP3",
"GOWP3", "GOWP3", "GOWP4", "GOWP4", "GOWP4", "GOWP4", "GOWP4",
"GOWP4", "GOWP5", "GOWP5", "GOWP5", "GOWP5", "GOWP5", "GOWP5",
"GOWP6", "GOWP6", "GOWP6", "GOWP6", "GOWP6", "GOWP6", "GOWP7",
"GOWP7", "GOWP7", "GOWP7", "GOWP7", "GOWP7", "GOWP7", "GOWP8",
"GOWP8", "GOWP8", "GOWP8", "GOWP8", "GOWP8", "GOWP8"), top = c(0L,
25L, 50L, 73L, 95L, 0L, 11L, 36L, 55L, 85L, 95L, 0L, 11L,
28L, 45L, 66L, 73L, 94L, 0L, 12L, 25L, 50L, 84L, 93L, 0L,
8L, 46L, 76L, 88L, 97L, 0L, 16L, 40L, 47L, 73L, 98L, 0L,
9L, 40L, 47L, 67L, 81L, 92L, 0L, 13L, 36L, 43L, 66L, 95L,
99L), bottom = c(25L, 50L, 73L, 95L, 100L, 11L, 36L, 55L,
85L, 95L, 100L, 11L, 28L, 45L, 66L, 73L, 94L, 100L, 12L,
25L, 50L, 84L, 93L, 100L, 8L, 46L, 76L, 88L, 97L, 100L, 16L,
40L, 47L, 73L, 98L, 100L, 9L, 40L, 47L, 67L, 81L, 92L, 100L,
13L, 36L, 43L, 66L, 95L, 99L, 100L), name = c("KV", "Ah",
"Bv ", "Bbv", "C", "KV", "Ah", "Al-Bv", "Btv", "Bbv", "KV",
"KV", "Ah", "Al-Bv", "Btv", "Bv ", "Bbv", "KV", "KV", "MAh",
"MAh-MBv", "Btv", "Btv", "KV", "KV", "MAH/Bv", "(M)Bv", "M?",
"Bbv", "KV", "KV", "Ah", "Al-Bv", "B(t)v", "Bbv", "KV", "KV",
"Ah", "Ah/Bv", "Bv ", "Bbv", "Cv", "KV", "KV", "Ah", "Ah/Bv",
"Bv ", "Bbv", "Cv", "KV"), p1 = c(0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L),
soil_color = c("FFFFFF", "482B12", "B06A32", "AB6D28",
"AB6D28", "FFFFFF", "A59181", "9A7359", "6E3B19", "9A7359",
"FFFFFF", "FFFFFF", "482B12", "B06A32", "6E3B19", "B06A32",
"9A7359", "FFFFFF", "FFFFFF", "482B12", "2B190D", "6E3B19",
"9A7359", "FFFFFF", "FFFFFF", "532404", "B06A32", "91765F",
"9A7359", "FFFFFF", "FFFFFF", "482B12", "A59181", "B06A32",
"A66E46", "FFFFFF", "FFFFFF", "482B12", "A59181", "B06A32",
"9A7359", "AB6D28", "FFFFFF", "FFFFFF", "482B12", "A59181",
"B06A32", "B06A32", "AB6D28", "FFFFFF"), lang = c("Kernverlust",
"dunkelbraun", "rostfärblich-braun", "hellbraun", "hellbraun",
"Kernverlust", "hellgrau-braun", "rostfärblich-hellbraun",
"dunkelrostfärblich-braun", "rostfärblich-hellbraun",
"Kernverlust", "Kernverlust", "dunkelbraun", "rostfärblich-braun",
"dunkelrostfärblich-braun", "rostfärblich-braun", "rostfärblich-hellbraun",
"Kernverlust", "Kernverlust", "dunkelbraun", "dunkelbraun-braun",
"dunkelrostfärblich-braun", "rostfärblich-hellbraun",
"Kernverlust", "Kernverlust", "dunkelbraun-rostfärblich-braun",
"rostfärblich-braun", "graubraun", "rostfärblich-hellbraun",
"Kernverlust", "Kernverlust", "dunkelbraun", "hellgrau-braun",
"rostfärblich-braun", "hellrostfärblich-braun", "Kernverlust",
"Kernverlust", "dunkelbraun", "hellgrau-braun", "rostfärblich-braun",
"rostfärblich-hellbraun", "hellbraun", "Kernverlust",
"Kernverlust", "dunklebraun", "hellgrau-braun", "rostfärblich-braun",
"rostfärblich-braun", "hellbraun", "Kernverlust"), kurz = c("",
"dbn", "rflibn", "hbn", "hbn", "", "hgrbn", "rflihbn",
"drflibn", "rflihbn", "", "", "dbn", "rflibn", "drflibn",
"rflibn", "rflihbn", "", "", "dbn", "dbnbn", "drflibn",
"rflihbn", "", "", "dbnrflibn", "rflibn", "grbn", "rflihbn",
"", "", "dbn", "hgrbn", "rflibn", "hrflibn", "", "",
"dbn", "hgrbn", "rflibn", "rflihbn", "hbn", "", "", "dbn",
"hgrbn", "rflibn", "rflibn", "hbn", ""), label = c("KV - ",
"Ah - dbn", "Bv - rflibn", "Bbv - hbn", "C - hbn", "KV - ",
"Ah - hgrbn", "Al-Bv - rflihbn", "Btv - drflibn", "Bbv - rflihbn",
"KV - ", "KV - ", "Ah - dbn", "Al-Bv - rflibn", "Btv - drflibn",
"Bv - rflibn", "Bbv - rflihbn", "KV - ", "KV - ", "MAh - dbn",
"MAh-MBv - dbnbn", "Btv - drflibn", "Btv - rflihbn",
"KV - ", "KV - ", "MAH/Bv - dbnrflibn", "(M)Bv - rflibn",
"M? - grbn", "Bbv - rflihbn", "KV - ", "KV - ", "Ah - dbn",
"Al-Bv - hgrbn", "B(t)v - rflibn", "Bbv - hrflibn", "KV - ",
"KV - ", "Ah - dbn", "Ah/Bv - hgrbn", "Bv - rflibn",
"Bbv - rflihbn", "Cv - hbn", "KV - ", "KV - ", "Ah - dbn",
"Ah/Bv - hgrbn", "Bv - rflibn", "Bbv - rflibn", "Cv - hbn",
"KV - "), hzID = 1:50), row.names = c(NA, 50L), class = "data.frame"),
site = structure(list(id = c("GOWP1", "GOWP2", "GOWP3", "GOWP4",
"GOWP5", "GOWP6", "GOWP7", "GOWP8")), class = "data.frame", row.names = c(NA,
-8L)), sp = new("SpatialPoints", coords = structure(0, .Dim = c(1L,
1L)), bbox = structure(NA, .Dim = c(1L, 1L)), proj4string = new("CRS",
projargs = NA_character_)), diagnostic = structure(list(), .Names = character(0), row.names = integer(0), class = "data.frame"),
restrictions = structure(list(), .Names = character(0), row.names = integer(0), class = "data.frame"))