a 3d image can be interpreted as a cellular complex. In this context, minimum-weight
surfaces have been used to segment flat structures in medical 3d images [15].
The segmentation of cracks in images of concrete is a broad field of research. A large
variety of segmentation methods have been studied for 2d and 3d images. Comprehensive
overviews and comparison studies in 2d and 3d can be found in [9] and [4], respectively.
In image segmentation, ground truths are necessary for an objective output evaluation.
Furthermore, they are prerequisite for training machine learning models such as convo-
lutional neural networks [8] and random forests [7]. In 2d, annotated data is abundantly
available, for example SDNET2018 [11]. For 3d images, manual annotation is not feas-
ible due to the large amount of data. Therefore, 3d ground truth images are scarce and
one has to rely on artificial data. Previously, the generation of synthetic cracks has been
realized via fractional Brownian surfaces [1, 4]. However, the Hurst index - a measure of
roughness - is the only parameter that can be used to control the shape of the surface.
Therefore, this model is rather limited.
Voronoi diagrams have been used previously for simulating crack propagation via finite
element methods [13, 16]. In this work, we propose a novel method to synthesize crack
structures in 3d images using Voronoi diagrams generated from random point processes.
The method includes two aspects:
First, we use minimum-weight surfaces as a tool to model the macrostructure of cracks.
Bounded 3d Voronoi diagrams serve as the underlying cell complexes. This leaves several
degrees of freedom such as the choice of the generating point processes, bounding cycles
and the facet weights.
Second, the computed surfaces are discretized to 3d binary images. The rough micro-
structure of cracks is modeled via a second Voronoi diagram on a finer scale. Then, the
cracks are embedded into patches of real computed tomography (CT) images of concrete.
This paper is structured as follows. In Section 2, we outline the concept of shortest
paths. Their extension to minimum-weight surfaces is described in Section 3. Section 4
comprises our crack modeling pipeline using Voronoi diagrams. It includes a description
of our approach for macro- and microstructure modeling. Section 5 serves as conclusion
and outlook to possible future research.
2 Shortest paths
Let G= (V, A) be a directed graph consisting of a set of vertices Vand directed arcs
A⊆V×V. We use the notation α(a) for the start vertex and ω(a) for the end vertex
of an arc a∈A.
Furthermore, let c:A→R>0be a function assigning a non-negative weight to every arc
in A. A path in Gis a finite sequence of vertices and arcs, P= (v0, a0, v1, a1, . . . , ak−1, vk),
k≥0, with vi∈Vand ai∈Awith α(ai) = viand ω(ai) = vi+1 for i= 0, . . . , k −1. It
is called a cycle (or closed contour) if no arc and no vertex is included more than once
except for v0=vk. The weight of a path Pis given as c(P) = Pk−1
i=0 c(ai).
Given two vertices s, t ∈V, the SPP is looking for a path of minimum weight from sto
t. Many algorithms for solving the SPP exist, for example Dijkstra’s [10] or Bellman’s
2