masqOLS: Extract raster pixel values within the boundaries of a...

In chrisvwn/Rnightlights: Satellite Nightlight Data Extraction

Description

Extract raster pixel values within the boundaries of a polygon

Usage

masqOLS(ctryPoly, ctryRast, idx, retVal = NULL, configName, extension)

Arguments

ctryPoly	the country Polygon layer as SpatialPolygon
ctryRast	the clipped country raster
idx	the index of the polygon in the country polygon layer (shp)
retVal	Whether to return the raster data as a vector, or data.frame with spatial context NULL returns a vector of all values, colrowval returns a data.frame with row, col and raster value while lonlatval returns a data.frame with lon,lat and val.
configName	character the config shortname of raster being processed
extension	character the extension of raster being processed

Value

numeric vector of radiances

Examples

## Not run: 
ctryPoly <- rgdal::readOGR(getPolyFnamePath("KEN"), getCtryShpLyrNames("KEN", 1))
ctryRaster <- raster::raster(getCtryRasterOutputFnamePath("KEN", "OLS", "1999"))
temp <- NULL
KenAdm1Sum <- NULL
for (i in 1:length(ctryPoly@polygons))
{
 temp$name <- as.character(ctryPoly@data$NAME_1[i])
 temp$sum <- sum(masqOLS(ctryPoly, ctryRaster, i), na.rm=T)

 KenAdm1Sum <- rbind(KenAdm1Sum)
}

## End(Not run)

chrisvwn/Rnightlights documentation built on Sept. 7, 2021, 1:44 a.m.