OrbiCG/
Workshop on Computational Geometry
INRIA Sophia Antipolis - Méditerranée, 8-10 December 2010
Room Euler violet
Co-organized by
ANR Triangles and
INRIA Associate team OrbiCG
This two days workshop will address various aspects of computational
geometry
and its applications. It will be followed on December 10 by the
PhD defense of
Manuel Caroli.
Print the practical informations
Participants list
Schedule
december 8
|
10:00 - 10:30
|
Welcome and coffee
|
10:30 - 12:00
|
Presentations
|
10:30 - 11:15
|
Curvature and combinatorics of triangulations
John Sullivan (Technische Universität Berlin)
There is a rich interplay between combinatorial and differential
geometry. For 2-manifolds, the Gauss-Bonnet theorem gives a
clear connection between combinatorics (say of a triangulation) and
curvature/topology. This picture also has new applications giving a more
detailed classification of surface triangulations. We will also examine
to what extent something similar can be done for 3-manifolds.
slides.
|
|
11:15 - 12:00
|
The Poincaré dodecahedral space
Gert Vegter (Johann
Bernoulli Institute for Mathematics and Computer Science, Groningen),
Rien van de Weijgaert
(Kapteyn Institute, Groningen)
For ages people have wondered what the shape of the universe could be.
Whereas the geometry of the Universe has, by now, been established by
Einstein's theory of relativity, this theory does not determine the
exact 3-manifold that is physical space. Enquiry into the global shape
of the Universe, more specifically its topology, has been ongoing for
about a century now, but recently some new evidence emerged. In 2003,
data from the WMAP experiment led Jean-Pierre Luminet to propose a
multiconnected spherical space model called the Poincaré dodecahedral
space (PDS).
We discuss the meaning of this claim from both a mathematical and a
physical perspective. Specifically, we will discuss the relevant
physical and cosmological evidence on the topological structure of the
Universe and investigate whether there is any possibility for a
PDS. Also, the PDS will be shown to be a spherical space with a
dodecahedral fundamental domain.
paper.
slides.
|
|
12:00 - 14:00
|
Lunch
|
14:00 - 15:30
|
Presentations
|
14:00 - 14:45
|
Shortest cut graph of a surface with prescribed vertex set
Éric Colin de
Verdière (ENS)
A cut graph of a surface is an embedded graph that splits the
surface into a topological disk. Erickson and Whittlesey [SODA 2005]
found a very nice greedy algorithm to compute the shortest one-vertex
cut graph of a combinatorial surface. This talk will show how their
algorithm extends to the computation of a shortest cut graph with a
given vertex set. Moreover, a simpler proof will be given, also
revealing that the algorithm actually computes a minimum-weight basis
of some matroid.
paper .
slides.
|
|
14:45 - 15:30
|
Sculpting shapes of arbitrary topology in real-time
Lucian Stanculescu (Liris)
I will present a method for sculpting a watertight and manifold
triangular mesh which can accommodate topology changes. The key idea
is to maintain a homogeneous sampling of the surface throughout the
deformation, thus allowing the detection of topology changes from
simple collisions between spheres. The example image shows an object
sculpted from a sphere.
slides.
|
|
15:30 - 16:00
|
Coffee Break
|
16:00 - 17:30
|
Presentations
|
16:00 - 16:45
|
Linear Cameras
Jean Ponce (ENS)
I will address in this talk the problem of characterizing a general
class of cameras under reasonable, "linear" assumptions. Concretely, I
will use the formalism and terminology of classical projective
geometry to model cameras by two-parameter linear families of straight
lines---that is, degenerate reguli (rank-3 families) and
non-degenerate linear congruences (rank-4 families). This model
captures both the general linear cameras of Yu and McMillan and the
linear oblique cameras of Pajdla. I will show that it affords a simple
classification of all possible geometric camera configurations. From
an analytical viewpoint, our model also provides a simple and unified
methodology for deriving general formulas for projection and inverse
projection, triangulation, and binocular and trinocular geometry. In
addition, it is possible to generalize Pajdla's model of oblique
cameras by certain projective maps to all linear cameras, and
introduce a natural notion of intrinsic parameters for these cameras,
further simplifying the corresponding single- and multi-view
formulas. I will conclude by showing a prototype of a parabolic
camera.
Joint work with Guillaume Batog and Xavier Goaoc at INRIA Nancy Grand Est.
|
|
16:45 - 17:30
|
Balls, sticks, triangles and molecules
Frédéric Cazals (INRIA)
``Do you think that you will succeed in modeling molecules using
balls, sticks and triangles?'' This is the admittedly provocative
question I was asked by a molecular biologist during a workshop at the
NIH in Bethesda in 2008, and this talk will provide two illustrations
of geometric constructions indeed improving our understanding of
correlations between biophysical and geometric (structural) properties
of macro-molecules.
In the first part, we shall consider the problem of describing and
comparing protein interfaces. We shall see that a sprinkle of affine
Voronoi diagrams mixed with elementary ideas form Morse theory and
baked with graph matching algorithms allows one to precisely
investigate macro-molecular recognition through the atomic-modeling of
interfaces.
In the second one, we shall switch scales and consider the biggest
protein assembly known to date in Eukaryotic cells, namely the Nuclear
Pore Complex. I shall explain how the alpha-shape associated to a
compoundly-weighted Voronoi diagram allows one to make a quantitative
assessment of qualitative models of the NPC proposed by structural
biologists.
slides.
|
|
december 9
|
09:30 - 10:15
|
Presentations
|
09:30 - 10:15
|
Asymptotic properties of the typical Poisson-Voronoi cell and other
random polytopes
Pierre Calka
(Université de Rouen)
Abstract. In this talk, we define several models of random polytopes in
the Euclidean space. Examples include the convex hull of independent and
uniformly distributed points in a fixed convex set and the typical
Poisson-Voronoi cell (chosen uniformly among all cells from the Voronoi
tessellation generated by a homogeneous Poisson point process). First we
show distributional properties of some of their characteristics. Then we
focus on the description of their asymptotic shape. In particular, we
develop a local scaling theory for the boundary itself of the random
polytope, seen as a random process. We obtain variance asymptotics,
invariance principles and extreme value results.
slides.
|
|
10:15 - 10:45
|
Coffee Break
|
10:45 - 12:15
|
Presentations
|
10:45 - 11:30
|
Smooth analysis of the convex hull of points
Olivier Devillers (INRIA)
Assume that Y is a noisy version of a point set X in convex
position. How many vertices does the convex hull of Y have, that
is, what is the number of extreme points of Y?
We consider the case where X is an (ε,κ)-sample of a
sphere in Rd and the noise is random and uniform:
Y is obtained
by replacing each point x in X by a point chosen uniformly at
random in some region S(x) of size a around x. We give
upper and lower bounds on the expected number of extreme points in
Y when S(x) is a ball (in arbitrary dimension) or an
axis-parallel square (in the plane). Our bounds depend on the size
n of X and a, and are tight up to a polylogarithmic
factor. These results naturally extend in various directions (more
general point sets, other regions S(x)...).
We also present experimental results, showing that our bounds for
random noise provide good estimators of the behavior of
snap-rounding, that is when Y is obtained by rounding each
point of X
to the nearest point on a grid of step a.
We describe how the technique generalizes
to obtain (up to polylog factor)
tight upper bound on size of the 3D Delaunay triangulation when the input
points are randomly perturbed.
paper .
slides.
|
|
11:30 - 12:15
|
Convex hulls of spheres and convex hulls of convex polytopes lying on parallel hyperplanes
Menelaos Karavelas
(University of Crete)
Given a set Σ of spheres in Ed, with
d ≥ 3
and d odd, having a fixed number of m distinct radii
r1,r2,...,rm, we show that the worst-case
combinatorial complexity of the convex hull CHd(Σ) of
Σ is
Θ(∑1≤i ≠ j ≤m
ninjfloor(d/2)),
where ni is the number of spheres in Σ with
radius
ri.
Our bound refines the currently best upper bound on the worst-case
combinatorial complexity of CHd(Σ) for all
odd
d ≥ 3.
To prove the lower bound, we construct a set of
n1+n2
spheres in Ed,
with d ≥ 3 odd, where ni spheres have radius
ri, i=1,2, and
r2 ≠ r1, such that their convex
hull has combinatorial complexity
Ω(n1n2floor(d/2)
+ n2n1floor(d/2) ) . Our construction is
then generalized to the case where the spheres have m
distinct radii.
For the upper bound, we reduce the sphere convex hull
problem to the problem of computing the worst-case combinatorial
complexity of the convex hull of a set of m d-dimensional
convex polytopes lying on m parallel hyperplanes in
Ed+1, where d ≥ 3 odd, a problem which is of
independent interest.
More precisely, we show that the worst-case combinatorial complexity
of the convex hull of a set
P1,P2,...,Pm
of m d-dimensional convex polytopes lying on m parallel
hyperplanes of Ed+1 is
O(∑1≤i≠j≤mninj
floor(d/2)),
where ni is the number of vertices of Pi.
This bound is an improvement over the worst-case
bound on the combinatorial complexity of the convex hull of a point
set where we impose no restriction on the points'
configuration; using the lower bound construction for the sphere
convex hull problem, it is also shown to be tight for all odd
d ≥ 3.
Finally, we: (1) briefly discuss how to compute convex hulls of
spheres with a fixed number of distinct radii, or convex hulls of a
fixed number of polytopes lying on parallel hyperplanes; (2)
how our tight bounds for the parallel polytope convex hull problem,
yield tight bounds on the combinatorial complexity of the Minkowski
sum of two convex polytopes in Ed; and (3) state some
open problems and directions for future work.
slides.
|
|
12:15 - 14:00
|
Lunch
|
14:00 - 15:30
|
Presentations
|
14:00 - 14:45
|
Exact Voronoi Diagram of Arbitrary Lines in Space - with fast point location
Michael
Hemmer (Tel Aviv University)
For a set of objects the Voronoi Diagram (VD) is the partition of
space into cells, where each cell is defined as the set of all points
that are closer to a particular object than to all other objects. The
VD is a fundamental data structure in Computational Geometry. We
introduce a new, efficient, and complete algorithm, and its exact
implementation, to compute the Voronoi diagram of lines in
three-dimensions. This is a major milestone towards the robust
construction of the Voronoi diagram of polyhedra. Following the exact
geometric computation paradigm, it is guaranteed that we always
computes the mathematically correct result. The algorithm is complete
in the sense that it can handle all configurations, in particular all
degenerate ones. The algorithm
requires O(n3+ε) time and space,
where n is the number of lines. The Voronoi diagram is
represented by a data structure that permits answering point location
queries in O(log2n) expected time. The
implementation employs the CGAL packages for constructing arrangements
and lower envelopes on parametric surfaces together with advanced
algebraic tools. The talk will also report on recent preliminary work
extending the algorithm towards polyhedra.
|
|
14:45 - 15:30
|
3D real time visualisation in a rich internet application
Fabien Cellier (Liris)
The GIS softwares are commonly used in cartography, urban planning and land surveying.
As in many other areas of computer science, new applications embedded in web browsers have emerged recently and we now have 3d clients for the web that can display models at resolutions superior to one triangle per pixel.
This talk will present the latest advances in visualization of digital terrain model and how these algorithms can be adapted to the context of new browsing devices (typically smartphones) and new CPU/GPU hardware architectures.
slides.
|
|
15:30 - 16:00
|
Coffee Break
|
16:00 - 16:45
|
Presentations
|
16:00 - 16:45
|
Persistent Alpha Shape Topology of the Cosmic Web
Rien van de Weijgaert
(Kapteyn Institute, Groningen)
Over the past decades a clear paradigm has emerged as large redshift
surveys opened the window onto the distribution of matter in our Local
Universe: galaxies, intergalactic gas and dark matter exist in a wispy
weblike spatial arrangement consisting of dense compact clusters,
elongated filaments, and sheetlike walls, amidst large near-empty void
regions. The Cosmic Web is the fundamental spatial organization of
matter on scales of a few up to a hundred Megaparsec, scales at which
the Universe still resides in a state of moderate dynamical evolution.
To understand the topological structure of the Cosmic Web, we study
the alpha shapes defined by the Megaparsec Cosmic Web. Alphashapes
were introduced by Edelsbrunner to formalize the concept of "shape"
for a spatial point dataset. A large value of alpha corresponds to the
convex hull of the dataset, while as alpha shrinks the alphashape
assumes cavities which may join to form tunnels and voids. A related
concept is the alpha complex, which for each alpha is the
corresponding subcomplex of the three-dimensional Delaunay
tessellation.
We have studied the alpha complex of the cosmic weblike point
patterns, in order to assess the signature of filaments, walls and
voids. The topology of the corresponding alpha complexes is
characterized by Betti numbers. The physical interpretation is
determined from a range of cellular point distributions. The findings
from the Voronoi clustering models is used to analyze the outcome of
cosmological N-body simulations and the SDSS galaxy redshift survey.
|
|
december 10
|
9:30 - 11:45
|
PhD defense
Triangulating Point Sets in Orbit Spaces
Manuel Caroli
Referees:
Kurt Mehlhorn and
John Sullivan
Defense Committee:
Éric Colin de Verdière,
Menelaos Karavelas,
Jean-Marc Schlenker,
John Sullivan,
Monique Teillaud (advisor),
Gert Vegter, and
Andreas Fabri (invited).
Abstract:
In this work we discuss triangulations of spaces of different topology defined by given point sets. We propose a general definition, valid for different classes of spaces, and an algorithm to compute Delaunay triangulations of these spaces. We provide an implementation for the specific case of the three-dimensional flat torus.
The work is originally motivated by the need for software computing three-dimensional periodic Delaunay triangulations in numerous domains including astronomy, material engineering, biomedical computing, fluid dynamics etc. Periodic triangulations can be understood as triangulations of the flat torus. We provide a definition and develop an efficient incremental algorithm to compute Delaunay triangulations of the flat torus. The algorithm is a modification of the incremental algorithm for computing Delaunay triangulations in Ed. Unlike previous work on periodic triangulations we avoid maintaining several periodic copies of the input point set whenever possible. Also the output of our algorithm is guaranteed to always be a triangulation of the flat torus.
We provide an implementation of our algorithm that has been made
available to a broad public as a part of the Computational Geometry
Algorithms Library CGAL (www.cgal.org).
We generalize the work to a more general class of flat quotient
spaces as well as quotient spaces of constant positive curvature. We
furthermore consider the double torus as a representative of the much
richer class of quotient spaces of constant negative curvature.
|
|
Registration and Accommodation
Print the practical informations
You can register for the
workshop by sending a mail to
the organizers:
Click here to register.
Registration is mandatory, do it as soon as possible and anyway before
December 1st.
There is no registration fee.
Participants have to pay for their lunch
at Inria (price is 7.06 euros (7.36 with a coffee), it is a good idea
if you can have exact change).
Suggested Hotels
The following hotels in Antibes are located at walking distance
from both a bus station to INRIA, and the scenic center of Antibes, as
well as the sea. We refer to the hotels web sites for more details.
Chrys Hotel***,
Le Collier**,
Hôtel de l'Étoile**,
Hôtel Josse***.
The latter is located in Juan les Pins,
in front of the sea, but further away from bus stations.
If you prefer an hotel in Sophia, (less fancy place, almost nothing to
do
at walking distance but going to INRIA) you can go to
Le Relais**.
How to get to INRIA Sophia Antipolis Méditerranée
(map)
From Nice and/or the airport take the bus
Line230
which provides regular shuttles between Nice, the airport
Nice-Côte d'Azur and Sophia Antipolis (40 mn). Ask the driver
for the INRIA bus stop.
From Antibes take the Envibus
(the name of the bus stop is Gare Sncf - Passerelle).
On Line 11
get off at the bus stop INRIA.
On Line 1
get off at the bus stop IUT.
On Line 100
get off at bus stop Templiers, then walk a little bit.
(Get out of Antibes train station, then cross the tracks by the
"passerelle" above the tracks to get the to bus stops "SNCF
Passerelle",
after crossing the passerelle, the place for the 1/11/100 is at left.
WARNING: if you use the 100 bus, after 9:46 it goes to Antibes center
first, and then goes to Sophia directly, so it may be written on the
bus "Antibes place de Gaulle" instead of "Sophia Gare routière".)
Look at the map of the bus stops
around INRIA.
Organization
This workshop is organized by
Olivier
Devillers,
Caroline French,
and
Monique Teillaud, INRIA.