A unified Pythagorean hodograph approach to the medial axis transform and offset approximation

AbstractAlgorithms based on Pythagorean hodographs (PH) in the Euclidean plane and in Minkowski space share common goals, the main one being rationality of offsets of planar domains. However, only separate interpolation techniques based on these curves can be found in the literature. It was recently revealed that rational PH curves in the Euclidean plane and in Minkowski space are very closely related. In this paper, we continue the discussion of the interplay between spatial MPH curves and their associated planar PH curves from the point of view of Hermite interpolation. On the basis of this approach we design a new, simple interpolation algorithm. The main advantage of the unifying method presented lies in the fact that it uses, after only some simple additional computations, an arbitrary algorithm for interpolation using planar PH curves also for interpolation using spatial MPH curves. We present the functionality of our method for G1 Hermite data; however, one could also obtain higher order algorithms

Continuous Medial Models in Two-Sample Statistics of Shape

In questions of statistical shape analysis, the foremost is how such shapes should be represented. The number of parameters required for a given accuracy and the types of deformation they can express directly influence the quality and type of statistical inferences one can make. One example is a medial model, which represents a solid object using a skeleton of a lower dimension and naturally expresses intuitive changes such as "bending", "twisting", and "thickening". In this dissertation I develop a new three-dimensional medial model that allows continuous interpolation of the medial surface and provides a map back and forth between the boundary and its medial axis. It is the first such model to support branching, allowing the representation of a much wider class of objects than previously possible using continuous medial methods. A measure defined on the medial surface then allows one to write integrals over the boundary and the object interior in medial coordinates, enabling the expression of important object properties in an object-relative coordinate system. I show how these properties can be used to optimize correspondence during model construction. This improved correspondence reduces variability due to how the model is parameterized which could potentially mask a true shape change effect. Finally, I develop a method for performing global and local hypothesis testing between two groups of shapes. This method is capable of handling the nonlinear spaces the shapes live in and is well defined even in the high-dimension, low-sample size case. It naturally reduces to several well-known statistical tests in the linear and univariate cases

Validity Determination for MAT Surface Representation

A vital issue to consider when exploiting the medial axis transform (MAT) as an object representation in its own right is whether the object recovered by the conversion process is a valid object. The general criterion for validity is that the recovered object should have no self-intersections. While any curve or surface could potentially be the medial axis (MA) of an object, the assignment of distances from the MA to the boundary can cause intersections in the resulting boundary. Boundary intersections can be caused by MAT points which are within a small "-neighborhood of each other, in which case the MAT is considered to be locally invalid. Additionally, two MAT points which are distant from each other may generate a self-intersection of the boundary, which makes the MAT globally invalid. While it is extremely difficult to determine global validity even after the boundary has been computed, it is conceivable that one could determine without first computing the boundary whether an MAT ..