horizonSearch: Horizon elevation angle from DEM RasterLayer

Description Usage Arguments Details Value See Also Examples

Description

Computes the horizon elevation angle, within a maximum search distance, for a given azimith angle for each pixel in a digital elevation model raster object.

Usage

1
2
horizonSearch(x, azimuth, maxDist = 1000, degrees = FALSE, 
              ll = TRUE, filename = "", blockSize = NULL)

Arguments

x

RasterLayer object.

azimuth

numeric. Azimuth angle.

maxDist

integer. Maximum search distance in meters.

degrees

logical. If TRUE, input (azimuth) and output (horizon elevation) angles are in degrees; else in radian.

ll

logical. If TRUE, coordinates should be in degrees; else planar.

filename

character. Output raster filename. Optional.

blockSize

integer>0. Number of raster rows processed simultaneously. Default is the number of rows in the input raster divided by 10.

Details

Azimuth is measured from North over East.

The 'blockSize' parameter may be reduced to avoid memory errors when processing large input RasterLayer objects.

Value

RasterLayer object

See Also

svf

Examples

1
2
3
r <- getData('alt', country='ALB')

h <- horizonSearch(r, 60, degrees=TRUE, maxDist=2000, ll=TRUE)


Search within the horizon package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.