Examples

Compute standard deviations in circles of diameter 10arcsec around grid pixels spaced every 4 x 5 pixels and anchored at pixel [30, 40], and use linear interpolation to compute values at non-grid-pixels:

Compute median of the absolute deviations from the median values using the z-score/Chauvenet algorithm, by fixing the maximum z-score to determine outliers to 5. Use cubic interpolation to compute values for non-grid-point pixels. Use a rectangular region of dimensions 5arcsec x 20arcsec centered on each grid point as the region in which to include pixels for the computation of stats for that grid point.