geometricRoadPadPropSearchTable: Calculate the areas of intersection of a series of nested annuli with an idealized access road and turbine pad.

Description

Calculate area of annulus bisected by 2 parallel lines (e.g. a road of a road/pad plot).

Usage

geometricRoadPadPropSearchTable(
  padRadius,
  roadWidth,
  maxSearchRadius,
  mastRadius,
  annulusWidth = 1,
  ...
)

Value

Data frame of proportion of area searched for each annulus. distanceFromTurbine column represents the outer radius of each annulus.

Arguments

padRadius: Integer, radius of turbine pad from the center of the turbine.
roadWidth: Integer, width of road leading to turbine pad.
maxSearchRadius: Integer, maximum search distance from the center of turbine.
mastRadius: Integer, radius of the turbine mast.
annulusWidth: Integer, width of annulus, default is 1.
...: Currently ignored.

Details

Searches are conducted on the road and turbine pad around wind turbines for bird and bat fatalities. This function creates a data frame of proportion of area searched within each annulus ring on an idealized road and pad. The turbine is assumed to be centered on a perfectly circular turbine pad with radius padRadius, and a perfectly straight access road of width roadWidth is oriented from the center of the circle away from the turbine. (The resulting road and pad looks like rather like a lollipop.)

The mastRadius argument is to account for the area taken up by the turbine mast.

The arguments padRadius, roadWidth, mastRadius, and annulusWidth are all rounded to the nearest integer. The maxSearchDistance is rounded up (ceiling function) to an integer. If half units are needed, then convert to a smaller unit. See examples.

Examples

Run this code


pad <- 10 #meters, turbine pad radius
road <- 4 #meters, width of the road to the turbine pad
maxDistance <- 100 #meters, max distance
mast <- 2 #meters, turbine mast radius

## proportion are area searched at each annulus
propSearch <- geometricRoadPadPropSearchTable(padRadius = pad,
                                              roadWidth = road,
                                              maxSearchRadius = maxDistance,
                                              mastRadius = mast)
head(propSearch, 20)

## if half meter annulus rings are desired:
convert <- 100 # meters * 100 = centimeters

## units in centimeters
propSearchHalfMeter <- geometricRoadPadPropSearchTable(padRadius = pad * convert,
                                              roadWidth = road*convert,
                                              maxSearchRadius = maxDistance * convert,
                                              mastRadius = mast * convert,
                                              annulusWidth = 50) ##50cm = half a meter
head(propSearchHalfMeter, 30)

## convert back to meters
propSearchHalfMeter$distanceFromTurbine <- propSearchHalfMeter$distanceFromTurbine/convert
head(propSearchHalfMeter, 30)