LOCI

class psfsubtraction.fitpsf.LOCI(psfbase, image=None)

Bases: psfsubtraction.fitpsf.BasePSFFitter

LOCI fitter (locally optimized combination of images)

The loci algorithm was introduced in the following paper Lafreniere et al. 2007, ApJ, 660, 770.

The default parameters in this fitter are chosen similar to the shape of the regions used in that paper.

Attributes Summary

dilation_region
sector_phi
sector_radius
sector_radius_inner
sector_radius_n

Methods Summary

findbase(region)	Return all bases that are not masked in any pixel in region
fitpsfcoeff(image1d, psfbase)	solve a linear algebra system for the best PSF
optregion(region, indpsf)
regions()	Generate a function that generates sector regions

Attributes Documentation

dilation_region = 10

sector_phi = 12

sector_radius

sector_radius_inner = 0

sector_radius_n = 10

Methods Documentation

findbase(region)

Return all bases that are not masked in any pixel in region

fitpsfcoeff(image1d, psfbase)

solve a linear algebra system for the best PSF

Parameters:

image1d : array in 1 dim

psfbase : array in [M,N]

M = number of pixels in flattened image N = number of images that form the space of potential PSFs

Returns:

psf_coeff : array in 1 dim

Coefficients for a linear combination of psfbase elements that that give the optimal PSF.

Raises:

ValueError : If given masked data, because numpy.linalg.solve would silently

use the “values behind the mask”.

optregion(region, indpsf)

regions()

Generate a function that generates sector regions

A pixel is included in a region, if the pixel center falls within the region boundaries.

This function makes use of the following fitter attributes, which have to be set to use this function:

fitter.sector_radius : np.array

boundaries for sector elements in pixels.

fitter.sector_phi : int or quantity

If this is an int it sets the number of sectors that make up a full circle. If this is an astropy.quantity it is interpreted as the boundaries of the angular bins. It should cover the range from 0 to 2 pi (or 360 deg, if units is degrees).

fitter.sector_center : tuple or None

x, y position of the center of all sectors (in pixel coordinates). None selects the center of the input image.

Returns:

regions : generator

sector regions