populR (version 0.2.1)

pp_estimate: Areal Interpolation of Population Data

Description

Areal Interpolation of Population Data

Usage

pp_estimate(
  target,
  source,
  sid,
  spop,
  volume = NULL,
  ancillary = NULL,
  point = FALSE,
  method
)

Value

An object of class sf including estimated population counts for target features using either awi or vwi

methods. The estimated population counts are stored in a new column called pp_est.

Arguments

target

An object of class sf that is used to interpolate data to. Usually, target may include polygon features representing building units

source

An object of class sf including data to be interpolated. Source may be a set of coarse polygon features such as city blocks or census tracts

sid

Source identification number

spop

Source population values to be interpolated

volume

Target feature volume information (height or number of floors). Required when method=vwi

ancillary

ancillary information

point

Whether to return point geometries (FALSE by default)

method

Two methods provided: awi (areal weighting interpolation) and vwi (volume weighting interpolation). awi proportionately interpolates the population values based on areal weights calculated by the area of intersection between the source and target zones. vwi proportionately interpolates the population values based on areal weights calculated by the area of intersection between the source and target zones multipled by the volume information (height or number of floors).

Examples

# read lib data
data('src')
data('trg')

# areal weighted interpolation - awi
pp_estimate(trg, src, sid = sid, spop = pop,
    method = awi)

# areal weighted interpolation - awi using point geometries
pp_estimate(trg, src, sid = sid, spop = pop,
    method = awi, point = TRUE)

# volume weighted interpolation - vwi
pp_estimate(trg, src, sid = sid, spop = pop,
    method = vwi, volume = floors)

# volume weighted interpolation - vwi using point geometries
pp_estimate(trg, src, sid = sid, spop = pop,
    method = vwi, volume = floors, point = TRUE)