mitigation Package

outlier_masking Module

Created on Jan 22, 2015

author: seehars

seek.mitigation.outlier_masking.getMask(tod, multiplicator)

Construct outlier mask based on user-specified rejection criterion.

Parameters:

• tod – TOD
• multiplicator – the number of standard deviation from the mean beyond which the data point is considered an outlier

Returns:

outlier mask

seek.mitigation.outlier_masking.rm_rfi(ctx)

Remove RFI by outlier rejection of a number of measurements in the same healpix pixel.

Parameters:ctx – context Returns:outlier mask

sum_threshold Module

Created on Jan 21, 2015

author: jakeret

seek.mitigation.sum_threshold.binary_mask_dilation(mask, struct_size_0, struct_size_1)

Dilates the mask.

Parameters:

• mask – original mask
• struct_size_0 – dilation parameter
• struct_size_1 – dilation parameter

Returns:

dilated mask

seek.mitigation.sum_threshold.get_di_kwrags(struct_size_0=3, struct_size_1=3)

Creates a dict with the dilation keywords.

Parameters:

• struct_size_0 – struct size in axis=0
• struct_size_1 – struct size in axis=1

Returns:

dictionary with the dilation keywords

seek.mitigation.sum_threshold.get_rfi_mask(tod, mask=None, chi_1=35000, eta_i=[0.5, 0.55, 0.62, 0.75, 1], normalize_standing_waves=True, suppress_dilation=False, plotting=True, sm_kwargs=None, di_kwargs=None)

Computes a mask to cover the RFI in a data set.

Parameters:

• data – array containing the signal and RFI
• mask – the initial mask
• chi_1 – First threshold
• eta_i – List of sensitivities
• normalize_standing_waves – whether to normalize standing waves
• suppress_dilation – if true, mask dilation is suppressed
• plotting – True if statistics plot should be displayed
• sm_kwargs – smoothing key words
• di_kwargs – dilation key words

Return mask:

the mask covering the identified RFI

seek.mitigation.sum_threshold.get_sm_kwargs(kernel_m=40, kernel_n=20, sigma_m=7.5, sigma_n=15)

Creates a dict with the smoothing keywords.

Parameters:

• kernel_m – kernel window size in axis=1
• kernel_n – kernel window size in axis=0
• sigma_m – kernel sigma in axis=1
• sigma_n – kernel sigma in axis=0

Returns:

dictionary with the smoothing keywords

seek.mitigation.sum_threshold.get_sumthreshold_kwargs(params)

Creates the smoothing and dilation kwargs from a params objects.

Parameters:params – the params object containing the configuration Returns:smoothing and dilation kwargs

seek.mitigation.sum_threshold.normalize(data, mask)

Simple normalization of standing waves: subtracting the median over time for each frequency.

Parameters:

• data – data
• mask – mask

Returns:

normalized data

seek.mitigation.sum_threshold.rm_rfi(ctx)

Call the main SumThreshold routine.

Parameters:ctx – context Returns:SumThreshold RFI mask.

sum_threshold_utils Module

Created on Dec 21, 2015

author: jakeret

seek.mitigation.sum_threshold_utils.get_stats(rfi, rfi_mask)

Returns the stats needed to compute a ROC curve.

Parameters:

• rfi – array containing the RFI pixels
• rfi_mask – boolean array that masks the RFI

Returns rl, ml, il:  

count of rfi pixels, count of masked pixels, count of intersecting pixels

seek.mitigation.sum_threshold_utils.plot_data(data, ax, title, vmin=None, vmax=None, cb=True, norm=None, extent=None, cmap=None)

Plot TOD.

seek.mitigation.sum_threshold_utils.plot_dilation(st_mask, mask, dilated_mask)

Plot mask and dilation.

seek.mitigation.sum_threshold_utils.plot_moments(data)

Plot standard divation and mean of data.

seek.mitigation.sum_threshold_utils.plot_steps(data, st_mask, smoothed_data, res, eta)

Plot individual steps of SumThreshold.