Reads in soil clay, silt, sand fractions and returns soil classification based on USDA or Norwegian definitions.
classify_soil(
input,
output,
wd = NULL,
version = "NOR",
plot_result = TRUE,
interactive = TRUE,
sheet = 1,
decimal = ",",
append = TRUE,
new_sums = FALSE,
write = FALSE
)path to input file (supported: .xlsx/.csv)
file name of output (select between csv or excel) (ie: .xlsx/.csv)
path to desired working directory. Default is current.
"NOR" for Norwegian, "USDA" for USDA.
TRUE/FALSE display diagnostic plot?
if plot is true, do you want an interactive or static plot?
If using excel, specify which sheet the data is in. 1 by default
Specify decimal separator ("," or "."). By default ","
append results to existing dataframe? (TRUE) or return single column (FALSE). default TRUE.
append new "rounded sums" to input dataframe? Default FALSE
optional flag to stop the function from writing any files. Default FALSE
Returns path to written file
Automatically Compatible with csv or xlsx.
Tries to automatically guess format of data, and acts accordingly.
Automatically rounds total fraction to 100%, if differences are under 1%.
Exports classification based on user settings
Creates an interactive diagnostic plot, as well as saving one to output directory.
Note: the script will automatically try to identify which columns contain the clay, silt, sand fractions and also try to identify which format they are in. If this does not work, you might get incorrect or flawed results, so make sure to READ what THE CONSOLE prints out, to make sure it did the correct things.
If the script cannot read your file automatically, then the following format is recommended:
Clay column: "clay" in percent format (%), decimal = "." Silt column: "silt" in percent format (%), decimal = "." Sand column: "sand" in percent format (%), decimal = "."
# for demonstration purposes, use path of package
example_file_path <- system.file(package = "miljotools", "/extdata/soil_classify/example.csv")
classify_soil(
input = example_file_path,
output = "output_file.csv",
version = "NOR",
plot_result = TRUE,
interactive = FALSE,
decimal = ".",
write = FALSE # FALSE just for demonstration
)
#> Identifying column: sand
#> sand column identified as ' sand ' index= 3
#> > mean value of sand column is 25
#> Assuming values are in PERCENT form. ie. 43.1% (instead of 0.431)
#> Identifying column: silt
#> silt column identified as ' silt ' index= 4
#> > mean value of silt column is 23.1
#> Assuming values are in PERCENT form. ie. 43.1% (instead of 0.431)
#> Identifying column: clay
#> clay column identified as ' clay ' index= 2
#> > mean value of clay column is 51.9
#> Assuming values are in PERCENT form. ie. 43.1% (instead of 0.431)
#>
#>
#> Checking for fraction sums over 100 %:
#> > 0 entries have a fractional sum of over 100.
#> > 0 entries have a fractional sum of under 100.
#> adjusting the sand fraction and rounding...
#>
#> > classifiying NORWEGIAN textures
#> Error in guides_train(scales = scales$non_position_scales(), theme = theme, guides = guides, labels = labels): could not find function "guides_train"