GSIF (version 0.5-5.1)

as.data.frame: Converts an object of class "SoilProfileCollection" to a data frame

Description

Converts an object of class "SoilProfileCollection" to an object of class "data.frame" with both site and horizon data sorted in one row. Each original column name in the horizons table receives a sufix *_A, B, …, Z where alphabetic letters represent horizon sequence.

Usage

# S4 method for SoilProfileCollection
as.data.frame(x, row.names = NULL, optional = FALSE, …)

Arguments

x

object of class "SoilProfileCollection"

row.names

character; giving the row names for the data frame (missing values are not allowed)

optional

logical; if `TRUE', setting row names and converting column names (to syntactic names: see `make.names') is optional

…

optional arguments

Details

The advantage of converting the "SoilProfileCollection" data to a single table is that, once both tables have been merged to a single data frame, it can be more easily exported and visualized in a GIS and/or imported into a data base.

See Also

as.geosamples, mpspline

Examples

# NOT RUN {
library(aqp)
library(plyr)
library(rgdal)
library(sp)
## sample profile from Nigeria:
lon = 3.90; lat = 7.50; id = "ISRIC:NG0017"; FAO1988 = "LXp" 
top = c(0, 18, 36, 65, 87, 127) 
bottom = c(18, 36, 65, 87, 127, 181)
ORCDRC = c(18.4, 4.4, 3.6, 3.6, 3.2, 1.2)
## prepare a SoilProfileCollection:
prof1 <- join(data.frame(id, top, bottom, ORCDRC), 
    data.frame(id, lon, lat, FAO1988), 
    type='inner')
depths(prof1) <- id ~ top + bottom
site(prof1) <- ~ lon + lat + FAO1988 
coordinates(prof1) <- ~ lon + lat
proj4string(prof1) <- CRS("+proj=longlat +datum=WGS84")
## convert to a simple table:
x <- as.data.frame(prof1)
str(x)
## horizons only
horizons <- getHorizons(x, idcol="id", sel=c("top", "bottom", "ORCDRC"))
horizons
# }