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.

Azimuthal jittered sampling of bandlimited functions in the two-dimensional Fourier transform and the Hankel transform domains

Controlling T 2 blurring in 3D RARE arterial spin labeling acquisition through optimal combination of variable flip angles and k-space filtering.

Sampling theorems associated with offset linear canonical transform by polar coordinates

Golden‐ratio Rotated Stack‐of‐stars Acquisition for Improved Volumetric MRI

SYMPLECTIC PSEUDOSPECTRAL TIME-DOMAIN SCHEME FOR SOLVING TIME-DEPENDENT SCHRODINGER EQUATION

Joint Angular-Frequency Distribution Estimation via Spatial-temporal Sparse Sampling and Low-rank Matrix Recovery

Mathematical Equivalence of Zero-Padding Interpolation and Circular Sampling Theorem Interpolation with Implications for Direct Fourier Image Reconstruction.

Orthogonal Array Sampling for Monte Carlo Rendering

Reduced Aliasing Artifacts Using Shaking Projection K-Space Sampling Trajectory

Stochastic Oscillation Undersampling to Reduce Frequency Aliasing Artifacts in the Magnetic Resonance Imaging

Diagnosing space telescope misalignment and jitter using stellar images

Spectral solver for Cauchy problems in polar coordinates using discrete Hankel transforms

Hankel Spectrum Analysis: A novel signal decomposition method and its geophysical applications

Utilizing Bochners Theorem for Constrained Evaluation of Missing Fourier Data

Anti-aliasing Weighting Functions for Single-Slice Helical CT

Aliasing Effects for Samples of Spin Random Fields on the Sphere

Uncertainty quantification for sparse Fourier recovery

Rational Approximation of Golden Angles: Accelerated Reconstructions for Radial MRI

Fourier-Based Approach To Interpolation In Single-Slice Helical Computed Tomography

Joint Angular-Frequency Distribution Estimation of Incoherently Distributed Wideband Sources via Low-Rank Matrix Recovery