A geometrical characterization of regions of uniqueness and applications to discrete tomography

Abstract

In the reconstruction problem of Discrete Tomography, projections are considered from a finite set S of lattice directions. Employing a limited number of projections implies that the injectivity of the Radon transform is lost, and, in general, images consistent with a given set of projections form a huge class. In order to lower the number of allowed solutions, one usually tries to include in the problem some a priori information. This suggests that modeling the tomographic reconstruction problem as a linear system of equations is preferable.In this paper we propose to restrict the usual notion of uniqueness, related to the solutions of the linear system, and to provide, for each set S, a geometrical characterization of the shape of a lattice subset, say region of uniqueness (ROU), forming a partial, fast reconstructible, solution.Any selected set S intrinsically determines its ROU inside an arbitrary lattice grid. For instance, trivially, if |S| = 1, the ROU is represented by two rectangles having sizes equal to the absolute values of the entries of the unique direction in S, and placed at two opposite corners of the chosen grid. Surprisingly, if |S| = 2, the problem becomes much more complicated.Our purpose is to provide a geometrical characterization of the ROU. This is based on a double Euclidean division algorithm (DEDA), which runs in polynomial time. It turns out that the ROU is delimited by a zigzag profile obtained by means of numerical relations among the entries of the employed directions. According to different inputs in DEDA, the shape of the ROU can change consistently, as it can be easily observed from the provided examples. Moreover, after selecting a region of interest (ROI) from a given phantom, we exploit DEDA to reconstruct the part of the ROI which falls in the ROU and, with a few further a priori knowledge, even parts of the ROI which are outside the ROU.

Geometrical Characterization of the Uniqueness Regions Under Special Sets of Three Directions in Discrete Tomography

“…For |S| = 2, the ROU has been completely characterized in [5]. We now aim to extend this result to sets of three directions.…”

“…In this case, quick means in linear time (we recall that the reconstruction of the whole grid is NP-hard for more than two directions). This question involves the definition of (geometric) region of uniqueness (geometric ROU), which has already been presented in [4,5]. Definition 1.…”

“…As a consequence, a suitable selection of the projections' directions may allow the reconstruction of the whole object. In [5] we have characterized the shape of subsets of pixels in a generic area, whose values are uniquely determined from two projections, without explicitly computing them. These regions are called Regions of Uniqueness, briefly ROU.…”

“…In Sect. 2 the notation is set, together with the main definitions and a summary of the algorithm proposed in [5]. In Sect.…”

Stability results for sets of uniqueness in binary tomography

“…Such images are known as ghosts. The presence of ghosts was the input to investigate which parts of the object could be anyway reconstructed; this problem has been developed in [2][3][4].…”