Barnes interpolation (Barnes objective analysis)

Smoothing interpolation of unstructured N-dimensional data using Barnes objective analysis
433 Downloads
Updated 20 Sep 2016

View License

BARNESN Barnes smoothing interpolation of unstructured data
Vq = BARNESN(X, V, Xv) returns the smoothing interpolation of
D-dimensional observations V(X) at query points Xq. Query points Xq are
created by meshing the vectors in the cell array Xv that define the
grid in each dimension. Smoothing interpolation is performed using
the Koch form of Barnes objective analysis [2]. Roughly, (in 2D) the
interpolated value (vq) at gridpoint (xq, yq) is determined as a
weighted-sum of the values (v) at data points (x, y), based on the
gaussian weighting function exp(-r^2 / s / g^j), where r is the
euclidian distance from (xq, yq) to (x, y), s is the Gaussian Variance,
and g is the Convergence Parameter.
-
Bibliography:
[1] Barnes, Stanley L. "Mesoscale objective map analysis using weighted
time-series observations." (1973)
[2] Koch, Steven E., Mary DesJardins, and Paul J. Kocin. "An
interactive Barnes objective map analysis scheme for use with
satellite and conventional data." Journal of Climate and Applied
Meteorology 22.9 (1983): 1487-1503.
[3] Daley, Roger. Atmospheric data anlysis. No. 2. Cambridge University
Press, 1993.
-
See documentation/help for more information.

Cite As

Lena Bartell (2024). Barnes interpolation (Barnes objective analysis) (https://www.mathworks.com/matlabcentral/fileexchange/58937-barnes-interpolation-barnes-objective-analysis), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2016a
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Time Series in Help Center and MATLAB Answers
Acknowledgements

Inspired by: Barnes objective analysis

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes
1.3.0.0

Update to appropriately handle the situation wherein some data points lie outside the grid (previously this just produced NaNs).

1.2.0.0

Fixed some errors in the algorithm implementation and implemented optimal parameters as default. Expanded functionality from 2D to N-dimensions.

1.1.0.0

Corrected typo in example implementation

1.0.0.0