... this page is part of the Web Site of George North ...

Homework Supplement

Point Location by: George North
439-68-5643

Computational Geometry
CSCI 6640
Fall 1995

Review -- Geophysical parametrization 3
Planar Point Location Bibliography 6
Geophysical parametrization and interpolation of irregular data using natural neighbours -- by: Malcom Sambridge, Jean Braun,
and Herbert McQueen 13
Geophysical parametrization and interpolation of irregular data using natural neighbours

by:
Malcom Sambridge, Jean Braun, and Herbert McQueen,

A SUMMARY by: George North

In searching the Internet for a paper on which to base this report, I found a site "Geometry in Action."
(http://www.ics.uci.edu/~eppstein/geom.html).
Its purpose is to present various areas in which ideas from computation geometry are used to solve real world problems. Information regarding many areas are presented (with links) including: architecture, astronomy, cartography, geographic information systems, computer aided design, computer aided manufacturing, constraint solving, data mining and multidimensional analysis, graph drawing, computer graphics, interpolation, medical imaging, mesh generation, metrology, molecular modeling and nanotechnology, origami, robot motion planning, signal processing, textile layout, typography, video game programming, computer vision, VLSI design, Voronoi diagrams, windowing systems. It was interesting browsing these many sites. The paper I chose, like most, is very technical in its presentation. Admittedly, I cannot follow many of its topics. Please note that the graphics included here are borrowed from the above mentioned web pages.

A goal of Geophysical parametrization is the estimation of mineral reserves in a deposit using information obtained from bore holes and modeling crack patterns in basalt due to contraction on cooling. Interpolation is to estimate a value of (a function or series) between two known values. Extrapolation is to estimate (a value of a variable outside a known range) from values within a known range by assuming that the estimated value follows logically from the known values. The natural neighbours (clearly British) of any node are those in the neighbouring Voronoi cells, or equivalently, those to which the node is connected by the sides of Delaunay triangles (see cover graphic).

A Voronoi diagram of a set of "sites" (points) is a collection of regions that divide up the plane. Each region corresponds to one of the sites, and all the points in one region are closer to the corresponding site than to any other site. All of the Voronoi regions are convex polygons. Some of them are infinite -- these correspond to the sites on the convex hull. The boundary between two adjacent regions is a line segment, and the line that contains it is the perpendicular bisector of the segment joining the two sites. Usually, Voronoi regions meet three at at time at Voronoi points. If three sites determine Voronoi regions that meet at a Voronoi point, the circle through those three sites is centered at that Voronoi point, and there are no other sites in the circle. The Delaunay triangulation is the dual of the Voronoi diagram. To build the Delaunay triangulation, draw a line segment between any two sites whose Voronoi regions share an edge. This procedure splits the convex hull of the sites into triangles. Please see cover graphic for example.

The interpolating property of the Earth's crust make the procedures presented in this paper ideally suited for gridding irregularly spaced geophysical data. Natural neighbour interpolation describes a unique set of nodes that define the 'neighourhood' of a point in the plane. A set of nautural neighbours reflect the distance between nodes when their distance is relatively large in some part and when the distribution is highly anisotropic (having properties that differ according to the direction of measurement). The size and shape of the region that can influence any point will adapt naturally to the local variation in node density. Natural neighbour interpolation uses the areas of Voronoi cells to determine these weights.

Two examples are presented to demonstrate the power of this technique: i) Combining Global and regional scale models on a sphere; and ii) From subduction zones to spherical harmonic models of the Mantle . With this approach is is possible to build 2-D and 3-D with an enormous variation in cell sizes and complex geometry. Fine detail may be imposed anywhere by simply adding more nodes locally and letting the Delaunay tessellation (form a mosaic pattern) produce the cells. Also presented is an algorithm, called "walking triangle," used to find the triangle containing a given point.

Planar Point Location -- Bibliography (1986 - 1995)

@inproceedings{acg-cdctd-87
author = M. J. Atallah and R. Cole and M. T. Goodrich"
title = Cascading divide-and-conquer: {A technique for designing parallel algorithms"
booktitle = Proc. 28th Annu. IEEE Sympos. Found. Comput. Sci."
year = 1987
pages = 151--160"
keywords = cascading divide-and-conquer algorithm parallel algorithms concurrent-read- exclusive-write parallel random-access machine"
precedes = acg-cdctd-89"
annote = Triangulation of a polygon using $O(\log n)$ time and
$O(n)$ processors."
abstract = Techniques are presented for parallel
divide-and-conquer resulting in improved parallel
algorithms for a number of problems including
intersection detection trapezoidal decomposition
(hence polygon triangulation) and planar point
location (hence Voronoi diagram construction).
Efficient parallel algorithms are also given for
fractional cascading three-dimensional maxima two-set dominance counting and visibility from a point. All of the algorithms run in O(log n) time with either a linear or sublinear number of processors in the
concurrent-read- exclusive-write parallel random-access machine model. 20 refs."

@article{cj-nrdpp-92
author = S. W. Cheng and R. Janardan"
title = New results on dynamic planar point location"
journal = SIAM J. Comput."
volume = 21
year = 1992
pages = 972--999"
keywords = point location dynamization"

author = Y.-J. Chiang and F. P. Preparata and R. Tamassia"
title = A Unified Approach to Dynamic Point Location Ray Shooting and Shortest Paths in Planar Maps"
booktitle = Proc. 4th ACM-SIAM Sympos. Discrete Algorithms"
year = 1993
pages = 44--53"
update = 93.05 smid"

@inproceedings{ct-dtmpp-91
author = Y.-J. Chiang and R. Tamassia"
title = Dynamization of the trapezoid method for planar point location"
booktitle = Proc. 7th Annu. ACM Sympos. Comput. Geom."
year = 1991
pages = 61--70"
keywords = point location subdivision dynamizing data structures"
precedes = ct-dtmpp-92"

@article{ct-dtmpp-92
author = Y.-J. Chiang and R. Tamassia"
title = Dynamization of the trapezoid method for planar point location in monotone subdivision"
journal = Internat. J. Comput. Geom. Appl."
volume = 2
number = 3
year = 1992
pages = 311--333"
keywords = planar point location dynamic data structures monotone subdivision analysis of algorithm"
succeeds = ct-dtmpp-91"

@article{cz-opabd-90
author = R. Cole and O. Zajicek"
title = An optimal parallel algorithm for building a data structure for planar point location"
journal = J. Parallel Distrib. Comput."
volume = 8
year = 1990
pages = 280--285"
update = 93.05 devillers"

@article{lp-pbppl-89
author = D. T. Lee and F. P. Preparata"
title = Parallel batched planar point location on the {CCC "
journal = Inform. Process. Lett."
volume = 33
year = 1989
pages = 175--179"
update = 93.05 devillers"

@article{p-pplr-90
author = F. P. Preparata"
title = Planar point location revisited"
journal = Internat. J. Found. Comput. Sci."
,@techreport{pt-fdppl-87
author = F. P. Preparata and R. Tamassia"
title = A fully dynamic planar point location technique"
type = Report"
number = UILU-ENG-87-2266"
institution = Coordinated Sci. Lab. Univ. Illinois"
year = 1987
keywords = point location data structuring on-line"
precedes = pt-dpplo-89"

@inproceedings{pt-dpplo-89
author = F. P. Preparata and R. Tamassia"
title = Dynamic planar point location with optimal query time"
booktitle = Proc. 6th Sympos. Theoret. Aspects Comput. Sci."
series = Lecture Notes in Computer Science"
volume = 349
publisher = Springer-Verlag"
year = 1989
pages = 84--95"
keywords = subdivisions point location on-line two-dimensional"
succeeds = pt-fdppl-87"
precedes = pt-dpplo-90"

@article{pt-dpplo-90
author = F. P. Preparata and R. Tamassia"
title = Dynamic planar point location with optimal query time"
journal = Theoret. Comput. Sci."
volume = 74
year = 1990
pages = 95--114"
succeeds = pt-dpplo-89"

@inproceedings{pt-fdtpl-88
author = F. P. Preparata and R. Tamassia"
title = Fully dynamic techniques for point location and transitive closure in planar structures"
booktitle = Proc. 29th Annu. IEEE Sympos. Found. Comput. Sci."
year = 1988
pages = 558--567"
comments = In two journal papers: \cite{tp-dmpda-90 ,
\cite{pt-fdplm-89 "

@inproceedings{rs-orpac-87
author = John H. Reif and Sandeep Sen"
title = Optimal Randomized Parallel Algorithms For Computational Geometry"
booktitle = Proceedings of the 1987 International Conference on Parallel Processing."
publisher = Pennsylvania State Univ Press"
year = 1987
pages = 270--277"
keywords = randomized parallel algorithms computational geometry concurrent-read exclusive write parallel random-access machines"
annote = Including polygon triangulation."
abstract = The authors present parallel algorithms for some
fundamental problems in computational geometry which
have optimal running time with very high probability
(approaching 1 as n yields infinity). These include
planar point location triangulation trapezoidal
decomposition 3-D maxima two- set dominance counting, and range-counting among others. Most of these algorithms run on the CREW PRAM model with O(n)
processors and execute in O(log n) time. Most of these
algorithms use a novel data structure called the nested plane sweep tree. Random splitting is used very
effectively to divide and conquer on the plane raising the possibility of extending this technique to higher dimensions. 20 refs."

@article{rs-orpac-92
author = John H. Reif and Sandeep Sen"
title = Optimal Randomized Parallel Algorithms For Computational Geometry"
journal = Algorithmica"
volume = 7
number = 1
year = 1992
pages = 91--117"
keywords = randomized algorithms parallel algorithms crew pram computational geometry"
update = 93.09 rote"
abstract = We present parallel algorithms for some fundamental
problems in computational geometry which have a running time of O(log n) using n processors with very high probability (approaching 1 as n yields infinity ). These include planar-point location triangulation and trapezoidal decomposition. We also present optimal
algorithms for three-dimensional maxima and two-set
dominance counting by an application of integer
sorting. Most of these algorithms run on a CREW PRAM
model and have optimal processor-time product which
improve on the previously best-known algorithms of
Atallah and Goodrich  for these problems. The crux
of these algorithms is a useful data structure which
emulates the plane- sweeping paradigm used for
sequential algorithms. We extend some of the techniques used by Reischuk  and Reif and Valiant  for flashso@article{st-pplps-86
author = N. Sarnak and R. E. Tarjan"
title = Planar point location using persistent search trees"
journal = Commun. ACM"
volume = 29
year = 1986
pages = 669--679"

@article{s-facdt-87
author = S. W. Sloan"
title = A fast algorithm for constructing {Delaunay triangulations in the plane"
journal = Advances in Engineering Software"
volume = 9
number = 1
month = jan
year = 1987
pages = 34--55"
annote = Your basic incremental switching algorithm. Points are
sorted into bins to increase speed of point location
(which walks the triangulation starting at the
previous triangle.)"
abstract = The primary purpose of the paper is to present a
FORTRAN 77 code for the efficient triangulation of
planar data points as might be required for the mesh
generation of a finite element analysis. The author
describes in detail the various features of the
algorithm and compares them to existing routines. The
efficiency of the algorithm is demonstrated for
randomly distributed data points where the number of
points varies between 100 and 10,000. It is shown that
the proposed scheme outperforms existing routines for
both large and small data sets. -G. E. Bell St.
Andrews U"

@article{t-irmdp-91
author = R. Tamassia"
title = An incremental reconstruction method for dynamic planar point location"
journal = Inform. Process. Lett."
volume = 37
year = 1991
pages = 79--83"
keywords = point location triangulation computational geometry planar point location"
abstract = We present a fully dynamic technique for point
location in triangulations that allows a tradeoff
between query and update time and can be used in
conjunction with any of the known static point location data structures. Let $S$ be a triangulation whose current number of vertices is $n$. We show that for any smooth nondecreasing integer function $b(n)$ with $2\le b(n) \le \sqrt n$ there exists a dynamic point location data structure for $S$ with space $O(n)$, query time $O((\log^2n)/\log b(n))$ and update time $O((\log^2n)b(n)/\log b(n))$. (Author abstract) 15 Refs."

@inproceedings{tv-opatc-89
author = R. Tamassia and J. S. Vitter"
title = Optimal Parallel Algorithms for Transitive Closure and Point Location in Planar Structures"
booktitle = Proc. ACM Sympos. Parallel Algorithms Architect."
year = 1989
pages = 399--408"
keywords = graph drawing"
update = 93.09 tamassia"