Angular Fast Smooth - NaN Tolerant

Fast smoothing algorithm for angles that is tolerant of NaNs in the data
Updated 26 Aug 2015

View License

nanfastsmoothAngle is a function of the form y = fastsmoothAngle(x,w,type,tol) where y becomes a smoothed version of x using a moving average function of width w where NaNs in the data are ignored. This function is specifically for angles and avoids the usual problem of trying to average angles such as 10 and 350 where a normal average would be 180 but the angular average is 0 or 360.
The argument "type" determines the smooth type:
If type=1, rectangular (sliding-average or boxcar)
If type=2, triangular (2 passes of sliding-average)
If type=3, pseudo-Gaussian (3 passes of sliding-average)
The argument "tol" controls the amount of tolerance to NaNs allowed between 0and 1. A value of zero means that if the window has any NaNs in it then the output is set as NaN. A value of 1 allows any number of NaNs in the window and will still give an answer for the smoothed signal. A value of 0.5 means that there must be at least half real values in the window for the output to be valid.

The start and end of the file are treated as if there are NaNs beyond the dataset. As such the behaviour depends on the value of 'tol' as described above. With 'tol' set at 0.5 the smoothed signal will start and end at the same time as the orgional signal. However it's accuracy will be reduced and the moving average will become more and more one-sided as the beginning and end is approached.

fastsmoothAngle(Y,w,type) smooths with tol = 0.5.
fastsmoothAngle(Y,w) smooths with type = 1 and tol = 0.5

Written using Matlab 2013b but should work with older releases

This file is heavily based on the fastsmooth function by T. C. O'Haver, May, 2008.

Cite As

Greg Pittam (2024). Angular Fast Smooth - NaN Tolerant (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2013b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Inspired by: Fast smoothing function

Community Treasure Hunt

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

Start Hunting!
Version Published Release Notes