estimate_dxdy

astropop.polarimetry.estimate_dxdy(x, y, steps=[100, 30, 5, 3], bins=30, dist_limit=100)

Estimate the displacement between the two beams.

To compute the displacement between the ordinary and extraordinary beams, this function computes the most common distances between the sources in image, using clipped histograms around the peak.

Parameters:

x, y: array_like

Arrays of x and y positions of the sources.

steps: list of int

Number of pixels, around the mode, to clip the histogram in each step.

bins: int

Number of bins to use in the histogram.

dist_limit: float

Maximum distance between the pairs of sources to consider.