On the conditions for the coincidence of two cubic Bézier curves

AbstractIn a recent article, Wang et al. [2] derive a necessary and sufficient condition for the coincidence of two cubic Bézier curves with non-collinear control points. The condition reads that their control points must be either coincident or in reverse order. We point out that this uniqueness of the control points for polynomial cubics is a straightforward consequence of a previous and more general result of Barry and Patterson, namely the uniqueness of the control points for rational Bézier curves. Moreover, this uniqueness applies to properly parameterized polynomial curves of arbitrary degree

On automatic tuning of basis functions in Bezier method

A transition from the fixed basis in Bezier's method to some class of base functions is proposed. A parameter vector of a basis function is introduced as additional information. This achieves a more universal form of presentation and analytical description of geometric objects as compared to the non-uniform rational B-splines (NURBS). This enables control of basis function parameters including control points, their weights and node vectors. This approach can be useful at the final stage of constructing and especially local modification of compound curves and surfaces with required differential and shape properties; it also simplifies solution of geometric problems. In particular, a simple elimination of discontinuities along local spline curves due to automatic tuning of basis functions is demonstrated

Gardener's spline curve

In the well-known gardener’s construction of the ellipse we replace the two foci by a finite set of points in the plane, that results in a G spline curve that consists of elliptic arcs, if the set contains at least three non-collinear points. An algorithm is provided for the specification of these elliptic arcs, along with their control point based representation

Robust multigrid methods for Isogeometric discretizations applied to poroelasticity problems

El análisis isogeométrico (IGA) elimina la barrera existente entre elementos finitos (FEA) y el diseño geométrico asistido por ordenador (CAD). Debido a esto, IGA es un método novedoso que está recibiendo una creciente atención en la literatura y recientemente se ha convertido en tendencia. Muchos esfuerzos están siendo puestos en el diseño de solvers eficientes y robustos para este tipo de discretizaciones. Dada la optimalidad de los métodos multimalla para elementos finitos, la aplicación de estosmétodos a discretizaciones isogeométricas no ha pasado desapercibida. Nosotros pensamos firmemente que los métodos multimalla son unos candidatos muy prometedores a ser solvers eficientes y robustos para IGA y por lo tanto en esta tesis apostamos por su aplicación. Para contar con un análisis teórico para el diseño de nuestros métodos multimalla, el análisis local de Fourier es propuesto como principal análisis cuantitativo. En esta tesis, a parte de considerar varios problemas escalares, prestamos especial atención al problema de poroelasticidad, concretamente al modelo cuasiestático de Biot para el proceso de consolidación del suelo. Actualmente, el diseño de métodos multimalla robustos para problemas poroelásticos respecto a parámetros físicos o el tamaño de la malla es un gran reto. Por ello, la principal contribución de esta tesis es la propuesta de métodos multimalla robustos para discretizaciones isogeométricas aplicadas al problema de poroelasticidad.La primera parte de esta tesis se centra en la construcción paramétrica de curvas y superficies dado que estas técnicas son la base de IGA. Así, la definición de los polinomios de Bernstein y curvas de Bézier se presenta como punto de partida. Después, introducimos los llamados B-splines y B-splines racionales no uniformes (NURBS) puesto que éstas serán las funciones base consideradas en nuestro estudio.La segunda parte trata sobre el análisis isogeométrico propiamente dicho. En esta parte, el método isoparamétrico es explicado al lector y se presenta el análisis isogeométrico de algunos problemas. Además, introducimos la formulación fuerte y débil de los problemas anteriores mediante el método de Galerkin y los espacios de aproximación isogeométricos. El siguiente punto de esta tesis se centra en los métodos multimalla. Se tratan las bases de los métodos multimalla y, además de introducir algunos métodos iterativos clásicos como suavizadores, también se introducen suavizadores por bloques como los métodos de Schwarz multiplicativos y aditivos. Llegados a esta parte, nos centramos en el LFA para el diseño de métodos multimalla robustos y eficientes. Además, se explican en detalle el análisis estándar y el análisis basado en ventanas junto al análisis de suavizadores por bloques y el análisis para sistemas de ecuaciones en derivadas parciales.Tras introducir las discretizaciones isogeométricas, los métodos multimalla y el LFA como análisis teórico, nuestro propósito es diseñar métodos multimalla eficientes y robustos respecto al grado polinomial de los splines para discretizaciones isogeométricas de algunos problemas escalares. Así, mostramos que el uso de métodos multimalla basados en suavizadores de tipo Schwarz multiplicativo o aditivo produce buenos resultados y factores de convergencia asintóticos robustos. La última parte de esta tesis está dedicada al análisis isogeométrico del problema de poroelasticidad. Para esta tarea, se introducen el modelo de Biot y su discretización isogeométrica. Además, presentamos una novedosa estabilización de masa para la formulación de dos campos de las ecuaciones de Biot que elimina todas las oscilaciones no físicas en la aproximación numérica de la presión. Después, nos centramos en dos tipos de solvers para estas ecuaciones poroelásticas: Solvers desacoplados y solvers monolíticos. En el primer grupo, le dedicamos una especial atención al método fixed-stress y a un método iterativo propuesto por nosotros que puede ser aplicado de forma automática a partir de la estabilización de masa ya mencionada.Por otro lado, realizamos un análisis de von Neumann para este método iterativo aplicado al problema de Terzaghi y demostramos su estabilidad y convergencia para los pares de elementos Q1 Q1, Q2 Q1 y Q3 Q2 (con suavidad global C1). Respecto al grupo de solvers monolíticos, nosotros proponemos métodos multimalla basados en suavizadores acoplados y desacoplados. En esta parte, métodosIsogeometric analysis (IGA) eliminates the gap between finite element analysis (FEA) and computer aided design (CAD). Due to this, IGA is an innovative approach that is receiving an increasing attention in the literature and it has recently become a trending topic. Many research efforts are being devoted to the design of efficient and robust solvers for this type of discretization. Given the optimality of multigrid methods for FEA, the application of these methods to IGA discretizations has not been unnoticed. We firmly think that they are a very promising approach as efficient and robust solvers for IGA and therefore in this thesis we are concerned about their application. In order to give a theoretical support to the design of multigrid solvers, local Fourier analysis (LFA) is proposed as the main quantitative analysis. Although different scalar problems are also considered along this thesis, we make a special focus on poroelasticity problems. More concretely, we focus on the quasi-static Biot's equations for the soil consolidation process. Nowadays, it is a very challenging task to achieve robust multigrid solvers for poroelasticity problems with respect physical parameters and/or the mesh size. Thus, the main contribution of this thesis is to propose robust multigrid methods for isogeometric discretizations applied to poroelasticity problems. The first part of this thesis is devoted to the introduction of the parametric construction of curves and surfaces since these techniques are the basis of IGA. Hence, with the definition of Bernstein polynomials and B\'ezier curves as a starting point, we introduce B-splines and non-uniform rational B-splines (NURBS) since these will be the basis functions considered for our numerical experiments. The second part deals with the isogeometric analysis. In this part, the isoparametric approach is explained to the reader and the isogeometric analysis of some scalar problems is presented. Hence, the strong and weak formulations by means of Galerkin's method are introduced and the isogeometric approximation spaces as well. The next point of this thesis consists of multigrid methods. The basics of multigrid methods are explained and, besides the presentation of some classical iterative methods as smoothers, block-wise smoothers such as multiplicative and additive Schwarz methods are also introduced. At this point, we introduce LFA for the design of efficient and robust multigrid methods. Furthermore, both standard and infinite subgrids local Fourier analysis are explained in detail together with the analysis for block-wise smoothers and the analysis for systems of partial differential equations. After the introduction of isogeometric discretizations, multigrid methods as our choice of solvers and LFA as theoretical analysis, our goal is to design efficient and robust multigrid methods with respect to the spline degree for IGA discretizations of some scalar problems. Hence, we show that the use of multigrid methods based on multiplicative or additive Schwarz methods provide a good performance and robust asymptotic convergence rates. The last part of this thesis is devoted to the isogeometric analysis of poroelasticity. For this task, Biot's model and its isogeometric discretization are introduced. Moreover, we present an innovative mass stabilization of the two-field formulation of Biot's equations that eliminates all the spurious oscillations in the numerical approximation of the pressure. Then, we deal with two types of solvers for these poroelastic equations: Decoupled and monolithic solvers. In the first group we devote special attention to the fixed-stress split method and a mass stabilized iterative scheme proposed by us that can be automatically applied from the mass stabilization formulation mentioned before. In addition, we perform a von Neumann analysis for this iterative decoupled solver applied to Terzaghi's problem and demonstrate that it is stable and convergent for pairs Q1-Q1, Q2-Q1 and Q3-Q2 (with global smoothness C1). Regarding the group of monolithic solvers, we propose multigrid methods based on coupled and decoupled smoothers. Coupled additive Schwarz methods are proposed as coupled smoothers for isogeometric Taylor-Hood elements. More concretely, we propose a 51-point additive Schwarz method for the pair Q2-Q1. In the last part, we also propose to use an inexact version of the fixed-stress split algorithm as decoupled smoother by applying iterations of different additive Schwarz methods for each variable. For the latter approach, we consider the pairs of elements Q2-Q1 and Q3-Q2 (with global smoothness C1). Finally, thanks to LFA we manage to design efficient and robust multigrid solvers for the Biot's equations and some numerical results are shown.<br /

A novel parallel algorithm for surface editing and its FPGA implementation

A thesis submitted to the University of Bedfordshire in partial fulfilment of the requirements for the degree of Doctor of PhilosophySurface modelling and editing is one of important subjects in computer graphics. Decades of research in computer graphics has been carried out on both low-level, hardware-related algorithms and high-level, abstract software. Success of computer graphics has been seen in many application areas, such as multimedia, visualisation, virtual reality and the Internet. However, the hardware realisation of OpenGL architecture based on FPGA (field programmable gate array) is beyond the scope of most of computer graphics researches. It is an uncultivated research area where the OpenGL pipeline, from hardware through the whole embedded system (ES) up to applications, is implemented in an FPGA chip. This research proposes a hybrid approach to investigating both software and hardware methods. It aims at bridging the gap between methods of software and hardware, and enhancing the overall performance for computer graphics. It consists of four parts, the construction of an FPGA-based ES, Mesa-OpenGL implementation for FPGA-based ESs, parallel processing, and a novel algorithm for surface modelling and editing. The FPGA-based ES is built up. In addition to the Nios II soft processor and DDR SDRAM memory, it consists of the LCD display device, frame buffers, video pipeline, and algorithm-specified module to support the graphics processing. Since there is no implementation of OpenGL ES available for FPGA-based ESs, a specific OpenGL implementation based on Mesa is carried out. Because of the limited FPGA resources, the implementation adopts the fixed-point arithmetic, which can offer faster computing and lower storage than the floating point arithmetic, and the accuracy satisfying the needs of 3D rendering. Moreover, the implementation includes Bézier-spline curve and surface algorithms to support surface modelling and editing. The pipelined parallelism and co-processors are used to accelerate graphics processing in this research. These two parallelism methods extend the traditional computation parallelism in fine-grained parallel tasks in the FPGA-base ESs. The novel algorithm for surface modelling and editing, called Progressive and Mixing Algorithm (PAMA), is proposed and implemented on FPGA-based ES’s. Compared with two main surface editing methods, subdivision and deformation, the PAMA can eliminate the large storage requirement and computing cost of intermediated processes. With four independent shape parameters, the PAMA can be used to model and edit freely the shape of an open or closed surface that keeps globally the zero-order geometric continuity. The PAMA can be applied independently not only FPGA-based ESs but also other platforms. With the parallel processing, small size, and low costs of computing, storage and power, the FPGA-based ES provides an effective hybrid solution to surface modelling and editing

CAD interface and framework for curve optimisation applications

Computer Aided Design is currently expanding its boundaries to include more design features in its processes. Design is identified as an iterative process converging to solutions satisfying a set of constraints. Its close relation with optimisation indicate that there is strong potential for the integration of optimisation and CAD. The problem addressed in this thesis lies in interfacing the geometric representation of design with other non-geometric aspects. The example of free-form curve modelling is taken to investigate such relationships. Assumptions are made that Optimisation is powered by Evolutionary Computing algorithms like Genetic Algorithms (GA). The geometric definition of curves is commonly supported by NURBS, whose construction constraints are defined locally at the data points. Here the NURBS formulation is used with GA in an attempt to provide complementary handles on the curves shape other than the usual data point coordinates and control points weights. Differential properties are used for optimising NURBS, Hermite interpolation allows for the definition of higher order constraints (tangent, normal, bi-normal) at data points. The assignment of parameter values at the data points, known as parameterisation also provides control of the curve’s shape. Curve optimisation is also performed at the geometric modelling level. Old mathematical theorems established by Frénet and further developed by other mathematicians provide means of defining a curve’s shape with it’s intrinsic equations. Such representation is possible by using Function Representation (F-rep) algebra available in the ACIS software. Frep allows more generic and exact means of interfacing with the curve’s geometry and new functionality for curve inspection and optimisation are proposed in this thesis. The integration of optimisation findings and CAD are documented in the definition of a framework. The framework architecture proposed reconstructs a new CAD environment from separate elements bolted together in a generic Application Programming Interface (API) named “Oli interface”. Functionality created to interface optimisation and CAD makes a requirement list of the work that both sides should undertake to achieve design optimisation in the CAD environment.EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Changing representation of curves and surfaces: exact and approximate methods

Το κύριο αντικείμενο μελέτης στην παρούσα διατριβή είναι η αλλαγή αναπαράστασης γεωμετρικών αντικειμένων από παραμετρική σε αλγεβρική (ή πεπλεγμένη) μορφή. Υπολογίζουμε την αλγεβρική εξίσωση παρεμβάλλοντας τους άγνωστους συντελεστές του πολυωνύμου δεδομένου ενός υπερσυνόλου των μονωνύμων του. Το τελευταίο υπολογίζεται απο το Newton πολύτοπο της αλγεβρικής εξίσωσης που υπολογίζεται από μια πρόσφατη μέθοδο πρόβλεψης του συνόλου στήριξης της εξίσωσης. H μέθοδος πρόβλεψης του συνόλου στήριξης βασίζεται στην αραιή (ή τορική) απαλοιφή: το πολύτοπο υπολογίζεται από το Newton πολύτοπο της αραιής απαλοίφουσας αν θεωρίσουμε την παραμετροποίηση ως πολυωνυμικό σύστημα. Στα μονώνυμα που αντιστοιχούν στα ακέραια σημεία του Newton πολυτόπου δίνονται τιμές ώστε να σχηματίσουν έναν αριθμητικό πίνακα. Ο πυρήνα του πίνακα αυτού, διάστασης 1 σε ιδανική περίπτωση, περιέχει τους συντελεστές των μονωνύμων στην αλγεβρική εξίσωση. Υπολογίζουμε τον πυρήνα του πίνακα είτε συμβολικά είτε αριθμητικά εφαρμόζοντας την μέθοδο του singular value decomposition (SVD). Προτείνουμε τεχνικές για να διαχειριστούμε την περίπτωση ενός πολυδιάστατου πυρήνα το οποίο εμφανίζεται όταν το προβλεπόμενο σύνολο στήριξης είναι ένα υπερσύνολο του πραγματικού. Αυτό δίνει έναν αποτελεσματικό ευαίσθητο-εξόδου αλγόριθμο υπολογισμού της αλγεβρικής εξίσωσης. Συγκρίνουμε διαφορετικές προσεγγίσεις κατασκευής του πίνακα μέσω των λογισμικών Maple και SAGE. Στα πειράματά μας χρησιμοποιήθηκαν ρητές καμπύλες και επιφάνειες καθώς και NURBS. Η μέθοδός μας μπορεί να εφαρμοστεί σε πολυώνυμα ή ρητές παραμετροποιήσεις επίπεδων καμπυλών ή (υπερ)επιφανειών οποιασδήποτε διάστασης συμπεριλαμβανομένων και των περιπτώσεων με παραμετροποίηση σεσημεία βάσης που εγείρουν σημαντικά ζητήματα για άλλες μεθόδους αλγεβρικοποίησης. Η μέθοδος έχει τον εξής περιορισμό: τα γεωμετρικά αντικείμενα πρέπει να αναπαριστώνται από βάσεις μονωνύμων που στην περίπτωση τριγωνομετρικών παραμετροποιήσεων θα πρέπει να μπορούν να μετασχηματιστούν σε ρητές συναρτήσεις. Επιπλέον η τεχνική που προτείνουμε μπορεί να εφαρμοστεί σε μη γεωμετρικά προβλήματα όπως ο υπολογισμόςτης διακρίνουσας ενός πολυωνύμου με πολλές μεταβλητές ή της απαλοίφουσας ενός συστήματος πολυωνύμων με πολλές μεταβλητές.The main object of study in our dissertation is the representation change of the geometric objects from the parametric form to implicit. We compute the implicit equation interpolating the unknown coefficients of the implicit polynomial given a superset of its monomials. The latter is derived from the Newton polytope of the implicit equation obtained by the recently developed method for support prediction. The support prediction method we use relies on sparse (or toric) elimination: the implicit polytope is obtained from the Newton polytope of the sparse resultant of the system in parametrization, represented as polynomials. The monomials that correspond to the lattice points of the Newton polytope are suitably evaluated to build a numeric matrix, ideally of corank 1. Its kernel contains their coefficients in the implicit equation. We compute kernel of the matrix either symbolically, or numerically, applying singular value decomposition (SVD). We propose techniques for handling the case of the multidimensional kernel space, caused by the predicted support being a superset of the actual. This yields an efficient, output-sensitive algorithm for computing the implicit equation. We compare different approaches for constructing the matrix in Maple and SAGE software. In our experiments we have used classical algebraic curves and surfaces as well as NURBS. Our method can be applied to polynomial or rational parametrizations of planar curves or (hyper)surfaces of any dimension including cases of parameterizations with base points which raise important issues for other implicitization methods. The method has its limits: geometric objects have to be presented using monomial basis; in the case of trigonometric parametrizations they have to be convertible to rational functions. Moreover, the proposed technique can be applied for nongeometric problems such as the computation of the discriminant of a multivariate polynomial or the resultant of a system of multivariate polynomials

Camera Calibration with Non-Central Local Camera Models

Kamerakalibrierung ist eine wichtige Grundvoraussetzung für viele Computer-Vision-Algorithmen wie Stereo-Vision und visuelle Odometrie. Das Ziel der Kamerakalibrierung besteht darin, sowohl die örtliche Lage der Kameras als auch deren Abbildungsmodell zu bestimmen. Das Abbildungsmodell einer Kamera beschreibt den Zusammenhang zwischen der 3D-Welt und der Bildebene. Aktuell werden häufig einfache globale Kamera-Modelle in einem Kalibrierprozess geschätzt, welcher mit vergleichsweise geringem Aufwand und einer großen Fehlertoleranz durchgeführt werden kann. Um das resultierende Kameramodell zu bewerten, wird in der Regel der Rückprojektionsfehler als Maß herangezogen. Jedoch können auch einfache Kameramodelle, die das Abbildungsverhalten von optischen Systemen nicht präzise beschreiben können, niedrige Rückprojektionsfehler erzielen. Dies führt dazu, dass immer wieder schlecht kalibrierte Kameramodelle nicht als solche identifiziert werden. Um dem entgegen zu wirken, wird in dieser Arbeit ein neues kontinuierliches nicht-zentrales Kameramodell basierend auf B-Splines vorgeschlagen. Dieses Abbildungsmodell ermöglicht es, verschiedene Objektive und nicht-zentrale Verschiebungen, die zum Beispiel durch eine Platzierung der Kamera hinter einer Windschutzscheibe entstehen, akkurat abzubilden. Trotz der allgemeinen Modellierung kann dieses Kameramodell durch einen einfach zu verwendenden Schachbrett-Kalibrierprozess geschätzt werden. Um Kalibrierergebnisse zu bewerten, wird anstelle des mittleren Rückprojektionsfehlers ein Kalibrier-Benchmark vorgeschlagen. Die Grundwahrheit des Kameramodells wird durch ein diskretes Sichtstrahlen-basiertes Modell beschrieben. Um dieses Modell zu schätzen, wird ein Kalibrierprozess vorgestellt, welches ein aktives Display als Ziel verwendet. Dabei wird eine lokale Parametrisierung für die Sichtstrahlen vorgestellt und ein Weg aufgezeigt, die Oberfläche des Displays zusammen mit den intrinsischen Kameraparametern zu schätzen. Durch die Schätzung der Oberfläche wird der mittlere Punkt-zu-Linien-Abstand um einen Faktor von mehr als 20 reduziert. Erst dadurch kann das so geschätzte Kameramodell als Grundwahrheit dienen. Das vorgeschlagene Kameramodell und die dazugehörigen Kalibrierprozesse werden durch eine ausführliche Auswertung in Simulation und in der echten Welt mithilfe des neuen Kalibrier-Benchmarks bewertet. Es wird gezeigt, dass selbst in dem vereinfachten Fall einer ebenen Glasscheibe, die vor der Kamera platziert ist, das vorgeschlagene Modell sowohl einem zentralen als auch einem nicht-zentralen globalen Kameramodell überlegen ist. Am Ende wird die Praxistauglichkeit des vorgeschlagenen Modells bewiesen, indem ein automatisches Fahrzeug kalibriert wird, das mit sechs Kameras ausgestattet ist, welche in unterschiedliche Richtungen zeigen. Der mittlere Rückprojektionsfehler verringert sich durch das neue Modell bei allen Kameras um den Faktor zwei bis drei. Der Kalibrier-Benchmark ermöglicht es in Zukunft, die Ergebnisse verschiedener Kalibrierverfahren miteinander zu vergleichen und die Genauigkeit des geschätzten Kameramodells mithilfe der Grundwahrheit akkurat zu bestimmen. Die Verringerung des Kalibrierfehlers durch das neue vorgeschlagene Kameramodell hilft die Genauigkeit weiterführender Algorithmen wie Stereo-Vision, visuelle Odometrie oder 3D-Rekonstruktion zu erhöhen

AutoGraff: towards a computational understanding of graffiti writing and related art forms.

The aim of this thesis is to develop a system that generates letters and pictures with a style that is immediately recognizable as graffiti art or calligraphy. The proposed system can be used similarly to, and in tight integration with, conventional computer-aided geometric design tools and can be used to generate synthetic graffiti content for urban environments in games and in movies, and to guide robotic or fabrication systems that can materialise the output of the system with physical drawing media. The thesis is divided into two main parts. The first part describes a set of stroke primitives, building blocks that can be combined to generate different designs that resemble graffiti or calligraphy. These primitives mimic the process typically used to design graffiti letters and exploit well known principles of motor control to model the way in which an artist moves when incrementally tracing stylised letter forms. The second part demonstrates how these stroke primitives can be automatically recovered from input geometry defined in vector form, such as the digitised traces of writing made by a user, or the glyph outlines in a font. This procedure converts the input geometry into a seed that can be transformed into a variety of calligraphic and graffiti stylisations, which depend on parametric variations of the strokes