interpol_view: Plot an interpolated viewshed

In windfarmGA: Genetic Algorithm for Wind Farm Layout Optimization

Description

Plot an interpolated view of the viewshed analysis

Usage

interpol_view(
  res,
  plot = TRUE,
  breakseq,
  breakform = NULL,
  plotDEM = FALSE,
  fun = mean,
  pal = NULL,
  ...
)

Arguments

res	The resulting list from viewshed
plot	Should the result be plotted? Default is TRUE
breakseq	The breaks for value plotting. By default, 5 equal intervals are generated.
breakform	If breakseq is missing, a sampling function to calculate the breaks, like quantile, fivenum, etc.
plotDEM	Plot the DEM? Default is FALSE
fun	Function used for rasterize. Default is mean
pal	A color palette
...	Arguments passed on to plot.

Value

An interpolated raster

See Also

Other Viewshed Analysis: cansee(), plot_viewshed(), rasterprofile(), viewTo(), viewshed()

Other Plotting Functions: plot_cloud(), plot_development(), plot_evolution(), plot_fitness_evolution(), plot_heatmap(), plot_parkfitness(), plot_result(), plot_viewshed(), plot_windfarmGA(), plot_windrose(), random_search_single()

Examples

## Not run: 
library(raster)
library(sf)
Polygon1 <- sf::st_as_sf(sf::st_sfc(
  sf::st_polygon(list(cbind(
    c(4496482, 4496482, 4499991, 4499991, 4496482),
    c(2666272, 2669343, 2669343, 2666272, 2666272)))),
  crs = 3035
))
DEM_meter <- getDEM(Polygon1)

turbloc = st_sample(DEM_meter[[2]], 10, type = "random");
res <- viewshed(r = DEM_meter[[1]], shape=DEM_meter[[2]], 
                turbine_locs = turbloc,  h1=1.8, h2=50)
interpol_view(res, plotDEM = TRUE)

interpol_view(res, breakseq = seq(0,max(colSums(res$Result)),1))
interpol_view(res, plotDEM = FALSE, breakform = quantile)
interpol_view(res, breakform = factor)

## Different color palettes
interpol_view(res, plotDEM = TRUE, pal=topo.colors)
interpol_view(res, plotDEM = TRUE, pal=colorRampPalette(c("white","purple")))

## ... Arguments are past on to the raster plot method
interpol_view(res, plotDEM = TRUE, alpha=0.5)
interpol_view(res, plotDEM = FALSE, breakseq = seq(0,10,1), colNA="black")


## End(Not run)

windfarmGA documentation built on May 5, 2021, 5:08 p.m.