create.pgrid
creates a grid of locations from the boundaries of domain and other information.
create.pgrid(
xmin,
xmax,
ymin,
ymax,
nx,
ny,
midpoints = FALSE,
poly.coords = NULL
)
Returns an object of class pgrid
with the following components:
An \(n \times 2\) matrix of locations (the midpoints of the pixelized grid).
The number of rows in pgrid.
A vector of 0s and 1s indicating whether the midpoint of each pixel is in the convex hull of poly.coords
. If poly.coords
is not provided, this is a vector of 1s.
The pixel boundaries in the x direction.
The pixel boundaries in the y direction.
The pixel midpoints in the x direction.
The pixel midpoints in the y direction.
The minimum value of the boundary of the x coordinates of the spatial domain.
The maximum value of the boundary of the x coordinates of the spatial domain.
The minimum value of the boundary of the y coordinates of the spatial domain.
The maximum value of the boundary of the y coordinates of the spatial domain.
The number of gridpoints/cells/pixels in the x direction.
The number of gridpoints/cells/pixels in the y direction.
A logical value (TRUE
or FALSE
) indicating whether the boundary values are for the midpoint of a pixel (midpoints = TRUE
) or for the boundary of the spatial domain in general (midpoints = FALSE
), in which case the midpoints are calculated internally). Default is FALSE
.
An \(n \times 2\) matrix with the coordinates specifying the polygon vertices of the true spatial domain of interest within the rectangular boundaries provided by xmin
, xmax
, ymin
, and ymax
. If this is provided, the pgrid
returned will be within the convex hull of poly.coords
.
Joshua French
The key argument in the function midpoints. If this is TRUE
, it is assumed that the boundaries of the spatial domain correspond to the midpoints of the cell/pixel in the grid. Otherwise, it is assumed that the boundaries correspond to the actual borders of the region of interest. If poly.coords
is supplied, the grid returned is the grid of midpoints contained in the convex hull of poly.coords
.
pgrida <- create.pgrid(0, 1, 0, 1, nx = 50, ny = 50, midpoints = FALSE)
pgridb <- create.pgrid(.01, .99, .01, .99, nx = 50, ny = 50, midpoints = TRUE)
Run the code above in your browser using DataLab