The silicate package provides tools and data models to work with decompositions of complex data types. The primary examples of use are with spatial data layers and most silicate workflows are geared to using the underlying structures and topology of spatial data. We can use these for a range of activities:
The functions in silicate are grouped into models and worker verbs. The models have capitalized names and the worker verbs are prefixed with ‘sc’.
The worker verbs work with all (there are some exceptions) model types and with several external formats, particularly sf, sp, trip, and other external package types.
Topology is the inherent underlying shape structure of data and silicate is clear on the distinction. There are two kinds of data structures that are abandoned by geospatial vector formats: surfaces and ordered linear paths. This might sound controversial but really is not, it’s a trivial phrasing of what is in the simple features standard. This is why silicate treats “spatial polygons” as PATH forms, because they are inherently linear structures and any surface interpretation is purely implicit (winding or even-odd rule). For an explicit surface, silicate requires a truly surface-based form, with triangular primitives.
Silicate itself only inclues ear-clipping triangulation, a
relatively cheap method suitable for implicit surfaces using closed
paths. A more sophisticated meshing triangulation is available in
the anglr package. The
anglr package is less strict than silicate, and will happily treat many
path-based data structures as implicit surfaces, and so we can
as.mesh3d() a polygon or line
layer and it will do its best to generate a true explicit surface. We
can otherwise use the silicate model
TRI0() or the anglr models
DEL0() to actually work with the data structure in use (not
just a visual side-effect).
It’s important to know what is the simplest vs more complex ways of
creating surfaces. In general,
TRI0() is cheapest, and
DEL() is most expensive - but the former requires
paths, and the latter works with edges. The cheapest
one, is TRI0(). TRI() is next expensive (adds visible property, and
native orientation, and arbitrary subsetting) … but cheap triangles,
next is DEL0() (control over area, angle, Steiner points, Delaunay
constraint, etc), then DEL().