powered by
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.
horizonSearch(x, azimuth, maxDist = 1000, degrees = FALSE, ll = TRUE, filename = "", blockSize = NULL)
RasterLayer object.
numeric. Azimuth angle.
integer. Maximum search distance in meters.
logical. If TRUE, input (azimuth) and output (horizon elevation) angles are in degrees; else in radian.
logical. If TRUE, coordinates should be in degrees; else planar.
character. Output raster filename. Optional.
integer>0. Number of raster rows processed simultaneously. Default is the number of rows in the input raster divided by 10.
RasterLayer object
Azimuth is measured from North over East.
The 'blockSize' parameter may be reduced to avoid memory errors when processing large input RasterLayer objects.
svf
# NOT RUN { r <- getData('alt', country='ALB') h <- horizonSearch(r, 60, degrees=TRUE, maxDist=2000, ll=TRUE) # }
Run the code above in your browser using DataLab