schlafly

import numpy as np
import scipy.ndimage
import scipy.interpolate
import scipy.special
from astropy.convolution import convolve, Box2DKernel


# some function we can oversample
def make_fake_image(oversample):
    pts = np.linspace(-10, 10, 10 * oversample + 1)