Box spline 2026-10-14
For a spanning direction matrix in , the box spline is the density of the pushforward of uniform measure on under . Equivalently,
Its support is the zonotope . It is piecewise polynomial of total degree at most . For integer directions, splitting each integration interval in half yields its binary refinement symbol
up to the monomial specifying the choice of origin. If is the smallest number of directions whose removal leaves a nonspanning set, the usual box spline smoothness criterion gives regularity. Repeated directions increase smoothness while enlarging support.
In a convenient cardinal parameter normalization, take the unit -cube and project by the sum of its coordinates. Its pushforward density is the unit-direction univariate box spline
Equivalently, with at its discontinuity points,
Here there are convolution factors; some conventions instead index the spline by its degree . The half-open convention for assigns knot values consistently and does not change the pushforward measure.
Support and positivity. A sum of numbers in lies in , so the closed support is finite. For , the interval of with has positive length. Induction in the convolution recurrence therefore gives in the interior. For , : strict positivity includes the open support interval, not its endpoints.
Partition of unity. The translates of cover the line exactly once under the half-open convention. Nonnegativity permits interchange of sum and integral, and induction gives
Thus , a locally finite sum. Each also has integral one, although partition of unity does not mean its maximum value is one.
For , repeated integration gives the explicit truncated-power form
On each unit knot interval this is a polynomial of degree at most , and on it is , so the degree is attained. Each truncated power has continuous derivatives; the derivative of order jumps at the first knot. Therefore the dependence on the projected box dimension is
For the pieces are constants and the function is discontinuous at the knots. If geometric width is measured under orthogonal projection onto the unit diagonal instead, the projected coordinate is , so the physical support width is ; its unit-mass density is . These are different scalings of the same cardinal spline, not different continuity or degree claims.
For a bivariate box spline, choose a spanning direction matrix in the plane and project the unit cube by . Define its normalization by
The support is the zonotope , and the polynomial pieces generally have total degree . The grid comes from the chosen direction families and their lattice translates. Two independent directions, with repetitions, give tensor-product B-splines on a parallelogram grid. Three families such as give a three-direction triangular grid; an affine change of basis can make it the usual equilateral triangular grid. The four families give the criss-cross square grid and the quadratic four-direction box spline. Changing the direction matrix and repeating directions therefore changes the grid, support and smoothness while retaining the box-projection construction.
In a regular region, moving one control point adds its displacement multiplied by a translate of the twice-smoothed four-direction box spline. Measure parameter distances in original grid-edge units, with the affected control point at the origin. The convex hull of the binary mask offsets is
After successive binary refinements the physical displacements are scaled by . Their possible accumulated offsets therefore form . The box spline description confirms that this zonotope is the actual support, not merely a loose enclosure. Equivalently,
It is an octagon with vertices . Its enclosing square is six original grid intervals wide in each direction. Cutting four right triangles of area two from this square gives area parameter-square units. Outside this octagon the control-point displacement has no effect; its closure is the support, although the basis function can vanish on its boundary.
The regular influence footprint is this width-six, area-28 parameter octagon. For a nonplanar control mesh, the affected part of the spatial surface is the image of the octagon, not necessarily a planar octagon of the same metric size. Near an extraordinary subdivision vertex, use the actual refinement neighborhood instead of assuming a globally regular lattice footprint.
Take two copies of each direction . The resulting box spline has total polynomial degree six and continuity: there are eight directions and a rank-deficient remaining set has at most two, so the removal number is six. With a centred origin its support is the zonotope
Its centred binary symbol is
This symbol factors as with , giving an efficient quincunx subdivision implementation. The box spline criterion is also described in the primary research paper 4–8 Subdivision by Velho and Zorin.
Zonotope 2026-10-14
A zonotope is a Minkowski sum of finitely many line segments, equivalently an affine map image of a cube. For vectors , the set is the support of the associated box spline. A centred version uses instead. In two dimensions, sums of differently directed segments give centrally symmetric polygons; repeated directions lengthen their corresponding edges without adding new edge directions.