tmaptools (version 1.2-4)

points_to_raster: Bin spatial points to a raster

Description

Bin spatial points to a raster. For each raster cell, the number of points are counted. Optionally, a factor variable can be specified by which the points are counts are split.

Usage

points_to_raster(shp, nrow = NA, ncol = NA, N = 250000, by = NULL,
  to.Raster = FALSE)

Arguments

shp

shape object. a SpatialPoints(DataFrame), a SpatialGrid(DataFrame), or an sf object that can be coerced as such.

nrow

number of raster rows. If NA, it is automatically determined by N and the aspect ratio of shp.

ncol

number of raster columns. If NA, it is automatically determined by N and the aspect ratio of shp.

N

preferred number of raster cells.

by

name of a data variable which should be a factor. The points are split and counted according to the levels of this factor.

to.Raster

logical; should the output be a Raster object (TRUE), or a SpatialGridDataFrame (FALSE). If TRUE, a RasterBrick is returned when by is specified, and a RasterLayer when by is unspecified.

Value

A SpatialGridDataFrame, or a Raster object when (to.Raster=TRUE)

Details

This function is a wrapper around rasterize.

References

Tennekes, M., 2018, tmap: Thematic Maps in R, Journal of Statistical Software, 84(6), 1-39, DOI

See Also

poly_to_raster

Examples

# NOT RUN {
if (require(tmap)) {
    data(NLD_muni, NLD_prov)

    # sample points (each point represents 1000 people)
    NLD_muni_points <- sample_dots(NLD_muni, vars = "population",
        w=1000, convert2density = TRUE)

    # dot map
    tm_shape(NLD_muni_points) + tm_dots()

    # convert points to raster
    NLD_rst <- points_to_raster(NLD_muni_points, N = 1e4)

    # plot raster
    tm_shape(NLD_rst) +
    	tm_raster() +
    tm_shape(NLD_prov) +
    	tm_borders() +
    tm_format_NLD() + tm_style_grey()
}
# }