A two-stage filter for smoothing multivariate noisy data on unstructured grids☆

Abstract Experimental data as well as numerical simulations are very often affected by “noise”, random fluctuations that distort the final output or the intermediate products of a numerical process. In this paper, a new method for smoothing data given on nonstructured grids is proposed. It takes advantage of the smoothing properties of kernel-weighted averaging and least-squares techniques. The weighted averaging is performed following a Shepard-like procedure, where Gaussian kernels are employed, while least-squares fitting is reduced to the use of very few basis functions so as to improve smoothness, though at the price of interpolation accuracy. Once we have defined an n -point grid and want to make a smoothed fit at a given grid point, this method reduces to consideration of all m -point stencils (for a given value of m where m n ) that include that point, and to make a least-squares fitting, for that particular point, in each of those stencils; finally, the various results, thus obtained, are weight-averaged, the weights being inversely proportional to the distance of the point to the middle of the stencil. Then this process is repeated for all the grid points so as to obtain the smoothing of the input function. Though this method is generalized for the multivariate case, one-, two-, and five-dimensional test cases are shown as examples of the performance of this method.

[1]  Tom Lyche,et al.  Mathematical methods in computer aided geometric design , 1989 .

[2]  T. A. Foley,et al.  Three-stage interpolation to scattered data , 1984 .

[3]  P. Lancaster.  Curve and surface fitting , 1986 .

[4]  A. F. de Vos,et al.  A synopsis of the smoothing formulae associated with the Kalman filter , 1993 .

[5]  D. Truhlar,et al.  A dual‐level Shepard interpolation method for generating potential energy surfaces for dynamics calculations , 1995 .

[6]  O. Talagrand,et al.  A new reduced-order adaptive filter for state estimation in high-dimensional systems , 1997, Autom..

[7]  I. K. Grain.  Computer interpolation and contouring of two-dimensional data: A review , 1970 .

[8]  Gung Feng,et al.  Data smoothing by cubic spline filters , 1998, IEEE Trans. Signal Process..

[9]  J. G. Kalbfleisch.  Probability and Statistical Inference , 1977 .

[10]  C. Chui,et al.  Approximation Theory II , 1976 .

[11]  P. Alfeld.  Scattered data interpolation in three or more variables , 1989 .

[12]  Paul Suetens,et al.  An algorithm for surface reconstruction from planar contours using smoothing splines , 1988 .

[13]  Lennart Ljung,et al.  Fast time-invariant implementations for linear least-squares smoothing filters , 1978 .

[14]  Surface smoothing with finite elements , 1985 .

[15]  W. J. Whitten.  The Use of Multi-Dimensional Cubic Spline Functions for Regression and Smoothing , 1971, Aust. Comput. J..

[16]  L. Schumaker,et al.  Multivariate scattered data fitting , 1996 .

[17]  H. Akima,et al.  On estimating partial derivatives for bivariate interpolation of scattered data , 1984 .

[18]  W. Cleveland,et al.  Locally Weighted Regression: An Approach to Regression Analysis by Local Fitting , 1988 .

[19]  M. A. Collins,et al.  Molecular potential energy surfaces by interpolation , 1994 .

[20]  P. Lancaster,et al.  Surfaces generated by moving least squares methods , 1981 .

[21]  Non-linear smoothing algorithm for multi-dimensional dynamic systems , 1989 .

[22]  Nobuhisa Kashiwagi,et al.  On use of the Kalman filter for spatial smoothing , 1993 .

[23]  L. Schumaker.  Fitting surfaces to scattered data , 1976 .

[24]  Robert E. Barnhill,et al.  Multistage trivariate surfaces , 1984 .

[25]  O. Talagrand,et al.  On the design of a stable adaptive filter for state estimation in high dimensional systems , 2001, Autom..

[26]  The whittaker-vondrak method of data smoothing as a numerical filter , 1981 .

[27]  G. Bierman.  Factorization methods for discrete sequential estimation , 1977 .

[28]  R. Wyatt,et al.  Quantum Wave Packet Dynamics with Trajectories , 1999 .

[29]  Surface smoothing with B-spline finite elements , 1992 .

[30]  Approximation and Smoothing of Surfaces In (P+1)-Dimensional Spaces , 1982 .

[31]  Eli Sternbers.  Proceedings of the first U.S. National Congress of Applied Mechanics , 1952 .

[32]  Paul Dierckx,et al.  Smoothing scattered data with a monotone Powell-Sabin spline surface , 2005, Numerical Algorithms.

[33]  G. Alexits.  Approximation theory , 1983 .

[34]  C. Chui,et al.  Approximation Theory VI , 1990 .

[35]  I. K Crain,et al.  Treatment of non-equispaced two-dimensional data with a digital computer , 1967 .

[36]  J. Meinguet.  Multivariate interpolation at arbitrary points made simple , 1979 .

[37]  Robert E. Wyatt.  Quantum wavepacket dynamics with trajectories: wavefunction synthesis along quantum paths , 1999 .

[38]  T. A. Foley,et al.  Shepard's Interpolation for Solution-Adaptive Methods , 1993 .

[39]  Attila Askar,et al.  Quantum fluid dynamics in the Lagrangian representation and applications to photodissociation problems , 1999 .

[40]  R. Wyatt,et al.  Quantum wave packet dynamics with trajectories: Implementation with adaptive Lagrangian grids , 2000 .

[41]  Hongmin Lu,et al.  Surface approximation by spline smoothing and generalized cross-validation , 1992 .

[42]  E. Dubois,et al.  Design of multidimensional finite-wordlength FIR and IIR filters by simulated annealing , 1995 .