create_coverage_map: Create a best server map and a coverage map.

View source: R/create_coverage_map.R

create_coverage_mapR Documentation

Create a best server map and a coverage map.

Description

The function create_coverage_map to create a coverage map. It shows the best signal strength per raster tile (either in dBm or in signal dominance). The function create_best_server_map is used to create a best (area) server map. It shows the best cell per raster tile.

Usage

create_coverage_map(strength, raster, type = c("dBm", "s"), cells = NULL)

create_best_server_map(llh, raster, cells = NULL)

Arguments

strength

a signal strength model object, which is the result of compute_sig_strength

raster

raster object that contains the raster tile index numbers (e.g. created with create_raster)

type

either "dBm" for absolute signal strength values and "s" for signal dominance values

cells

selection of cells

llh

a likelihood object, e.g. the result of create_voronoi_llh and create_strength_llh

See Also

vignette("mobloc")

Examples

## Not run: 
# set parameters
ZL_param <- mobloc_param()

# load data
data("ZL_cellplan", "ZL_muni", "ZL_elevation", "ZL_landuse")

# create environment layer (needed to calculate path loss exponent (ple))
ZL_envir <- combine_raster_layers(ZL_landuse, weights = c(1, 1, 1, 0, 0))

# validate cellplan
ZL_cellplan <- validate_cellplan(ZL_cellplan, param = ZL_param, region = ZL_muni,
                                 envir = ZL_envir, elevation = ZL_elevation)

# create raster
ZL_bbox <- sf::st_bbox(c(xmin = 4012000, ymin = 3077000, xmax = 4048000, ymax = 3117000),
                       crs = sf::st_crs(3035))
ZL_raster <- create_raster(ZL_bbox)

# compute the signal strength model
ZL_strength <- compute_sig_strength(cp = ZL_cellplan, raster = ZL_raster,
                                    elevation = ZL_elevation, param = ZL_param)

# create likelihoods
ZL_strength_llh <- create_strength_llh(ZL_strength, param = ZL_param)
ZL_voronoi_llh <- create_voronoi_llh(ZL_cellplan, ZL_raster)

bsm_strength <- create_best_server_map(ZL_strength_llh, raster = ZL_raster)
bsm_voronoi <- create_best_server_map(ZL_voronoi_llh, raster = ZL_raster)

bsms <- brick(bsm_strength, bsm_voronoi)
names(bsms) <- c("strength", "voronoi")

if (require(tmap)) {
    tmap_mode("view")
    tm_shape(bsms) +
        tm_raster(stretch.palette = FALSE) +
    qtm(move_cells_into_prop_direction(ZL_cellplan, 100))
}

## End(Not run)

MobilePhoneESSnetBigData/mobloc documentation built on Feb. 18, 2024, 3:41 a.m.