Robust smoothing of gridded data in one and higher dimensions with missing values

Robust smoothing of gridded data in one and higher dimensions with missing values