Split a spatial object into quadrants

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

Other geo: bbox_scale(), gclip(), geo_bb_matrix(), geo_bb(), mapshape_available(), mapshape(), reproject()

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))