Practical Box Splines for Reconstruction on the Body Centered Cubic Lattice

We introduce a family of box splines for efficient, accurate, and smooth reconstruction of volumetric data sampled on the body-centered cubic (BCC) lattice, which is the favorable volumetric sampling pattern due to its optimal spectral sphere packing property. First, we construct a box spline based on the four principal directions of the BCC lattice that allows for a linear $C ^{ 0 }$ reconstruction. Then, the design is extended for higher degrees of continuity. We derive the explicit piecewise polynomial representations of the $C ^{ 0 }$ and $C ^{ 2 }$ box splines that are useful for practical reconstruction applications. We further demonstrate that approximation in the shift-invariant space—generated by BCC-lattice shifts of these box splines—is twice as efficient as using the tensor-product B-spline solutions on the Cartesian lattice (with comparable smoothness and approximation order and with the same sampling density). Practical evidence is provided demonstrating that the BCC lattice not only is generally a more accurate sampling pattern, but also allows for extremely efficient reconstructions that outperform tensor-product Cartesian reconstructions

Toward high-quality gradient estimation on regular lattices

Abstract—In this paper, we present two methods for accurate gradient estimation from scalar field data sampled on regular lattices. The first method is based on the multidimensional Taylor series expansion of the convolution sum and allows us to specify design criteria such as compactness and approximation power. The second method is based on a Hilbert space framework and provides a minimum error solution in the form of an orthogonal projection operating between two approximation spaces. Both methods lead to discrete filters, which can be combined with continuous reconstruction kernels to yield highly accurate estimators as compared to the current state of the art. We demonstrate the advantages of our methods in the context of volume rendering of data sampled on Cartesian and Body-Centered Cubic lattices. Our results show significant qualitative and quantitative improvements for both synthetic and real data, while incurring a moderate preprocessing and storage overhead. Index Terms—Approximation theory, Taylor series expansion, normal reconstruction, orthogonal projection, body-centered cubic lattice, box splines. Ç

A Practical Box Spline Compendium

Box splines provide smooth spline spaces as shifts of a single generating function on a lattice and so generalize tensor-product splines. Their elegant theory is laid out in classical papers and a summarizing book. This compendium aims to succinctly but exhaustively survey symmetric low-degree box splines with special focus on two and three variables. Tables contrast the lattices, supports, analytic and reconstruction properties, and list available implementations and code.Comment: 15 pages, 10 figures, 8 table

Symmetric box-splines on root lattices

AbstractRoot lattices are efficient sampling lattices for reconstructing isotropic signals in arbitrary dimensions, due to their highly symmetric structure. One root lattice, the Cartesian grid, is almost exclusively used since it matches the coordinate grid; but it is less efficient than other root lattices. Box-splines, on the other hand, generalize tensor-product B-splines by allowing non-Cartesian directions. They provide, in any number of dimensions, higher-order reconstructions of fields, often of higher efficiency than tensored B-splines. But on non-Cartesian lattices, such as the BCC (Body-Centered Cubic) or the FCC (Face-Centered Cubic) lattice, only some box-splines and then only up to dimension three have been investigated.This paper derives and completely characterizes efficient symmetric box-spline reconstruction filters on all irreducible root lattices that exist in any number of dimensions n≥2 (n≥3 for Dn and Dn∗ lattices). In all cases, box-splines are constructed by convolution using the lattice directions, generalizing the known constructions in two and three variables. For each box-spline, we document the basic properties for computational use: the polynomial degree, the continuity, the linear independence of shifts on the lattice and optimal quasi-interpolants for fast approximation of fields

Fast space-variant elliptical filtering using box splines

The efficient realization of linear space-variant (non-convolution) filters is a challenging computational problem in image processing. In this paper, we demonstrate that it is possible to filter an image with a Gaussian-like elliptic window of varying size, elongation and orientation using a fixed number of computations per pixel. The associated algorithm, which is based on a family of smooth compactly supported piecewise polynomials, the radially-uniform box splines, is realized using pre-integration and local finite-differences. The radially-uniform box splines are constructed through the repeated convolution of a fixed number of box distributions, which have been suitably scaled and distributed radially in an uniform fashion. The attractive features of these box splines are their asymptotic behavior, their simple covariance structure, and their quasi-separability. They converge to Gaussians with the increase of their order, and are used to approximate anisotropic Gaussians of varying covariance simply by controlling the scales of the constituent box distributions. Based on the second feature, we develop a technique for continuously controlling the size, elongation and orientation of these Gaussian-like functions. Finally, the quasi-separable structure, along with a certain scaling property of box distributions, is used to efficiently realize the associated space-variant elliptical filtering, which requires O(1) computations per pixel irrespective of the shape and size of the filter.Comment: 12 figures; IEEE Transactions on Image Processing, vol. 19, 201

On the Treatment of Field Quantities and Elemental Continuity in FEM Solutions

As the finite element method (FEM) and the finite volume method (FVM), both traditional and high-order variants, continue their proliferation into various applied engineering disciplines, it is important that the visualization techniques and corresponding data analysis tools that act on the results produced by these methods faithfully represent the underlying data. To state this in another way: the interpretation of data generated by simulation needs to be consistent with the numerical schemes that underpin the specific solver technology. As the verifiable visualization literature has demonstrated: visual artifacts produced by the introduction of either explicit or implicit data transformations, such as data resampling, can sometimes distort or even obfuscate key scientific features in the data. In this paper, we focus on the handling of elemental continuity, which is often only C0 continuous or piecewise discontinuous, when visualizing primary or derived fields from FEM or FVM simulations. We demonstrate that traditional data handling and visualization of these fields introduce visual errors. In addition, we show how the use of the recently proposed line-SIAC filter provides a way of handling elemental continuity issues in an accuracy-conserving manner with the added benefit of casting the data in a smooth context even if the representation is element discontinuous