Description Usage Arguments Value Author(s) Examples
Wrapper for the over
method to find triangle centroids
or vertices inside sp
polygon objects
1 2 | 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
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | # 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"))
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")), 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")), ignore.CRS = TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.