diff --git a/articles/Introduction-to-SoilProfileCollection-Objects.html b/articles/Introduction-to-SoilProfileCollection-Objects.html index 44d8bb01..f511528a 100644 --- a/articles/Introduction-to-SoilProfileCollection-Objects.html +++ b/articles/Introduction-to-SoilProfileCollection-Objects.html @@ -373,8 +373,8 @@
#> [1] 21 27 32 55 25 34 3 15 27 32 25 31 33 13 21 23 15 17 12 19 14 14 22 25 40 51 67 24 25 32
sp4$elevation # vector of simulated elevation (site data)
#> [1] 1009.9889 1006.1053 944.9123 932.4596 1065.9364 1022.1107 1368.4942 1113.8545 1122.1764
-#> [10] 1212.8237
+#> [1] 1079.8920 895.2652 1087.5804 829.4244 743.0426 1042.8912 755.5324 1150.2872 1117.3215
+#> [10] 789.5267
# unit-length value explicitly targeting site data
site(sp4)$collection_id <- 1
@@ -430,12 +430,12 @@ Horizon and Site Data#>
#> ----- Sites (6 / 10 rows | 5 / 5 columns) -----
#> id elevation collection_id constant group
-#> colusa 1009.9889 1 1 A
-#> glenn 1006.1053 1 1 B
-#> kings 944.9123 1 1 A
-#> mariposa 932.4596 1 1 B
-#> mendocino 1065.9364 1 1 A
-#> napa 1022.1107 1 1 B
+#> colusa 1079.8920 1 1 A
+#> glenn 895.2652 1 1 B
+#> kings 1087.5804 1 1 A
+#> mariposa 829.4244 1 1 B
+#> mendocino 743.0426 1 1 A
+#> napa 1042.8912 1 1 B
#> [... more sites ...]
#>
#> Spatial Data:
@@ -538,14 +538,13 @@ Object Metadata
# metadata structure
str(metadata(sp4))
#> List of 7
+#> List of 6
#> $ aqp_df_class : chr "data.frame"
#> $ aqp_group_by : chr ""
#> $ aqp_hzdesgn : chr "name"
#> $ aqp_hztexcl : chr ""
#> $ depth_units : chr "cm"
-#> $ stringsAsFactors: logi FALSE
-#> $ original.order : int [1:30] 1 2 3 4 5 6 7 8 9 10 ...
+#> $ stringsAsFactors: logi FALSE
# alter the depth unit metadata attribute
depth_units(sp4) <- 'inches' # units are really 'cm'
@@ -557,14 +556,13 @@ Object Metadata
# check new values have been added
str(metadata(sp4))
#> List of 10
+#> List of 9
#> $ aqp_df_class : chr "data.frame"
#> $ aqp_group_by : chr ""
#> $ aqp_hzdesgn : chr "name"
#> $ aqp_hztexcl : chr ""
#> $ depth_units : chr "inches"
#> $ stringsAsFactors: logi FALSE
-#> $ original.order : int [1:30] 1 2 3 4 5 6 7 8 9 10 ...
#> $ describer : chr "DGM"
#> $ date : Date[1:1], format: "2009-01-01"
#> $ citation : chr "McGahan, D.G., Southard, R.J, Claassen, V.P. 2009. Plant-Available Calcium Varies Widely in Soils on Serpentini"| __truncated__
@@ -594,16 +592,16 @@ Spatial Data# extract coordinates as matrix
getSpatial(sp4)
#> x y
-#> [1,] 353783.1 4109357
-#> [2,] 354163.2 4109556
-#> [3,] 354088.0 4109566
-#> [4,] 353944.1 4109435
-#> [5,] 353877.0 4109613
-#> [6,] 353957.6 4109556
-#> [7,] 353994.9 4109491
-#> [8,] 354048.2 4109584
-#> [9,] 353947.6 4109520
-#> [10,] 353937.6 4109442
+#> [1,] 353980.0 4109481
+#> [2,] 353915.3 4109323
+#> [3,] 354097.9 4109695
+#> [4,] 354060.2 4109444
+#> [5,] 354070.8 4109407
+#> [6,] 354016.3 4109710
+#> [7,] 354174.9 4109467
+#> [8,] 354082.5 4109446
+#> [9,] 354098.5 4109491
+#> [10,] 353971.3 4109562
# get/set spatial reference system using prj()<-
prj(sp4) <- '+proj=utm +zone=11 +datum=NAD83'
@@ -1660,13 +1658,13 @@ Aggregation over “slabs”# note: horizon names are lost due to aggregation
head(d.gsm, 7)
#> variable id value contributing_fraction top bottom
-#> 1 p1 1 19.77679 1.0 0 5
-#> 2 p1 1 19.79667 1.0 5 15
-#> 3 p1 1 19.84657 1.0 15 30
-#> 4 p1 1 21.41035 1.0 30 60
-#> 5 p1 1 23.46578 0.7 60 100
-#> 6 p1 1 NaN 0.0 100 200
-#> 7 p1 2 24.10446 1.0 0 5
+#> 1 p1 1 15.06565 1.0000000 0 5
+#> 2 p1 1 15.11775 1.0000000 5 15
+#> 3 p1 1 15.58663 1.0000000 15 30
+#> 4 p1 1 18.64598 1.0000000 30 60
+#> 5 p1 1 23.29671 1.0000000 60 100
+#> 6 p1 1 23.03729 0.2982456 100 200
+#> 7 p1 2 23.49896 1.0000000 0 5
A simple graphical comparison of the original and re-aligned soil
profile data, after converting slab()
result from long
-> wide format with {data.table} dcast()
:
simulateColor()
gains new method mvnorm
for simulating plausible colors
profileInformationIndex()
+profileInformationIndex()
, dice()
, slab()
, spc2mpspline()
, fillHzGaps()
, and flagOverlappingHz()
huePositionCircle()