MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2016/7351861 7351861 Research Article The Projections of Convex Lattice Sets of Points in E2 Gu Yu 1 http://orcid.org/0000-0003-0003-0763 Si Lin 1 Soldovieri Francesco College of Science Beijing Forestry University Beijing 100083 China bjfu.edu.cn 2016 15112016 2016 10 05 2016 17 10 2016 25 10 2016 2016 Copyright © 2016 Yu Gu and Lin Si. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Can one determine a centrally symmetric lattice polygon by its projections? In 2005, Gardner et al. proposed the above discrete version of Aleksandrov’s projection theorem. In this paper, we define a coordinate matrix for a centrally symmetric convex lattice set and suggest an algorithm to study this problem.

Young Talent Program of Beijing YETP0770 State Forestry Administration 948 Project 2013-4-66 Fundamental Research Funds for the Central Universities 2015ZCQ-LY-01 National Natural Science Foundation of China 11371224 11001014
1. Introduction

According to Gardner , geometric tomography deals with the retrieval of information about a geometric object from data concerning its projections (shadows) on planes or cross sections by planes. When the geometric object is a convex body, there are many results from convex geometry. A famous example is the celebrated Aleksandrov’s projection theorem, that is, for an origin symmetric convex body (compact, convex subset with nonempty interior in the n-dimensional Euclidean space En); it can be determined by its projections. For more details, see . When the object is a convex lattice set, it models the atoms in a crystal which will lead to a range of topics concerned with electron microscopy in discrete tomography [2, 3].

If A is a set in En, we denote the convex hull of A by convA. The dimension of A is the dimension of its affine hull and is denoted by dimA. The notation for the usual orthogonal projection of A on a subspace H is AH. A convex lattice set F is a finite subset of the n-dimensional integer lattice Zn such that F=convFZn and we denote the cardinality of F by {F}. Let C(Zn) denote the set of all centrally symmetric convex lattice sets about origin of Zn.

In 2005, Gardner et al.  proposed the following discrete version of Aleksandrov’s projection theorem.

Question 1.

Let C1 and C2 belong to C(Zn) with dimC1=dimC2. If {C1H}={C2H} holds for every hyperplane H with a norm uZn{o}, do C1 and C2 coincide up to a translation?

In , Gardner et al. discovered a counterexample to this problem in Z2 with cardinality 11 (Figure 1). In , Zhou proved that the counterexample discovered by Gardner et al. was the only counterexample in Z2, up to unimodular transformations and with cardinality not larger than 17. In , Xiong proved that an origin symmetric convex lattice set can be uniquely determined by its lattice projection counts if y-coordinate’s absolute value is not bigger than 2 for every point in this origin symmetric convex lattice set and its cardinality is not 11.

Counterexample to Question 1 in Z2 with cardinality 11.

All the results above  show that the discrete Aleksandrov’s projection theorem must have additional hypotheses. Motivated by the work of Gardner, Gronchi, Zong, Zhou, and Xiong, in the present paper, we introduced the coordinate matrix for convex lattice set and present some results about coordinate matrix. Considering that counting of projection of a convex lattice set in the discrete Aleksandrov’s projection theorem is important, we propose an algorithm based on the coordinate matrix and its properties.

2. Coordinate Matrix

Let P belong to C(Zn). If P has ai points on the line x=i(-mim), where am0 and am+1=am+2==0, then we denote P by P[a-m,,a-1,a0,a1,,am]. Clearly, since P is an origin symmetric convex lattice set, we have that ai=a-i(1im), {P}=i=-mmai.

Definition 1.

For P[a-m,,a-1,a0,a1,,am], denote the minimum y-coordinate by yi0 and the maximum y-coordinate by yi1 on the line x=i(-mim), respectively. Define the coordinate matrix of P[a-m,,a-1,a0,a1,,am] by the following: (1)x0x1xixmy01y11yi1ym1y00y10yi0ym0.

Obviously, a convex lattice set and its coordinate matrix can be uniquely determined by each other.

Example 2.

Considering the coordinate matrix of C1 in Figure 1, we get (2)012100-1-20.

We summarize some basic properties of coordinate matrix as the following propositions.

Theorem 3.

Let(3)PCZn,M=x0x1xixmy01y11yi1ym1y00y10yi0ym0be its coordinate matrix. Then M has the following properties:

{P}=2i=1myi1-yi0+y01-y00+2m+1.

{Py=0}=2m+1.

{Px=0}=2max1,,m{y01,,ym1;y00,,ym0}+1.

Proof.

Let(4)PCZn,M=x0x1xixmy01y11yi1ym1y00y10yi0ym0be its coordinate matrix.

By the definition of the coordinate matrix and the symmetry of P, we have (5)P=2i=1myi1-yi0+1+y01-y00+1=2i=1myi1-yi0+y01-y00+2m+1.

From the the definition of the coordinate matrix P, we can get it.

Note that P is symmetric about the origin.

For PC(Zn) and its coordinate matrix(6)M=x0x1xixmy01y11yi1ym1y00y10yi0ym0,we can denote by (±xi,±yir), (±xj,±yjs), and (±xk,±ykt) any three points in P, where 0i,j,km, and yi0yiryi1, yj0yjsyj1, yk0yktyk1.

Remark 4.

Here (±xi,±yir), (±xj,±yjs), and (±xk,±ykt) denote one of the following cases: (7)-1u1xi,-1u2yir,-1v1xj,-1v2yjs,-1w1xk,-1w2ykt,for u1,u2, v1,v2, w1, and w2{1,2}.

If(8)±xi±yir1±xj±yjs1±xk±ykt10,we know that (±xi,±yir), (±xj,±yjs), and (±xk,±ykt) can be regarded as three vertexes of a nondegenerate triangle.

Then, with 0λ1,λ2,λ31 and λ1+λ2+λ3=1, we have that (9)xp=λ1±xi+λ2±xj+λ3±xk,yp=λ1±yir+λ2±yjs+λ3±yktis a convex combination of (±xi,±yir), (±xj,±yjs), and (±xk,±ykt). If (xp,yp)Z2, by the convexity of P, (xp,yp) belongs to P. Thus we obtain the following theorem.

Theorem 5.

Let(10)M=x0x1xixmy01y11yi1ym1y00y10yi0ym0with xi=i and yi1>yi0, yi1,yi0Z,i=0,1,2,,m. Then M is the coordinate matrix of some PC(Zn) if and only if for any λ1,λ2,λ3 with 0λ1,λ2,λ31 and λ1+λ2+λ3=1 there exist xq,yq1,yq0 in M such that (11)xq=±xp,yq1±ypyq0.

3. Algorithm for Counting of Projection of a Convex Lattice Set

Let PC(Zn). The primitive directions of P, DS(P) is defined as (12)DSP=u=a,ba,bZ,  gcda,b=1,  p,qP,  such  that  lpqu.Here, gcd represents the greatest common divisor.

Let(13)M=x0x1xixmy01y11yi1ym1y00y10yi0ym0be the coordinate matrix of P, and u=(a,b)DS(P). Based on the line segment intersection algorithm , the computation of planar convex hulls algorithm , and the coordinate matrix of P, an algorithm named as plane sweep algorithm about {Pu},uDS(P) is given in Algorithm 1.

<bold>Algorithm 1: </bold>Plane sweep algorithm.

Require:

coordinate matrix M and primitive direction u.

Ensure:

{Pu}.

( 1 )    M¬. Let M correspond to points pir=(xi,yir),xiZ,yirZ.

( 2 )    Initialize an empty event queue D. Next, insert the point pir into D. (Rename pir to

pk [Remark 7] by the insert-order [Remark 6] of lattice points).

( 3 )    Initialize an empty stack S and an empty linked list T.

( 4 )    while (D is not empty) do

( 5 )     Determine the next event point pk in D and delete it.

( 6 )      HANDLEEVENTPOINT

( 7 )      ORIGINSYMMETRY

( 8 )    end while

( 9 )    Initialize an empty variable v [Remark 9].

( 10 ) if (xi=0,xiT) then

( 11 )   v=v+1.

( 12 ) else

( 13 )   v=v+0.

( 14 ) end if

( 15 )  Compute the value of (s+2t-w-v) [Remark 8], and assigned to {Pu}.

Remark 6.

Moving a horizontal sweep line l downwards over the plane, the sweep line halts at certain event points. If sweep line reaches an event point, then a lattice point must be added to event queue D. This point is tested against the ones already have been reached. If two points have the same y-coordinate, then the one with smaller x-coordinate will be added to event queue D first.

Remark 7.

Based on Theorem 3, M is a 3×(m+1) matrix and it presents i=0m(yi1-yi0+1) points. We denote i=0m(yi1-yi0+1) by N. Thus D could be recorded as {pk,1kN}.

Remark 8.

s , t, and w are obtained after the While-loop is completed. At the end of Algorithm 2, s and t represent the number of points in S and T, respectively. w is defined in Algorithm 3.

<bold>Algorithm 2: </bold>HANDLEEVENTPOINT.

( 1 )    Insert pk into T.

( 2 )    if (T) then

( 3 )     Keep pk in T.

( 4 )     if (S) then

( 5 )      Keep pk in T.

( 6 )     else

( 7 )      if (lkju) then

( 8 )       Delete pk from T. (Suppose pjS, 1jk and jZ. The projection of pk and

pj in the primitive direction u is the same point.)

( 9 )      else

( 10 )       Keep pk in T.

( 11 )      end if

( 12 )     end if

( 13 )   else

( 14 )    if (lkiu) then

( 15 )   Insert pi into S, and delete pk from T at the same time. (Suppose piT, 1ik

and iZ.)

( 16 )   else

( 17 )   if (S) then

( 18 )    Keep pk in T.

( 19 )   else

( 20 )    if (lkju) then

( 21 )     Delete pk from T. (Suppose pjS, 1jk and jZ)

( 22 )    else

( 23 )     Keep pk in T.

( 24 )    end if

( 25 )   end if

( 26 ) end if

( 27 ) end if

<bold>Algorithm 3: </bold>ORIGINSYMMETRY.

( 1 )   The origin symmetric point of pi except those on the y-coordinate is denoted by qi.

( 2 )    l q i p i should be used to express the line though piT and qi.

( 3 )   Initialize an empty variable w.

( 4 )   if (lqipiu) then

( 5 )    w=w+1.

( 6 )   else

( 7 )    w=w+0.

( 8 )   end if

Remark 9.

The points on the y-axis should be considered separately because their symmetric points about the origin are still on the y-axis. The If-loop at the end of Algorithm 1 is used to do this. v represents the number of points whose x-coordinates are zero in T.

The procedures for finding the lattice points of the projection overlapped as follows: we simply test if the direction of a straight line that is determined by two lattice points is parallel to the primitive directions. We could delete the lattice points which conform to the conditions in order to simplify the computation. Hence, the procedure HANDLEEVENTPOINT is defined as in Algorithm 2.

By the origin symmetry, we have less calculation of nearly half of the points. The next procedure ORIGINSYMMETRY, defined as in Algorithm 3, is used to handle the missing points which are caused by the origin symmetry. In this procedure the way for process the lattice points whose projection will be overlap is the same as the procedure HANDLEEVENTPOINT.

Remark 10.

If we do not consider the origin symmetry, we should calculate 2N-1-y01+y00 points. This will take O((2N-1-y01+y00)×(2N-2-y01+y00)) time.

The time required for the parts of HANDLEEVENTPOINT and ORIGINSYMMETRY dominates the total running time. In HANDLEEVENTPOINT, considering the worst case, we check every point with other N-1 points whether they are in the same straight line with primitive direction. This will take O((1/2)N×(N-1)) time. In ORIGINSYMMETRY, we check N points with other N-(y01-y00+1) points. This will take O((1/2)N×(N-1-y01+y00)) time.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the Young Talent Program of Beijing (YETP0770), the State Forestry Administration 948 Project (no. 2013-4-66), the Fundamental Research Funds for the Central Universities (no. 2015ZCQ-LY-01), and the National Natural Science Foundation of China (nos. 11371224 and 11001014).

Gardner R. J. Geometric Tomography 1995 58 New York, NY, USA Cambridge University Press Encyclopedia of Mathematics and Its Applications MR1356221 Gardner R. J. Gritzmann P. Discrete tomography: determination of finite sets by X-rays Transactions of the American Mathematical Society 1997 349 6 2271 2295 10.1090/s0002-9947-97-01741-8 MR1376547 2-s2.0-21744459240 Herman G. T. Kuba A. Discrete Tomography: Foundations, Algorithms, and Applications 1999 Boston, Mass, USA Birkhäuser Gardner R. J. Gronchi P. Zong C. Sums, projections, and sections of lattice sets, and the discrete covariogram Discrete & Computational Geometry 2005 34 3 391 409 10.1007/s00454-005-1169-z MR2160045 2-s2.0-23944469960 Zhou J. On the projections of convex lattice sets Acta Mathematica Sinica 2010 26 10 1969 1980 10.1007/s10114-010-7164-1 MR2718094 2-s2.0-77956615220 Xiong H. On a discrete version of Alexandrov's projection theorem Acta Mathematica Sinica (English Series) 2013 29 8 1597 1606 10.1007/s10114-013-1354-6 MR3073577 2-s2.0-84880064468 de Berg M. Cheong O. van Kreveld M. Overmars M. Computational Geometry Algorithms and Applicationsy 2008 Berlin, Germany Springer 10.1007/978-3-540-77974-2 MR2723879