Skip to contents

Number of rows, columns, depths, cells, and layers

Source: R/dim.r
dim.Rd

For GRegions: Number of rows, columns, depths, and cells:

  • dim(): Rows and columns

  • dim3d(): Rows, columns, and depths

  • nrow(): Rows

  • ncol(): Columns

  • ndepth(): Depths (for 3-dimensional rasters only)

  • ncell(): Number of cells (2 dimensions)

  • ncell3d(): Number of cells (3 dimensions)

For GRasters: As above, plus number of layers:

  • nlyr(): Layers (number of "stacked" rasters–different from depths of a raster).

For GVectors: Number of geometries and fields (columns):

  • dim(): Number of geometries and number of columns in data table

  • nrow(): Number of geometries

  • ncol(): Number of columns in data table

Usage

# S4 method for class 'GRegion'
dim(x)

# S4 method for class 'missing'
dim3d(x)

# S4 method for class 'GRegion'
dim3d(x)

# S4 method for class 'missing'
nrow(x)

# S4 method for class 'GRegion'
nrow(x)

# S4 method for class 'missing'
ncol(x)

# S4 method for class 'GRegion'
ncol(x)

# S4 method for class 'missing'
ndepth(x)

# S4 method for class 'GRegion'
ndepth(x)

# S4 method for class 'missing'
ncell(x)

# S4 method for class 'GRegion'
ncell(x)

# S4 method for class 'missing'
ncell3d(x)

# S4 method for class 'GRegion'
ncell3d(x)

# S4 method for class 'GVector'
dim(x)

# S4 method for class 'GVector'
nrow(x)

# S4 method for class 'GVector'
ncol(x)

# S4 method for class 'missing'
nlyr(x)

# S4 method for class 'GRaster'
nlyr(x)

Arguments

x

A GRegion, GRaster, GVector, or missing. If missing, then the dimensions of the currently active "region" are returned (see vignette("regions", package = "fasterRaster")).

Value

A numeric value or vector.

See also

ngeom(), nsubgeom(), nacell(), nonnacell(), terra::dim()

Examples

if (grassStarted()) {

# Setup
library(sf)
library(terra)

# Example data
madElev <- fastData("madElev")
madForest2000 <- fastData("madForest2000")
madCoast0 <- fastData("madCoast0")
madRivers <- fastData("madRivers")
madDypsis <- fastData("madDypsis")

### GRaster properties

# convert SpatRasters to GRasters
elev <- fast(madElev)
forest <- fast(madForest2000)

# plot
plot(elev)

dim(elev) # rows, columns, depths, layers
nrow(elev) # rows
ncol(elev) # columns
ndepth(elev) # depths
nlyr(elev) # layers

res(elev) # resolution

ncell(elev) # cells
ncell3d(elev) # cells (3D rasters only)

topology(elev) # number of dimensions
is.2d(elev) # is it 2D?
is.3d(elev) # is it 3D?

minmax(elev) # min/max values

# name of object in GRASS
sources(elev)

# "names" of the object
names(elev)

# coordinate reference system
crs(elev)

# extent (bounding box)
ext(elev)

# data type
datatype(elev)

# assigning
copy <- elev
copy[] <- pi # assign all cells to the value of pi
copy

# concatenating multiple GRasters
rasts <- c(elev, forest)
rasts

# adding a raster "in place"
add(rasts) <- ln(elev)
rasts

# subsetting
rasts[[1]]
rasts[["madForest2000"]]

# assigning
rasts[[4]] <- elev > 500

# number of layers
nlyr(rasts)

# names
names(rasts)
names(rasts) <- c("elev_meters", "forest", "ln_elev", "high_elevation")
rasts

### GVector properties

# convert sf vectors to GVectors
coast <- fast(madCoast4)
rivers <- fast(madRivers)
dypsis <- fast(madDypsis)

# extent
ext(rivers)

W(rivers) # western extent
E(rivers) # eastern extent
S(rivers) # southern extent
N(rivers) # northern extent
top(rivers) # top extent (NA for 2D rasters like this one)
bottom(rivers) # bottom extent (NA for 2D rasters like this one)

# coordinate reference system
crs(rivers)
st_crs(rivers)

# column names and data types
names(coast)
datatype(coast)

# name of object in GRASS
sources(rivers)

# points, lines, or polygons?
geomtype(dypsis)
geomtype(rivers)
geomtype(coast)

is.points(dypsis)
is.points(coast)

is.lines(rivers)
is.lines(dypsis)

is.polygons(coast)
is.polygons(dypsis)

# dimensions
nrow(rivers) # how many spatial features
ncol(rivers) # hay many columns in the data frame

# number of geometries and sub-geometries
ngeom(coast)
nsubgeom(coast)

# 2- or 3D
topology(rivers) # dimensionality
is.2d(elev) # is it 2D?
is.3d(elev) # is it 3D?

# Update values from GRASS
# (Reads values from GRASS... will not appear to do anything in this case)
coast <- update(coast)

### operations on GVectors

# convert to data frame
as.data.frame(rivers)
as.data.table(rivers)

# subsetting
rivers[c(1:2, 5)] # select 3 rows/geometries
rivers[-5:-11] # remove rows/geometries 5 through 11
rivers[ , 1] # column 1
rivers[ , "NAM"] # select column
rivers[["NAM"]] # select column
rivers[1, 2:3] # row/geometry 1 and column 2 and 3
rivers[c(TRUE, FALSE)] # select every other geometry (T/F vector is recycled)
rivers[ , c(TRUE, FALSE)] # select every other column (T/F vector is recycled)

# removing data table
noTable <- dropTable(rivers)
noTable
nrow(rivers)
nrow(noTable)

# Refresh values from GRASS
# (Reads values from GRASS... will not appear to do anything in this case
# since the rivers object is up-to-date):
rivers <- update(rivers)

# Concatenating multiple vectors
rivers2 <- rbind(rivers, rivers)
dim(rivers)
dim(rivers2)

}