Obtain the corners and bounding box of the MODFLOW grid
rmf_extent(dis, prj = rmf_get_prj(dis))
dis |
|
---|---|
prj | optional |
a list with (1) the coordinates of the grid corners as a data.frame and (2) the bounding box of the grid as
a sf bbox
object
dis <- rmf_create_dis() prj <- rmf_create_prj(origin = c(152082, 168000.2), rotation = -12, crs = 31370) rmf_extent(dis, prj)#> $corners #> x y #> ll 152082.0 168000.2 #> ul 152289.9 168978.3 #> ur 153268.1 168770.4 #> lr 153060.1 167792.3 #> #> $bbox #> xmin ymin xmax ymax #> 152082.0 167792.3 153268.1 168978.3 #>