Localized bases for kernel spaces on the unit sphere

Approximation/interpolation from spaces of positive definite or conditionally positive definite kernels is an increasingly popular tool for the analysis and synthesis of scattered data, and is central to many meshless methods. For a set of $N$ scattered sites, the standard basis for such a space utilizes $N$ \emph{globally} supported kernels; computing with it is prohibitively expensive for large $N$. Easily computable, well-localized bases, with "small-footprint" basis elements - i.e., elements using only a small number of kernels -- have been unavailable. Working on \sphere, with focus on the restricted surface spline kernels (e.g. the thin-plate splines restricted to the sphere), we construct easily computable, spatially well-localized, small-footprint, robust bases for the associated kernel spaces. Our theory predicts that each element of the local basis is constructed by using a combination of only $\mathcal{O}((\log N)^2)$ kernels, which makes the construction computationally cheap. We prove that the new basis is $L_p$ stable and satisfies polynomial decay estimates that are stationary with respect to the density of the data sites, and we present a quasi-interpolation scheme that provides optimal $L_p$ approximation orders. Although our focus is on $\mathbb{S}^2$, much of the theory applies to other manifolds - $\mathbb{S}^d$, the rotation group, and so on. Finally, we construct algorithms to implement these schemes and use them to conduct numerical experiments, which validate our theory for interpolation problems on $\mathbb{S}^2$ involving over one hundred fifty thousand data sites.Comment: This article supersedes arXiv:1111.1013 "Better bases for kernel spaces," which proved existence of better bases for various kernel spaces. This article treats a smaller class of kernels, but presents an algorithm for constructing better bases and demonstrates its effectiveness with more elaborate examples. A quasi-interpolation scheme is introduced that provides optimal linear convergence rate

Computational Modeling for 3D Data Reconstruction of Solar Coronal Magnetic Fields

The solar corona -- the atmosphere of the Sun -- supports a plethora of magnetic structures that constitute major drivers of space weather. Such a rope of magnetic field can twist and break to launch the plasma locked inside it in a Coronal Mass Ejection, or CME. The analysis and prediction of events of this nature is a major goal of the heliophysics community, and reconstructing the current state of the coronal magnetic field is a central component. Determining the shape and structure of the magnetic fields arcing through the solar corona is a form of nonlinear inverse problem. The solution strategy for reconstructing a given magnetic structure is composed of three separate routines: a parameterized forward model, which can generate the sorts of structures native to the corona; a methodology by which realistic observations of such models can be synthesized; and an optimization framework that searches many computed models for the one whose synthetic observations best match true observations of the real coronal structure. The goal of this work is to provide a solid foundation for the construction of coronal inversion frameworks. First, a full working version of such a framework is presented, using existing forms of all three components. For a variety of reasons, the parameterized forward model is deemed the component of this full inverse model most needing improvement. Subsequent work lays out a brand new methodology for computation of the forward model. This new forward model is then inspected for its ability to be implemented in an improved inversion framework. The result is a construction that addresses many of the traditional issues around coronal field modeling.</p

Self-Similar Vector Fields

We propose statistically self-similar and rotation-invariant models for vector fields, study some of the more significant properties of these models, and suggest algorithms and methods for reconstructing vector fields from numerical observations, using the same notions of self-similarity and invariance that give rise to our stochastic models. We illustrate the efficacy of the proposed schemes by applying them to the problems of denoising synthetic flow phantoms and enhancing flow-sensitive magnetic resonance imaging (MRI) of blood flow in the aorta. In constructing our models and devising our applied schemes and algorithms, we rely on two fundamental notions. The first of these, referred to as "innovation modelling" in the thesis, is the principle —applicable both analytically and synthetically— of reducing complex phenomena to combinations of simple independent components or "innovations". The second fundamental idea is that of "invariance", which indicates that in the absence of any distinguishing factor, two equally valid models or solutions should be given equal consideration

Compression, pose tracking, and halftoning

In this thesis, we discuss image compression, pose tracking, and halftoning. Although these areas seem to be unrelated at first glance, they can be connected through video coding as application scenario. Our first contribution is an image compression algorithm based on a rectangular subdivision scheme which stores only a small subsets of the image points. From these points, the remained of the image is reconstructed using partial differential equations. Afterwards, we present a pose tracking algorithm that is able to follow the 3-D position and orientation of multiple objects simultaneously. The algorithm can deal with noisy sequences, and naturally handles both occlusions between different objects, as well as occlusions occurring in kinematic chains. Our third contribution is a halftoning algorithm based on electrostatic principles, which can easily be adjusted to different settings through a number of extensions. Examples include modifications to handle varying dot sizes or hatching. In the final part of the thesis, we show how to combine our image compression, pose tracking, and halftoning algorithms to novel video compression codecs. In each of these four topics, our algorithms yield excellent results that outperform those of other state-of-the-art algorithms.In dieser Arbeit werden die auf den ersten Blick vollkommen voneinander unabhängig erscheinenden Bereiche Bildkompression, 3D-Posenschätzung und Halbtonverfahren behandelt und im Bereich der Videokompression sinnvoll zusammengeführt. Unser erster Beitrag ist ein Bildkompressionsalgorithmus, der auf einem rechteckigen Unterteilungsschema basiert. Dieser Algorithmus speichert nur eine kleine Teilmenge der im Bild vorhandenen Punkte, während die restlichen Punkte mittels partieller Differentialgleichungen rekonstruiert werden. Danach stellen wir ein Posenschätzverfahren vor, welches die 3D-Position und Ausrichtung von mehreren Objekten anhand von Bilddaten gleichzeitig verfolgen kann. Unser Verfahren funktioniert bei verrauschten Videos und im Falle von Objektüberlagerungen. Auch Verdeckungen innerhalb einer kinematischen Kette werden natürlich behandelt. Unser dritter Beitrag ist ein Halbtonverfahren, das auf elektrostatischen Prinzipien beruht. Durch eine Reihe von Erweiterungen kann dieses Verfahren flexibel an verschiedene Szenarien angepasst werden. So ist es beispielsweise möglich, verschiedene Punktgrößen zu verwenden oder Schraffuren zu erzeugen. Der letzte Teil der Arbeit zeigt, wie man unseren Bildkompressionsalgorithmus, unser Posenschätzverfahren und unser Halbtonverfahren zu neuen Videokompressionsalgorithmen kombinieren kann. Die für jeden der vier Themenbereiche entwickelten Verfahren erzielen hervorragende Resultate, welche die Ergebnisse anderer moderner Verfahren übertreffen

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

Mathematical and Numerical Aspects of Dynamical System Analysis

From Preface: This is the fourteenth time when the conference “Dynamical Systems: Theory and Applications” gathers a numerous group of outstanding scientists and engineers, who deal with widely understood problems of theoretical and applied dynamics. Organization of the conference would not have been possible without a great effort of the staff of the Department of Automation, Biomechanics and Mechatronics. The patronage over the conference has been taken by the Committee of Mechanics of the Polish Academy of Sciences and Ministry of Science and Higher Education of Poland. It is a great pleasure that our invitation has been accepted by recording in the history of our conference number of people, including good colleagues and friends as well as a large group of researchers and scientists, who decided to participate in the conference for the first time. With proud and satisfaction we welcomed over 180 persons from 31 countries all over the world. They decided to share the results of their research and many years experiences in a discipline of dynamical systems by submitting many very interesting papers. This year, the DSTA Conference Proceedings were split into three volumes entitled “Dynamical Systems” with respective subtitles: Vibration, Control and Stability of Dynamical Systems; Mathematical and Numerical Aspects of Dynamical System Analysis and Engineering Dynamics and Life Sciences. Additionally, there will be also published two volumes of Springer Proceedings in Mathematics and Statistics entitled “Dynamical Systems in Theoretical Perspective” and “Dynamical Systems in Applications”

Social work with airports passengers

Social work at the airport is in to offer to passengers social services. The main methodological position is that people are under stress, which characterized by a particular set of characteristics in appearance and behavior. In such circumstances passenger attracts in his actions some attention. Only person whom he trusts can help him with the documents or psychologically

Commonwealth of Independent States aerospace science and technology, 1992: A bibliography with indexes

This bibliography contains 1237 annotated references to reports and journal articles of Commonwealth of Independent States (CIS) intellectual origin entered into the NASA Scientific and Technical Information System during 1992. Representative subject areas include the following: aeronautics, astronautics, chemistry and materials, engineering, geosciences, life sciences, mathematical and computer sciences, physics, social sciences, and space sciences