Get a data.frame with raster cell values, or coerce SpatialPolygons, Lines, or Points to a data.frame

as.matrix returns all values of a Raster* object as a matrix. For RasterLayers, rows and columns in the matrix represent rows and columns in the RasterLayer object. For other Raster* objects, the matrix returned by as.matrix has columns for each layer and rows for each cell.

as.array returns an array of matrices that are like those returned by as.matrix for a RasterLayer

If there is insufficient memory to load all values, you can use getValues or getValuesBlock to read chunks of the file. You could also first use sampleRegular

The methods for Spatial* objects allow for easy creation of a data.frame with the coordinates and attributes; the default method only returns the attributes data.frame

# S4 method for class 'Raster'
as.data.frame(x, row.names=NULL, optional=FALSE, xy=FALSE, 
              na.rm=FALSE, long=FALSE, ...)

# S4 method for class 'SpatialPolygons'
as.data.frame(x, row.names=NULL, optional=FALSE,
              xy=FALSE, centroids=TRUE, sepNA=FALSE, ...)

# S4 method for class 'SpatialLines'
as.data.frame(x, row.names=NULL, optional=FALSE, 
              xy=FALSE, sepNA=FALSE, ...)

Arguments

x

Raster* object

row.names

NULL or a character vector 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

xy

logical. If TRUE, also return the spatial coordinates

na.rm

logical. If TRUE, remove rows with NA values. This can be particularly useful for very large datasets with many NA values

long

logical. If TRUE, values are reshaped from a wide to a long format

centroids

logical. If TRUE return the centroids instead of all spatial coordinates (only relevant if xy=TRUE)

sepNA

logical. If TRUE the parts of the spatial objects are separated by lines that are NA (only if xy=TRUE and, for polygons, if centroids=FALSE

...

Additional arguments (none)

Value

data.frame

Examples

r <- raster(ncol=3, nrow=3)
values(r) <- sqrt(1:ncell(r))
r[3:5] <- NA
as.data.frame(r)
#>      layer
#> 1 1.000000
#> 2 1.414214
#> 3       NA
#> 4       NA
#> 5       NA
#> 6 2.449490
#> 7 2.645751
#> 8 2.828427
#> 9 3.000000
s <- stack(r, r*2)
as.data.frame(s)
#>    layer.1  layer.2
#> 1 1.000000 2.000000
#> 2 1.414214 2.828427
#> 3       NA       NA
#> 4       NA       NA
#> 5       NA       NA
#> 6 2.449490 4.898979
#> 7 2.645751 5.291503
#> 8 2.828427 5.656854
#> 9 3.000000 6.000000
as.data.frame(s, na.rm=TRUE)
#>    layer.1  layer.2
#> 1 1.000000 2.000000
#> 2 1.414214 2.828427
#> 6 2.449490 4.898979
#> 7 2.645751 5.291503
#> 8 2.828427 5.656854
#> 9 3.000000 6.000000