Abstract:We generalize N‐rooks, jittered, and (correlated) multi‐jittered sampling to higher dimensions by importing and improving upon a class of techniques called orthogonal arrays from the statistics literature. Renderers typically combine or "pad" a collection of lower‐dimensional (e.g. 2D and 1D) stratified patterns to form higher‐dimensional samples for integration. This maintains stratification in the original dimension pairs, but looses it for all other dimension pairs. For truly multi‐dimensional integrands like those in rendering, this increases variance and deteriorates its rate of convergence to that of pure random sampling. Care must therefore be taken to assign the primary dimension pairs to the dimensions with most integrand variation, but this complicates implementations. We tackle this problem by developing a collection of practical, in‐place multi‐dimensional sample generation routines that stratify points on all t‐dimensional and 1‐dimensional projections simultaneously. For instance, when t=2, any 2D projection of our samples is a (correlated) multi‐jittered point set. This property not only reduces variance, but also simplifies implementations since sample dimensions can now be assigned to integrand dimensions arbitrarily while maintaining the same level of stratification. Our techniques reduce variance compared to traditional 2D padding approaches like PBRT's (0,2) and Stratified samplers, and provide quality nearly equal to state‐of‐the‐art QMC samplers like Sobol and Halton while avoiding their structured artifacts as commonly seen when using a single sample set to cover an entire image. While in this work we focus on constructing finite sampling point sets, we also discuss potential avenues for extending our work to progressive sequences (more suitable for incremental rendering) in the future.

An Improvement of Normal Approximation of Randomized Orthogonal Array Sampling Designs

A Uniform Bound on a Combinatorial Central Limit Theorem for Randomized Orthogonal Array Sampling Designs

Approximation of Bounds on Mixed Level Orthogonal Arrays

A multivariate central limit theorem for randomized orthogonal array sampling designs in computer experiments

Approximate Truthful Mechanism Design For Two-Dimensional Orthogonal Knapsack Problem

Nested Lattice Sampling: A New Sampling Scheme Derived by Randomizing Nested Orthogonal Arrays

Approximation Algorithms for D-optimal Design

Unbiased orthogonal designs

Efficient Bounds and Estimates for Canonical Angles in Randomized Subspace Approximations

A GENERAL THEORY FOR ORTHOGONAL ARRAY BASED LATIN HYPERCUBE SAMPLING

Randomized low-rank approximations beyond Gaussian random matrices

Using a Non-Commutative Bernstein Bound to Approximate Some Matrix Algorithms in the Spectral Norm

Orthogonal Array Sampling for Monte Carlo Rendering

Improved NP-Hardness of Approximation for Orthogonality Dimension and Minrank

Generalized resolution for orthogonal arrays

Sampling projections in the uniform norm

Randomized least-squares with minimal oversampling and interpolation in general spaces

Exact discretization, tight frames and recovery via D-optimal designs

Adaptive and Oblivious Randomized Subspace Methods for High-Dimensional Optimization: Sharp Analysis and Lower Bounds

Three-Parameter Approximations of Sums of Locally Dependent Random Variables via Stein's Method