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

In Rnightlights: Satellite Nightlight Data Extraction

Description

Extract raster pixel values within the boundaries of a polygon

Usage

masqOLS(shp, rast, i, retVal = "colrowval", configName)

Arguments

shp	the country Polygon layer as SpatialPolygon
rast	the clipped country raster
i	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 type 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)

Rnightlights documentation built on Aug. 29, 2019, 5:02 p.m.