Next Article in Journal
Thermo-Elastodifusive Waves in Semiconductor Excitation Medium with Laser Pulses under Two Temperature Photo-Thermoelasticity Theory
Previous Article in Journal
A Novel Optical-Based Methodology for Improving Nonlinear Fourier Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Equivalence, Partial Order and Lattice of Neighborhood Sequences on the Triangular Grid

Department of Mathematics, Faculty of Arts and Sciences, Eastern Mediterranean University, North Cyprus, Mersin-10, Famagusta 99450, Turkey
Some parts of the paper were written during a research leave of absence at Faculty of Informatics, Eszterházy Károly Catholic University, 3300 Eger, Hungary.
Submission received: 24 October 2022 / Revised: 23 November 2022 / Accepted: 24 November 2022 / Published: 29 November 2022

Abstract

:
In (digital) grids, neighbor relation is a crucial concept; digital distances are based on paths through neighbor points. Digital distances are significant, e.g., in digital image processing for giving an approximation of the Euclidean distance and allowing incremental algorithms on images. Neighborhood sequences (i.e., infinite sequences of the possible types of neighbors) are defining digital distances with a lower rotational dependency than the distances based only on a sole neighborhood. They allow one to change the used neighborhood condition in every step along a path. They are defined in various grids, and they can be periodic. Generalized neighborhood sequences do not need to be periodic. In this paper, the triangular grid is studied. An equivalence and two partial order relations on the set of generalized and periodic neighborhood sequences are shown on this grid. The first partial order, the “faster” relation, is based on distances defined by neighborhood sequences, and it does not provide a lattice but gives a relatively complex relation for neighborhood sequences with a short period. The other partial order, the relation “componentwise dominate”, defines a complete distributive lattice on the set of generalized neighborhood sequences. Finally, a relation of the above-mentioned relations is established. Important differences regarding the cases of the square and triangular grids are also highlighted.

1. Introduction

The present paper has various roots in mathematics. On the one hand, digital geometry (a field that is closely connected to graph theory, geometry and combinatorics and that has applications in crystallography, image processing and computer graphics [1,2]), including the theory of digital (path-based) distances on various grids, was started by the seminal works of Rosenfeld and Pfaltz in the 1960s. They defined the two basic movements (steps) on the square grid (see [3,4]): the cityblock (also called the Manhattan taxicab) motion that allows horizontal and vertical steps only, while the chessboard motion allows diagonal steps also as when a king moves on the chessboard.
In digital spaces, e.g., on a computer screen, the geometry is not the same as in the continuous spaces: instead of the Euclidean geometry, digital geometry is applied. There is an underlying grid, and pixels or gridpoints are arranged accordingly; the neighborhood relation is crucial, while no neighbor points exist in Euclidean geometry. Consequently, digital distances are based on the steps/moves to neighbor pixels/points. However, the digital distances based on the cityblock or on the chessboard moves are a rough approximation of the Euclidean distance and have a large rotational dependency, i.e., to go in different directions, e.g., horizontal and diagonal (0 and 45 , respectively), the same number of steps leads to connected points with a Euclidean distance of a 2 factor.
Therefore, Rosenfeld and Pflatz previously coined the idea of alternating use of these two motions to have a digital distance with a lower rotational dependency [4]. These octagonal distances are obtained by the mixed use of the two basic motions. The theory was further developed when arbitrarily long periodic sequences of cityblock and chessboard steps, named neighborhood sequences, were studied in [5,6] (to distinguish those neighborhood sequences. Some more general ones were introduced later on, and we refer to these as periodic neighborhood sequences throughout this paper).
Digital distances are useful for stepwise approximations and simulations with incremental algorithms, e.g., computing distance transform in digital image processing and various algorithms in computer graphics. Moreover, as these distances always have integer values, it is easy to use them without rounding errors. Das et al. also provided a formula for calculating the distance between two arbitrary gridpoints based on any periodic neighborhood sequence. Based on this formula, a partial order relation on the set of periodic neighborhood sequences was investigated in [7].
This result connects our topic to other fields of mathematics. When there is a set of objects, e.g., the set of neighborhood sequences (and distance functions based on them), one may be interested to know how the objects of the set are related to each other. Various relations and algebraic concepts are known to help to describe the inner structure of a set. Now, we come back to the history of neighborhood sequences in the next section with a brief literature review.
In this paper, we use both periodic and generalized neighborhood sequences. Our motivation here is to show that similar algebraic concepts can be studied on sets of neighborhood sequences on a triangular grid and a square grid; thus, our results are somewhat analogous to the results of [7,8] with some interesting and remarkable differences. Moreover, a better understanding of the structure of the sets of neighborhood sequences allows us to know more about the distances based on them and also helps to find adequate neighborhood sequences and/or distances based on them in various applications.
On the triangular grid, there is also a natural equivalence relation on the sets of neighborhood sequences based on the distance function that they define; thus, after recalling the necessary definitions and preliminaries, we start our results with this equivalence relation. Although this result appears to be an easy consequence of some previous studies (see, e.g., [9]), we present it for the following two reasons.
On the one hand, it highlights an important difference between the square and triangular grids, which has various interesting consequences later on. On the other hand, the partial order that we consider first is a relation on the above-mentioned equivalence sets. This partial order and equivalence relation have a close connection to the distances defined by the neighborhood sequences. Further, another partial order relation is also studied. By this second partial order relation, lattice structures can also be obtained with some additional properties.
We also highlight the main differences between the cases of square and triangular grids in view of the mentioned relations and their analogous relations on the other grid, as the latter one is, in fact, not a point lattice (i.e., the triangular grid is not translation invariant to all grid vectors). Our last two theorems establish some connections among the considered relations. More precisely, our results are as follows. (For a better explanation, see Section 4).
The first two results are more or less the rephrasing or formal stating of known things:
  • On the set T of generalized neighborhood sequences of the triangular grid, the relation of defining the same distance function is an equivalence relation (Theorem 1).
  • The faster relation is a partial order both on the equivalence classes of T and of the set T of periodic neighborhood sequences on the triangular grid (Theorem 2).
The following results are our main results:
  • The set T does not form a lattice under the faster relation (Theorem 3).
  • The set T does not form a lattice under the faster relation (Theorem 4).
  • Incomparability of two generalized neighborhood sequences of the square grid under the componentwise dominate relation implies that the given neighborhood sequences cannot have the same speed, i.e., it cannot happen that both of them are faster than the other (Proposition 5).
  • The set S of generalized neighborhood sequences of the square grid is a Boolean lattice under the componentwise dominate relation (Theorem 5).
  • The set T is a complete distributive lattice under the componentwise dominate relation (Theorem 6); however, it is not a Boolean lattice (Proposition 7).
  • The set T is a distributive lattice under the componentwise dominate relation (Theorem 7); however, it may not be complete (Theorem 10).
  • The maximal equivalent neighborhood sequence of any generalized neighborhood sequence of T is defined and computed (in the proof of Theorem 8).
  • The relations faster and componentwise dominate have more complex and sophisticated structures on T than do their analogous relations on S (Theorems 11 and 12).
In the next section, a brief literature review is given allowing the reader to better understand and place our results. Then, in Section 3, necessary definitions and preliminaries are recalled about the triangular grid itself, about the neighborhood sequences on the triangular grid and about lattice theory. These parts before starting with the real results in Section 4 may seem lengthy; however, we would like to have the paper self-consistent such that readers without much background in these fields can also read and understand how the results are connected to these fields. Section 5 shows how our results can be applied, and our concluding section summarizes our results and also give links to future works.

2. Literature Review

Here, we recall parts of the history of digital distances. As we previously mentioned, Rosenfeld and Pfaltz defined the first digital distances. Around the same time when Das and their co-authors conducted their fundamental studies regarding periodic neighborhood sequences in the 1980s, a similar type of generalization of the basic digital distances was introduced under the same “neighborhood sequence” name [10,11] using some more general predefined motions (but still restricting the work to periodic sequences).
In [8,12], the concept of neighborhood sequences appeared with another type of generalization that also allowed non-periodic infinite neighborhood sequences. In [8], partial order relations were introduced, and the lattice properties of the set of these neighborhood sequences were analyzed extending the results of [7]. In [13,14], a greedy shortest path algorithm, a necessary and sufficient condition for neighborhood sequences to define metrical distances and a formula to compute distances between any pair of points was proven.
Higher dimensional extensions of these digital distances were also studied, see, e.g., [15]. Approximation of the Euclidean distance has been studied in [6,16,17,18] in various dimensions. Various image processing algorithms, e.g., distance transforms [19,20,21], have been developed; moreover, the distances based on neighborhood sequences have a close connection to integer sequences [22].
A generalization giving various real-valued weights to the diagonal steps was also investigated to give perfect approximation of the Euclidean distance on a perimeter of a square [23,24], while another type of generalization using a larger set of neighbor relations was developed under the name “broadcasting sequences” [25,26,27]. These papers show that neighborhood sequences have already been efficiently applied on a square grid.
Another way to have a digital distance with a lower rotational dependency is to use a non-traditional grid instead of the traditional (i.e., the square/rectangular) grid. The hexagonal and the triangular grids also provide regular tiling of the plane, and moreover they have better symmetric properties compared with the square grid [28,29,30].
In [31], Deutsch described the basic three types of neighborhoods in a triangular grid. Stojmenovic and Nagy, in [32] and in [9,33], respectively, used three coordinate values to represent a honeycomb grid and triangular tiling of the plane. These are dual notations of the same structure. In this paper, we use a triangular grid with triangle pixels (trixels for short). For defining digital distances on this grid, the concept of a generalized neighborhood sequence was used in [9,33]. The applicability of the neighborhood sequences and of the distances based on them is shown in [34,35], where applications in image processing and in networking are presented.

3. Basic Definitions and Preliminaries

To allow us to formally write our results, in this section, we provide the basic definitions and notations. In this paper, N , N 0 and Z denote the set of positive integers, nonnegative integers and integers, respectively. We also recall some basic facts about the neighborhood sequences and the distances based on them.
First, we recall the description of the triangular grid using three coordinate values to address the cells of the triangular grid, i.e., trixels.

3.1. Description of the Triangular Grid

Figure 1 (right) shows a part of the grid with coordinate axes and coordinate triplets for the points of the grid (for details, we refer to [9,33,36]). The coordinate system reflects the symmetry of the grid. The points are represented by trixels; however, throughout the paper, we usually refer to them as points.
If the sum of the coordinates of a point p is 0, then we call this point even. If the sum is 1, then we call this point odd. Graphically, the even points (as triangles) have the same orientation as the origin, and the odd points (triangles) have the opposite orientation.
Since we have two types of points, the triangular grid is not a point lattice, e.g., not every grid vector transforms the grid to itself.
Definition 1.
Let p and q be two points in the triangular grid (plane). The i-th coordinate of the point p is indicated by p ( i ) . Let m be 1, 2 or 3. The points p and q are m-neighbors, if the following two conditions hold:
  • | p ( i ) q ( i ) | 1 for 1 i 3 ,
  • i = 1 3 | p ( i ) q ( i ) | m .
In the case of equality at the second condition, we use the term strict m-neighbors as well.
The coordinate-values of the points and the definition of neighborhood relation harmonize (see Figure 1) and return the same neighborhood relations as the ones used in [31].
On the square grid, the description of neighborhood relations is similar to Definition 1 using the Cartesian coordinate system and addressing the points of grid by the elements of Z 2 . Accordingly, the one-neighbor relation corresponds to the cityblock neighborhood, while the two-neighbor relation corresponds to the chessboard neighborhood.
Now, we recall some known parts of the theory of neighborhood sequences concentrating on the triangular grid.

3.2. Neighborhood Sequences

We start this subsection with the definition of neighborhood sequences for both square (see, e.g., [5,8,12]) and triangular grids ([9,33,34]).
Definition 2.
A neighborhood sequence on the triangular grid (on the square grid) is an infinite sequence B = ( b ( i ) ) i = 1 , where b ( i ) { 1 , 2 , 3 } ( b ( i ) { 1 , 2 } , resp.) for all i N . If there exists l N such that i N , b ( i ) = b ( i + l ) , then B is periodic (with period l). In this case, the abbreviation B = ( b ( 1 ) , , b ( l ) ) is used.
Let T (S) and T ( S ) denote the sets of general and periodic neighborhood sequences on the triangular (square) grid, respectively.
In some papers, to distinguish the general scenario, when non-periodic neighborhood sequences can also be used, the term “generalized neighborhood sequence” is used [8,12]; however, here and on the triangular grid in general, if we do not restrict ourselves to only periodic ones, we may simply use the term “neighborhood sequence”. We use neighborhood sequences to define path-based (digital) distances as follows.
Definition 3.
Let p and q be two points (in the triangular grid or in the square grid), and let B = ( b ( i ) ) i = 1 be a neighborhood sequence (respectively). The point sequence Π ( p , q ; B ) of the form ( p = p 0 ) , p 1 , , ( p m = q ) , with b ( i ) -neighbor points p i 1 and p i for each i { 1 , , m } , is called a B-path, i.e., a path from p to q with respect the neighborhood sequence B. Further, the length | Π ( p , q ; B ) | of the path Π ( p , q ; B ) is m. Let a B-path from p to q with minimum length (with respect to B) be denoted by Π * ( p , q ; B ) . The distance from p to q is defined as the length of such a shortest path and is written as d ( p , q ; B ) = | Π * ( p , q ; B ) | . Since the distance is defined by B, it is called the B-distance.
The above-defined distance functions may not satisfy the metric properties. On the rectangular grid, there are already neighborhood sequences that define distances without satisfying the triangular inequality; however, these distances are always symmetric on that grid. Here, we show that the case of the triangular grid is more interesting and challenging.
Example 1.
Let B = ( 3 , 1 ) , p = ( 0 , 0 , 0 ) , q = ( 1 , 1 , 1 ) and r = ( 2 , 2 , 1 ) . In this case, d ( p , q ; B ) = 1 , d ( q , r ; B ) = 1 but d ( p , r ; B ) = 3 . Thus, this distance does not meet the triangular inequality. Furthermore, there is a property of B-distances of a triangular grid that does not occur in the rectangular case; it is possible that the B-distance is non-symmetric. Let s = ( 0 , 1 , 0 ) , t = ( 1 , 1 , 1 ) , then d ( s , t ; B ) = 3 opposite to d ( t , s ; B ) = 2 .
To use the neighborhood sequences in various applications, a/the shortest B-path and the B-distance are computed from any point to any other and necessary and sufficient conditions are proven for a neighborhood sequence to define a metric (see [12,13,14] for these results on the rectangular grids and [9,33,37] on the triangular grid).
The following fact about B-distances on the square grid can be shown, e.g., based on the algorithm providing a short/the shortest B-path. This is connected to one of the important differences between B-distances on the two grids.
Proposition 1.
On the square grid, each neighborhood sequence defines a unique distance function injectively, i.e., if B 1 , B 2 are two different neighborhood sequences ( B 1 B 2 ), then there are points p , q Z 2 such that d ( p , q ; B 1 ) d ( p , q ; B 2 ) .
This is not the case on the triangular grid; thus, we require some further notions. The next concept, the minimal equivalent neighborhood sequence, and the related result are recalled from [9].
Definition 4.
Let B = ( b ( i ) ) i = 1 and B m = ( b m ( i ) ) i = 1 be two neighborhood sequences. B m is the minimal equivalent neighborhood sequence of B if the following conditions hold:
  • d ( p , q ; B ) = d ( p , q ; B m ) for all points p , q ; and
  • for each neighborhood sequence B , if d ( p , q ; B ) = d ( p , q ; B ) for all points p , q , then b m ( i ) b ( i ) for all i N .
Lemma 1
(Nagy [9]). The minimal equivalent neighborhood sequence B m of B is uniquely determined and is given by
  • b m ( i ) = b ( i ) , if b ( i ) < 3 ;
  • b m ( i ) = 3 , if b ( i ) = 3 and if there is no j < i , such that b m ( j ) = 3 ;
  • b m ( i ) = 3 , if b ( i ) = 3 and there are some b m ( l ) = 3 with l < i and k = j + 1 i 1 b m ( k ) is odd, where j = max l l < i , b m ( l ) = 3 ;
  • b m ( i ) = 2 , otherwise.
As we shall see later, this concept and result has a good connection with the “componentwise dominate” relation (that we recall formally in Definition 9). By that terminology, one can say that any neighborhood sequence B componentwise dominates its minimal equivalent neighborhood sequence.
Now, we recall the main theorem of [37]. First, we need some further notions and notations.
For a neighborhood sequence B, let B m be the minimal equivalent neighborhood sequence of B, and let B m be the neighborhood sequence that is obtained from B m by replacing its first value 3 by 2 (if it has any value 3). Further, the sequence obtained from B by replacing all its values of 3 with values of 2 is denoted by B ( 2 ) = b 2 ( i ) i = 1 . Thus, B ( 2 ) = B m ( 2 ) = ( B m ) ( 2 ) = min { b m ( i ) , 2 } i = 1 .
For two points p = ( p ( 1 ) , p ( 2 ) , p ( 3 ) ) and q = ( q ( 1 ) , q ( 2 ) , q ( 3 ) ) , let x ( 1 ) = p ( 1 ) q ( 1 ) , x ( 2 ) = p ( 2 ) q ( 2 ) , x ( 3 ) = p ( 3 ) q ( 3 ) and w = ( w ( 1 ) , w ( 2 ) , w ( 3 ) ) be defined as the sorted difference vector of p and q: let w ( 1 ) be the value with highest absolute value among { x ( 1 ) , x ( 2 ) , x ( 3 ) } , let w ( 3 ) be the value with the smallest absolute value among these values, and let w ( 2 ) be the third value. (For disambiguation, if the largest absolute value appears twice, let the positive one be w ( 1 ) and the negative one be w ( 2 ) .)
Now, we are ready to recall the formula for computing the B-distance.
Lemma 2
(Nagy [37]). Let two points p = ( p ( 1 ) , p ( 2 ) , p ( 3 ) ) and q = ( q ( 1 ) , q ( 2 ) , q ( 3 ) ) and a neighborhood sequence B = ( b ( i ) ) i = 1 be given. Let the sequences B m , B m and B ( 2 ) be given accordingly. Further, let Z = i = 1 k 1 b ( i ) where k is the index of the first element 3 of B (if there is no value 3 in B, then Z = ).
Now, let
D ( B ) = max | w ( 1 ) | , max i | w ( 1 ) | + | w ( 2 ) | > j = 1 i 1 b ( 2 ) ( j ) , max i j = 1 3 | w ( j ) | > j = 1 i 1 b ( j ) .
The B-distance from p to q can be computed by the following formula
d ( p , q ; B ) = D ( B m ) , if p is even , Z is or even and w contains two negative values ; D ( B m ) , if p is odd , Z is or odd and w contains two negative values ; D ( B m ) , if p is odd , Z is or even and w contains two positive values ; D ( B m ) , if p is even , Z is or odd and w contains two positive values ; D ( B m ) , otherwise .
The concept of digital disks based on neighborhood sequences is also recalled (see, e.g., [34]).
Definition 5.
Let B be a neighborhood sequence on the triangular grid. Then, for any point p of the grid and for any natural number r N ,
Δ p ( B , r ) = { q | d ( p , q ; B ) r }
denotes the digital disk (based on the B-distance) with center p and radius r. In the case where p is the origin, the notation can be simplified to Δ ( B , r ) .
We note here that the digital disk with a radius r depends only on the first r elements of B.

3.3. Lattice Theory

We assume that the reader is familiar with the basic concepts of lattice theory, otherwise they aer referred to review the following concepts in [38,39] or another standard textbook on the topic. Here, we list them to show the notation we use:
  • Equivalence relation ( P , ) with equivalence classes x ¯ = { y | x y } .
  • Partially ordered set ( P , ) .
  • Supremum (least upper bound, x y for any x , y P ) of a set.
  • Infimum (greatest lower bound x y ) of a set.
  • Lattice, distributive lattice.
  • Complete lattice (with a supremum X and an infimum X for any X P ).
  • Complete distributive lattice and Boolean lattice (Boolean algebra).

4. Relations on the Set of Neighborhood Sequences in the Triangular Grid

The results of the paper are shown in this section. First, the equivalence relation and the partial order relation “faster” on the equivalence classes of neighborhood sequences are shown, and then some algebraic properties of this partial order relation are presented. Finally, the “componentwise dominate” relation will be investigated with details.

4.1. An Equivalence Relation

In this section, we formally establish the equivalence relation among the neighborhood sequences that has been previously mentioned.
Definition 6.
Let B 1 and B 2 be two neighborhood sequences in T. We say that they are equivalent, i.e., B 1 B 2 , if, for every pair of points p , q , d ( p , q ; B 1 ) = d ( p , q ; B 2 ).
Based on the definition of minimal equivalent neighborhood sequences, we may establish the following statement.
Lemma 3.
For any two neighborhood sequences B 1 and B 2 , B 1 B 2 , if and only if their minimal equivalent sequences coincide.
Based on that, we can conclude the following result.
Theorem 1.
The relation ∼ is an equivalence relation on the set T of generalized neighborhood sequences.
Proof. 
Using Definition 6, ∼ is reflexive, symmetric and transitive, and therefore it is an equivalence relation. □
By comparing this result to the square and nD rectangular grids, it is important to note that there are no two different neighborhood sequences that generate the same distance function on those grids. By the mirror of this fact, the forthcoming relations on the neighborhood sequences on the triangular grid also have some interesting properties in which the case of the triangular grid differs from the case of rectangular grids.
For further purposes, we fix the following notation.
Definition 7.
Let B ¯ denote the equivalence class of neighborhood sequences defined by ∼:
B ¯ = { B | B B } .
We say that B is a representative of the set B ¯ or that B represents this class.
There exists neighborhood sequence B for which the class B ¯ contains only the neighborhood sequence B itself and no other neighborhood sequences. One can obtain neighborhood sequences with this property based on Lemma 1, e.g., every neighborhood sequence without containing any value of 3 is satisfying this property.

4.2. A Partial Order Defined by B-Distances

Since there are several, there are actually infinitely many distance functions that are defined by neighborhood sequences; thus, the question arises of how they are related to each other. The next definition is recalled from [9].
Definition 8.
For any B 1 , B 2 T , we define the relation “faster” denoted by * in the following way:
B 1 * B 2 d ( p , q ; B 1 ) d ( p , q ; B 2 )
for all points p , q .
The relation * can be understood on the set of equivalence classes of neighborhood sequences as well; however, for simplicity, we use the above notations, where B 1 and B 2 represent their equivalence sets, respectively. When we are talking about periodic neighborhood sequences, then only the elements of B ¯ T are counted.
Notice that equality is allowed by the relation * ; thus, the term “is faster than or has the same speed” would match better to it; however, for traditional reasons, we use the short term “faster”.
Similarly, the relation * is defined on the sets of S and S , and it was studied in [8,13] as it also plays an important role in the metricity condition.
Proposition 2
(Fazekas et. al. [8]). On the square grid S , * forms a complete distributive lattice. On the other hand, S , * is not a lattice.
The equivalence classes for the neighborhood sequences on the square grid are singleton sets, i.e., each of them has exactly one element, and therefore this element can be used to represent this class or used instead of the class itself.
The next examples demonstrate the non-triviality of our partial order for the case of periodic neighborhood sequences.
Example 2.
Let B 1 , B 2 T such that B 1 = ( 1 , 2 ) , B 2 = ( 1 , 1 , 2 , 2 , 2 ) . Further, let the points o = ( 0 , 0 , 0 ) , p = ( 2 , 1 , 0 ) , q = ( 1 , 0 , 0 ) and r = ( 4 , 4 , 0 ) . Then, d ( o , p ; B 1 ) = 2 < 3 = d ( o , p ; B 2 ) , d ( o , q ; B 1 ) = 1 = d ( o , q ; B 2 ) but d ( o , r ; B 1 ) = 6 > 5 = d ( o , r ; B 2 ) . Therefore, the functions B 1 -distance and B 2 -distance are not comparable, in general.
We should note here that similar examples can also be constructed on the square and other grids using B-distances. Our next example, however, does not work on point lattices but does work on the triangular grid.
Example 3.
Let B 1 = ( 1 , 3 ) , B 2 = ( 3 , 1 ) and the points p = ( 0 , 0 , 0 ) , q = ( 1 , 1 , 2 ) and r = ( 1 , 1 , 2 ) . Then, one may see that d ( p , q ; B 1 ) = 3 > 2 = d ( p , q ; B 2 ) but d ( p , r ; B 1 ) = 2 < 3 = d ( p , r ; B 2 ) .
This example shows that, although we have the same sum of the values in the sequences, (i.e., b 1 ( 1 ) + b 1 ( 2 ) = b 2 ( 1 ) + b 2 ( 2 ) = 4 ), their order makes a difference on the distance: the steps of a shortest path may not commute [33]. We note here that, on the square- and nD grids the B-distance of the point-pair p , q depends only on the multiset of used elements of the neighborhood sequence and does not depend on the order of the first d ( p , q ; B ) 1 elements of the neighborhood sequence [13,14].
Now, after these introductory examples and thoughts, let us see how and which neighborhood sequences are comparable on the triangular grid.
Theorem 2.
The relation * is a partial order on the sets of equivalence classes of T and T .
Proof. 
By definition, the relation * is reflexive and transitive. In defining it for the equivalence sets, it is clear that, if B 1 * B 2 and B 2 * B 1 , then d ( p , q ; B 1 ) = d ( p , q ; B 2 ) , and thus B 1 ¯ = B 2 ¯ ; therefore, they represent the same class. □
As we have previously seen in Example 3, the neighborhood sequences ( 1 , 3 ) and ( 3 , 1 ) are not comparable under the relation * . In terms of digital disks, this means that Δ ( ( 1 , 3 ) , 2 ) Δ ( ( 3 , 1 ) , 2 ) and Δ ( ( 1 , 3 ) , 2 ) Δ ( ( 3 , 1 ) , 2 ) ; see also Figure 2, where digital disks with radii up to 3 are presented with the (periodic) neighborhood sequences ( 3 , 1 ) , ( 2 ) and ( 1 , 3 ) .
Based on Lemma 2, and by observing the digital disks obtained by neighborhood sequences (see [34,35,40] for details), we may consider the following important observations.
Proposition 3.
In the triangular grid, from any point
(a) 
By a 2-step, exactly the same region (i.e., set of points) is reached as by two consecutive 1-steps.
(b) 
By two consecutive 2-steps, a larger region can be reached than by using a 1-step and a 3-step in any order.
We can transform part (b) of the above proposition into the next lemma about the partial order * .
Lemma 4.
The following statements hold for the specific periodic neighborhood sequences mentioned below.
( 2 ) * ( 1 , 3 ) . ( 1 , 2 , 2 , 2 ) * ( 1 , 1 , 2 , 3 , 1 ) .
However, ( 2 ) * ( 3 , 1 ) and ( 3 , 1 ) * ( 2 ) .
Proof. 
The first and the second statements follow from the formula of Lemma 2.
The first part of the last statement of the lemma is proven, for instance, by looking the points that can be reached in one step: in one step, a larger set of points is reached by ( 3 , 1 ) than by ( 2 ) , i.e., Δ ( ( 2 ) , 1 ) Δ ( ( 3 , 1 ) , 1 ) . To prove the other part, consider the B-distances of the points p = ( 0 , 0 , 0 ) and q = ( 2 , 1 , 1 ) . By the formula of Lemma 2, d ( p , q ; ( 2 ) ) = 2 and d ( p , q ; ( 3 , 1 ) ) = 3 , proving the incomparability of the neighborhood sequences ( 3 , 1 ) and ( 2 ) under * on the triangular grid. □
Now, we establish our results on the structure of T and T under * . Some of the results may seem simple or easy to prove; however, they give important information about the sets T and T of neighborhood sequences. They also underline the differences between the distances defined on the square and on the triangular grids. Figure 3 shows the partial order * on the sets of periodic neighborhood sequences with periods 1, 2 and 3 for both of the square and triangular grids. As one can observe, the diagram of T with period 2 already appears more complicated than the diagram of S with period 3 (the diagrams on the top are recalled from [7]). The Hasse diagrams shown here compress a number of comparisons of B-distances.
One can observe on Figure 3, bottom-right, that, for instance, the neighborhood sequence ( 2 ) = ( 2 , 2 , 2 ) is faster than ( 1 , 3 , 2 ) , ( 1 , 3 , 3 ) , ( 2 , 1 , 3 ) , ( 2 , 2 , 1 ) and all the neighborhood sequences below them, i.e., ( 1 , 3 , 1 ) , ( 1 , 2 , 3 ) , ( 1 , 2 , 2 ) , ( 2 , 1 , 2 ) , ( 2 , 1 , 1 ) , ( 1 , 2 , 1 ) , ( 1 , 1 , 3 ) , ( 1 , 1 , 2 ) and ( 1 ) = ( 1 , 1 , 1 ) . By, e.g., the formula of Lemma 2, one can see that there is no point that can be reached in one step by any of the listed neighborhood sequences B but not with the neighborhood sequence ( 2 ) . That could be the base of the induction. The inheritance can be proven for each particular case (applying also, e.g., Proposition 3) by checking the sums of the prefix of the (minimal equivalent) neighborhood sequence B and comparing the same length prefix of ( 2 ) , i.e., to the number that is twice the length of the prefix.
On the other hand, ( 3 , 3 , 3 ) , ( 3 , 3 , 2 ) , ( 3 , 2 , 3 ) , ( 3 , 2 , 2 ) , ( 2 , 3 , 3 ) , ( 2 , 3 , 2 ) and ( 2 , 2 , 3 ) are faster than ( 2 ) . All those relations can be seen since any path defined by ( 2 ) is also a path with any of the listed neighborhood sequences; thus, the distance determined by ( 2 ) cannot be smaller than the B-distance defined by another neighborhood sequence in this list. (Trivially, ( 2 ) is faster than itself. We mention this only for the sake of completeness).
Further, there are neighborhood sequences such that they are incomparable with ( 2 ) under the relation faster. Each of those neighborhood sequences B are listed below with two points such that the B-distance of exactly one of the points, p, is smaller from the origin than its ( 2 ) -distance, while, for the second point, q, its B-distance is larger than its ( 2 ) -distance from the origin.
  • B = ( 3 , 1 , 1 ) , p = ( 1 , 1 , 1 ) and q = ( 2 , 3 , 1 ) ;
  • B = ( 3 , 1 , 2 ) , p = ( 1 , 1 , 1 ) and q = ( 3 , 3 , 0 ) ;
  • B = ( 2 , 3 , 1 ) , p = ( 2 , 2 , 1 ) and q = ( 3 , 3 , 0 ) ;
  • B = ( 3 , 3 , 1 ) , p = ( 1 , 1 , 1 ) and q = ( 3 , 3 , 0 ) ;
  • B = ( 3 , 2 , 1 ) , same points as above;
  • B = ( 3 , 1 , 3 ) , p = ( 1 , 1 , 1 ) and q = ( 3 , 0 , 3 ) .
The relations among other neighborhood sequences can be shown in a similar manner.
Now, as the first of our main theoretical results about * , we show the following.
Theorem 3.
The set of periodic neighborhood sequences T , * is not a lattice.
Proof. 
Let B 1 = ( 2 , 1 , 1 ) and B 2 = ( 1 , 3 , 2 ) be periodic neighborhood sequences. Notice that B 1 ¯ = { B 1 } , while
B 2 ¯ = ( b y ( i ) ) i = 1 b y ( i ) = 1 , i f   i 1 ( mod 3 ) b y ( i ) = 3 , i f   i 2 ( mod 3 ) b y ( i ) 2 , i f   i 0 ( mod 3 ) .
The least upper bound of these sequences in T , * can be written as
B = ( 2 , 2 , 2 , 1 , 3 , 2 , 1 , 3 , 2 , 1 , 3 , 2 , ) , which belongs to the equivalence set
B ¯ = ( b ( i ) ) i = 1 b ( i ) = 2 , i f   i { 1 , 2 , 3 } b ( i ) = 1 , i f   i > 3 , i 1 ( mod 3 ) b ( i ) = 3 , i f   i > 3 , i 2 ( mod 3 ) b ( i ) 2 , i f   i > 3 , i 0 ( mod 3 ) .
Notice that B M = ( 2 , 2 , 2 , 1 , 3 , 3 , 1 , 3 , 3 , 1 , 3 , 3 , 1 , 3 , 3 , ) B ¯ contains only three values of 2. Clearly, there is no element of B ¯ that is in T . Now, let us suppose, on the contrary, that there is a neighborhood sequence B = B 1 B 2 in T and B = ( b ( 1 ) , , b ( l ) ) for some value l N . Then, B * B in T but B B ¯ ; therefore, there must be i N such that b ( i ) > b M ( i ) holds.
  • B may contain an element larger than 2 among its first three elements. (Clearly, it cannot contain values smaller than 2 at these positions without contradicting the fact that B * B .) However, in this case, for the periodic sequence B = ( 2 ) , we have B * B (and B * B ), but (as one can easily check) B * B 1 and B * B 2 . That leads to a contradiction.
  • In the other case, it means that there is a value i = 3 j + 1 (for a value j N ) such that b M ( i ) = 1 ; however, b ( i ) > 1 . Without loss of generality, suppose that, for all j < i : b ( j ) = b ( j ) holds. Then, defining b ( k ) = b ( k ) for 1 k i + 5 and B = ( b ( 1 ) , , b ( i + 5 ) ) as a periodic neighborhood sequence with period i + 5 , clearly B * B 1 and B * B 2 . However, in this case, B * B , which leads to a contradiction.
Therefore, B = B 1 B 2 cannot be in T . □
By the previous proof, we demonstrated the following fact.
Corollary 1.
Even if the supremum of two periodic neighborhood sequences exists, it may not be periodic.
By the following theorem, we show the case with the set of general neighborhood sequences.
Theorem 4.
The set of neighborhood sequences on the triangular grid T , * is not a lattice.
Proof. 
Let
B 1 = { 3 , 2 , 1 , 1 , 1 , 3 , 1 , 3 , 1 , 3 , 1 , 3 , 1 , } ,
B 2 = { 2 , 2 , 1 , 2 , 1 , 2 , 2 , 2 , 2 , 2 , 2 , 2 , 2 , } ,
B 3 = { 2 , 2 , 1 , 1 , 1 , 3 , 1 , 3 , 1 , 3 , 1 , 3 , 1 , }
and
B 4 = { 2 , 2 , 1 , 1 , 1 , 1 , 3 , 1 , 3 , 1 , 3 , 1 , 3 , } .
(The last two values are repeated for each of these sequences up to infinity). Observe that B 1 , B 2 , B 3 and B 4 form a singleton equivalence set under ∼. Clearly, B 1 * B 3 and B 2 * B 3 . Let B 5 T such that B 1 * B 5 and B 2 * B 5 . Then, one can easily infer that B 5 * B 3 . Further, the first two elements of B 5 can be at most 2. Let us assume that the first difference of B 3 = b 3 ( k ) k = 1 and B 5 = b 5 ( k ) k = 1 occurs at their ( i + 1 ) st position, i.e., b 3 ( k ) = b 5 ( k ) for all k i ; however, b 3 ( i + 1 ) b 5 ( i + 1 ) . Then, B 5 * B 3 would imply that b 5 ( i + 1 ) > b 3 ( i + 1 ) . It can only be if b 3 ( i + 1 ) = 1 and b 5 ( i + 1 ) 2 , which contradicts B 1 * B 5 .
Clearly, B 1 * B 4 , B 2 * B 4 but B 3 * B 4 (since the order of elements 1 and 3 matters). Therefore, B 1 and B 2 have no greatest lower bound, completing the proof that T is not a lattice. □
In the previous proof, we used non-periodic sequences. Thus, it seems that none of T and T form a lattice under the relation * . This also highlights an important difference between the square and the triangular grid. In the next subsection, another relation will be investigated on T and T .

4.3. Lattices

The results of the previous subsection recommend that, instead of the relation * , we may use another relation to find some structure in T and in T . Now, we adopt another partial order relation that was investigated in [8] for S and S and is called “componentwise dominates” in [5] for some subsets of S . We have also seen that it plays a natural role when we defined and used minimal equivalent neighborhood sequences in T.
Definition 9.
For any two neighborhood sequences B 1 = ( b 1 ( i ) ) i = 1 , B 2 = ( b 2 ( i ) ) i = 1 , the relation is defined in the following way:
B 1 B 2 b 1 ( i ) b 2 ( i ) , for every i N .
Let us have a more detailed view of the structure that this relation gives to the sets S and S on the square grid. It should be noted that * is a proper refinement of ⊒ in S and S ; more formally, we can state the following results (based, e.g., on the formula to compute B-distances on the square grid [14] or on the properties of a faster relation from [13]).
Proposition 4.
Let B 1 , B 2 S . If B 1 B 2 , then B 1 * B 2 . Furthermore, if B 1 B 2 and B 1 B 2 , then B 2 B 1 .
Proof. 
It is a corollary of the formula to compute B-distances on the square grid; see Proposition 1 and Definition 9. □
Moreover, we state the following.
Proposition 5.
Let B 1 , B 2 S . If B 1 B 2 and B 2 B 1 , then B 1 * B 2 and B 2 * B 1 cannot be at the same time, i.e., at most, one of them may hold.
Proof. 
Generally, for B 1 , B 2 S , there could be three cases, either both B 1 * B 2 and B 2 * B 1 hold, or maybe exactly one of them, or maybe none of them.
If the condition B 1 B 2 and B 2 B 1 holds:
  • We may have that B 1 and B 2 are incomparable under * . Let us consider, e.g., B 1 = ( 1 , 1 , 2 , 2 , 2 ) and B 2 = ( 1 , 2 , 1 , 1 , 1 ) . Then, considering the points ( 0 , 0 ) and ( 1 , 2 ) , their B 2 -distance is 2; however, their B 1 -distance is 3. On the other hand, by considering ( 0 , 0 ) and ( 4 , 2 ) , their B 1 -distance is 4, while their B 2 -distance is 5.
  • However, it may also happen that, e.g., B 1 * B 2 , and we check the case, e.g., with B 1 = ( 1 , 2 , 2 , 1 , 2 ) and B 2 = ( 1 , 1 , 1 , 2 , 1 ) . (By changing the role of B 1 and B 2 , the opposite relation can be fulfilled.)
Now, to complete the proof, let us assume that B 1 B 2 and B 2 B 1 for B 1 = ( b 1 ( i ) ) i = 1 , B 2 = ( b 2 ( i ) ) i = 1 S . It follows that B 1 B 2 . Thus, there is a position i where the first mismatch occurs in these two sequences, i.e., b 1 ( j ) = b 2 ( j ) for all j < i but b 1 ( i ) b 2 ( i ) . Then, considering the points ( 0 , 0 ) and ( i , k + 1 ) with k = { j | b 1 ( j ) = 2 , j < i } , i.e., k is the number of 2s among the first i 1 elements of B 1 (and also of B 2 since their i 1 long prefix sequence are identical). Then, one of the B 1 - and B 2 -distances of these points is i, the other is i + 1 , and hence it cannot be that both B 1 * B 2 and B 2 * B 1 hold. □
Moreover, the following results are known.
Proposition 6
(Fazekas et al. [8]). S , is a complete distributive lattice with the greatest lower bound S = ( 1 ) and least upper bound S = ( 2 ) .
S , is a distributive lattice with the greatest lower bound S = ( 1 ) and least upper bound S = ( 2 ) ; however, it is not a complete lattice.
Our first result in this subsection is a kind of strengthening of (the first part of) the previous result on the square grid.
Theorem 5.
( S , ) is a Boolean lattice (Boolean algebra).
Proof. 
A Boolean algebra, by definition, is a complemented distributive lattice. By Proposition 6, ( S , ) is a complete distributive lattice, and we need to show only that it is also complemented. Having exactly two types of neighborhood, the complementary relation can easily be defined on them: with 1 and 2 as complements of each other. This complementarity relation can also be extended to the neighborhood sequences: the complement of the neighborhood sequence B = ( b ( i ) ) i = 1 , B c = ( 3 b ( i ) ) i = 1 contains the opposite element in every of its positions. In this way, both the axioms of complements: B B c = ( 2 ) and B B c = ( 1 ) are satisfied for every neighborhood sequence B on the square grid. Thus, the statement is proven. □
Now, we turn to analyze the componentwise dominate relation on neighborhood sequences of the triangular grid. In the following part, the sets T and T are analyzed under ⊒.
Theorem 6.
T , is a complete distributive lattice with the greatest lower bound T = ( 1 ) and least upper bound T = ( 3 ) .
Proof. 
By Definition 9, ⊒ is a partial order on T. Moreover, for any B 1 = ( b 1 ( i ) ) i = 1 , B 2 = ( b 2 ( i ) ) i = 1 T , we can define B 1 B 2 as min { b 1 ( i ) , b 2 ( i ) } i = 1 and B 1 B 2 as max { b 1 ( i ) , b 2 ( i ) } i = 1 . Furthermore, for any set P of neighborhood sequences B 1 , B 2 , , their greatest lower bound
P = min { b j ( i ) | B j P , B j = b j ( k ) k = 1 } i = 1 and least upper bound
P = max { b j ( i ) | B j P , B j = b j ( k ) k = 1 } i = 1 is determined; thus, the completeness of the lattice is shown. This implies also that T = ( 1 ) and T = ( 3 ) . Finally, the distributivity of the lattice is straightforward. □
However, there is a remarkable difference between the T and S under ⊒.
Proposition 7.
T , is not a Boolean lattice with B 1 B 2 = min { b 1 ( i ) , b 2 ( i ) } i = 1 and B 1 B 2 = max { b 1 ( i ) , b 2 ( i ) } i = 1 where B 1 = ( b 1 ( i ) ) i = 1 and B 2 = ( b 2 ( i ) ) i = 1 .
Proof. 
In the proof of the previous theorem, we previously proved that T , is a distributive lattice with the given operators. We need to show that the complement operation cannot be introduced on T fulfilling the required axioms. Let B = ( 2 ) , if there is its complement B c T , then both B B c = ( 1 ) and B B c = ( 3 ) would be fulfilled. However, then B c = ( 1 ) should be to make B B c = ( 1 ) true; however, B c = ( 3 ) would be required to fulfill B B c = ( 3 ) . This contradiction proves the result. □
Now, we continue with a result on T .
Theorem 7.
T , is a distributive lattice with the greatest lower bound T = ( 1 ) and least upper bound T = ( 3 ) .
Proof. 
The proof is similar to the proof of the previous theorem except the completeness by noting that the constructions mentioned there provide periodic neighborhood sequences with period of the least common multiplier of the periods of the given neighborhood sequences. □
In the next part of the paper, we establish connections on the relations studied in this paper. Let us see how the ordering ⊒ and the sets T and T can be connected to the B-distances. First, we use the relation ⊒ on specific subsets of T.
Theorem 8.
Let B be any neighborhood sequence in T. The relation ⊒ defines a complete distributive lattice on B ¯ .
Proof. 
The reflexivity, the antisymmetry and the transitivity follows from the the definition of ⊒. Therefore, B ¯ , is a partially ordered set for any B T .
It is clear that, for every B 1 , B 2 B ¯ , B 1 B 2 and B 1 B 2 exist in T and can be defined as B 1 B 2 = ( min { b 1 ( i ) , b 2 ( i ) } ) i = 1 , and B 1 B 2 = ( max { b 1 ( i ) , b 2 ( i ) } ) i = 1 .
It can be checked (e.g., by using Lemma 1) that both B 1 B 2 and B 1 B 2 are in B ¯ . Thus, B ¯ , is a lattice.
Let X = { B k | B k B ¯ , k K } be a subset of B ¯ with some index set K. Let now, b ( i ) = max k K b k ( i ) . Clearly, B = ( b ( i ) ) i = 1 is the least upper bound of the subset X, and thus ( B ¯ , ) is a complete lattice.
The statement that this lattice is distributive is trivial. Moreover, B ¯ = B m , the minimal equivalent neighborhood sequence for B, and
B ¯ =
b M ( i ) = b ( i ) , if   b ( i ) = 1   or   b ( i ) = 3 , 2 , if   b ( i ) = 2   and   there   is   no   j < i , such   that   b M ( j ) = 3 , 3 , if   b ( i ) = 2   and   there   are   some   b M ( l ) = 3 with   l < i , and   | { k | b M ( k ) = 1 , j < k < i } | is   even , where   j = max l l < i , b M ( l ) = 3 , 2 , otherwise . i = 1
Thus, the proof of the theorem is completed. □
The neighborhood sequence defined by the last formula of the previous proof can be called the maximal equivalent neighborhood sequence  B M of B. Analogously to Lemma 1, it can be shown that it is uniquely defined for every B T . The case of periodic neighborhood sequences is more interesting due to the fact that there are periodic neighborhood sequences such that their minimal/maximal equivalent neighborhood sequences are not periodic and vice versa.
Example 4.
B = ( 3 ) is periodic. It is the same as its maximal equivalent neighborhood sequence. However, its minimal equivalent neighborhood sequence is not periodic, it is ( 3 , 2 , 2 , 2 , 2 , 2 , 2 , ) (having only 2’s after the first 3). The next non-periodic neighborhood sequence belongs to the same equivalence class: ( 3 , 2 , 2 , 3 , 2 , 2 , 2 , 2 , 3 , 2 , 2 , ) having 3 at every place i where i = j 2 for some j N .
Example 5.
Let B = ( 2 , 1 , 3 ) , and then its minimal equivalent neighborhood sequence is the same; however, its maximal equivalent neighborhood sequence,
B M = ( 2 , 1 , 3 , 3 , 1 , 3 , 3 , 1 , 3 , 3 , 1 , 3 , 3 , 1 , 3 , ) , is not periodic (having exactly one value 2).
Therefore, it is interesting to analyze what we can say about subsets of T containing only periodic neighborhood sequences.
Theorem 9.
Let B be any neighborhood sequence in T . The relation ⊒ defines a distributive lattice on B ¯ T that may not be complete.
Proof. 
Since ⊒ is a partial order on any set B ¯ (see (the proof of) Theorem 8), it is a partial order on the subset B ¯ T . Now, we show that, for any two representatives B 1 and B 2 of the class B ¯ T have B 1 B 2 and B 1 B 2 elements in B ¯ T . Let l 1 and l 2 be the periods of B 1 and B 2 , respectively. For simplicity, let l = l 1 · l 2 . (The least common multiplier of l 1 and l 2 suffices.) Then, let us construct B 1 B 2 with period l: ( min { b 1 ( i ) , b 2 ( i ) } ) i = 1 l . Similarly, the construction B 1 B 2 = ( max { b 1 ( i ) , b 2 ( i ) } ) i = 1 l .
Therefore, ( B ¯ T , ) form a lattice for any B T . By the constructions above, it is clear that this lattice is distributive.
Now, we show an example for ( B ¯ T , ) that proves that this lattice is not complete.
Let B @ = ( 3 ) . Then, B @ ¯ contains every neighborhood sequence B that have the following properties: b ( 1 ) = 3 , and b ( i ) 2 for every i N .
Now, a monotonously increasing and a monotonously decreasing sequence in B @ ¯ will be shown, such that they have a common “limit sequence” B l i m but B l i m is not in B @ ¯ T . In this way, we show that B @ ¯ , is not complete. For technical purposes, we introduce the following notation. Any periodic neighborhood sequence B = ( b ( 1 ) , , b ( l ) ) T can also be written in the form having period 2 l : ( b ( 1 ) , , b ( l ) , b ( 1 ) , , b ( l ) ) , and this form is denoted by B 2 ( l ) .
Now, for any neighborhood sequence B B @ ¯ that has period l, i.e., B = ( b ( 1 ) , , b ( l ) ) , and for any i N with 1 i l let B i ( l ) = ( b ( 1 ) , , b ( l ) ) defined with the elements b ( j ) = b ( j ) for j i and b ( i ) = 5 b ( i ) . In this way, only the i-th (modulo l) element of B is changed from 2 to 3 or vice versa, each other element remains the same.
Further, for any k N 0 , we say that k is even-square-bounded (ESB) if there is an even j N 0 such that j 2 k < ( j + 1 ) 2 , else k is odd-square-bounded (OSB). Clearly, the even-square-bounded and the odd-square-bounded numbers form longer and longer intervals in N 0 as they are growing. Let B 1 = ( 3 ) , B 1 = ( 2 ) , and now we define the sequences B k and B k in T as follows. For an ESB k N 0 , use B k = ( B k 1 ) 2 ( 2 k ) and B k = ( B k 1 ) 2 ( 2 k ) 2 k ( 2 k + 1 ) , else (in case k is OSB) let B k = ( B k 1 ) 2 ( 2 k ) 2 k ( 2 k + 1 ) and B k = ( B k 1 ) 2 ( 2 k ) .
In this way, to describe B k and B k , we give the 2 k long periods of B k 1 and B k 1 after B k 1 and B k 1 , respectively, such that exactly one of the 2 k + 1 -th long periods is modified at the 2 k -th position; where the modified period depends on the ESBness of k. For the sake of a better view, we show some of the first constructed values of these sequences:
B 0 = ( 3 , 3 ) ,
B 1 = ( 3 , 2 , 3 , 3 ) ,
B 2 = ( 3 , 2 , 3 , 2 , 3 , 2 , 3 , 3 ) ,
B 3 = ( 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 3 )
B 4 = ( 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 3 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 3 )
and
B 0 = ( 3 , 2 ) ,
B 1 = ( 3 , 2 , 3 , 2 ) ,
B 2 = ( 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 ) ,
B 3 = ( 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 ) ,
B 4 = ( 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 3 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 , 3 , 2 ) .
As one may observe, and as can also be seen from the construction, every element of both sequences for k 0 is in B @ ¯ , i.e., B k B @ ¯ and B k B @ ¯ for all non-negative values of k. It is also clear that B k B k 1 and B k 1 B k for every non-negative k, moreover, in both sequences, there are strict inequalities infinitely often. However, by comparing the elements of the two sequences with the same index k, from the construction, we have the following. They both have 2 k + 1 long periods, i.e., B k = ( b k ( 1 ) , , b k ( 2 k + 1 ) ) and B k = ( b k ( 1 ) , , b k ( 2 k + 1 ) ) .
Moreover, their elements b k ( i ) = b k ( i ) for every 1 i < 2 k + 1 and b k ( 2 k + 1 ) < b k ( 2 k + 1 ) for every k N 0 { 1 } (the mismatch is at the last position of the 2 k + 1 long period). Then, for each i N , let us choose a value k N with 2 k + 1 > i and let b l i m ( i ) = b k ( i ) and, further, B l i m = ( b l i m ( i ) ) i = 1 . By the previously analyzed property of the sequences B k and B k , the sequence B l i m is well-defined. It is clear that B l i m B @ ¯ holds. Now, we prove that B l i m is not periodic.
The definitions of B k and B k imply that, for each non-negative j b l i m ( i ) = b l i m ( i + n t ) with t = 2 ( j 2 ) , 1 i t and 0 n 2 2 j + 1 1 , but b l i m ( 2 ( j 2 ) ) b l i m ( 2 ( j + 1 ) 2 ) . However, there is no periodic sequence with this property.
Finally, let X = { B k | k 0 } . It is clear that X B @ ¯ . Now, suppose that B = X in B @ ¯ T for some B = ( b ( 1 ) , , b ( l ) ) . However, in this case, in B @ ¯ (such as in T) B B l i m but B B l i m ; therefore, there is i N such that b ( i ) > b l i m ( i ) . However, this would yield B k B in B @ ¯ T , if 2 k + 1 > i .
However, by B k B j , for each j 0 , which contradicts B = X in T . Thus, there is no least upper bound for X in B @ ¯ T ; therefore, B @ ¯ T , is not complete. □
Now, we are ready to complement the results of Theorem 7.
Theorem 10.
T , is not a complete lattice.
Proof. 
Observe that, in the previous proof, we demonstrated that ( B ¯ T ) is not complete as all the neighborhood sequences constructed were in T ; however, as B l i m T , the mentioned property is also inherited to T itself. □
Thus, the sets T and T have some similar properties under the partial order ⊒ as S and S with the differences mentioned.
We are closing in on results that connect the orderings * and ⊒ on the triangular grid. These results, with respect to Propositions 4 and 5, show that these relations give a more complex (theoretical) structure on the triangular grid compared with on the square grid.
Theorem 11.
Let B 1 , B 2 T . The following statements hold:
(a) 
If B 1 B 2 , then B 1 * B 2 .
(b) 
Moreover, B 1 * B 2 and B 2 * B 1 both may hold, even if B 1 B 2 and B 1 B 2 .
(c) 
Further, if none of B 1 B 2 and B 2 B 1 hold, all the following possibilities may occur:
  • both B 1 * B 2 and B 2 * B 1 hold, i.e., B 1 and B 2 have the same “speed” according to the faster relation;
  • exactly one of B 1 * B 2 and B 2 * B 1 holds, i.e., one of them is (really) faster (and not equally fast) compared with the other; or
  • none of B 1 * B 2 and B 2 * B 1 hold: B 1 and B 2 are also incomparable under the faster relation.
Proof. 
From Part (a) from B 1 B 2 , it follows that every B 2 -path is also a B 1 -path, and thus B 1 * B 2 is immediately proven.
If the class B ¯ contains more than one element, then one can easily pick two elements from there, e.g., the minimal and the maximal equivalent neighborhood sequences B m and B M , respectively. Then, clearly B m B M , but B M B m . On the other hand, since these neighborhood sequences are equivalent ( B m B M ; they provide the same distance function), both B m * B M and B m * B M hold.
Finally, the proof of part (c) is by examples.
Considering B 1 = ( 3 , 2 , 3 ) and B 2 = ( 3 , 3 , 2 ) , it is clear that none of B 1 B 2 and B 2 B 1 hold; however, they provide the same distance function with their minimal equivalent neighbor sequence B m = ( 3 , 2 , 2 , 2 , 2 , ) with exactly one value of 3 and all other values are 2. Thus, both B 1 * B 2 and B 2 * B 1 hold in this case.
Now, consider B 1 = ( 1 , 3 , 3 , 3 , 3 ) and B 2 = ( 2 , 3 , 2 , 2 , 2 ) . Clearly, B 1 B 2 and B 2 B 1 . On the other hand, the minimal equivalent neighborhood sequence of B 1 is B m = ( 1 , 3 , 2 , 2 , 2 , 2 , 2 , ) contains exactly one value of 1 and exactly one value of 3 and then only 2s. The minimal equivalent neighborhood sequence of B 2 is itself. These facts imply that B 2 is faster than B 1 but not vice versa, i.e., B 2 * B 1 and B 1 * B 2 .
To show the last possibility, let B 1 = ( 1 , 2 , 2 ) and B 2 = ( 2 , 1 , 1 ) . One may observe, e.g., on the bottom right of Figure 3 that these two neighborhood sequences are incomparable under the faster relation, i.e., B 1 * B 2 and B 2 * B 1 . □
In the next result, the equivalence relation ∼ is also involved.
Theorem 12.
Let every class B ¯ be represented by the minimal (maximal) equivalent neighborhood sequence B m ( B M , respectively) of the class. Then, the relation ⊒ provides a proper refinement of * defined on the sets of these classes of T, i.e., for any two neighborhood sequences B 1 , B 2 that are minimal (resp. maximal) equivalent neighborhood sequences of themselves B 1 B 2 implies that B 1 * B 2 . However, it may occur that B 1 * B 2 holds without B 1 B 2 .
Proof. 
Let B 1 and B 2 be two neighborhood sequences such that both of them are equivalent to their minimal (maximal, resp.) neighborhood sequences. Assume that B 1 B 2 holds. Let p and q be any two points of the grid. Then, let d ( p , q ; B 2 ) be determined by a short/the shortest B 2 -path Π . Observe that, by the definition of B 1 B 2 , the B 2 -path Π is also a B 1 path from p to q, but then d ( p , q ; B 1 ) d ( p , q ; B 2 ) clearly holds.
It is also clear from the definitions that B 1 B 2 and B 2 B 1 if and only if B 1 = B 2 , and then both B 1 * B 2 and B 2 * B 1 trivially hold.
The refinement is proper, since there are neighborhood sequences B 1 and B 2 , such that B 1 * B 2 but B 1 B 2 . Let B 1 = ( 2 , 2 , 2 , 1 ) and B 2 = ( 1 , 1 , 1 , 2 ) . Notice that each of these sequences defines a singleton equivalence class, therefore B 1 ¯ = { B 1 } and B 2 ¯ = { B 2 } , and thus B 1 , m = B 1 , M = B 1 and B 2 , m = B 2 , M = B 2 . □
We need to note that the previous proof works for any representatives of the classes. On the other hand, this last result shows also one of the good applicability features of our result, and this will lead us to the next section.

5. Applications

This paper is of theoretical nature; however, in this section, our aim is to briefly show how our results can be applied.
In [35], neighborhood sequences were used in communication networks in various scenarios. The case of periodic neighborhood sequences were studied as special cases, as it is sufficient to attach a period-long part of the neighborhood sequence along the signal. Based on our equivalence relation described in this paper, one may easily find a periodic equivalent neighborhood sequence to the original sequence (if such equivalent exists) even if the original sequence is not a periodic one.
Example 6.
Let B 1 = ( 1 , 3 , 1 , 1 , 3 , 1 , 1 , 2 , 1 , 1 , 2 , 1 , 1 , 3 , 1 , ) be a neighborhood sequence defined by b ( i ) = 1 , if i 1 ( mod 3 ) ; 3 , if i = 3 ( n 2 ) + 2 , for n N 0 ; 2 , if i 2 ( mod 3 ) but i { 3 ( n 2 ) + 2 | n N 0 } ; and 1 , if i 0 ( mod 3 ) . It is easy to see that B 1 is not periodic. However, one can see that B 1 is equivalent to B 2 = ( 1 , 3 , 1 ) , and this latter is periodic with period 3, which makes it possible to use the distance function defined by B 1 in communication schemes that require periodic neighborhood sequences. The minimal equivalent neighborhood sequence of B 1 is not periodic, but B 2 is its maximal equivalent neighborhood sequence (which we defined in this paper).
Another possible application can be in the approximation of Euclidean disks on the triangular grid: We recall here the fact that there are digital disks with the same radius that are incomparable under set theoretic inclusion on the triangular grid also yields to advanced types of approximation of the Euclidean disks based on neighborhood sequences that is impossible on the square grid: it uses the intersection of two disks; see, e.g., in [41].
To make such approximations allowed and efficient, one requires two such neighborhood sequences, let us say B 1 and B 2 , that are not in any of the relations that we have studied, i.e., none of the following is true for them: B 1 B 2 , B 1 * B 2 , B 2 * B 1 , B 1 B 2 , and B 2 B 1 . If any of these relations were true, then the intersection of the disks after each step will be the same as one of the disks based on one of the neighborhood sequences, and one will not gain anything by using the other neighborhood sequence.
From an application point of view, it is essential to be able to decide if two given neighborhood sequences are in any of the studied relations. Regarding our equivalence relation studied in Section 4.1, one can use Lemma 1 and the fact that equivalent neighborhood sequences have the same minimal equivalent neighborhood sequence.
For the faster relation, the decision is based on the comparisons of formulae of Lemma 2; finally, it is clear to see which neighborhood sequence componentwise dominates another neighborhood sequence by the definition of the relation.
On the other hand, based on our results, we may also able to decide about the faster relation based on Theorem 12: using an arbitrary representative of an equivalent class, it could happen that B 1 B 2 does not hold, even if it could hold for some of their other representatives. This causes the possibilities listed in Theorem 11. We highlight again the difference between neighborhood sequences on the square and on the triangular grid with the following example: B 1 = ( 2 , 3 , 2 , 2 , 2 ) and B 2 = ( 1 , 3 , 3 , 3 , 1 ) . If there are representatives B 1 B 1 ¯ and B 2 B 2 ¯ of these classes for which B 1 B 2 holds, then it is easy to see that it also holds for their minimal/maximal equivalent neighborhood sequences, and thus B 1 * B 2 can be more easily deducted by using these special representatives.

6. Conclusions

We investigated various relations among the neighborhood sequences on a triangular grid with various properties. There are specific equivalence and partial order relations, and we also proved a distributive lattice structure. More precisely, we characterized and described the equivalent neighborhood sequences that define the same distance function and characterized the partial order relation “faster” on the aforementioned equivalence classes.
Further, among other results, we showed that the relation “componentwise dominate” gives a lattice structure to the sets T and T of general and periodic neighborhood sequences. Based on our results, we can say that the structure of the set of neighborhood sequences on the triangular grid is proven to be more interesting than the structure of the set of neighborhood sequences on the square grid. This is based, e.g., on the fact that there are various neighborhood sequences that define the same distance function. Thus, the minimal and the newly investigated maximal equivalent neighborhood sequences play crucial roles not only in the computation of B-distances and in various applications of them but also in the description of various relations among the neighborhood sequences.
Further, by viewing the digital disks (circles) of the same radius, they do not form a well-ordered set on the triangular grid for a fixed radius [34,35,40], and the order that the elements of the neighborhood sequence takes also matters. As we also presented, our results could be applied in communication networks with a honeycomb structure in a similar way as to the examples shown in [35]. Our theoretical results can help to find appropriate neighborhood sequences with some desired properties and relations.
On the one hand, it is a future task to develop further applications in the field of image processing, computer imaging and/or in networking where the appropriate neighborhood sequences are selected based on the algebraic relations that we investigated. On the other hand, in [35], a possible generalization of the faster relation was also defined where some parameters were also added to make theoretical investigations on these relations a possible future work.

Funding

This research was partly funded by the Hungarian National Foundation for Scientific Research: OTKA F043090.

Data Availability Statement

Not applicable.

Acknowledgments

The author is thankful for the discussions on the topic with their OTKA F043090 research mates, A. Fazekas, A. Hajdu and L. Hajdu as well as for their support.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Klette, R.; Rosenfeld, A. Digital Geometry-Geometric Methods for Digital Picture Analysis; Morgan Kaufmann, Elsevier Science B.V.: San Franciso, CA, USA, 2004. [Google Scholar]
  2. Kiselman, C.O. Elements of Digital Geometry, Mathematical Morphology, and Discrete Optimization; Published by World Scientific: Singapore, 2022; 488p. [Google Scholar]
  3. Rosenfeld, A.; Pfaltz, J.L. Sequential operations in digital picture processing. J. ACM 1966, 13, 471–494. [Google Scholar] [CrossRef]
  4. Rosenfeld, A.; Pfaltz, J.L. Distance functions on digital pictures. Pattern Recognit. 1968, 1, 33–61. [Google Scholar] [CrossRef]
  5. Das, P.P.; Chakrabarti, P.P.; Chatterji, B.N. Distance functions in digital geometry. Inf. Sci. 1987, 42, 113–136. [Google Scholar] [CrossRef]
  6. Das, P.P.; Chatterji, B.N. Octagonal distances for digital pictures. Inf. Sci. 1990, 50, 123–150. [Google Scholar] [CrossRef]
  7. Das, P.P. Lattice of octagonal distances in digital geometry. Pattern Recognit. Lett. 1990, 11, 663–667. [Google Scholar] [CrossRef]
  8. Fazekas, A.; Hajdu, A.; Hajdu, L. Lattice of generalized neighbourhood sequences in nD and D. Publ. Math. Debr. 2002, 60, 405–427. [Google Scholar] [CrossRef]
  9. Nagy, B. Metrics based on neighbourhood sequences in triangular grids. Pure Math. Appl. 2002, 13, 259–274. [Google Scholar]
  10. Yamashita, M.; Honda, N. Distance functions defined by variable neighborhood sequences. Pattern Recognit. 1984, 17, 509–513. [Google Scholar] [CrossRef]
  11. Yamashita, M.; Ibaraki, T. Distances defined by neighborhood sequences. Pattern Recognit. 1986, 19, 237–246. [Google Scholar] [CrossRef]
  12. Nagy, B. Distance functions based on generalized neighbourhood sequences in finite and infinite dimensional space. In Proceedings of the ICAI’01: Fifth International Conference on Applied Informatics, Eger, Hungary, 28 January–3 February 2001; pp. 183–190. [Google Scholar]
  13. Nagy, B. Distance functions based on neighbourhood sequences. Publ. Math. Debr. 2003, 63, 483–493. [Google Scholar] [CrossRef]
  14. Nagy, B. Distance with generalized neighbourhood sequences in nD and D. Discret. Appl. Math. 2008, 156, 2344–2351. [Google Scholar] [CrossRef]
  15. Danielsson, P.E. 3D octagonal metrics. In Proceedings of the Eighth Scandinavian Conference on Image Analysis (SCIA 1993), Tromsø, Norway, 25–28 May 1993; pp. 727–736. [Google Scholar]
  16. Das, P.P. Best simple octagonal distances in digital geometry. J. Approx. Theory 1992, 68, 155–174. [Google Scholar] [CrossRef] [Green Version]
  17. Mukherjee, J. Error analysis of octagonal distances defined by periodic neighborhood sequences for approximating Euclidean metrics in arbitrary dimension. Pattern Recognit. Lett. 2016, 75, 16–23. [Google Scholar] [CrossRef]
  18. Mukherjee, J.; Das, P.P.; Kumar, M.A.; Chatterji, B.N. On approximating Euclidean metrics by digital distances in 2D and 3D. Pattern Recognit. Lett. 2000, 21, 573–582. [Google Scholar] [CrossRef]
  19. Strand, R.; Nagy, B.; Fouard, C.; Borgefors, G. Generating Distance Maps with Neighbourhood Sequences. In Discrete Geometry for Computer Imagery, DGCI 2006, LNCS 4245; Springer: Berlin/Heidelberg, Germany, 2006; pp. 295–307. [Google Scholar]
  20. Normand, N.; Strand, R.; Évenou, P.; Arlicot, A. Minimal-delay distance transform for neighborhood-sequence distances in 2D and 3D. Comput. Vis. Image Underst. 2013, 117, 409–417. [Google Scholar] [CrossRef] [Green Version]
  21. Normand, N.; Strand, R.; Évenou, P.; Arlicot, A. A Streaming Distance Transform Algorithm for Neighborhood-Sequence Distances. Image Process. Line 2014, 4, 196–203. [Google Scholar] [CrossRef] [Green Version]
  22. Normand, N.; Strand, R.; Évenou, P. Digital Distances and Integer Sequences. In Discrete Geometry for Computer Imagery. DGCI 2013; Lecture Notes in Computer Science, LNCS 7749; Gonzalez-Diaz, R., Jimenez, M.J., Medrano, B., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 169–179. [Google Scholar]
  23. Nagy, B.; Strand, R.; Normand, N. Distance Functions Based on Multiple Types of Weighted Steps Combined with Neighborhood Sequences. J. Math. Imaging Vis. 2018, 60, 1209–1219. [Google Scholar] [CrossRef]
  24. Nagy, B.; Strand, R.; Normand, N. Distance Transform Based on Weight Sequences. In Discrete Geometry for Computer Imagery-21st IAPR International Conference, DGCI 2019, LNCS 11414; Springer: Cham, Switzerland, 2019; pp. 62–74. [Google Scholar]
  25. Martin, R.; Nickson, T.; Potapov, I. Geometric computations by broadcasting automata. Nat. Comput. 2012, 11, 623–635. [Google Scholar] [CrossRef]
  26. Potapov, I. Pattern formations with broadcasting automata model. In Proceedings of the Invited Talk at Eighth Workshop on Non-Classical Models of Automata and Applications, NCMA 2016, Debrecen, Hungary, 29–30 August 2016; [email protected] 321. Österreichische Computer Gesellschaft/Austrian Computer Society: Vienna, Austria, 2016; pp. 61–73. [Google Scholar]
  27. Song, H.; Potapov, I. Polygon Approximations of the Euclidean Circles on the Square Grid by Broadcasting Sequences. In Discrete Geometry for Computer Imagery-21st IAPR International Conference, DGCI 2019, LNCS 11414; Springer: Cham, Switzerland, 2019; pp. 444–456. [Google Scholar]
  28. Conway, J.H.; Burgiel, H.; Goodman-Strass, C. The Symmetries of Things; AK Peters: Natick, MA, USA, 2008. [Google Scholar]
  29. Grünbaum, B.; Shephard, G.C. Tilings by regular polygons. Math. Mag. 1977, 50, 227–247. [Google Scholar] [CrossRef]
  30. Nagy, B. Non-traditional 2D Grids in Combinatorial Imaging–Advances and Challenges. In Keynote Talk at 21st International Workshops on Combinatorial Image Analysis, IWCIA22, to Appear in LNCS 13348; Springer: Cham, Switzerland, 2022. [Google Scholar]
  31. Deutsch, E.S. Thinning algorithms on rectangular, hexagonal and triangular arrays. Commun. ACM 1972, 15, 827–837. [Google Scholar] [CrossRef] [Green Version]
  32. Stojmenovic, I. Honeycomb Networks: Topological Properties and Communication Algorithms. IEEE Trans. Parallel Distrib. Syst. 1997, 8, 1036–1042. [Google Scholar] [CrossRef]
  33. Nagy, B. Finding shortest path with neighborhood sequences in triangular grids. In Proceedings of the ITI-ISPA 2001: Second IEEE R8-EURASIP International Symposium on Image and Signal Processing and Analysis, Pula, Croatia, 19–21 June 2001; pp. 55–60. [Google Scholar]
  34. Nagy, B. Characterization of digital circles in triangular grid. Pattern Recognit. Lett. 2004, 25, 1231–1242. [Google Scholar] [CrossRef]
  35. Nagy, B. Application of neighborhood sequences in communication of hexagonal networks. Discret. Appl. Math. 2017, 216, 424–440. [Google Scholar] [CrossRef]
  36. Nagy, B. Isometric transformations of the dual of the hexagonal lattice. In Proceedings of the sixth IEEE International Symposium on Image and Signal Processing and Analysis, ISPA 2009, Salzburg, Austria, 16–18 September 2009; pp. 432–437. [Google Scholar]
  37. Nagy, B. Distances with neighbourhood sequences in cubic and triangular grids. Pattern Recognit. Lett. 2007, 28, 99–109. [Google Scholar] [CrossRef]
  38. Grätzer, G. Lattice Theory: First Concepts and Distributive Lattices; Dover Books on Mathematics; Dover Publications: Dover, UK, 2009. [Google Scholar]
  39. Grätzer, G. Lattice Theory: Foundation; Birkhäuser: Basel, Switzerland, 2011. [Google Scholar]
  40. Nagy, B. Number of Words Characterizing Digital Balls on the Triangular Tiling. In Discrete Geometry for Computer Imagery-19th IAPR International Conference, DGCI 2016, LNCS 9647; Springer: Berlin/Heidelberg, Germany, 2016; pp. 31–44. [Google Scholar]
  41. Nagy, B.; Strand, R. Approximating Euclidean circles by neighbourhood sequences in a hexagonal grid. Theor. Comput. Sci. 2011, 412, 1364–1377. [Google Scholar] [CrossRef]
Figure 1. Various (strict) neighborhood relations (left) and a part of the grid with the symmetric coordinate system (right).
Figure 1. Various (strict) neighborhood relations (left) and a part of the grid with the symmetric coordinate system (right).
Mathematics 10 04514 g001
Figure 2. Digital disks with neighborhood sequence ( 3 , 1 ) on the left, with ( 2 ) on the right and with ( 1 , 3 ) at the bottom. The origin is marked with gray color, disks with radius 1 are dark shaded, disks with radius 2 are (dark and) medium color, while for disks with radius 3, all shaded pixels are counted.
Figure 2. Digital disks with neighborhood sequence ( 3 , 1 ) on the left, with ( 2 ) on the right and with ( 1 , 3 ) at the bottom. The origin is marked with gray color, disks with radius 1 are dark shaded, disks with radius 2 are (dark and) medium color, while for disks with radius 3, all shaded pixels are counted.
Mathematics 10 04514 g002
Figure 3. The Hasse diagrams of the partial orders ( S , * ) (top) and ( T , * ) (bottom) restricted to neighborhood sequences with periods 1, 2 and 3, from left to right, respectively.
Figure 3. The Hasse diagrams of the partial orders ( S , * ) (top) and ( T , * ) (bottom) restricted to neighborhood sequences with periods 1, 2 and 3, from left to right, respectively.
Mathematics 10 04514 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nagy, B. Equivalence, Partial Order and Lattice of Neighborhood Sequences on the Triangular Grid. Mathematics 2022, 10, 4514. https://0-doi-org.brum.beds.ac.uk/10.3390/math10234514

AMA Style

Nagy B. Equivalence, Partial Order and Lattice of Neighborhood Sequences on the Triangular Grid. Mathematics. 2022; 10(23):4514. https://0-doi-org.brum.beds.ac.uk/10.3390/math10234514

Chicago/Turabian Style

Nagy, Benedek. 2022. "Equivalence, Partial Order and Lattice of Neighborhood Sequences on the Triangular Grid" Mathematics 10, no. 23: 4514. https://0-doi-org.brum.beds.ac.uk/10.3390/math10234514

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop