I try to fill NA cells in a SpatialPixelsDataFrame based on mean values from neighbouring cells. I wonder if there is any automatized approach equivalent to the focal function in terra (see here).