Create a square or hexagonal grid covering the bounding box of the geometry of an sf or sfc object

st_make_grid(
x,
cellsize = c(diff(st_bbox(x)[c(1, 3)]), diff(st_bbox(x)[c(2, 4)]))/n,
offset = st_bbox(x)[c("xmin", "ymin")],
n = c(10, 10),
crs = if (missing(x)) NA_crs_ else st_crs(x),
what = "polygons",
square = TRUE,
flat_topped = FALSE
)

## Arguments

x |
object of class sf or sfc |

cellsize |
target cellsize |

offset |
numeric of length 2; lower left corner coordinates (x, y) of the grid |

n |
integer of length 1 or 2, number of grid cells in x and y direction (columns, rows) |

crs |
object of class `crs` ; coordinate reference system of the target of the target grid in case argument `x` is missing, if `x` is not missing, its crs is inherited. |

what |
character; one of: `"polygons"` , `"corners"` , or `"centers"` |

square |
logical; if `FALSE` , create hexagonal grid |

flat_topped |
logical; if `TRUE` generate flat topped hexagons, else generate pointy topped |

## Value

Object of class `sfc`

(simple feature geometry list column) with, depending on `what`

and `square`

,
square or hexagonal polygons, corner points of these polygons, or center points of these polygons.

## Examples

#> although coordinates are longitude/latitude, st_intersects assumes that they are planar