Journal of Universal Computer Science, vol. 6, no. 4 (2000), 447-459
submitted: 30/4/99, accepted: 28/1/00, appeared: 28/4/00 Springer Pub. Co.
Region-based Discrete Geometry
M. B. Smyth
Department of Computing
Imp erial College
London SW7 2BZ
mbs@do c.ic.ac.uk
Abstract: This pap er is an essay in axiomatic foundations for discrete geometry in-
tended, in principle, to b e suitable for digital image pro cessing and more sp eculatively
for spatial reasoning and description as in AI and GIS. Only the geometry of convex-
ity and linearity is treated here. A digital image is considered as a nite collection
of regions; regions are primitiveentities they are not sets of p oints. The main re-
sult Theorem 20 shows that nite spaces are sucient. The theory draws on b oth
\region-based top ology" also known as mereotop ology and abstract convexity theory.
Key Words: Discrete geometry, regions, mereotop ology, convexity.
1 Intro duction
Digital top ology, as a study of the connectivity prop erties of digital images and
the spaces in which they lie, is by now a fairly well-develop ed discipline. The
approach is either graph-theoretic, as in [14], or else top ological in the strict
sense, as in [8]. See also [9] and references given there.
In our own work in this area, the emphasis has b een on an amalgamation
of digital top ology with ordinary top ology, achieved by working with a single
category in which all the spaces including the graphs exist as ob jects. In this
category TopGr, for top ological graphs the usual \continuous" spaces arise as
inverse limits of the digital spaces [15,16,17].
In the present pap er we take up the task of extending this approach to
geometry prop er. We are not aware of any previous axiomatic theory of digital
geometry, although it has o ccasionally b een prop osed that such a theory should
b e develop ed for example by Zeeman [18]. Most often what is studied is some
n
version of the \grid" mo del: in e ect, Z taken with some graph adjacency
structure, and with n restricted to b e 2 or 3.
What we shall undertake here is an axiomatic approach to the geometry of
convexity and linearity: no attempt will b e made to deal with parallels and con-
gruence for the time b eing. Abstract convexity theory, esp ecially as develop ed
by W. Prenowitz [12], is our exemplar; the problem is to adapt this material
to \discrete" geometry. See also Copp el [4] for a recent treatment of abstract
convexity theory. Both Prenowitz and Copp el provide non-trivial nite mo dels
of parts of the theory. These, however, like the nite geometries of combinatorial
theory [5], cannot b e said to resemble the spaces traditionally studied Euclidean
or non-Euclidean, nor to b e capable of b eing viewed as discrete or approximate
versions of these.
In order to obtain discrete spaces that are suitable for digital image pro cess-
ing, we shall adopt the region-based approach of what is known as mereotop ol-
ogy [2, 6, 13, 7]. Pixels or voxels and other simple convex regions are not
448 Smyth M.B.: Region-based Discrete Geometry
considered as sets of p oints, but as primitiveentities which are sub ject to a re-
lation of partial order \part of " and a symmetric binary relation expressing
closeness of two regions the \connection predicate". Our theory maythus b e
considered as a sort of fusion of mereotop ology with abstract convexity theory.
The primitives we adopt are not to o far removed from those whichhave b een
adopted in some previous work on spatial reasoning [3] and references given
there. The main di erences b etween what we are attempting here and what was
done in those works are as follows. First, we aim for a mathematically adequate
account of convexity. A simple criterion for this, often adopted in studies of
abstract convexity, is that it p ermits the derivation of abstract versions of three
famous theorems in convexity, namely those of Radon, Helly and Caratheo dory;
see for example [1]. The axioms concerning convexityhavetobechosen so as to
p ermit this. A second di erence is that we are willing to admit some p oints as
regions. We are mainly interested in discrete geometries, and it do es not greatly
matter whether the nitely many p oints which can b e found in a b ounded region
of a discrete space are treated as primitive or as nite lters of regions. For
the general theory, it would be desirable to resolve this question, however. A
further imp ortant technical di erence from many region-based theories [13,10]
is that we do not require our structures to b e Bo olean algebras; in particular,
there is no op eration corresp onding to the intersection of regions. Despite these
di erences, we think it p ossible that some development of the theory presented
here will prove to b e useful for studies in spatial reasoning.
n
Note : As a reminder, Radon's Theorem for < says that any set of n + 2 or
more p oints can b e partitioned into two disjoint subsets whose convex closures
intersect. It is enough to lo ok at the case n = 2 to grasp the meaning of the
theorem; while the case n = 1 is already of interest, as we shall see later. The
Helly and Carathe o dory theorems are closely related to the Radon theorem, but
they will not b e considered explicitly here.
2 Axiomatics
Our basic structure is a triple Q,1,, where Q is a set of regions , 1 is
a symmetric binary relation connection , and is a commutative asso ciative
op eration product, satisfying the following Axioms:
A 1 is \almost re exive":
9X:A 1 X A 1 A
The e ect of this will b e that a non-null region is connected with itself.
W
B Fusion For any collection B Q, there exists a unique region B such
that
_
X 1 B , 9B 2B:X 1 B:
W
C Distr. distributes over .
D Extension For any A; B 2 Q, there exists a region B=A such that
X 1 B=A , AX 1 B
Notice that o ccurrences of the op erator are generally omitted.
Smyth M.B.: Region-based Discrete Geometry 449
E A 1 C & B 1 D _ B=A 1 D=C AD 1 BC
E A stronger form of Axiom E
A 1 C & B 1 D B=A 1 D=C AD 1 BC
Some immediate consequences of these axioms:
Prop osition 1. The relation , de nedonQ by
*
A B 8X:A 1 X B 1 X
W
isapartial order, with respect to which is the sup operation.
Proof. Trivially, is a pre-order. It is a partial order since, if 8X:A 1 X ,
W
B 1 X , then, by uniqueness of fusion, A = B = fAg.Itis then a trivial
veri cation that fusion is the join for this partial order.
Prop osition 2. Extension distributes over join.
Proof.
_ _ _ _
X 1 A= B,X B1 A
_ _
, X B 1 A
B 2B
,9A 2A; B 2B:X B 1 A
,9A 2A; B 2B:X 1 A=B:
There follow some comments on the axioms:
W
{ Taking Q together with and ,wehave a non-unital quantale
W
{ Since wehave ,wehave a complete lattice. However, this is not required
V
to b e distributive, and has little signi cance compared with and 1.
{ It is easy to see that for any A 2 Q we have the \complement"
W
fX j:X 1 Ag, which satis es the usual join and meet conditions for
Bo olean complement. By the preceding remark, however, we do not thereby
obtain a Bo olean algebra. Moreover the complement so de ned is not in
general an involution thus it is not an ortho complement.
{ An immediate consequence of E but not of E is that, if A; B are non-
null regions, then A=B is non-null. A mo del based on a closed b ounded
subset of Euclidean space will fail Axiom E see the interpretation 3 to
be given in a moment, but may nevertheless be of interest for discrete
geometry. In the pro ofs b elowwehave taken care to use E only when E is
not sucient. Prenowitz works with p oint-based axioms corresp onding to
Axiom E , whereas Bryant& Webster have, in e ect, only the weaker form
E.
{ It may sometimes be convenient to drop the uniqueness requirement from
Axiom B. The join is then a sp eci ed op eration satisfying the prop erty given
in B, and the order Prop. 1 can only b e claimed to b e a pre-order. In the
conclusions of one or two of the subsequent Prop ositions, equality should b e
replaced by equivalence.
450 Smyth M.B.: Region-based Discrete Geometry
{ The main remaining axiom to b e considered is an Axiom of Order see Sec-
tion 5.
We next consider some of the main intended mo dels of the axiom system.
n
1. Euclidean Q is: the subsets of < . A 1 B means: A meets B . The pro duct
is de ned in the rst instance for p oints:
x; y if x 6= y
x y =
x if x = y
By distributivity, the de nition extends to arbitrary regions. Taking x y as
the set of p oints strictly b etween x and y the op en interval, in case x 6= y ,
is Prenowitz' preferred interpretation. Notice that if, in apparent conformity
with this, wewere to try to take x x as empty,we should lose asso ciativity
consider x x y . A variant corresp onding to what is usually done in
de ning \interval convexities" is to de ne x y as the closed interval [x; y ]
in all cases.
2. More in the spirit of mereotop ology would b e to take Q as the set of regular
n
op en subsets of < , and to take A 1 B as meaning that the closures of
A,B meet, with the pro duct de ned exactly as b efore it clearly makes no
di erence here whether intervals are taken as op en or closed. Of course with
this interpretation wewould obtain a complete Bo olean algebra.
3. The rst interpretation may be generalized by starting with an arbitrary
convex subset C of the Euclidean n-space while de ning 1 and exactly
as b efore. It is easily seen that, if C is op en, all the axioms are satis ed. If
C is not op en, then we do not in general haveE , but only the weak form
Axiom E . This is b ecause, if x is an extreme p ointof C , x=y may b e empty
assuming Prenowitz' interpretation of . In any case, let us denote by GC
the geometry de ned in this way. Supp ose next that S is a nite subset of
C . Then wemay consider the geometry \generated" by S : the collection of
all the regions that may b e obtained from S by rep eatedly taking pro ducts
and extensions and union, with 1 de ned as in GC . From the p ointof
view of discrete geometry, the question that is now of interest is, whether
the geometry so obtained is itself nite that is, whether only nitely many
distinct regions arise in the construction. We shall return to this later for
n
the case C = < .
We may remark that, implicit in the third interpretation just given, is the
0 0 0
notion of a subgeometry. Given two geometries Q; 1; , Q ; 1 ; , the rst
0
is said to b e a subgeometry of the second if Q Q and the connection, pro duct,
0
extension and join fusion of Q coincide with the restrictions of those of Q to
n
Q.Thus, if C is a prop er convex subset of < , then GC is not a subgeometry
n
of G< , since the extension op eration di ers. On the other hand, the question
just raised is a question as to whether a certain nitely generated subgeometry
of GC is nite.
Note that the subgeometry construction is one which can lead to geome-
tries lacking uniqueness of fusion \pre-ordered geometries", since there may
be to o few elements remaining to make via 1 and Axiom B all the distinc-
tions b etween regions which hold in the larger geometry. See comments following Theorem 20.
Smyth M.B.: Region-based Discrete Geometry 451
3 Development of geometry I
In developing geometry on the basis of an axiomatization of the ab ovetyp e, a
key role is played by a series of basic lemmas, rst develop ed by Prenowitz [12].
A similar development is given by Copp el [4]. With their aid, much geometric
reasoning is reduced to simple algebraic manipulations though in a completely
di erent fashion from the way in which this happ ens in ordinary analytic geome-
try. The lemmas as will b e seen resemble the formulas found in the elementary
calculus of fractions, though with equality usually replaced by inequality .
The development of these lemmas given by the authors just cited is, of course,
entirely p oint-based regions are assumed to b e sets of p oints. Wehave to ensure
that, by adopting a region-based approach, wehave not lost anything essential.
Lemma 3. A=B =C = A=B C
Proof. For any region X
X 1 A=B =C , XC 1 A=B
, XCB 1 A
, X 1 A=B C
all by the Extension axiom. Due to the Fusion axiom, this proves the result.
Lemma 4. Assume A 6=0. Then
a B AB=A
b B AB =A
c B A=A=B
Proof. Supp ose that X 1 B . Then X=A 1 B=A, AX 1 AB , and A=X 1 A=B
all by Axiom E . Invoking Extension, we immediately obtain a,b,c resp ec-
tively.
Lemma 5. AB=C AB =C
Proof. Supp ose X 1 AB=C . Then
X=A 1 B=C Extension
XC 1 AB Axiom E
X 1 AB =C
Lemma 6. A=B C=D AC =B D
Proof. A=B C=D A=B C=D Lemma 5
C A=B =D Lemma 5, monotonicity
AC =B D Lemma 3.
The monotonicityof and =, whichwehave not troubled to sp ell out, is of
W
course a consequence of the distributivityof over these op erators.
The following lemma is proved similarly to the preceding one:
452 Smyth M.B.: Region-based Discrete Geometry
Lemma 7. A=B =C=D AD =B C
The pattern of these pro ofs should be clear. Where, in p oint-based work,
an inclusion A B is proved via the assumption \x 2 A", the region-based
argument will yield A B , via the assumption \X 1 A".
4 Convex and Linear Regions
A region A should b e considered convex if, whenever U ,V A, U V A. This
simpli es to:
De nition 8. Region A is convex if AA A.
Theorem 9. The col lection of al l convex regions is closed under arbitrary meets
and directed joins.
Proof. Let X b e the meet of an arbitrary collection A of convex regions. Then
V
XX AA, for every A 2A. Hence XX AA X . Next, let Y b e the
A2A
W
BC . join of a directed collection B of convex regions. By Distr., YY =
B;C 2B
W
Since B is directed, YY DD Y .
D 2B
This means that convex closure or hull can be de ned, with the usual
general prop erties, including Scott-continuity.
As for the other op erators, it is immediate that the pro duct of nitely many
convex regions is convex. Only a little less trivial is:
Prop osition 10. A,B conv ex A=B convex.
Proof. AA A and BB B . Hence by Lemma 6,
A=B A=B AA=B B A=B :
A linear region should contain, with any parts U ,V , not only what lies b e-
tween them, but also what lies in the extension V=U. On simplifying, we get:
De nition 11. Region A is linearly closed provided that A is convex and
A=A A.
In deference to analytic geometry, it might be preferable to dub these re-
gions ane rather than linear. Prenowitz, and also Bryant and Webster, have
\linear", Copp el has \ane". Wehave the same general prop erties for linear-
ity, including the existence of the linear hull of any region, as indicated ab ove
for convexity.For the sp ecial op erators, however, matters are a little di erent.
Evidently, the pro duct of two linear regions need not b e linear. For the exten-
sion, wehave a p ositive result, corresp onding to Theorem 6.5 of Prenowitz [12].
Prenowitz' pro of including the pro of of his Lemma 6.5 cannot b e adapted to
Smyth M.B.: Region-based Discrete Geometry 453
our setting, as it dep ends essentially on the use of p oints, and sp eci cally on the
assumption that, if two regions are connected, they have a p oint in common. So
we provide the region-based pro of here in detail:
Lemma 12. Let A,L beregions such that A 1 L.IfL is convex, then L L=A.
If L is linear, then also L A=L.
Proof. Supp ose that L is convex, and X 1 L. Then Axiom E XA 1 LL L.
So Extension X 1 L=A. This shows that L L=A. Next, assume L linear,
and supp ose that X 1 L. Then Axiom E A=X 1 L=L L. Hence A 1 XL,
and so X 1 A=L.
Theorem 13. Suppose that A,B are linear, and A 1 B . Then A=B is linear.
Proof.
A=B =A=B A=B =B=A=A=B Lemma 12
A=B =BB=AA Lemma 7
A=B =B=A
A=B Lemma 7
This shows that A=B is linear.
An easy application of Lemma 7 also yields:
Prop osition 14. If A is convex, then A=A is linear.
Prop osition 15. The linear hul l of a nite set fA ,...,A g of convex regions is
1 n
given by R=R, where R = A ::: A .
1 n
Proof. The region R /R is linear by the preceding Prop osition, and is evidently
contained in any linear region which contains all the A . It remains only to show
i
that A R=R for each i. Wlg, let us show this for i = 1. Supp ose X 1 A .
i 1
Then for any P in particular, for P = A ::: A , XA P 1 A A P . Since
2 n 1 1 1
A A A , XA P A P . Hence X 1 A P /A P = R /R .
1 1 1 1 1 1 1
A simple result which will b e useful in Section 6 is:
Prop osition 16. If A is a convex region, and B; C are any regions, then A=B
A=C A=B C and B=A C=A BC=A.
Proof. By two uses of Lemma 5 and one use of Lemma 3:
A A A=B A AA=B A
B C C C BC The other part is similar.
454 Smyth M.B.: Region-based Discrete Geometry
As already mentioned, we do not at present see how to eliminate p oints alto-
gether from the development. In particular, we seem to need them for formulating
the Axiom of Order to b e considered shortly, and thus for establishing the key
theorems which dep end on it. It is also true that, in describing some of the ge-
ometries in which we are mainly interested in this pap er, namely grid mo dels
for digital image pro cessng, we nd it convenient to admit, b esides the pixels,
also the edges and vertices of pixels as primitive regions. We shall regard p oints
as linear regions which are as \small as p ossible" without b eing null. There
are at least two notions of smallness which suggest themselves here. The rst
is complete primality: if a p oint is covered by a collection of regions, it must
b e contained in one or more of those regions. The second is that, if a p ointis
connected with a region, it must b e contained in that region. Fortunately, these
are equivalent:
Prop osition 17. Aregion A is a complete prime if and only if it satis es:
8X:A 1 X A X 1
W
Proof. ONLY IF: Supp ose A B9B 2B. A B and A 1 X .Take
B as
fX g[ fY jY 1 A & :Y 1 X g
W
Clearly, A B . But X is the only member of B which can contain A .
IF: Assume 1. Then
_
A B 9B 2B:A 1 B
9B 2B:A B
De nition 18. A p oint is a non-null completely prime linear region.
Prop osition 19. If P isapoint, then PP = P=P = P .
5 Development of geometry II
This section on further development of geometry could of course b e rather exten-
sive, but will on the contrary b e very brief. This is b ecause we shall have recourse
to p oints, after whichwe need do little more than showhow the developmentin
[12]or[4]may b e imitated.
The Axiom of Order, as usually formulated, concerns the order of p oints on a
line \Given any three distinct collinear p oints, exactly one of them lies b etween
the other two", say. This is problematic for us in a number of ways. Ignore, for
the moment, scruples ab out the emphasis on p oints. A line is presumably to b e
de ned as the linear hull of two distinct p oints. Now, it mayvery well happ en,
in a non-trivial mo del of our theory, that lines - even, all lines - have only two
p oints. In that case,the Axiom of Order as formulated ab ovewould b e vacuous.
We can try to reformulate the Axiom so that it refers to a \linear arrangement"
of regions that need not be p oints; but this seems hard to do, in a way that
preserves the p ower of the axiom in p oint-based work.
Smyth M.B.: Region-based Discrete Geometry 455
A partial solution maybe found by switching attention from lines to rays.
Aray is a half-line. Mention of lines, however, is inessential in treating rays. If
p,q are any two p oints,the ray from p in the direction of q may be expressed
as p q=p. Notice that the case p = q is not excluded. If p = q , we get the
\degenerate" ray p p=p = p. In the expression for a ray, the p oint q may with
advantage b e replaced by an arbitrary region A. The p oint p may b e replaced,
for reasons whichwe cannot go into here, by an arbitrary linear region, yielding
the notion of a generalizedray or cone as a region of the form MA/M , where
M is linear. The b eautiful theory of rays and spherical geometries of Prenowitz
[11,12] can now, if desired, b e adapted to the region-based style. But here we
haveintro duced rays only in order to address the Axiom of Order.
One of several principles shown by Prenowitz to b e equivalent in p oint-based
theory to the Axiom of Order is the following: for any p oint P and region for
Prenowitz: set of p oints R ,
PR=P = PR _ R _ R=P
In words: the ray PR/P consists of what lies b etween P and R , R itself, and the
extension of R away from P . Notice that the p oint P cannot b e replaced even
by a general linear region. Taking P as a line, the resulting statementwould b e
false already in Euclidean geometry. This \ray principle" is taken as the basic
form of the Order Axiom by Bryant&Webster [1], and this is what we shall do
to o. But whereas for these authors the rayversion is simply a variant of the line
version, for us it is essentially stronger: in a region geometry, the ray principle
mayvery well b e non-vacuous even when every line contains only two p oints.
We shall pro ceed by giving an indication little more than a hint of how
the Order Axiom, in its rayversion, underlies p owerful theorems on convexity.
Notice that in the equation emb o dying the ray principle, the \quotient" on the
left has the factor P o ccurring in b oth numerator and denominator, whereas
no such rep eated factor o ccurs on the right hand side. Indeed, in the geometric
calculus, the principle functions mainly as an aid in eliminating rep etitions of
factors of this kind. An imp ortant application is in connection with the formula
for linear hull Prop osition 15. In case the A are p oints, rep eated use of the
i
principle enables an expansion of the linear hull formulas to b e given, in which
no rep eated factors o ccur. For the details of this, we refer to [12, 1]. Here we just
consider the case n=2. It is an easy exercise to show, using the ray principle,
that, if A ,A are p oints, wehave:
1 2
A A =A A = A =A _ A _ A A _ A _ A =A
1 2 1 2 1 2 1 1 2 2 2 1
This clearly implies the following : if A is any third p oint on the line through
3
A ,A , then one of the three p oints lies b etween the other two. This shows the
1 2
connection with the line version of the Order Axiom. More interesting in the
present context, however, is that it shows that the ray principle, applied to the
linear hull formula for two p oints, yields what is essentially the 1-dimensional
Radon theorem. In a similar way, the expansion for n p oints enables the full
Radon theorem to be obtained. With a little more work, general versions of
the Helly and Caratheo dory theorems may b e obtained as well: see [12], or the
extremely concise [1], for details.
456 Smyth M.B.: Region-based Discrete Geometry
6 Finite Mo dels
It might b e asked: whydowe not admit intersection of regions as a primitive
op eration, enabling the whole development to b e simpli ed? The answer is that
wewould thereby lose the p ossibilty of nite mo dels of the kind in whichwe
are interested.
Supp ose that we construct a geometry on a nite grid 1 x 3 suces, with
unit squares. Let A,B ,C ,D b e the p oints 0,0, 3,0, 0,1, 3,1 resp ectively.
1 1
The segments AD , BC meet the ordinates through 0,1, 0,2 at, say, U ,V .
1 1
The line through U ,V meets AC , BD at p oints which we shall call C ,D .
1 1 2 2
Taking C ,D in place of C , D , and rep eating the construction, we obtain
2 2 1 1
sequences of distinct p oints C ,D , which we are forced to admit into the
i i
geometry if intersections of regions are accepted as regions.
Now, it must b e observed that, in certain circumstances, the extension op-
erator is capable of generating in nite sequences in a similar way.To see this,
let us start with the 2 x 3 grid, and let the p oints A,B ,C ,D ,U ,V b e the p oints
1
0,1, 3,1,3,1.5, 0,0, 1,1,2,1 resp ectively. The p oint region C do es not
1
b elong to the mo del, but the region R =VBC do es, as it is V /AD V . It is easily
1 1
seen that the region R =V /U /R is VBC , where C lies strictly between B
2 1 2 2
and C . Continuing in this manner, we obtain a strictly decreasing sequence of
1
regions R .
i
In terms of the discussion in Section 2 third interpretation, what this argu-
ment has shown is that, if K is the rectangle with vertices 0,0,3,0,3,2,0,2,
then the subgeometry of GK generated by the unit grid p oints of K is in -
nite.
To approach the construction of nite mo dels, we b egin by considering a
prop erty enjoyed by some, but not all, geometries:
F A=B C A=B A=C & BC=A B=A C=A
n
It is easy to see that F is satis ed by the geometry G< it suces to verify
the prop ertyin the case that A,B ,C are p oints. By Prop osition 16, we know
that, in case A is convex, the inequalities in prop erty F maybe replaced by
equalities. The signi cance of this is that any term built up, using * and /, from
constants A ,...,A representing convex regions can b e reduced to a pro duct of
1 n
terms in A ,...,A built using / alone.
1 n
The preceding reduction drives o ccurrences of / inwards. A further reduction
will drive them to the left. Indeed, wehave
A=B =C = A=B C Lemma3 = A=B A=C F
this b eing a pro duct of terms of the form A/X , eachhaving fewer o ccurrences of /
than A/B /C . Let an expression A /A /.../A be understo o d as
1 2 n
A /A /.../A ... asso ciation to right. Then by an easy ordinal induction
1 2 n
based on the ab ove reductions, wehave:
Theorem 20 Normal form. Let T be any term built, using * and /, from
constants representing convex regions. Then, assuming that Property F holds,
T can be expressed asaproduct of terms of the form A /A /.../A A ,...,A
1 2 n 1 n
constants
Smyth M.B.: Region-based Discrete Geometry 457
We can now obtain our main result, namely the existence of nite geometries
n
whichmay b e considered to approximate < .
n
Theorem 21. Let S be a nite set of points in < . Then the subgeometry of
n
G< generatedbyS is nite.
Proof. In view of the Normal Form Theorem, wehave only to show that there are
no more than nitely many distinct regions A /A /.../A , where A ,...,A 2 S .
1 2 n 1 n
We adopt the following notation: A is the ray with end-p oint A in the di-
i
ijk