Methods

Alard & Lupton

Alard and Lupton [alard1998] introduced a method to simultaneously fit a background \(B_{kg}\) and a convolution kernel \(K\) that will minimize the difference between a reference image \(R\) and an another science image \(I\).

\[I \approx R \otimes K + B_{kg}\]

This method assumes that the convolution kernel can be approximated by a linear combination of fixed Gaussians, where the linear coefficients are left free to vary. In addition, each Gaussian is modulated with a polynomial of an arbitrary degree.

\[\begin{split}K &= \sum_i a_i B_i \\ &= \sum_n a_n \times \left[ \exp \left(- \frac{u^2 + v^2}{2 \sigma_n^2} \right) \sum_{d_n^x} \sum_{d_n^y} u^{d_n^x} v^{d_n^y} \right]\end{split}\]

Note

The centre, width and orientation of the Gaussians are fixed beforehand as well as the number of Gaussians to use in the expansion.

In OIS this is specified with a list of dictionaries, one for each gaussian we want to use. Below is an example of a basis with 3 Gaussians:

gausslist=[{center: (5, 5), sx: 2., sy: 2., modPolyDeg: 3},
           {sx: 1.0, sy: 2.5, modPolyDeg: 1},
           {sx: 3.0, sy: 1.0},]

Bramich

The method developed in [bramich2008] modifies Alard-Lupton making each pixel of the kernel an independent parameter to fit. This is equivalent as having a vector basis consisting of Delta kernels. It can also simultaneously fit a polynomial background.

\[K = \sum_n a_n \times \delta_{nn'}\]

This method does not make assumptions on the kernel shape and can thus model completely arbitrary kernels. It can also correct for small translations between the images.

While more accurate, this method is computationally more expensive than Alard & Lupton’s.

Warning

Since each pixel is treated independently, a 11 by 11 kernel will have 121 free parameters just for the kernel. It grows quadratically with the kernel side. This needs to be taken in consideration for large kernels.

Adaptive Bramich

Like Bramich, this method also treats each pixel independently, but it will also multiply each pixel by a polynomial on the coordinates of the image.

This requires a special type of convolution where the kernel varies point to point in the image.

It is especially suited for situations where the PSF varies significatively across the image.

The method is described in more detail in [miller2008].

Warning

Just like Bramich, the number of free parameters scales quadratically with the kernel side. Furthermore, the degree of the polynomial multiplies the number of parameters by (deg + 1) * (deg + 2) / 2. This needs to be taken in consideration for large kernels.