Conformal geometry processing

This thesis introduces fundamental equations and numerical methods for manipulating surfaces in three dimensions via conformal transformations. Conformal transformations are valuable in applications because they naturally preserve the integrity of geometric data. To date, however, there has been no clearly stated and consistent theory of conformal transformations that can be used to develop general-purpose geometry processing algorithms: previous methods for computing conformal maps have been restricted to the flat two-dimensional plane, or other spaces of constant curvature. In contrast, our formulation can be used to produce---for the first time---general surface deformations that are perfectly conformal in the limit of refinement. It is for this reason that we commandeer the title Conformal Geometry Processing.

The main contribution of this thesis is analysis and discretization of a certain time-independent Dirac equation, which plays a central role in our theory. Given an immersed surface, we wish to construct new immersions that (i) induce a conformally equivalent metric and (ii) exhibit a prescribed change in extrinsic curvature. Curvature determines the potential in the Dirac equation; the solution of this equation determines the geometry of the new surface. We derive the precise conditions under which curvature is allowed to evolve, and develop efficient numerical algorithms for solving the Dirac equation on triangulated surfaces.

From a practical perspective, this theory has a variety of benefits: conformal maps are desirable in geometry processing because they do not exhibit shear, and therefore preserve textures as well as the quality of the mesh itself. Our discretization yields a sparse linear system that is simple to build and can be used to efficiently edit surfaces by manipulating curvature and boundary data, as demonstrated via several mesh processing applications. We also present a formulation of Willmore flow for triangulated surfaces that permits extraordinarily large time steps and apply this algorithm to surface fairing, geometric modeling, and construction of constant mean curvature (CMC) surfaces.

Conformal structures and twistors in the paravector model of spacetime

Some properties of the Clifford algebras Cl-3,Cl-0, Cl-1,Cl-3, Cl-4,Cl-1 similar or equal to C circle times Cl-1,Cl-3 and Cl-2,Cl-4 are presented, and three isomorphisms between the Dirac-Clifford algebra C circle times Cl-1,Cl-3 and Cl-4,Cl-1 are exhibited, in order to construct conformal maps and twistors, using the paravector model of spacetime. The isomorphism between the twistor space inner product isometry group SU( 2,2) and the group $pin(+)(2,4) is also investigated, in the light of a suitable isomorphism between C circle times Cl-1,Cl-3 and Cl-4,Cl-1. After reviewing the conformal spacetime structure, conformal maps are described in Minkowski spacetime as the twisted adjoint representation of $ pin(+)(2,4), acting on paravectors. Twistors are then presented via the paravector model of Clifford algebras and related to conformal maps in the Clifford algebra over the Lorentzian R-4,(1) spacetime.We construct twistors in Minkowski spacetime as algebraic spinors associated with the Dirac-Clifford algebra C circle times Cl-1,Cl-3 using one lower spacetime dimension than standard Clifford algebra formulations, since for this purpose, the Clifford algebra over R-4,R-1 is also used to describe conformal maps, instead of R-2,(4). Our formalism sheds some new light on the use of the paravector model and generalizations.

Hausdorff dimension of non-hyperbolic repellers. I: Maps with holes

This is the first paper in a two-part series devoted to studying the Hausdorff dimension of invariant sets of non-uniformly hyperbolic, non-conformal maps. Here we consider a general abstract model, that we call piecewise smooth maps with holes. We show that the Hausdorff dimension of the repeller is strictly less than the dimension of the ambient manifold. Our approach also provides information on escape rates and dynamical dimension of the repeller.

Broadband electromagnetic cloaking with smart metamaterials.

The ability to render objects invisible with a cloak that fits all objects and sizes is a long-standing goal for optical devices. Invisibility devices demonstrated so far typically comprise a rigid structure wrapped around an object to which it is fitted. Here we demonstrate smart metamaterial cloaking, wherein the metamaterial device not only transforms electromagnetic fields to make an object invisible, but also acquires its properties automatically from its own elastic deformation. The demonstrated device is a ground-plane microwave cloak composed of an elastic metamaterial with a broad operational band (10-12 GHz) and nearly lossless electromagnetic properties. The metamaterial is uniform, or perfectly periodic, in its undeformed state and acquires the necessary gradient-index profile, mimicking a quasi-conformal transformation, naturally from a boundary load. This easy-to-fabricate hybrid elasto-electromagnetic metamaterial opens the door to implementations of a variety of transformation optics devices based on quasi-conformal maps.

Exam 1998

Exam and solutions in LaTex

Exam 1998 pdf

Exam and solutions in PDF

Radiobiological model comparison of 3D conformal radiotherapy and IMRT plans for the treatment of prostate cancer

The main aim of radiotherapy is to deliver a dose of radiation that is high enough to destroy the tumour cells while at the same time minimising the damage to normal healthy tissues. Clinically, this has been achieved by assigning a prescription dose to the tumour volume and a set of dose constraints on critical structures. Once an optimal treatment plan has been achieved the dosimetry is assessed using the physical parameters of dose and volume. There has been an interest in using radiobiological parameters to evaluate and predict the outcome of a treatment plan in terms of both a tumour control probability (TCP) and a normal tissue complication probability (NTCP). In this study, simple radiobiological models that are available in a commercial treatment planning system were used to compare three dimensional conformal radiotherapy treatments (3D-CRT) and intensity modulated radiotherapy (IMRT) treatments of the prostate. Initially both 3D-CRT and IMRT were planned for 2 Gy/fraction to a total dose of 60 Gy to the prostate. The sensitivity of the TCP and the NTCP to both conventional dose escalation and hypo-fractionation was investigated. The biological responses were calculated using the Källman S-model. The complication free tumour control probability (P+) is generated from the combined NTCP and TCP response values. It has been suggested that the alpha/beta ratio for prostate carcinoma cells may be lower than for most other tumour cell types. The effect of this on the modelled biological response for the different fractionation schedules was also investigated.