Split a spatial object (initially tested on SpatialPolygons) into quadrants.

quadrant(sp_obj, number_out = FALSE)

Arguments

sp_obj

Spatial object

number_out

Should the output be numbers from 1:4 (FALSE by default)

Details

Returns a character vector of NE, SE, SW, NW corresponding to north-east, south-east quadrants respectively. If number_out is TRUE, returns numbers from 1:4, respectively.

See also

Examples

data(zones)
sp_obj <- zones
(quads <- quadrant(sp_obj))
#> [1] "SE" "SW" "SE" "NE" "NE" "SW" "NW" "SE"
plot(sp_obj, col = factor(quads))
points(rgeos::gCentroid(sp_obj), col = "white")

# edge cases (e.g. when using rasters) lead to NAs
sp_obj <- raster::rasterToPolygons(raster::raster(ncol = 3, nrow = 3))
(quads <- quadrant(sp_obj))
#> [1] "SW" NA   "NE" NA   NA   NA   "NW" NA   "SE"
plot(sp_obj, col = factor(quads))