inla.over_sp_mesh {INLA} | R Documentation |
Wrapper for the sp::over()
method to find triangle centroids
or vertices inside sp
polygon objects
inla.over_sp_mesh(x, y, type = c("centroid", "vertex"), ignore.CRS = FALSE)
x |
geometry (typically a |
y |
an |
type |
the query type; either |
ignore.CRS |
logical; whether to ignore the coordinate system information in |
A vector of triangle indices (when type
is 'centroid'
) or
vertex indices (when type
is 'vertex'
)
Haakon Bakka, bakka@r-inla.org, and Finn Lindgren finn.lindgren@gmail.com
# Create a polygon and a mesh obj <- sp::SpatialPolygons(list(Polygons(list(Polygon(rbind( c(0, 0), c(50, 0), c(50, 50), c(0, 50) ))), ID = 1 )), proj4string = inla.CRS("longlat_globe") ) mesh <- inla.mesh.create(globe = 2, crs = inla.CRS("sphere")) ## 3 vertices found in the polygon inla.over_sp_mesh(obj, mesh, type = "vertex") ## 3 triangles found in the polygon inla.over_sp_mesh(obj, mesh) ## Multiple transformations can lead to slightly different results due to edge cases ## 4 triangles found in the polygon inla.over_sp_mesh(obj, inla.spTransform(mesh, CRSobj = inla.CRS("mollweide_norm")), ignore.CRS = FALSE) ## Ignoring mismatching coordinate systems is rarely useful ## 20 triangles "found in" the polygon inla.over_sp_mesh(obj, inla.spTransform(mesh, CRSobj = inla.CRS("mollweide_norm")), ignore.CRS = TRUE)