so for those who are interested, I found a partial solution in how VoronoiDelaunay plots meshes. Basically, one can interrupt lines and polygons using NaNs. For example:
points = Tuple{Float64, Float64}[] for n = 1:size(T, 2) p = [X[:, T[:, n]] X[:, T[1,n]]] for m = 1:size(p, 2) push!(points, tuple(p[1, m], p[2, m])) end push!(points, tuple(NaN, NaN)) end compose(context(units=ub), polygon(points), fill(elcol), stroke(linecol), linewidth(lwidth) ) This works well for me - for now. Though I am sure one can improve this further. Christoph