“…It should be pointed out that the problem of analyzing data spatially distributed over irregularly shaped two-dimensional domains has recently attracted an increasing interest, and other regularized least-square smoothers have been proposed that can tackle this issue, such as bivariate splines over triangulations [see, e.g., Lai and Schumaker, 2007, Guillas and Lai, 2010, Ettinger et al, 2012, Lai and Wang, 2013, soap film smoothing [Wood et al, 2008], and low-rank thin-plate spline approximations Ranalli, 2007, Scott-Hayward et al, 2014]. All these methods have isotropic and stationary regularizing terms; bivariate splines over triangulations can include high order derivatives.…”