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\).
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.
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.
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.