Abstract
A finite subset K of ℤd is said to be lattice-convex if K is the intersection of ℤd with a convex set. The covariogram g K of K⊆ℤd is the function associating to each u∈ℤd the cardinality of K∩(K+u). Daurat, Gérard, and Nivat and independently Gardner, Gronchi, and Zong raised the problem of the reconstruction of lattice-convex sets K from g K . We provide a partial positive answer to this problem by showing that for d=2 and under mild extra assumptions, g K determines K up to translations and reflections. As a complement to the theorem on reconstruction we also extend the known counterexamples (i.e., planar lattice-convex sets which are not reconstructible, up to translations and reflections) to an infinite family of counterexamples.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In the research area of (geometric) tomography one is usually concerned with the reconstruction of an object from information on the projections and/or sections of that object. Discrete tomography is concerned with the reconstruction of discrete objects and relies on discrete tools (see [18–20]) while analytic tomography is concerned with the reconstruction of objects determined by analytic data (see [16]). Both research directions are naturally unified within the realm of convex geometry, since convex sets can be of both discrete and analytic nature. Geometric tomography has numerous applications outside mathematics. E.g., one of the many links between tomography and physics is the so-called phase retrieval problem, that is, the problem of retrieving an unknown object from its diffraction data (see [22]). The diffraction data can be described in mathematical language as an autocorrelation or a covariogram.
In this paper we present results on a discrete version of the covariogram problem under convexity assumptions. For presenting the problem we first start with the (standard) covariogram problem in the continuous case. Let d∈ℕ. By \(\operatorname{vol}\) we denote the d-dimensional volume. Let K⊆ℝd be a body, i.e., a bounded set coinciding with the closure of its interior. The (continuous) covariogram of K is the function g K :ℝd→ℝ that sends u to \(\operatorname{vol}(K \cap(K+u))\). The covariogram problem is the problem of determining K from g K within a given class of bodies. Clearly, g K remains unchanged if we translate K or reflect K with respect to a point (in what follows, the term ‘reflection’ will always mean a point reflection). Thus, the determination of K up to translations and reflections is considered the unique determination. It is not hard to show that g K is the autocorrelation of the characteristic function of K. By this, the covariogram problem is a special case of the phase retrieval problem. Recently there has been much progress on the retrieval of K from g K within the class of convex bodies (this version of the covariogram problem is usually called Matheron’s problem); see [1, 3, 7, 8, 11, 12, 26].
If K is a finite subset of ℤd we define the (discrete) covariogram of K to be g K :ℤd→ℤ, u↦|K∩(K+u)|; here |⋅| stands for the cardinality. The problem of retrieving K from g K is called the discrete covariogram problem. The discrete covariogram problem (for d=1 also known as the partial digest problem) was studied by many authors; see for example [4, 24, 27]. It turns out that the discrete covariogram problem is a special case of the continuous covariogram problem. In fact, for K⊆ℤd, g K determines \(g_{K+[0,1]^{d}}\) and vice versa.
A subset K of ℤd is said to be lattice-convex if K is the intersection of ℤd with a convex set. By \(\mathcal{K}(\mathbb{Z}^{d})\) we denote the class of finite, lattice-convex subsets of ℤd which affinely span ℝd. To the best of our knowledge, the problem of determination of K from g K within \(\mathcal{K}(\mathbb {Z}^{d})\) was first addressed in [14] and [18]. Both papers also give counterexamples to the uniqueness of reconstruction. In our paper we contribute to the study of this problem. In view of the previous observations, on the one hand, this problem is a special case of the continuous covariogram problem and, on the other hand, a discrete version of Matheron’s problem. We emphasize that recently the covariogram problem for lattice-convex sets attracted attention of researchers working in X-ray crystallography and the theory of quasicrystals; see [6]. (We refer to [21] and [25] for further information on the latter topics.) The discrete covariogram problem is also related to the discrete X-ray problem (i.e., the problem of reconstruction from discrete X-ray pictures). In [18, Theorem 4.1] it was proved that sets \(K, L \in\mathcal{K}(\mathbb{Z}^{2})\) have equal covariograms if and only if for every direction the discrete X-ray pictures of K and L are (so-called) rearrangements of each other. Strong positive results on the discrete X-ray problem for the class \(\mathcal{K}(\mathbb{Z}^{2})\) were obtained in [17].
In this paper we address the covariogram problem for \(\mathcal {K}(\mathbb{Z}^{2})\) and provide both positive and negative results. One can see that for \(K \in \mathcal{K}(\mathbb{Z}^{2})\) the directions of the (discrete) edges of K are determined by g K . Our positive result (Theorem 2.1) implies that, for a given collection of prescribed edge directions, a set \(K \in\mathcal{K}(\mathbb{Z}^{2})\) is determined by g K if all edges of K have sufficiently large cardinality and the difference between cardinalities of parallel edges is either zero or sufficiently large. The proof of the positive result suggests that, under the given sufficient conditions, the retrieval of K from g K is analogous to the retrieval of polygons from the continuous covariogram (the problem solved in [26], see also [8] for a short solution). Thus, the discrete nature of the problem is only revealed when the above sufficient conditions are violated. As examples in [18, p. 402] and [14, Sect. 4] show, \(K \in\mathcal{K}(\mathbb{Z}^{2})\) is not always determined by g K . Therefore, as a complement to the positive result, we also study those \(K \in \mathcal{K}(\mathbb{Z}^{2})\) which are not determined by g K . All known \(K \in\mathcal{K}(\mathbb{Z}^{2})\) not determined by g K can be represented as a direct sum K=S⊕T, where S⊆ℤ2 and \(T \in\mathcal {K}(\mathbb{Z}^{2})\) are not centrally symmetric and T can be covered by two distinct parallel lines. Thus, we are concerned with the characterization of all K having the above form. It turns out that the class of all such K is quite narrow (the convex hull of S is a certain polygon with at most six edges and T is also of a certain particular form); this is our second main result which is contained in Theorem 2.5. The complete classification of all \(K\in\mathcal{K}(\mathbb{Z}^{2})\) determined by g K remains an open question.
The paper is organized as follows. In Sect. 2 we introduce necessary notation, present the main results (Theorems 2.1 and 2.5), and state consequences and open questions. Theorems 2.1 and 2.5 are proved in Sects. 3 and 4, respectively.
2 Results and Open Questions
For formulating our main results we need further notations and terminology. The origin in ℝd will be denoted by o. If u∈ℤd∖{o}, then gcd(u) stands for the greatest common divisor of the components of u. For u 1,u 2∈ℤ2 we denote by det(u 1,u 2) the determinant of the matrix with columns u 1 and u 2 (in this order). The interior, boundary, affine hull, and conical hull operations will be denoted by \(\operatorname {int}\), \(\operatorname{bd}\), \(\operatorname{aff}\), and \(\operatorname{cone}\), respectively. The convex hull of K⊆ℝd will be denoted by \(\overline{K}\).Footnote 1 The closure of the set {u∈ℝd:g K (u)>0} is called the support of g K and is denoted by \(\operatorname{supp}g_{K}\). As usual, for S,T⊆ℝd the (Minkowski) sum S+T is defined to be {s+t:s∈S ∧ t∈T}; for t∈ℝd we also write S+t instead of S+{t}. Further, we put −T:={−t:t∈T}. The notations t−S and S−t are introduced in a similar fashion. For R⊆ℝ and S⊆ℝd we use the notation RS={rs:r∈R ∧ s∈S}.
For K⊆ℝd we put [K]:={K+t:t∈ℝd}∪{−K+t:t∈ℝd}, which is the set of translations and reflections of K. A pair (K,L) of subsets of ℤd is said to be homometric if g K =g L and nontrivially homometric if additionally [K]≠[L].
Let us introduce the notion of lattice-convexity with respect to an arbitrary lattice. A subset K of a lattice \(\mathbb{L}\subseteq \mathbb{R}^{d}\) is said to be lattice-convex (with respect to \(\mathbb{L}\)) if K is the intersection of \(\mathbb{L}\) and a convex set. If \(\mathbb{L}\) is clear from the context, we will drop the qualifier ‘with respect to \(\mathbb {L}\)’. By \(\mathcal{K}(\mathbb{L})\) we denote the class of finite, lattice-convex subsets K of \(\mathbb{L}\) with \(\operatorname{aff}K= \mathbb{R}^{d}\).
We use a few standard notations and definitions of convex geometry (in some of the cases, carried over to the case of lattice-convex sets); for details see [28]. If K is a subset of ℝd, the support function h(K,⋅):ℝd→ℝ∪{∞} of K is defined by h(K,u):=sup{〈x,u〉:x∈K}, where 〈⋅,⋅〉 is the standard scalar product in ℝd. If K⊆ℝd is compact, then the support set of K in direction u∈ℝd is defined by
Clearly, if \(K \in\mathcal{K}(\mathbb{Z}^{d})\) and u∈ℝd, then F(K,u) is lattice-convex. For support sets we have
for bounded compact sets K,L⊆ℝd and each u∈ℝd.
Some terminology from classical convexity theory can be directly carried over to the setting of lattice-convex sets. If \(K \in\mathcal {K}(\mathbb{Z}^{2})\), then a support set F(K,u) with u∈ℝ2∖{o} and that contains exactly one element will be called a vertex of K; a support set F(K,u) with u∈ℝ2∖{o} and with more than one element will be called an edge of K with outer normal u. We introduce
Vectors u∈ℤ2∖{o} with gcd(u)=1 are called primitive.
To measure the number of lattice points on the edges and the difference of parallel edges of \(K\in\mathcal{K}(\mathbb{Z}^{2})\) we introduce
where as usual min∅=∞. Observe that m″(K)=∞ if and only if K is centrally symmetric, and that m(K)=m′(K) if K has no parallel edges of different lengths.
If U is a finite set of vectors in ℝ2 which linearly spans ℝ2, then the discrepancy δ(U) of U is
where D(U):={|det(u 1,u 2)|:u 1,u 2∈U}∖{0}. For \(K\in\mathcal{K}(\mathbb{Z}^{2})\) we define δ(K):=δ(U(K)). The geometric meaning of the discrepancy of a system of vectors can be easily visualized in terms of cardinalities of lattice-convex sets. For two primitive vectors u,v∈ℤ2 the number |det(u,v)| is equal to the number of lattice points in the half-open parallelogram [0,1)u+[0,1)v; this can be derived, e.g., from Pick’s formula. In this sense, for primitive u and v, |det(u,v)| is a measure for ‘how densely’ \(\operatorname{cone}\{u,v\}\) is filled with lattice points, relatively to ‘how densely’ the border of this cone is filled. The discrepancy δ(U) then can be understood as a number measuring the ‘unbalancedness’ of this kind of density of the cones that can be formed by pairs of nonparallel vectors from U.
Now we are ready to formulate our first main result.
Theorem 2.1
Let \(K, L \in\mathcal{K}(\mathbb{Z}^{2})\) be sets with g K =g L . Then m′(K), m″(K), m(K), U(K)∪U(−K) and δ(K) are determined by g K . Moreover, if m(K)≥δ(K)2+δ(K)+1, then [K]=[L]; hence g K determines K within the class \(\mathcal {K}(\mathbb{Z}^{2})\) up to translations and reflections.
The proof of Theorem 2.1 will be given in Sect. 3. It is based on the proof ideas of results for continuous covariograms presented in [1, 2, 8, 13]. Most parts of our proof of Theorem 2.1 can be rendered into an algorithmic form. A detailed analysis of algorithmic issues would certainly be interesting and valuable, but it would go beyond the scope of this paper.
Remark 2.2
Note that δ(K)≤2n 2 whenever U(K)⊆{−n,…,n}2 and n∈ℕ. Thus, the assumption on m(K) in Theorem 2.1 can be replaced by the less technical but somewhat stronger assumption m(K)≥4n 4+2n 2+1.
Since examples in [14] and [18] show that \(K \in\mathcal{K}(\mathbb {Z}^{2})\) is not always reconstructible from g K , we also study nonreconstructible sets \(K\in\mathcal{K}(\mathbb{Z}^{2})\). That is, we work towards a description of nontrivially homometric pairs of sets lying in \(\mathcal{K}(\mathbb{Z}^{2})\).
We first present known information on the related question on the structure of homometric pairs of finite subsets of ℤd. In [27] it is shown that homometry for subsets of ℤd can be explained by considering the reconstruction from the autocorrelation within the group ring \(\mathbb{F}[\mathbb{Z}^{d}]\), where \(\mathbb{F}\) is ℤ, ℝ, or ℂ. Let us give more details. Let Z:=Z 1,…,Z d be a sequence of independent commuting indeterminates. Given u=(u 1,…,u d )∈ℤd, by Z u we denote the expression \(Z_{1}^{u_{1}} \cdots Z_{d}^{u_{d}}\). Then \(\mathbb{F}[\mathbb{Z}^{d}]\) is the ring of all expressions \(f :=\sum_{u \in \mathbb{Z}^{d}} a_{u} Z^{u}\) such that \(a_{u} \in\mathbb{F}\) is equal to zero for all but finitely many u; the ring addition and multiplication are introduced in the same way as in polynomial rings. For f as above we define \(f^{\ast}:= \sum_{u \in\mathbb{Z}^{d}} \overline{a_{u}} Z^{-u}\), where \(\overline{a_{u}}\) is the complex conjugate of a u . The autocorrelation of f is then defined to be ff ∗. In [27, Theorem 2.3] it is shown that \(f_{1}, f_{2} \in\mathbb {F}[\mathbb{Z}^{d}]\) are homometric (i.e., have equal autocorrelations) if and only if f 1=Z u h 1 h 2 and \(f_{2} = c Z^{v} h_{1} h_{2}^{\ast}\) for some u,v∈ℤd, \(c \in\mathbb{F}\) with |c|=1 and \(h_{1}, h_{2} \in \mathbb{F}[\mathbb{Z}^{d}]\). (The proof is based on the fact that \(\mathbb {F}[\mathbb{Z}^{d}]\) is a unique factorization domain.) Since every finite subset K of ℤd can be identified with its generating function f K :=∑ u∈K Z u and since \(f_{K} f_{K}^{\ast}\) provides the same information as g K , the latter result can also be used for homometric pairs of finite subsets of ℤd.
As the problem of reconstruction of a finite K⊆ℤd from g K is more geometric than the analogous problem for \(\mathbb{F}[\mathbb{Z}^{d}]\), one might expect a more geometric reason of homometry. This is sometimes the case. The explanation will require the notion of direct (Minkowski) sums. We call the Minkowski sum K=S+T of two lattice sets S,T∈ℤd direct, if for each s,s′∈S and t,t′∈T the equality s+t=s′+t′ implies s=s′ and t=t′. In this case we write K=S⊕T. It is easy to see that if the sum of S and T is direct, then the sum of S and −T is also direct. One also readily verifies that if K,L∈ℤd satisfy K=S⊕T and L=S⊕(−T) for suitable S,T∈ℤd, then g K =g L . Moreover, in this case one has [K]≠[L] if and only if both S and T are not centrally symmetric. So direct Minkowski sums provide a way to construct examples of pairs of nontrivially homometric lattice sets. We say that the homometric pairs (K,L) obtained this way are geometrically constructible. For K,L,S,T as above the relations K=S⊕T and L=S⊕(−T) carry over to f K =f S⊕T =f S f T and \(f_{S \oplus (-T)} = f_{S} f_{-T} = f_{S} f_{T}^{\ast}\). This provides a link to the result for \(\mathbb{F}[\mathbb{Z}^{d}]\). In [27] a homometric pair of subsets of ℤd (for d=1) was given which is not geometrically constructible. This justifies the extension from the class of finite subsets of ℤd to \(\mathbb{F}[\mathbb{Z}^{d}]\).
We wish to study the retrieval of K from g K for \(K \in\mathcal {K}(\mathbb{Z}^{d})\). Taking into account the presented information on \(\mathbb{F}[\mathbb{Z}^{d}]\), the following problem can be formulated:
Problem 2.3
Let d≥2 and let \(\mathbb{F}\) be ℤ, ℝ or ℂ. Describe all \(h \in\mathbb{F}[\mathbb{Z}^{d}]\) such that h is a factor of f K for some \(K \in\mathcal{K}(\mathbb{Z}^{d})\).
Let us briefly discuss the analogous problem for d=1. The problem resembles to the study of the factors of the polynomial \(\sum _{k=0}^{n-1}Z^{k}=(Z^{n}-1)/(Z-1)\) for some n∈ℕ. For \(\mathbb{F}=\mathbb{C}\), the irreducible factors correspond to just the roots of unity (except 1), and for \(\mathbb{F}=\mathbb{R}\) one can describe the factors with the results for the complex case. For \(\mathbb{F}=\mathbb{Z}\) the problem is well-studied, see [23, Ch. VI, §3]; the corresponding factors can be determined with the help of the so-called cyclotomic polynomials.
Also, the following question arises naturally:
Question 2.4
Do there exist nontrivially homometric pairs (K,L) of sets in \(\mathcal{K}(\mathbb{Z}^{2})\) which are not geometrically constructible?
We performed an exhaustive computer search which showed that all homometric pairs of sets \(\mathcal{K}(\mathbb{Z}^{2})\) which are contained in {1,…,6}×{1,…,5} are geometrically constructible. Our computer search also covered the homometric pairs presented in [14] and [18], up to affine automorphisms of ℤ2. The fact that the example in [18] is geometrically constructible was noticed in [5, pp. 280–281].
Another natural research direction is the study of geometrically constructible homometric pairs from \(\mathcal{K}(\mathbb{Z}^{2})\). A partial information is provided by Theorem 2.5 given below. A set \(T \in\mathcal{K}(\mathbb{Z}^{2})\) is said to have lattice width 1 if T can be covered by two distinct parallel lines. Theorem 2.5 characterizes lattice-convex direct sums S⊕T with S⊆ℤ2 and with \(T\in \mathcal{K}(\mathbb{Z}^{2})\) having lattice width 1.
Theorem 2.5
Let k,ℓ be integers with k>ℓ≥0. We define T:=({0,…,k}×{0})∪({0,…,ℓ}×{1}), the vectors w 1:=(−k−1,1), w 2:=(ℓ+1,1), and the lattice \(\mathbb{L}:=\mathbb{Z}w_{1}+\mathbb{Z}w_{2}\). Let S be a set with o∈S⊆ℤ2. Then the following conditions are equivalent:
-
(i)
The sum of S and T is direct and lattice-convex.
-
(ii)
S is lattice-convex with respect to \(\mathbb{L}\) and \(\overline{S}\) is a polyhedron (in ℝ2) such that
-
every edge of \(\overline{S}\) is parallel to w 1 or w 2 if k>ℓ+1 and
-
every edge of \(\overline{S}\) is parallel to w 1, w 2, or w 1+w 2 if k=ℓ+1.
-
The assumption o∈S is not really restrictive, since Condition (i) is invariant under translations of S. Thus, assuming o∈S, Condition (ii) can be formulated in a less technical form. Condition (ii) is invariant under the change of S by −S. Thus, if (ii) holds, then S⊕T and S⊕(−T)=−((−S)⊕T) are both lattice-convex. If k>ℓ+1 and S is a finite set satisfying Condition (ii), then \(\overline{S}\) is a parallelogram or a line segment (and by this, centrally symmetric). Thus, as a direct consequence of Theorem 2.5 we obtain:
Corollary 2.6
Let S, T, k and ℓ be as in the formulation of Theorem 2.5, let S be finite and the sum of S and T be direct. Then the following conditions are equivalent:
-
(i)
(S⊕T,S⊕(−T)) is a nontrivially homometric pair of sets in \(\mathcal{K}(\mathbb{Z}^{2})\).
-
(ii)
All of the following properties are satisfied:
-
k=ℓ+1,
-
S is lattice-convex with respect to \(\mathbb{L}\),
-
\(\overline{S}\) is a polyhedron (in ℝ2) with every edge parallel to w 1, w 2, or w 1+w 2,
-
S is not centrally symmetric.
-
Corollary 2.6 describes (up to affine automorphisms of ℤd) all pairs of nontrivially homometric sets that can be written as (S⊕T,S⊕(−T)) with a lattice-convex summand T of lattice width 1. Thus, on the one hand, Corollary 2.6 shows that there are infinitely many nontrivially homometric pairs of sets in \(\mathcal{K}(\mathbb{Z}^{2})\) (see Fig. 1 for an illustration). On the other hand, under the given restrictions on T, these homometric pairs are ‘sporadic’ in the sense that \(\overline{S}\) has at most six edges and very restricted normal vectors and the shape of T is also very specific (see the condition k=ℓ+1 in Corollary 2.6).
We note that, up to affine automorphisms of ℤ2, all homometric pairs found by our computer search are covered by Corollary 2.6. Thus, we ask the following question:
Question 2.7
Do there exist nontrivially homometric pairs which are geometrically constructible and do not coincide (up to affine automorphisms of ℤ2) with the homometric pairs in Corollary 2.6?
Remark 2.8
In [6] Baake and Grimm used the example from [18] to construct a pair of ‘distinct homometric model sets’. Loosely speaking, these are mathematical descriptions of quasicrystals that differ from each other in many significant ways, but that cannot be distinguished by their diffraction data. Our Corollary 2.6 can be used to extend the result of [6] to an infinite series of pairs of distinct homometric model sets.
3 Proof of Theorem 2.1
The proof of Theorem 2.1 is organized as follows. First we show that g K contains (up to translations and reflections) the ‘local boundary information’ on K, namely edges (Lemma 3.1) and normal cones (Lemma 3.3). These arguments are discrete equivalents of some of the arguments given in [8]. The key arguments contained in Lemma 3.5 show that the ‘local boundary information’ can be assembled in the unique way to reproduce the boundary of \(\overline {K}\) up to translations and reflections. One the one hand, the proof of Lemma 3.5 is related to arguments in [8] and, on the other hand, to the arguments on the so-called capturing arcs used in [13] and [1] (see also [2] which is a full-length version of [1]). Somewhat more precisely, if we have a pair of arcs of \(\overline{K}\) which is symmetric (with respect to reflection) and this pair is adjacent to a nonsymmetric pair of arcs of \(\overline{K}\), then, using certain values of g K , the symmetric and the nonsymmetric pair can be ‘glued together’. This ‘gluing procedure’ is analogous to a continuous procedure used in [1, Lemma 3 and the proof of Theorem 1] (see also [2, the proofs of Lemma 9 and Theorem 1]). The mentioned arguments from [1] are refined versions of the arguments from [13, Definition 4.1, Lemma 4.2, Proposition 4.3]. Of course, in our setting, we have to adjust the arguments of [8, 13] and [1] (or [2]) to match the discrete case.
If K⊆ℝ2, then DK:=K+(−K) is said to be the difference set of K. If K is a convex body, then DK is called the difference body of K. It is known and easy to show that for each finite K⊆ℤd one has \(\operatorname {supp}g_{K} = DK\) and \(\operatorname{supp}g_{\overline{K}} = D\overline{K}=\overline{DK}\). If v is a vertex of K, we define the normal cone N(K,v):={u∈ℝd:v∈F(K,u)} and the supporting cone \(S(K,v) :=\operatorname{cone}(K-v)\) of K at v (see also [28, p. 70]).
Lemma 3.1
Let \(K,L \in\mathcal{K}(\mathbb{Z}^{2})\) with g K =g L . Then for each u∈ℤ2∖{o} there exists L′∈[L] satisfying
In particular, m′(K)=m′(L), m″(K)=m″(L), m(K)=m(L), U(K)∪U(−K)=U(L)∪U(−L), and δ(K)=δ(L).
Proof
Since \(DK = \operatorname{supp}g_{K}\), the set DK is determined by g K . We reconstruct F(K,u)∪F(K,−u) by analyzing the values of g K on \(\operatorname{aff}F(DK,u) \cap\mathbb{Z}^{2}\).
Case 1: F(DK,u) is an edge of DK. Then at least one of F(K,u) and F(K,−u) is an edge of K, see (1). We consider a nonsingular affine transformation τ mapping ℤ onto \(\operatorname{aff}F(DK,u) \cap\mathbb{Z}^{2}\) (that is, we enumerate the elements of \(\operatorname{aff}F(DK,u) \cap\mathbb{Z}^{2}\) by elements of ℤ in a consecutive way). Then there exist k 1,…,k 4∈ℤ with k 1≤k 2≤k 3≤k 4 such that g K (τ(k))=0 for k<k 1 and k>k 4 and g K (τ(k)) is the maximum of g K on \(\operatorname{aff}F(DK,u) \cap\mathbb{Z}^{2}\) for k 2≤k≤k 3. It can be seen that F(K,u)∪F(K,−u) is a translation or a reflection of {τ(k)−τ(k 2):k 1≤k≤k 2}∪{τ(k):k 1≤k≤k 3}, see Fig. 2 for an illustration.
Case 2: F(DK,u) is a vertex of DK. Then both F(K,u) and F(K,−u) are vertices of K. In analogy to the argumentation above one shows that F(K,u)∪F(K,−u) is a translation of {o}∪F(DK,u). □
In view of Lemma 3.1 we see that for \(K \in \mathcal{K}(\mathbb{Z}^{2})\) the parameters m′(K), m″(K), m(K) and set U(K)∪U(−K) are determined by g K , which proves the first statement of Theorem 2.1.
Lemma 3.2 given below will be used to ensure that certain triangles are contained in homometric lattice-convex sets.
Lemma 3.2
Let \(K \in\mathcal{K}(\mathbb{Z}^{2})\) and let m:=m(K) and δ:=δ(K) satisfy m≥δ 2+δ+1. Then there exists k∈ℕ such that δ≤k≤(m−1)/δ. Moreover, if u,v,w∈U(K)∪U(−K) are mutually nonparallel, then ku=αv+βw for some α,β∈ℝ that satisfy 1≤|α|≤m−1 and 1≤|β|≤m−1. The same holds if u,v,w∈ℤ2∖{o} are mutually nonparallel, primitive vectors that are parallel to edges of \(\overline{K}\).
Proof
The first assertion of Lemma 3.2 follows from (m−1)/δ−δ≥1. Let u,v,w∈U(K)∪U(−K) be mutually nonparallel. There exist α,β∈ℝ∖{0} with ku=αv+βw. Employing Cramer’s rule, we get
Using the definition of δ and the choice of k, this implies that 1≤k/δ≤|α|≤kδ≤m−1 and, similarly, 1≤k/δ≤|β|≤kδ≤m−1. If U is the set of primitive vectors parallel to edges of \(\overline {K}\), then δ(K)=δ(U). This implies the second assertion. □
Lemma 3.3
Let \(K, L \in\mathcal{K}(\mathbb{Z}^{2})\) be such that g K =g L . Let m(K)≥δ(K)2+δ(K)+1. Let v be a vertex of K and let u∈ℤ2∖{o} satisfy F(K,u)={v}. Then there exists L′∈[L] such that F(K,u)=F(L′,u), F(K,−u)=F(L′,−u), and S(K,v)=S(L′,v). Moreover, if F(K,−u)={v′} for some vertex v′ of K, then L′ can be chosen to additionally satisfy S(K,v′)=S(L′,v′).
Proof
Put m:=m(K) and δ:=δ(K).
Case 1: F(K,−u) is an edge of K. Due to Lemma 3.1 we find L′∈[L] such that F(K,u)=F(L′,u), F(K,−u)=F(L′,−u). It remains to show that S(K,v)=S(L′,v). Let w 1,…,w 4∈ℤ2 be vectors with gcd(w i )=1 for each i∈{1,…,4} such that if E is an edge of K or L′ with v∈E, then for some i∈{1,…,4} the point v+w i is the element of E∖{v} closest to v. We argue by contradiction and assume that S(K,v)≠S(L′,v). After possibly interchanging K with L′ and reordering w 1,…,w 4, we assume that \(S(K,v)= \operatorname{cone}\{w_{1},w_{2}\}\) and w 1∉S(L′,v). Let w∈ℤ2 with gcd(w)=1 be parallel to F(K,−u) and directed in such a way that \(w_{1}\in\operatorname{cone}\{w_{i},-w\}\) for i∈{2,3,4}. See Fig. 3 for an illustration of the situation and the following arguments.
By Lemma 3.2 there exists k∈ℕ with k≤m−1 and there exist real numbers α 2,α 3,α 4,β 2,β 3,β 4 such that kw 1=α i w i −β i w and 1≤α i ≤m−1 and 1≤β i ≤m−1 for each i∈{2,3,4}. Let p be the vertex of F(K,−u)=F(L′,−u) such that p+w belongs to F(K,−u)=F(L′,−u). Since α i ,β i ≤m−1 for each i∈{2,3,4} and by definition of m=m(K), the sets \(\overline{K} \cap(\overline{K}+ v - p + k w_{1})\) and \(\overline{L'}\cap(\overline{L'}+ v - p + k w_{1})\) are triangles whose edges are parallel to the vectors from {w,w 1,…,w 4}. The fact that β 3,β 4≥1 implies that v+kw 1+w is not contained in the interior of \(\overline{L'}\). So, by construction, one has
It follows that g K (v−p+kw 1+w)<g K (v−p+kw 1) and g L′(v−p+kw 1+w)=g L′(v−p+kw 1) and by this g K ≠g L′, a contradiction.
Case 2: F(K,−u)={v′} for some vertex v′ of K. The main proof idea for this case is essentially borrowed from [8, Lemma 3.1 (Case 2)]. Let us first show that S(K,v) or −S(K,v) is a supporting cone of L′ for some L′∈[L]. We do this by proving an equivalent statement, namely that N(K,v) or −N(K,v) is a normal cone of L (with respect to an appropriate vertex of L), because the arguments for the normal cones will turn out to be more natural. We distinguish the following subcases:
Case 2.1: There exists u′ such that F(K,−u′)={v} and F(K,u′) is an edge.
Case 2.2: For each u′ such that F(K,−u′)={v}, the set F(K,u′) is a vertex of K.
For Case 2.1 we see that N(K,v) or −N(K,v) is a normal cone of L (this follows directly from Case 1). For Case 2.2 observe that F(K,u′) is the same vertex for each u′ with F(K,−u′)={v}. In fact, otherwise we could find two directions u 1,u 2 such that F(K,−u i )={v} for i∈{1,2} and F(K,u i ) are consecutive vertices of K. But then the outer normal u′ of the edge joining F(K,u 1) with F(K,u 2) satisfies F(K,−u′)={v}, a contradiction to the assumption for Case 2.2. Thus, we have {v′}=F(K,u′) for each u′ with F(K,−u′)={v}. It follows that −N(K,v)⊆N(K,v′). If −N(K,v)=N(K,v′), we apply Lemma 3.1 and infer that −N(K,v) and N(K,v′) are normal cones of L, as well. Otherwise −N(K,v)⊊N(K,v′), and thus there exists a direction u′ such that F(K,−u′) is an edge containing v and F(K,u′)={v′}. From Case 2.1, it follows that N(K,v′) is determined by g K , up to a sign. That is, for some L′∈[L] the point v′ is a vertex of L′ and N(K,v′)=N(L′,v′). Then Lemma 3.1 implies that the set of those edge normals of K and L′ that lie in −N(K,v′) coincide. Consequently N(K,v)=N(L′,v).
Summarizing the arguments of the subcases, we conclude that N(K,v) or −N(K,v) is a normal cone of L. The arguments of the subcases are still valid if we replace u′ by −u′ and exchange v and v′. Thus, we also conclude that N(K,v′) or −N(K,v′) is a normal cone of L. The latter implies that either N(K,v) and N(K,v′) or −N(K,v) and −N(K,v′) are both normal cones of L (e.g., N(K,v) and −N(K,v′) cannot be simultaneously normal cones of L since their interiors intersect). Taking into account that F(K,u)∪F(K,−u) and F(L,u)∪F(L,−u) are translations of each other (see Lemma 3.1), we arrive at the assertion of Lemma 3.3 within Case 2. □
Remark 3.4
Lemma 3.1 and Lemma 3.3 already imply that if \(K, L \in\mathcal{K}(\mathbb{Z}^{2})\) satisfy g K =g L and if K (and thus L) does not contain parallel edges of different lengths, then m′(K)≥δ(K)2+δ(K)+1 ensures that [K]=[L]. This is a special case of Theorem 2.1.
If K contains parallel edges of different lengths, the information gained so far by Lemmas 3.1 and 3.3 is in general not enough to guarantee that K can be reassembled from g K up to translations and reflections. But g K provides more information. In fact, Lemma 3.5 (which is essentially the same as [8, Lemma 4.1]) tells us that if we could not reconstruct [K] from g K in a unique fashion, then for another solution L there must occur certain ‘symmetries’ in the parts where the boundaries of K and L (or translations or reflections of L) overlap. This will be the key argument in the proof of Theorem 2.1.
By \(\mathbb{S}^{1}\) we denote the unit circle {u∈ℝ2:〈u,u〉=1} in ℝ2.
Lemma 3.5
Let \(K, L\in\mathcal{K}(\mathbb{Z}^{2})\) be such that g K =g L and m(K)≥δ(K)2+δ(K)+1. Assume that there is a nonempty arc \(U\subseteq \mathbb{S}^{1}\) that satisfies
and that is maximal (with respect to inclusion) with the property (3). Let A be the inclusion-maximal arc satisfying \(G(U) \subseteq A \subseteq\operatorname{bd}\overline{K}\cap\operatorname{bd}\overline{L}\) and B be the inclusion-maximal arc satisfying \(G(-U)\subseteq B \subseteq\operatorname{bd}\overline{K} \cap \operatorname{bd}\overline{L}\). Let [K]≠[L]. Then one of the following holds: |A|=1 or |B|=1, or A and B are parallel line segments, or A is a reflection of B.
Proof
Put m:=m(K) and δ:=δ(K). We follow the lines of the proof of [8, Lemma 4.1], but replace a ‘continuous’ argument that makes use of derivatives with a discrete one. First, observe that [K]≠[L] implies that neither A nor B coincides with \(\operatorname{bd}\overline{K}\). Observe that U is closed.
If |A|=1 or |B|=1 we are done. So assume that both A and B contain more than one point. Let a 1, a 2 be the endpoints of A, and let b 1, b 2 be the endpoints of B. Choose the labeling in such a way that a 1, a 2, b 1, and b 2 are in counterclockwise order on \(\operatorname{bd}\overline{K}\). Further on, let u 1, u 2 be the first and second endpoint of U in counterclockwise order, respectively. We claim that
Thus, we need to show u i ∈N(K,a i ), u i ∈N(L,a i ), −u i ∈N(K,b i ), −u i ∈N(L,b i ) for i∈{1,2}. We prove u 1∈N(K,a 1) by contradiction; the remaining assertions can be settled in complete analogy. By construction we have F(K,u 1)∩A≠∅. Let U′ be the set of vectors \(u'\in\mathbb{S}^{1}\) which, in counterclockwise order, strictly precede u 1 and strictly follow every v∈N(K,a 1). If a 1∉F(K,u 1), then U′ is nonempty. For u′∈U′∪{u 1} we have u′∉N(K,a 1) and for u′∈U′ the set F(K,u′) is contained in the relative interior of A. By definition of A and due to (2) we have F(K,−u′)=F(L,−u′) for each u′∈U′∪{u 1}. This contradicts the maximality of U and concludes the proof of the claim.
Now, for i∈{1,2}, let c i =(a i +b i )/2 and let \(A'_{i}\) be the maximal subarc of A that starts in a i and in which each point is also present in B after reflection in c i . It is possible that \(A'_{i}\) contains just one point, in which case we say that \(A'_{i}\) is degenerate. Further on, let \(B'_{i}\) be the reflection of \(A'_{i}\) in c i . By construction we have \(a_{i}\in A'_{i}\) and \(b_{i}\in B'_{i}\). We also observe that a i and b i are vertices of K or L for i∈{1,2}, but if parallel edges of different lengths exist in K and L, they do not have to be vertices of both K and L.
Case 1: For each i∈{1,2} we have \(A'_{i}=A\) or \(B'_{i}=B\). If for some i∈{1,2} we have \(A'_{i}=A\) and \(B'_{i}=B\), then A is a reflection of B and we are done. So assume that for each i∈{1,2} we either have \(A'_{i}=A\) or \(B'_{i}=B\). Then, in view of the symmetries arising from the exchanging of A and B and/or replacing i by 3−i, it is sufficient to distinguish the following subcases:
Case 1.1: \(A'_{1} = A \ne A'_{2}\) and \(B'_{1} \ne B = B'_{2}\).
Case 1.2: \(A'_{1} = A = A'_{2}\) and \(B'_{1} \ne B \ne B'_{2}\).
In Case 1.1 we have
But this implies that A is properly contained in a translation of A, which yields a contradiction. In Case 1.2, \(B'_{1}\) and \(B'_{2}\) are two distinct translations of each other contained in B and each sharing an endpoint with B. This is only possible if \(B'_{1}, B'_{2}\) and B are parallel line segments.
Case 2: There exists i∈{1,2} such that \(A'_{i}\neq A\) and \(B'_{i}\neq B\). We prove that this case cannot occur by showing that it implies g K ≠g L , a contradiction. We show the arguments for i=1; the case i=2 can be proved in complete analogy. We refer to Fig. 4 for an illustration of the following arguments.
Let \(a'_{1}\) and \(b'_{1}\) be chosen so that \(\{a_{1},a'_{1}\}\) is the set of endpoints of \(A'_{1}\) and \(\{b_{1},b'_{1}\}\) is the set of endpoints of \(B'_{1}\). (If \(A'_{1}\) and thus \(B'_{1}\) is degenerate, then \(a_{1}=a'_{1}\) and \(b_{1}=b'_{1}\).) The point \(a'_{1}\) or \(b'_{1}\) is a vertex of both K and L, but not necessarily both of them. (Again, the latter can happen if K and L have parallel edges of different lengths.) But there is an edge E of \(\overline{K}\) or \(\overline{L}\) that contains \(a'_{1}\), that is not completely contained in \(A'_{1}\), and that is a subset of both \(\operatorname{bd}\overline{K}\) and \(\operatorname{bd}\overline {L}\). Similarly, there is an edge F of \(\overline{K}\) or \(\overline{L}\) that contains \(b'_{1}\), that is not completely contained in \(B'_{1}\), and that is a subset of both \(\operatorname{bd}\overline{K}\) and \(\operatorname{bd}\overline {L}\). By the choice of \(A'_{1}\), the edges E and F are not parallel. Observe that, by the definition of m, the sets \(E\setminus A'_{1}\) and \(F\setminus B'_{1}\) contain at least m−1 lattice points.
By definition of U there exist nonparallel edges E′ and E″ of K and L, respectively, that contain a 1 (not necessarily as a vertex) and are not completely contained in \(A'_{1}\). By the definition of m, the sets \(E'\setminus A'_{1}\) and \(E'\setminus A'_{1}\) contain at least m−1 lattice points.
Choose w 1,w 2,w 3,w∈ℤ2 such that gcd(w i )=1 for each i∈{1,…,4} and such the following conditions are fulfilled:
-
w 1, w 2, w 3, w are parallel to E, E′, E″, F, respectively,
-
\(a'_{1}+w_{1}\in\operatorname{bd}\overline{K}\setminus A_{1}'\), \(a_{1}+w_{2}\in\operatorname{bd}\overline{K}\setminus A_{1}'\), \(a_{1}+w_{3}\in\operatorname{bd}\overline {L}\setminus A_{1}'\), \(b'_{1}+w\in\operatorname{bd}\overline{K}\setminus B_{1}'\).
After possibly exchanging U with −U and/or interchanging the roles of K and L, the following inequalities are ensured:Footnote 2
This is shown as follows. By assumption for Case 2 we have det(w 1,w)≠0. Thus, possibly interchanging the roles of U and −U (and by this also the roles of A and B) we achieve det(w 1,w)>0. By construction det(w 1,w i )≥0 for each i∈{1,2}, and these inequalities are even strict. In fact, if we assume that, say, for i=2 one has det(w 1,w i )=0, then w 1=w 2 and \(\overline{K}= \overline{A_{1}' \cup B_{1}'}\). Hence \(\overline{K}\) is centrally symmetric and using Lemma 3.1 we deduce K=L, a contradiction. Let us show the nonstrict inequalities det(w i ,w)≥0 for i∈{1,2}. If \(A_{1}'\) is nondegenerate, the inequalities follow directly by construction (taking into account that \(a_{1},a_{1}',b_{1},b_{1}'\) lie in counterclockwise order on \(\operatorname{bd}\overline{K}\)). If \(A_{1}'\) is degenerate, then b 1∈F. Hence, using a 1∈E′∩E″, the relations in (4), and the definition of the vectors w,w 2,w 3, we see that 〈u 1,w i 〉≤0 for each i∈{2,3} and 〈u 1,w〉≥0. Taking into account the fact that a 1,a 2,b 1,b 2 are in counterclockwise order on \(\operatorname{bd}\overline{K}\) we obtain the inequalities det(w i ,w)≥0 for i∈{2,3}. The inequalities are strict. In fact, if we assume that, say, for i=2 one has det(w 2,w)=0, then w is parallel to \(a_{1}-b_{1}'\). The latter implies det(w 1,w)≤0, a contradiction to det(w 1,w)>0. Possibly interchanging the roles of K and L we can assume that det(w 2,w 3)>0, which can also be stated as w 3∉S(K,a 1).
By Lemma 3.2 there exists k∈ℕ with k≤m−1 and there exist real numbers α 2,α 3,β 2,β 3 such that kw 1=α i w i −β i w and 1≤|α 1|≤m−1 and 1≤|β i |≤m−1 for each i∈{2,3}. By Cramer’s rule and (6)–(8), we see that both α i and β i are positive for each i∈{1,2}. By construction \(a_{1}' + [o,k] w_{1} \in\operatorname {bd}\overline{K} \cap \operatorname{bd}\overline{L}\). Due to convexity of \(\overline{K}\) and \(\overline {L}\) we also have
Using the definition of m we conclude that the triangle with vertices a 1,a 1+kw 1,a 1+α 2 w 2 and the triangle with vertices a 1,a 1+kw 1,a 1+α 3 w 3 are both contained in \(\overline {K} \cap\overline{L}\). Now, using α 3≥1 and (9), we have
giving \(g_{K}(a'_{1}-b_{1}+k w_{1})<g_{L}(a'_{1}-b_{1}+k w_{1})\), a contradiction. □
Remark 3.6
The argument of Lemma 3.5 can be re-translated for the case of covariograms of convex polygons. This way one can establish a more discrete version of the proof of [8, Lemma 4.1] and, by this, of the main result in [8].
Proof of Theorem 2.1
This proof is essentially the same as the proof of [8, Theorem 1.1]. We repeat it here as a service to the reader.
Let \(K, L \in\mathcal{K}(\mathbb{Z}^{2})\) be sets with g K =g L and put m:=m(K) and δ:=δ(K). It remains to show that if m≥δ 2+δ+1, then [K]=[L]. To this end assume [K]≠[L]. We will prove that K and L are centrally symmetric. This implies that \(\overline{K}\) is a translation of \(\frac{1}{2} \overline{DK}\); the same is true for L since \(\overline{DK}=\overline{DL}\). This yields the result that K equals L up to translations, a contradiction.
To prove central symmetry of K and L, let v and v′ be opposite vertices of K, i.e.,
By Lemma 3.3 there exists L′∈[L] such that v and v′ are vertices of L′, N(K,v)=N(L′,v), and N(K,v′)=N(L′,v′). Now we apply Lemma 3.5 with U taken so that all the normalized vectors in N(K,v) are in U. Let A and B be defined as in the statement of Lemma 3.5. Observe that A and B both are neither singletons nor line segments, because \(\operatorname{bd}\overline{K}\) and \(\operatorname{bd}\overline{L'}\) coincide in a neighborhood of both the vertices v and v′. So Lemma 3.5 implies that A is a reflection of B.
We have shown that the normal cones of opposite edges of K and L are reflections of each other. Consequently, each edge of K (resp. L) is parallel to another edge of K (resp. L). It remains to show that all pairs of parallel edges of K (and thus L) have equal length to prove that K and L are centrally symmetric. So let \(E=\overline{\{v_{1},v_{2}\}}\) and \(F=\overline{\{v_{1}',v_{2}'\}}\) be two parallel edges of K, where v 1, v 2, \(v_{1}'\), and \(v_{2}'\) are in counterclockwise order of \(\operatorname{bd}\overline{K}\). By Lemma 3.1 we can choose L′∈[L] such that E and F are also edges of L′. Now both pairs \(v_{1},v_{1}'\) and \(v_{2},v_{2}'\) are opposite vertices of K as well as of L′. This gives
So the boundaries of \(\overline{K}\) and \(\overline{L}\) coincide also in a neighborhood of E and F. Lemma 3.5 then shows that E must be a reflection of F, and so they have the same length. □
Remark 3.7
Extending Theorem 2.1 to higher dimensions seems to be a nontrivial task. Recent results for the continuous covariogram problem such as those presented in [10] or [7] suggest that an analogue of Theorem 2.1 for higher dimensions might exist. However, at least for dimension d≥4, no straightforward generalization of Theorem 2.1 seems to be possible. In fact, for d≥4 one can borrow the construction of [9, Theorem 1.2] to obtain a vast class of nontrivially homometric pairs of lattice convex sets: If K⊆ℤℓ and L⊆ℤm (ℓ,m≥2) are both lattice-convex and not centrally symmetric, then the lattice convex sets K×L and K×(−L) in ℤℓ+m form a nontrivially homometric pair. So an essentially different approach seems to be needed to cope with higher dimensions.
4 Proof of Theorem 2.5
For the remainder of this section we assume that the situation of Theorem 2.5 is given, that is, k and ℓ are integers with k>ℓ≥0, and T=T 1∪T 2 with
Further we fix the vectors w 1:=(−k−1,1), w 2:=(ℓ+1,1) and the lattice \(\mathbb{L}:= \mathbb{Z}w_{1} + \mathbb{Z}w_{2}\).
Lemma 4.1
With \(\mathbb{L}\) and T defined as above one has \(\mathbb{Z}^{2} =\mathbb{L}\oplus T\).
Proof
It can be verified directly that the sum of T and \(\mathbb{L}\) is direct if and only if \((DT) \cap\mathbb{L}= \{o\}\). We have
Thus, we see that \((DT) \cap\mathbb{L}= \{o\}\) is in fact fulfilled. For showing \(\mathbb{Z}^{2} = \mathbb{L}\oplus T\) we compare ℤ2 and \(\mathbb{L}\oplus T\) modulo \(\mathbb{L}\). We have
The above equality \(|(T \oplus\mathbb{L}) / \mathbb{L}| = |T|\) holds since the sum of T and \(\mathbb{L}\) is direct, and the equality \(\det\mathbb{L}=|\mathbb{Z}^{2}/ \mathbb{L}|\) is a standard fact. Taking into account the inclusion \((T \oplus\mathbb{L}) / \mathbb{L}\subseteq\mathbb{Z}^{2}/ \mathbb{L}\), we obtain \((T \oplus\mathbb{L}) / \mathbb{L}= \mathbb{Z}^{2} / \mathbb {L}\). The latter implies \(T \oplus\mathbb{L}= \mathbb{Z}^{2}\). □
In this section we give a proof of Theorem 2.5. We shall use several standard graph-theoretic notions (see, e.g., [15] for a comprehensive account on graph theory). The graphs we use will all be undirected and in some cases infinite. As usual, V(G) stands for the node set of a graph G. Under a grid graph we shall understand a Cartesian product of two graph-theoretic paths (where a graph-theoretic path is allowed to be finite, infinite in one direction, or infinite in two directions). Given a set S⊆ℤ2 we define an undirected graph G S as follows. We declare two nodes s,s′∈S to be adjacent if and only if s−s′∈{±w 1,±w 2} for k−ℓ>1 and if and only if s−s′∈{±w 1,±w 2,±(w 1+w 2)} for k−ℓ=1. Thus, for k−ℓ>1, \(G_{\mathbb{L}}\) is an infinite grid graph and G S is the subgraph of \(G_{\mathbb{L}}\) induced by the node set S (where \(S \subseteq\mathbb{L}\)).
Let us first give a sketch of the proof of Theorem 2.5. The lengthy part of the proof is to provide enough auxiliary information for the implication (i) ⇒ (ii) of the theorem. This is done in a number of lemmas. In a first step we prove that (i) implies that the graph G S is connected (Lemma 4.5). Having shown this, we derive that the lattice-convexity of S⊕T implies certain combinatorial conditions on G S (namely, the presence of certain pairs in S force ‘neighboring’ pairs to be also elements of S; see Lemmas 4.7 and 4.8). These conditions together with the connectedness of G S imply (ii). From this we can establish (i) ⇒ (ii); the reverse implication will follow quite straightforwardly.
When dealing with a set S⊆ℤ2 such that the sum of S and T is direct and K:=S⊕T is lattice-convex we introduce the following notations. For i∈{1,2} by \(\mathcal {T}_{i}\) we denote the set of all s+T i with s∈S. We set \(\mathcal {T}:=\mathcal{T}_{1} \cup \mathcal{T}_{2}\). Then K is a disjoint union of the elements of \(\mathcal{T}\). For h∈ℤ, the set I h :=(ℤ×{h})∩K is the horizontal section of K at the height h. The set I h is empty or a (probably infinite) interval of integer points. In the latter case it decomposes into the disjoint union of the elements of \(\mathcal{T}(h):=\{J \in\mathcal{T}:J \subseteq I_{h}\}\). We introduce a consecutive order on \(\mathcal{T}(h)\) by ordering elements from left to right. In particular, for distinct \(J_{1}, J_{2} \in\mathcal{T}(h)\), J 1 strictly precedes J 2 if J 1∪J 2 is lattice-convex and the first coordinates of the points of J 1 are smaller than the first coordinates of the points of J 2.
Lemma 4.2
Let S⊆ℤ2 be such that the sum of S and T is direct and lattice-convex. Then for each h∈ℤ, the elements of \(\mathcal{T}_{1}\) and \(\mathcal{T}_{2}\) alternate in \(\mathcal{T}(h)\) with respect to the consecutive order of \(\mathcal{T}(h)\).
Proof
We show by contradiction that neither two elements of \(\mathcal{T}_{1}\) nor two elements of \(\mathcal{T}_{2}\) can be consecutive in \(\mathcal{T}(h)\). Let s,s′∈S.
If the set s+T 2 strictly precedes s′+T 2 in \(\mathcal{T}(h)\), then (s+T)∩(s′+T)≠∅, a contradiction to the fact that the sum of S and T is direct. Thus, for each h∈ℤ the set \(\mathcal{T}(h)\) does not contain two consecutive elements of \(\mathcal{T}_{2}\).
If s+T 1 strictly precedes s′+T 1, then there are precisely k−ℓ integer points between s+T 2 and s′+T 2 within I h+1. Since no two elements of \(\mathcal{T}_{2}\) are consecutive in a horizontal section of K, we see that there must be an element of \(\mathcal{T}_{1}\) between s+T 2 and s′+T 2. Elements of \(\mathcal{T}_{1}\) have cardinality k+1, hence k−ℓ≥k+1, a contradiction to ℓ≥0. □
Lemma 4.3
Let S⊆ℤ2 be such that the sum of S and T is direct and lattice-convex. Let h∈ℤ be such that \(\mathcal{T}(h)\) consists of precisely one element s+T 2 with s∈S and \(\mathcal{T}(h+1)\) consists of precisely one element s′+T 1 with s′∈S. Then k−ℓ=1 and s′−s=(−1,2)=w 1+w 2.
Proof
Let P be the convex hull of (s+T 1)∪(s′+T 1). Then P is a parallelogram with horizontal edges of length k. Therefore, the section of P by ℝ×{h} is also of length k; see Fig. 5(a) for an illustration. This section then contains at least k integer points. By lattice-convexity of S⊕T, all those points belong to S⊕T. By the assumption, I h contains precisely ℓ+1 points. Hence ℓ+1≥k. Since we assume k>ℓ, we deduce k−ℓ=1.
Clearly, s′−s=(t,2) for some t∈ℤ. If t≥0, then s+w 2 belongs to the convex hull of (s+T)∪(s′+T), a contradiction to the fact that s+T 2 is the intersection of K with ℤ×{h}. If t<−1, then s+(−1,1) belongs to the convex hull of (s+T)∪(s′+T), a contradiction the fact that s+T 2 is the intersection of K with ℤ×{h}. □
Lemma 4.4
Let S⊆ℤ2 be nonempty and such that the sum of S and T is direct and K:=S⊕T is lattice-convex. Then {h∈ℤ:I h ≠∅} is an interval of integers, i.e., a set of the form {i∈ℤ:h 1≤i≤h 2} for some h 1,h 2∈ℤ∪{±∞} satisfying h 1≤h 2.
Proof
We consider arbitrary h 1,h 2∈ℤ satisfying \(I_{h_{1}}\ne \emptyset\ne I_{h_{2}}\). It suffices to show that I h ≠0 for every h with h 1≤h≤h 2. We may assume without loss of generality that \(\mathcal{T}(h_{i})\) contains an element J i of \(\mathcal{T}_{1}\) for each i∈{1,2}; otherwise we decrease h 1 or h 2 or both by 1. If h 2≤h 1 we are done. Otherwise we consider the parallelogram \(P:=\overline{J_{1} \cup J_{2}}\), see Fig. 5(b) for an illustration. For each h∈ℤ with h 1≤h≤h 2 the set (ℝ×{h})∩P is a section of length k. Thus, (ℤ×{h})∩P contains a least k>0 integer points. Since in view of convexity we have (ℤ×{h})∩P⊆I h , we see that I h ≠∅. □
Lemma 4.5
Let S⊆ℤ2 be nonempty and such that the sum of S and T is direct and lattice-convex. Then the graph G S is connected.
Proof
If |S|=1 we are done. So let s,s′∈S be distinct. We show that they are connected by a path in G S . Choose i,i′∈{1,2} and h,h′∈ℤ such that s+T i ∈I h and s′+T i′∈I h′. Without loss of generality let h′≥h.
Case 1: h=h′. If s+T i and s′+T i′ (in this order) are consecutive in \(\mathcal{T}(h)\), then for (i,i′)=(1,2) one has s−s′=w 1 and for (i,i′)=(2,1) one has s′−s=w 2. See Fig. 5(c) for an illustration of the two cases. If s+T i and s′+T i′ are not consecutive, then using Lemma 4.2 and the previous observation for the consecutive elements of \(\mathcal{T}(h)\) lying between s+T i and s′+T i′ we see that s and s′ are connected in S.
Case 2: h′=h+1. If there exists s″∈S such that s″+T 1∈I h and s″+T 2∈I h+1, then using Case 1 we see that s and s″ are connected in S and s″ and s′ are connected in S. But then also s and s′ are connected in S. Otherwise, \(\mathcal{T}(h) \cap\mathcal{T}_{1} = \emptyset\) and \(\mathcal{T}(h+1) \cap\mathcal{T}_{2} = \emptyset\). In view of Lemma 4.2, we deduce that \(\mathcal{T}(h+1)\) consists of precisely one element of \(\mathcal {T}_{1}\) and \(\mathcal{T}(h)\) consists of precisely one element of \(\mathcal{T}_{2}\). Then, by Lemma 4.3, we see that k−ℓ=1 and s′−s=w 1+w 2; so s and s′ are connected in S.
Case 3: h′>h+1. By Lemma 4.4, I h″≠∅ for each h″ that satisfies h≤h″≤h′. Thus, we can employ the conclusion of Case 2 for those consecutive sections which are above I h and below I h′. □
Lemma 4.6
Let S⊆ℤ2 be nonempty and such that the sum of S and T is lattice-convex. Let s,s′∈S. Then the following holds.
-
I.
If s−s′=w 2−w 1, then s′−w 1,s′+w 2∈S.
-
II.
If k−ℓ>1 and s−s′=w 1+w 2, then s′+w 1,s′+w 2∈S.
-
III.
If k−ℓ=1 and s−s′=w 1+2w 2, then s′+w 2,s′+w 1+w 2∈S.
-
IV.
If k−ℓ=1 and s−s′=2w 1+w 2, then s′+w 1,s′+w 1+w 2∈S.
Proof
For an illustration of Parts I, II, and III and the according proof we refer to Fig. 5(d), 5(e), and 5(f), respectively.
Part I. After possibly translating S we assume without loss of generality that s,s′∈ℝ×{0}. Then s+T 1 and s′+T 1 are two elements of \(\mathcal{T}(0)\), and, in view of Lemma 4.2, there exists s″+T 2 with s″∈S, lying between s+T 1 and s′+T 2. Since between s+T 1 and s′+T 1 there is space for only one copy of T 2, s″+T 2 and by this s″ is determined uniquely. It follows that s″=s′−w 1. The proof that s′+w 2 uses analogous arguments for sets in \(\mathcal{T}(1)\).
Part II. Without loss of generality let s=o, so s′=(k−ℓ,−2). Let us show that (k−ℓ−1,−1)∈S⊕T. We have \((k-\ell-1,-1) = \frac{1}{2} (k- \ell-2,0 ) +\frac{1}{2} s'\), where (k−ℓ−2,0)∈s+T 1⊆s+T and s′∈s′+T. Thus, by the lattice-convexity of S⊕T, we have (k−ℓ−1,−1)∈S⊕T. The lattice point (k−ℓ−1,−1) strictly precedes s′+T 2 in I −1. Consequently, in view of Lemma 4.2, we see that s−w 2=s′+w 1∈S.
We show that (k+1,−1)∈S⊕T. We have \((k+1,-1) = \frac {1}{2} (k,0) + \frac{1}{2} (k+2,-2)\), where (k,0)∈s+T 1⊆s+T and (k+2,−2)∈s′+T 1 (since k>ℓ+1). Thus, by lattice-convexity of S⊕T, (k+1,−1)∈S⊕T. The point (k+1,−1) directly follows s′+T 2 in I −1. Thus, in view of Lemma 4.2, we see that s−w 1=s′+w 2∈S.
Part III. Without loss of generality let s′=o, so s=(k−1,3). Let P be the convex hull of (s+T 1)∪(s′+T 1). The point s″:=w 2=(k,1)=s′+w 2 lies in P, and by this, belongs to S⊕T. Since the point s″ strictly follows s′+T 2, in view of Lemma 4.2 applied for elements of \(\mathcal{T}(1)\), we deduce \(s''+T_{2} \in\mathcal{T}(1)\). Now, analogously, we use the point (k−1,2), which strictly precedes s″+T 2 and lies in P to show that s′+w 1+w 2 belongs to S.
Part IV. The proof is analogous to the proof of Part III. □
The assertions of Lemma 4.6 together with the connectedness of G S presented in Lemma 4.5 impose strong restrictions on S. This is the topic of Lemmas 4.7 and 4.8.
Lemma 4.7
Let k−ℓ>1 (and so \(G_{\mathbb{L}}\) is an infinite grid graph). Let \(S\subseteq\mathbb{L}\) be such that G S is connected and for all s,s′∈S one has
Then G S is a grid graph, i.e.,
for some α 1,α 2,β 1,β 2∈ℤ∪{−∞,∞}.
Proof
In what follows we are only interested in the case that S is finite. We give a proof for this, while the proof for the general situation is essentially the same. The following arguments are illustrated in Fig. 6(a). Let G′ be an inclusion-maximal grid graph contained in G S . If G′≠G S , in view of connectedness of G S , we can find two nodes s,s′∈ℤ2 adjacent in G S such that s∈V(G′) and s′∈S∖V(G′). Possibly translating or reflecting S, we assume that V(G′)={0,…,α}w 1+{0,…,β}w 2 for integers α,β≥0 and s′=iw 1−w 2 for some i∈{0,…,α} or s′=jw 2−w 1 for some j∈{1,…,β}. We consider the case s′=iw 1−w 2, for the other case is handled analogously. Consecutively using (10) and (11), we get i′w 1−w 2∈V(G′) for each i′∈{0,…,α}, a contradiction to the maximality of G′.
□
Lemma 4.8
Let k−ℓ=1. Let \(S \subseteq\mathbb{L}\) be such that G S is connected and for all s,s′∈S one has
Then
for some α 1,α 2,β 1,β 2,γ 1,γ 2∈ℤ∪{−∞,+∞}.
Proof
The proof of Lemma 4.8 follows along the same lines as the proof of Lemma 4.7, we give a sketch. Again one considers be the inclusion-maximal subgraph G′ of G whose node set can be given as in (15). If G′≠G S , then by Lemma 4.5 there exist s∈V(G′) and s′∈S∖V(G′) such that s and s′ are adjacent in G S . After possibly translating or reflecting S, we may assume that s−s′∈{w 1,w 2,w 1+w 2}. Successively applying (12), (13), and (14) one concludes that G′ was not maximal, a contradiction. Here one has to perform more case distinctions than in the proof of Lemma 4.7. See Fig. 6(b) for an illustration. □
Proof of Theorem 2.5
(i) ⇒ (ii). If the sum of S and T is direct and S⊕T is lattice-convex, then Lemma 4.5 implies that G S is connected. Furthermore, in view Lemma 4.6, we can apply Lemma 4.7 or Lemma 4.8, depending on whether k−ℓ>1 or k−ℓ=1. This yields the necessity.
(ii) ⇒ (i). Assume that S is given as in (ii). Since \(S \subseteq\mathbb{L}\) and the sum of \(\mathbb{L}\) and T is direct (see Lemma 4.1), we see that the sum of S and T is direct. It remains to prove that S⊕T is lattice-convex. For this we show that \(\mathbb{Z}^{2} \cap\overline{S\oplus T} \subseteq S \oplus T\). Taking into account \(\overline {S+T}=\overline{S}+\overline{T}\) and using Lemma 4.1 we obtain
Case 1: k−ℓ=1. We represent \(\overline{S}\) by
where α 1,α 2,β 1,β 2,γ 1,γ 2∈ℤ.
If we represent \(p \in\overline{DT}\) as p=x p w 1+y p w 2, then it can be verified directly that
In view of (16), each element of \(\mathbb{Z}^{2}\cap\overline{S \oplus T}\) can be given as s+t=s′+t′ with \(s\in\mathbb{L}\), t∈T, \(s' \in\overline{S}\), and \(t' \in \overline {T}\). Let s=iw 1+jw 2 for i,j∈ℤ. We show that s∈S arguing by contradiction. Assume that s∉S. Then, since \(S = \mathbb{L}\cap\overline{S}\), one of the inequalities on the right hand side of (17) is violated for x=i and y=j. But then, in view of (18), the inequality which is violated for x=i and y=j is also violated for x=i+x p and j+y p with p=t−t′. This contradicts the fact that \(s + t-t' = s'\in\overline{S}\). It follows that s∈S. Consequently, \(\mathbb{Z}^{2} \cap\overline{S \oplus T} \subseteq S + T\).
Case 2: k−ℓ>1. The proof of this case follows along the same lines. For x p ,y p as for the previous case we have −1<x p ,y p <1; the remaining arguments are analogous. □
Notes
We decided to use the notation \(\overline{K}\) instead of the standard notation \(\operatorname{conv}K\), because it allows us to write formulas in a more compact way.
Recall that, for x,y∈ℝ2, one has det(x,y)>0 if and only if (x,y) is a positively oriented basis of ℝ2.
References
Averkov, G., Bianchi, G.: Retrieving convex bodies from restricted covariogram functions. Adv. Appl. Probab. 39(3), 613–629 (2007)
Averkov, G., Bianchi, G.: Retrieving convex bodies from restricted covariogram functions, 21 pp., preprint, available at http://arxiv.org/abs/math/0702892, 2007
Averkov, G., Bianchi, G.: Confirmation of Matheron’s conjecture on the covariogram of a planar convex body. J. Eur. Math. Soc. 11(6), 1187–1202 (2009)
Averkov, G.: Detecting and reconstructing centrally symmetric sets from the autocorrelation: two discrete cases. Appl. Math. Lett. 22(9), 1476–1478 (2009)
Benassi, A., Bianchi, G., D’Ercole, G.: Covariogram of non-convex sets. Mathematika 56, 267–284 (2010)
Baake, M., Grimm, U.: Homometric model sets and window covariograms. Z. Kristallogr. 222, 54–58 (2007)
Bianchi, G., Gardner, R.J., Kiderlen, M.: Phase retrieval for characteristic functions of convex bodies and reconstruction from covariograms. J. Am. Math. Soc. 24(2), 293–343 (2011)
Bianchi, G.: Determining convex polygons from their covariograms. Adv. Appl. Probab. 34(2), 261–266 (2002)
Bianchi, G.: Matheron’s conjecture for the covariogram problem. J. Lond. Math. Soc. 71(1), 203–220 (2005)
Bianchi, G.: The covariogram determines three-dimensional convex polytopes. Adv. Math. 220(6), 1771–1808 (2009)
Bianchi, G.: The cross covariogram of a pair of polygons determines both polygons, with a few exceptions. Adv. Appl. Math. 42(4), 519–544 (2009)
Bianchi, G.: Geometric tomography of convex cones. Discrete Comput. Geom. 41(1), 61–76 (2009)
Bianchi, G., Segala, F., Volčič, A.: The solution of the covariogram problem for plane \(\mathcal {C}^{2}_{+}\) convex bodies. J. Differ. Geom. 60(2), 177–198 (2002)
Daurat, A., Gérard, Y., Nivat, M.: Some necessary clarifications about the chords’ problem and the partial digest problem. Theor. Comput. Sci. 347(1–2), 432–436 (2005)
Diestel, R.: Graph Theory, 3rd edn. Graduate Texts in Mathematics, vol. 173. Springer, Berlin, (2005)
Gardner, R.J.: Geometric Tomography, 2nd edn. Encyclopedia of Mathematics and its Applications, vol. 58. Cambridge University Press, Cambridge (2006)
Gardner, R.J., Gritzmann, P.: Discrete tomography: determination of finite sets by X-rays. Trans. Am. Math. Soc. 349(6), 2271–2295 (1997)
Gardner, R.J., Gronchi, P., Zong, Ch.: Sums, projections, and sections of lattice sets, and the discrete covariogram. Discrete Comput. Geom. 34(3), 391–409 (2005)
Herman, G.T., Kuba, A. (eds.): Discrete Tomography. Applied and Numerical Harmonic Analysis. Foundations, Algorithms, and Applications. Birkhäuser, Boston (1999)
Herman, G.T., Kuba, A. (eds.): Advances in Discrete Tomography and its Applications. Applied and Numerical Harmonic Analysis. Birkhäuser, Boston (2007)
Janot, C.: Quasicrystals: A Primer. Oxford University Press, Oxford (1997)
Klibanov, M.V., Sacks, P.E., Tikhonravov, A.V.: The phase retrieval problem. Inverse Probl. 11(1), 1–28 (1995)
Lang, S.: Algebra, 3rd edn. Graduate Texts in Mathematics, vol. 211. Springer, New York (2002)
Lemke, P., Skiena, S.S., Smith, W.D.: Reconstructing Sets from Interpoint Distances. Discrete and Computational Geometry. Algorithms Combin., vol. 25, pp. 507–631. Springer, Berlin (2003)
Moody, R.V.: Model sets: A survey, 28 pp., preprint, available at http://arxiv.org/pdf/math/0002020 (2000)
Nagel, W.: Orientation-dependent chord length distributions characterize convex polygons. J. Appl. Probab. 30(3), 730–736 (1993)
Rosenblatt, J., Seymour, P.D.: The structure of homometric sets. SIAM J. Algebr. Discrete Methods 3(3), 343–350 (1982)
Schneider, R.: Convex Bodies: The Brunn-Minkowski Theory. Encyclopedia of Mathematics and Its Applications, vol. 44. Cambridge University Press, Cambridge (1993)
Acknowledgements
We thank the anonymous referee for valuable suggestions regarding a previous version of this paper.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Averkov, G., Langfeld, B. On the Reconstruction of Planar Lattice-Convex Sets from the Covariogram. Discrete Comput Geom 48, 216–238 (2012). https://doi.org/10.1007/s00454-012-9416-6
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00454-012-9416-6