protect_neighborhood: protects raster by summing over the neighborhood

In sdcSpatial: Statistical Disclosure Control for Spatial Data

protects raster by summing over the neighborhood

Description

protects raster by summing over the neighborhood

Usage

protect_neighborhood(x, radius = 10 * raster::res(x$value)[1], ...)

Arguments

x	sdc_raster() object to be protected
radius	of the neighborhood to take
...	not used at the moment

Value

sdc_raster object

Examples

data(enterprises)

# create a sdc_raster from point data with raster with
# a resolution of 200m
production <- sdc_raster(enterprises, variable = "production"
                        , r = 200, min_count = 3)

print(production)

# plot the raster
zlim <- c(0, 3e4)
# show which raster cells are sensitive
plot(production, zlim=zlim)

# let's smooth to reduce the sensitivity
smoothed <- protect_smooth(production, bw = 400)
plot(smoothed)

neighborhood <- protect_neighborhood(production, radius=1000)
plot(neighborhood)

# what is the sensitivy fraction?
sensitivity_score(neighborhood)

sdcSpatial documentation built on March 24, 2022, 5:05 p.m.